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

Identification of Novel Genetic Loci Affecting Age at Onset of Parkinson's Disease: A Genome-wide Association Study

BACKGROUND

The age at onset (AAO) of Parkinson's disease (PD) varies widely among individuals and significantly influences disease progression and prognosis. However, few genome-wide association studies (GWASs) have investigated genetic variants determining AAO, particularly in East Asian populations.

OBJECTIVES

To identify single-nucleotide polymorphisms (SNPs) affecting AAO of PD in Korean patients.

METHODS

We conducted a GWAS on AAO of PD in 1048 Korean patients using sex-adjusted linear regression models. Additionally, we conducted downstream analyses of our primary GWAS results.

RESULTS

rs2134545 demonstrated genome-wide significance (β = -2.459; standard error [SE] = 0.851; P = 1.898 × 10-8) and is an intergenic SNP near the ALCAM gene associated with an average AAO reduction of 3.47 years. Additionally, rs4366309 (LYST; MIR1537) demonstrated suggestive significance (β = 2.949; SE = 1.072; P = 8.68 × 10-8) and was associated with an average delay of 3.05 years. The polygenic risk score based on known PD risk loci also affected the AAO for European and Korean PD risk loci, respectively (β = -0.149; P < 0.001 and β = -0.096; P = 0.002). However, the proportion of variance was small (r2 = 0.022 and 0.009, respectively).

CONCLUSION

We identified a novel SNP associated with the AAO of PD near the ALCAM gene, distinct from previously reported PD risk loci. These findings need further functional validation; however, they suggest unique genetic pathways influencing the AAO of PD and highlight the need for further research in diverse populations. © 2024 The Author(s). Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Enhancing Familial Hypercholesterolemia Detection in South Korea: A Targeted Screening Approach Integrating National Program and Genetic Cascade Screening

BACKGROUND AND OBJECTIVES

Familial hypercholesterolemia (FH) increases the risk of premature cardiovascular disease through disrupted low-density lipoprotein cholesterol (LDL-C) metabolism. Although FH is a severe condition, it remains widely underdiagnosed, which can be attributed to barriers in genetic testing and a lack of awareness. This study aims to propose and evaluate a targeted screening program for FH in South Korea by integrating the General Health Screening Program (GHSP) with cascade genetic screening.

METHODS

The study included individuals with LDL-C levels ≥190 mg/dL identified during the 2021 GHSP (primary participants). Data on demographics, lifestyle, medical history, and family history were collected through questionnaires. Targeted next-generation sequencing was used to identify pathogenic mutations in the PCSK9, APOB, LDLRAP1, and LDLR genes associated with FH. Pathogenic mutations found in primary participants were confirmed in their relatives (secondary participants) using Sanger sequencing. Participant characteristics were analyzed based on the presence of pathogenic mutations.

RESULTS

Among 83 individuals with severe hypercholesterolemia identified through the GHSP, 7 primary participants (8.4%) carried pathogenic mutations in the LDLR and PCSK9 genes. In secondary participants, pathogenic mutations were identified in 61.1% of the relatives of 4 patients with pathogenic mutations. The prevalence of pathogenic mutations was significantly higher in primary participants compared to secondary participants.

CONCLUSIONS

Integrating community resources with FH screening can enhance the early detection and treatment of FH. By utilizing GHSP data and adding genetic screening, the proposed model provides a strategy to reduce the cardiovascular risks associated with FH, supporting its wider adoption at the national level.

Effects of polygenic risk score and sodium and potassium intake on hypertension in Asians: A nationwide prospective cohort study

Genetic factors, lifestyle, and diet have been shown to play important roles in the development of hypertension. Increased salt intake is an important risk factor for hypertension. However, research on the involvement of genetic factors in the relationship between salt intake and hypertension in Asians is lacking. We aimed to investigate the risk of hypertension in relation to sodium and potassium intake and the effects of genetic factors on their interactions. We used Korean Genome and Epidemiology Study data and calculated the polygenic risk score (PRS) for the effect of systolic and diastolic blood pressure (SBP and DBP). We also conducted multivariable logistic modeling to evaluate associations among incident hypertension, PRSSBP, PRSDBP, and sodium and potassium intake. In total, 41,351 subjects were included in the test set. The top 10% PRSSBP group was the youngest of the three groups (bottom 10%, middle, top 10%), had the highest proportion of women, and had the highest body mass index, baseline BP, red meat intake, and alcohol consumption. The multivariable logistic regression model revealed the risk of hypertension was significantly associated with higher PRSSBP, higher sodium intake, and lower potassium intake. There was significant interaction between sodium intake and PRSSBP for incident hypertension especially in sodium intake ≥2.0 g/day and PRSSBP top 10% group (OR 1.27 (1.07-1.51), P = 0.007). Among patients at a high risk of incident hypertension due to sodium intake, lifestyle modifications and sodium restriction were especially important to prevent hypertension.

Gut microbial and human genetic signatures of inflammatory bowel disease increase risk of comorbid mental disorders

The high prevalence of comorbid mental disorders (CMDs) in patients with inflammatory bowel disease (IBD) is well-documented. This study delves into the intricate CMD-IBD relationship through comprehensive analyses using human variants, gut microbiome, and anxiety/depression estimates from a cohort of 507 IBD patients and 75 controls. Notably, patients with IBD, especially those with CMD, exhibited lower diversity than controls. We identified 106 differentially abundant taxa (DATs) in IBD patients compared to controls and 21 DATs distinguishing CMD-affected from CMD-free IBD patients. Microbial IBD-risk scores, reflecting an individual's microbial burden for IBD, revealed a significant enrichment of IBD-risk signatures in CMD-affected patients compared to CMD-free patients. Additionally, there was an IBD-risk variant potentially regulating the abundance of an IBD/CMD-associated DAT, suggesting an interplay between IBD-risk variants and dysbiosis in CMD. Our investigation underscores the pivotal role of IBD-associated gut dysbiosis in predisposing IBD patients to CMD, partially through genetic variant-mediated mechanisms.

Cardiovascular risk according to genetic predisposition to gout, lifestyle and metabolic health across prospective European and Korean cohorts

OBJECTIVE

Recent studies have reported that gout is associated with a risk of cardiovascular disease (CVD) later in life. However, the predictive value of genetic predisposition to gout combined with lifestyle habits for CVD risk remains unclear. This study aimed to examine the association between genetic predisposition to gout and lifestyle habits and the risk of developing CVD in two diverse prospective cohorts from different ancestries.

METHODS

A total of 224 689 participants of European descent from the UK Biobank and 50 364 participants of East Asian descent from the Korean Genome and Epidemiology Study were included. The genetic risk for gout was assessed using a polygenic risk score (PRS) derived from a meta-genome-wide association study (n=444 533). The incident CVD risk was evaluated according to genetic risk, lifestyle and metabolic syndrome (MetS).

RESULTS

Individuals at high genetic risk for gout had a higher risk of incident CVD than those with low genetic risk across ancestry. Notably, a reduction in CVD risk by up to 62% (HR 0.38; 95% CI 0.31 to 0.46; p <0.001) was observed in individuals at both low and high genetic risk for gout when they maintained ideal MetS and favourable lifestyle habits.

CONCLUSIONS

Our findings indicate that a higher genetic risk of gout is significantly associated with an increased risk of CVD. Moreover, adherence to a favourable lifestyle can significantly reduce CVD risk, particularly in individuals with high genetic risk. These results underscore the potential of PRS-based risk assessment to improve clinical outcomes through tailored preventative strategies.

A combination of red and processed meat intake and polygenic risk score influences the incidence of hyperuricemia in middle-aged Korean adults

BACKGROUND/OBJECTIVES

The high consumption of purine-rich meat is associated with hyperuricemia. However, there is limited evidence linking the consumption of red and processed meat to the genetic risk of hyperuricemia. We investigated the relationship between various combinations of red and processed meat consumption and the polygenic risk scores (PRSs) and the incidence of hyperuricemia in middle-aged Koreans.

SUBJECTS/METHODS

We analyzed the data from 44,053 participants aged ≥40 years sourced from the Health Examinees (HEXA) cohort of the Korean Genome and Epidemiology Study (KoGES). Information regarding red and processed meat intake was obtained using a semiquantitative food frequency questionnaire (SQ-FFQ). We identified 69 independent single-nucleotide polymorphisms (SNPs) at uric acid-related loci using genome-wide association studies (GWASs) and clumping analyses. The individual PRS, which is the weighted sum of the effect size of each allele at the SNP, was calculated. We used multivariable Cox proportional hazards models adjusted for covariates to determine the relationship between red and processed meat intake and the PRS in the incidence of hyperuricemia.

RESULTS

During an average follow-up period of 5 years, 2,556 patients with hyperuricemia were identified. For both men and women, the group with the highest red and processed meat intake and the highest PRS was positively associated with the development of hyperuricemia when compared with the group with the lowest red and processed meat intake and the lowest PRS (hazard ratio [HR], 2.72; 95% confidence interval [CI], 2.10-3.53; P < 0.0001; HR, 3.28; 95% CI, 2.45-4.40; P < 0.0001).

CONCLUSION

Individuals at a high genetic risk for uric acid levels should moderate their consumption of red and processed meat to prevent hyperuricemia.

Assessment of the Impact of Preanalytical DNA Integrity on the Genome Data Quality

Many molecular approaches have been employed for the quality control (QC) of biobanked DNA samples. Since 2003, the National Biobank of Korea (NBK) has provided various studies with over half a million quality-controlled genomic DNA samples using conventional agarose gel electrophoresis and spectrophotometry. We assessed the postanalytical genomic data quality of DNA samples (n = 41) with a different range of the DNA quality index such as genomic quality number (GQN) for developing an evidence-based best practice for DNA quality criteria. We examined the quality indices of three different platforms, including single nucleotide polymorphism arrays, methylation arrays, and next-generation sequencing, using the same DNA samples (n = 41) of different quality, ranging from 4.0 to 10.0 values of the GQN. Our data analysis revealed that higher GQN value and/or double-stranded DNA concentration resulted in higher quality genomic data. In addition, all the analyzed DNA samples successfully generated good-quality genomic data. This study provides a guide for the QC of biobanked DNA samples for genomic analysis platforms.

NPC1L1 rs217434 A > G as a Novel Single Nucleotide Polymorphism Related to Dyslipidemia in a Korean Population

A relationship between cholesterol levels and Niemann-Pick C1-Like 1 (NPC1L1) polymorphisms in diverse populations was found in previous studies. However, relevant research on this association in the Korean population is relatively scarce. Therefore, the current study sought to examine the correlation between the NPC1L1 rs217434 A > G polymorphism and clinical as well as biochemical variables pertaining to dyslipidemia in the Korean population. This cross-sectional single-center study included 1404 Korean subjects aged 20-86 years, grouped based on dyslipidemia presence (normal and dyslipidemia) and genotype (AA or AG). After adjusting for sex and age, it was discovered that the dyslipidemia group's BMI, diastolic blood pressure, glucose-related indicators, lipid profile, high-sensitivity C-reactive protein (hs-CRP), and parameters of oxidative stress were considerably different from the normal group's values. When grouped according to genotype, individuals in the AG group exhibited greater total cholesterol, low-density lipoprotein cholesterol, hs-CRP, and 8-epi-prostaglandin F2α in comparison to those in the AA group. Moreover, individuals with dyslipidemia and the AG genotype exhibited unfavorable outcomes for lipid profiles, markers related to glucose and inflammation, and markers of oxidative stress. This study provided evidence for a relationship between the NPC1L1 rs217434 A > G genotype and dyslipidemia in the Korean population, which highlights the potential of the NPC1L1 rs217434 A > G genotype as an early predictor of dyslipidemia.

Fine-mapping genomic loci refines bipolar disorder risk genes

Bipolar disorder (BD) is a heritable mental illness with complex etiology. While the largest published genome-wide association study identified 64 BD risk loci, the causal SNPs and genes within these loci remain unknown. We applied a suite of statistical and functional fine-mapping methods to these loci, and prioritized 17 likely causal SNPs for BD. We mapped these SNPs to genes, and investigated their likely functional consequences by integrating variant annotations, brain cell-type epigenomic annotations, brain quantitative trait loci, and results from rare variant exome sequencing in BD. Convergent lines of evidence supported the roles of genes involved in neurotransmission and neurodevelopment including SCN2A, TRANK1, DCLK3, INSYN2B, SYNE1, THSD7A, CACNA1B, TUBBP5, PLCB3, PRDX5, KCNK4, CRTC3, AP001453.3, TRPT1, FKBP2, DNAJC4, RASGRP1, FURIN, FES, DPH1, GSDMB, MED24 and THRA in BD. These represent promising candidates for functional experiments to understand biological mechanisms and therapeutic potential. Additionally, we demonstrated that fine-mapping effect sizes can improve performance of BD polygenic risk scores across diverse populations, and present a high-throughput fine-mapping pipeline (https://github.com/mkoromina/SAFFARI).

Predicting structure-targeted food bioactive compounds for middle-aged and elderly Asians with myocardial infarction: insights from genetic variations and bioinformatics-integrated deep learning analysis

Myocardial infarction (MI) is a significant global health issue. Despite the advances in genome-wide association studies, a complete genetic and molecular understanding of MI is elusive and needs to be fully explored. This study aimed to elucidate the genetic framework of MI and explore the potential health benefits of natural compounds (NCs). The genetic architecture of MI was explored using data from the Korean Genome and Epidemiology Study. We pinpointed crucial protein-coding genes related to MI by multi-marker analysis of genomic annotation for gene-based analysis. The bioinformatics-integrated deep neural analysis of NCs (BioDeepNat), a novel disease discovery application, was utilized to assess the influence of NCs on MI-related target proteins and validated with molecular docking analysis. The BioDeepNat application revealed significant NCs on MI-related target proteins, such as E-resveratrol, epicatechin 3-gallate, and kaempferol. The E3 region of RNF213 protein with a point mutation (Arg4810Lys) had different binding energies with NCs, such as ursolic acid and olean-12-en-28-oic acid, compared to the wild type. However, ginsenosides, eleutheroside, oleanolic acid, and hederagenin showed similar binding energies to wild and mutated types of RNF213 protein. The predicted NCs were primarily sourced from foods such as common grapes and teas. Aromatic hydrocarbons are frequently observed as the prevalent functional groups with high binding affinity for NCs in a molecular docking analysis. In conclusion, the proteins encoded by these genes identified by gene-based analysis interacted with several NCs with health promotion found in day-to-day foods, particularly E-resveratrol and kaempferol. This understanding offers promising directions for precision nutrition strategies in MI.

Genome-wide association study and polygenic risk score analysis for schizophrenia in a Korean population

Although large-scale genome-wide association studies (GWASs) have revealed the genetic architecture of schizophrenia, these studies have mainly focused on populations of European ancestry. This study aimed to identify common genetic variants associated with schizophrenia in the Korean population and evaluate the performance of polygenic risk scores (PRSs) derived from large-scale GWASs across ancestries. In the Korean psychiatric GWAS project (KPGP), seven academic institutes and their affiliated hospitals across South Korea recruited a cohort of 1670 patients with DSM-IV-defined schizophrenia and 2271 healthy controls. A total of 6690,822 SNPs were tested for association with schizophrenia. We identified one previously unreported genome-wide significant locus rs2423464 (P = 2.83 × 10-11; odds ratio = 1.65; 95 % confidence interval = 1.43-1.91, minor allele frequency = 0.126). This variant was also associated with increased lysosomal-associated membrane protein family member 5 (LAMP5) gene expression. The PRS derived from the meta-analysis results of East Asian and European GWASs explained a larger proportion of the phenotypic variance in the Korean schizophrenia sample than the PRS of an East Asian or European GWAS. (R2 = 0.073 for meta-analysis; 0.028 for East Asian GWAS; 0.037 for European GWAS). GWASs involving diverse ethnic groups will expand our understanding of the genetic architecture of schizophrenia.

Prediction of metabolic syndrome using machine learning approaches based on genetic and nutritional factors: a 14-year prospective-based cohort study

INTRODUCTION

Metabolic syndrome is a chronic disease associated with multiple comorbidities. Over the last few years, machine learning techniques have been used to predict metabolic syndrome. However, studies incorporating demographic, clinical, laboratory, dietary, and genetic factors to predict the incidence of metabolic syndrome in Koreans are limited. In the present study, we propose a genome-wide polygenic risk score for the prediction of metabolic syndrome, along with other factors, to improve the prediction accuracy of metabolic syndrome.

METHODS

We developed 7 machine learning-based models and used Cox multivariable regression, deep neural network (DNN), support vector machine (SVM), stochastic gradient descent (SGD), random forest (RAF), Naïve Bayes (NBA) classifier, and AdaBoost (ADB) to predict the incidence of metabolic syndrome at year 14 using the dataset from the Korean Genome and Epidemiology Study (KoGES) Ansan and Ansung.

RESULTS

Of the 5440 patients, 2,120 were considered to have new-onset metabolic syndrome. The AUC values of model, which included sex, age, alcohol intake, energy intake, marital status, education status, income status, smoking status, dried laver intake, and genome-wide polygenic risk score (gPRS) Z-score based on 344,447 SNPs (p-value < 1.0), were the highest for RAF (0.994 [95% CI 0.985, 1.000]) and ADB (0.994 [95% CI 0.986, 1.000]).

CONCLUSIONS

Incorporating both gPRS and demographic, clinical, laboratory, and seaweed data led to enhanced metabolic syndrome risk prediction by capturing the distinct etiologies of metabolic syndrome development. The RAF- and ADB-based models predicted metabolic syndrome more accurately than the NBA-based model for the Korean population.

Big data in sarcoidosis

PURPOSE OF REVIEW

This review provides an overview of recent advancements in sarcoidosis research, focusing on collaborative networks, phenotype characterization, and molecular studies. It highlights the importance of collaborative efforts, phenotype characterization, and the integration of multilevel molecular data for advancing sarcoidosis research and paving the way toward personalized medicine.

RECENT FINDINGS

Sarcoidosis exhibits heterogeneous clinical manifestations influenced by various factors. Efforts to define sarcoidosis endophenotypes show promise, while technological advancements enable extensive molecular data generation. Collaborative networks and biobanks facilitate large-scale studies, enhancing biomarker discovery and therapeutic protocols.

SUMMARY

Sarcoidosis presents a complex challenge due to its unknown cause and heterogeneous clinical manifestations. Collaborative networks, comprehensive phenotype delineation, and the utilization of cutting-edge technologies are essential for advancing our understanding of sarcoidosis biology and developing personalized medicine approaches. Leveraging large-scale epidemiological resources and biobanks and integrating multilevel molecular data offer promising avenues for unraveling the disease's heterogeneity and improving patient outcomes.

Genome-Wide Polygenic Risk Score Predicts Incident Type 2 Diabetes in Women With History of Gestational Diabetes

OBJECTIVE

Women with a history of gestational diabetes mellitus (GDM) are at increased risk of developing type 2 diabetes (T2D). It remains unclear whether genetic information improves prediction of incident T2D in these women.

RESEARCH DESIGN AND METHODS

Using five independent cohorts representing four different ancestries (n = 1,895), we investigated whether a genome-wide T2D polygenic risk score (PRS) is associated with increased risk of incident T2D. We also calculated the area under the receiver operating characteristics curve (AUROC) and continuous net reclassification improvement (NRI) following the incorporation of T2D PRS into clinical risk models to assess the diagnostic utility.

RESULTS

Among 1,895 women with previous history of GDM, 363 (19.2%) developed T2D in a range of 2 to 30 years. T2D PRS was higher in those who developed T2D (-0.08 vs. 0.31, P = 2.3 × 10-11) and was associated with an increased risk of incident T2D (odds ratio 1.52 per 1-SD increase, 95% CI 1.05-2.21, P = 0.03). In a model that includes age, family history of diabetes, systolic blood pressure, and BMI, the incorporation of PRS led to an increase in AUROC for T2D from 0.71 to 0.74 and an intermediate improvement of NRI (0.32, 95% CI 0.15-0.49, P = 3.0 × 10-4). Although there was variation, a similar trend was observed across study cohorts.

CONCLUSIONS

In cohorts of GDM women with diverse ancestry, T2D PRS was significantly associated with future development of T2D. A significant but small improvement was observed in AUROC when T2D PRS was integrated into clinical risk models to predict incident T2D.

Genome-wide association study for metabolic syndrome reveals APOA5 single nucleotide polymorphisms with multilayered effects in Koreans

BACKGROUND AND PURPOSE

Genome-wide association studies (GWAS) of metabolic syndrome (MetS) have predominantly focused on non-Asian populations, with limited representation from East Asian cohorts. Moreover, previous GWAS analyses have primarily emphasized the significance of top single nucleotide polymorphisms (SNPs), poorly explaining other SNP signals in linkage disequilibrium. This study aimed to reveal the interaction between rs651821 and rs2266788, the principal variants of apolipoprotein A5 (APOA5), within the most significant loci identified through GWAS on MetS.

METHODS

GWAS on MetS and its components was conducted using the data from the Korean Genome and Epidemiology Study (KoGES) city cohort comprising 58,600 individuals with available biochemical, demographic, lifestyle factors, and the most significant APOA5 locus was analyzed further in depth.

RESULTS

According to GWAS of MetS and its diagnostic components, a significant association between the APOA5 SNPs rs651821/rs2266788 and MetS/triglycerides/high-density lipoprotein phenotypes was revealed. However, a conditional analysis employing rs651821 unveiled a reversal in the odds ratio for rs2266788. Therefore, rs651821 and rs2266788 emerged as independent and opposing signals in the extended GWAS analysis, i.e., the multilayered effects. Further gene-environment interaction analyses regarding lifestyle factors such as smoking, alcohol consumption, and physical activity underscored these multilayered effects.

CONCLUSION

This study unveils the intricate interplay between rs651821 and rs2266788 derived from MetS GWAS. Removing the influence of lead SNP reveals an independent protective signal associated with rs2266788, suggesting a multilayered effect between these SNPs. These findings underline the need for novel perspectives in future MetS GWAS.

Forensic height estimation using polygenic score in Korean population

Height is known to be a classically heritable trait controlled by complex polygenic factors. Numerous height-associated genetic variants across the genome have been identified so far. It is also a representative of externally visible characteristics (EVC) for predicting appearance in forensic science. When biological evidence at a crime scene is deficient in identifying an individual, the examination of forensic DNA phenotyping using some genetic variants could be considered. In this study, we aimed to predict 'height', a representative forensic phenotype, by using a small number of genetic variants when short tandem repeat (STR) analysis is hard with insufficient biological samples. Our results not only replicated previous genetic signals but also indicated an upward trend in polygenic score (PGS) with increasing height in the validation and replication stages for both genders. These results demonstrate that the established SNP sets in this study could be used for height estimation in the Korean population. Specifically, since the PGS model constructed in this study targets only a small number of SNPs, it contributes to enabling forensic DNA phenotyping even at crime scenes with a minimal amount of biological evidence. To the best of our knowledge, this was the first study to evaluate a PGS model for height estimation in the Korean population using GWAS signals. Our study offers insight into the polygenic effect of height in East Asians, incorporating genetic variants from non-Asian populations.

A comprehensive multi-task deep learning approach for predicting metabolic syndrome with genetic, nutritional, and clinical data

Metabolic syndrome (MetS) is a complex disorder characterized by a cluster of metabolic abnormalities, including abdominal obesity, hypertension, elevated triglycerides, reduced high-density lipoprotein cholesterol, and impaired glucose tolerance. It poses a significant public health concern, as individuals with MetS are at an increased risk of developing cardiovascular diseases and type 2 diabetes. Early and accurate identification of individuals at risk for MetS is essential. Various machine learning approaches have been employed to predict MetS, such as logistic regression, support vector machines, and several boosting techniques. However, these methods use MetS as a binary status and do not consider that MetS comprises five components. Therefore, a method that focuses on these characteristics of MetS is needed. In this study, we propose a multi-task deep learning model designed to predict MetS and its five components simultaneously. The benefit of multi-task learning is that it can manage multiple tasks with a single model, and learning related tasks may enhance the model's predictive performance. To assess the efficacy of our proposed method, we compared its performance with that of several single-task approaches, including logistic regression, support vector machine, CatBoost, LightGBM, XGBoost and one-dimensional convolutional neural network. For the construction of our multi-task deep learning model, we utilized data from the Korean Association Resource (KARE) project, which includes 352,228 single nucleotide polymorphisms (SNPs) from 7729 individuals. We also considered lifestyle, dietary, and socio-economic factors that affect chronic diseases, in addition to genomic data. By evaluating metrics such as accuracy, precision, F1-score, and the area under the receiver operating characteristic curve, we demonstrate that our multi-task learning model surpasses traditional single-task machine learning models in predicting MetS.

Genetic determinants of obesity in Korean populations: exploring genome-wide associations and polygenic risk scores

East Asian populations exhibit a genetic predisposition to obesity, yet comprehensive research on these traits is limited. We conducted a genome-wide association study (GWAS) with 93,673 Korean subjects to uncover novel genetic loci linked to obesity, examining metrics such as body mass index, waist circumference, body fat ratio, and abdominal fat ratio. Participants were categorized into non-obese, metabolically healthy obese (MHO), and metabolically unhealthy obese (MUO) groups. Using advanced computational methods, we developed a multifaceted polygenic risk scores (PRS) model to predict obesity. Our GWAS identified significant genetic effects with distinct sizes and directions within the MHO and MUO groups compared with the non-obese group. Gene-based and gene-set analyses, along with cluster analysis, revealed heterogeneous patterns of significant genes on chromosomes 3 (MUO group) and 11 (MHO group). In analyses targeting genetic predisposition differences based on metabolic health, odds ratios of high PRS compared with medium PRS showed significant differences between non-obese and MUO, and non-obese and MHO. Similar patterns were seen for low PRS compared with medium PRS. These findings were supported by the estimated genetic correlation (0.89 from bivariate GREML). Regional analyses highlighted significant local genetic correlations on chromosome 11, while single variant approaches suggested widespread pleiotropic effects, especially on chromosome 11. In conclusion, our study identifies specific genetic loci and risks associated with obesity in the Korean population, emphasizing the heterogeneous genetic factors contributing to MHO and MUO.

Epigenetic link between Agent Orange exposure and type 2 diabetes in Korean veterans

Conflicting findings have been reported regarding the association between Agent Orange (AO) exposure and type 2 diabetes. This study aimed to examine whether AO exposure is associated with the development of type 2 diabetes and to verify the causal relationship between AO exposure and type 2 diabetes by combining DNA methylation with DNA genotype analyses. An epigenome-wide association study and DNA genotype analyses of the blood of AO-exposed and AO-unexposed individuals with type 2 diabetes and that of healthy controls were performed. Methylation quantitative trait locus and Mendelian randomisation analyses were performed to evaluate the causal effect of AO-exposure-identified CpGs on type 2 diabetes. AO-exposed individuals with type 2 diabetes were associated with six hypermethylated CpG sites (cg20075319, cg21757266, cg05203217, cg20102280, cg26081717, and cg21878650) and one hypo-methylated CpG site (cg07553761). Methylation quantitative trait locus analysis showed the methylation levels of some CpG sites (cg20075319, cg20102280, and cg26081717) to be significantly different. Mendelian randomisation analysis showed that CpG sites that were differentially methylated in AO-exposed individuals were causally associated with type 2 diabetes; the reverse causal effect was not significant. These findings reflect the need for further epigenetic studies on the causal relationship between AO exposure and type 2 diabetes.

Identification of shared genetic risks underlying metabolic syndrome and its related traits in the Korean population

Introduction: Observational studies have demonstrated strong correlations between metabolic syndrome (MetS) and its related traits. To gain insight into the genetic architecture and molecular mechanism of MetS, we investigated the shared genetic basis of MetS and its related traits and further tested their causal relationships. Methods: Using summary statistics from genome-wide association analyses of about 72,000 subjects from the Korean Genome and Epidemiological Study (KoGES), we conducted genome-wide multi-trait analyses to quantify the overall genetic correlation and Mendelian randomization analyses to infer the causal relationships between traits of interest. Results: Genetic correlation analyses revealed a significant correlation of MetS with its related traits, such as obesity traits (body mass index and waist circumference), lipid traits (triglyceride and high-density lipoprotein cholesterol), glycemic traits (fasting plasma glucose and hemoglobin A1C), and blood pressure (systolic and diastolic). Mendelian randomization analyses further demonstrated that the MetS-related traits showing significant overall genetic correlation with MetS could be genetically determined risk factors for MetS. Discussion: Our study suggests a shared genetic basis of MetS and its related traits and provides novel insights into the biological mechanisms underlying these complex traits. Our findings further inform public health interventions by supporting the important role of the management of metabolic risk factors such as obesity, unhealthy lipid profiles, diabetes, and high blood pressure in the prevention of MetS.

Sex-specific survival gene mutations are discovered as clinical predictors of clear cell renal cell carcinoma

Although sex differences have been reported in patients with clear cell renal cell carcinoma (ccRCC), biological sex has not received clinical attention and genetic differences between sexes are poorly understood. This study aims to identify sex-specific gene mutations and explore their clinical significance in ccRCC. We used data from The Cancer Genome Atlas-Kidney Renal Clear Cell Carcinoma (TCGA-KIRC), The Renal Cell Cancer-European Union (RECA-EU) and Korean-KIRC. A total of 68 sex-related genes were selected from TCGA-KIRC through machine learning, and 23 sex-specific genes were identified through verification using the three databases. Survival differences according to sex were identified in nine genes (ACSS3, ALG13, ASXL3, BAP1, JADE3, KDM5C, KDM6A, NCOR1P1, and ZNF449). Female-specific survival differences were found in BAP1 in overall survival (OS) (TCGA-KIRC, p = 0.004; RECA-EU, p = 0.002; and Korean-KIRC, p = 0.003) and disease-free survival (DFS) (TCGA-KIRC, p = 0.001 and Korean-KIRC, p = 0.000004), and NCOR1P1 in DFS (TCGA-KIRC, p = 0.046 and RECA-EU, p = 0.00003). Male-specific survival differences were found in ASXL3 (OS, p = 0.017 in TCGA-KIRC; and OS, p = 0.005 in RECA-EU) and KDM5C (OS, p = 0.009 in RECA-EU; and DFS, p = 0.016 in Korean-KIRC). These results suggest that biological sex may be an important predictor and sex-specific tailored treatment may improve patient care in ccRCC.

Transcriptomic and epigenomic landscapes of muscle growth during the postnatal period of broilers

BACKGROUND

Broilers stand out as one of the fastest-growing livestock globally, making a substantial contribution to animal meat production. However, the molecular and epigenetic mechanisms underlying the rapid growth and development of broiler chickens are still unclear. This study aims to explore muscle development patterns and regulatory networks during the postnatal rapid growth phase of fast-growing broilers. We measured the growth performance of Cornish (CC) and White Plymouth Rock (RR) over a 42-d period. Pectoral muscle samples from both CC and RR were randomly collected at day 21 after hatching (D21) and D42 for RNA-seq and ATAC-seq library construction.

RESULTS

The consistent increase in body weight and pectoral muscle weight across both breeds was observed as they matured, with CC outpacing RR in terms of weight at each stage of development. Differential expression analysis identified 398 and 1,129 genes in the two dimensions of breeds and ages, respectively. A total of 75,149 ATAC-seq peaks were annotated in promoter, exon, intron and intergenic regions, with a higher number of peaks in the promoter and intronic regions. The age-biased genes and breed-biased genes of RNA-seq were combined with the ATAC-seq data for subsequent analysis. The results spotlighted the upregulation of ACTC1 and FDPS at D21, which were primarily associated with muscle structure development by gene cluster enrichment. Additionally, a noteworthy upregulation of MUSTN1, FOS and TGFB3 was spotted in broiler chickens at D42, which were involved in cell differentiation and muscle regeneration after injury, suggesting a regulatory role of muscle growth and repair.

CONCLUSIONS

This work provided a regulatory network of postnatal broiler chickens and revealed ACTC1 and MUSTN1 as the key responsible for muscle development and regeneration. Our findings highlight that rapid growth in broiler chickens triggers ongoing muscle damage and subsequent regeneration. These findings provide a foundation for future research to investigate the functional aspects of muscle development.

Interplay of Mediterranean-diet adherence, genetic factors, and metabolic dysfunction-associated steatotic liver disease risk in Korea

BACKGROUNDS

Metabolic dysfunction-associated steatotic liver disease (MASLD) has gained attention owing to its severe complications. This study aimed to explore the interaction between Mediterranean-diet (MD) adherence, genetic factors, and MASLD risk in a Korean population.

METHODS

In total, 33,133 individuals aged 40 years and older from the Korean Genome and Epidemiology Study (KoGES) were analyzed. Participants were assessed for MASLD based on criteria and MD adherence measured by the Korean version of the Mediterranean-Diet Adherence Screener (K-MEDAS). Individuals were categorized into two groups based on their MD adherence: high adherence (K-MEDAS > 6) and low adherence (K-MEDAS < 5). Single nucleotide polymorphism (SNP) genotypes were obtained using the Korea Biobank array. Logistic regression was used to examine the single-marker variants for genetic associations with MASLD prevalence.

RESULTS

Individuals were categorized into MASLD (10,018 [30.2%]) and non-MASLD (23,115 [69.8%]) groups. A significant interaction was observed between the rs780094 glucokinase regulatory protein (GCKR) gene and K-MEDAS on MASLD (p < 10 - 2 ). Of individuals with K-MEDAS > 6, those carrying the minor allele (C) of the GCKR gene rs780094 exhibited a lower risk of MASLD compared to those without the allele (odds ratio [OR] = 0.88 [0.85-0.91], p-value = 5.54e-13).

CONCLUSION

The study identified a significant interaction involving the rs780094 variant near the GCKR gene, with carriers of the minor allele exhibiting a lower MASLD risk among those adhering well to the MD. Dietary habits influence the MASLD risk associated with the rs780094 allele, emphasizing the need for personalized nutrition recommendations.

Exploring intra-diagnosis heterogeneity and inter-diagnosis commonality in genetic architectures of bipolar disorders: association of polygenic risks of major psychiatric illnesses and lifetime phenotype dimensions

BACKGROUND

Bipolar disorder (BD) shows heterogeneous illness presentation both cross-sectionally and longitudinally. This phenotypic heterogeneity might reflect underlying genetic heterogeneity. At the same time, overlapping characteristics between BD and other psychiatric illnesses are observed at clinical and biomarker levels, which implies a shared biological mechanism between them. Incorporating these two issues in a single study design, this study investigated whether phenotypically heterogeneous subtypes of BD have a distinct polygenic basis shared with other psychiatric illnesses.

METHODS

Six lifetime phenotype dimensions of BD identified in our previous study were used as target phenotypes. Associations between these phenotype dimensions and polygenic risk scores (PRSs) of major psychiatric illnesses from East Asian (EA) and other available populations were analyzed.

RESULTS

Each phenotype dimension showed a different association pattern with PRSs of mental illnesses. PRS for EA schizophrenia showed a significant negative association with the cyclicity dimension (p = 0.044) but a significant positive association with the psychotic/irritable mania dimension (p = 0.001). PRS of EA major depressive disorder demonstrated a significant negative association with the elation dimension (p = 0.003) but a significant positive association with the comorbidity dimension (p = 0.028).

CONCLUSION

This study demonstrates that well-defined phenotype dimensions of lifetime-basis in BD have distinct genetic risks shared with other major mental illnesses. This finding supports genetic heterogeneity in BD and suggests a pleiotropy among BD subtypes and other psychiatric disorders beyond BD. Further genomic analyses adopting deep phenotyping across mental illnesses in ancestrally diverse populations are warranted to clarify intra-diagnosis heterogeneity and inter-diagnoses commonality issues in psychiatry.

Discovery and Validation of Survival-Specific Genes in Papillary Renal Cell Carcinoma Using a Customized Next-Generation Sequencing Gene Panel

PURPOSE

Papillary renal cell carcinoma (PRCC), the second most common kidney cancer, is morphologically, genetically, and molecularly heterogeneous with diverse clinical manifestations. Genetic variations of PRCC and their association with survival are not yet well-understood. This study aimed to identify and validate survival-specific genes in PRCC and explore their clinical utility.

MATERIALS AND METHODS

Using machine learning, 293 patients from the Cancer Genome Atlas-Kidney Renal Papillary Cell Carcinoma (TCGA-KIRP) database were analyzed to derive genes associated with survival. To validate these genes, DNAs were extracted from the tissues of 60 Korean PRCC patients. Next generation sequencing was conducted using a customized PRCC gene panel of 202 genes, including 171 survival-specific genes. Kaplan-Meier and Log-rank tests were used for survival analysis. Fisher's exact test was performed to assess the clinical utility of variant genes.

RESULTS

A total of 40 survival-specific genes were identified in the TCGA-KIRP database through machine learning and statistical analysis. Of them, 10 (BAP1, BRAF, CFDP1, EGFR, ITM2B, JAK1, NODAL, PCSK2, SPATA13, and SYT5) were validated in the Korean-KIRP database. Among these survival gene signatures, three genes (BAP1, PCSK2, and SPATA13) showed survival specificity in both overall survival (OS) (p = 0.00004, p = 1.38 × 10-7, and p = 0.026, respectively) and disease-free survival (DFS) (p = 0.00002, p = 1.21 × 10-7, and p = 0.036, respectively). Notably, the PCSK2 mutation demonstrated survival specificity uniquely in both the TCGA-KIRP (OS: p = 0.010 and DFS: p = 0.301) and Korean-KIRP (OS: p = 1.38 × 10-7 and DFS: p = 1.21 × 10-7) databases.

CONCLUSIONS

We discovered and verified genes specific for the survival of PRCC patients in the TCGA-KIRP and Korean-KIRP databases. The survival gene signature, including PCSK2 commonly obtained from the 40 gene signature of TCGA and the 10 gene signature of the Korean database, is expected to provide insight into predicting the survival of PRCC patients and developing new treatment.

Genome-wide association study and fine-mapping on Korean biobank to discover renal trait-associated variants

BACKGROUND

Chronic kidney disease is a significant health burden worldwide, with increasing incidence. Although several genome- wide association studies (GWAS) have investigated single nucleotide polymorphisms (SNP) associated with kidney trait, most studies were focused on European ancestry.

METHODS

We utilized clinical and genetic information collected from the Korean Genome and Epidemiology Study (KoGES).

RESULTS

More than five million SNPs from 58,406 participants were analyzed. After meta-GWAS, 1,360 loci associated with estimated glomerular filtration rate (eGFR) at a genome-wide significant level (p = 5 × 10-8) were identified. Among them, 399 loci were validated with at least one other biomarker (blood urea nitrogen [BUN] or eGFRcysC) and 149 loci were validated using both markers. Among them, 18 SNPs (nine known ones and nine novel ones) with 20 putative genes were found. The aggregated effect of genes estimated by MAGMA gene analysis showed that these significant genes were enriched in kidney-associated pathways, with the kidney and liver being the most enriched tissues.

CONCLUSION

In this study, we conducted GWAS for more than 50,000 Korean individuals and identified several variants associated with kidney traits, including eGFR, BUN, and eGFRcysC. We also investigated functions of relevant genes using computational methods to define putative causal variants.

A novel approach of kinship determination based on the physical length of genetically shared regions of chromosomes

BACKGROUND

Determination of genetic relatedness between individuals plays a crucial role in resolving numerous civil cases involving familial relationships and in forensic investigation concerning missing persons. Short tandem repeats (STRs), known for their high degree of DNA polymorphism, have traditionally been the primary choice of DNA markers in genetic testing, but their application for kinships testing is limited to cases involving close kinship. SNPs have emerged as promising supplementary markers for kinship determination. Nevertheless, the challenging remains in discriminating between third-degree or more distant relatives, such as first cousins, using SNPs.

OBJECTIVE

To investigate a kinship analysis method for distant degree of familial relationships using high-density SNP data.

METHODS

A high-density SNP data from 337 individuals of Korean families using Affymetrix Axiom KORV1.0-96 Array was obtained for this study. SNPs were aligned by chromosomal positions, and identity-by-state (IBS) was determined, and then shared regions as consecutive SNPs with IBS of 1 or 2 were investigated. The physical lengths of these IBS segments were measured and summed them to create an Index, as a measure of kinship.

RESULTS

The kinship was determined by the physical length of shared chromosomal regions that are distinguished by each kinship. Using this method, the relationship was able be distinguished up to the fourth degree of kinship, and non-relatives were clearly distinguished from true relatives. We also found a potential for this approach to be used universally, regardless of microarray platforms for SNP genotyping and populations.

CONCLUSION

This method has a potential to determine the different degree of kinship between individuals and to distinguish non-relatives from true relatives, which can be of great help for practical applications in kinship determination.

Associations between plasma metabolites and heavy metal exposure in residents of environmentally polluted areas

Heavy metals are commonly released into the environment through industrial processes such as mining and refining. The rapid industrialization that occurred in South Korea during the 1960s and 1970s contributed significantly to the economy of the country; however, the associated mining and refining led to considerable environmental pollution, and although mining is now in decline in South Korea, the detrimental effects on residents inhabiting the surrounding areas remain. The bioaccumulation of toxic heavy metals leads to metabolic alterations in human homeostasis, with disruptions in this balance leading to various health issues. This study used metabolomics to explore metabolomic alterations in the plasma samples of residents living in mining and refining areas. The results showed significant increases in metabolites involved in glycolysis and the surrounding metabolic pathways, such as glucose-6-phosphate, phosphoenolpyruvate, lactate, and inosine monophosphate, in those inhabiting polluted areas. An investigation of the associations between metabolites and blood clinical parameters through meet-in-the-middle analysis indicated that female residents were more affected by heavy metal exposure, resulting in more metabolomic alterations. For women, inhabiting the abandoned mine area, metabolites in the glycolysis and pentose phosphate pathways, such as ribose-5-phosphate and 3-phosphoglycerate, have shown a negative correlation with albumin and calcium. Finally, Mendelian randomization(MR) was used to determine the causal effects of these heavy metal exposure-related metabolites on heavy metal exposure-related clinical parameters. Metabolite biomarkers could provide insights into altered metabolic pathways related to exposure to toxic heavy metals and improve our understanding of the molecular mechanisms underlying the health effects of toxic heavy metal exposure.

Placental expression quantitative trait loci in an East Asian population

Expression quantitative trait loci (eQTL) analysis measures the contribution of genetic variation in gene expression on complex traits. Although this methodology has been used to examine gene regulation in numerous human tissues, eQTL research in solid tissues is relatively lacking. We conducted eQTL analysis on placentas collected from an East Asian population in an effort to identify gene regulatory mechanisms in this tissue. Placentas (n = 102) were collected at the time of cesarean delivery. mRNA was extracted, sequenced with NGS, and compared with matched maternal and fetal DNA arrays performed using maternal and neonatal cord blood. Linear regression modeling was performed using tensorQTL. Fine-mapping along with epigenomic annotation was used to select putative functional variants. We identified 2,703 coding genes that contained at least one eQTL with statistical significance (false discovery rate <0.05). After fine-mapping, we found 108 previously unreported eQTL variants with posterior inclusion probability >0.1. Of these, 19% were located in genomic regions with evidence from public placental epigenome suggesting that they may be functionally relevant. For example, variant rs28379289 located in the placenta-specific regulatory region changes the binding affinity of transcription factor leading to higher expression of LGALS3, which is known to affect placental function. This study expands the knowledge base of regulatory elements within the human placenta and identifies 108 previously unreported placenta eQTL signals, which are listed in our publicly available GMI eQTL database. Further studies are needed to identify and characterize genetic regulatory mechanisms that affect placental function in normal pregnancy and placenta-related diseases.

Evaluating cardiovascular disease risk stratification using multiple-polygenic risk scores and pooled cohort equations: insights from a 17-year longitudinal Korean cohort study

Cardiovascular disease (CVD) remains the leading cause of mortality worldwide, caused by a complex interplay of genetic and environmental factors. This study aimed to evaluate the combined efficacy of multi-polygenic risk scores and pooled cohort equations (PCE) for predicting future CVD risks in the Korean population. In this longitudinal study, 7,612 individuals from the Ansan and Ansung cohorts were analyzed over a 17-year follow-up period. The participants were genotyped using the Korea Biobank Array, and quality-controlled genetic data were subjected to imputation analysis. The weighted sum of the PRSs (wPRSsum) was calculated using PRS-CS with summary statistics from myocardial infarction, ischemic stroke, coronary artery disease, and hypertension genome-wide association studies. The recalibrated PCE was used to assess clinical risk, and the participants were stratified into risk groups based on the wPRSsum and PCE. Associations between these risk scores and incident CVD were evaluated using Cox proportional hazards models and Kaplan-Meier analysis. The wPRSsum approach showed a significant association with incident CVD (HR = 1.15, p = 7.49 × 10-5), and the top 20% high-risk genetic group had an HR of 1.50 (p = 5.04 × 10-4). The recalibrated PCE effectively differentiated between the low and high 10-year CVD risk groups, with a marked difference in survival rates. The predictive models constructed using the wPRSsum and PCE demonstrated a slight improvement in prediction accuracy, particularly among males aged <55 years (C-index = 0.640). We demonstrated that while the integration of wPRSsum with PCE did not significantly outperform the PCE-only model (C-index: 0.703 for combined and 0.704 for PCE-only), it provided enhanced stratification of CVD risk. The highest risk group, identified through the combination of high wPRSsum and PCE scores, exhibited an HR of 4.99 for incident CVD (p = 1.45 × 10-15). These findings highlight the potential of integrating genetic risk assessments with traditional clinical tools for effective CVD risk stratification. Although the addition of wPRSsum to the PCE provided a marginal predictive improvement, it proved valuable in identifying high-risk individuals and supporting personalized treatment strategies. This study reinforces the utility of multi-PRS in conjunction with clinical risk assessment tools, paving the way for more tailored approaches for CVD prevention and management in diverse populations.

Risk prediction and interaction analysis using polygenic risk score of type 2 diabetes in a Korean population

Joint modelling of genetic and environmental risk factors can provide important information to predict the risk of type 2 diabetes (T2D). Therefore, to predict the genetic risk of T2D, we constructed a polygenic risk score (PRS) using genotype data of one Korean cohort, KARE (745 cases and 2549 controls), and the genome-wide association study summary statistics of Biobank Japan. We evaluated the performance of PRS in an independent Korean cohort, HEXA (5684 cases and 35,703 controls). Individuals with T2D had a significantly higher mean PRS than controls (0.492 vs. - 0.078, p ≈ 0 ). PRS predicted the risk of T2D with an AUC of 0.658 (95% CI 0.651-0.666). We also evaluated interaction between PRS and waist circumference (WC) in the HEXA cohort. PRS exhibited a significant sub-multiplicative interaction with WC (ORinteraction 0.991, 95% CI 0.987-0.995, pinteraction = 4.93 × 10-6) in T2D. The effect of WC on T2D decreased as PRS increased. The sex-specific analyses produced similar interaction results, revealing a decreased WC effect on T2D as the PRS increased. In conclusion, the risk of WC for T2D may differ depending on PRS and those with a high PRS might develop T2D with a lower WC threshold. Our findings are expected to improve risk prediction for T2D and facilitate the identification of individuals at an increased risk of T2D.

The clinical relevance of a polygenic risk score for type 2 diabetes mellitus in the Korean population

The clinical utility of a type 2 diabetes mellitus (T2DM) polygenic risk score (PRS) in the East Asian population remains underexplored. We aimed to examine the potential prognostic value of a T2DM PRS and assess its viability as a clinical instrument. We first established a T2DM PRS for 5490 Korean individuals using East Asian Biobank data (269,487 samples). Subsequently, we assessed the predictive capability of this T2DM PRS in a prospective longitudinal study with baseline data and data from seven additional follow-ups. Our analysis showed that the T2DM PRS could predict the transition of glucose tolerance stages from normal glucose tolerance to prediabetes and from prediabetes to T2DM. Moreover, T2DM patients in the top-decile PRS group were more likely to be treated with insulin (hazard ratio = 1.69, p value = 2.31E-02) than were those in the remaining PRS groups. T2DM PRS values were significantly high in the severe diabetes subgroup, characterized by insulin resistance and β -cell dysfunction (p value = 0.0012). The prediction models with the T2DM PRS had significantly greater Harrel's C-indices than did corresponding models without it. By utilizing prospective longitudinal study data and extensive clinical risk factor information, our analysis provides valuable insights into the multifaceted clinical utility of the T2DM PRS.

Influence of blood pressure polygenic risk scores and environmental factors on coronary artery disease in the Korean Genome and Epidemiology Study

The present study aimed to investigate the association of blood pressure polygenic risk scores (BP PRSs) with coronary artery disease (CAD) in a Korean population and the interaction effects between PRSs and environmental factors on CAD. Data were derived from the Cardiovascular Disease Association Study (CAVAS; N = 5100) and the Health Examinee Study (HEXA; N = 58,623) within the Korean Genome and Epidemiology Study. PRSs for systolic and diastolic BP were calculated with the weighted allele sum of >200 single-nucleotide polymorphisms. Multivariable logistic regression models were used. BP PRSs were strongly associated with systolic BP (SBP), diastolic BP (DBP), and hypertension in both CAVAS and HEXA (p < 0.0001). PRSSBP was significantly associated with CAD in CAVAS, while PRSSBP and PRSDBP were significantly associated with CAD in HEXA. There was an interaction effect between the BP PRSs and environmental factors on CAD. The odds ratios (ORs) for CAD were 1.036 (95% confidence interval [CI], 1.016-1.055) for obesity, 1.028 (95% CI, 1.011-1.045) for abdominal obesity, 1.030 (95% CI, 1.009-1.050) for triglyceride, 1.024 (95% CI, 1.008-1.041) for high-density lipoprotein cholesterol, and 1.039 for smoking (95% CI, 1.003-1.077) in CAVAS. There was no significant interaction in HEXA, except between PRSDBP and triglyceride (OR, 1.012; 95% CI, 1.001-1.024). BP PRS was associated with an increased risk of hypertension and CAD. The interactions among PRSs and environmental risk factors increased the risk of CAD. Multi-component interventions to lower BP in the population via healthy behaviors are needed to prevent CAD regardless of genetic predisposition.

Shared genetic architectures of educational attainment in East Asian and European populations

Educational attainment (EduYears), a heritable trait often used as a proxy for cognitive ability, is associated with various health and social outcomes. Previous genome-wide association studies (GWASs) on EduYears have been focused on samples of European (EUR) genetic ancestries. Here we present the first large-scale GWAS of EduYears in people of East Asian (EAS) ancestry (n = 176,400) and conduct a cross-ancestry meta-analysis with EduYears GWAS in people of EUR ancestry (n = 766,345). EduYears showed a high genetic correlation and power-adjusted transferability ratio between EAS and EUR. We also found similar functional enrichment, gene expression enrichment and cross-trait genetic correlations between two populations. Cross-ancestry fine-mapping identified refined credible sets with a higher posterior inclusion probability than single population fine-mapping. Polygenic prediction analysis in four independent EAS and EUR cohorts demonstrated transferability between populations. Our study supports the need for further research on diverse ancestries to increase our understanding of the genetic basis of educational attainment.

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.

Genetic Variants Linked to Myocardial Infarction in Individuals with Non-Alcoholic Fatty Liver Disease and Their Potential Interaction with Dietary Patterns

In recent studies, non-alcoholic fatty liver disease (NAFLD) has been associated with a high risk of ischemic heart disease. This study aimed to investigate a genetic variant within a specific gene associated with myocardial infarction (MI) among patients with NAFLD. We included 57,205 participants from a Korean genome and epidemiology study. The baseline population consisted of 45,400 individuals, with 11,805 identified as patients with NAFLD. Genome-wide association studies were conducted for three groups: the entire sample, the healthy population, and patients with NAFLD. We defined the p-value < 1 × 10-5 as the nominal significance and the p-value < 5 × 10-2 as statistically significant for the gene-by-nutrient interaction. Among the significant single-nucleotide polymorphisms (SNPs), the lead SNP of each locus was further analyzed. In this cross-sectional study, a total of 1529 participants (2.8%) had experienced MI. Multivariable logistic regression was performed to evaluate the association of 102 SNPs across nine loci. Nine SNPs (rs11891202, rs2278549, rs13146480, rs17293047, rs184257317, rs183081683, rs1887427, rs146939423, and rs76662689) demonstrated an association with MI in the group with NAFLD Notably, the MI-associated SNP, rs134146480, located within the SORCS2 gene, known for its role in secreting insulin in islet cells, showed the most significant association with MI (p-value = 2.55 × 10-7). Our study identifies candidate genetic polymorphisms associated with NAFLD-related MI. These findings may serve as valuable indicators for estimating MI risk and for conducting future investigations into the underlying mechanisms of NAFLD-related MI.

Genome-wide association study implicates the role of TBXAS1 in the pathogenesis of depressive symptoms among the Korean population

Although depression is an emerging disorder affecting many people worldwide, most genetic studies have been performed in European descent populations. Herein, a genome-wide association study (GWAS) was conducted in Korean population to elucidate the genomic loci associated with depressive symptoms. Two independent cohorts were used as discovery datasets, which consisted of 6474 (1484 cases and 4990 controls) and 1654 (557 cases and 1097 controls) Korean participants, respectively. The participants were divided into case and control groups based on the Beck Depression Inventory (BDI). Meta-analysis using the two cohorts revealed that rs6945590 was significantly associated with the risk of depressive symptoms [P = 2.83 × 10-8; odds ratio (OR) = 1.23; 95% confidence interval (CI): 1.15-1.33]. This association was validated in other independent cohorts which were another Korean cohort (258 cases and 1757 controls) and the East Asian study of the Psychiatric Genomics Consortium (PGC) (12,455 cases and 85,548 controls). The predicted expression levels of thromboxane A synthase 1 gene (TBXAS1), which encodes the enzyme thromboxane A synthase 1 and participates in the arachidonic acid (AA) cascade, was significantly decreased in the whole blood tissues of the participants with depressive symptoms. Furthermore, Mendelian randomization (MR) analysis showed a causal association between TBXAS1 expression and the risk of depressive symptoms. In conclusion, as the number of risk alleles (A) of rs6945590 increased, TBXAS1 expression decreased, which subsequently caused an increase in the risk of depressive symptoms.

A Prospective Sociocentric Study of 2 Entire Traditional Korean Villages: The Korean Social Life, Health, and Aging Project (KSHAP)

The Korean Social Life, Health, and Aging Project (KSHAP) was a multidisciplinary prospective study conducted in South Korea that measured various health biomarkers from blood, hair, and brain magnetic resonance imaging, and we examined their associations with sociocentric (global) social network data of older adults in 2 entire villages (or cohorts). Cohort K included participants aged 60 years or older, and cohort L included participants aged 65 years or older. We performed a baseline survey involving 814 of the 860 individuals (94.7% response rate) in cohort K in 2012 and 947 of the 1,043 individuals (90.8% response rate) in cohort L in 2017. We gathered longitudinal data for 5 waves in cohort K from 2011 to 2019 and 2 waves in cohort L from 2017 to 2022. Here, we describe for the first time the follow-up design of the KSHAP, the changes in social networks, and various biomarkers over a number of years. The data for cohort K are publicly available via the Korean Social Science Data Archive as well as the project website, and the data for cohort L will be shared soon.

Evaluation of penalized and machine learning methods for asthma disease prediction in the Korean Genome and Epidemiology Study (KoGES)

BACKGROUND

Genome-wide association studies have successfully identified genetic variants associated with human disease. Various statistical approaches based on penalized and machine learning methods have recently been proposed for disease prediction. In this study, we evaluated the performance of several such methods for predicting asthma using the Korean Chip (KORV1.1) from the Korean Genome and Epidemiology Study (KoGES).

RESULTS

First, single-nucleotide polymorphisms were selected via single-variant tests using logistic regression with the adjustment of several epidemiological factors. Next, we evaluated the following methods for disease prediction: ridge, least absolute shrinkage and selection operator, elastic net, smoothly clipped absolute deviation, support vector machine, random forest, boosting, bagging, naïve Bayes, and k-nearest neighbor. Finally, we compared their predictive performance based on the area under the curve of the receiver operating characteristic curves, precision, recall, F1-score, Cohen's Kappa, balanced accuracy, error rate, Matthews correlation coefficient, and area under the precision-recall curve. Additionally, three oversampling algorithms are used to deal with imbalance problems.

CONCLUSIONS

Our results show that penalized methods exhibit better predictive performance for asthma than that achieved via machine learning methods. On the other hand, in the oversampling study, randomforest and boosting methods overall showed better prediction performance than penalized methods.

Coffee consumption affects kidney function based on GCKR polymorphism in a Korean population

Kidney function can be preserved through pharmacological interventions and nonpharmacological strategies, such as lifestyle and dietary adjustments. Among these, coffee has been linked to protective effects on kidney function. However, few studies have investigated the effect of coffee consumption on kidney function according to specific genes. We hypothesized that the impact of coffee consumption on kidney function might vary depending on GCKR polymorphism. GCKR rs1260326 polymorphism was examined using the Korean genome and epidemiology data from 656 chronic kidney disease (CKD) cases and 38,540 individuals without CKD (non-CKD). GCKR polymorphism has been previously associated with both coffee consumption and kidney function in Europeans. We replicated the associations between GCKR rs1260326 and coffee consumption and kidney function in Korean individuals. We also explored the effect of coffee consumption on kidney function by multivariate logistic regression analysis. Individuals with the rs1260326 (TC/CC) genotype did not experience significant changes in CKD risk based on their coffee consumption habits. In contrast, individuals with the TT genotype exhibited a significantly lower risk of CKD based on coffee consumption. Interestingly, in the non-CKD group, a beneficial effect on estimated glomerular filtration rate was observed in individuals with the T allele as coffee consumption increased. Our findings supported the hypothesis and revealed that the impact of coffee consumption habits on kidney function may vary based on the GCKR rs1260326 genotype of Korean individuals.

Regulatory Variants on the Leukocyte Immunoglobulin-Like Receptor Gene Cluster are Associated with Crohn's Disease and Interact with Regulatory Variants for TAP2

BACKGROUND AND AIMS

Crohn's disease [CD] has a complex polygenic aetiology with high heritability. There is ongoing effort to identify novel variants associated with susceptibility to CD through a genome-wide association study [GWAS] in large Korean populations.

METHODS

Genome-wide variant data from 902 Korean patients with CD and 72 179 controls were used to assess the genetic associations in a meta-analysis with previous Korean GWAS results from 1621 patients with CD and 4419 controls. Epistatic interactions between CD-risk variants of interest were tested using a multivariate logistic regression model with an interaction term.

RESULTS

We identified two novel genetic associations with the risk of CD near ZBTB38 and within the leukocyte immunoglobulin-like receptor [LILR] gene cluster [p < 5 × 10-8], with highly consistent effect sizes between the two independent Korean cohorts. CD-risk variants in the LILR locus are known quantitative trait loci [QTL] for multiple LILR genes, of which LILRB2 directly interacts with various ligands including MHC class I molecules. The LILR lead variant exhibited a significant epistatic interaction with CD-associated regulatory variants for TAP2 involved in the antigen presentation of MHC class I molecules [p = 4.11 × 10-4], showing higher CD-risk effects of the TAP2 variant in individuals carrying more risk alleles of the LILR lead variant (odds ratio [OR] = 0.941, p = 0.686 in non-carriers; OR = 1.45, p = 2.51 × 10-4 in single-copy carriers; OR = 2.38, p = 2.76 × 10-6 in two-copy carriers).

CONCLUSIONS

This study demonstrated that genetic variants at two novel susceptibility loci and the epistatic interaction between variants in LILR and TAP2 loci confer a risk of CD.

Polygenic risk score model for renal cell carcinoma in the Korean population and relationship with lifestyle-associated factors

BACKGROUND

The polygenic risk score (PRS) is used to predict the risk of developing common complex diseases or cancers using genetic markers. Although PRS is used in clinical practice to predict breast cancer risk, it is more accurate for Europeans than for non-Europeans because of the sample size of training genome-wide association studies (GWAS). To address this disparity, we constructed a PRS model for predicting the risk of renal cell carcinoma (RCC) in the Korean population.

RESULTS

Using GWAS analysis, we identified 43 Korean-specific variants and calculated the PRS. Subsequent to plotting receiver operating characteristic (ROC) curves, we selected the 31 best-performing variants to construct an optimal PRS model. The resultant PRS model with 31 variants demonstrated a prediction rate of 77.4%. The pathway analysis indicated that the identified non-coding variants are involved in regulating the expression of genes related to cancer initiation and progression. Notably, favorable lifestyle habits, such as avoiding tobacco and alcohol, mitigated the risk of RCC across PRS strata expressing genetic risk.

CONCLUSION

A Korean-specific PRS model was established to predict the risk of RCC in the underrepresented Korean population. Our findings suggest that lifestyle-associated factors influencing RCC risk are associated with acquired risk factors indirectly through epigenetic modification, even among individuals in the higher PRS category.

Advances in RNA Interference for Plant Functional Genomics: Unveiling Traits, Mechanisms, and Future Directions

RNA interference (RNAi) is a conserved molecular mechanism that plays a critical role in post-transcriptional gene silencing across diverse organisms. This review delves into the role of RNAi in plant functional genomics and its applications in crop improvement, highlighting its mechanistic insights and practical implications. The review begins with the foundational discovery of RNAi's mechanism, tracing its origins from petunias to its widespread presence in various organisms. Various classes of regulatory non-coding small RNAs, including siRNAs, miRNAs, and phasiRNAs, have been uncovered, expanding the scope of RNAi-mediated gene regulation beyond conventional understanding. These RNA classes participate in intricate post-transcriptional and epigenetic processes that influence gene expression. In the context of crop enhancement, RNAi has emerged as a powerful tool for understanding gene functions. It has proven effective in deciphering gene roles related to stress resistance, metabolic pathways, and more. Additionally, RNAi-based approaches hold promise for integrated pest management and sustainable agriculture, contributing to global efforts in food security. This review discusses RNAi's diverse applications, such as modifying plant architecture, extending shelf life, and enhancing nutritional content in crops. The challenges and future prospects of RNAi technology, including delivery methods and biosafety concerns, are also explored. The global landscape of RNAi research is highlighted, with significant contributions from regions such as China, Europe, and North America. In conclusion, RNAi remains a versatile and pivotal tool in modern plant research, offering novel avenues for understanding gene functions and improving crop traits. Its integration with other biotechnological approaches such as gene editing holds the potential to shape the future of agriculture and sustainable food production.

Ethnic-specific associations between body mass index and gastric cancer: a Mendelian randomization study in European and Korean populations

BACKGROUND

Given the uncertainties surrounding the associations in previous epidemiological studies, we conducted linear and nonlinear Mendelian randomization (MR) studies to evaluate whether body mass index (BMI) associated with gastric cancer (GC) risk in European and Korean.

METHODS

Genome-wide association study-summary statistics were used from the Pan-UK Biobank, the Genetic Investigation of Anthropometric Traits consortium, the K-CHIP consortium, and BioBank Japan. BMI-associated single nucleotide polymorphisms (SNPs) were used as instrumental variables (IVs) in MR to identify the association between BMI and GC. Both linear and nonlinear MR analyses were performed. Sensitivity analyses were also conducted for individuals below or above a BMI of 24 kg/m2.

RESULTS

The study used 22 and 55 SNPs as IVs for BMI in European and Korean populations, respectively. Genetically predicted BMI was positively associated with GC risk in the European population (Odds ratio per 1 kg/m2 increase; 95% CI = 1.17; 1.01-1.36 using simple median method), but no significant association was observed in the Korean population. However, the nonlinear MR identified a U-shaped association between BMI and GC in the Korean population, with both low and high BMIs associated with increased GC risk. A BMI of 24 kg/m2 presented the lowest risk. Sensitivity analyses did not yield any genome-wide significant SNPs.

CONCLUSION

While MR analysis suggests a linear association between BMI and GC in those of European ancestry, nonlinear MR hints at a U-shaped association in Koreans. This suggests the association between BMI and GC risk may vary according to ethnic ancestry.

Single-nucleotide polymorphisms link gout with health-related lifestyle factors in Korean cohorts

Gout-a very painful inflammatory arthritis caused by the deposition of monosodium urate crystals in the joints-is influenced by several factors. We identified the association of single- nucleotide polymorphisms (SNPs) that link gout with health-related lifestyle factors using genomic data from the Korean Genome and Epidemiology Study. We conducted a genome-wide association study (GWAS) on 18,927 samples of 438 Korean patients with gout and 18,489 controls for the discovery stage. For the replication stage, another batch containing samples of 326 patients with gout and 2,737 controls were analyzed. Lastly, a meta-analysis was performed using these two cohorts. We analyzed the effects of health-related lifestyle factors, including eating habits, physical activity, drinking behavior, and smoking behavior, on gout. After identifying the association between GWAS-derived SNPs and health-related lifestyle factors, we confirmed the interaction between the polygenic risk score (PRS) and health-related lifestyle factors. We identified 15 SNPs related to gout, among which rs1481012 of ABCG2 located on chromosome 4 has been newly discovered (P = 2.46e-11). On examining the interaction between SNPs and health-related lifestyles, rs3109823-located in ABCG2-was found to be associated with smoking status. In addition, rs11936395-located in SLC2A9-was significantly associated with the average momentum of exercise per session, whereas rs11066325 located in PTPN11, showed a significant association with the number of exercise sessions per week, smoking status, drinking status, and amount of soju drink per session. rs9421589-located in FAM35A-was significantly associated with the duration of smoking. In addition, we verified that the association between PRS and duration of smoking affects gout. Thus, in this study, we identified novel SNPs that link gout with health-related lifestyle factors in the Korean population.

Genetic variant rs9939609 in FTO is associated with body composition and obesity risk in Korean females

INTRODUCTION

The fat mass and obesity-associated (FTO) gene is a significant locus in obesity. However, the association between FTO genetic variants and body composition has not been fully elucidated.

RESEARCH DESIGN AND METHODS

This observational study examined the associations of FTO rs9939609 T>A with obesity and body composition markers in Koreans. A total of 6474 participants from the Korean Genome and Epidemiology Study were analyzed for their general characteristics, body composition and FTO genotype with a sex-stratified approach.

RESULTS

Females with the obesity risk A allele showed significantly greater body weight, hip circumference, and body mass index and were at a 1.28-fold higher risk of obesity (95% CI=1.088 to 1.507) than those with the TT genotype. Analyses of body composition also showed that females with the A allele had a greater body fat mass and percentage, abdominal fat percentage, and degree of obesity, and this association and FTO genetic variation and adiposity was observed in females, especially aged under 50 years. However, the effect of the variant allele on non-fat tissue markers was not evident in females and was not associated with any parameters examined in males.

CONCLUSIONS

The FTO rs9939609 variant is associated with body composition in Koreans, especially body fat markers in females. These results support that the FTO rs9939609 variant is a genetic risk factor in the etiology of obesity.

Evaluation of optimal methods and ancestries for calculating polygenic risk scores in East Asian population

Polygenic risk scores (PRSs) have been studied for predicting human diseases, and various methods for PRS calculation have been developed. Most PRS studies to date have focused on European ancestry, and the performance of PRS has not been sufficiently assessed in East Asia. Herein, we evaluated the predictive performance of PRSs for East Asian populations under various conditions. Simulation studies using data from the Korean cohort, Health Examinees (HEXA), demonstrated that SBayesRC and PRS-CS outperformed other PRS methods (lassosum, LDpred-funct, and PRSice) in high fixed heritability (0.3 and 0.7). In addition, we generated PRSs using real-world data from HEXA for ten diseases: asthma, breast cancer, cataract, coronary artery disease, gastric cancer, glaucoma, hyperthyroidism, hypothyroidism, osteoporosis, and type 2 diabetes (T2D). We utilized the five previous PRS methods and genome-wide association study (GWAS) data from two biobank-scale datasets [European (UK Biobank) and East Asian (BioBank Japan) ancestry]. Additionally, we employed PRS-CSx, a PRS method that combines GWAS data from both ancestries, to generate a total of 110 PRS for ten diseases. Similar to the simulation results, SBayesRC showed better predictive performance for disease risk than the other methods. Furthermore, the East Asian GWAS data outperformed those from European ancestry for breast cancer, cataract, gastric cancer, and T2D, but neither of the two GWAS ancestries showed a significant advantage on PRS performance for the remaining six diseases. Based on simulation data and real data studies, it is expected that SBayesRC will offer superior performance for East Asian populations, and PRS generated using GWAS from non-East Asian may also yield good results.

A genome-wide association scan reveals novel loci for facial traits of Koreans

DNA-based prediction of externally visible characteristics (EVC) with SNPs is one of the research areas of interest in the forensic field. Based on a previous study performing GWAS on facial traits in a Korean population, herein, we present results stemming from GWA analysis with KoreanChip and novel genetic loci satisfying genome-wide significant level. We discovered a total of 20 signals and 12 loci were found to have novel associations with facial traits, including six loci located in intergenic regions and six loci located at UBE2O, HECTD2, CCDC108, TPK1, FCN2, and FRMPD1. Additionally, we performed a polygenic score analysis for 33 distance-related traits in facial phenotyping and determined genetic relationships between facial traits and SNPs using the GCTA program. The results of the current study offer an understanding of how facial morphology is influenced by complex genetic structures and provide insights into forensic investigation and population genetics.

Identification of asthma-related genes using asthmatic blood eQTLs of Korean patients

BACKGROUND

More than 200 asthma-associated genetic variants have been identified in genome-wide association studies (GWASs). Expression quantitative trait loci (eQTL) data resources can help identify causal genes of the GWAS signals, but it can be difficult to find an eQTL that reflects the disease state because most eQTL data are obtained from normal healthy subjects.

METHODS

We performed a blood eQTL analysis using transcriptomic and genotypic data from 433 Korean asthma patients. To identify asthma-related genes, we carried out colocalization, Summary-based Mendelian Randomization (SMR) analysis, and Transcriptome-Wide Association Study (TWAS) using the results of asthma GWASs and eQTL data. In addition, we compared the results of disease eQTL data and asthma-related genes with two normal blood eQTL data from Genotype-Tissue Expression (GTEx) project and a Japanese study.

RESULTS

We identified 340,274 cis-eQTL and 2,875 eGenes from asthmatic eQTL analysis. We compared the disease eQTL results with GTEx and a Japanese study and found that 64.1% of the 2,875 eGenes overlapped with the GTEx eGenes and 39.0% with the Japanese eGenes. Following the integrated analysis of the asthmatic eQTL data with asthma GWASs, using colocalization and SMR methods, we identified 15 asthma-related genes specific to the Korean asthmatic eQTL data.

CONCLUSIONS

We provided Korean asthmatic cis-eQTL data and identified asthma-related genes by integrating them with GWAS data. In addition, we suggested these asthma-related genes as therapeutic targets for asthma. We envisage that our findings will contribute to understanding the etiological mechanisms of asthma and provide novel therapeutic targets.

Associations between the Genetic Heritability of Dyslipidemia and Dietary Patterns in Korean Adults Based on Sex Differences

Dyslipidemia can be defined as an abnormality in serum lipid levels that is substantially linked to genetic variations and lifestyle factors, such as diet patterns, and has distinct sex-specific characteristics. We aimed to elucidate the genetic impact of dyslipidemia according to sex and explore the associations between genetic variants and dietary patterns in large-scale population-based cohorts. After performing genome-wide association studies (GWASs) in male, female, and entire cohorts, significant single nucleotide polymorphisms (SNPs) were identified in the three groups, and genetic risk scores (GRSs) were calculated by summing the risk alleles from the selected SNPs. After adjusting for confounding variables, the risk of dyslipidemia was 2.013-fold and 2.535-fold higher in the 3rd quartile GRS group in the male and female cohorts, respectively, than in the 1st quartile GRS group. While instant noodle and soft drink intake were significantly associated with GRS related to hyperlipidemia in male cohorts, coffee consumption was substantially related to GRS related to hyperlipidemia in female cohorts. Considering the influence of genetic factors and dietary patterns, the findings of this study suggest the potential for implementing sex-specific strategic interventions to avoid dyslipidemia.