Genetic profiles of 103,106 individuals in the Taiwan Biobank provide insights into the health and history of Han Chinese

Multiple genetic polymorphisms are associated with the risk of metabolic syndrome, fatty liver, and airflow limitation: A Taiwan Biobank study

BACKGROUND

Links have been reported between the airflow limitation and both metabolic syndrome (MetS) and fatty liver (FL). Additionally, associations between genetic factors and risks of MetS, FL, and airflow limitation have been identified separately in different studies. Our study aims to simultaneously explore the association between specific single nucleotide polymorphisms (SNPs) of certain genes and the risk of the three associated diseases.

METHODS

In this retrospective cross-sectional nationwide study, 150,709 participants from the Taiwan Biobank (TWB) were enrolled. We conducted a genotype-phenotype association analysis of nine SNPs on seven genes (ApoE-rs429358, MBOAT7-rs641738, LEPR-rs1805096, APOC3-rs2854116, APOC3-rs2854117, PPP1R3B-rs4240624, PPP1R3B-rs4841132, TM6SF2-rs58542926, and IFNL4-rs368234815) using data from the TWB1.0 and TWB2.0 genotype dataset. Participants underwent a series of assessments including questionnaires, blood examinations, abdominal ultrasounds, and spirometry examinations.

RESULTS

MetS was associated with FL and airflow limitation. ApoE-rs429358, LEPR-rs1805096, APOC3-rs2854116, APOC3-rs2854117, PPP1R3B-rs4240624, PPP1R3B-rs4841132, and TM6SF2-rs58542926 were significantly associated with the risk of MetS. The cumulative impact of T alleles of ApoE-rs429358 and TM6SF2-rs58542926 on the risk of FL was observed (p-value for trend < 0.001). Individuals without MetS and airflow limitation carrying LEPR-rs1805096 G_G genotype exhibited a reduction in the forced expiratory volume in 1 s percentage prediction (Coefficient -35, 95 % confidence interval (CI) -69.7- -0.4), low forced vital capacity percentage prediction (Coefficient -41.6, 95 % CI -82.6- -0.6), and low vital capacity percentage prediction (Coefficient -42.2, 95 % CI -84.2- -0.1).

CONCLUSIONS

MetS significantly correlated with FL and airflow limitation. Multiple SNPs were notably associated with MetS. Specifically, T alleles of ApoE-rs429358 and TM6SF2-rs58542926 cumulatively increased the risk of FL. LEPR-rs1805096 shows a trend-wise association with pulmonary function, which is significant in patients without MetS or airflow limitation.

Unsupervised clustering identified clinically relevant metabolic syndrome endotypes in UK and Taiwan Biobanks

Metabolic syndrome (MetS) is a collection of cardiovascular risk factors; however, the high prevalence and heterogeneity impede effective clinical management. We conducted unsupervised clustering on individuals from UK Biobank to reveal endotypes. Five MetS subgroups were identified: Cluster 1 (C1): non-descriptive, Cluster 2 (C2): hypertensive, Cluster 3 (C3): obese, Cluster 4 (C4): lipodystrophy-like, and Cluster 5 (C5): hyperglycemic. For all of the endotypes, we identified the corresponding cardiometabolic traits and their associations with clinical outcomes. Genome-wide association studies (GWASs) were conducted to identify associated genotypic traits. We then determined endotype-specific genotypic traits and constructed polygenic risk score (PRS) models specific to each endotype. GWAS of each MetS clusters revealed different genotypic traits. C1 GWAS revealed novel findings of TRIM63, MYBPC3, MYLPF, and RAPSN. Intriguingly, C1, C3, and C4 were associated with genes highly expressed in brain tissues. MetS clusters with comparable phenotypic and genotypic traits were identified in Taiwan Biobank.

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.

Detecting gene-environment interactions from multiple continuous traits

MOTIVATION

Genetic variants present differential effects on humans according to various environmental exposures, the so-called "gene-environment interactions" (GxE). Many diseases can be diagnosed with multiple traits, such as obesity, diabetes, and dyslipidemia. I developed a multivariate scale test (MST) for detecting the GxE of a disease with several continuous traits. Given a significant MST result, I continued to search for which trait and which E enriched the GxE signals. Simulation studies were performed to compare MST with the univariate scale test (UST).

RESULTS

MST can gain more power than UST because of (1) integrating more traits with GxE information and (2) the less harsh penalty on multiple testing. However, if only few traits account for GxE, MST may lose power due to aggregating non-informative traits into the test statistic. As an example, MST was applied to a discovery set of 93 708 Taiwan Biobank (TWB) individuals and a replication set of 25 200 TWB individuals. From among 2 570 487 SNPs with minor allele frequencies ≥5%, MST identified 18 independent variance quantitative trait loci (P < 2.4E-9 in the discovery cohort and P < 2.8E-5 in the replication cohort) and 41 GxE signals (P < .00027) based on eight trait domains (including 29 traits).

AVAILABILITY AND IMPLEMENTATION

https://github.com/WanYuLin/Multivariate-scale-test-MST.

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.

Screening embryos for polygenic disease risk: a review of epidemiological, clinical, and ethical considerations

BACKGROUND

The genetic composition of embryos generated by in vitro fertilization (IVF) can be examined with preimplantation genetic testing (PGT). Until recently, PGT was limited to detecting single-gene, high-risk pathogenic variants, large structural variants, and aneuploidy. Recent advances have made genome-wide genotyping of IVF embryos feasible and affordable, raising the possibility of screening embryos for their risk of polygenic diseases such as breast cancer, hypertension, diabetes, or schizophrenia. Despite a heated debate around this new technology, called polygenic embryo screening (PES; also PGT-P), it is already available to IVF patients in some countries. Several articles have studied epidemiological, clinical, and ethical perspectives on PES; however, a comprehensive, principled review of this emerging field is missing.

OBJECTIVE AND RATIONALE

This review has four main goals. First, given the interdisciplinary nature of PES studies, we aim to provide a self-contained educational background about PES to reproductive specialists interested in the subject. Second, we provide a comprehensive and critical review of arguments for and against the introduction of PES, crystallizing and prioritizing the key issues. We also cover the attitudes of IVF patients, clinicians, and the public towards PES. Third, we distinguish between possible future groups of PES patients, highlighting the benefits and harms pertaining to each group. Finally, our review, which is supported by ESHRE, is intended to aid healthcare professionals and policymakers in decision-making regarding whether to introduce PES in the clinic, and if so, how, and to whom.

SEARCH METHODS

We searched for PubMed-indexed articles published between 1/1/2003 and 1/3/2024 using the terms 'polygenic embryo screening', 'polygenic preimplantation', and 'PGT-P'. We limited the review to primary research papers in English whose main focus was PES for medical conditions. We also included papers that did not appear in the search but were deemed relevant.

OUTCOMES

The main theoretical benefit of PES is a reduction in lifetime polygenic disease risk for children born after screening. The magnitude of the risk reduction has been predicted based on statistical modelling, simulations, and sibling pair analyses. Results based on all methods suggest that under the best-case scenario, large relative risk reductions are possible for one or more diseases. However, as these models abstract several practical limitations, the realized benefits may be smaller, particularly due to a limited number of embryos and unclear future accuracy of the risk estimates. PES may negatively impact patients and their future children, as well as society. The main personal harms are an unindicated IVF treatment, a possible reduction in IVF success rates, and patient confusion, incomplete counselling, and choice overload. The main possible societal harms include discarded embryos, an increasing demand for 'designer babies', overemphasis of the genetic determinants of disease, unequal access, and lower utility in people of non-European ancestries. Benefits and harms will vary across the main potential patient groups, comprising patients already requiring IVF, fertile people with a history of a severe polygenic disease, and fertile healthy people. In the United States, the attitudes of IVF patients and the public towards PES seem positive, while healthcare professionals are cautious, sceptical about clinical utility, and concerned about patient counselling.

WIDER IMPLICATIONS

The theoretical potential of PES to reduce risk across multiple polygenic diseases requires further research into its benefits and harms. Given the large number of practical limitations and possible harms, particularly unnecessary IVF treatments and discarded viable embryos, PES should be offered only within a research context before further clarity is achieved regarding its balance of benefits and harms. The gap in attitudes between healthcare professionals and the public needs to be narrowed by expanding public and patient education and providing resources for informative and unbiased genetic counselling.

Impact of polygenic risk score for triglyceride trajectory and diabetic complications in subjects with type 2 diabetes based on large electronic medical record data from Taiwan: a case control study

BACKGROUND

The prevalence of diabetic dyslipidemia has gradually increased worldwide and individuals with hypertriglyceridemia often have a high polygenic burden of triglyceride (TG)-increasing variants. However, the contribution of genetic variants to dyslipidemia in patients with type 2 diabetes (T2D) remains limited. Therefore, in this study, we aimed to investigate the genetic characteristics of longitudinal changes in TG levels among patients with T2D and summarize the genetic effects of polygenic risk score (PRS) on TG trajectory and risk of diabetic complications.

METHODS

We conducted a case-control study. A total of 11,312 patients with T2D with longitudinal TG and genetic data were identified from a large hospital database in Taiwan. We then performed a genome-wide association study and calculated the relative PRS.

RESULTS

In total, 21 single-nucleotide polymorphisms (SNPs) related to TG trajectory were identified and yielded an area under the receiver operating characteristic curve (ROC) of 0.712 for high TG trajectory risk among Taiwanese patients with T2D. A cumulative genetic effect was observed for high TG trajectory, even when considering the adherence of a lipid-lowering agent in stratified analysis. An increased PRS increases high TG trajectory risk in a logistic regression model (odds ratio = 1.55; 95% confidence interval [CI] = 1.31-1.83 in the validation cohort). The TG-specific PRS was associated with the risk of diabetic microvascular complications, including diabetic retinopathy and nephropathy (with hazard ratios of 1.11 [95% CI = 1.01-1.21, P = 0.027] and 1.05 [95% CI = 1.01-1.1, P = 0.018], respectively).

CONCLUSIONS

This study may contribute to the identification of patients with T2D who are at risk of abnormal TG levels and diabetic microvascular complications using polygenic information.

Polygenic risk score predicting susceptibility and outcome of benign prostatic hyperplasia in the Han Chinese

BACKGROUND

Given the high prevalence of BPH among elderly men, pinpointing those at elevated risk can aid in early intervention and effective management. This study aimed to explore that polygenic risk score (PRS) is effective in predicting benign prostatic hyperplasia (BPH) incidence, prognosis and risk of operation in Han Chinese.

METHODS

A retrospective cohort study included 12,474 male participants (6,237 with BPH and 6,237 non-BPH controls) from the Taiwan Precision Medicine Initiative (TPMI). Genotyping was performed using the Affymetrix Genome-Wide TWB 2.0 SNP Array. PRS was calculated using PGS001865, comprising 1,712 single nucleotide polymorphisms. Logistic regression models assessed the association between PRS and BPH incidence, adjusting for age and prostate-specific antigen (PSA) levels. The study also examined the relationship between PSA, prostate volume, and response to 5-α-reductase inhibitor (5ARI) treatment, as well as the association between PRS and the risk of TURP.

RESULTS

Individuals in the highest PRS quartile (Q4) had a significantly higher risk of BPH compared to the lowest quartile (Q1) (OR = 1.51, 95% CI = 1.274-1.783, p < 0.0001), after adjusting for PSA level. The Q4 group exhibited larger prostate volumes and a smaller volume reduction after 5ARI treatment. The Q1 group had a lower cumulative TURP probability at 3, 5, and 10 years compared to the Q4 group. PRS Q4 was an independent risk factor for TURP.

CONCLUSIONS

In this Han Chinese cohort, higher PRS was associated with an increased susceptibility to BPH, larger prostate volumes, poorer response to 5ARI treatment, and a higher risk of TURP. Larger prospective studies with longer follow-up are warranted to further validate these findings.

Genetics in neovascular age-related macular degeneration susceptibility and treatment response to anti-VEGF intravitreal injection: A case series study

BACKGROUND

To identify genotypes associated with neovascular age-related macular degeneration (nAMD) and investigate the associations between genotype variations and anti-vascular endothelial growth factor (VEGF) treatment response.

METHODS

This observational, retrospective, case series study enrolled patients diagnosed with nAMD who received anti-VEGF treatment in National Taiwan University Hospital with at least one-year follow-up between 2012 and 2020. A genome-wide association study (GWAS) was conducted on enrolled patients and controls. Correlations between the genotypes identified from GWAS and the treatment response of functional/anatomical biomarkers, including visual acuity (VA), presence of intraretinal or subretinal fluid (SRF), serous or fibrovascular pigmented epithelium detachment (PED), and disruption of the ellipsoid zone (EZ), were analysed.

RESULTS

In total, 182 patients with nAMD and 1748 controls were enrolled. GWAS revealed 16 single nucleotide polymorphisms (SNPs) as risk loci for nAMD, including seven loci in CFH and ARMS2/HTRA1 and nine novel loci, including rs117517872 and rs79835234(COPB2-DT), rs7525578(RAP1A), rs2123738(LOC105376755), rs1374879(CNTN3), rs3812692(SAR1A), rs117501587(PRKCA), rs9965945(CNDP1), and rs189769231(MATK). Our study revealed rs800292(CFH), rs11200638(HTRA1), and rs2123738(LOC105376755) correlated with poor treatment response in VA (P = 0.005), SRF (P = 0.044), and fibrovascular PED (P = 0.007), respectively. Rs9965945(CNDP1) was correlated with poor response in disruption of EZ (P = 0.046) and serous PED (P = 0.049).

CONCLUSIONS

Among the 16 SNPs found in the GWAS, four loci-CFH, ARMS2/HTRA1, and two novel loci-were correlated with the susceptibility of nAMD and anatomical/functional responses after anti-VEGF treatment.

Phenomic landscape and pharmacogenomic implications for HLA region in a Taiwan Han Chinese population

BACKGROUND

The human leukocyte antigen (HLA) genes, exhibiting significant genetic diversity, are associated with susceptibility to various clinical diseases and diverse in drug responses. High costs of HLA sequencing and the population-specific architecture of this genetic region necessitate the establishment of a population-specific HLA imputation reference panel. Moreover, there is a lack of understanding about the genetic and phenotypic landscape of HLA variations within the Taiwanese population.

METHODS

We created models for a Taiwanese-specific HLA imputation reference panel. These models were trained with the array genotype data and HLA sequencing data from 845 Taiwanese subjects. HLA imputation was applied for 59,448 Taiwanese subjects to characterize the HLA allele and haplotype frequencies. Additionally, a phenome-wide association study (PheWAS) was conducted to identify the phenotypes associated with HLA variations. The association of the biallelic HLA variants with the binary and quantitative traits were evaluated with additive logistic and linear regression models, respectively. Furthermore, an omnibus test with likelihood-ratio test was applied for each HLA amino acid position in the multiallelic HLA amino acid polymorphisms to compare the difference between a fitted model and a null model following a χ2 distribution of n-1 degree of freedom at a position with n residues. Finally, we estimated the prevalence of adverse drug reactions (ADR)-related HLA alleles in the Taiwanese population.

RESULTS

In this study, the reference panel models displayed remarkable accuracy, with averages of 99.3%, 98.9%, and 99.1% for 2-, 4-, 6-digit alleles of the eight classical HLA genes, respectively. For PheWAS, a total of 18,136 significant associations with HLA variants across 26 phenotypes are identified (p < 5×10-8), highlighting the pleiotropy feature of the HLA region. Among the independent signals, 15 are novel, including the association of HLA-B pos 138 variation with ankylosing spondylitis (AS), and rs9266290 and rs9266292 with allergy. Through an analysis spanning the entire HLA region, we identified clusters of phenotype correlations. Finally, the carriers of pharmacogenomic related HLA alleles, including HLA-C*01:02 (35.86%), HLA-B*58:01 (20.9%), and HLA-B*15:02 (8.38%), were characterized in the Taiwanese general population.

CONCLUSIONS

We successfully delivered the HLA imputation for 59,448 Taiwanese subjects and characterized the genetic and phenotypic landscapes of the HLA variations. In addition, we quantified the estimated prevalence of the ADR-related HLA alleles in the Taiwanese population. The developed HLA imputation reference panel could be used for estimation of population HLA allele frequencies, which can facilitate further studies in the role of HLA variants in a wider range of phenotypes in the population.

Imputation accuracy across global human populations

Genotype imputation is now fundamental for genome-wide association studies but lacks fairness due to the underrepresentation of references from non-European ancestries. The state-of-the-art imputation reference panel released by the Trans-Omics for Precision Medicine (TOPMed) initiative improved the imputation of admixed African-ancestry and Hispanic/Latino samples, but imputation for populations primarily residing outside of North America may still fall short in performance due to persisting underrepresentation. To illustrate this point, we imputed the genotypes of over 43,000 individuals across 123 populations around the world and identified numerous populations where imputation accuracy paled in comparison to that of European-ancestry populations. For instance, the mean imputation r-squared (Rsq) for variants with minor allele frequencies between 1% and 5% in Saudi Arabians (n = 1,061), Vietnamese (n = 1,264), Thai (n = 2,435), and Papua New Guineans (n = 776) were 0.79, 0.78, 0.76, and 0.62, respectively, compared to 0.90-0.93 for comparable European populations matched in sample size and SNP array content. Outside of Africa and Latin America, Rsq appeared to decrease as genetic distances to European-ancestry reference increased, as predicted. Using sequencing data as ground truth, we also showed that Rsq may over-estimate imputation accuracy for non-European populations more than European populations, suggesting further disparity in accuracy between populations. Using 1,496 sequenced individuals from Taiwan Biobank as a second reference panel to TOPMed, we also assessed a strategy to improve imputation for non-European populations with meta-imputation, but this design did not improve accuracy across frequency spectra. Taken together, our analyses suggest that we must ultimately strive to increase diversity and size to promote equity within genetics research.

Impact of the trans-ancestry polygenic risk score on type 2 diabetes risk, onset age and progression among population in Taiwan

Type 2 diabetes (T2D) prevalence in adults at a younger age has increased but the disease status may go unnoticed. This study aimed to determine whether the onset age and subsequent diabetic complications can be attributed to the polygenic architecture of T2D in the Taiwan Han population. A total of 9,627 cases with T2D and 85,606 controls from the Taiwan Biobank were enrolled. Three diabetic polygenic risk scores (PRSs), PRS_EAS and PRS_EUR, and a trans-ancestry PRS (PRS_META), calculated using summary statistic from East Asian and European populations. The onset age was identified by linking to the National Taiwan Insurance Research Database, and the incidence of different diabetic complications during follow-up was recorded. PRS_META (7.4%) explained a higher variation for T2D status. And the higher percentile of PRS is also correlated with higher percentage of T2D family history and prediabetes status. More, the PRS was negatively associated with onset age (β = -0.91 yr), and this was more evident among males (β = -1.11 vs. -0.76 for males and females, respectively). The hazard ratio of diabetic retinopathy (DR) and diabetic foot were significantly associated with PRS_EAS and PRS_META, respectively. However, the PRS was not associated with other diabetic complications, including diabetic nephropathy, cardiovascular disease, and hypertension. Our findings indicated that diabetic PRS which combined susceptibility variants from cross-population could be used as a tool for early screening of T2D, especially for high-risk populations, such as individuals with high genetic risk, and may be associated with the risk of complications in subjects with T2D. NEW & NOTEWORTHY Our findings indicated that diabetic polygenic risk score (PRS) which combined susceptibility variants from Asian and European population affect the onset age of type 2 diabetes (T2D) and could be used as a tool for early screening of T2D, especially for individuals with high genetic risk, and may be associated with the risk of diabetic complications among people in Taiwan.

Predicting Risks of Dry Eye Disease Development Using a Genome-Wide Polygenic Risk Score Model

Purpose

The purpose of this study was to conduct a large-scale genome-wide association study (GWAS) and construct a polygenic risk score (PRS) for risk stratification in patients with dry eye disease (DED) using the Taiwan Biobank (TWB) databases.

Methods

This retrospective case-control study involved 40,112 subjects of Han Chinese ancestry, sourced from the publicly available TWB. Cases were patients with DED (n = 14,185), and controls were individuals without DED (n = 25,927). The patients with DED were further divided into 8072 young (<60 years old) and 6113 old participants (≥60 years old). Using PLINK (version 1.9) software, quality control was carried out, followed by logistic regression analysis with adjustments for sex, age, body mass index, depression, and manic episodes as covariates. We also built PRS prediction models using the standard clumping and thresholding method and evaluated their performance (area under the curve [AUC]) through five-fold cross-validation.

Results

Eleven independent risk loci were identified for these patients with DED at the genome-wide significance levels, including DNAJB6, MAML3, LINC02267, DCHS1, SIRPB3P, HULC, MUC16, GAS2L3, and ZFPM2. Among these, MUC16 encodes mucin family protein. The PRS model incorporated 932 and 740 genetic loci for young and old populations, respectively. A higher PRS score indicated a greater DED risk, with the top 5% of PRS individuals having a 10-fold higher risk. After integrating these covariates into the PRS model, the area under the receiver operating curve (AUROC) increased from 0.509 and 0.537 to 0.600 and 0.648 for young and old populations, respectively, demonstrating the genetic-environmental interaction.

Conclusions

Our study prompts potential candidates for the mechanism of DED and paves the way for more personalized medication in the future.

Translational Relevance

Our study identified genes related to DED and constructed a PRS model to improve DED prediction.

GBA1 as a risk gene for osteoporosis in the specific populations and its role in the development of Gaucher disease

BACKGROUND

Osteoporosis and its primary complication, fragility fractures, contribute to substantial global morbidity and mortality. Gaucher disease (GD) is caused by glucocerebrosidase (GBA1) deficiency, leading to skeletal complications. This study aimed to investigate the impact of the GBA1 gene on osteoporosis progression in GD patients and the specific populations.

METHODS

We selected 8115 patients with osteoporosis (T-score ≤ - 2.5) and 55,942 healthy individuals (T-score > - 1) from a clinical database (N = 95,223). Monocytes from GD patients were evaluated in relation to endoplasmic reticulum (ER) stress, inflammasome activation, and osteoclastogenesis. An in vitro model of GD patient's cells treated with adeno-associated virus 9 (AAV9)-GBA1 to assess GBA1 enzyme activity, chitotriosidase activity, ER stress, and osteoclast differentiation. Longitudinal dual-energy X-ray absorptiometry (DXA) data tracking bone density in patients with Gaucher disease (GD) undergoing enzyme replacement therapy (ERT) over an extended period.

RESULTS

The GBA1 gene variant rs11264345 was significantly associated [P < 0.002, Odds Ratio (OR) = 1.06] with an increased risk of bone disease. Upregulation of Calnexin, NOD-, LRR- and pyrin domain-containing protein 3 (NLRP3) and Apoptosis-associated speck-like protein containing a C-terminal caspase recruitment domain (ASC) was positively associated with osteoclastogenesis in patients with GD. In vitro AAV9-GBA1 treatment of GD patient cells led to enhanced GBA1 enzyme activity, reduced chitotriosidase activity, diminished ER stress, and decreased osteoclast differentiation. Long-term bone density data suggests that initiating ERT earlier in GD leads to greater improvements in bone density.

CONCLUSIONS

Elevated ER stress and inflammasome activation are indicative of osteoporosis development, suggesting the need for clinical monitoring of patients with GD. Furthermore, disease-associated variant in the GBA1 gene may constitute a risk factor predisposing specific populations to osteoporosis.

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.

Secondhand smoke increases the risk of developing chronic obstructive pulmonary disease

Smoking is the most important risk factor for chronic obstructive pulmonary disease (COPD), however evidence from large-scale studies on whether secondhand smoke (SHS) increases the risk of COPD is still lacking. We conducted this large longitudinal study to investigate the association between SHS and the development of COPD. This is a longitudinal study. Data on 6519 subjects who were never-smokers, had no history of COPD, and had complete lung function records were extracted from the Taiwan Biobank. They were divided into two groups according to SHS exposure: no exposure and exposure groups. Data were collected when participants enrolled in the study and during regular follow-up. Cox proportional hazards regression models were used to estimate the relative risk (RR) and 95% confidence interval (CI) for the association between SHS and the risk of developing COPD. At 48 months of follow-up, 260 (4%) participants in the no exposure group and 34 (7%) participants in the exposure group developed COPD. The RR of incident COPD development was significantly higher in the exposure group than that in the no exposure group after adjusting for confounders (RR = 1.49; 95% CI 1.04 to 2.14; P value = 0.031). There is a dose-response relationship between the duration of exposure to SHS and the risk of incident COPD, which demonstrates that an additional hour of exposure to SHS per week was associated with a 1.03-fold increased likelihood of developing COPD after adjusting for confounders (RR = 1.03; 95% CI 1.00 to 1.05; P value = 0.027). SHS exposure contributes to the development of COPD. This finding can help raise awareness of the harms of SHS and provide a reference for formulating anti-smoking policies.

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.

Twelve-year review of galactosemia newborn screening in Taiwan: Evolving methods and insights

Background

Galactosemia was introduced into Taiwan's routine newborn screening (NBS) program in 1985. This study presents a 12-year experience, emphasizing disease diagnosis and screening performance.

Method

NBS for galactosemia utilized dried blood spot samples taken 48-72 h post-delivery, with total galactose (TGal) level as the primary marker. Newborns with critical TGal levels were referred immediately, while those with borderline TGal underwent a recall test. GALT activity measurement was applied simultaneously as the second-tier marker. Further confirmatory tests, such as whole exome sequencing (WES), were conducted upon referral.

Results

From January 1st, 2011, to December 31st, 2022, 51 cases were identified from 817,906 newborns. Of these, nine individuals had persistently elevated TGal. Diagnoses included one case of GALT deficiency, one of GALM deficiency, and seven of GALE deficiencies. Notably, the classic galactosemia patient (GALT deficiency) presented with extreme high TGal and was referred to the hospital for diet management immediately. All affected patients were instructed to adopt a galactose-restricted diet. By the median age of 2.5 years, all exhibited normal development and liver function.

Conclusion

The incidence of classical galactosemia and its variants is extremely low in Taiwan. Incorporating WES into NBS has improved our ability to detect various galactosemia forms, enriching our understanding of the genetic underpinnings. While these newly discovered forms often present with milder initial elevations in TGal, specific biochemical investigations and regular monitoring are essential to understanding the long-term implications and outcomes.

Ethnically unique disease burden and limitations of current expanded carrier screening panels

OBJECTIVES

The purpose of the study is to identify the recessive diseases currently affecting real-world pediatric patients in Taiwan, and whether current extended carrier screening panels have the coverage and detective power to identify the pathogenic variants in the carrier parents.

METHODS

A total of 132 trio-samples were collected from May 2017 to March 2022. The participants were parents of pediatric intensive care unit patients who were critically ill or infants with abnormal newborn screening results. A retrospective carrier screening scheme was applied to analyze only the carrier status of pathogenic or likely pathogenic recessive variants resulting in diseases in their children. The recessive disorders diagnosed in our cohort were compared with the gene content in commercial panels.

RESULTS

Mutations in COQ4, PEX1, OTC, and IKBKG were the most frequently identified. In the parents of 44 children with confirmed diagnoses of recessive diseases, 47 (53.40%) screened positive for being the carriers of the same recessive disorders diagnosed in their children. The commercial panels covered 35.13% to 54.05% of the disorders diagnosed in this cohort.

CONCLUSION

Clinicians and genetic counselors should be aware of the limitations of current extended carrier screening and interpret negative screening results with caution. Future panels should also consider genes with ethnically unique mutations such as pathogenic variants of the COQ4 gene in the East Asian population.

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.

The characterization and comorbidities of heterozygous Bardet-Biedl syndrome carriers

Introduction: Bardet-Biedl syndrome (BBS) is a rare autosomal recessive disorder with clinical features of retinal dystrophy, obesity, postaxial polydactyly, renal anomalies, learning disabilities, hypogonadism, and genitourinary abnormalities. Nevertheless, previous studies on the phenotypic traits of BBS heterozygous carriers have generated inconclusive results. The aim of our study was to investigate the impact of BBS heterozygosity on carriers when compared to non-carriers within the Taiwanese population. Materials and Methods: This study follows a hospital-based case-control design. We employed the Taiwan Biobank version 2 (TWBv2) array to identify three specific loci associated with BBS (rs773862084, rs567573386, and rs199910690). In total, 716 patients were included in the case group, and they were compared to a control group of 2,864 patients who lacked BBS alleles. The control group was selected through gender and age matching at a ratio of 1:4. The association between BBS-related loci and comorbidity was assessed using logistic regression models. Results: We found that BBS heterozygous carriers exhibited a significant association with elevated BMI levels, especially the variant rs199910690 in MKS1 (p=0.0037). The prevalence of comorbidities in the carriers' group was not higher than that in the non-carriers' group. Besides, the average values of the biochemistry data showed no significant differences, except for creatinine level. Furthermore, we conducted a BMI-based analysis to identify specific risk factors for chronic kidney disease (CKD). Our findings revealed that individuals carrying the CA/AA genotype of the BBS2 rs773862084 variant or the CT/TT genotype of the MKS1 rs199910690 variant showed a reduced risk of developing CKD, irrespective of their BMI levels. When stratified by BMI level, obese males with the MKS1 rs199910690 variant and obese females with the BBS2 rs773862084 variant exhibited a negative association with CKD development. Conclusion: We found that aside from the association with overweight and obesity, heterozygous BBS mutations did not appear to increase the predisposition of individuals to comorbidities and metabolic diseases. To gain a more comprehensive understanding of the genetic susceptibility associated with Bardet-Biedl Syndrome (BBS), further research is warranted.

Blastocyst telomere length predicts successful implantation after frozen-thawed embryo transfer

STUDY QUESTION

Do embryos with longer telomere length (TL) at the blastocyst stage have a higher capacity to survive after frozen-thawed embryo transfer (FET)?

SUMMARY ANSWER

Digitally estimated TL using low-pass whole genome sequencing (WGS) data from the preimplantation genetic testing for aneuploidy (PGT-A) process demonstrates that blastocyst TL is the most essential factor associated with likelihood of implantation.

WHAT IS KNOWN ALREADY

The lifetime TL is established in the early cleavage cycles following fertilization through a recombination-based lengthening mechanism and starts erosion beyond the blastocyst stage. In addition, a telomerase-mediated slow erosion of TL in human fetuses has been observed from a gestational age of 6-11 weeks. Finally, an abnormal shortening of telomeres is likely involved in embryo loss during early development.

STUDY DESIGN SIZE DURATION

Blastocyst samples were obtained from patients who underwent PGT-A and FET in an IVF center from March 2015 to May 2018. Digitally estimated mitochondrial copy number (mtCN) and TL were used to study associations with the implantation potential of each embryo.

PARTICIPANTS/MATERIALS SETTING AND METHODS

In total, 965 blastocysts from 232 cycles (164 patients) were available to investigate the biological and clinical relevance of TL. A WGS-based workflow was applied to determine the ploidy of each embryo. Data from low-pass WGS-PGT-A were used to estimate the mtCN and TL for each embryo. Single-variant and multi-variant logistic regression, decision tree, and random forest models were applied to study various factors in association with the implantation potential of each embryo.

MAIN RESULTS AND THE ROLE OF CHANCE

Of the 965 blastocysts originally available, only 216 underwent FET. While mtCN from the transferred embryos is significantly associated with the ploidy call of each embryo, mtCN has no role in impacting IVF outcomes after an embryo transfer in these women. The results indicate that mtCN is a marker of embryo aneuploidy. On the other hand, digitally estimated TL is the most prominent univariant factor and showed a significant positive association with pregnancy outcomes (P < 0.01, odds ratio 79.1). We combined several maternal and embryo parameters to study the joint effects on successful implantation. The machine learning models, namely decision tree and random forest, were trained and yielded classification accuracy of 0.82 and 0.91, respectively. Taken together, these results support the vital role of TL in governing implantation potential, perhaps through the ability to control embryo survival after transfer.

LIMITATIONS REASONS FOR CAUTION

The small sample size limits our study as only 216 blastocysts were transferred. The number was further reduced to 153 blastocysts, where pregnancy outcomes could be accurately traced. The other limitation of this study is that all data were collected from a single IVF center. The uniform and controlled operation of IVF cycles in a single center may cause selection bias.

WIDER IMPLICATIONS OF THE FINDINGS

We present novel findings to show that digitally estimated TL at the blastocyst stage is a predictor of pregnancy capacity after a FET cycle. As elective single-embryo transfer has become the mainstream direction in reproductive medicine, prioritizing embryos based on their implantation potential is crucial for clinical infertility treatment in order to reduce twin pregnancy rate and the time to pregnancy in an IVF center. The AI-powered, random forest prediction model established in this study thus provides a way to improve clinical practice and optimize the chances for people with fertility problems to achieve parenthood.

STUDY FUNDING/COMPETING INTERESTS

This study was supported by a grant from the National Science and Technology Council, Taiwan (MOST 108-2321-B-006-013 -). There were no competing interests.

TRIAL REGISTRATION NUMBER

N/A.

Host-Pathogen Interactions in K. pneumoniae Urinary Tract Infections: Investigating Genetic Risk Factors in the Taiwanese Population

BACKGROUND

Klebsiella pneumoniae (K. pneumoniae) urinary tract infections pose a significant challenge in Taiwan. The significance of this issue arises because of the growing concerns about the antibiotic resistance of K. pneumoniae. Therefore, this study aimed to uncover potential genomic risk factors in Taiwanese patients with K. pneumoniae urinary tract infections through genome-wide association studies (GWAS).

METHODS

Genotyping data are obtained from participants with a history of urinary tract infections enrolled at the Tri-Service General Hospital as part of the Taiwan Precision Medicine Initiative (TPMI). A case-control study employing GWAS is designed to detect potential susceptibility single-nucleotide polymorphisms (SNPs) in patients with K. pneumoniae-related urinary tract infections. The associated genes are determined using a genome browser, and their expression profiles are validated via the GTEx database. The GO, Reactome, DisGeNET, and MalaCards databases are also consulted to determine further connections between biological functions, molecular pathways, and associated diseases between these genes.

RESULTS

The results identified 11 genetic variants with higher odds ratios compared to controls. These variants are implicated in processes such as adhesion, protein depolymerization, Ca2+-activated potassium channels, SUMOylation, and protein ubiquitination, which could potentially influence the host immune response.

CONCLUSIONS

This study implies that certain risk variants may be linked to K. pneumoniae infections by affecting diverse molecular functions that can potentially impact host immunity. Additional research and follow-up studies are necessary to elucidate the influence of these risk variants on infectious diseases and develop targeted interventions for mitigating the spread of K. pneumoniae urinary tract infections.

Polygenic Risk Score in Predicting Esophageal, Oropharyngeal, and Hypopharynx Cancer Risk among Taiwanese Population

Esophageal cancer shares strong associations with oropharyngeal and hypopharyngeal cancers, primarily due to shared risk factors like excessive tobacco and alcohol use. This retrospective study at Taichung Veterans General Hospital involved 54,692 participants, including 385 with squamous cell carcinoma (SCC) of the esophagus, oropharynx, or hypopharynx. Using a polygenic risk score (PRS) derived from 8353 single-nucleotide polymorphisms, researchers aimed to assess its correlation with cancer incidence and prognosis. The study found a 1.83-fold higher risk of esophageal, oropharyngeal, and hypopharyngeal SCCs in participants with a high PRS (Q4) compared to the low-PRS group (Q1). Esophageal cancer risk demonstrated a significant positive association with the PRS, as did hypopharyngeal cancer. Clinical parameters and staging showed limited associations with PRS quartiles, and the PRS did not significantly impact recurrence or mortality rates. The research highlighted that a higher PRS is linked to increased susceptibility to esophageal and hypopharyngeal cancer. Notably, a specific polygenic risk score, PGS001087, exhibited a discernible association with SCC risk, particularly in specific subtypes and advanced disease stages. However, it was not significantly linked to clinical cancer staging, emphasizing the multifactorial nature of cancer development. This hospital study reveals that a higher PRS correlates with increased susceptibility to esophageal and hypopharyngeal cancers. Notably, PGS001087 shows a discernible association with SCC risk in specific subtypes and advanced stages, although not significantly linked to clinical cancer staging. These findings enhance our understanding of genetic factors in upper aerodigestive tract cancers, particularly esophageal SCC, guiding future research and risk assessment strategies.

More than a decade of genetic research on the Denisovans

Denisovans, a group of now extinct humans who lived in Eastern Eurasia in the Middle and Late Pleistocene, were first identified from DNA sequences just over a decade ago. Only ten fragmentary remains from two sites have been attributed to Denisovans based entirely on molecular information. Nevertheless, there has been great interest in using genetic data to understand Denisovans and their place in human history. From the reconstruction of a single high-quality genome, it has been possible to infer their population history, including events of admixture with other human groups. Additionally, the identification of Denisovan DNA in the genomes of present-day individuals has provided insights into the timing and routes of dispersal of ancient modern humans into Asia and Oceania, as well as the contributions of archaic DNA to the physiology of present-day people. In this Review, we synthesize more than a decade of research on Denisovans, reconcile controversies and summarize insights into their population history and phenotype. We also highlight how our growing knowledge about Denisovans has provided insights into our own evolutionary history.

Prosaposin variants in sporadic, familial, and early-onset Parkinson's disease: a Taiwanese case-control study and meta-analysis

Polymorphisms in the PSAP gene, which encodes prosaposin and is involved in the lysosomal function, yielded conflicting results regarding the association with Parkinson's disease (PD). Therefore, this study aims to investigate the role of PSAP in familial PD (FPD), early onset PD (EOPD) with age at onset before 50 years old, and sporadic PD (SPD) among Taiwanese population, and summarize relevant studies via meta-analysis. By sequencing exon 1 to 14 in 183 FPD and 219 EOPD, two novel exonic variants were found in EOPD, including p.A146E (c.437C > A) on exon 5 and p.Y248C (c.743A > G) on exon 7. Furthermore, four previously reported intronic variants (rs142614739/rs74733861), rs749823, rs4747203 and rs885828) in intron 11 and 12 were analyzed in 485 SPD and 712 in-hospital controls, in addition to the aforementioned FPD and EOPD groups. The adjusted odd ratios (ORs) by age and sex, only rs142614739 was significantly associated with higher risk of EOPD (OR = 1.85, 95% CI = 1.33-2.58). The risk effect was further confirmed by the meta-analysis of the association between rs142614739 and the risk of PD in both common effect (OR = 1.29, 95% CI = 1.11-1.50) and random effect (OR = 1.29, 95% CI = 1.11-1.50). Our findings suggest that the PSAP rs142614739 variant is associated with the risk of EOPD. Further functional studies are warranted to elucidate the biochemical mechanisms.

Machine learning approaches to identify systemic lupus erythematosus in anti-nuclear antibody-positive patients using genomic data and electronic health records

BACKGROUND

Although the 2019 EULAR/ACR classification criteria for systemic lupus erythematosus (SLE) has required at least a positive anti-nuclear antibody (ANA) titer (≥ 1:80), it remains challenging for clinicians to identify patients with SLE. This study aimed to develop a machine learning (ML) approach to assist in the detection of SLE patients using genomic data and electronic health records.

METHODS

Participants with a positive ANA (≥ 1:80) were enrolled from the Taiwan Precision Medicine Initiative cohort. The Taiwan Biobank version 2 array was used to detect single nucleotide polymorphism (SNP) data. Six ML models, Logistic Regression, Random Forest (RF), Support Vector Machine, Light Gradient Boosting Machine, Gradient Tree Boosting, and Extreme Gradient Boosting (XGB), were used to identify SLE patients. The importance of the clinical and genetic features was determined by Shapley Additive Explanation (SHAP) values. A logistic regression model was applied to identify genetic variations associated with SLE in the subset of patients with an ANA equal to or exceeding 1:640.

RESULTS

A total of 946 SLE and 1,892 non-SLE controls were included in this analysis. Among the six ML models, RF and XGB demonstrated superior performance in the differentiation of SLE from non-SLE. The leading features in the SHAP diagram were anti-double strand DNA antibodies, ANA titers, AC4 ANA pattern, polygenic risk scores, complement levels, and SNPs. Additionally, in the subgroup with a high ANA titer (≥ 1:640), six SNPs positively associated with SLE and five SNPs negatively correlated with SLE were discovered.

CONCLUSIONS

ML approaches offer the potential to assist in diagnosing SLE and uncovering novel SNPs in a group of patients with autoimmunity.

The largest genome-wide association study for breast cancer in Taiwanese Han population

PURPOSE

Breast cancer is a molecularly heterogeneous disease, and multiple genetic variants contribute to its development and prognosis. Most of previous genome-wide association studies (GWASs) and polygenic risk scores (PRSs) analyses focused on studying breast cancers of Caucasian populations, which may not be applicable to other population. Therefore, we conducted the largest breast cancer cohort of Taiwanese population to fill in the knowledge gap.

METHODS

A total of 152,534 Participants recruited by China Medical University Hospital between 2003 and 2019 were filtered by several patient selection criteria and GWAS quality control steps, resulting in the inclusion of 2496 cases and 9984 controls for this study. We then conducted GWAS for all breast cancers and PRS analyses for all breast cancers and the four breast cancer subtypes, including luminal A, luminal B, basal-like, and HER2-enriched.

RESULTS

The GWAS analyses identified 113 SNPs, 50 of which were novel. The PRS models for all breast cancers and the luminal A subtype showed positively correlated trends between the PRS and the risk of developing breast cancer. The odds ratios (95% confidence intervals) for the groups with the highest PRS in all breast cancers and the luminal A subtype were 5.33 (3.79-7.66) and 3.55 (2.13-6.14), respectively.

CONCLUSION

In summary, we explored the association of genetic variants with breast cancer in the largest Taiwanese cohort and developed two PRS models that can predict the risk of developing any breast cancer and the luminal A subtype in Taiwanese women.

Modifiable vascular risk factors contribute to stroke in 1080 NOTCH3 R544C carriers in Taiwan Biobank

BACKGROUND AND AIM

Previous studies have suggested cardiovascular risk factors increase the risk of not only common sporadic stroke but also of stroke in patients with monogenic stroke disorders including CADASIL. We investigated the effects of the NOTCH3 Arg544Cys (R544C) variant and associated vascular risk factors on stroke in the Taiwanese population.

METHODS

This study was conducted using data from the Taiwan Biobank, consisting of at least 130,000 Han Chinese participants. The genotype was derived from customized genome-wide arrays for 650,000 to 750,000 single-nucleotide polymorphisms (SNPs). Individuals with NOTCH3 R544C were subsequently matched with noncarriers based on the propensity score at a 1:10 ratio by demographic and cardiovascular risk factors. The odds ratio (OR) for stroke or other phenotypes in NOTCH3 R544C carriers and matched noncarriers was then calculated. Univariate and multivariate regression analyses were performed on cardiovascular risk factors in NOTCH3 R544C carriers with and without stroke. The polygenic risk score (PRS) model, adopted from the UK Biobank, was then applied to evaluate the role of NOTCH3 R544C in stroke.

RESULTS

From the 114,282 participants with both genotype and questionnaire results, 1080 (0.95%) harbored the pathogenic NOTCH3 R544C variant. When compared to the matched controls (n = 10,800), the carriers presented with a history of stroke (OR: 2.52, 95% confidence interval (CI) (1.45, 4.37)), dementia (OR: 30.1, 95% CI (3.13, 289.43)), and sibling history of stroke (OR: 2.48, 95% CI (1.85, 3.34)) phenotypes. The risk of stroke increased with every 10-year increase in age (p = 0.006, Cochran-Mantel-Haenszel test). Among NOTCH3 R544C carriers, 16 (1.3%) of the 1080 carriers with a stroke history were older, male, and more likely to have hypertension, diabetes, dyslipidemia, and a family history of stroke. In the stepwise multivariate analysis, hypertension (OR: 11.28, 95% CI (3.54, 43.3)) and diabetes mellitus (OR: 4.10, 95% CI (1.31, 12.4)) were independently associated with stroke. Harboring the NOTCH3 R544C variant in the Taiwan Biobank is comparable with a 6.74 standard deviations increase in individual's polygenic risk score for stroke.

CONCLUSION

While the NOTCH3 R544C variant alone increased the risk of stroke, modifiable vascular risk factors also played a role in the occurrence of stroke in Taiwanese community-dwelling individuals carrying the NOTCH3 variant.

Complete genomic profiles of 1496 Taiwanese reveal curated medical insights

INTRODUCTION

The population of Taiwan has a long history of ethno-cultural evolution. The Taiwanese population was isolated from other large populations such as the European, Han Chinese, and Japanese population. The Taiwan Biobank (TWB) project has built a nationwide database, particularly for personal whole-genome sequence (WGS) to facilitate basic and clinical collaboration nationally and internationally, making it one of the most valuable public datasets of the East Asian population.

OBJECTIVES

This study provides comprehensive medical genomic findings from TWB WGS data, for better characterization of disease susceptibility and the choice of ideal treatment regimens in Taiwanese population.

METHODS

We reanalyzed 1496 WGS using a PrecisionFDA Truth challenge winner method Sentieon DNAscope. Single nucleotide variants (SNV) and small insertions/deletions (INDEL) were benchmarked. We also analyzed pharmacogenomic (PGx) drug-associated alleles, and copy number variants (CNV). Multiple practicing clinicians reviewed and curated the clinically significant variants. Variant annotations can be browsed at TaiwanGenomes (https://genomes.tw).

RESULTS

We found that each participant had an average of 6,870.7 globally novel variants and 75.3% (831/1103) of the participants harbored at least one PharmGKB-selected high evidence level human leukocyte antigen (HLA) risk allele. 54 PharmGKB-reported high-level instances of evidence of Cytochrome P450 variant-drug pairs, with a population frequency of over 13.2%. We also identified 23 variants in the ACMG secondary finding V3 gene list from 25 participants, suggesting that 1.67% (25/1496) of the population is harboring at least one medical actionable variant. Our carrier status analyses suggest that one in 25 couples (3.94%) would risk having offspring with at least one pathogenic variant, which is in line with rates found in Japan and Singapore. For pathogenic CNV, we detected 6.88% and 2.02% carrier rates for alpha thalassemia and spinal muscular atrophy, respectively.

CONCLUSION

Our study highlights the overall medical insights of a complete Taiwanese genomic profile.

A Genome-Wide Association Study of Metabolic Syndrome in the Taiwanese Population

The purpose of this study was to investigate genetic factors associated with metabolic syndrome (MetS) by conducting a large-scale genome-wide association study (GWAS) in Taiwan, addressing the limited data on Asian populations compared to Western populations. Using data from the Taiwan Biobank, comprehensive clinical and genetic information from 107,230 Taiwanese individuals was analyzed. Genotyping data from the TWB1.0 and TWB2.0 chips, including over 650,000 single nucleotide polymorphisms (SNPs), were utilized. Genotype imputation using the 1000 Genomes Project was performed, resulting in more than 9 million SNPs. MetS was defined based on a modified version of the Adult Treatment Panel III criteria. Among all participants (mean age: 50 years), 23% met the MetS definition. GWAS analysis identified 549 SNPs significantly associated with MetS, collectively mapping to 10 genomic risk loci. Notable risk loci included rs1004558, rs3812316, rs326, rs4486200, rs2954038, rs10830963, rs662799, rs62033400, rs183130, and rs34342646. Gene-set analysis revealed 22 associated genes: CETP, LPL, APOA5, SIK3, ZPR1, APOC1, BUD13, MLXIPL, TOMM40, GCK, YKT6, RPS6KB1, FTO, VMP1, TUBD1, BCL7B, C19orf80 (ANGPTL8), SIDT2, SENP7, PAFAH1B2, DOCK6, and FOXA2. This study identified genomic risk loci for MetS in a large Taiwanese population through a comprehensive GWAS approach. These associations provide novel insights into the genetic basis of MetS and hold promise for the potential discovery of clinical biomarkers.

Analysis across Taiwan Biobank, Biobank Japan, and UK Biobank identifies hundreds of novel loci for 36 quantitative traits

Genome-wide association studies (GWASs) have identified tens of thousands of genetic loci associated with human complex traits. However, the majority of GWASs were conducted in individuals of European ancestries. Failure to capture global genetic diversity has limited genomic discovery and has impeded equitable delivery of genomic knowledge to diverse populations. Here we report findings from 102,900 individuals across 36 human quantitative traits in the Taiwan Biobank (TWB), a major biobank effort that broadens the population diversity of genetic studies in East Asia. We identified 968 novel genetic loci, pinpointed novel causal variants through statistical fine-mapping, compared the genetic architecture across TWB, Biobank Japan, and UK Biobank, and evaluated the utility of cross-phenotype, cross-population polygenic risk scores in disease risk prediction. These results demonstrated the potential to advance discovery through diversifying GWAS populations and provided insights into the common genetic basis of human complex traits in East Asia.

Utilize polygenic risk score to enhance fracture risk estimation and improve the performance of FRAX in patients with osteoporosis

This study examined the use of polygenic risk scores (PGS) in combination with the Fracture Risk Assessment Tool (FRAX) to enhance fragility fractures risk estimation in osteoporosis patients. Analyzing data from over 57,000 participants, PGS improved fracture risk estimation, especially for individuals with intermediate to low risks, allowing personalized preventive strategies.

INTRODUCTION

Osteoporosis and fragility fractures are multifactorial, with contributions from both clinical and genetic determinants. However, whether using polygenic risk scores (PGS) may enhance the risk estimation of osteoporotic fracture in addition to Fracture Risk Assessment Tool (FRAX) remains unknown. This study investigated the collective association of PGS and FRAX with fragility fracture.

METHODS

We conducted a cohort study from the Taiwan Precision Medicine Initiative (TPMI) at Taichung Veterans General Hospital, Taiwan. Genotyping was performed to compute PGS associated with bone mineral density (BMD). Phenome-wide association studies were executed to pinpoint phenotypes correlated with the PGS. Logistic regression analysis was conducted to ascertain factors associated with osteoporotic fractures.

RESULTS

Among all 57,257 TPMI participants, 3744 (904 men and 2840 women, with a mean age of 66.7) individuals had BMD testing, with 540 (14.42%) presenting with fractures. The 3744 individuals who underwent BMD testing were categorized into four quartiles (Q1-Q4) based on PGS; 540 (14.42%) presented with fractures. Individuals with PGS-Q1 exhibited lower BMD, a higher prevalence of major fractures, and elevated FRAX-major and FRAX-hip than those with PGS-Q4. PGS was associated with major fractures after adjusting age, sex, and FRAX scores. Notably, the risk of major fractures (PGS-Q1 vs. Q4) was significantly higher in the subgroups of FRAX-major scores < 10% and 10-20%, but not in participants with a FRAX-major score ≧ 20%.

CONCLUSIONS

Our study highlights the potential of PGS to augment fracture risk estimation in conjunction with FRAX, particularly in individuals with middle to low risks. Incorporating genetic testing could empower physicians to tailor personalized preventive strategies for osteoporosis.

Populational pan-ethnic screening panel enabled by deep whole genome sequencing

Birth defect is a global threat to the public health systems. Mitigating neonatal anomalies is hampered by elusive molecular mechanisms of pathogenic mutations and poor subsequent translation into preventative measures. Applying appropriate strategies in China to promote reproductive health is particularly challenging, as the Chinese population compromises complex genomic diversity due to the inclusion of many ethnic groups with distinct genetic backgrounds. To investigate and evaluate the feasibility of implementing a pan-ethnic screening strategy, and guide future reproductive counselling, high-quality variants associated with autosome recessive (AR) diseases derived from the largest publicly available cohort of the Chinese population were re-analysed using a bottom-up approach. The analyses of gene carrier rates (GCRs) across distinct ethnic groups revealed that substantial heterogeneity existed potentially due to diverse evolutionary selection. The sampling population, sequencing coverage and underlying population structure contributed to the differential variants observed between ChinaMAP and the East Asian group in gnomAD. Beyond characteristics of GCR, potential druggable targets were additionally explored according to genomic features and functional roles of investigated genes, demonstrating that phase separation could be a therapeutic target for autosomal recessive diseases. A further examination of estimated GCR across ethnic groups indicated that most genes shared by at least two populations could be utilised to direct the design of a pan-ethnic screening application once sequencing and interpreting costs become negligible. To this end, a list of autosomal recessive disease genes is proposed based on the prioritised rank of GCR to formulate a tiered screening strategy.

Genotype-phenotype correlation in Taiwanese children with diazoxide-unresponsive congenital hyperinsulinism

Objective

Congenital hyperinsulinism (CHI) is a group of clinically and genetically heterogeneous disorders characterized by dysregulated insulin secretion. The aim of the study was to elucidate genetic etiologies of Taiwanese children with the most severe diazoxide-unresponsive CHI and analyze their genotype-phenotype correlations.

Methods

We combined Sanger with whole exome sequencing (WES) to analyze CHI-related genes. The allele frequency of the most common variant was estimated by single-nucleotide polymorphism haplotype analysis. The functional effects of the ATP-sensitive potassium (KATP) channel variants were assessed using patch clamp recording and Western blot.

Results

Nine of 13 (69%) patients with ten different pathogenic variants (7 in ABCC8, 2 in KCNJ11 and 1 in GCK) were identified by the combined sequencing. The variant ABCC8 p.T1042QfsX75 identified in three probands was located in a specific haplotype. Functional study revealed the human SUR1 (hSUR1)-L366F KATP channels failed to respond to intracellular MgADP and diazoxide while hSUR1-R797Q and hSUR1-R1393C KATP channels were defective in trafficking. One patient had a de novo dominant mutation in the GCK gene (p.I211F), and WES revealed mosaicism of this variant from another patient.

Conclusion

Pathogenic variants in KATP channels are the most common underlying cause of diazoxide-unresponsive CHI in the Taiwanese cohort. The p.T1042QfsX75 variant in the ABCC8 gene is highly suggestive of a founder effect. The I211F mutation in the GCK gene and three rare SUR1 variants associated with defective gating (p.L366F) or traffic (p.R797Q and p.R1393C) KATP channels are also associated with the diazoxide-unresponsive phenotype.

The cysteine-altering p.R544C variant in the NOTCH3 gene is a probable candidate for blood pressure and relevant traits in the Taiwan Biobank

BACKGROUND

The cysteine-altering variants in NOTCH3 have been suggested to be associated with stroke, dementia, and cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL), where aberrant blood pressure levels represent the characteristics of these diseases. We aimed to assess whether the cysteine-altering p.Arg544Cys (p.R544C; rs201118034) variant and common single nucleotide variants (SNVs) in NOTCH3 could contribute to systolic and diastolic blood pressure and related phenotypes in the Taiwan Biobank.

METHODS

We employed a discovery sample of 68,925 individuals, an independent replication sample of 45,676 individuals, and a combined/total sample of 114,601 individuals; all from the Taiwan Biobank. Blood pressure, such as systolic and diastolic blood pressure, was measured for all participants. Association was evaluated using a general linear model, where results were considered statistically significant if the P value < 0.05 divided by the number of independent tests per model.

RESULTS

From our analysis, we identified and replicated three novel candidates for blood pressure that have not previously been reported: the cysteine-altering p.R544C variant for systolic blood pressure, the common SNV rs11669950 for diastolic blood pressure, and the common SNV rs4808235 for diastolic blood pressure. We also generalized two previously identified SNVs (i.e., rs10418305 and rs7408868) in NOTCH3 for blood pressure in European and non-Taiwanese East Asian populations to the Taiwanese population. Moreover, the participants with NOTCH3 p.R544C had an increased stroke frequency (P < 1.0 × 10-5) and a higher dementia frequency (P = 2.0 × 10-4) compared with the whole Taiwan Biobank population in the combined/total sample.

CONCLUSION

NOTCH3 is a strong candidate for a role in stroke, dementia, and CADASIL, which has previously been linked to blood pressure changes. While our preliminary study suggests that NOTCH3 p.R544C may influence blood pressure, stroke, and dementia in the Taiwan Biobank, replication in a well-powered external sample is required. This study also underlines considerable prospects of detecting novel genetic biomarkers in underrepresented worldwide populations.

GJB2 p.V37I Mutation Associated With Moderate Nonsyndromic Hearing Loss in an Adult Taiwanese Population

BACKGROUND

Gap junction protein beta 2 ( GJB2 ) p.V37I mutations are the most important hereditary cause of sensorineural hearing loss (SNHL) in Taiwan. Hearing outcomes are associated with hearing levels at baseline and the duration of follow-up. However, the audiological features of GJB2 p.V37I mutations in the adult population are unknown. The objectives of the present study were to investigate the audiological features, progression rate, and allele frequency of GJB2 p.V37I mutations among an adult Taiwanese population.

METHODS

Subjects of this case-control study were chosen from 13,580 participants of the Taiwan Precision Medicine Initiative. The genetic variations of GJB2 p.V37I were determined by polymerase chain reaction. We analyzed existing pure-tone threshold data from 38 individuals who were homozygous or compound heterozygotes for GJB2 p.V37I, 129 who were heterozygotes, and 602 individuals who were wild-type. Phenome-wide association studies (PheWAS) analysis was also performed to identify phenotypes associated with GJB2 p.V37I.

RESULTS

The minor allele frequency of GJB2 p.V37I was 0.92% in our study population. The mean hearing level of participants with a p.V37I mutation indicated moderate to severe hearing loss with 38.2% ± 22.3% binaural hearing impairment. GJB2 p.V37I was associated with an increased risk of hearing disability (odds ratio: 21.46, 95% confidence interval: 8.62 to 53.44, p < 0.001) in an autosomal recessive pattern. In addition, PheWAS discovered a significant association between GJB2 p.V37I and fracture of the humerus. GJB2 p.V37I is a pathogenic and prevalent variant of SNHL among the adult population.

CONCLUSIONS

The present study recommends patients with known GJB2 p.V37I mutations receive regular audiometric evaluation and genetic counseling. Early assistive listening device intervention is suggested to improve the quality of hearing.

Genetic Factors Contribute to the Phenotypic Variability in GJB2-Related Hearing Impairment

Recessive variants in GJB2 are the most important genetic cause of sensorineural hearing impairment (SNHI) worldwide. Phenotypes vary significantly in GJB2-related SNHI, even in patients with identical variants. For instance, patients homozygous for the GJB2 p.V37I variant, which is highly prevalent in the Asian populations, usually present with mild-to-moderate SNHI; yet severe-to-profound SNHI is occasionally observed in approximately 10% of p.V37I homozygotes. To investigate the genomic underpinnings of the phenotypic variability, we performed next-generation sequencing of GJB2 and other deafness genes in 63 p.V37I homozygotes with extreme phenotypic severities. Additional pathogenic variants of other deafness genes were identified in five of the 35 patients with severe-to-profound SNHI. Furthermore, case-control association analyses were conducted for 30 unrelated p.V37I homozygotes with severe-to-profound SNHI against 28 p.V37I homozygotes with mild-to-moderate SNHI, and 120 population controls from the Taiwan Biobank. The severe-to-profound group exhibited a higher frequency of the crystallin lambda 1 (CRYL1) variant (rs14236), located upstream of GJB2, than the mild-to-moderate and Taiwan Biobank groups. Our results demonstrated that pathogenic variants in other deafness genes and a possible modifier, the CRYL1 rs14236 variant, may contribute to phenotypic variability in GJB2-realted SNHI, highlighting the importance of comprehensive genomic surveys to delineate the genotype-phenotype correlations.

CYP2C19 genotypes and osteoporotic fractures in long-term users of proton pump inhibitors: A hospital-based study

Proton pump inhibitors (PPIs) are commonly prescribed medications. The existing data suggest that individuals at a high risk of fractures have been exposed to high doses of PPIs for prolonged durations. CYP2C19 plays a pivotal role in metabolism of PPIs and thereby influences their pharmacokinetic profile. Hence, we hypothesize that CYP2C19 genotypes may be associated with fragility fracture among PPIs users due to PPI exposure. This study aimed to investigate the association between CYP2C19 genotypes, bone mineral density (BMD), and osteoporotic fracture in a hospital-based population. This retrospective cohort study enrolled patients who were prescribed long-term PPIs at Taichung Veterans General Hospital using data extracted from the Taiwan Precision Medicine Initiative between January 2010 and April 2021. Associations between CYP2C19 phenotypes, comorbidities, and fractures in PPI users were analyzed. We enrolled 1518 long-term PPI users; 571 (38%), 727 (48%), and 220 (14%) CYP2C19 normal metabolizers (NMs), intermediate metabolizers (IMs), and poor metabolizers (PMs), respectively. Among them, 49 (3.2%) patients developed fractures within the 1-year follow-up period; 20 (3.5%) fractures in NMs, 24 (3.3%) in IMs, and 5 (2.3%) in PMs, respectively. No significant difference was observed among CYP2C19 genotypes and fracture. Additionally, BMD measurements during the 1-year follow-up period were made available among 75 participants. No significant difference in BMD between CYP2C19 PMs and non-PMs was found. This real-world, hospital-based study concludes that CYP2C19 PMs/IMs are not associated with an increased risk for fractures or reduced BMD in individuals on long-term PPI therapy.

CYP2C19 Loss-of-Function Variants Associated With Long-Term Ischemic Stroke Events During Clopidogrel Treatment in the Chinese Population

This study aims to determine whether CYP2C19 loss-of-function (LoF) variants were associated with long-term ischemic stroke risk in Chinese primary care patients treated with clopidogrel. Patients treated with clopidogrel were ascertained from Chinese electronic medical records linked with a biobank for a retrospective cohort study. Their medical information was examined for the period from January 2018 to December 2021. Two CYP2C19 major loss of function variants (*2:rs4244285 and *3: rs4986893) were genotyped. The clinical outcome was ischemic stroke event. Cox regression analysis was used to evaluate the association between the occurrence of ischemic stroke events and CYP2C19 LoF variants. Covariates included age, gender, body mass index, prior ischemic stroke, transient ischemic attack, hypertension, diabetes mellitus, hyperlipoidemia, smoke status, aspirin use, proton-pump inhibitor use, and statin use. Of the 1,141 patients included in the clopidogrel therapy cohort, 61.9% carried at least one CYP2C19 LoF variant. During a median follow-up period of 12 months, 103 patients (9.0%) had an ischemic stroke. After adjusting for other risk factors, carriers of CYP2C19 LoF variants had significantly higher risk of ischemic stroke compared with non-carriers (hazard ratio: 1.64, 95% confidence interval: 1.06-2.53, P = 0.025). This pharmacogenetic study of clopidogrel provides novel insights into the association between the CYP2C19 LoF variant and long-term stroke risk. We established that there is still a need for CYP2C19 genotype-guided personalized antiplatelet therapy in those who have returned to the primary care setting for clopidogrel prescription.

Imputation Accuracy Across Global Human Populations

Genotype imputation is now fundamental for genome-wide association studies but lacks fairness due to the underrepresentation of populations with non-European ancestries. The state-of-the-art imputation reference panel released by the Trans-Omics for Precision Medicine (TOPMed) initiative contains a substantial number of admixed African-ancestry and Hispanic/Latino samples to impute these populations with nearly the same accuracy as European-ancestry cohorts. However, imputation for populations primarily residing outside of North America may still fall short in performance due to persisting underrepresentation. To illustrate this point, we curated genome-wide array data from 23 publications published between 2008 to 2021. In total, we imputed over 43k individuals across 123 populations around the world. We identified a number of populations where imputation accuracy paled in comparison to that of European-ancestry populations. For instance, the mean imputation r-squared (Rsq) for 1-5% alleles in Saudi Arabians (N=1061), Vietnamese (N=1264), Thai (N=2435), and Papua New Guineans (N=776) were 0.79, 0.78, 0.76, and 0.62, respectively. In contrast, the mean Rsq ranged from 0.90 to 0.93 for comparable European populations matched in sample size and SNP content. Outside of Africa and Latin America, Rsq appeared to decrease as genetic distances to European reference increased, as predicted. Further analysis using sequencing data as ground truth suggested that imputation software may over-estimate imputation accuracy for non-European populations than European populations, suggesting further disparity between populations. Using 1496 whole genome sequenced individuals from Taiwan Biobank as a reference, we also assessed a strategy to improve imputation for non-European populations with meta-imputation, which can combine results from TOPMed with smaller population-specific reference panels. We found that meta-imputation in this design did not improve Rsq genome-wide. Taken together, our analysis suggests that with the current size of alternative reference panels, meta-imputation alone cannot improve imputation efficacy for underrepresented cohorts and we must ultimately strive to increase diversity and size to promote equity within genetics research.

Identification of Genetic Variants Associated with Severe Myocardial Bridging through Whole-Exome Sequencing

Myocardial bridging (MB) is a congenital coronary artery anomaly and an important cause of angina. The genetic basis of MB is currently unknown. This study used a whole-exome sequencing technique and analyzed genotypic differences. Eight coronary angiography-confirmed cases of severe MB and eight age- and sex-matched control patients were investigated. In total, 139 rare variants that are potentially pathogenic for severe MB were identified in 132 genes. Genes with multiple rare variants or co-predicted by ClinVar and CADD/REVEL for severe MB were collected, from which heart-specific genes were selected under the guidance of tissue expression levels. Functional annotation indicated significant genetic associations with abnormal skeletal muscle mass, cardiomyopathies, and transmembrane ion channels. Candidate genes were reviewed regarding the functions and locations of each individual gene product. Among the gene candidates for severe MB, rare variants in DMD, SGCA, and TTN were determined to be the most crucial. The results suggest that altered anchoring proteins on the cell membrane and intracellular sarcomere unit of cardiomyocytes play a role in the development of the missed trajectory of coronary vessels. Additional studies are required to support the diagnostic application of cardiac sarcoglycan and dystroglycan complexes in patients with severe MB.

Risk Prediction of Chronic Rhinosinusitis with or without Nasal Polyps in Taiwanese Population Using Polygenic Risk Score for Nasal Polyps

The association between single nucleotide polymorphisms and chronic rhinosinusitis (CRS) has been determined. However, it was not known whether the polygenic risk score (PRS) for nasal polyps (NP) could predict CRS with NP (CRSwNP) or without NP (CRSsNP). The aim of this study was to investigate the association between PRSs for NP and the risk of CRS with or without NP. Data from 535 individuals with CRS and 5350 control subjects in the Taiwan Precision Medicine Initiative project were collected. Four PRSs for NP, including PGS000933, PGS000934, PGS001848, and PGS002060 from UK Biobank, were tested in these participants. They were divided into four groups according to quartiles of PRSs. The logistic regression model was performed to evaluate CRSwNP and CRSsNP risk according to PRSs for NP. The PGS002060 had the highest area under the curve at 0.534 for CRSsNP prediction and at 0.588 for CRSwNP prediction. Compared to subjects in the lowest PRS category, the PGS002060 significantly increased the odds for CRSsNP by 1.48 at the highest quintile (p = 0.003) and by 2.32 at the highest quintile for CRSwNP (p = 0.002). In addition, the odds for CRSwNP increased by 3.01 times in female CRSwNP patients (p = 0.009) at the highest quintile compared with those in the lowest PRS category. The PRSs for NP developed from European populations could be applied to the Taiwanese population to predict CRS risk, especially for female CRSwNP.

Curated incidence of lysosomal storage diseases from the Taiwan Biobank

Lysosomal storage diseases (LSDs) are a group of metabolic disorders resulting from a deficiency in one of the lysosomal hydrolases. Most LSDs are inherited in an autosomal or X-linked recessive manner. As LSDs are rare, their true incidence in Taiwan remains unknown. In this study, we used high-coverage whole-genome sequencing data from 1,495 Taiwanese individuals obtained from the Taiwan Biobank. We found 3826 variants in 71 genes responsible for autosomal recessive LSDs. We first excluded benign variants by allele frequency and other criteria. As a result, 270 variants were considered disease-causing. We curated these variants using published guidelines from the American College of Medical Genetics and Genomics (ACMG). Our results revealed a combined incidence rate of 13 per 100,000 (conservative estimation by pathologic and likely pathogenic variants; 95% CI 6.92-22.23) to 94 per 100,000 (extended estimation by the inclusion of variants of unknown significance; 95% CI 75.96-115.03) among 71 autosomal recessive disease-associated genes. The conservative estimations were similar to those in published clinical data. No disease-causing mutations were found for 18 other diseases; thus, these diseases are likely extremely rare in Taiwan. The study results are important for designing screening and treatment methods for LSDs in Taiwan and demonstrate the importance of mutation curation to avoid overestimating disease incidences from genomic data.

Epigenetic regulation of Parkinson's disease risk variant GPNMB cg17274742 methylation by sex and exercise from Taiwan Biobank

Background

Parkinson's disease (PD) is a complex neurodegenerative disease with an elusive etiology that involves the interaction between genetic, behavioral, and environmental factors. Recently, epigenetic modifications, particularly DNA methylation, have been recognized to play an important role in the onset of PD. Glycoprotein non-metastatic melanoma protein B (GPNMB), a type I transmembrane protein crucial for immune cell activation and maturation, has emerged as a potential biomarker for the risk of PD. This research aims to investigate the influence of exercise and gender on the regulation of methylation levels of GPNMB cg17274742 in individuals.

Methods

We analyze data from 2,474 participants in the Taiwan Biobank, collected from 2008 and 2016. Methylation levels at the GPNMB cg17274742 CpG site were measured using Illumina Infinium MethylationEPIC beads. After excluding individuals with incomplete data or missing information on possible risk factors, our final analysis included 1,442 participants. We used multiple linear regression models to assess the association between sex and exercise with adjusted levels of GPNMB cg17274742 for age, BMI, smoking, drinking, coffee consumption, serum uric acid levels, and hypertension.

Results

Our results demonstrated that exercise significantly influenced the methylation levels of GPNMB cg17274742 in males (β = -0.00242; p = 0.0026), but not in females (β = -0.00002362; p = 0.9785). Furthermore, male participants who exercised showed significantly lower levels of methylation compared to the reference groups of the female and non-exercising reference groups (β = -0.00357; p = 0.0079). The effect of the interaction between gender and exercise on the methylation of GPNMB cg17274742 was statistically significant (p = 0.0078).

Conclusion

This study suggests that gender and exercise can modulate GPNMB cg17274742, with hypomethylation observed in exercise men. More research is needed to understand the underlying mechanisms and implications of these epigenetic changes in the context of risk and prevention strategies.

A Rapid Review on the Value of Biobanks Containing Genetic Information

OBJECTIVES

Increasing access to health data through biobanks containing genetic information has the potential to expand the knowledge base and thereby improve screening, diagnosis, and treatment options for many diseases. Nevertheless, although privacy concerns and risks surrounding genetic data sharing are well documented, direct evidence in favor of the hypothesized benefits of data integration is scarce, which complicates decision making in this area. Therefore, the objective of this study is to summarize the available evidence on the research and clinical impacts of biobanks containing genetic information, so as to better understand how to quantify the value of expanding genomic data access.

METHODS

Using a rapid review methodology, we performed a search of MEDLINE/PubMed and Embase databases; and websites of biobanks and genomic initiatives published from 2010 to 2022. We classified findings into 11 indicators including outputs (a direct product of the biobank activities) and outcomes (changes in scientific and clinical capacity).

RESULTS

Of 8479 abstracts and 101 gray literature sources were reviewed, 96 records were included. Although most records did not report key indicators systematically, the available evidence concentrated on research indicators such as publications and gene-disorder association discoveries (63% of studies), followed by research infrastructure (26%), and clinical indicators (11%) such as supporting the diagnosis of individual patients.

CONCLUSIONS

Existing evidence on the benefits of biobanks is skewed toward easily quantifiable research outputs. Measuring a comprehensive set of outputs and outcomes inspired by value frameworks is necessary to generate better evidence on the benefits of genomic data sharing.

Age and sex differences in the association between APOE genotype and Alzheimer's disease in a Taiwan Chinese population

Introduction

The Apolipoprotein E (APOE) epsilon (ε) 4 allele is a well-established risk factor for late-onset Alzheimer's disease (AD). Reports on white ancestry populations have showed that age, sex, and ethnicity have different effects on the association between APOE genotype and AD. However, studies on Asian populations such as Taiwan Chinese populations are limited. This study aimed to evaluate the association between APOE genotype and AD in a Taiwan Chinese population, and to explore if the association varies by age and sex.

Methods

We conducted a case-control study in 725 patients with AD and 1,067 age- and sex- matched controls without dementia from a Taiwan Chinese population. Logistic regression models were used to test the association between AD and APOE genotypes. Secondary analyses considered age (<75 or ≥75 years old), and sex stratified models.

Results

The risk of AD was significantly increased for people with at least one copy of APOE ε4 (OR = 2.52, 95% CI = 2.01-3.17, p < 0.001) and in a dose-dependent manner. Our results did not show an statistically significance different in AD risk when women and men carrying APOEε4 were compared. Despite not reaching statistical significance, the risk of APOE ε4 for AD was higher among younger participants (OR = 3.21, 95% CI = 2.26-4.56, p < 0.001) compared to older ones (OR = 2.13, 95% CI = 1.53-2.97, p < 0.001). When considering both sex and age, the risk of AD was higher among older men carrying APOE ε4 (OR = 2.64, 95% CI = 1.51-4.60 in men; OR = 1.90, 95% CI = 1.26-2.86 in women), while women carrying APOE ε4 appeared to have an increased risk at a younger age (OR = 3.29, 95% CI = 2.20-4.93 in women; OR = 2.91, 95% CI = 1.40-6.05 in men).

Discussion

The APOE ε4 allele represents a major risk factor for AD in the Taiwanese population. The effect of APOE ε4 allele on AD risk appeared to be stronger among men aged 75 years or more and among younger women.

Hypertension as a Novel Link for Shared Heritability in Age at Menarche and Cardiometabolic Traits

CONTEXT

Extremely early age at menarche, also called precocious puberty, has been associated with various cardiometabolic traits, but their shared heritability remains unclear.

OBJECTIVES

This work aimed to identify new shared genetic variants and their pathways for age at menarche and cardiometabolic traits and to investigate the influence of central precocious puberty on childhood cardiometabolic traits.

METHODS

Using the conjunction false discovery rate method, this study analyzed genome-wide association study data from the menarche-cardiometabolic traits among 59 655 females of Taiwanese ancestry and systemically investigated pleiotropy between age at menarche and cardiometabolic traits. To support the novel hypertension link, we used the Taiwan Puberty Longitudinal Study (TPLS) to investigate the influence of precocious puberty on childhood cardiometabolic traits.

RESULTS

We discovered 27 novel loci, with an overlap between age at menarche and cardiometabolic traits, including body fat and blood pressure. Among the novel genes discovered, SEC16B, CSK, CYP1A1, FTO, and USB1 are within a protein interaction network with known cardiometabolic genes, including traits for obesity and hypertension. These loci were confirmed through demonstration of significant changes in the methylation or expression levels of neighboring genes. Moreover, the TPLS provided evidence regarding a 2-fold higher risk of early-onset hypertension that occurred in girls with central precocious puberty.

CONCLUSION

Our study highlights the usefulness of cross-trait analyses for identifying shared etiology between age at menarche and cardiometabolic traits, especially early-onset hypertension. The menarche-related loci may contribute to early-onset hypertension through endocrinological pathways.

Efficacy of Warfarin Therapy Guided by Pharmacogenetics: A Real-world Investigation Among Han Taiwanese

PURPOSE

The anticoagulation activity of warfarin in populations with CYP2C9, VKORC1, and CYP4F2 variants differs between individuals and is correlated with poor international normalized ratio (INR) control. Pharmacogenetics-guided warfarin dosing has been successfully developed for patients with genetic variations in recent years. However, few real-world data have been used to investigate the INR and warfarin dosage and the time to target INR. This study examined the largest collection of genetic and clinical real-world data related to warfarin to provide further evidence supporting the benefits of pharmacogenetics in clinical outcomes.

METHODS

We retrieved a total of 69,610 INR-warfarin records after the index date from 2,613 patients in the China Medical University Hospital database between January 2003 and December 2019. Each INR reading was obtained from the latest laboratory data after the hospital visit date. Patients with a history of malignant neoplasms or pregnancy before the index date were excluded, as were patients without data on INR measurements after the fifth day of prescription, genetic information, or gender variables. The primary outcomes were the INR and warfarin dosage during days 7, 14, 28, 56, and 84 after prescription. The secondary outcome was the time required to reach the INR ranges of 1.5 to 3.0 and >4.0.

FINDINGS

A total of 59,643 INR-warfarin records from 2188 patients were retrieved. The average INR was higher for homozygous carriers of the minor allele at CYP2C9 and VKORC1 during the first 7 days (1.83 [1.03] [CYP2C9*1] and 2.46 [1.44] [CYP2C9*3], P < 0.001; 1.39 [0.36] [rs9923231 G/G], 1.55 [0.79] [rs9923231 G/A], and 1.96 [1.13] [rs9923231 A/A], P < 0.001) than for the wild-type allele. These patients with variants required lower warfarin doses than those with the wild-type allele during the first 28 days. CYP4F2 variant patients seemed to require higher doses of warfarin than those in the wild-type group; however, no significant difference in the average INR was observed (1.95 [1.14] [homozygous V433 carriers], 1.78 [0.98] [heterozygous V433M carriers], and 1.66 [0.91] [homozygous M433 carriers], P = 0.016).

IMPLICATIONS

Our study indicates that genetic variants in the Han population may enhance warfarin responsiveness, which holds clinical relevance. An increased warfarin dosage was not linked to a shorter time to therapeutic INR between CYP4F2 variant patients and those with a wild-type allele. Assessing CYP2C9 and VKORC1 genetic polymorphisms before initiating warfarin treatment in real-world practice is essential for potentially vulnerable patients and is likely to optimize therapeutic dosing.

South Asian medical cohorts reveal strong founder effects and high rates of homozygosity

The benefits of large-scale genetic studies for healthcare of the populations studied are well documented, but these genetic studies have traditionally ignored people from some parts of the world, such as South Asia. Here we describe whole genome sequence (WGS) data from 4806 individuals recruited from the healthcare delivery systems of Pakistan, India and Bangladesh, combined with WGS from 927 individuals from isolated South Asian populations. We characterize population structure in South Asia and describe a genotyping array (SARGAM) and imputation reference panel that are optimized for South Asian genomes. We find evidence for high rates of reproductive isolation, endogamy and consanguinity that vary across the subcontinent and that lead to levels of rare homozygotes that reach 100 times that seen in outbred populations. Founder effects increase the power to associate functional variants with disease processes and make South Asia a uniquely powerful place for population-scale genetic studies.