We aimed to investigate whether maternal and fetal genetic predispositions to insulin deficiency and resistance affect offspring fetal growth through distinct pathways in multi-ethnic populations. In Show more
We aimed to investigate whether maternal and fetal genetic predispositions to insulin deficiency and resistance affect offspring fetal growth through distinct pathways in multi-ethnic populations. In 5065 multi-ethnic mother-infant pairs, we examined the conditional associations of maternal and fetal partitioned polygenic risk scores (pPRSs) for type 2 diabetes-related pathways with fetal growth outcomes, including birthweight, sum of skinfold thicknesses (SSF), large-for-gestational-age (LGA) births and small-for-gestational-age (SGA) births. Two-sample Mendelian randomisation (2SMR) in Europeans was performed for triangulation. Exposures were eight type 2 diabetes-related pathways (n=1,812,017), eight beta cell function indices (n=26,356) and two insulin sensitivity indices (n=53,657). Outcomes were maternal and fetal genetically determined birthweight (n=406,063). Mediation analysis was used to assess the mediation effects of maternal glucose levels and BMI on maternal genetic effects and of cord blood C-peptide on fetal genetic effects. Co-localisation analyses were performed to test for shared causal variants. Fetal type 2 diabetes polygenic risk score (PRS) and pPRSs for lipodystrophy-related insulin resistance and impaired fasting glucose (IFG)-related insulin deficiency were associated with lower birthweight and SSF, while maternal type 2 diabetes PRS and pPRSs for IFG-related insulin deficiency and obesity-related insulin resistance were associated with higher offspring birthweight, SSF and LGA. These associations were consistent across five ethnic groups. Maternal post-load hyperglycaemia mediated 44.2% and 34.2% of the effects of type 2 diabetes PRS and IFG pPRS, respectively, while maternal BMI mediated 43.4% of the effect of Obesity pPRS. 2SMR found consistent results in Europeans and further revealed that fetal insulin sensitivity index and corrected insulin response were associated with higher birthweight. Some loci with shared causal variants acted through multiple pathways, including CDKAL1, TCF7L2, ADCY5 and MACF1. Reduced fetal growth may be driven by lipodystrophy-related insulin resistance and IFG-related insulin deficiency pathways. Targeting pregnant women with high type 2 diabetes PRS/pPRS and prescribing interventions to reduce their post-load hyperglycaemia and BMI may help reduce offspring risk of LGA. Show less
Daytime napping is a common, heritable behavior, but its genetic basis and causal relationship with cardiometabolic health remain unclear. Here, we perform a genome-wide association study of self-repo Show more
Daytime napping is a common, heritable behavior, but its genetic basis and causal relationship with cardiometabolic health remain unclear. Here, we perform a genome-wide association study of self-reported daytime napping in the UK Biobank (n = 452,633) and identify 123 loci of which 61 replicate in the 23andMe research cohort (n = 541,333). Findings include missense variants in established drug targets for sleep disorders (HCRTR1, HCRTR2), genes with roles in arousal (TRPC6, PNOC), and genes suggesting an obesity-hypersomnolence pathway (PNOC, PATJ). Association signals are concordant with accelerometer-measured daytime inactivity duration and 33 loci colocalize with loci for other sleep phenotypes. Cluster analysis identifies three distinct clusters of nap-promoting mechanisms with heterogeneous associations with cardiometabolic outcomes. Mendelian randomization shows potential causal links between more frequent daytime napping and higher blood pressure and waist circumference. Show less
Lean body mass (LM) plays an important role in mobility and metabolic function. We previously identified five loci associated with LM adjusted for fat mass in kilograms. Such an adjustment may reduce Show more
Lean body mass (LM) plays an important role in mobility and metabolic function. We previously identified five loci associated with LM adjusted for fat mass in kilograms. Such an adjustment may reduce the power to identify genetic signals having an association with both lean mass and fat mass. To determine the impact of different fat mass adjustments on genetic architecture of LM and identify additional LM loci. We performed genome-wide association analyses for whole-body LM (20 cohorts of European ancestry with n = 38,292) measured using dual-energy X-ray absorptiometry) or bioelectrical impedance analysis, adjusted for sex, age, age2, and height with or without fat mass adjustments (Model 1 no fat adjustment; Model 2 adjustment for fat mass as a percentage of body mass; Model 3 adjustment for fat mass in kilograms). Seven single-nucleotide polymorphisms (SNPs) in separate loci, including one novel LM locus (TNRC6B), were successfully replicated in an additional 47,227 individuals from 29 cohorts. Based on the strengths of the associations in Model 1 vs Model 3, we divided the LM loci into those with an effect on both lean mass and fat mass in the same direction and refer to those as "sumo wrestler" loci (FTO and MC4R). In contrast, loci with an impact specifically on LM were termed "body builder" loci (VCAN and ADAMTSL3). Using existing available genome-wide association study databases, LM increasing alleles of SNPs in sumo wrestler loci were associated with an adverse metabolic profile, whereas LM increasing alleles of SNPs in "body builder" loci were associated with metabolic protection. In conclusion, we identified one novel LM locus (TNRC6B). Our results suggest that a genetically determined increase in lean mass might exert either harmful or protective effects on metabolic traits, depending on its relation to fat mass. Show less
Most genetic variants identified in genome-wide association studies (GWASs) of complex traits are thought to act by affecting gene regulation rather than directly altering the protein product. As a co Show more
Most genetic variants identified in genome-wide association studies (GWASs) of complex traits are thought to act by affecting gene regulation rather than directly altering the protein product. As a consequence, the actual genes involved in disease are not necessarily the most proximal to the associated variants. By integrating data from GWAS analyses with those from genetic studies of regulatory variation, it is possible to identify variants pleiotropically associated with both a complex trait and measures of gene regulation. In this study, we used summary-data-based Mendelian randomization (SMR), a method developed to identify variants pleiotropically associated with both complex traits and gene expression, to identify variants associated with complex traits and DNA methylation. We used large DNA methylation quantitative trait locus (mQTL) datasets generated from two different tissues (blood and fetal brain) to prioritize genes for >40 complex traits with robust GWAS data and found considerable overlap with the results of SMR analyses performed with expression QTL (eQTL) data. We identified multiple examples of variable DNA methylation associated with GWAS variants for a range of complex traits, demonstrating the utility of this approach for refining genetic association signals. Show less
Circulating metabolites associated with insulin sensitivity may represent useful biomarkers, but their causal role in insulin sensitivity and diabetes is less certain. We previously identified novel m Show more
Circulating metabolites associated with insulin sensitivity may represent useful biomarkers, but their causal role in insulin sensitivity and diabetes is less certain. We previously identified novel metabolites correlated with insulin sensitivity measured by the hyperinsulinemic-euglycemic clamp. The top-ranking metabolites were in the glutathione and glycine biosynthesis pathways. We aimed to identify common genetic variants associated with metabolites in these pathways and test their role in insulin sensitivity and type 2 diabetes. With 1,004 nondiabetic individuals from the RISC study, we performed a genome-wide association study (GWAS) of 14 insulin sensitivity-related metabolites and one metabolite ratio. We replicated our results in the Botnia study (n = 342). We assessed the association of these variants with diabetes-related traits in GWAS meta-analyses (GENESIS [including RISC, EUGENE2, and Stanford], MAGIC, and DIAGRAM). We identified four associations with three metabolites-glycine (rs715 at CPS1), serine (rs478093 at PHGDH), and betaine (rs499368 at SLC6A12; rs17823642 at BHMT)-and one association signal with glycine-to-serine ratio (rs1107366 at ALDH1L1). There was no robust evidence for association between these variants and insulin resistance or diabetes. Genetic variants associated with genes in the glycine biosynthesis pathways do not provide consistent evidence for a role of glycine in diabetes-related traits. Show less
Proinsulin is a precursor of mature insulin and C-peptide. Higher circulating proinsulin levels are associated with impaired β-cell function, raised glucose levels, insulin resistance, and type 2 diab Show more
Proinsulin is a precursor of mature insulin and C-peptide. Higher circulating proinsulin levels are associated with impaired β-cell function, raised glucose levels, insulin resistance, and type 2 diabetes (T2D). Studies of the insulin processing pathway could provide new insights about T2D pathophysiology. We have conducted a meta-analysis of genome-wide association tests of ∼2.5 million genotyped or imputed single nucleotide polymorphisms (SNPs) and fasting proinsulin levels in 10,701 nondiabetic adults of European ancestry, with follow-up of 23 loci in up to 16,378 individuals, using additive genetic models adjusted for age, sex, fasting insulin, and study-specific covariates. Nine SNPs at eight loci were associated with proinsulin levels (P < 5 × 10(-8)). Two loci (LARP6 and SGSM2) have not been previously related to metabolic traits, one (MADD) has been associated with fasting glucose, one (PCSK1) has been implicated in obesity, and four (TCF7L2, SLC30A8, VPS13C/C2CD4A/B, and ARAP1, formerly CENTD2) increase T2D risk. The proinsulin-raising allele of ARAP1 was associated with a lower fasting glucose (P = 1.7 × 10(-4)), improved β-cell function (P = 1.1 × 10(-5)), and lower risk of T2D (odds ratio 0.88; P = 7.8 × 10(-6)). Notably, PCSK1 encodes the protein prohormone convertase 1/3, the first enzyme in the insulin processing pathway. A genotype score composed of the nine proinsulin-raising alleles was not associated with coronary disease in two large case-control datasets. We have identified nine genetic variants associated with fasting proinsulin. Our findings illuminate the biology underlying glucose homeostasis and T2D development in humans and argue against a direct role of proinsulin in coronary artery disease pathogenesis. Show less
Obesity is globally prevalent and highly heritable, but its underlying genetic factors remain largely elusive. To identify genetic loci for obesity susceptibility, we examined associations between bod Show more
Obesity is globally prevalent and highly heritable, but its underlying genetic factors remain largely elusive. To identify genetic loci for obesity susceptibility, we examined associations between body mass index and ∼ 2.8 million SNPs in up to 123,865 individuals with targeted follow up of 42 SNPs in up to 125,931 additional individuals. We confirmed 14 known obesity susceptibility loci and identified 18 new loci associated with body mass index (P < 5 × 10⁻⁸), one of which includes a copy number variant near GPRC5B. Some loci (at MC4R, POMC, SH2B1 and BDNF) map near key hypothalamic regulators of energy balance, and one of these loci is near GIPR, an incretin receptor. Furthermore, genes in other newly associated loci may provide new insights into human body weight regulation. Show less
Glucose levels 2 h after an oral glucose challenge are a clinical measure of glucose tolerance used in the diagnosis of type 2 diabetes. We report a meta-analysis of nine genome-wide association studi Show more
Glucose levels 2 h after an oral glucose challenge are a clinical measure of glucose tolerance used in the diagnosis of type 2 diabetes. We report a meta-analysis of nine genome-wide association studies (n = 15,234 nondiabetic individuals) and a follow-up of 29 independent loci (n = 6,958-30,620). We identify variants at the GIPR locus associated with 2-h glucose level (rs10423928, beta (s.e.m.) = 0.09 (0.01) mmol/l per A allele, P = 2.0 x 10(-15)). The GIPR A-allele carriers also showed decreased insulin secretion (n = 22,492; insulinogenic index, P = 1.0 x 10(-17); ratio of insulin to glucose area under the curve, P = 1.3 x 10(-16)) and diminished incretin effect (n = 804; P = 4.3 x 10(-4)). We also identified variants at ADCY5 (rs2877716, P = 4.2 x 10(-16)), VPS13C (rs17271305, P = 4.1 x 10(-8)), GCKR (rs1260326, P = 7.1 x 10(-11)) and TCF7L2 (rs7903146, P = 4.2 x 10(-10)) associated with 2-h glucose. Of the three newly implicated loci (GIPR, ADCY5 and VPS13C), only ADCY5 was found to be associated with type 2 diabetes in collaborating studies (n = 35,869 cases, 89,798 controls, OR = 1.12, 95% CI 1.09-1.15, P = 4.8 x 10(-18)). Show less
Levels of circulating glucose are tightly regulated. To identify new loci influencing glycemic traits, we performed meta-analyses of 21 genome-wide association studies informative for fasting glucose, Show more
Levels of circulating glucose are tightly regulated. To identify new loci influencing glycemic traits, we performed meta-analyses of 21 genome-wide association studies informative for fasting glucose, fasting insulin and indices of beta-cell function (HOMA-B) and insulin resistance (HOMA-IR) in up to 46,186 nondiabetic participants. Follow-up of 25 loci in up to 76,558 additional subjects identified 16 loci associated with fasting glucose and HOMA-B and two loci associated with fasting insulin and HOMA-IR. These include nine loci newly associated with fasting glucose (in or near ADCY5, MADD, ADRA2A, CRY2, FADS1, GLIS3, SLC2A2, PROX1 and C2CD4B) and one influencing fasting insulin and HOMA-IR (near IGF1). We also demonstrated association of ADCY5, PROX1, GCK, GCKR and DGKB-TMEM195 with type 2 diabetes. Within these loci, likely biological candidate genes influence signal transduction, cell proliferation, development, glucose-sensing and circadian regulation. Our results demonstrate that genetic studies of glycemic traits can identify type 2 diabetes risk loci, as well as loci containing gene variants that are associated with a modest elevation in glucose levels but are not associated with overt diabetes. Show less
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 Show more
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⁻⁶⁰) and 9q31.2 (P = 2.2 × 10⁻³³), we identified 30 new menarche loci (all P < 5 × 10⁻⁸) and found suggestive evidence for a further 10 loci (P < 1.9 × 10⁻⁶). The new loci included four previously associated with body mass index (in or near FTO, SEC16B, TRA2B and TMEM18), three in or near other genes implicated in energy homeostasis (BSX, CRTC1 and MCHR2) and three in or near genes implicated in hormonal regulation (INHBA, PCSK2 and RXRG). Ingenuity and gene-set enrichment pathway analyses identified coenzyme A and fatty acid biosynthesis as biological processes related to menarche timing. Show less
Common genetic variants influence plasma triacylglycerol, HDL-cholesterol (HDL-C) and glucose levels in cross-sectional studies. However, the longitudinal effects of these established variants have no Show more
Common genetic variants influence plasma triacylglycerol, HDL-cholesterol (HDL-C) and glucose levels in cross-sectional studies. However, the longitudinal effects of these established variants have not been studied. Our aim was to examine the longitudinal associations of four such variants in the apolipoprotein A-V (APOA5), lipoprotein lipase (LPL), and glucokinase (GCK) genes with fasting glucose or lipid levels. The individuals analysed were participants in the Busselton Health Survey (n = 4,554). Cross-sectional analyses of family data used the total association test. Longitudinal association analyses of unrelated participant data (n = 2,864) used linear mixed-effects models. The findings of cross-sectional association analyses replicated those of previous studies. We observed associations of the G and C alleles at the APOA5 single nucleotide polymorphisms (SNPs) rs662799 and rs3135506 with raised triacylglycerol levels (p = 0.0003 and p < 0.0001, respectively), the 447X allele at the LPL SNP rs328 with reduced triacylglycerol levels (p = 0.0004) and raised HDL-C levels (p = 0.0004), and the A allele of the GCK SNP rs1799884 with raised fasting glucose level (p = 0.015). Longitudinal association analyses showed that most of these associations did not change in the same individuals over an average follow-up time of 17.4 years, though there was some evidence that the association of the 447X allele of rs328 with raised HDL-C level significantly increased with age (p = 0.01), and that the association of the C allele of rs3135506 with raised triacylglycerol level significantly increased over time (p = 0.0007). The current study suggests that the effects of established gene variants on lipid and glucose traits do not tend to alter with age during adulthood or over time. Show less
The APOA5 gene variants, -1131T>C and S19W, are associated with altered triglyceride concentrations in studies of subjects of Caucasian and East Asian descent. There are few studies of these variants Show more
The APOA5 gene variants, -1131T>C and S19W, are associated with altered triglyceride concentrations in studies of subjects of Caucasian and East Asian descent. There are few studies of these variants in South Asians. We investigated whether the two APOA5 variants also show similar association with various lipid parameters in Indian population as in the UK white subjects. We genotyped 557 Indian adults from Pune, India, and 237 UK white adults for -1131T>C and S19W variants in the APOA5 gene, compared their allelic and genotype frequency and determined their association with fasting serum triglycerides, total cholesterol, HDL and LDL cholesterol levels using univariate general linear analysis. APOC3 SstI polymorphism was also analyzed in 175 Pune Indian subjects for analysis of linkage disequilibrium with the APOA5 variants. The APOA5 -1131C allele was more prevalent in Indians from Pune (Pune Indians) compared to UK white subjects (allele frequency 20% vs. 4%, p = 0.00001), whereas the 19W allele was less prevalent (3% vs. 6% p = 0.0015). Patterns of linkage disequilibrium between the two variants were similar between the two populations and confirmed that they occur on two different haplotypes. In Pune Indians, the presence of -1131C allele and the 19W allele was associated with a 19% and 15% increase respectively in triglyceride concentrations although only -1131C was significant (p = 0.0003). This effect size was similar to that seen in the UK white subjects. Analysis of the APOC3 SstI polymorphism in 175 Pune Indian subjects showed that this variant is not in appreciable linkage disequilibrium with the APOA5 -1131T>C variant (r2 = 0.07). This is the first study to look at the role of APOA5 in Asian Indian subjects that reside in India. The -1131C allele is more prevalent and the 19W allele is less prevalent in Pune Indians compared to UK Caucasians. We confirm that the APOA5 variants are associated with triglyceride levels independent of ethnicity and that this association is similar in magnitude in Asian Indians and Caucasians. The -1131C allele is present in 36% of the Pune Indian population making it a powerful marker for looking at the role of elevated triglycerides in important conditions such as pancreatitis, diabetes and coronary heart disease. Show less