Genetic drivers of heterogeneity in type 2 diabetes pathophysiology

Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes1,2 and molecular mechanisms that are often specific to cell type3,4. Here, to characterize the genetic contribution to these processes across ancestry groups, we aggregate genome-wide association study data from 2,535,601 individuals (39.7% not of European ancestry), including 428,452 cases of T2D. We identify 1,289 independent association signals at genome-wide significance (P < 5 × 10-8) that map to 611 loci, of which 145 loci are, to our knowledge, previously unreported. We define eight non-overlapping clusters of T2D signals that are characterized by distinct profiles of cardiometabolic trait associations. These clusters are differentially enriched for cell-type-specific regions of open chromatin, including pancreatic islets, adipocytes, endothelial cells and enteroendocrine cells. We build cluster-specific partitioned polygenic scores5 in a further 279,552 individuals of diverse ancestry, including 30,288 cases of T2D, and test their association with T2D-related vascular outcomes. Cluster-specific partitioned polygenic scores are associated with coronary artery disease, peripheral artery disease and end-stage diabetic nephropathy across ancestry groups, highlighting the importance of obesity-related processes in the development of vascular outcomes. Our findings show the value of integrating multi-ancestry genome-wide association study data with single-cell epigenomics to disentangle the aetiological heterogeneity that drives the development and progression of T2D. This might offer a route to optimize global access to genetically informed diabetes care.

Rapid switching from levetiracetam to brivaracetam in pharmaco-resistant epilepsy in people with and without intellectual disabilities: a naturalistic case control study

BACKGROUND

Approximately one quarter of people with an intellectual disability (PwID) have epilepsy of whom nearly three-quarters are pharmaco-resistant. There are higher reported neuropsychiatric side-effects to anti-seizure medication (ASM) in this group. Levetiracetam (LEV) is a first-line ASM with a stronger association with neuropsychiatric symptoms for PwID than other ASMs. Brivaracetam (BRV) is a newer ASM. Recent studies suggest a beneficial effect of swapping people who experience neuropsychiatric events with LEV to BRV. However, there is limited evidence of this for PwID. This evaluation analyses real world outcomes of LEV to BRV swap for PwID compared to those without ID.

METHODS

We performed a multicentre, retrospective review of clinical records. Demographic, clinical characteristics and reported adverse events of patients switched from LEV to BRV (2016-2020) were recorded at 3 months pre and 6- and 12-month post-BRV initiation. Outcomes were compared between PwID and those without and summarised using cross-tabulations and logistic regression models. A Bonferroni correction was applied.

RESULTS

Of 77 participants, 46 had ID and 52% had a past psychiatric illness. 71% participants switched overnight from LEV to BRV. Seizure reduction of > 50% was seen in 40% patients. Psychiatric illness history was predictive of having neuropsychiatric side-effects with LEV but not BRV (p = 0.001). There was no significant difference for any primary outcomes between PwID versus without ID.

CONCLUSIONS

Switching from LEV to BRV appears as well tolerated and efficacious in PwID as those without ID with over 90% still on BRV after 12 months.

Multi-ancestry genome-wide study in >2.5 million individuals reveals heterogeneity in mechanistic pathways of type 2 diabetes and complications

Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes. To characterise the genetic contribution to these processes across ancestry groups, we aggregate genome-wide association study (GWAS) data from 2,535,601 individuals (39.7% non-European ancestry), including 428,452 T2D cases. We identify 1,289 independent association signals at genome-wide significance (P<5×10-8) that map to 611 loci, of which 145 loci are previously unreported. We define eight non-overlapping clusters of T2D signals characterised by distinct profiles of cardiometabolic trait associations. These clusters are differentially enriched for cell-type specific regions of open chromatin, including pancreatic islets, adipocytes, endothelial, and enteroendocrine cells. We build cluster-specific partitioned genetic risk scores (GRS) in an additional 137,559 individuals of diverse ancestry, including 10,159 T2D cases, and test their association with T2D-related vascular outcomes. Cluster-specific partitioned GRS are more strongly associated with coronary artery disease and end-stage diabetic nephropathy than an overall T2D GRS across ancestry groups, highlighting the importance of obesity-related processes in the development of vascular outcomes. Our findings demonstrate the value of integrating multi-ancestry GWAS with single-cell epigenomics to disentangle the aetiological heterogeneity driving the development and progression of T2D, which may offer a route to optimise global access to genetically-informed diabetes care.

Entropic comparison of Landau-Zener and Demkov interactions in the phase space of a quadrupole billiard

We investigate two types of avoided crossings in a chaotic billiard within the framework of information theory. The Shannon entropy in the phase space for the Landau-Zener interaction increases as the center of the avoided crossing is approached, whereas for the Demkov interaction, the Shannon entropy decreases as the center of avoided crossing is passed by with an increase in the deformation parameter. This feature can provide a new indicator for scar formation. In addition, it is found that the Fisher information of the Landau-Zener interaction is significantly larger than that of the Demkov interaction.

Multi-ancestry genetic study of type 2 diabetes highlights the power of diverse populations for discovery and translation

We assembled an ancestrally diverse collection of genome-wide association studies (GWAS) of type 2 diabetes (T2D) in 180,834 affected individuals and 1,159,055 controls (48.9% non-European descent) through the Diabetes Meta-Analysis of Trans-Ethnic association studies (DIAMANTE) Consortium. Multi-ancestry GWAS meta-analysis identified 237 loci attaining stringent genome-wide significance (P < 5 × 10-9), which were delineated to 338 distinct association signals. Fine-mapping of these signals was enhanced by the increased sample size and expanded population diversity of the multi-ancestry meta-analysis, which localized 54.4% of T2D associations to a single variant with >50% posterior probability. This improved fine-mapping enabled systematic assessment of candidate causal genes and molecular mechanisms through which T2D associations are mediated, laying the foundations for functional investigations. Multi-ancestry genetic risk scores enhanced transferability of T2D prediction across diverse populations. Our study provides a step toward more effective clinical translation of T2D GWAS to improve global health for all, irrespective of genetic background.

Unraveling the Genomic Architecture of the CYP3A Locus and ADME Genes for Personalized Tacrolimus Dosing

BACKGROUND

Tacrolimus (TAC) is an immunosuppressant widely prescribed following an allogenic organ transplant. Due to wide interindividual pharmacokinetic (PK) variability, optimizing TAC dosing based on genetic factors is required to minimize nephrotoxicity and acute rejections.

METHODS

We enrolled 1133 participants receiving TAC from 4 cohorts, consisting of 3 with kidney transplant recipients and 1 with healthy males from clinical trials. The effects of clinical factors were estimated to appropriately control confounding variables. A genome-wide association study, haplotype analysis, and a gene-based association test were conducted using the Korea Biobank Array or targeted sequencing for 114 pharmacogenes.

RESULTS

Genome-wide association study verified that CYP3A5*3 is the only common variant associated with TAC PK variability in Koreans. We detected several CYP3A5 and CYP3A4 rare variants that could potentially affect TAC metabolism. The haplotype structure of CYP3A5 stratified by CYP3A5*3 was a significant factor for CYP3A5 rare variant interpretation. CYP3A4 rare variant carriers among CYP3A5 intermediate metabolizers displayed higher TAC trough levels. Gene-based association tests in the 61 absorption, distribution, metabolism, and excretion genes revealed that CYP1A1 are associated with additional TAC PK variability: CYP1A1 rare variant carriers among CYP3A5 poor metabolizers showed lower TAC trough levels than the noncarrier controls.

CONCLUSIONS

Our study demonstrates that rare variant profiling of CYP3A5 and CYP3A4, combined with the haplotype structures of CYP3A locus, provide additive value for personalized TAC dosing. We also identified a novel association between CYP1A1 rare variants and TAC PK variability in the CYP3A5 nonexpressers that needs to be further investigated.

The trans-ancestral genomic architecture of glycemic traits

Glycemic traits are used to diagnose and monitor type 2 diabetes and cardiometabolic health. To date, most genetic studies of glycemic traits have focused on individuals of European ancestry. Here we aggregated genome-wide association studies comprising up to 281,416 individuals without diabetes (30% non-European ancestry) for whom fasting glucose, 2-h glucose after an oral glucose challenge, glycated hemoglobin and fasting insulin data were available. Trans-ancestry and single-ancestry meta-analyses identified 242 loci (99 novel; P < 5 × 10-8), 80% of which had no significant evidence of between-ancestry heterogeneity. Analyses restricted to individuals of European ancestry with equivalent sample size would have led to 24 fewer new loci. Compared with single-ancestry analyses, equivalent-sized trans-ancestry fine-mapping reduced the number of estimated variants in 99% credible sets by a median of 37.5%. Genomic-feature, gene-expression and gene-set analyses revealed distinct biological signatures for each trait, highlighting different underlying biological pathways. Our results increase our understanding of diabetes pathophysiology by using trans-ancestry studies for improved power and resolution.

Genome-wide association studies identify two novel loci conferring susceptibility to diabetic retinopathy in Japanese patients with type 2 diabetes

Several reports have suggested that genetic susceptibility contributes to the development and progression of diabetic retinopathy. We aimed to identify genetic loci that confer susceptibility to diabetic retinopathy in Japanese patients with type 2 diabetes. We analysed 5 790 508 single nucleotide polymorphisms (SNPs) in 8880 Japanese patients with type 2 diabetes, 4839 retinopathy cases and 4041 controls, as well as 2217 independent Japanese patients with type 2 diabetes, 693 retinopathy cases and 1524 controls. The results of these two genome-wide association studies (GWAS) were combined with an inverse variance meta-analysis (Stage-1), followed by de novo genotyping for the candidate SNP loci (P < 1.0 × 10-4) in an independent case-control study (Stage-2, 2260 cases and 723 controls). After combining the association data (Stages 1 and 2) using meta-analysis, the associations of two loci reached a genome-wide significance level: rs12630354 near STT3B on chromosome 3, P = 1.62 × 10-9, odds ratio (OR) = 1.17, 95% confidence interval (CI) 1.11-1.23, and rs140508424 within PALM2 on chromosome 9, P = 4.19 × 10-8, OR = 1.61, 95% CI 1.36-1.91. However, the association of these two loci was not replicated in Korean, European or African American populations. Gene-based analysis using Stage-1 GWAS data identified a gene-level association of EHD3 with susceptibility to diabetic retinopathy (P = 2.17 × 10-6). In conclusion, we identified two novel SNP loci, STT3B and PALM2, and a novel gene, EHD3, that confers susceptibility to diabetic retinopathy; however, further replication studies are required to validate these associations.

The upgraded TOMAS device: A toroidal plasma facility for wall conditioning, plasma production, and plasma-surface interaction studies

The Toroidal Magnetized System device has been significantly upgraded to enable development of various wall conditioning techniques, including methods based on ion and electron cyclotron (IC/EC) range of frequency plasmas, and to complement plasma-wall interaction research in tokamaks and stellarators. The toroidal magnetic field generated by 16 coils can reach its maximum of 125 mT on the toroidal axis. The EC system is operated at 2.45 GHz with up to 6 kW forward power. The IC system can couple up to 6 kW in the frequency range of 10 MHz-50 MHz. The direct current glow discharge system is based on a graphite anode with a maximum voltage of 1.5 kV and a current of 6 A. A load-lock system with a vertical manipulator allows exposure of material samples. A number of diagnostics have been installed: single- and triple-pin Langmuir probes for radial plasma profiles, a time-of-flight neutral particle analyzer capable of detecting neutrals in the energy range of 10 eV-1000 eV, and a quadrupole mass spectrometer and video systems for plasma imaging. The majority of systems and diagnostics are controlled by the Siemens SIMATIC S7 system, which also provides safety interlocks.

Extracting hypertrophic cardiomyopathy features from cardiac magnetic resonance reports by natural language processing

Objective

Determine if information regarding hypertrophic cardiomyopathy (HCM) can be accurately retrieved from cardiac magnetic resonance (CMR) reports using natural language processing (NLP).

Background

CMR imaging is used for diagnosis and risk stratification of HCM. Manual annotation of information from CMR is time-consuming. NLP is an artificial intelligence method for automating extraction of information from narrative text.

Methods

We identified 200 HCM patients who had CMR reports from 1998 to 2018. These patients were randomly allocated into training (100 patients with 185 CMR reports) and testing sets (100 patients with 206 reports). An NLP system with 2 tiers was developed; the first extracted information regarding HCM diagnosis while second extracted categorical or numeric concepts for HCM classification. NLP performance was compared with gold-standard manual annotation.

Results

NLP algorithms achieved very high performance across all concepts with mean positive predictive value (PPV) = 0.96. An outlier was the performance for abstracting the presence of an apical pouch from CMR reports, which had noticeably lower PPV= 0.78 which be attributed to the low number of cases with this finding.

Conclusions

The algorithms developed can be translated to clinical decision support systems to increase efficiency and contribute to improved quality of care.

Funding Acknowledgement

Type of funding source: Other. Main funding source(s): Study supported by the National Heart, Lung and Blood Institute of National Institutes of Health (K01HL124045), the Mayo Clinic Center for Clinical and Translational Science (CCaTS), and the Mayo Clinic K2R award. Content is solely the responsibility of authors and does not necessarily represent official views of the National Institutes of Health.

Obesity susceptible novel DNA methylation marker on regulatory region of inflammation gene: results from the Korea Epigenome Study (KES)

INTRODUCTION

Obesity is growing global health concern and highly associated with increased risk of metabolic diseases including type 2 diabetes. We aimed to discover new differential DNA methylation patterns predisposing obesity and prioritize surrogate epigenetic markers in Koreans.

RESEARCH DESIGN AND METHODS

We performed multistage epigenome-wide analyses to identify differentially expressed CpGs in obesity using the Illumina HumanMethylationEPIC array (EPIC). Forty-eight CpGs showed significant differences across three phases: 902 whole blood DNAs from two cohorts (phase 1: n=450, phase 2: n=377) and a hospital-based sample (phase 3: n=75). Samples from phase III participants were used to examine whether the 48 CpGs are significant in the fat tissue and influenced gene expression. Furthermore, we investigated the epigenetic effect of CpG loci in childhood obesity (n=94).

RESULTS

Seven of the 48 CpGs exhibited similar changes in the fat tissue along with gene expression changes. In particular, hypomethylated CpG (cg13424229) on the GATA1 transcription factor cluster of CPA3 promoter was related to its increased gene expression and showed consistent effect in childhood obesity. Interestingly, subsequent analysis using RNA sequencing data from 21 preadipocytes and 26 adipocytes suggested CPA3 as a potential obesity-related gene. Moreover, expression patterns from RNA sequencing and public Gene Expression Omnibus showed the correlation between CPA3 and type 2 diabetes (T2D) and asthma.

CONCLUSIONS

Our finding prioritizes influential genes in obesity and provides new evidence for the role of CPA3 linking obesity, T2D, and asthma.

Integrated Analysis of Tissue-Specific Promoter Methylation and Gene Expression Profile in Complex Diseases

This study investigated whether the promoter region of DNA methylation positively or negatively regulates tissue-specific genes (TSGs) and if it correlates with disease pathophysiology. We assessed tissue specificity metrics in five human tissues, using sequencing-based approaches, including 52 whole genome bisulfite sequencing (WGBS), 52 RNA-seq, and 144 chromatin immunoprecipitation sequencing (ChIP-seq) data. A correlation analysis was performed between the gene expression and DNA methylation levels of the TSG promoter region. The TSG enrichment analyses were conducted in the gene-disease association network (DisGeNET). The epigenomic association analyses of CpGs in enriched TSG promoters were performed using 1986 Infinium MethylationEPIC array data. A correlation analysis showed significant associations between the promoter methylation and 449 TSGs' expression. A disease enrichment analysis showed that diabetes- and obesity-related diseases were high-ranked. In an epigenomic association analysis based on obesity, 62 CpGs showed statistical significance. Among them, three obesity-related CpGs were newly identified and replicated with statistical significance in independent data. In particular, a CpG (cg17075888 of PDK4), considered as potential therapeutic targets, were associated with complex diseases, including obesity and type 2 diabetes. The methylation changes in a substantial number of the TSG promoters showed a significant association with metabolic diseases. Collectively, our findings provided strong evidence of the relationship between tissue-specific patterns of epigenetic changes and metabolic diseases.

Identification of type 2 diabetes loci in 433,540 East Asian individuals

Meta-analyses of genome-wide association studies (GWAS) have identified more than 240 loci that are associated with type 2 diabetes (T2D)1,2; however, most of these loci have been identified in analyses of individuals with European ancestry. Here, to examine T2D risk in East Asian individuals, we carried out a meta-analysis of GWAS data from 77,418 individuals with T2D and 356,122 healthy control individuals. In the main analysis, we identified 301 distinct association signals at 183 loci, and across T2D association models with and without consideration of body mass index and sex, we identified 61 loci that are newly implicated in predisposition to T2D. Common variants associated with T2D in both East Asian and European populations exhibited strongly correlated effect sizes. Previously undescribed associations include signals in or near GDAP1, PTF1A, SIX3, ALDH2, a microRNA cluster, and genes that affect the differentiation of muscle and adipose cells3. At another locus, expression quantitative trait loci at two overlapping T2D signals affect two genes-NKX6-3 and ANK1-in different tissues4-6. Association studies in diverse populations identify additional loci and elucidate disease-associated genes, biology, and pathways.

Timing of antibiotics in septic patients: a prospective cohort study

OBJECTIVES

To evaluate the effect of timing and appropriateness of antibiotics administration on mortality in patients diagnosed with sepsis according to the Sepsis-3 definition.

METHODS

This prospective cohort study was conducted in patients diagnosed with sepsis according to the Sepsis-3 definition at the emergency department of Korea University Ansan Hospital from January 2016 to January 2019. The time to antibiotics was defined as the time in hours from emergency department arrival to the first antibiotic administration. Cox proportional hazards regression analysis was used to estimate the association between time to antibiotics and 7-, 14- and 28-day mortality.

RESULTS

Of 482 patients enrolled onto this study, 203 (42.1%) of 482 and 312 (64.7%) of 482 were diagnosed with septic shock and high-grade infection respectively. The median time to receipt of antibiotic therapy was 115 minutes. Antibiotics were administered within 3 and 6 hours in 340 (70.4%) of 482 and 450 (93.2%) of 482 patients respectively. Initial appropriate empirical antibiotics were administered in 375 (77.8%) of 482 patients. The time to and appropriateness of the initial antibiotics were not associated with 7-, 14- and 28-day mortality in multivariate analysis. The Sequential Organ Failure Assessment (SOFA) score (adjusted hazard ratio (aHR) 1.229, 95% confidence interval (CI) 1.093-1.381, p 0.001) and initial lactate levels (aHR 1.128, 95% CI 1.034-1.230, p 0.007), Charlson comorbidity index (aHR 1.115, 95% CI 1.027-1.210, p 0.014), 2-hour lactate level (aHR 1.115, 95% CI 1.027-1.210, p 0.009) and SOFA score (aHR 1.077, 95% CI 1.013-1.144, p 0.018) affected 7-, 14- and 28-day mortality respectively. Subgroup analysis with septic shock, bacteraemia and high-grade infection did not affect mortality rates.

CONCLUSIONS

Time to receipt of antibiotics may not affect the prognosis of patients with sepsis if a rapid and well-trained resuscitation is combined with appropriate antibiotic administration within a reasonable time.

Programmable Nuclease-Based Integration into Novel Extragenic Genomic Safe Harbor Identified from Korean Population-Based CNV Analysis

Here, we found two genomic safe harbor (GSH) candidates from chromosomes 3 and 8, based on large-scale population-based cohort data from 4,694 Koreans by CNV analysis. Furthermore, estimated genotype of these CNVRs was validated by quantitative real-time PCR, and epidemiological data examined no significant genetic association between diseases or traits and two CNVRs. After screening the GSH candidates by in silico approaches, we designed TALEN pairs to integrate EGFP expression cassette into human cell lines in order to confirm the functionality of GSH candidates in an in vitro setting. As a result, transgene insertion into one of the two loci using TALEN showed robust transgene expression comparable to that with an AAVS1 site without significantly perturbing neighboring genes. Changing the promoter or cell type did not noticeably disturb this trend. Thus, we could validate two CNVRs as a site for effective and safe transgene insertion in human cells.

Exome chip-driven association study of lipidemia in >14,000 Koreans and evaluation of genetic effect on identified variants between different ethnic groups

Lipid levels in blood are widely used to diagnose and monitor chronic diseases. It is essential to identify the genetic traits involved in lipid metabolism for understanding chronic diseases. However, the influence of genetic traits varies depending on race, sex, age, and ethnicity. Therefore, research focusing on populations of individual countries is required, and the results can be used as a basis for comparison of results of other studies at the cross-racial and cross-country levels. In the present study, we selected lipid-related variants and evaluated their effects on lipid-related diseases in more than 14,000 subjects of three cohorts using the Illumina Human Exome Beadchip. A genome-wide association study was conducted using EPACTs after adjusting for age, sex, and recruitment area. A genome-wide significance cutoff was defined as p < 5E-08 in all the three cohorts. Sixteen variants represented the lipid traits and were classified as vulnerable to borderline hypertriglyceridemia, hyper-LDL-cholesterolemia, or hypo-HDL-cholesterolemia. Moreover, we compared the genetic effects of the 16 variants between ethnic groups and identified the missense variants in apolipoprotein A-V, cholesterol ester transfer protein, and apolipoprotein E as Asian-specific. Our study provides candidate genes as markers for chronic diseases through the evaluation of genetic effects.

Abstract PD1-07: Withdrawn

This abstract was withdrawn by the authors.

Citation Format: Howell SJ, Waters S, Twelves C, Joffe J, Moon S, Bale C, Venkitaraman R, Bezecny P, Casbard A, Wilhelm-Benartzi C, Carucci M, Butler R, Alchami F, Jones R. Withdrawn [abstract]. In: Proceedings of the 2018 San Antonio Breast Cancer Symposium; 2018 Dec 4-8; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2019;79(4 Suppl):Abstract nr PD1-07.

The Korea Biobank Array: Design and Identification of Coding Variants Associated with Blood Biochemical Traits

We introduce the design and implementation of a new array, the Korea Biobank Array (referred to as KoreanChip), optimized for the Korean population and demonstrate findings from GWAS of blood biochemical traits. KoreanChip comprised >833,000 markers including >247,000 rare-frequency or functional variants estimated from >2,500 sequencing data in Koreans. Of the 833 K markers, 208 K functional markers were directly genotyped. Particularly, >89 K markers were presented in East Asians. KoreanChip achieved higher imputation performance owing to the excellent genomic coverage of 95.38% for common and 73.65% for low-frequency variants. From GWAS (Genome-wide association study) using 6,949 individuals, 28 associations were successfully recapitulated. Moreover, 9 missense variants were newly identified, of which we identified new associations between a common population-specific missense variant, rs671 (p.Glu457Lys) of ALDH2, and two traits including aspartate aminotransferase (P = 5.20 × 10⁻¹³) and alanine aminotransferase (P = 4.98 × 10⁻⁸). Furthermore, two novel missense variants of GPT with rare frequency in East Asians but extreme rarity in other populations were associated with alanine aminotransferase (rs200088103; p.Arg133Trp, P = 2.02 × 10⁻⁹ and rs748547625; p.Arg143Cys, P = 1.41 × 10⁻⁶). These variants were successfully replicated in 6,000 individuals (P = 5.30 × 10⁻⁸ and P = 1.24 × 10⁻⁶). GWAS results suggest the promising utility of KoreanChip with a substantial number of damaging variants to identify new population-specific disease-associated rare/functional variants.

Multiple genotype-phenotype association study reveals intronic variant pair on SIDT2 associated with metabolic syndrome in a Korean population

Background

Metabolic syndrome is a risk factor for type 2 diabetes and cardiovascular disease. We identified common genetic variants that alter the risk for metabolic syndrome in the Korean population. To isolate these variants, we conducted a multiple-genotype and multiple-phenotype genome-wide association analysis using the family-based quasi-likelihood score (MFQLS) test. For this analysis, we used 7211 and 2838 genotyped study subjects for discovery and replication, respectively. We also performed a multiple-genotype and multiple-phenotype analysis of a gene-based single-nucleotide polymorphism (SNP) set.

Results

We found an association between metabolic syndrome and an intronic SNP pair, rs7107152 and rs1242229, in SIDT2 gene at 11q23.3. Both SNPs correlate with the expression of SIDT2 and TAGLN, whose products promote insulin secretion and lipid metabolism, respectively. This SNP pair showed statistical significance at the replication stage.

Conclusions

Our findings provide insight into an underlying mechanism that contributes to metabolic syndrome.

Electronic supplementary material

The online version of this article (10.1186/s40246-018-0180-4) contains supplementary material, which is available to authorized users.

Nonsynonymous Variants in PAX4 and GLP1R Are Associated With Type 2 Diabetes in an East Asian Population

We investigated ethnicity-specific exonic variants of type 2 diabetes (T2D) and its related clinical phenotypes in an East Asian population. We performed whole-exome sequencing in 917 T2D case and control subjects, and the findings were validated by exome array genotyping in 3,026 participants. In silico replication was conducted for seven nonsynonymous variants in an additional 13,122 participants. Single-variant and gene-based association tests for T2D were analyzed. A total of 728,838 variants were identified by whole-exome sequencing. Among nonsynonymous variants, PAX4 Arg192His increased risk of T2D and GLP1R Arg131Gln decreased risk of T2D in genome-wide significance (odds ratio [OR] 1.48, P = 4.47 × 10^-16 and OR 0.84, P = 3.55 × 10^-8, respectively). Another variant at PAX4 192 codon Arg192Ser was nominally associated with T2D (OR 1.62, P = 5.18 × 10^-4). In T2D patients, PAX4 Arg192His was associated with earlier age at diagnosis, and GLP1R Arg131Gln was associated with decreased risk of cardiovascular disease. In control subjects without diabetes, the PAX4 Arg192His was associated with higher fasting glucose and GLP1R Arg131Gln was associated with lower fasting glucose and HbA_1c level. Gene-based analysis revealed that SLC30A8 was most significantly associated with decreased risk of T2D (P = 1.0 × 10^-4). In summary, we have identified nonsynonymous variants associated with risk of T2D and related phenotypes in Koreans.

Association analyses of East Asian individuals and trans-ancestry analyses with European individuals reveal new loci associated with cholesterol and triglyceride levels

Large-scale meta-analyses of genome-wide association studies (GWAS) have identified >175 loci associated with fasting cholesterol levels, including total cholesterol (TC), high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), and triglycerides (TG). With differences in linkage disequilibrium (LD) structure and allele frequencies between ancestry groups, studies in additional large samples may detect new associations. We conducted staged GWAS meta-analyses in up to 69,414 East Asian individuals from 24 studies with participants from Japan, the Philippines, Korea, China, Singapore, and Taiwan. These meta-analyses identified (P < 5 × 10-8) three novel loci associated with HDL-C near CD163-APOBEC1 (P = 7.4 × 10-9), NCOA2 (P = 1.6 × 10-8), and NID2-PTGDR (P = 4.2 × 10-8), and one novel locus associated with TG near WDR11-FGFR2 (P = 2.7 × 10-10). Conditional analyses identified a second signal near CD163-APOBEC1. We then combined results from the East Asian meta-analysis with association results from up to 187,365 European individuals from the Global Lipids Genetics Consortium in a trans-ancestry meta-analysis. This analysis identified (log10Bayes Factor ≥6.1) eight additional novel lipid loci. Among the twelve total loci identified, the index variants at eight loci have demonstrated at least nominal significance with other metabolic traits in prior studies, and two loci exhibited coincident eQTLs (P < 1 × 10-5) in subcutaneous adipose tissue for BPTF and PDGFC. Taken together, these analyses identified multiple novel lipid loci, providing new potential therapeutic targets.

Impact of common genetic determinants of Hemoglobin A1c on type 2 diabetes risk and diagnosis in ancestrally diverse populations: A transethnic genome-wide meta-analysis

BACKGROUND

Glycated hemoglobin (HbA1c) is used to diagnose type 2 diabetes (T2D) and assess glycemic control in patients with diabetes. Previous genome-wide association studies (GWAS) have identified 18 HbA1c-associated genetic variants. These variants proved to be classifiable by their likely biological action as erythrocytic (also associated with erythrocyte traits) or glycemic (associated with other glucose-related traits). In this study, we tested the hypotheses that, in a very large scale GWAS, we would identify more genetic variants associated with HbA1c and that HbA1c variants implicated in erythrocytic biology would affect the diagnostic accuracy of HbA1c. We therefore expanded the number of HbA1c-associated loci and tested the effect of genetic risk-scores comprised of erythrocytic or glycemic variants on incident diabetes prediction and on prevalent diabetes screening performance. Throughout this multiancestry study, we kept a focus on interancestry differences in HbA1c genetics performance that might influence race-ancestry differences in health outcomes.

METHODS & FINDINGS

Using genome-wide association meta-analyses in up to 159,940 individuals from 82 cohorts of European, African, East Asian, and South Asian ancestry, we identified 60 common genetic variants associated with HbA1c. We classified variants as implicated in glycemic, erythrocytic, or unclassified biology and tested whether additive genetic scores of erythrocytic variants (GS-E) or glycemic variants (GS-G) were associated with higher T2D incidence in multiethnic longitudinal cohorts (N = 33,241). Nineteen glycemic and 22 erythrocytic variants were associated with HbA1c at genome-wide significance. GS-G was associated with higher T2D risk (incidence OR = 1.05, 95% CI 1.04-1.06, per HbA1c-raising allele, p = 3 × 10-29); whereas GS-E was not (OR = 1.00, 95% CI 0.99-1.01, p = 0.60). In Europeans and Asians, erythrocytic variants in aggregate had only modest effects on the diagnostic accuracy of HbA1c. Yet, in African Americans, the X-linked G6PD G202A variant (T-allele frequency 11%) was associated with an absolute decrease in HbA1c of 0.81%-units (95% CI 0.66-0.96) per allele in hemizygous men, and 0.68%-units (95% CI 0.38-0.97) in homozygous women. The G6PD variant may cause approximately 2% (N = 0.65 million, 95% CI 0.55-0.74) of African American adults with T2D to remain undiagnosed when screened with HbA1c. Limitations include the smaller sample sizes for non-European ancestries and the inability to classify approximately one-third of the variants. Further studies in large multiethnic cohorts with HbA1c, glycemic, and erythrocytic traits are required to better determine the biological action of the unclassified variants.

CONCLUSIONS

As G6PD deficiency can be clinically silent until illness strikes, we recommend investigation of the possible benefits of screening for the G6PD genotype along with using HbA1c to diagnose T2D in populations of African ancestry or groups where G6PD deficiency is common. Screening with direct glucose measurements, or genetically-informed HbA1c diagnostic thresholds in people with G6PD deficiency, may be required to avoid missed or delayed diagnoses.

On the association analysis of CNV data: a fast and robust family-based association method

BACKGROUND

Copy number variation (CNV) is known to play an important role in the genetics of complex diseases and several methods have been proposed to detect association of CNV with phenotypes of interest. Statistical methods for CNV association analysis can be categorized into two different strategies. First, the copy number is estimated by maximum likelihood and association of the expected copy number with the phenotype is tested. Second, the observed probe intensity measurements can be directly used to detect association of CNV with the phenotypes of interest.

RESULTS

For each strategy we provide a statistic that can be applied to extended families. The computational efficiency of the proposed methods enables genome-wide association analysis and we show with simulation studies that the proposed methods outperform other existing approaches. In particular, we found that the first strategy is always more efficient than the second strategy no matter whether copy numbers for each individual are well identified or not. With the proposed methods, we performed genome-wide CNV association analyses of hematological trait, hematocrit, on 521 Korean family samples.

CONCLUSIONS

We found that statistical analysis with the expected copy number is more powerful than the statistic with the probe intensity measurements regardless of the accuracy of the estimation of copy numbers.

Genome-Wide Association Study Meta-Analysis of Long-Term Average Blood Pressure in East Asians

BACKGROUND

Genome-wide single marker and gene-based meta-analyses of long-term average (LTA) blood pressure (BP) phenotypes may reveal novel findings for BP.

METHODS AND RESULTS

We conducted genome-wide analysis among 18 422 East Asian participants (stage 1) followed by replication study of ≤46 629 participants of European ancestry (stage 2). Significant single-nucleotide polymorphisms and genes were determined by a P<5.0×10-8 and 2.5×10-6, respectively, in joint analyses of stage-1 and stage-2 data. We identified 1 novel ARL3 variant, rs4919669 at 10q24.32, influencing LTA systolic BP (stage-1 P=5.03×10-8, stage-2 P=8.64×10-3, joint P=2.63×10-8) and mean arterial pressure (stage-1 P=3.59×10-9, stage-2 P=2.35×10-2, joint P=2.64×10-8). Three previously reported BP loci (WBP1L, NT5C2, and ATP2B1) were also identified for all BP phenotypes. Gene-based analysis provided the first robust evidence for association of KCNJ11 with LTA systolic BP (stage-1 P=8.55×10-6, stage-2 P=1.62×10-5, joint P=3.28×10-9) and mean arterial pressure (stage-1 P=9.19×10-7, stage-2 P=9.69×10-5, joint P=2.15×10-9) phenotypes. Fourteen genes (TMEM180, ACTR1A, SUFU, ARL3, SFXN2, WBP1L, CYP17A1, C10orf32, C10orf32-ASMT, AS3MT, CNNM2, and NT5C2 at 10q24.32; ATP2B1 at 12q21.33; and NCR3LG1 at 11p15.1) implicated by previous genome-wide association study meta-analyses were also identified. Among the loci identified by the previous genome-wide association study meta-analysis of LTA BP, we transethnically replicated associations of the KCNK3 marker rs1275988 at 2p23.3 with LTA systolic BP and mean arterial pressure phenotypes (P=1.27×10-4 and 3.30×10-4, respectively).

CONCLUSIONS

We identified 1 novel variant and 1 novel gene and present the first direct evidence of relevance of the KCNK3 locus for LTA BP among East Asians.

Lithium suppression of tau induces brain iron accumulation and neurodegeneration

Lithium is a first-line therapy for bipolar affective disorder. However, various adverse effects, including a Parkinson-like hand tremor, often limit its use. The understanding of the neurobiological basis of these side effects is still very limited. Nigral iron elevation is also a feature of Parkinsonian degeneration that may be related to soluble tau reduction. We found that magnetic resonance imaging T₂ relaxation time changes in subjects commenced on lithium therapy were consistent with iron elevation. In mice, lithium treatment lowers brain tau levels and increases nigral and cortical iron elevation that is closely associated with neurodegeneration, cognitive loss and parkinsonian features. In neuronal cultures lithium attenuates iron efflux by lowering tau protein that traffics amyloid precursor protein to facilitate iron efflux. Thus, tau- and amyloid protein precursor-knockout mice were protected against lithium-induced iron elevation and neurotoxicity. These findings challenge the appropriateness of lithium as a potential treatment for disorders where brain iron is elevated (for example, Alzheimer's disease), and may explain lithium-associated motor symptoms in susceptible patients.

Identification of a Copy Number Variation on Chromosome 20q13.12 Associated with Osteoporotic Fractures in the Korean Population

Osteoporotic fractures (OFs) are critical hard outcomes of osteoporosis and are characterized by decreased bone strength induced by low bone density and microarchitectural deterioration in bone tissue. Most OFs cause acute pain, hospitalization, immobilization, and slow recovery in patients and are associated with increased mortality. A variety of genetic studies have suggested associations of genetic variants with the risk of OF. Genome-wide association studies have reported various single-nucleotide polymorphisms and copy number variations (CNVs) in European and Asian populations. To identify CNV regions associated with OF risk, we conducted a genome-wide CNV study in a Korean population. We performed logistic regression analyses in 1,537 Korean subjects (299 OF cases and 1,238 healthy controls) and identified a total of 8 CNV regions significantly associated with OF (p < 0.05). Then, one CNV region located on chromosome 20q13.12 was selected for experimental validation. The selected CNV region was experimentally validated by quantitative polymerase chain reaction. The CNV region of chromosome 20q13.12 is positioned upstream of a family of long non-coding RNAs, LINC01260. Our findings could provide new information on the genetic factors associated with the risk of OF.