Andrew Bjonnes

Genome-Wide Association Study Identifies a Novel Locus Contributing to Type 2 Diabetes Susceptibility in Sikhs of Punjabi Origin From India

We performed a genome-wide association study (GWAS) and a multistage meta-analysis of type 2 diabetes (T2D) in Punjabi Sikhs from India. Our discovery GWAS in 1,616 individuals (842 case subjects) was followed by in silico replication of the top 513 independent single nucleotide polymorphisms (SNPs) (P < 10−3) in Punjabi Sikhs (n = 2,819; 801 case subjects). We further replicated 66 SNPs (P < 10−4) through genotyping in a Punjabi Sikh sample (n = 2,894; 1,711 case subjects). On combined meta-analysis in Sikh populations (n = 7,329; 3,354 case subjects), we identified a novel locus in association with T2D at 13q12 represented by a directly genotyped intronic SNP (rs9552911, P = 1.82 × 10−8) in the SGCG gene. Next, we undertook in silico replication (stage 2b) of the top 513 signals (P < 10−3) in 29,157 non-Sikh South Asians (10,971 case subjects) and de novo genotyping of up to 31 top signals (P < 10−4) in 10,817 South Asians (5,157 case subjects) (stage 3b). In combined South Asian meta-analysis, we observed six suggestive associations (P < 10−5 to < 10−7), including SNPs at HMG1L1/CTCFL, PLXNA4, SCAP, and chr5p11. Further evaluation of 31 top SNPs in 33,707 East Asians (16,746 case subjects) (stage 3c) and 47,117 Europeans (8,130 case subjects) (stage 3d), and joint meta-analysis of 128,127 individuals (44,358 case subjects) from 27 multiethnic studies, did not reveal any additional loci nor was there any evidence of replication for the new variant. Our findings provide new evidence on the presence of a population-specific signal in relation to T2D, which may provide additional insights into T2D pathogenesis.

Causal mechanisms and balancing selection inferred from genetic associations with polycystic ovary syndrome

Polycystic ovary syndrome (PCOS) is the most common reproductive disorder in women, yet there is little consensus regarding its aetiology. Here we perform a genome-wide association study of PCOS in up to 5,184 self-reported cases of White European ancestry and 82,759 controls, with follow-up in a further ∼2,000 clinically validated cases and ∼100,000 controls. We identify six signals for PCOS at genome-wide statistical significance (P<5 × 10−8), in/near genes ERBB4/HER4, YAP1, THADA, FSHB, RAD50 and KRR1. Variants in/near three of the four epidermal growth factor receptor genes (ERBB2/HER2, ERBB3/HER3 and ERBB4/HER4) are associated with PCOS at or near genome-wide significance. Mendelian randomization analyses indicate causal roles in PCOS aetiology for higher BMI (P=2.5 × 10−9), higher insulin resistance (P=6 × 10−4) and lower serum sex hormone binding globulin concentrations (P=5 × 10−4). Furthermore, genetic susceptibility to later menopause is associated with higher PCOS risk (P=1.6 × 10−8) and PCOS-susceptibility alleles are associated with higher serum anti-Müllerian hormone concentrations in girls (P=8.9 × 10−5). This large-scale study implicates an aetiological role of the epidermal growth factor receptors, infers causal mechanisms relevant to clinical management and prevention, and suggests balancing selection mechanisms involved in PCOS risk.

Impact of common diabetes risk variant in MTNR1B on sleep, circadian and melatonin physiology

The risk of type 2 diabetes (T2D) is increased by abnormalities in sleep quantity and quality, circadian alignment, and melatonin regulation. A common genetic variant in a receptor for the circadian-regulated hormone melatonin (MTNR1B) is associated with increased fasting blood glucose and risk of T2D, but whether sleep or circadian disruption mediates this risk is unknown. We aimed to test if MTNR1B diabetes risk variant rs10830963 associates with measures of sleep or circadian physiology in intensive in-laboratory protocols (n = 58–96) or cross-sectional studies with sleep quantity and quality and timing measures from self-report (n = 4,307–10,332), actigraphy (n = 1,513), or polysomnography (n = 3,021). In the in-laboratory studies, we found a significant association with a substantially longer duration of elevated melatonin levels (41 min) and delayed circadian phase of dim-light melatonin offset (1.37 h), partially mediated through delayed offset of melatonin synthesis. Furthermore, increased T2D risk in MTNR1B risk allele carriers was more pronounced in early risers versus late risers as determined by 7 days of actigraphy. Our results provide the surprising insight that the MTNR1B risk allele influences dynamics of melatonin secretion, generating a novel hypothesis that the MTNR1B risk allele may extend the duration of endogenous melatonin production later into the morning and that early waking may magnify the diabetes risk conferred by the risk allele.

Genome-Wide Association Study Identifies Variants in Casein Kinase II (CSNK2A2) to be Associated With Leukocyte Telomere Length in a Punjabi Sikh Diabetic Cohort

Background—Telomere length is a heritable trait, and short telomere length has been associated with multiple chronic diseases. We investigated the relationship of relative leukocyte telomere length with cardiometabolic risk and performed the first genome-wide association study and meta-analysis to identify variants influencing relative telomere length in a population of Sikhs from South Asia. Methods and Results—Our results revealed a significant independent association of shorter relative telomere length with type 2 diabetes mellitus and heart disease. Our discovery genome-wide association study (n=1616) was followed by stage 1 replication of 25 top signals (P<10–6) in an additional Sikhs (n=2397). On combined discovery and stage 1 meta-analysis (n= 4013), we identified a novel relative telomere length locus at chromosome 16q21 represented by an intronic variant (rs74019828) in the CSNK2A2 gene (&bgr;=−0.38; P=4.5×10−8). We further tested 3 top variants by genotyping in UK cardiovascular disease (UKCVD) (whites n=2952) for stage 2. Next, we performed in silico replication of 139 top signals (P<10–5) in UK Twin, Nurses Heart Study, Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial, and MD Anderson Cancer Controls (n=10 033) and joint meta-analysis (n=16 998). The observed signal in CSNK2A2 was confined to South Asians and could not be replicated in whites because of significant difference in allele frequencies (P<0.001). CSNK2A2 phosphorylates telomeric repeat binding factor 1 and plays an important role for regulation of telomere length homoeostasis. Conclusions—By identification of a novel signal in telomere pathway genes, our study provides new molecular insight into the underlying mechanism that may regulate telomere length and its association with human aging and cardiometabolic pathophysiology.

Han Chinese polycystic ovary syndrome risk variants in women of European ancestry: relationship to FSH levels and glucose tolerance

STUDY QUESTION Are PCOS risk variants identified in women of Han Chinese ethnicity also associated with risk of PCOS or the phenotypic features of PCOS in European women? SUMMARY ANSWER One variant, rs2268361-T, in the intron of FSHR was associated with PCOS and lower FSH levels, while another variant rs705702-G near the RAB5B and SUOX genes was associated with insulin and glucose levels after oral glucose testing in women with PCOS of European ethnicity. WHAT IS KNOWN ALREADY Three of the eleven variants associated with PCOS in the Han Chinese genome-wide association studies were also associated with PCOS in at least one European population when corrected for multiple testing (DENND1A, THADA and YAP1). However, additional replication is needed to establish the importance of these variants in European women and to determine the relationship to PCOS phenotypic traits. STUDY DESIGN, SIZE, DURATION The study was a case-control examination in a discovery cohort of women with PCOS (n = 485) and controls (n = 407) from Boston (Boston 1). Replication was performed in women from Greece (cases n = 884 and controls n = 311) and an additional cohort from Boston (Boston electronic medical record (EMR); n = 350 cases and n = 1258 controls). PARTICIPANTS/MATERIALS, SETTINGS, METHODS Women had PCOS defined by the National Institutes of Health criteria in Boston 1 and Greece (n = 783), with additional subjects fulfilling the Rotterdam criteria (hyperandrogenism, polycystic ovary morphology and regular menses) in Greece (n = 101). Controls in Boston and Greece had regular menstrual cycles and no hyperandrogenism. The second cohort from Boston was defined using the EMR and natural language processing. Allele frequencies for variants associated with PCOS in Han Chinese women were examined in PCOS cases and controls, along with the relationship to quantitative traits. MAIN RESULTS AND THE ROLE OF CHANCE A variant rs2268361-T in an intron of FSHR was associated with PCOS (0.84 [0.76-0.93], OR [95% CI]; P = 0.002). The rs2268361-T was associated with lower FSH levels (-0.15 ± 0.05; P = 0.0029). A variant rs705702-G near RAB5B and SUOX was associated with insulin (-0.16 ± 0.05, P = 0.0029) and glucose levels (-0.20 ± 0.05, P = 0.0002) 120 min after an oral glucose test. LIMITATIONS, REASONS FOR CAUTION The study was large and contained replication cohorts, but was limited by a small number of controls in the Greek cohort and a small number of cases in the second Boston cohort. The second Boston group was identified using electronic medical record review, but was validated for the cardinal features of PCOS. WIDER IMPLICATIONS OF THE FINDINGS This study demonstrates a cross-ethnic PCOS risk locus in FSHR in women of European ancestry with PCOS. The variant may influence FSH receptor responsiveness as suggested by the associated change in FSH levels. The relationship between a variant near RAB5B and SUOX and glucose stimulated insulin and glucose levels suggests an influence of one of these genes on glucose tolerance, but the absence of a relationship with PCOS points to potential differences in the international PCOS patient populations. STUDY FUNDING/COMPETING INTERESTS The project was supported by Award Number R01HD065029 from the Eunice Kennedy Shriver National Institute Of Child Health & Human Development, Award Number 1 UL1 RR025758, Harvard Clinical and Translational Science Center, from the National Center for Research Resources, award 1-10-CT-57 from the American Diabetes Association and the Partners Healthcare Center for Personalized Genetics Project Grant. C.K.W. is a consultant for Takeda Pharmaceuticals. TRIAL REGISTRATION NUMBER NCT00166569.

Genome-wide association study of 25(OH) Vitamin D concentrations in Punjabi Sikhs: Results of the Asian Indian diabetic heart study.

Vitamin D deficiency is implicated in multiple disease conditions and accumulating evidence supports that the variation in serum vitamin D (25(OH)D) levels, including deficiency, is under strong genetic control. However, the underlying genetic mechanism associated with vitamin 25(OH)D concentrations is poorly understood. We earlier reported a very high prevalence of vitamin D deficiency associated with an increased risk for type 2 diabetes and obesity in a Punjabi Sikh diabetic cohort as part of the Asian Indian diabetic heart study (AIDHS). Here we have performed the first genome-wide association study (GWAS) of serum 25(OH)D on 3538 individuals from this Punjabi Sikh population. Our discovery GWAS comprised of 1387 subjects followed by validation of 24 putative SNPs (P<10(-4)) using an independent replication sample (n=2151) from the same population by direct genotyping. A novel locus at chromosome 20p11.21 represented by rs2207173 with minor allele frequency (MAF) 0.29, [β=-0.13, p=4.47×10(-9)] between FOXA2 and SSTR4 was identified to be associated with 25(OH)D levels. Another suggestive association signal at rs11586313 (MAF 0.54) [β=0.90; p=1.36×10(-6)] was found within the regulatory region of the IVL gene on chromosome 1q21.3. Additionally, our study replicated 3 of 5 known GWAS genes associated with 25(OH)D concentrations including GC (p=0.007) and CYP2R1 (p=0.019) reported in Europeans and the DAB1 (p=0.003), reported in Hispanics. Identification of novel association signals in biologically plausible regions with 25(OH)D metabolism will provide new molecular insights on genetic drivers of vitamin D status and its implications in health disparities.

Gene variants associated with age at menopause are also associated with polycystic ovary syndrome, gonadotrophins and ovarian volume

STUDY QUESTION Is there a relationship between the genetic risk for polycystic ovary syndrome (PCOS) and genetic variants that influence timing of menopause? SUMMARY ANSWER The genetic risk score, which sums the contribution of variants at all menopause loci, was associated with PCOS. WHAT IS ALREADY KNOWN Ovarian parameters and anti-Mullerian hormone levels suggest that women with PCOS should have a later age at menopause. STUDY DESIGN, SIZE, DURATION The study was a case-control examination of genetic variants associated with age at menopause in a discovery cohort of women with PCOS (n = 485) and controls (n = 407) from Boston recruited from 2003 to 2012. Replication was performed in women from Greece (cases, n = 884 and controls, n = 311). PARTICIPANTS/MATERIALS, SETTINGS, METHODS PCOS was defined by the National Institutes of Health criteria in Boston and Greece (n = 783), with additional subjects fulfilling the Rotterdam criteria (hyperandrogenism, polycystic ovary morphology and regular menses) in Greece (n = 101). Controls in Boston and Greece had regular menstrual cycles and no hyperandrogenism. Allele frequencies for variants previously associated with age at menopause were examined in PCOS cases and controls, along with the relationship to quantitative traits. MAIN RESULTS AND ROLE OF CHANCE The variant rs11668344-G was associated with decreased risk of PCOS (odds ratio: 0.77 [0.59-0.93]; P = 0.004). There was a strong relationship between the late menopause allele rs12294104-T and increased LH levels (β ± SE; 0.26 ± 0.06; P = 5.2 × 10(-5)) and the LH:FSH ratio (0.28 ± 0.06; P = 2.7 × 10(-6)). The minor allele at rs10852344-T was associated with smaller ovarian volume (-0.16 ± 0.05; P = 0.0012). A genetic risk score calculated from 16 independent variants associated with age at menopause was also associated with PCOS (P < 0.02), LH and the LH:FSH ratio (both P < 0.05). LIMITATIONS, REASONS FOR CAUTION The variant rs11668344 was not associated with PCOS in the Greek cohort, but results exhibited the same direction of effect as the Boston cohort. However, it is possible that the individual association was a false positive in the Boston cohort. WIDER IMPLICATIONS OF THE FINDINGS The study demonstrates that gene variants known to influence age at menopause are also associated with risk for PCOS. Further, our data suggest that the relationship between age at menopause and PCOS may be explained, at least in part, by effects on LH levels and follicle number. The data point to opposing influences of the genetic variants on both menopausal age and PCOS. STUDY FUNDING/COMPETING INTERESTS The project was supported by award number R01HD065029 from the Eunice Kennedy Shriver National Institute of Child Health & Human Development, award number 1 UL1 RR025758, Harvard Clinical and Translational Science Center, from the National Center for Research Resources and award 1-10-CT-57 from the American Diabetes Association. C.K.W. is a consultant for Takeda Pharmaceuticals. TRIAL REGISTRATION NUMBER NCT00166569.

No Association of Coronary Artery Disease with X-Chromosomal Variants in Comprehensive International Meta-Analysis.

In recent years, genome-wide association studies have identified 58 independent risk loci for coronary artery disease (CAD) on the autosome. However, due to the sex-specific data structure of the X chromosome, it has been excluded from most of these analyses. While females have 2 copies of chromosome X, males have only one. Also, one of the female X chromosomes may be inactivated. Therefore, special test statistics and quality control procedures are required. Thus, little is known about the role of X-chromosomal variants in CAD. To fill this gap, we conducted a comprehensive X-chromosome-wide meta-analysis including more than 43,000 CAD cases and 58,000 controls from 35 international study cohorts. For quality control, sex-specific filters were used to adequately take the special structure of X-chromosomal data into account. For single study analyses, several logistic regression models were calculated allowing for inactivation of one female X-chromosome, adjusting for sex and investigating interactions between sex and genetic variants. Then, meta-analyses including all 35 studies were conducted using random effects models. None of the investigated models revealed genome-wide significant associations for any variant. Although we analyzed the largest-to-date sample, currently available methods were not able to detect any associations of X-chromosomal variants with CAD.

Circadian gene variants influence sleep and the sleep electroencephalogram in humans.

ABSTRACT The sleep electroencephalogram (EEG) is highly heritable in humans and yet little is known about the genetic basis of inter-individual differences in sleep architecture. The aim of this study was to identify associations between candidate circadian gene variants and the polysomnogram, recorded under highly controlled laboratory conditions during a baseline, overnight, 8 h sleep opportunity. A candidate gene approach was employed to analyze single-nucleotide polymorphisms from five circadian-related genes in a two-phase analysis of 84 healthy young adults (28 F; 23.21 ± 2.97 years) of European ancestry. A common variant in Period2 (PER2) was associated with 20 min less slow-wave sleep (SWS) in carriers of the minor allele than in noncarriers, representing a 22% reduction in SWS duration. Moreover, spectral analysis in a subset of participants (n = 37) showed the same PER2 polymorphism was associated with reduced EEG power density in the low delta range (0.25–1.0 Hz) during non-REM sleep and lower slow-wave activity (0.75–4.5 Hz) in the early part of the sleep episode. These results indicate the involvement of PER2 in the homeostatic process of sleep. Additionally, a rare variant in Melatonin Receptor 1B was associated with longer REM sleep latency, with minor allele carriers exhibiting an average of 65 min (87%) longer latency from sleep onset to REM sleep, compared to noncarriers. These findings suggest that circadian-related genes can modulate sleep architecture and the sleep EEG, including specific parameters previously implicated in the homeostatic regulation of sleep.

Genetic Predisposition to Dyslipidemia and Risk of Preeclampsia

BACKGROUND Large epidemiologic studies support the role of dyslipidemia in preeclampsia; however, the etiology of preeclampsia or whether dyslipidemia plays a causal role remains unclear. We examined the association between the genetic predisposition to dyslipidemia and risk of preeclampsia using validated genetic markers of dyslipidemia. METHODS Preeclampsia cases (n = 164) and normotensive controls (n = 110) were selected from live birth certificates to nulliparous Iowa women during the period August 2002 to May 2005. Disease status was verified by medical chart review. Genetic predisposition to dyslipidemia was estimated by 4 genetic risk scores (GRS) (total cholesterol (TC), LDL cholesterol (LDL-C), HDL cholesterol (HDL-C), and triglycerides) on the basis of established loci for blood lipids. Logistic regression analyses were used to evaluate the relationships between each of the 4 genotype scores and preeclampsia. Replication analyses were performed in an independent, US population of preeclampsia cases (n = 516) and controls (n = 1,097) of European ancestry. RESULTS The GRS related to higher levels of TC, LDL-C, and triglycerides demonstrated no association with the risk of preeclampsia in either the Iowa or replication population. The GRS related to lower HDL-C was marginally associated with an increased risk for preeclampsia (odds ratio (OR) = 1.03, 95% confidence interval (CI) = 0.99-1.07; P = 0.10). In the independent replication population, the association with the HDL-C GRS was also marginally significant (OR = 1.03, 95% CI: 1.00-1.06; P = 0.04). CONCLUSIONS Our data suggest a potential effect between the genetic predisposition to dyslipidemic levels of HDL-C and an increased risk of preeclampsia, and, as such, suggest that dyslipidemia may be a component along the causal pathway to preeclampsia.