Transethnic insight into the genetics of glycaemic traits: fine-mapping results from the Population Architecture using Genomics and Epidemiology (PAGE) consortium

Diabetologia, Sep 2017

Aims/hypothesis Elevated levels of fasting glucose and fasting insulin in non-diabetic individuals are markers of dysregulation of glucose metabolism and are strong risk factors for type 2 diabetes. Genome-wide association studies have discovered over 50 SNPs associated with these traits. Most of these loci were discovered in European populations and have not been tested in a well-powered multi-ethnic study. We hypothesised that a large, ancestrally diverse, fine-mapping genetic study of glycaemic traits would identify novel and population-specific associations that were previously undetectable by European-centric studies. Methods A multiethnic study of up to 26,760 unrelated individuals without diabetes, of predominantly Hispanic/Latino and African ancestries, were genotyped using the Metabochip. Transethnic meta-analysis of racial/ethnic-specific linear regression analyses were performed for fasting glucose and fasting insulin. We attempted to replicate 39 fasting glucose and 17 fasting insulin loci. Genetic fine-mapping was performed through sequential conditional analyses in 15 regions that included both the initially reported SNP association(s) and denser coverage of SNP markers. In addition, Metabochip-wide analyses were performed to discover novel fasting glucose and fasting insulin loci. The most significant SNP associations were further examined using bioinformatic functional annotation. Results Previously reported SNP associations were significantly replicated (p ≤ 0.05) in 31/39 fasting glucose loci and 14/17 fasting insulin loci. Eleven glycaemic trait loci were refined to a smaller list of potentially causal variants through transethnic meta-analysis. Stepwise conditional analysis identified two loci with independent secondary signals (G6PC2-rs477224 and GCK-rs2908290), which had not previously been reported. Population-specific conditional analyses identified an independent signal in G6PC2 tagged by the rare variant rs77719485 in African ancestry. Further Metabochip-wide analysis uncovered one novel fasting insulin locus at SLC17A2-rs75862513. Conclusions/interpretation These findings suggest that while glycaemic trait loci often have generalisable effects across the studied populations, transethnic genetic studies help to prioritise likely functional SNPs, identify novel associations that may be population-specific and in turn have the potential to influence screening efforts or therapeutic discoveries. Data availability The summary statistics from each of the ancestry-specific and transethnic (combined ancestry) results can be found under the PAGE study on dbGaP here: https://​www.​ncbi.​nlm.​nih.​gov/​projects/​gap/​cgi-bin/​study.​cgi?​study_​id=​phs000356.​v1.​p1

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://link.springer.com/content/pdf/10.1007%2Fs00125-017-4405-1.pdf

Transethnic insight into the genetics of glycaemic traits: fine-mapping results from the Population Architecture using Genomics and Epidemiology (PAGE) consortium

C E = Transethnic insight into the genetics of glycaemic traits: fine-mapping results from the Population Architecture using Genomics and Epidemiology (PAGE) consortium Stephanie A. Bien 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 James S. Pankow 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Jeffrey Haessler 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Yinchang N. Lu 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Nathan Pankratz 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Rebecca R. Rohde 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Alfred Tamuno 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Christopher S. Carlson 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Fredrick R. Schumacher 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Petra Bůžková 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Martha L. Daviglus 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Unhee Lim 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Steven Buyske 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Myron D. Gross 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Mariaelisa Graff 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Carmen R. Isasi 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Lewis H. Kuller 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 JoAnn E. Manson 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Tara C. Matise 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Ross L. Prentice 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Lynne R. Wilkens 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Sachiko Yoneyama 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Ruth J. F. Loos 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Lucia A. Hindorff 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Loic Le Marchand 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Kari E. North 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Christopher A. Haiman 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Ulrike Peters 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 Charles Kooperberg 0 1 2 3 4 5 6 8 9 10 11 12 13 14 15 16 17 18 19 20 21 0 Division of Cardiovascular Sciences, National Heart, Lung, and Blood Institute, National Institutes of Health , Bethesda, MD , USA 1 Department of Biological Sciences, Vanderbilt University , Nashville, TN , USA 2 Division of Epidemiology and Community Health, University of Minnesota , Minneapolis, MN , USA 3 Human Genetics Center, University of Texas Health Science Center at Houston , Houston, TX , USA 4 Department of Statistics, Rutgers University , Newark, NJ , USA 5 Division of Public Health Sciences, Fred Hutchinson Cancer Research Center , 1100 Fairview Ave N., Seattle, WA 98109-1024 , USA 6 Department of Genetics, Rutgers University , Piscataway, NJ , USA 7 Stephanie A. Bien 8 Epidemiology Program, University of Hawaii Cancer Center , Honolulu, HI , USA 9 Department of Medicine, Institute for Minority Health Research, University of Illinois at Chicago , Chicago, IL , USA 10 Department of Biostatistics, University of Washington , Seattle, WA , USA 11 Department of Epidemiology and Biostatistics, Case Western Reserve University , Cleveland, OH , USA 12 Department of Epidemiology & Population Health, Albert Einstein College of Medicine , Bronx, NY , USA 13 The Department of Preventive Medicine, The Icahn School of Medicine at Mount Sinai , New York, NY , USA 14 Department of Epidemiology, School of Public Health, University of North Carolina at Chapel Hill , Chapel Hill, NC , USA 15 Department of Laboratory Medicine and Pathology, University of Minnesota , Minneapolis, MN , USA 16 Department of Ophthalmology and Visual Sciences, University of Michigan , Ann Arbor, MI , USA 17 Department of Epidemiology, Harvard T.H. Chan School of Public Health , Boston, MA , USA 18 Department of Epidemiology, University of Pittsburgh , Pittsburgh, PA , USA 19 Department of Preventive Medicine, Keck School of Medicine, University of Southern California/Norris Comprehensive Cancer Center , Los Angeles, CA , USA 20 MRC Epidemiology Unit, Institute of Metabolic Science, University of Cambridge , Cambridge , UK 21 Department of Epidemiology, University of Michigan , Ann Arbor, MI , USA Aims/hypothesis Elevated levels of fasting glucose and fasting insulin in non-diabetic individuals are markers of dysregulation of glucose metabolism and are strong risk factors for type 2 diabetes. Genome-wide association studies have discovered over 50 SNPs associated with these traits. Most of these loci were discovered in European populations and have not been tested in a well-powered multi-ethnic study. We hypothesised that a large, ancestrally diverse, fine-mapping genetic study of glycaemic traits would identify novel and population-specific associations that were previously undetectable by Europeancentric studies. - Methods A multiethnic study of up to 26,760 unrelated individuals without diabetes, of predominantly Hispanic/Latino a n d A f r i c a n a n c e s t r i e s , w e r e g e n o t y p e d u s i n g t h e Metabochip. Transethnic meta-analysis of racial/ethnicspecific linear regression analyses were performed for fasting glucose and fasting insulin. We attempted to replicate 39 fasting glucose and 17 fasting insulin loci. Genetic finemapping was performed through sequential conditional analyses in 15 regions that included both the initially reported SNP association(s) and denser coverage of SNP markers. In addition, Metabochip-wide analyses were performed to discover novel fasting glucose and fasting insulin loci. The most significant SNP associations were further examined using bioinformatic functional annotation. Results Previously reported SNP associations were significantly replicated (p ≤ 0.05) in 31/39 fasting glucose loci and 14/17 fasting insulin loci. Eleven glycaemic trait loci were refined to a smaller list of potentially causal variants through transethnic meta-analysis. Stepwise conditional analysis identified two loci with independent secondary signals (G6PC2rs477224 and GCK-rs2908290), which had not previously been reported. Population-specific conditional analyses identified an independent signal in G6PC2 tagged by the rare variant rs77719485 in African ancestry. Further Metabochipwide analysis uncovered one novel fasting insulin locus at SLC17A2-rs75862513. Conclusions/interpretation These findings suggest that while glycaemic trait loci often have generalisable effects across the studied populations, transethnic genetic studies help to prioritise likely functional SNPs, identify novel associations that may be population-specific and in turn have the potential to influence screening efforts or therapeutic discoveries. 16 17 18 19 20 24 25 21 The Charles Bronfman Institute for Personalized Medicine, The Icahn School of Medicine at Mount Sinai, New York, NY, USA 22 The Icahn School of Medicine at Mount Sinai, New York, NY, USA 23 National Human Genome Research Institute, National Institutes of Health, Bethesda, MD, USA Carolina Center for Genome Sciences, University of North Carolina at Chapel Hill, Chapel Hill, NC, USA Data availability The summary statistics from each of the ancestry-specific and transethnic (combined ancestry) results can be found under the PAGE study on dbGaP here: https:// www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_ id=phs000356.v1.p1 Abbreviations AA AFR AI/AN AMR ARIC ASN CARDIA CEU EUR GWAS HCHS/SOL H/L MAF MAGIC MEC NHGRI PAGE SHARe WHI Introduction African ancestry African ancestry (1000 Genomes Super Population Code) American Indian/Alaskan Native Admixed American ancestry (1000 Genomes Super Population Code) Atherosclerosis Risk in Communities Asian and Pacific Islander Coronary Artery Risk Development in Young Adults Utah Residents (CEPH) with Northern and Western European Ancestry (HapMap Population Code) European ancestry (1000 Genomes Super Population Code) Genome-wide association studies Hispanic Community Health Study/Study of Latinos Hispanic/Latino Minor allele frequency Meta-Analyses of Glucose and Insulinrelated traits The Multiethnic Cohort National Human Genome Research Institute Population Architecture using Genetic Epidemiology WHI SNP Health Association Resource Women’s Health Initiative Type 2 diabetes is a growing epidemic that disproportionally burdens US minority populations [ 1 ]. Elevated levels of fasting glucose and fasting insulin in individuals without diabetes are markers of dysregulated glucose metabolism and are strong risk factors for type 2 diabetes [ 2 ]. Although twin and family studies provide heritability estimates of 10–50% for these traits [ 3, 4 ], family-based linkage studies have been largely unsuccessful in identifying specific contributing loci. Genome-wide association studies (GWAS) greatly accelerated the pace of discovery of genetic variants contributing to glycaemic traits. For example, the Meta-Analyses of Glucose and Insulin-related traits (MAGIC) consortium performed a large-scale investigation of glycaemic traits in individuals of European descent without diabetes and identified 24 fasting glucose loci and eight fasting insulin loci, three of which were associated with both traits [ 5, 6 ]. These findings have implicated genes and pathways known to be related to glucose metabolism (e.g. GCK/G6PC2 and glucose dephosphorylation), as well as novel pathways (e.g. MTNR1B and circadian rhythmicity). However, in some instances, the interpretation of GWAS findings has been challenging. For instance, many of the known loci are positioned in non-coding, putative regulatory regions of the genome, which in turn makes it difficult to identify the gene target(s). Additionally, the most significant variant is often not the causal variant but is a correlated variant in linkage disequilibrium with the functional variant(s). While early GWAS efforts were focused on populations of European descent, initial attempts to generalise GWAS findings to more diverse populations have had limited success [ 7–9 ]. Importantly, these studies tended to be small and only included the initial most significant GWAS variant (index SNP). However, it is critical that transethnic investigation of GWAS loci include both the index variant and all correlated variants, given that patterns of linkage disequilibrium vary by ancestry and the functional SNP(s) are rarely known. On average, European populations have more highly correlated SNPs and extended haplotypes in comparison with populations of African ancestry (AA). Hispanic/Latino (H/L) populations, on the other hand, are more admixed with highly variable contributions of African, European and New World ancestry. Due in part to reduction in linkage disequilibrium with neighbouring SNPs, transethnic studies can utilise these differences across and within admixed populations to localise causal variants, and discover novel population-specific associations that were undetectable in genetically homogeneous studies. Thus, transethnic studies may provide insight into the underlying biology of complex traits, which may differ among groups. The Metabochip was developed to fine-map GWAS loci for metabolic and cardiovascular traits, as well as replicate promising loci with suggestive, but not genome-wide, significant p values [ 10 ]. Among the 196,725 Metabochip variants selected for fine-mapping metabolic and cardiovascularrelated loci, approximately 40,000 were selected for type 2 diabetes and related biomarkers. Among the 39 fasting glucose loci and 17 fasting insulin loci [ 5, 6 ] that were available for replication, 15 loci included not only the index SNP but also denser coverage of SNPs on the Metabochip that could be utilised for fine-mapping. Importantly, despite very large sample sizes, attempted Metabochip fine-mapping in a population of European descent generally did not yield stronger associations than the original GWAS index SNP and did not reduce the number of SNPs reaching similar levels of significance [11]. As such, this effort was unable to narrow in on functional candidate SNP(s). This study examined the association of Metabochip SNPs with fasting glucose and fasting insulin in a multiethnic study of up to 26,760 participants: 14,953 H/L, 10,380 AA, 998 Asian and Pacific Islander (ASN) and 429 American Indian/Alaskan Native (AI/AN) populations from the Population Architecture using Genetic Epidemiology (PAGE) consortium. Specifically, we carried out the following procedures: (1) tested the association of index SNPs previously reported for 39 fasting glucose and 17 fasting insulin loci from studies of individuals of European descent; (2) used transethnic meta-analysis to refine known glycaemic trait loci in 15 loci which were densely covered with SNPs on the Metabochip; (3) investigated remaining metabolic and cardiovascular trait loci on the Metabochip for association with these glycaemic traits and (4) performed bioinformatic functional annotation of the most significant (lead) SNPs to further prioritise likely causal variants. Methods Ethics statement This study was performed in accordance with the tenets of the Declaration of Helsinki and approved by the Institutional Review Boards of each participating study. All study participants provided written informed consent. Study population and trait measurement The PAGE con sortium was funded by the National Human Genome Research Institute (NHGRI) to investigate the epidemiological architecture of well-replicated genetic variants associated with human diseases or traits [ 12 ]. This analysis includes selfreported H/L, AA, ASN and AI/AN individuals without diabetes, aged 18 years or over, from the Multiethnic Cohort Study (MEC), the Women’s Health Initiative (WHI), Atherosclerosis Risk in Communities (ARIC), Coronary Artery Risk Development in Young Adults (CARDIA), the Hispanic Community Health Study/Study of Latinos (HCHS/SOL) and the Mount Sinai School of Medicine’s (MSSM) DNA biobank (BioMe). Further details about each cohort can be found in the electronic supplementary materials (ESM) Methods (study population and trait measurement section). Fasting glucose and fasting insulin concentrations were measured using standard assays, at laboratories specific to each PAGE site (ESM Table 1). Individuals self-reporting that they had ever been diagnosed with diabetes or taken diabetes medications or who had fasting blood glucose levels ≥ 6.99 mmol/l (≥ 126 mg/dl) were excluded from analyses. Individuals with BMI < 16.5 kg/m2 or BMI > 70 kg/m2 were also excluded on the assumption that these extremes could be attributable to data coding errors or underlying illness or could reflect a familial syndrome. Prior to analyses, each study removed race/ethnicity outliers using ancestry informative principal components. After exclusions, fasting glucose analyses consisted of 14,953 H/L, 10,380 AA, 998 ASN and 429 AI/AN individuals. Fasting insulin analyses involved fewer individuals: 12,895 H/L, 8361 AA, 998 ASN and 420 AI/AN. Fasting insulin was not available for BioMe. Race/ethnicity was self-reported. Descriptive characteristics of PAGE study participants by cohort can be found in ESM Table 2. While ASN and AI/AN were included for transethnic meta-analysis, populationspecific analyses were underpowered due to small sample sizes. As such, ASN and AI/AN population-specific analyses were used as a comparison for consistency in the direction of effect. Genotyping and quality control Genotyping was performed using the Metabochip, the design of which has been described elsewhere [ 10 ]. In brief, the 200K Metabochip is designed to cost effectively analyse putative association signals identified through GWAS of many glucose- and insulin-related metabolic and cardiovascular traits and to fine-map established loci [ 10 ]. More than 122,000 SNPs were included to fine-map 257 GWAS loci for 23 traits [ 10 ]. Fine-mapping loci were defined as the GWAS index SNP and all correlated SNPs (r2 ≥ 0.5) that were within 0.02 cM of the index and having a minor allele frequency (MAF) > 1% in at least one HapMap Phase I population. SNPs were excluded if the Illumina design score was < 0.5 or there were SNPs within 15 bp of the SNP of interest with MAF of > 2% among Europeans (CEU [HapMap Population Code for Utah residents (CEPH) with Northern and Western European ancestry]). Metobochip genotyping was performed for MEC, ARIC, CARDIA, HCHS/SOL and WHI [ 13 ] individuals. Standard quality control filters were applied for samples and SNPs, including missing rate and Hardy–Weinberg equilibrium (p < 1 × 10−7). A portion of WHI individuals of AA had both Metabochip and the Affymetrix 6.0 genotype data available from the SNP Health Association Resource (SHARe); this was used to impute Metabochip SNPs in the remaining SHARe participants with only Affymetrix 6.0 GWAS [ 8 ] and only dosages with imputation R2 > 0.3 were included in the analyses. In BioMe, genotypes from the Illumina HumanOmniExpress array were imputed to 1000 Genome Phase I haplotype panels (March 2012) [ 14 ]. Metabochip SNPs with ‘proper info’ score ≥ 0.4 were included in the analysis. Principal components were determined within each study using the Eigensoft software [ 15 ]. We excluded SNPs with a minor allele count less than 5 within each study by racial/ethnic population. The sample success rate and concordance rate for duplicate pairs across all studies was ≥ 95% and ≥ 99%, respectively. Further genotyping and analytical characteristics of the participating studies are further summarised in ESM Methods (genotyping and quality control section) and ESM Table 1. Replication and fine-mapping approach The overall study design for replication, fine-mapping and discovery of novel loci is summarised in Fig. 1. For replication of known loci, unconditional association analyses were performed for previously reported index SNPs listed in ESM Table 3. A nominal significance level (α = 0.05) was used to define replication of a locus. Next, unconditional association analyses were performed for all SNPs in a locus by race/ethnicity and by transethnic meta-analysis. A locus-specific p value threshold was defined as 0.05 divided by the number of SNPs passing quality control in each region (ranging from α = 1.4 × 10−5 to α = 4.1 × 10−4, Table 1). Locus-specific significance was used to conservatively adjust for multiple testing, while also acknowledging that genetic variation is known to influence glycaemic traits in these regions. Linkage disequilibrium was calculated for PAGE H/L, AA and Asian samples with 500 kb sliding windows using PLINK [ 16 ]. Metabochip linkage disequilibrium and frequency information in Europeans was provided by the 1000 Genomes Phase 3 population. These linkage disequilibrium p a t t e r n s w e r e us e d t o e v a l u a te l o c u s r e f i n e m e n t . Additionally, LocusZoom plots [ 17 ] were used to graphically display the fine-mapping results and linkage disequilibrium for these plots used 1000 Genomes Phase I Super Populations (European ancestry [EUR], admixed American ancestry [AMR], African ancestry [AFR]). After identifying the most significant lead SNP in each region, we searched for additional independent association signals by including the lead SNP in the conditional model and then testing each of the remaining SNPs in a region. These conditional analyses were repeated, adding in the lead SNP and conditional lead SNP(s), until no SNP in the model had a conditional p value less than the locus-specific significance. Sequential conditional analyses were performed for each race/ethnicity and transethnic meta-analysis. Further details on our approach to locus refinement are provided in ESM Methods (replication and fine-mapping of known glycaemic trait loci section). Discovery of novel loci Metabochip-wide analyses were performed to identify novel associations with fasting glucose and fasting insulin. Statistical significance for the Metabochipwide analysis was set at 0.05 divided by the number of Metabochip SNPs passing quality control (α = 2.7 × 10−7). Results were examined through qq plots and Manhattan plots for each model, highlighting known regions defined in ESM Table 4. Further details are provided in ESM Methods (strategy for selecting novel associations section). ilitcaeponR 42FF65aassSSttiiNNnnggPP ssignlsuucloinse Statistical analysis First, in each study with unrelated individuals we performed race/ethnic-specific analyses for fasting glucose and natural log-transformed fasting insulin, excluding ancestry outliers and first-degree relatives. In HCHS/SOL, a weighted version of generalised estimation equations was used to account for unequal inclusion probabilities and complex family-based sampling designs [ 18 ]. Models adjusted for age, sex (except WHI), study site (as applicable), smoking status (current vs former/never), continuous BMI and ancestry principal components. Like previous studies [ 11 ], primary analyses adjusted for BMI because it is a major risk factor for type 2 diabetes and is correlated with glycaemic traits. For comparison, all models were also run without adjustment for BMI. Next, fixed-effect models with inverse-variance weighting were used to pool the study-specific SNP effect estimates and their standard errors by race/ethnicity as implemented in METAL [ 19 ]. Finally, summary statistics from METAL for H/L, AA, NA/AI and ASN were combined using inverse-variance weighted fixed effects meta-analysis in METAL. Q statistics and I2 were used to evaluate heterogeneity across studies and race/ethnicity. Further details are provided in ESM Methods (statistical analysis section). Functional annotation Detailed information on the functional annotation methods and various datasets used is provided in ESM Methods (functional annotation section). In brief, it is expected that the lead SNPs are more likely to be functional or to be in stronger linkage disequilibrium with underlying functional variant(s). Therefore, lead SNPs and all correlated SNPs (r2 > 0.2 in 1000 Genomes Phase 3 AFR/AMR populations) were annotated using publicly available functional datasets. Potential functional effects were assessed using PolyPhen2 [ 20 ] (http://genetics.bwh.harvard.edu/pph2/, accessed 24 August 2016) for non-synonymous variants, SPANR (http:// tools.genes.toronto.edu/) [ 21 ] for variants near splice sites, TargetScan miRNA Regulatory Sites for 3′-UTR regions [ 22 ], ENCODE/NIH Roadmap data [ 23–25 ] and GTEx (https://www.gtexportal.org/home/) [ 26 ] to identify noncoding variants positioned in predicted regulatory elements. Results Demographics We included a total of 26,760 participants (14,953 H/L, 10,380 AA, 998 ASN, and 429 AI/AN) in fasting glucose analyses. The sample sizes for fasting insulin analyses were slightly smaller, with a total of 22,674 participants (12,895 H/L, 8361 AA, 998 ASN and 420 AI/AN). The mean age across the five cohorts was 55 years for men and 59 years for women (range 18–93 years). Study-specific descriptive characteristics are shown in ESM Table 2. Particularly due to the inclusion of the WHI cohort, the proportion of women in the total study population was high, with the highest fraction observed in AA (82.6% for fasting glucose and 97.1% for fasting insulin). Glycaemic trait distributions were similar across studies and ethnicities, with average fasting glucose levels ranging from 4.7 ± 0.7 mmol/l to 5.5 ± 0.6 mmol/l and average fasting insulin levels ranging from 43.3 ± 23.6 pmol/l to 75.9 ± 38.8 pmol/l. Generalisation of European glycaemic trait loci We found that 31/39 (79.5%) fasting glucose loci and 14/17 (82.3%) fasting insulin loci had a p value smaller than 0.05. Index SNP associations were directionally consistent in our transethnic PAGE meta-analysis and only four SNPs had heterogeneity p values less than 0.05 (Table 2). The effect estimates (βs) of index SNPs in the transethnic meta-analysis were very similar to those published in Metabochip analysis of individuals of European descent (Pearson’s r2 = 0.86, 95% CI 0.78, 0.91; p < 2.2 × 10−16; ESM Fig. 1). At three loci (WARS, GIPR and DPYSLS) we observed replication in only H/L and not the transethnic meta-analysis. Interestingly, while the sample sizes were much smaller for Asian individuals than for H/L and AA individuals, the transethnic meta-analysis of the PROX1 index (rs340874) was only nominally significant and directionally consistent in the Asian samples. In the remaining loci that did not replicate in transethnic meta-analysis or the race/ethnic-specific analyses, the effects were generally similar or at least in the same direction. Analyses without inclusion of BMI as a covariate were generally similar, with slightly lower significance at some loci. Full summary statistics for models with and without BMI covariate are reported in ESM Table 5 and ESM Table 6, respectively. Fine-mapping of European glycaemic trait loci Among the 15 glycaemic trait loci for which fine-mapping was attempted on the Metabochip, ten fasting glucose loci and two fasting insulin loci had one or more SNPs that reached locus-specific significance (α = 0.05/number of SNPs in the locus) in the transethnic meta-analysis. The p values ranged from 1.0 × 10−29 at G6PC2-rs560887 to 1.5 × 10−4 at PROX1rs10494973 (Table 3). Although AI/AN ancestries were included in the transethnic meta-analysis, the AI/AN results are not shown because the small sample size was underpowered for population-specific analysis. At four fasting glucose loci, the most significant lead SNP in PAGE transethnic meta-analysis was the same as the European index SNP from prior Metabochip evaluation (G6PC2, ADCY5, MTNR1B and FADS2). For six fasting glucose loci (PROX1, GCKR, SLC2A2, DGKB, GCK and GLIS3) and the one fasting insulin locus (GCKR), the lead SNP in PAGE transethnic meta-analysis was in moderate or weak linkage disequilibrium with the index SNP in 1000 Genomes Population EUR (r2 > 0.2). At these six fasting glucose loci and one fasting insulin locus, the PAGE lead SNP and EUR index SNP were not independent of each other as only one of the two SNP associations maintained nominal significance in transethnic conditional meta-analysis where both lead and index variants were included in the model. This was further supported by investigation of potential fine-mapping through locus zoom plots. For each of the 11 glycaemic trait loci with potential transethnic fine-mapping (fasting glucose loci–PROX1, G6PC2, ADCY5, MTNR1B, FADS2, GCKR, SLC2A2, DGKB, GCK and GLIS3; fasting insulin locus–GCKR), we found that the number of SNPs in linkage disequilibrium with the most significant marker in the transethnic results (r2 ≥ 0.2 in the 1KG super populations AFR and AMR) were less than the number of SNPs tagged by the EUR marker (r2 ≥ 0.2 in EUR). Visual inspection of locus zoom plots indicated that transethnic meta-analysis refined each of these loci by reducing the number of highly correlated SNPs reaching the same level of significance and/or narrowing the genomic region containing putative causal SNPs (ESM Fig. 2). On average, the number of variants in high linkage disequilibrium was reduced by 72.5% with the number of linkage disequilibrium SNPs ranging from one at MTNR1B to 162 at SLC2A2 in the PAGE transethnic metaanalysis results. Refinement was most evident at the SLC2A2 locus (Fig. 2). Bioinformatic functional follow-up was performed for each of the eleven glycaemic trait loci with one or more variants passing the region-specific significance threshold in our transethnic meta-analysis. We observed an overlap of promoter and enhancer sequences at each locus and identified potential target genes. These data not only provided further support for the fine-mapping results but also revealed additional insights into the aetiology of glycaemic traits. UCSC Genome Browser images of each locus are provided ,N /G /G /G /G /G /C /G /G /C /G /T /C /G /G /G /G /G /G /G /C /G 0 A A A A A A A A A A A A A A C A A A A A A n /N lle a e C a p 6 7 , 6 2 o r itfacnouoE eadELRU i(c=NTE rs807434 rs009478 rs808756 rs07806711 rs0812 rs1934921 rs409773 rs55847111 rs18491610 rs91550211 rs6056545 rs06592411 rs30918211 rs545017 rs38096310 rs2015645 rs18166110 rs9627248 rs16931911 rs383099 rs4931536 il lo p e e s R e co 8 B A B 3B le2 /segnu litgugn .323q 1XRO .323p RKC .11q3 26PC .121q 5YCD .262q 22LAC .212p BKG 3p1 KC .114q2 30LAC .224p I3SL .5220q 2AARD .522q0 72LFC .121p1 2YR .121p1 ADD .221q1 2SAD .413q1 1TRN .222q5 42CD .321p 2NKD 51q 1SKC .212q3 1XD .132p 1RPP .1p12 01RB b c s 1 P 2 G 2 G 3 A 3 S 7 D 7 G 8 S 9 G 1 A 1 T 1 C 1 M1 F 1 M1 C 9 C 5 P 1 P 8 P 7 G a o a T L F ) pleauvTE .tt()eaeHM 5−.()006191× 6−.()055012× ..()480150 4−.()091015× ..()400800 ..()550500 5−.()094013× ..()880210 3−.()053015× ..()900700 4−.()008107× ..()560950 ..()970200 ..)(74040 ..()110340 ..()050300 ..()580360 ..()620950 ..()304100 7−.()041102× 10−.(029101× E w T , E T , E T , E T , E E T T , , C le /N lle /G /G /C /G /C /G /C /G /G /G /G /G /G /G /G /G /G /G N /G /G /G C a A A A A A C A A A A A A A A A A A A ,4 A A A 7 6 * , )eud eadLERU rs11336034 rs6172132 rs16369193 rs3810299 rs3734837 rs2359023 rs9322682 rs10707483 rs6027725 rs7609510 rs577466 rs11597115 rs17247654 rs7728085 rs2687579 rs1361714 rs12044695 rs1431821 22=i(cNET rs8466545 rs709804 rs19501252 U E n l o i t n n il le2(co /seegun .1q143 1ARP .1p0211 2XAO .13q3 BAPK .433 LZN .34q22 SRA .1q9332 IRP .22p3 1LAKD .2q2433 22RX q120 1PO .72q2 22FBP .1q31 .1231 TM.42p3 1ERB .13q3 3EBD .31q23 S2L .3p23 S5YLP .25q22 24BCD .11p12 S14R itisnngu 4q1 1YLLPA .2p33 RKC .4q23 41BR b c 1 A 2 F 9 IK 9q D 1 W1 G 6 C 1 P 2 T 3 IG 13 LK p3 A 6 R 5 Z 1 G 2 D 1 C 1 O sa 1 L 2 G 2 G a o F T L ) E .t T e H e ( u l a a t v e p M 2 /g e s l u b c a o T L E E E E E E E E ,/TL ,/TL ,TA ,/TL eno ,TA ll eno ,/TL /L ll ,TA eno ,/TL A p H H A H N A A N H H A A N H ) ) ) ) ) ) ) ) 6 ) ) 6 ) ) ) ) 6 6 5 6 7 6 0 6 6 0 6 5 6 6 0 0 0 tae .(000 .(000 .(000 .(005 .(000 .(000 .(020 .(000 .(000 .(000 .(000 .(090 .(040 .(020 EM .016 .022 .012 .001 .009 .020 .002 .001 .011 .011 .026 .001 .000 .001 T 0 0 0 − 0 0 − 0 0 0 0 − − − ) ) ) ) ) ) ) ) ) ) 7 ) ) 6 7 8 3 ) 2 6 6 8 8 2 9 6 2 2 2 0 7 3 2 .(004 .(006 .(002 .(030 .(002 .(003 .(000 .(090 .(040 .(021 .(005 .(000 .(070 2 6 4 2 9 6 4 3 0 8 3 5 3 SN .60 .30 .02 .00 .00 .00 .00 .00 .00 .00 .02 .00 .00 A 0 0 0 − 0 0 − … − − − 0 − − a t 8 5 7 4 5 7 6 1 9 3 5 7 6 1 E in lo p n i l u ) s SLND (ET% )(571 ()6709 )(299 ()7581 (24961 ()9121 ()2881 )(787 ()7044 )(991 ()6709 ()8046 iitsgnn a f d e AA .010< .402 1 1 .304 .003 .502 .008 1 1 .402 .01< frsom 0 0 /L .10 .19 o ) t 4 l − a e × ic .1 n 4 = 0 1 s i s y l a n a b r t-eam leau teaM 4−01 9−01 29−01 8−01 7−01 5−01 25−01 5−01 7−01 27−01 8−01 5−01 /llfoo c v E × × × × × × × × × × × × i itlssceanggou lll)aee SAAAN ..010010 ..150025 ..070030 ..840079 ..580032 ..270005 ..180002 ..860036 ..910055 ..930095 ..160250 ..740560 liicnudgnBM frandoSAN saeengohPm trfanoSPNm2r.lacc0u2> /ti-scecehnp tttrseaednh ,ttrseanhn a C f ( e c E l c e T itsennP fceoodd /ERUH L ..480003 ..410034 ..310017 ..820075 ..290034 ..500038 ..18002 ..620071 ..660052 ..780079 ..410035 ..830048 itraaeeodm ,//ILHAAN li-ssaany rfo1000m tiiifsscogn iirsuungm rfseoPSN .frra0001o irrf-rrecono ,iliiruqbum ltaeadSN rceenuyqF a/ETC N ./300CG ./920AG ./410AG ./870AG ./440AG ./430TA ./910AG ./670CA ./060AG ./180CG ./920AG ./60CA ltiitieuvvm itiitsecenh ittecanhm tllacecaud ifsaeednm iliiseuqbd teunbhm p***d< iifceconB iseeadgk n n le3tiiifsscognM ieadnoL PEAG SPN illiscecoouggn .rs49497301332q 1XRO rs2603261.323p†RKC .rs60887511q3†26PC .rs70806711121q†5YCD .rs604038126q2†22LAC .rs244861862p12†BKG rs0828629p13†CK .rs097443812p24†I3SL .rs745471221q1†2SAD .rs08309631413q1†1TRBN liiliisconngun .rs2603261332p†RKC .rs860845103222q 1FG itfrffsezceeaaonddm rrsseaeaacdgovFA ttlrrfseaeaehnouvm liiiiirseaequbudgknm ireaendxndpouPSN .liifsaennkgooPSN tiitrreeecaecduongnn p.,.a*001*005n<< ittiiir-fsaacegonnngp ,,,lirseanpouLERD abT egR tsaF 1 P 2 G 2 G 3 A 3 S 7 D 7 G 9 G 1 F 1 MtsaF 2 G 1 I :βaMbp cL dEeNfP p* †S EU Fig. 2 SLC2A2 regional plot. Regional plots of SNP associations (−log10(p value)) with fasting glucose are shown for the MAGIC European (a) and the PAGE transethnic (b) meta-analyses. Not all SNPs used in the transethnic meta-analysis were present in the available MAGIC data (www.magicinvestigators.org/downloads/, accessed 26 June 2017) because of mapping issues [ 11 ]. SNPs not passing quality control or outside the fine-mapping region were removed from the transethnic plots. The colour scale indicates linkage disequilibrium (r2) between each fine-mapping SNP and the GWAS index SNP (rs1280, purple diamond), which was calculated using 1000 Genomes Populations (CEU for MAGIC and AMR for PAGE). The population chosen for linkage disequilibrium colouring in the transethnic metaanalysis was based on population-specific analysis results (choosing the one with strongest underlying SNP associations). The most significant SNPs in MAGIC fine-mapping (rs11709140) and PAGE (rs1604038) are labelled in ESM Fig. 3. The results of our in silico functional annotations are summarised in ESM Table 7. Secondary associations at known glycaemic trait loci To identify additional independent association signals at significant loci, conditional analyses were performed. Results of these analyses and population-specific associations are shown in Table 4. For transethnic conditional meta-analyses, ten fasting glucose loci and two fasting insulin loci were analysed. Independent secondary associations were identified at two fasting glucose loci (G6PC2-rs477224 and GCK-rs2908286). The second round of conditional analyses did not identify significant tertiary signals. Bioinformatic follow-up of rs477224 suggested that the variant is positioned within a pancreatic islet enhancer. The rs2908290 variant was in weak linkage disequilibrium (AMR r2 = 0.26, AFR r2 = 0.23) with a variant, rs2971677, predicted to alter splicing efficiency of GCK. To identify population-specific loci, we conducted separate conditional analyses for significant loci in the primary H/L (GCKR-rs1260326, G6PC2-rs560887, SLC2A2-rs1280, DGKB-rs1005256, GCK-rs1799884, FADS3-rs12577276, MTNR1B-rs10830963, C2CD4A-rs7167881), AA (G6PC2rs77719485, GCK-rs2908286, CRY2-rs117493014, MADDrs77082299, ADCY5-rs11708067, MTNR1B-rs10830963) and Asian populations (GLIS3-rs4395942). A populationspecific variant was detected in the AA analysis of the G6PC2 locus. The lead fasting glucose SNP, rs77719485, is less frequent in AA population (MAF 2.4%) and rare or monomorphic in the other populations (MAF 0.4% in H/L). Like the transethnic lead SNP, rs560887, bioinformatic follow-up suggested that rs77719485 may affect splicing efficiency for exon 4 for G6PC2. Association testing outside of glycaemic trait fine-mapping regions to identify potential novel variants In secondary analyses, we conducted a Metabochip-wide scan to identify potential novel or pleiotropic variants, given that the chip included variants with suggestive signals in established loci for many known metabolic traits. Models were run with and without BMI as a covariate (ESM Table 8, ESM Figs 4,5). Using the Bonferroni significance threshold (0.05/182,055 = 2.7 × 10−7), we identified one novel association for fasting insulin (rs75862513, p = 4.3 × 10−8, Fig. 3) at the SLC17A2 locus previously associated with height and blood pressure [ 27, 28 ]. After BMI adjustment (ESM Fig. 5), the association was attenuated suggesting that the effects may be mediated by BMI. Discussion In this large multiethnic study population of close to 30,000 participants, we used transethnic fine-mapping to narrow the list of putative causal variants for eleven glycaemic trait loci. On average, we observed a 72% reduction in the number of candidate SNPs, before bioinformatic follow-up. We further demonstrated that many of the genetic variants associated with glycaemic traits likely exert their effects through regulatory mechanisms (splicing or enhancer activity), and provide detailed annotations for subsequent laboratory follow-up. These regulatory annotations provide putative targets for laboratory follow-up (e.g. genome editing) and important insights into strong targets for future therapeutic interventions. For example, this study found that most of the implicated enhancer elements were binding sites for the transcription factor FOXA2 in pancreatic islets, and previous studies have suggested that differential expression of FOXA2 is a genetic determinant of fasting glucose levels, as well as type 2 diabetes risk [ 29, 30 ]. Like the previous European Metabochip analysis, we found that rs6113722, which is positioned within a lncRNA adjacent to FOXA2, was e u l a v ./ e) p. nd ray d o n c m o se ir C ( p d 2 r D L b e u l a v p y r a c m P i r N P S 6 6 2 1 7 − − − 0 0 0 1 1 1 2 2 2 7 8 P y r a d n o c e s i r c fo o l e l iln lle u a s ) n i C ( g n d it e s d a o f c d f n o a t e c s fe o f c E * 8 * * * 0 * * * 6 )7 )* )* ic rs5 0 7 2 .(004 .(000 .(002 illno and 5 y F a C N litiltf-ssssaeecaaaougngyn ../rs44057506842277AG ../rs500540439028092AG ltiiiffssceecacogugnAA ../rs7906790375489177CA lilitirfrsssceaeaaadopdynnonowm ,illtirfsssssaaaecn1777y489nuogg ililitirrfsscaeaaacynnghnondooPm jtttrrssaeaceeudondydnohoPSNm liitli)(rrrfacaaanodnonunypoPmm itrrrseeaceaanydnodnypPSNwm illitirfssaaacynnodnoom p.it/rr-fsaeecacnhop1000***dn< iiliirseeubuqdm f e m 2 sp 2 tisecahnn 6PCG CKG lit-anoup 6PCG itleanuq tehAA aedSN laeuv aedSN 2rbD lsaeuv .005< ,ilankg r o e n L p L L p p D T P S I a b c d e * L Fig. 3 Fasting insulin association p values for each Metabochip variant from the transethnic meta-analysis in model without BMI. The –log10 of p values for each SNP on the Metabochip is plotted against chromosomal positions. Grey and black circles, SNPs alternating by chromosome; red squares, SNPs in previously reported glycaemic trait loci (within 1 Mb of index SNP n = 28,580); blue diamonds, novel SNP associations reaching Metabochip-wide significance (all are in the SLC17A2 locus); solid line, threshold for Metabochip-wide significance (0.05/174,898 = 2.9 × 10−7); dashed line, threshold for genome-wide significance α = 5.0 × 10−8 associated (p = 3.2 × 10−8) with fasting glucose. As such, expression levels of FOXA2 could be a particularly important regulator of glucose homeostasis and a putative target for genome editing. Although the clinical application of genome editing is in its infancy, in vivo studies have already demonstrated the utility of the CRISPR/Cas9 technique. For example, to mimic observations of the naturally occurring loss-of-function mutation in the gene encoding LDL receptor antagonist PCSK9, a previous study in mice used CRISPR/ Cas9 vectors to decrease PCSK9 protein levels, which resulted in increased hepatic LDL receptor levels, and a subsequent decrease in blood cholesterol levels [ 31 ]. Identification of key targets, such as FOXA2, and potential regulatory elements of these targets for laboratory follow-up is a critical first step in the translation of GWAS findings. Analysis of known glycaemic trait loci in this diverse population study suggests the genetic determinants of glycaemic trait levels are likely to be similar across populations. In comparison with previous glycaemic trait studies conducted in diverse populations [ 7, 32 ], the replication of effects across populations is more extensive, likely due to the size of this study population. Although most of the loci in the European study were generalisable across populations, this study exemplifies the notion that analysis in diverse populations can refine known loci as well as help in the discovery of novel, sometimes population-specific, associations. For instance, in addition to the well-established splice variant rs560887 that has been robustly associated with fasting glucose, transethnic metaanalysis of the G6PC2 locus identified an additional signal that may implicate regulatory functionality in glycaemia-related tissues. At this same locus, an AA-specific variant, rs77719485, was found to be strongly associated with fasting glucose and, like rs560887 [33], is predicted to affect splicing efficiency. By expanding our analysis to the entire Metabochip, we discovered strong associations with SLC17A2, that were not previously reported by the Metabochip analysis carried out by Scott et al [ 11 ] in Europeans. rs75862513 is a relatively rare variant that appears to be monomorphic in Europeans and was most frequent in the Asian (MAF = 0.04-A) and H/L (MAF = 0.001-A) populations in this study. If replicated in an independent dataset, this finding may represent a new locus not previously detected in European- or AA-specific analyses. These examples illustrate the power of transethnic analysis for locus refinement and novel discovery. Strengths of this study include the large study size, highdensity genotyping and representation of multiple diverse populations. In light of the heavy burden of hyperglycaemia in H/L and AA populations, this study begins to address the major gap in knowledge related to the genetic architecture of glycaemic traits in understudied American minority populations. The large study population, combined with new annotation resources, allowed transethnic fine-mapping and prediction of regulatory elements. However, there were several limitations that should be noted. Although this study included populations from four major racial/ethnic groups, the greatest proportions of participants were H/L and AA. As such, this study was limited in its ability to detect associations with more prominent effects in Asian populations [ 34, 35 ]. We also acknowledge that fine-mapping approaches only serve as an initial step in determining the underlying causal variant(s) driving association signals by prioritising likely causal candidates for more onerous laboratory follow-up. To further meet this objective, functional elements and variants were identified using bioinformatics databases. However, given that the functional evidence detected by these datasets is incomplete, future functional studies are critical in determining the underlying causal variants. That being said, the combination of finemapping with bioinformatics data is particularly useful for reducing both the physical genomic regions of interest and prioritising candidates for molecular characterisation. Furthermore, the in silico approaches help to provide richer inferences regarding the biological mechanisms modulating fasting glucose and insulin levels. As such, fine-mapping is an essential step in functional interpretation of GWAS signals because laboratory follow-up of all possible variants in GWAS loci is prohibitively expensive and time-intensive. This transethnic study comprehensively fine-mapped known common variants associated with concentrations of fasting glucose and insulin. Genomic regions harbouring known risk variants were refined, novel functional candidates were proposed, new independent signals in previously fasting glucose-implicated genes were identified and one novel locus was discovered. Thus, these results suggest that transethnic meta-analysis can help in transforming GWAS results into new biological insight. Acknowledgements The PAGE programme is supported by Genetic Epidemiology of Causal Variants Across the Life Course (CALiCo), MEC, WHI and the Coordinating Center. The contents of this paper are solely the responsibility of the authors and do not necessarily represent the official views of the NIH. The complete list of PAGE members can be found at www.pagestudy.org, accessed 29 April 2016. The data and materials included in this report result from a collaboration between the following studies: (1) The MEC characterisation of epidemiological architecture; (2) The Mount Sinai BioMe Biobank; (3) ‘Epidemiology of putative genetic variants: The Women’s Health Initiative’ study. The authors thank the WHI investigators and staff for their dedication and the study participants for making the program possible. Full listing of WHI investigators can be found at www.whi.org/researchers/Documents% 20%20Write%20a%20Paper/WHI%20Investigator%20Short%20List. pdf, accessed 2 June 2016; (4) The CALiCo programme and the ARIC, CARDIA and HCHS/SOL studies contributed to this manuscript. The authors thank the staff and participants of the ARIC study for their important contributions. Assistance with phenotype harmonisation, SNP selection and annotation, data cleaning, data management, integration and dissemination and general study coordination was provided by the PAGE Coordinating Center. The authors gratefully acknowledge B. Voight for sharing the Metabochip SNP linkage disequilibrium and MAF statistics estimated in the Malmö Diet and Cancer Study. The PAGE consortium thanks the staff and participants of all PAGE studies for their important contributions. Funding The PAGE programme is funded by the NHGRI, supported by U01HG004803 (CALiCo), U01HG004802 (MEC), U01HG004790 (WHI) and U01HG004801 (Coordinating Center), and their respective NHGRI ARRA supplements. The MEC characterisation of epidemiological architecture is funded through the NHGRI PAGE programme (U01HG004802 and its NHGRI ARRA supplement). The MEC study is funded by the National Cancer Institute (R37CA54281, R01 CA63, P01CA33619, U01CA136792 and U01CA98758). Funding support for the ‘Epidemiology of putative genetic variants: The Women’s Health Initiative’ study is provided through the NHGRI PAGE programme (U01HG004790 and its NHGRI ARRA supplement). The WHI programme is funded by the National Heart, Lung, and Blood Institute (NHLBI), the National Institutes of Health (NIH) and the US Department of Health and Human Services through contracts N01WH22110, 24152, 32100-2, 32105-6, 32108-9, 32111-13, 32115, 32118-32119, 32122, 42107-26, 42129-32 and 44221. Funding support for the CALiCo programme was provided by the NHGRI PAGE programme (U01HG004803 and its NHGRI ARRA supplement). The following studies are funded as follows: The ARIC study is carried out as a collaborative study supported by NHLBI contracts ( H H S N 2 6 8 2 0 1 1 0 0 0 0 5 C , H H S N 2 6 8 2 0 1 1 0 0 0 0 6 C , HHSN268201100007C, HHSN268201100008C, HSN268201100009C, H H S N 2 6 8 2 0 1 1 0 0 0 1 0 C , H H S N 2 6 8 2 0 1 1 0 0 0 1 1 C a n d HHSN26 820 11 000 12C), R01 HL087 641, R0 1HL593 67 and R01HL086694; National Human Genome Research Institute contract U 0 1H G 0 04 4 0 2 ; a n d NI H c o n t r a c t H H S N 2 6 8 20 0 6 2 5 2 2 6 C . Infrastructure was partly supported by grant no. UL1RR025005, a component of the NIH and NIH Roadmap for Medical Research. The CARDIA study is supported by contracts HHSN268201300025C, H H S N 2 6 8 2 0 1 3 0 0 0 2 6 C , H H S N 2 6 8 2 0 1 3 0 0 0 2 7 C , H H S N 2 6 8 2 0 1 3 0 0 0 2 8 C , H H S N 2 6 8 2 0 1 3 0 0 0 2 9 C a n d HHSN268200900041C from the NHLBI, the Intramural Research Program of the National Institute on Aging (NIA) and an intra-agency agreement between NIA and NHLBI (AG0005). The HCHS/SOL was carried out as a collaborative study supported by contracts from the NHLBI to the University of North Carolina (N01HC65233), University of Miami (N01-HC65234), Albert Einstein College of Medicine (N01-HC65235), Northwestern University (N01HC65236) and San Diego State University (N01-HC65237). Additional support was provided by 1R01DK101855-01 and 13GRNT16490017. The following Institutes/Centres/Offices contribute to the HCHS/SOL through a transfer of funds to the NHLBI: National Center on Minority Health and Health Disparities, the National Institute of Deafness and Other Communications Disorders, the National Institute of Dental and Craniofacial Research, the National Institute of Diabetes and Digestive and Kidney Diseases, the National Institute of Neurological Disorders and Stroke and the Office of Dietary Supplements. The Mount Sinai BioMe Biobank is supported by The Andrea and Charles Bronfman Philanthropies. Funding support for the PAGE Coordinating Center is provided through the NHGRI PAGE programme (U01HG004801-01 and its NHGRI ARRA supplement). The National Institute of Mental Health also contributes to the support for the Coordinating Center. Duality of interest The authors declare that there is no duality of interest associated with this manuscript. Contribution statement Each of the co-authors made substantial contributions in each of the three following areas: (1) conception and design, acquisition of data or analysis and interpretation of data; (2) drafting the article or revising it critically for important intellectual content and (3) final approval of the version to be published. CK and SAB are responsible for the integrity of the work as a whole. Open Access This article is distributed under the terms of the Creative C o m m o n s A t t r i b u t i o n 4 . 0 I n t e r n a t i o n a l L i c e n s e ( h t t p : / / creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 1. Guariguata L , Whiting DR , Hambleton I et al ( 2014 ) Global estimates of diabetes prevalence for 2013 and projections for 2035 . Diabetes Res Clin Pract 103 : 137 - 149 2. Yki-Järvinen H ( 1994 ) Pathogenesis of non-insulin-dependent diabetes mellitus . Lancet 343 : 91 - 95 3. Almgren P , Lehtovirta M , Isomaa B et al ( 2011 ) Heritability and familiality of type 2 diabetes and related quantitative traits in the Botnia study . Diabetologia 54 : 2811 - 2819 4. Schousboe K , Visscher PM , Henriksen JE et al ( 2003 ) Twin study of genetic and environmental influences on glucose tolerance and indices of insulin sensitivity and secretion . Diabetologia 46 : 1276 - 1283 5. Dupuis J , Langenberg C , Prokopenko I et al ( 2010 ) New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk . Nat Genet 42 : 105 - 116 6. Manning AK , Hivert M-F , Scott RA et al ( 2012 ) A genome-wide approach accounting for body mass index identifies genetic variants influencing fasting glycemic traits and insulin resistance . Nat Genet 44 : 659 - 669 7. Fesinmeyer MD , Meigs JB , North KE et al ( 2013 ) Genetic variants associated with fasting glucose and insulin concentrations in an ethnically diverse population: results from the population architecture using genomics and epidemiology (PAGE) study . BMC Med Genet 14 : 98 8. Liu EY , Buyske S , Aragaki AK et al ( 2012 ) Genotype imputation of Metabochip SNPs using a study-specific reference panel of ~4,000 haplotypes in African Americans from the Women's Health Initiative . Genet Epidemiol 36 : 107 - 117 9. Rasmussen-Torvik LJ , Guo X , Bowden DW et al ( 2012 ) Fasting glucose GWAS candidate region analysis across ethnic groups in the Multiethnic Study of Atherosclerosis (MESA) . Genet Epidemiol 36 : 384 - 391 10. Voight BF , Kang HM , Ding J et al ( 2012 ) The metabochip, a custom genotyping array for genetic studies of metabolic, cardiovascular, and anthropometric traits . PLoS Genet 8 : e1002793 11. Scott RA , Lagou V , Welch RP et al ( 2012 ) Large-scale association analyses identify new loci influencing glycemic traits and provide insight into the underlying biological pathways . Nat Genet 44 : 991 - 1005 12. Matise TC , Ambite JL , Buyske S et al ( 2011 ) The Next PAGE in understanding complex traits: design for the analysis of Population Architecture Using Genetics and Epidemiology (PAGE) Study . Am J Epidemiol 174 : 849 - 859 13. Buyske S , Wu Y , Carty CL et al ( 2012 ) Evaluation of the metabochip genotyping array in African Americans and implications for fine mapping of GWAS-identified loci: the PAGE study . PLoS One 7 : e35651 14. Genomes Project Consortium , Abecasis GR , Auton A et al ( 2012 ) An integrated map of genetic variation from 1,092 human genomes . Nature 491 : 56 - 65 15. Peterson LE ( 2003 ) Partitioning large-sample microarray-based gene expression profiles using principal components analysis . Comput Methods Prog Biomed 70 : 107 - 119 16. Purcell S , Neale B , Todd-Brown K et al ( 2007 ) PLINK: a tool set for whole-genome association and population-based linkage analyses . Am J Hum Genet 81 : 559 - 575 17. Pruim RJ , Welch RP , Sanna S et al ( 2010 ) LocusZoom: regional visualization of genome-wide association scan results . Bioinformatics 26 : 2336 - 2337 18. Lin D-Y , Tao R , Kalsbeek WD et al ( 2014 ) Genetic association analysis under complex survey sampling: the Hispanic Community Health Study/Study of Latinos . Am J Hum Genet 95 : 675 - 688 19. Willer CJ , Li Y , Abecasis GR ( 2010 ) METAL: fast and efficient meta-analysis of genomewide association scans . Bioinformatics 26 : 2190 - 2191 20. Adzhubei IA , Schmidt S , Peshkin L et al ( 2010 ) A method and server for predicting damaging missense mutations . Nat Methods 7 : 248 - 249 21. Xiong HY , Alipanahi B , Lee LJ et al ( 2015 ) RNA splicing . The human splicing code reveals new insights into the genetic determinants of disease . Science 347 : 1254806 22. Friedman RC , Farh KK -H, Burge CB , Bartel DP ( 2009 ) Most mammalian mRNAs are conserved targets of microRNAs . Genome Res 19 : 92 - 105 23. Roadmap Epigenomics Consortium, Kundaje A , Meuleman W et al ( 2015 ) Integrative analysis of 111 reference human epigenomes . Nature 518 : 317 - 330 24. ENCODE Project Consortium ( 2012 ) An integrated encyclopedia of DNA elements in the human genome . Nature 489 : 57 - 74 25. Raney BJ , Dreszer TR , Barber GP et al ( 2014 ) Track data hubs enable visualization of user-defined genome-wide annotations on the UCSC Genome Browser . Bioinformatics 30 : 1003 - 1005 26. GTEx Consortium ( 2013 ) The genotype-tissue expression (GTEx) project . Nat Genet 45 : 580 - 585 27. He M , Xu M , Zhang B et al ( 2015 ) Meta-analysis of genome-wide association studies of adult height in east Asians identifies 17 novel loci . Hum Mol Genet 24 : 1791 - 1800 28. Simino J , Shi G , Bis JC et al ( 2014 ) Gene-age interactions in blood pressure regulation: a large-scale investigation with the CHARGE , Global BPgen , and ICBP Consortia. Am J Hum Genet 95 : 24 - 38 29. Pasquali L , Gaulton KJ , Rodríguez-Seguí SA et al ( 2014 ) Pancreatic islet enhancer clusters enriched in type 2 diabetes riskassociated variants . Nat Genet 46 : 136 - 143 30. Ng MCY , Saxena R , Li J et al ( 2013 ) Transferability and fine mapping of type 2 diabetes loci in African Americans: the Candidate Gene Association Resource Plus Study . Diabetes 62 : 965 - 976 31. Ding Q , Strong A , Patel KM et al ( 2014 ) Permanent alteration of PCSK9 with in vivo CRISPR-Cas9 genome editing . Circ Res 115 : 488 - 492 32. Liu C-T, Ng MCY , Rybin D et al ( 2012 ) Transferability and finemapping of glucose and insulin quantitative trait loci across populations: CARe, the Candidate Gene Association Resource . Diabetologia 55 : 2970 - 2984 33. Baerenwald DA , Bonnefond A , Bouatia-Naji N et al ( 2013 ) Multiple functional polymorphisms in the G6PC2 gene contribute to the association with higher fasting plasma glucose levels . Diabetologia 56 : 1306 - 1316 34. Wang H , Liu L , Zhao J et al ( 2013 ) Large scale meta-analyses of fasting plasma glucose raising variants in GCK, GCKR, MTNR1B and G6PC2 and their impacts on type 2 diabetes mellitus risk . PLoS One 8 : e67665 35. Sakai K , Imamura M , Tanaka Y et al ( 2013 ) Replication study for the association of 9 East Asian GWAS-derived loci with susceptibility to type 2 diabetes in a Japanese population . PLoS One 8 : e76317


This is a preview of a remote PDF: https://link.springer.com/content/pdf/10.1007%2Fs00125-017-4405-1.pdf

Stephanie A. Bien, James S. Pankow, Jeffrey Haessler, Yinchang N. Lu, Nathan Pankratz, Rebecca R. Rohde, Alfred Tamuno, Christopher S. Carlson, Fredrick R. Schumacher, Petra Bůžková, Martha L. Daviglus, Unhee Lim, Myriam Fornage, Lindsay Fernandez-Rhodes, Larissa Avilés-Santa, Steven Buyske, Myron D. Gross, Mariaelisa Graff, Carmen R. Isasi, Lewis H. Kuller, JoAnn E. Manson, Tara C. Matise, Ross L. Prentice, Lynne R. Wilkens, Sachiko Yoneyama, Ruth J. F. Loos, Lucia A. Hindorff, Loic Le Marchand, Kari E. North, Christopher A. Haiman, Ulrike Peters, Charles Kooperberg. Transethnic insight into the genetics of glycaemic traits: fine-mapping results from the Population Architecture using Genomics and Epidemiology (PAGE) consortium, Diabetologia, 2017, 1-15, DOI: 10.1007/s00125-017-4405-1