Genome-wide association study identifies novel susceptible loci and evaluation of polygenic risk score for chronic obstructive pulmonary disease in a Taiwanese population.

BACKGROUND: Chronic Obstructive Pulmonary Disease (COPD) describes a group of progressive lung diseases causing breathing difficulties. While COPD development typically involves a complex interplay between genetic and environmental factors, genetics play a role in disease susceptibility. This study used genome-wide association studies (GWAS) and polygenic risk score (PRS) to elucidate the genetic basis for COPD in Taiwanese patients. RESULTS: GWAS was performed on a Taiwanese COPD case-control cohort with a sample size of 5,442 cases and 17,681 controls. Additionally, the PRS was calculated and assessed in our target groups. GWAS results indicate that although there were no single nucleotide polymorphisms (SNPs) of genome-wide significance, prominent COPD susceptibility loci on or nearby genes such as WWTR1, EXT1, INTU, MAP3K7CL, MAMDC2, BZW1/CLK1, LINC01197, LINC01894, and CFAP95 (C9orf135) were identified, which had not been reported in previous studies. Thirteen susceptibility loci, such as CHRNA4, AFAP1, and DTWD1, previously reported in other populations were replicated and confirmed to be associated with COPD in Taiwanese populations. The PRS was determined in the target groups using the summary statistics from our base group, yielding an effective association with COPD (odds ratio [OR] 1.09, 95% confidence interval [CI] 1.02-1.17, p = 0.011). Furthermore, replication a previous lung function trait PRS model in our target group, showed a significant association of COPD susceptibility with PRS of Forced Expiratory Volume in one second (FEV1)/Forced Vital Capacity (FCV) (OR 0.89, 95% CI 0.83-0.95, p = 0.001). CONCLUSIONS: Novel COPD-related genes were identified in the studied Taiwanese population. The PRS model, based on COPD or lung function traits, enables disease risk estimation and enhances prediction before suffering. These results offer new perspectives on the genetics of COPD and serve as a basis for future research.

The genome-wide association study of serum IgE levels demonstrated a shared genetic background in allergic diseases.

Immunoglobulin E (IgE) synthessis is highly related to a variety of atopic diseases, and several genome-wide association studies (GWASs) have demonstrated the association between genes and IgE level. In this study, we conducted the largest genome-wide association study of IgE involving a Taiwanese Han population. Eight independent variants exhibited genome-wide significance. Among them, an intronic SNP of CD28, rs1181388, and an intergenic SNP, rs1002957030, on 11q23.2 were identified as novel signals for IgE. Seven of the loci were replicated successfully in a meta-analysis using data on Japanese population. Among all the human leukocyte antigen (HLA) regions, HLA-DQA1*03:02 - HLA-DQB1*03:03 was the most significant haplotype (OR = 1.25, SE = 0.02, FDR = 1.6 × 10-14), corresponding to HLA-DQA1 Asp160 and HLA-DQB1 Leu87 amino acid residues. The genetic correlation showed significance between IgE and allergic diseases including asthma, atopic dermatitis, and pollinosis. IgE PRS was significantly correlated with total IgE levels. Furthermore, the top decile IgE polygenic risk score (PRS) group had the highest risk of asthma for the Taiwan Biobank and Biobank Japan cohorts. IgE PRS may be used to aid in predicting the occurrence of allergic reactions before symptoms occur and biomarkers are detectable. Our study provided a more comprehensive understanding of the impact of genomic variants, including complex HLA alleles, on serum IgE levels.

Genetic susceptibility to prostate cancer in Taiwan: A genome-wide association study.

We conducted the first genome-wide association study (GWAS) of prostate cancer (PCa) in Taiwan with 1844 cases and 80,709 controls. Thirteen independent single-nucleotide polymorphisms (SNPs) reached genome-wide significance (p < 5 × 10-8 ). Among these, three were distinct from previously identified loci: rs76072851 in CORO2B gene (15q23), odds ratio (OR) = 1.54, 95% confidence interval (CI), 1.36-1.76, p = 5.30 × 10-11 ; rs7837051, near two long noncoding RNA (lncRNA) genes, PRNCR1 and PCAT2 (8q24.21), OR = 1.41 (95% CI, 1.31-1.51), p = 8.77 × 10-21 ; and rs56339048, near an lncRNA gene, CASC8 (8q24.21), OR = 1.25 (95% CI, 1.16-1.35), p = 2.14 × 10-8 . We refined the lead SNPs for two previously identified SNPs in Taiwanese: rs13255059 (near CASC8), p = 9.02 × 10-43 , and rs1456315 (inside PRNCR1), p = 4.33 × 10-42 . We confirmed 35 out of 49 GWAS-identified East Asian PCa susceptibility SNPs. In addition, we identified two SNPs more specific to Taiwanese than East Asians: rs34295433 in LAMC1 (1q25.3) and rs6853490 in PDLIM5 (4q22.3). A weighted genetic risk score (GRS) was developed using the 40 validated SNPs and the area under the receiver-operating characteristic curve for the GRS to predict PCa was 0.67 (95% CI, 0.63-0.71). These identified SNPs provide valuable insights into the molecular mechanisms of prostate carcinogenesis in Taiwan and underscore the significant role of genetic susceptibility in regional differences in PCa incidence.

Ethnic-specific genetic susceptibility loci for endometriosis in Taiwanese-Han population: a genome-wide association study.

Endometriosis is a common gynecological disorder affecting around 10% of reproductive-age women. Although many hypotheses were proposed, genetic alteration has been considered as one of the key factors promoting pathogenesis. Due to racial/ethnic disparities in the process of hormone regulation and nutrition metabolism, a genome-wide association study (GWAS) with 2794 cases and 27,940 controls was conducted in a Taiwanese-Han population. Our study identified five significant susceptibility loci for endometriosis, and three of them, WNT4 (on the 1p36.12), RMND1 (6q25.1), and CCDC170 (6q25.1), have been previously associated with endometriosis across different populations, including European and Japanese descent cohorts. Other two including C5orf66/C5orf66-AS2 (5q31.1) and STN1 (10q24.33) are newly identified ones. Functional network analysis of potent risk genes revealed the involvement of cancer susceptibility and neurodevelopmental disorders in endometriosis development. In addition, long non-coding RNAs (lncRNAs) C5orf66 and C5orf66-AS2 can interact with many RNA-binding proteins (RBPs) which can influence RNA metabolic process, mRNA stabilization, and mRNA splicing, leading to dysregulation in tumor-promoting gene expression. Those findings support clinical observations of differences in the presentation of endometriosis in Taiwanese-Han population with higher risks of developing deeply infiltrating/invasive lesions and the associated malignancies.

Menarche-a journey into womanhood: age at menarche and health-related outcomes in East Asians.

STUDY QUESTION: Are there associations of age at menarche (AAM) with health-related outcomes in East Asians? SUMMARY ANSWER: AAM is associated with osteoporosis, Type 2 diabetes (T2D), glaucoma, and uterine fibroids, as demonstrated through observational studies, polygenic risk scores, genetic correlations, and Mendelian randomization (MR), with additional findings indicating a causal effect of BMI and T2D on earlier AAM. WHAT IS KNOWN ALREADY: Puberty timing is linked to adult disease risk, but research predominantly focuses on European populations, with limited studies in other groups. STUDY DESIGN, SIZE, DURATION: We performed an AAM genome-wide association study (GWAS) with 57 890 Han Taiwanese females and examined the association between AAM and 154 disease outcomes using the Taiwanese database. Additionally, we examined genetic correlations between AAM and 113 diseases and 67 phenotypes using Japanese GWAS summary statistics. PARTICIPANTS/MATERIALS, SETTING, METHODS: We performed AAM GWAS and gene-based GWAS studies to obtain summary statistics and identify potential AAM-related genes. We applied phenotype, polygenic risk scores, and genetic correlation analyses of AAM to explore health-related outcomes, using multivariate regression and linkage disequilibrium score regression analyses. We also explored potential bidirectional causal relationships between AAM and related outcomes through univariable and multivariable MR analyses. MAIN RESULTS AND THE ROLE OF CHANCE: Fifteen lead single-nucleotide polymorphisms and 24 distinct genes were associated with AAM in Taiwan. AAM was genetically associated with later menarche and menopause, greater height, increased osteoporosis risk, but lower BMI, and reduced risks of T2D, glaucoma, and uterine fibroids in East Asians. Bidirectional MR analyses indicated that higher BMI/T2D causally leads to earlier AAM. LIMITATIONS, REASONS FOR CAUTION: Our findings were specific to Han Taiwanese individuals, with genetic correlation analyses conducted in East Asians. Further research in other ethnic groups is necessary. WIDER IMPLICATIONS OF THE FINDINGS: Our study provides insights into the genetic architecture of AAM and its health-related outcomes in East Asians, highlighting causal links between BMI/T2D and earlier AAM, which may suggest potential prevention strategies for early puberty. STUDY FUNDING/COMPETING INTEREST(S): The work was supported by China Medical University, Taiwan (CMU110-S-17, CMU110-S-24, CMU110-MF-49, CMU111-SR-158, CMU111-MF-105, CMU111-MF-21, CMU111-S-35, CMU112-SR-30, and CMU112-MF-101), the China Medical University Hospital, Taiwan (DMR-111-062, DMR-111-153, DMR-112-042, DMR-113-038, and DMR-113-103), and the Ministry of Science and Technology, Taiwan (MOST 111-2314-B-039-063-MY3, MOST 111-2314-B-039-064-MY3, MOST 111-2410-H-039-002-MY3, and NSTC 112-2813-C-039-036-B). The funders had no influence on the data collection, analyses, or conclusions of the study. No conflict of interests to declare. TRIAL REGISTRATION NUMBER: N/A.

Identification of Genetic Variants for Risk Prediction and Early Diagnosis of Age-Related Macular Degeneration in the Taiwanese Population.

Age-related macular degeneration (AMD) is the leading cause of blindness in the elderly worldwide. The prevalence and phenotypes of AMD differ among populations, including between people in Taiwan and other regions. We performed a genome-wide association study to identify genetic variants and to develop genetic models to predict the risk of AMD development and progression in the Taiwanese population. In total, 4039 patients with AMD and 16,488 non-AMD controls (aged ≥ 65 years) were included. We identified 31 AMD-associated variants (p < 5 × 10-8) on chromosome 10q26, surrounding PLEKHA1-ARMS2-HTRA1. Two genetic models were constructed using the clump and threshold method. Model 1 included the single nucleotide polymorphism rs11200630 and showed a 1.31-fold increase in the risk of AMD per risk allele (95% confidence interval (CI) = 1.20-1.43, p < 0.001). In model 2, 1412 single-nucleotide polymorphisms were selected to construct a polygenic risk score (PRS). Individuals with the top 5% PRS had a 1.40-fold higher AMD risk compared with that of individuals with a PRS in the bottom quartile (95% CI = 1.04-1.89, p = 0.025). Moreover, the PRS in the upper quartile was related to a decreased age at AMD diagnosis by 0.62 years (95% CI = -1.15, -0.09, p = 0.023). Both genetic models provide useful predictive power for populations at high risk of AMD, affording a basis for identifying patients requiring close follow-up and early intervention.

Identification of the efficacy of parentage testing based on bi-allelic autosomal single nucleotide polymorphism markers in Taiwanese population.

Parentage testing is crucial for forensic DNA analysis, using short tandem repeats (STRs). Single nucleotide polymorphisms (SNPs) with high minor allele frequency (MAF) are promising for human identification. This study aimed to develop SNP markers for parentage testing in the Taiwanese population and compare their accuracy with STRs. The TPMv1 SNP microarray (714,457 SNPs) was used to screen 180,000 Taiwanese individuals and analyze the SNP data using PLINK. After quality control, allelic distribution, and MAF considerations, a set of SNPs with significant inheritance information was selected. Parentage testing was conducted on 355 single parent-child pairs using both STRs and SNPs, employing three kinship algorithms: identity by descent, kinship-based inference for genome-wide association studies, and the combined paternity index/probability of paternity (CPI/PP). An Affymetrix signature probe for kinship testing (ASP) was also used. Based on the quality control and selection criteria, 176 SNPs with MAF > 0.4995 were selected from the Taiwanese population. The CPI/PP results calculated using SNPs were consistent with the STR results. The accuracy of the SNPs used in the single-parent-child parentage testing was > 99.99%. The set of 176 SNPs had a higher identification rate in the single parent-child parentage test than in the ASP. The CPI/PP value calculated using 176 SNPs was also more accurate than that calculated using ASP. Our findings suggest that these 176 SNPs could be used for single-parent-child parentage identification in the Taiwanese population.

Genome-wide association study identifies DRAM1 associated with Tourette syndrome in Taiwan.

BACKGROUND: Tourette syndrome (TS) is a neurodevelopmental disorder characterized by motor and vocal tics. Several susceptibility loci associated with TS have been identified previously in populations of European descent using genome-wide association studies (GWAS). However, the exact pathogenic mechanism underlying TS is unknown; additionally, the results of previous GWAS for TS were based on Western populations, which may not translate to other populations. Therefore, we conducted a GWAS in Taiwanese patients with TS and chronic tic disorders (CTDs), with an aim to elucidate the genetic basis and potential risk factors for TS in this population. METHODS: GWAS was performed on a Taiwanese TS/CTDs cohort with a sample size of 1,007 patients with TS and 25,522 ancestry-matched controls. Additionally, polygenic risk score was calculated and assessed. RESULTS: Genome-wide significant locus, rs12313062 (p=1.43 × 10-8) and other 9 single nucleotide polymorphisms, were identified in chromosomes 12q23.2, associated with DRAM1 and was a novel susceptibility locus identified in TS/CTDs group. DRAM1, a lysosomal transmembrane protein regulated by p53, modulates autophagy and apoptosis, with potential implications for neuropsychiatric conditions associated with autophagy disruption. CONCLUSIONS: This study conducted the first GWAS for TS in a Taiwanese population, identifying a significant locus on chromosome 12q23.2 associated with DRAM1. These findings provide novel insights into the neurobiology of TS and potential directions for future research in this area.

Genome­wide association study and polygenic risk scores predict psoriasis and its shared phenotypes in Taiwan.

Psoriasis is a chronic inflammatory dermatological disease, and there is a lack of understanding of the genetic factors involved in psoriasis in Taiwan. To establish associations between genetic variations and psoriasis, a genome­wide association study was performed in a cohort of 2,248 individuals with psoriasis and 67,440 individuals without psoriasis. Using the ingenuity pathway analysis software, biological networks were constructed. Human leukocyte antigen (HLA) diplotypes and haplotypes were analyzed using Attribute Bagging (HIBAG)­R software and chi­square analysis. The present study aimed to assess the potential risks associated with psoriasis using a polygenic risk score (PRS) analysis. The genetic association between single nucleotide polymorphisms (SNPs) in psoriasis and various human diseases was assessed by phenome­wide association study. METAL software was used to analyze datasets from China Medical University Hospital (CMUH) and BioBank Japan (BBJ). The results of the present study revealed 8,585 SNPs with a significance threshold of P<5x10­8, located within 153 genes strongly associated with the psoriasis phenotype, particularly on chromosomes 5 and 6. This specific genomic region has been identified by analyzing the biological networks associated with numerous pathways, including immune responses and inflammatory signaling. HLA genotype analysis indicated a strong association between HLA­A*02:07 and HLA­C*06:02 in a Taiwanese population. Based on our PRS analysis, the risk of psoriasis associated with the SNPs identified in the present study was quantified. These SNPs are associated with various dermatological, circulatory, endocrine, metabolic, musculoskeletal, hematopoietic and infectious diseases. The meta­analysis results indicated successful replication of a study conducted on psoriasis in the BBJ. Several genetic loci are significantly associated with susceptibility to psoriasis in Taiwanese individuals. The present study contributes to our understanding of the genetic determinants that play a role in susceptibility to psoriasis. Furthermore, it provides valuable insights into the underlying etiology of psoriasis in the Taiwanese community.

Increased risk of hearing loss associated with MT-RNR1 gene mutations: a real-world investigation among Han Taiwanese Population.

BACKGROUND: Previous studies have implicated inherited mutations in mitochondrial DNA (mtDNA) in sensorineural hearing loss (SNHL). However, the definitive association between mitochondrial 12S rRNA (MT-RNR1) variants and hearing loss in the population has not been well established, particularly in Asia. The objective of this retrospective cohort study was to assess the association between MT-RNR1 variants and the risk of SNHL in patients in Taiwan. METHODS: The cohort included 306,068 participants from Taiwan between January 2003 and December 2020. Participants were classified based on genetic variants, particularly mitochondrial mutations (rs267606618, rs267606619, rs267606617). MT-RNR1 variant cases were matched 1:10 with non-mutant patients by age, gender, and visit year, excluding those with pre-existing hearing loss. The primary endpoint was SNHL, identified using specific ICD-TM codes with a 90% positive predictive value. Medication exposure history was determined via self-report or electronic medical records in the hospital. Cox proportional hazard regression models were used to assess the association between MT-RNR1 variants and hearing loss, adjusting for various covariates. Kaplan-Meier survival curves and log-rank tests compared hearing loss incidence between groups. RESULTS: The mean age of the mtDNA variants group is 32.4 years, with a standard deviation of 19.2 years. The incidence density of hearing loss for the mutation group was 36.42 per 10,000 person-years (95% Confidence Interval [CI], 27.21-47.73), which was higher than the 23.77per 10,000 person-years (95% CI, 21.32-26.42) in the wild-type group (p = 0.0036). Additionally, diabetes mellitus was associated with an increased risk of developing SNHL in individuals with MT-RNR1 variants (adjusted hazard ratio = 1.76 [95% CI, 1.00-3.09], p < 0.05). CONCLUSION: This study highlights the increased risk of hearing loss in patients carrying MT-RNR1 variants, particularly those with diabetes mellitus. Future research that integrates genetic and clinical data is crucial for developing more precise interventions to monitor and treat hearing loss in this vulnerable population.

Identification of novel genetic susceptibility loci for calcium-containing kidney stone disease by genome-wide association study and polygenic risk score in a Taiwanese population.

Approximately 80% of kidney stone diseases contain calcium. Inherited genetic factors are among the variables that influence the development of calcium-containing kidney stone diseases (CKSD). Previous genome-wide association studies (GWAS) on stone diseases have been reported worldwide; however, these are not focused on calcium-containing stones. We conducted a GWAS to identify germline genetic polymorphisms associated with CKSD in a Medical Center in Taiwan; hence, this study was based primarily on a hospital-based database. CKSD was diagnosed using the chart records. Patients infected with urea-splitting-microorganisms and those with at least two urinary pH value below 5.5 were excluded. None of the patients had cystic stones based on stone analysis. Those over 40 years of age with no history of CKSD and no microscopic hematuria on urinalysis were considered as controls. The DNA isolated from the blood of 14,934 patients (63.7% male and 36.3% female) with CKSD and 29,868 controls (10,830 men and 19,038 women) at a medical center was genotyped for approximately 714,457 single nucleotide polymorphisms (SNPs) with minor allele frequency of ≥ 0.05. We used PLINK 1.9 to calculate the polygenic risk score (PRS) to investigate the association between CKSD and controls. The accuracy of the PRS was verified by dividing it into the training and testing groups. The statistical analyses were calculated with the area under the curve (AUC) using IBM SPSS version 22. We identified 432 susceptibility loci that reached a genome-wide threshold of P < 1.0 × 10- 5. A total of 132 SNPs reached a threshold of P < 5 × 10- 8 using a stricter definition of significance on chromosomes 4, 13, 16, 17, and 18. At the top locus of our study, SNPs in DGKH, PDILT, BCAS3, and ABCG2 have been previously reported. RN7SKP27, HDAC4, PCDH15, AP003068.2, and NFATC1 were novel findings in this study. PRS was adjusted for sex and age, resulting in an AUC of 0.65. The number of patients in the top quartile of PRS was 1.39 folds in the risk of CKSD than patients in the bottom quartile. Our data identified the significance of GWAS for patients with CKSD in a hospital-based study. The PRS also had a high AUC for discriminating patients with CKSD from controls. A total of 132 SNP loci of SNPs significantly associated with the development of CKSD. This first survey, which focused on patients with CKSD, will provide novel insights specific to CKSD and its potential clinical biomarkers.

Multiple polygenic risk scores can improve the prediction of systemic lupus erythematosus in Taiwan.

OBJECTIVE: To identify new genetic variants associated with SLE in Taiwan and establish polygenic risk score (PRS) models to improve the early diagnostic accuracy of SLE. METHODS: The study enrolled 2429 patients with SLE and 48 580 controls from China Medical University Hospital in Taiwan. A genome-wide association study (GWAS) and PRS analyses of SLE and other three SLE markers, namely ANA, anti-double-stranded DNA antibody (dsDNA) and anti-Smith antibody (Sm), were conducted. RESULTS: Genetic variants associated with SLE were identified through GWAS. Some novel genes, which have been previously reported, such as RCC1L and EGLN3, were revealed to be associated with SLE in Taiwan. Multiple PRS models were established, and optimal cut-off points for each PRS were determined using the Youden Index. Combining the PRSs for SLE, ANA, dsDNA and Sm yielded an area under the curve of 0.64 for the optimal cut-off points. An analysis of human leucocyte antigen (HLA) haplotypes in SLE indicated that individuals with HLA-DQA1*01:01 and HLA-DQB1*05:01 were at a higher risk of being classified into the SLE group. CONCLUSIONS: The use of PRSs to predict SLE enables the identification of high-risk patients before abnormal laboratory data were obtained or symptoms were manifested. Our findings underscore the potential of using PRSs and GWAS in identifying SLE markers, offering promise for early diagnosis and prediction of SLE.

Utility of polygenic scores across diverse diseases in a hospital cohort for predictive modeling.

Polygenic scores estimate genetic susceptibility to diseases. We systematically calculated polygenic scores across 457 phenotypes using genotyping array data from China Medical University Hospital. Logistic regression models assessed polygenic scores' ability to predict disease traits. The polygenic score model with the highest accuracy, based on maximal area under the receiver operating characteristic curve (AUC), is provided on the GeneAnaBase website of the hospital. Our findings indicate 49 phenotypes with AUC greater than 0.6, predominantly linked to endocrine and metabolic diseases. Notably, hyperplasia of the prostate exhibited the highest disease prediction ability (P value = 1.01 × 10-19, AUC = 0.874), highlighting the potential of these polygenic scores in preventive medicine and diagnosis. This study offers a comprehensive evaluation of polygenic scores performance across diverse human traits, identifying promising applications for precision medicine and personalized healthcare, thereby inspiring further research and development in this field.

Novel Genes Associated with Atrial Fibrillation and the Predictive Models for AF Incorporating Polygenic Risk Scores and PheWAS-Derived Risk Factors.

BACKGROUND: Atrial fibrillation (AF), the most common atrial arrhythmia, presents varied clinical manifestations. Despite the identification of genetic loci associated with AF, particularly in specific populations, research within Asian ethnicities remains limited. This study aimed to develop predictive models for AF using AF-associated single-nucleotide polymorphisms (SNPs) from a Genome-Wide Association Study (GWAS) on a substantial cohort of Taiwanese individuals, evaluating the models' predictive efficacy. METHODS: Involving 75,121 subjects, including 5,694 AF patients and 69,427 normal controls with GWAS data, the study merged polygenic risk scores (PRS) from AF-associated SNPs with Phenome-wide association study (PheWAS)-derived risk factors. Advanced statistical and machine learning techniques were employed to develop and evaluate AF predictive models for discrimination and calibration. RESULTS: The study identified the top 30 significant SNPs associated with AF, predominantly on chromosomes 10 and 16, implicating genes like NEURL1, SH3PXD2A, INA, NT5C2, STN1, and ZFHX3. Notably, INA, NT5C2, and STN1 were newly linked to AF. The GWAS predictive power using PRS-CS analysis for AF exhibited an area under the curve (AUC) of 0.600 (P < 0.001), improving to 0.855 (P < 0.001) after adjusting for age and gender. PheWAS analysis showed the top 10 diseases associated with these genes were circulatory system diseases. CONCLUSIONS: Integrating genetic and phenotypic data enhanced the accuracy and clinical relevance of AF predictive models. The findings suggest promise for refining AF risk assessment, enabling personalized interventions, and reducing AF-related morbidity and mortality burdens.

Genetic predisposition to bone mineral density and their health conditions in East Asians.

Osteoporosis, a condition defined by low BMD (typically < -2.5 SD), causes a higher fracture risk and leads to significant economic, social, and clinical impacts. Genome-wide studies mainly in Caucasians have found many genetic links to osteoporosis, fractures, and BMD, with limited research in East Asians (EAS). We investigated the genetic aspects of BMD in 86 716 individuals from the Taiwan Biobank and their causal links to health conditions within EAS. A genome-wide association study (GWAS) was conducted, followed by observational studies, polygenic risk score assessments, and genetic correlation analyses to identify associated health conditions linked to BMD. GWAS and gene-based GWAS studies identified 78 significant SNPs and 75 genes related to BMD, highlighting pathways like Hedgehog, WNT-mediated, and TGF-ß. Our cross-trait linkage disequilibrium score regression analyses for BMD and osteoporosis consistently validated their genetic correlations with BMI and type 2 diabetes (T2D) in EAS. Higher BMD was linked to lower osteoporosis risk but increased BMI and T2D, whereas osteoporosis linked to lower BMI, waist circumference, hemoglobinA1c, and reduced T2D risk. Bidirectional Mendelian randomization analyses revealed that a higher BMI causally increases BMD in EAS. However, no direct causal relationships were found between BMD and T2D, or between osteoporosis and either BMI or T2D. This study identified key genetic factors for bone health in Taiwan, and revealed significant health conditions in EAS, particularly highlighting the genetic interplay between bone health and metabolic traits like T2D and BMI.

We investigated how genetics affect bone health and related conditions like diabetes and obesity in 86 716 EAS. Previously, most studies focused on Caucasian populations, but our work helps to understand these issues in EAS. Our findings show that stronger bones are linked to a lower chance of osteoporosis but a higher risk of obesity and T2D. On the other hand, those with osteoporosis tend to have lower body weight and a decreased risk of diabetes, illustrating a complex relationship between bone health and body metabolism. Future research will focus on deeper genetic interactions and developing targeted interventions for bone health and related metabolic disorders in EAS.