It has recently been suggested that the low-frequency c.136–14_136–13insC variant in high-mobility group A1 (HMGA1) may strongly contribute to insulin resistance and type 2 diabetes risk. In our study, we attempted to confirm that HMGA1 is a novel type 2 diabetes locus in French Caucasians. The gene was sequenced in 368 type 2 diabetic case subjects with a family history of type 2 diabetes and 372 normoglycemic control subjects without a family history of type 2 diabetes. None of the 41 genetic variations identified were associated with type 2 diabetes. The lack of association between the c.136–14_136–13insC variant and type 2 diabetes was confirmed in an independent French group of 4,538 case subjects and 4,015 control subjects and in a large meta-analysis of 16,605 case subjects and 46,179 control subjects. Finally, this variant had no effects on metabolic traits and was not involved in variations of HMGA1 and insulin receptor (INSR) expressions. The c.136–14_136–13insC variant was not associated with type 2 diabetes in individuals of European descent. Our study emphasizes the need to analyze a large number of subjects to reliably assess the association of low-frequency variants with the disease.

Early menopause (EM) affects up to 10% of the female population, reducing reproductive lifespan considerably. Currently, it constitutes the leading cause of infertility in the western world, affecting mainly those women who postpone their first pregnancy beyond the age of 30 years. The genetic aetiology of EM is largely unknown in the majority of cases. We have undertaken a meta-analysis of genome-wide association studies (GWASs) in 3493 EM cases and 13 598 controls from 10 independent studies. No novel genetic variants were discovered, but the 17 variants previously associated with normal age at natural menopause as a quantitative trait (QT) were also associated with EM and primary ovarian insufficiency (POI). Thus, EM has a genetic aetiology which overlaps variation in normal age at menopause and is at least partly explained by the additive effects of the same polygenic variants. The combined effect of the common variants captured by the single nucleotide polymorphism arrays was estimated to account for ∼30% of the variance in EM. The association between the combined 17 variants and the risk of EM was greater than the best validated non-genetic risk factor, smoking.

The male-to-female sex ratio at birth is constant across world populations with an average of 1.06 (106 male to 100 female live births) for populations of European descent. The sex ratio is considered to be affected by numerous biological and environmental factors and to have a heritable component. The aim of this study was to investigate the presence of common allele modest effects at autosomal and chromosome X variants that could explain the observed sex ratio at birth. We conducted a large-scale genome-wide association scan (GWAS) meta-analysis across 51 studies, comprising overall 114 863 individuals (61 094 women and 53 769 men) of European ancestry and 2 623 828 common (minor allele frequency >0.05) single-nucleotide polymorphisms (SNPs). Allele frequencies were compared between men and women for directly-typed and imputed variants within each study. Forward-time simulations for unlinked, neutral, autosomal, common loci were performed under the demographic model for European populations with a fixed sex ratio and a random mating scheme to assess the probability of detecting significant allele frequency differences. We do not detect any genome-wide significant (P < 5 × 10−8) common SNP differences between men and women in this well-powered meta-analysis. The simulated data provided results entirely consistent with these findings. This large-scale investigation across ∼115 000 individuals shows no detectable contribution from common genetic variants to the observed skew in the sex ratio. The absence of sex-specific differences is useful in guiding genetic association study design, for example when using mixed controls for sex-biased traits.

Common diseases such as type 2 diabetes are phenotypically heterogeneous. Obesity is a major risk factor for type 2 diabetes, but patients vary appreciably in body mass index. We hypothesized that the genetic predisposition to the disease may be different in lean (BMI<25 Kg/m2) compared to obese cases (BMI≥30 Kg/m2). We performed two case-control genome-wide studies using two accepted cut-offs for defining individuals as overweight or obese. We used 2,112 lean type 2 diabetes cases (BMI<25 kg/m2) or 4,123 obese cases (BMI≥30 kg/m2), and 54,412 un-stratified controls. Replication was performed in 2,881 lean cases or 8,702 obese cases, and 18,957 un-stratified controls. To assess the effects of known signals, we tested the individual and combined effects of SNPs representing 36 type 2 diabetes loci. After combining data from discovery and replication datasets, we identified two signals not previously reported in Europeans. A variant (rs8090011) in the LAMA1 gene was associated with type 2 diabetes in lean cases (P = 8.4×10−9, OR = 1.13 [95% CI 1.09–1.18]), and this association was stronger than that in obese cases (P = 0.04, OR = 1.03 [95% CI 1.00–1.06]). A variant in HMG20A—previously identified in South Asians but not Europeans—was associated with type 2 diabetes in obese cases (P = 1.3×10−8, OR = 1.11 [95% CI 1.07–1.15]), although this association was not significantly stronger than that in lean cases (P = 0.02, OR = 1.09 [95% CI 1.02–1.17]). For 36 known type 2 diabetes loci, 29 had a larger odds ratio in the lean compared to obese (binomial P = 0.0002). In the lean analysis, we observed a weighted per-risk allele OR = 1.13 [95% CI 1.10–1.17], P = 3.2×10−14. This was larger than the same model fitted in the obese analysis where the OR = 1.06 [95% CI 1.05–1.08], P = 2.2×10−16. This study provides evidence that stratification of type 2 diabetes cases by BMI may help identify additional risk variants and that lean cases may have a stronger genetic predisposition to type 2 diabetes.

Author Summary

Individuals with Type 2 diabetes (T2D) can present with variable clinical characteristics. It is well known that obesity is a major risk factor for type 2 diabetes, yet patients can vary considerably—there are many lean diabetes patients and many overweight people without diabetes. We hypothesized that the genetic predisposition to the disease may be different in lean (BMI<25 Kg/m2) compared to obese cases (BMI≥30 Kg/m2). Specifically, as lean T2D patients had lower risk than obese patients, they must have been more genetically susceptible. Using genetic data from multiple genome-wide association studies, we tested genetic markers across the genome in 2,112 lean type 2 diabetes cases (BMI<25 kg/m2), 4,123 obese cases (BMI≥30 kg/m2), and 54,412 healthy controls. We confirmed our results in an additional 2,881 lean cases, 8,702 obese cases, and 18,957 healthy controls. Using these data we found differences in genetic enrichment between lean and obese cases, supporting our original hypothesis. We also searched for genetic variants that may be risk factors only in lean or obese patients and found two novel gene regions not previously reported in European individuals. These findings may influence future study design for type 2 diabetes and provide further insight into the biology of the disease.

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−8–1.2×10−43). 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−4). 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−3, n = 22,044), increased triglycerides (p = 2.6×10−14, n = 93,440), increased waist-to-hip ratio (p = 1.8×10−5, n = 77,167), increased glucose two hours post oral glucose tolerance testing (p = 4.4×10−3, n = 15,234), increased fasting insulin (p = 0.015, n = 48,238), but with lower in HDL-cholesterol concentrations (p = 4.5×10−13, n = 96,748) and decreased BMI (p = 1.4×10−4, n = 121,335). These findings identify novel genetic determinants of adiponectin levels, which, taken together, influence risk of T2D and markers of insulin resistance.

Author Summary

Serum adiponectin levels are highly heritable and are inversely correlated with the risk of type 2 diabetes (T2D), coronary artery disease, stroke, and several metabolic traits. To identify common genetic variants associated with adiponectin levels and risk of T2D and metabolic traits, we conducted a meta-analysis of genome-wide association studies of 45,891 multi-ethnic individuals. In addition to confirming that variants at the ADIPOQ and CDH13 loci influence adiponectin levels, our analyses revealed that 10 new loci also affecting circulating adiponectin levels. We demonstrated that expression levels of several genes in these candidate regions are associated with serum adiponectin levels. Using a powerful novel method to assess the contribution of the identified variants with other traits using summary-level results from large-scale GWAS consortia, we provide evidence that the risk alleles for adiponectin are associated with deleterious changes in T2D risk and metabolic syndrome traits (triglycerides, HDL, post-prandial glucose, insulin, and waist-to-hip ratio), demonstrating that the identified loci, taken together, impact upon metabolic disease.

To identify loci for age at menarche, we performed a meta-analysis of 32 genome-wide association studies in 87,802 women of European descent, with replication in up to 14,731 women. In addition to the known loci at LIN28B (P=5.4×10−60) and 9q31.2 (P=2.2×10−33), we identified 30 novel menarche loci (all P<5×10−8) and found suggestive evidence for a further 10 loci (P<1.9×10−6). New loci included four previously associated with BMI (in/near FTO, SEC16B, TRA2B and TMEM18), three in/near other genes implicated in energy homeostasis (BSX, CRTC1, and MCHR2), and three in/near genes implicated in hormonal regulation (INHBA, PCSK2 and RXRG). Ingenuity and MAGENTA pathway analyses identified coenzyme A and fatty acid biosynthesis as biological processes related to menarche timing.

Testosterone concentrations in men are associated with cardiovascular morbidity, osteoporosis, and mortality and are affected by age, smoking, and obesity. Because of serum testosterone's high heritability, we performed a meta-analysis of genome-wide association data in 8,938 men from seven cohorts and followed up the genome-wide significant findings in one in silico (n = 871) and two de novo replication cohorts (n = 4,620) to identify genetic loci significantly associated with serum testosterone concentration in men. All these loci were also associated with low serum testosterone concentration defined as <300 ng/dl. Two single-nucleotide polymorphisms at the sex hormone-binding globulin (SHBG) locus (17p13-p12) were identified as independently associated with serum testosterone concentration (rs12150660, p = 1.2×10−41 and rs6258, p = 2.3×10−22). Subjects with ≥3 risk alleles of these variants had 6.5-fold higher risk of having low serum testosterone than subjects with no risk allele. The rs5934505 polymorphism near FAM9B on the X chromosome was also associated with testosterone concentrations (p = 5.6×10−16). The rs6258 polymorphism in exon 4 of SHBG affected SHBG's affinity for binding testosterone and the measured free testosterone fraction (p<0.01). Genetic variants in the SHBG locus and on the X chromosome are associated with a substantial variation in testosterone concentrations and increased risk of low testosterone. rs6258 is the first reported SHBG polymorphism, which affects testosterone binding to SHBG and the free testosterone fraction and could therefore influence the calculation of free testosterone using law-of-mass-action equation.

Author Summary

Testosterone is the most important testicular androgen in men. Low serum testosterone concentrations are associated with cardiovascular morbidity, metabolic syndrome, type 2 diabetes mellitus, atherosclerosis, osteoporosis, sarcopenia, and increased mortality risk. Thus, there is growing evidence that serum testosterone is a valuable biomarker of men's overall health status. Studies in male twins indicate that there is a strong heritability of serum testosterone. Here we perform a large-scale genome-wide association study to examine the effects of common genetic variants on serum testosterone concentrations. By examining 14,429 men, we show that genetic variants in the sex hormone-binding globulin (SHBG) locus and on the X chromosome are associated with a substantial variation in serum testosterone concentrations and increased risk of low testosterone. The reported associations may now be used in order to better understand the functional background of recently identified disease associations related to low testosterone. Importantly, we identified the first known genetic variant, which affects SHBG's affinity for binding testosterone and the free testosterone fraction and could therefore influence the calculation of free testosterone. This finding suggests that individual-based SHBG-testosterone affinity constants are required depending on the genotype of this single-nucleotide polymorphism.

Human height is a classic, highly heritable quantitative trait. To begin to identify genetic variants influencing height, we examined genome-wide association data from 4,921 individuals. Common variants in the HMGA2 oncogene, exemplified by rs1042725, were associated with height (P = 4 × 10−8). HMGA2 is also a strong biological candidate for height, as rare, severe mutations in this gene alter body size in mice and humans, so we tested rs1042725 in additional samples. We confirmed the association in 19,064 adults from four further studies (P = 3 × 10−11, overall P = 4 × 10−16, including the genome-wide association data). We also observed the association in children (P = 1 × 10−6, N = 6,827) and a tall/short case-control study (P = 4 × 10−6, N = 3,207). We estimate that rs1042725 explains ~0.3% of population variation in height (~0.4 cm increased adult height per C allele). There are few examples of common genetic variants reproducibly associated with human quantitative traits; these results represent, to our knowledge, the first consistently replicated association with adult and childhood height.

Dehydroepiandrosterone sulphate (DHEAS) is the most abundant circulating steroid secreted by adrenal glands—yet its function is unknown. Its serum concentration declines significantly with increasing age, which has led to speculation that a relative DHEAS deficiency may contribute to the development of common age-related diseases or diminished longevity. We conducted a meta-analysis of genome-wide association data with 14,846 individuals and identified eight independent common SNPs associated with serum DHEAS concentrations. Genes at or near the identified loci include ZKSCAN5 (rs11761528; p = 3.15×10−36), SULT2A1 (rs2637125; p = 2.61×10−19), ARPC1A (rs740160; p = 1.56×10−16), TRIM4 (rs17277546; p = 4.50×10−11), BMF (rs7181230; p = 5.44×10−11), HHEX (rs2497306; p = 4.64×10−9), BCL2L11 (rs6738028; p = 1.72×10−8), and CYP2C9 (rs2185570; p = 2.29×10−8). These genes are associated with type 2 diabetes, lymphoma, actin filament assembly, drug and xenobiotic metabolism, and zinc finger proteins. Several SNPs were associated with changes in gene expression levels, and the related genes are connected to biological pathways linking DHEAS with ageing. This study provides much needed insight into the function of DHEAS.

Author Summary

Dehydroepiandrosterone sulphate (DHEAS), mainly secreted by the adrenal gland, is the most abundant circulating steroid in humans. It shows a significant physiological decline after the age of 25 and diminishes about 95% by the age of 85 years, which has led to speculation that a relative DHEAS deficiency may contribute to the development of common age-related diseases or diminished longevity. Twin- and family-based studies have shown that there is a substantial genetic effect with heritability estimate of 60%, but no specific genes regulating serum DHEAS concentration have been identified to date. Here we take advantage of recent technical and methodological advances to examine the effects of common genetic variants on serum DHEAS concentrations. By examining 14,846 Caucasian individuals, we show that eight common genetic variants are associated with serum DHEAS concentrations. Genes at or near these genetic variants include BCL2L11, ARPC1A, ZKSCAN5, TRIM4, HHEX, CYP2C9, BMF, and SULT2A1. These genes have various associations with steroid hormone metabolism—co-morbidities of ageing including type 2 diabetes, lymphoma, actin filament assembly, drug and xenobiotic metabolism, and zinc finger proteins—suggesting a wider functional role for DHEAS than previously thought.

Objective

Plasma adiponectin is strongly associated with various components of metabolic syndrome, type 2 diabetes and cardiovascular outcomes. Concentrations are highly heritable and differ between men and women. We therefore aimed to investigate the genetics of plasma adiponectin in men and women.

Methods

We combined genome-wide association scans of three population-based studies including 4659 persons. For the replication stage in 13795 subjects, we selected the 20 top signals of the combined analysis, as well as the 10 top signals with p-values less than 1.0*10-4 for each the men- and the women-specific analyses. We further selected 73 SNPs that were consistently associated with metabolic syndrome parameters in previous genome-wide association studies to check for their association with plasma adiponectin.

Results

The ADIPOQ locus showed genome-wide significant p-values in the combined (p=4.3*10-24) as well as in both women- and men-specific analyses (p=8.7*10-17 and p=2.5*10-11, respectively). None of the other 39 top signal SNPs showed evidence for association in the replication analysis. None of 73 SNPs from metabolic syndrome loci exhibited association with plasma adiponectin (p>0.01).

Conclusions

We demonstrated the ADIPOQ gene as the only major gene for plasma adiponectin, which explains 6.7% of the phenotypic variance. We further found that neither this gene nor any of the metabolic syndrome loci explained the sex differences observed for plasma adiponectin. Larger studies are needed to identify more moderate genetic determinants of plasma adiponectin.

Women become infertile approximately 10 years before menopause, and as more women delay childbirth into their 30s, the number of women who experience infertility is likely to increase. Tests that predict the timing of menopause would allow women to make informed reproductive decisions. Current predictors are only effective just prior to menopause, and there are no long-range indicators. Age at menopause and early menopause (EM) are highly heritable, suggesting a genetic aetiology. Recent genome-wide scans have identified four loci associated with variation in the age of normal menopause (40–60 years). We aimed to determine whether theses loci are also risk factors for EM. We tested the four menopause-associated genetic variants in a cohort of approximately 2000 women with menopause ≤45 years from the Breakthrough Generations Study (BGS). All four variants significantly increased the odds of having EM. Comparing the 4.5% of individuals with the lowest number of risk alleles (two or three) with the 3.0% with the highest number (eight risk alleles), the odds ratio was 4.1 (95% CI 2.4–7.1, P = 4.0 × 10−7). In combination, the four variants discriminated EM cases with a receiver operator characteristic area under the curve of 0.6. Four common genetic variants identified by genome-wide association studies, had a significant impact on the odds of having EM in an independent cohort from the BGS. The discriminative power is still limited, but as more variants are discovered they may be useful for predicting reproductive lifespan.

Aims and Hypothesis

Circulating β-carotene levels are inversely associated with type 2 diabetes risk, but the causal direction of this association is not certain. In this study we used a Mendelian Randomization approach to provide evidence for or against the causal role of the anti-oxidant vitamin β-carotene in type 2 diabetes.

Methods

We used a common polymorphism (rs6564851) near the β-carotene 15,15'-Monooxygenase 1 (BCMO1) gene that is strongly associated with circulating β-carotene levels (P = 2×10−24) - each G allele is associated with a 0.27 standard deviation increase in levels. We used data from the InCHIANTI study and the ULSAM study to estimate the association between β-carotene levels and type 2 diabetes. We next used a triangulation approach to estimate the expected effect of rs6564851 on type 2 diabetes risk, and compared this to the observed effect using data from 4549 type 2 diabetes cases and 5579 controls from the DIAGRAM consortium.

Results

A 0.27 standard deviation increase in β-carotene levels is associated with an odds ratio of 0.90 (0.86–0.95) for type 2 diabetes in the InCHIANTI study. This association is similar to that of the ULSAM study, OR (0.90 (0.84–0.97)). In contrast there was no association between rs6564851 and type 2 diabetes (OR 0.98 (0.93–1.04, P = 0.58), and this effect size was smaller than that expected given the known associations between rs6564851 and β-carotene levels and the associations between β-carotene levels and type 2 diabetes.

Conclusion

Our Mendelian Randomization studies are in keeping with randomized controlled trials that suggest β-carotene is not causally protective against type 2 diabetes.

We conducted a meta-analysis of genome-wide association data to detect genes influencing age at menarche in 17,510 women. The strongest signal was at 9q31.2 (P = 1.7 × 10−9), where the nearest genes include TMEM38B, FKTN, FSD1L, TAL2 and ZNF462. The next best signal was near the LIN28B gene (rs7759938; P = 7.0 × 10−9), which also influences adult height. We provide the first evidence for common genetic variants influencing female sexual maturation.

OBJECTIVE

Low birth weight is associated with an increased risk of type 2 diabetes. The mechanisms underlying this association are unknown and may represent intrauterine programming or two phenotypes of one genotype. The fetal insulin hypothesis proposes that common genetic variants that reduce insulin secretion or action may predispose to type 2 diabetes and also reduce birth weight, since insulin is a key fetal growth factor. We tested whether common genetic variants that predispose to type 2 diabetes also reduce birth weight.

RESEARCH DESIGN AND METHODS

We genotyped single-nucleotide polymorphisms (SNPs) at five recently identified type 2 diabetes loci (CDKAL1, CDKN2A/B, HHEX-IDE, IGF2BP2, and SLC30A8) in 7,986 mothers and 19,200 offspring from four studies of white Europeans. We tested the association between maternal or fetal genotype at each locus and birth weight of the offspring.

RESULTS

We found that type 2 diabetes risk alleles at the CDKAL1 and HHEX-IDE loci were associated with reduced birth weight when inherited by the fetus (21 g [95% CI 11–31], P = 2 × 10−5, and 14 g [4–23], P = 0.004, lower birth weight per risk allele, respectively). The 4% of offspring carrying four risk alleles at these two loci were 80 g (95% CI 39–120) lighter at birth than the 8% carrying none (Ptrend = 5 × 10−7). There were no associations between birth weight and fetal genotypes at the three other loci or maternal genotypes at any locus.

CONCLUSIONS

Our results are in keeping with the fetal insulin hypothesis and provide robust evidence that common disease-associated variants can alter size at birth directly through the fetal genotype.

OBJECTIVE

Recent genome-wide association studies have resulted in a dramatic increase in our knowledge of the genetic loci involved in type 2 diabetes. In a complementary approach to these single-marker studies, we attempted to identify biological pathways associated with type 2 diabetes. This approach could allow us to identify additional risk loci.

RESEARCH DESIGN AND METHODS

We used individual level genotype data generated from the Wellcome Trust Case Control Consortium (WTCCC) type 2 diabetes study, consisting of 393,143 autosomal SNPs, genotyped across 1,924 case subjects and 2,938 control subjects. We sought additional evidence from summary level data available from the Diabetes Genetics Initiative (DGI) and the Finland-United States Investigation of NIDDM Genetics (FUSION) studies. Statistical analysis of pathways was performed using a modification of the Gene Set Enrichment Algorithm (GSEA). A total of 439 pathways were analyzed from the Kyoto Encyclopedia of Genes and Genomes, Gene Ontology, and BioCarta databases.

RESULTS

After correcting for the number of pathways tested, we found no strong evidence for any pathway showing association with type 2 diabetes (top Padj = 0.31). The candidate WNT-signaling pathway ranked top (nominal P = 0.0007, excluding TCF7L2; P = 0.002), containing a number of promising single gene associations. These include CCND2 (rs11833537; P = 0.003), SMAD3 (rs7178347; P = 0.0006), and PRICKLE1 (rs1796390; P = 0.001), all expressed in the pancreas.

CONCLUSIONS

Common variants involved in type 2 diabetes risk are likely to occur in or near genes in multiple pathways. Pathway-based approaches to genome-wide association data may be more successful for some complex traits than others, depending on the nature of the underlying disease physiology.

The adipocyte-derived protein adiponectin is highly heritable and inversely associated with risk of type 2 diabetes mellitus (T2D) and coronary heart disease (CHD). We meta-analyzed 3 genome-wide association studies for circulating adiponectin levels (n = 8,531) and sought validation of the lead single nucleotide polymorphisms (SNPs) in 5 additional cohorts (n = 6,202). Five SNPs were genome-wide significant in their relationship with adiponectin (P≤5×10−8). We then tested whether these 5 SNPs were associated with risk of T2D and CHD using a Bonferroni-corrected threshold of P≤0.011 to declare statistical significance for these disease associations. SNPs at the adiponectin-encoding ADIPOQ locus demonstrated the strongest associations with adiponectin levels (P-combined = 9.2×10−19 for lead SNP, rs266717, n = 14,733). A novel variant in the ARL15 (ADP-ribosylation factor-like 15) gene was associated with lower circulating levels of adiponectin (rs4311394-G, P-combined = 2.9×10−8, n = 14,733). This same risk allele at ARL15 was also associated with a higher risk of CHD (odds ratio [OR] = 1.12, P = 8.5×10−6, n = 22,421) more nominally, an increased risk of T2D (OR = 1.11, P = 3.2×10−3, n = 10,128), and several metabolic traits. Expression studies in humans indicated that ARL15 is well-expressed in skeletal muscle. These findings identify a novel protein, ARL15, which influences circulating adiponectin levels and may impact upon CHD risk.

Author Summary

Through a meta-analysis of genome-wide association studies of 14,733 individuals, we identified common base-pair variants in the genome which influence circulating adiponectin levels. Since adiponectin is an adipocyte-derived circulating protein which has been inversely associated with risk of obesity-related diseases such as type 2 diabetes (T2D) and coronary heart disease (CHD), we next sought to understand if the identified variants influencing adiponectin levels also influence risk of T2D, CHD, and several metabolic traits. In addition to confirming that variation at the ADIPOQ locus influences adiponectin levels, our analyses point to a variant in the ARL15 (ADP-ribosylation factor-like 15) locus which decreases adiponectin levels and increases risk of CHD and T2D. Further, this same variant was associated with increased fasting insulin levels and glycated hemoglobin. While the function of ARL15 is not known, we provide insight into the tissue specificity of ARL15 expression. These results thus provide novel insights into the physiology of the adiponectin pathway and obesity-related diseases.

Objective

Low birth weight is associated with an increased risk of type 2 diabetes. The mechanisms underlying this association are unknown and may represent intrauterine programming or two phenotypes of one genotype. The fetal insulin hypothesis proposes that common genetic variants that reduce insulin secretion or action may predispose to type 2 diabetes and also reduce birth weight, since insulin is a key fetal growth factor. We tested whether common genetic variants that predispose to type 2 diabetes also reduce birth weight.

Research design and methods

We genotyped single nucleotide polymorphisms (SNPs) at five recently identified type 2 diabetes loci (CDKAL1, CDKN2A/B, HHEX-IDE, IGF2BP2 and SLC30A8) in 7986 mothers and 19200 offspring from four studies of white Europeans. We tested the association between maternal or fetal genotype at each locus and birth weight of the offspring.

Results

We found that type 2 diabetes risk alleles at the CDKAL1 and HHEX-IDE loci were associated with reduced birth weight when inherited by the fetus: 21g [95%CI:11-31g], P=2×10-5 and 14g [4-23g], P=0.004 lower birth weight per risk allele, respectively. The 4% of offspring carrying four risk alleles at these two loci were 80g [39-120g] lighter at birth than the 8% carrying none (Ptrend =5×10-7). There were no associations between birth weight and fetal genotypes at the three other loci, or maternal genotypes at any locus.

Conclusions

Our results are in keeping with the fetal insulin hypothesis and provide robust evidence that common disease-associated variants can alter size at birth directly through the fetal genotype.

Epidemiological studies consistently show that circulating sex hormone binding globulin (SHBG) levels are lower in type 2 diabetes patients than non-diabetic individuals, but the causal nature of this association is controversial. Genetic studies can help dissect causal directions of epidemiological associations because genotypes are much less likely to be confounded, biased or influenced by disease processes. Using this Mendelian randomization principle, we selected a common single nucleotide polymorphism (SNP) near the SHBG gene, rs1799941, that is strongly associated with SHBG levels. We used data from this SNP, or closely correlated SNPs, in 27 657 type 2 diabetes patients and 58 481 controls from 15 studies. We then used data from additional studies to estimate the difference in SHBG levels between type 2 diabetes patients and controls. The SHBG SNP rs1799941 was associated with type 2 diabetes [odds ratio (OR) 0.94, 95% CI: 0.91, 0.97; P = 2 × 10−5], with the SHBG raising allele associated with reduced risk of type 2 diabetes. This effect was very similar to that expected (OR 0.92, 95% CI: 0.88, 0.96), given the SHBG-SNP versus SHBG levels association (SHBG levels are 0.2 standard deviations higher per copy of the A allele) and the SHBG levels versus type 2 diabetes association (SHBG levels are 0.23 standard deviations lower in type 2 diabetic patients compared to controls). Results were very similar in men and women. There was no evidence that this variant is associated with diabetes-related intermediate traits, including several measures of insulin secretion and resistance. Our results, together with those from another recent genetic study, strengthen evidence that SHBG and sex hormones are involved in the aetiology of type 2 diabetes.

Adult height is a model polygenic trait, but there has been limited success in identifying the genes underlying its normal variation. To identify genetic variants influencing adult human height, we used genome-wide association data from 13,665 individuals and genotyped 39 variants in an additional 16,482 samples. We identified 20 variants associated with adult height (P < 5 × 10−7, with 10 reaching P < 1 × 10−10). Combined, the 20 SNPs explain ~3% of height variation, with a ~5 cm difference between the 6.2% of people with 17 or fewer ‘tall’ alleles compared to the 5.5% with 27 or more ‘tall’ alleles. The loci we identified implicate genes in Hedgehog signaling (IHH, HHIP, PTCH1), extracellular matrix (EFEMP1, ADAMTSL3, ACAN) and cancer (CDK6, HMGA2, DLEU7) pathways, and provide new insights into human growth and developmental processes. Finally, our results provide insights into the genetic architecture of a classic quantitative trait.

Obesity is a serious international health problem that increases the risk of several common diseases. The genetic factors predisposing to obesity are poorly understood. A genome-wide search for type 2 diabetes–susceptibility genes identified a common variant in the FTO (fat mass and obesity associated) gene that predisposes to diabetes through an effect on body mass index (BMI). An additive association of the variant with BMI was replicated in 13 cohorts with 38,759 participants. The 16% of adults who are homozygous for the risk allele weighed about 3 kilograms more and had 1.67-fold increased odds of obesity when compared with those not inheriting a risk allele. This association was observed from age 7 years upward and reflects a specific increase in fat mass.

There is considerable evidence that human genetic variation influences gene expression. Genome-wide studies have revealed that mRNA levels are associated with genetic variation in or close to the gene coding for those mRNA transcripts – cis effects, and elsewhere in the genome – trans effects. The role of genetic variation in determining protein levels has not been systematically assessed. Using a genome-wide association approach we show that common genetic variation influences levels of clinically relevant proteins in human serum and plasma. We evaluated the role of 496,032 polymorphisms on levels of 42 proteins measured in 1200 fasting individuals from the population based InCHIANTI study. Proteins included insulin, several interleukins, adipokines, chemokines, and liver function markers that are implicated in many common diseases including metabolic, inflammatory, and infectious conditions. We identified eight Cis effects, including variants in or near the IL6R (p = 1.8×10−57), CCL4L1 (p = 3.9×10−21), IL18 (p = 6.8×10−13), LPA (p = 4.4×10−10), GGT1 (p = 1.5×10−7), SHBG (p = 3.1×10−7), CRP (p = 6.4×10−6) and IL1RN (p = 7.3×10−6) genes, all associated with their respective protein products with effect sizes ranging from 0.19 to 0.69 standard deviations per allele. Mechanisms implicated include altered rates of cleavage of bound to unbound soluble receptor (IL6R), altered secretion rates of different sized proteins (LPA), variation in gene copy number (CCL4L1) and altered transcription (GGT1). We identified one novel trans effect that was an association between ABO blood group and tumour necrosis factor alpha (TNF-alpha) levels (p = 6.8×10−40), but this finding was not present when TNF-alpha was measured using a different assay , or in a second study, suggesting an assay-specific association. Our results show that protein levels share some of the features of the genetics of gene expression. These include the presence of strong genetic effects in cis locations. The identification of protein quantitative trait loci (pQTLs) may be a powerful complementary method of improving our understanding of disease pathways.

Author Summary

One of the central dogmas of molecular genetics is that DNA is transcribed to RNA which is translated to protein and alterations to proteins can influence human diseases. Genome-wide association studies have recently revealed many new DNA variants that influence human diseases. To complement these efforts, several genome-wide studies have established that DNA variation influences mRNA expression levels. Loci influencing mRNA levels have been termed “eQTLs”. In this study we have performed the first genome-wide association study of the third piece in this jigsaw – the role of DNA variation in relation to protein levels, or “pQTLs”. We analysed 42 proteins measured in blood fractions from the InCHIANTI study. We identified eight cis effects including common variants in or near the IL6R, CCL4, IL18, LPA, GGT1, SHBG, CRP and IL1RN genes, all associated with blood levels of their respective protein products. Mechanisms implicated included altered transcription (GGT1) but also rates of cleavage of bound to unbound soluble receptor (IL6R), altered secretion rates of different sized proteins (LPA) and variation in gene copy number (CCL4). Blood levels of many of these proteins are correlated with human diseases and the identification of “pQTLs” may in turn help our understanding of disease.

Sarcoidosis is a systemic inflammatory disease characterized by the formation of granulomas in affected organs. Genome-wide association studies (GWASs) of this disease have been conducted only in European population. We present the first sarcoidosis GWAS in African Americans (AAs, 818 cases and 1,088 related controls) followed by replication in independent sets of AAs (455 cases and 557 controls) and European Americans (EAs, 442 cases and 2,284 controls). We evaluated >6 million SNPs either genotyped using the Illumina Omni1-Quad array or imputed from the 1000 Genomes Project data. We identified a novel sarcoidosis-associated locus, NOTCH4, that reached genome-wide significance in the combined AA samples (rs715299, PAA-meta = 6.51×10−10) and demonstrated the independence of this locus from others in the MHC region in the same sample. We replicated previous European GWAS associations within HLA-DRA, HLA-DRB5, HLA-DRB1, BTNL2, and ANXA11 in both our AA and EA datasets. We also confirmed significant associations to the previously reported HLA-C and HLA-B regions in the EA but not AA samples. We further identified suggestive associations with several other genes previously reported in lung or inflammatory diseases.

Background

Recent genome-wide association studies (GWAS) have identified common variants in the UMOD region associated with kidney function and disease in the general population. To identify novel rare variants as well as common variants that may account for this GWAS signal, the exons and 4 kb upstream region of UMOD were sequenced.

Methodology/Principal Findings

Individuals (n = 485) were selected based on presence of the GWAS risk haplotype and chronic kidney disease (CKD) in the ARIC Study and on the extremes of of the UMOD gene product, uromodulin, in urine (Tamm Horsfall protein, THP) in the Framingham Heart Study (FHS). Targeted sequencing was conducted using capillary based Sanger sequencing (3730 DNA Analyzer). Variants were tested for association with THP concentrations and estimated glomerular filtration rate (eGFR), and identified non-synonymous coding variants were genotyped in up to 22,546 follow-up samples. Twenty-four and 63 variants were identified in the 285 ARIC and 200 FHS participants, respectively. In both studies combined, there were 33 common and 54 rare (MAF<0.05) variants. Five non-synonymous rare variants were identified in FHS; borderline enrichment of rare variants was found in the extremes of THP (SKAT p-value = 0.08). Only V458L was associated with THP in the FHS general-population validation sample (p = 9*10−3, n = 2,522), but did not show direction-consistent and significant association with eGFR in both the ARIC (n = 14,635) and FHS (n = 7,520) validation samples. Pooling all non-synonymous rare variants except V458L together showed non-significant associations with THP and eGFR in the FHS validation sample. Functional studies of V458L revealed no alternations in protein trafficking.

Conclusions/Significance

Multiple novel rare variants in the UMOD region were identified, but none were consistently associated with eGFR in two independent study samples. Only V458L had modest association with THP levels in the general population and thus could not account for the observed GWAS signal.

Objective

Low serum 25(OH)D levels are associated with cardiovascular risk factors, and also predict future myocardial infarction (MI), type 2 diabetes (T2DM), cancer and all-cause mortality. Recently several single nucleotide polymorphisms (SNPs) associated with serum 25-hydroxyvitamin D (25(OH)D) level have been identified. If these relations are causal one would expect a similar association between these SNPs and health.

Methods

DNA was prepared from subjects who participated in the fourth survey of the Tromsø Study in 1994–1995 and who were registered with the endpoints MI, T2DM, cancer or death as well as a randomly selected control group. The endpoint registers were complete up to 2007–2010. Genotyping was performed for 17 SNPs related to the serum 25(OH)D level.

Results

A total of 9528 subjects were selected for genetic analyses which were successfully performed for at least one SNP in 9471 subjects. Among these, 2025 were registered with MI, 1092 with T2DM, 2924 with cancer and 3828 had died. The mean differences in serum 25(OH)D levels between SNP genotypes with the lowest and highest serum 25(OH)D levels varied from 0.1 to 7.8 nmol/L. A genotype score based on weighted risk alleles regarding low serum 25(OH)D levels was established. There was no consistent association between the genotype score or individuals SNPs and MI, T2DM, cancer, mortality or risk factors for disease. However, for rs6013897 genotypes (located at the 24-hydroxylase gene (CYP24A1)) there was a significant association with breast cancer (P<0.05).

Conclusion

Our results do not support nor exclude a causal relationship between serum 25(OH)D levels and MI, T2DM, cancer or mortality, and our observation on breast cancer needs confirmation. Further genetic studies are warranted, particularly in populations with vitamin D deficiency.

Trial Registration

ClinicalTrials.gov NCT01395303

The Metabochip is a custom genotyping array designed for replication and fine mapping of metabolic, cardiovascular, and anthropometric trait loci and includes low frequency variation content identified from the 1000 Genomes Project. It has 196,725 SNPs concentrated in 257 genomic regions. We evaluated the Metabochip in 5,863 African Americans; 89% of all SNPs passed rigorous quality control with a call rate of 99.9%. Two examples illustrate the value of fine mapping with the Metabochip in African-ancestry populations. At CELSR2/PSRC1/SORT1, we found the strongest associated SNP for LDL-C to be rs12740374 (p = 3.5×10−11), a SNP indistinguishable from multiple SNPs in European ancestry samples due to high correlation. Its distinct signal supports functional studies elsewhere suggesting a causal role in LDL-C. At CETP we found rs17231520, with risk allele frequency 0.07 in African Americans, to be associated with HDL-C (p = 7.2×10−36). This variant is very rare in Europeans and not tagged in common GWAS arrays, but was identified as associated with HDL-C in African Americans in a single-gene study. Our results, one narrowing the risk interval and the other revealing an associated variant not found in Europeans, demonstrate the advantages of high-density genotyping of common and rare variation for fine mapping of trait loci in African American samples.