Comparative analysis of genetic risk scores for predicting biochemical recurrence in prostate cancer patients after radical prostatectomy.

BACKGROUND: In recent years, Genome-Wide Association Studies (GWAS) has identified risk variants related to complex diseases, but most genetic variants have less impact on phenotypes. To solve the above problems, methods that can use variants with low genetic effects, such as genetic risk score (GRS), have been developed to predict disease risk. METHODS: As the GRS model with the most incredible prediction power for complex diseases has not been determined, our study used simulation data and prostate cancer data to explore the disease prediction power of three GRS models, including the simple count genetic risk score (SC-GRS), the direct logistic regression genetic risk score (DL-GRS), and the explained variance weighted GRS based on directed logistic regression (EVDL-GRS). RESULTS AND CONCLUSIONS: We used 26 SNPs to establish GRS models to predict the risk of biochemical recurrence (BCR) after radical prostatectomy. Combining clinical variables such as age at diagnosis, body mass index, prostate-specific antigen, Gleason score, pathologic T stage, and surgical margin and GRS models has better predictive power for BCR. The results of simulation data (statistical power = 0.707) and prostate cancer data (area under curve = 0.8462) show that DL-GRS has the best prediction performance. The rs455192 was the most relevant locus for BCR (p = 2.496 × 10-6) in our study.

Genetic risk score constructed from common genetic variants is associated with cardiovascular disease risk in type 2 diabetes mellitus.

BACKGROUND: Patients with type 2 diabetes mellitus (T2DM) experience a two-fold increased risk of cardiovascular diseases. Genome-wide association studies (GWAS) have identified T2DM susceptibility genetic variants. Interestingly, the genetic variants associated with cardiovascular disease risk in T2DM Han Chinese remain to be elucidated. The present study aimed to investigate the genetic variants associated with cardiovascular disease risk in T2DM. METHODS: We performed bootstrapping, GWAS and an investigation of genetic variants associated with cardiovascular disease risk in a discovery T2DM cohort and in a replication cohort. The discovery cohort included 326 cardiovascular disease patients and 1209 noncardiovascular disease patients. The replication cohort included 68 cardiovascular disease patients and 317 noncardiovascular disease patients. The main outcome measures were genetic variants for genetic risk score (GRS) in cardiovascular disease risk in T2DM. RESULTS: In total, 35 genetic variants were associated with cardiovascular disease risk. A GRS was generated by combining risk alleles from these variants weighted by their estimated effect sizes (log odds ratio [OR]). T2DM patients with weighted GRS ≥ 12.63 had an approximately 15-fold increase in cardiovascular disease risk (odds ratio = 15.67, 95% confidence interval [CI] = 10.33-24.00) compared to patients with weighted GRS < 10.39. With the addition of weighted GRS, receiver-operating characteristic curves showed that area under the curve with conventional risk factors was improved from 0.719 (95% CI = 0.689-0.750) to 0.888 (95% CI = 0.866-0.910). CONCLUSIONS: These 35 genetic variants are associated with cardiovascular disease risk in T2DM, alone and cumulatively. T2DM patients with higher levels of weighted genetic risk score have higher cardiovascular disease risks.

Association of Plasma Thrombospondin-1 Level with Cardiovascular Disease and Mortality in Hemodialysis Patients.

BACKGROUND: Thrombospondin-1 (TSP-1) is known to be involved in the regulation of angiogenesis, inflammation, and vascular function. Clinical studies have demonstrated its correlation with peripheral artery disease, coronary artery disease, and pulmonary hypertension. In this study, we explored its potential roles in the background of end-stage renal disease (ESRD). METHODS: A total of 140 ESRD outpatients (ages 61.0 ± 12.4 years) were prospectively followed for 34 ± 7 months. Their TSP-1 levels were analyzed from pre-hemodialysis blood sample. Cardiovascular survey included ankle- brachial index (ABI), echocardiography and Tl-201 dipyridamole single-photon emission computed tomography (SPECT). RESULTS: Plasma TSP-1 levels were higher in those patients with preexisting clinical evidence of cardiovascular disease (CVD) than those without (p = 0.002). TSP-1 concentrations were also correlated with ABI, left ventricular ejection fraction, and scar burden in SPECT. Stepwise logistic regression analysis revealed that TSP-1 level was independently associated with the presence of CVD, with an odds ratio of 1.38 [95% confidence interval (CI), 1.09-1.75, p = 0.008]. In survival analyses, 31 patients (22%) died during the follow-up, 16 (52%) arising from cardiovascular causes. Cox hazards analysis revealed that the patients with TSP-1 levels in the highest tertile had a 5.32- and 6.75-fold higher risk for all-cause and cardiovascular mortality than those in the lowest tertile. This predictive value for all-cause mortality still persisted after multivariate adjustment (hazard ratio, 8.71; 95% CI, 1.36-55.68; p = 0.02). CONCLUSIONS: This study hallmarks the association of elevated TSP-1 level with CVD and adverse outcome among hemodialysis patients. KEY WORDS: Thrombospondin-1; End-stage renal disease; Cardiovascular disease; Mortality.

Predicting HLA genotypes using unphased and flanking single-nucleotide polymorphisms in Han Chinese population.

BACKGROUND: Genetic variation associated with human leukocyte antigen (HLA) genes has immunological functions and is associated with autoimmune diseases. To date, large-scale studies involving classical HLA genes have been limited by time-consuming and expensive HLA-typing technologies. To reduce these costs, single-nucleotide polymorphisms (SNPs) have been used to predict HLA-allele types. Although HLA allelic distributions differ among populations, most prediction model of HLA genes are based on Caucasian samples, with few reported studies involving non-Caucasians. RESULTS: Our sample consisted of 437 Han Chinese with Affymetrix 5.0 and Illumina 550 K SNPs, of whom 214 also had data on Affymetrix 6.0 SNPs. All individuals had HLA typings at a 4-digit resolution. Using these data, we have built prediction model of HLA genes that are specific for a Han Chinese population. To optimize our prediction model of HLA genes, we analyzed a number of critical parameters, including flanking-region size, genotyping platform, and imputation. Predictive accuracies generally increased both with sample size and SNP density. CONCLUSIONS: SNP data from the HapMap Project are about five times more dense than commercially available genotype chip data. Using chips to genotype our samples, however, only reduced the accuracy of our HLA predictions by only ~3%, while saving a great deal of time and expense. We demonstrated that classical HLA alleles can be predicted from SNP genotype data with a high level of accuracy (80.37% (HLA-B) ~95.79% (HLA-DQB1)) in a Han Chinese population. This finding offers new opportunities for researchers in obtaining HLA genotypes via prediction using their already existing chip datasets. Since the genetic variation structure (e.g. SNP, HLA, Linkage disequilibrium) is different between Han Chinese and Caucasians, and has strong impact in building prediction models for HLA genes, our findings emphasize the importance of building ethnic-specific models when analyzing human populations.

Identifying rare and common disease associated variants in genomic data using Parkinson's disease as a model.

BACKGROUND: Genome-wide association studies have been successful in identifying common genetic variants for human diseases. However, much of the heritable variation associated with diseases such as Parkinson's disease remains unknown suggesting that many more risk loci are yet to be identified. Rare variants have become important in disease association studies for explaining missing heritability. Methods for detecting this type of association require prior knowledge on candidate genes and combining variants within the region. These methods may suffer from power loss in situations with many neutral variants or causal variants with opposite effects. RESULTS: We propose a method capable of scanning genetic variants to identify the region most likely harbouring disease gene with rare and/or common causal variants. Our method assigns a score at each individual variant based on our scoring system. It uses aggregate scores to identify the region with disease association. We evaluate performance by simulation based on 1000 Genomes sequencing data and compare with three commonly used methods. We use a Parkinson's disease case-control dataset as a model to demonstrate the application of our method. Our method has better power than CMC and WSS and similar power to SKAT-O with well-controlled type I error under simulation based on 1000 Genomes sequencing data. In real data analysis, we confirm the association of α-synuclein gene (SNCA) with Parkinson's disease (p = 0.005). We further identify association with hyaluronan synthase 2 (HAS2, p = 0.028) and kringle containing transmembrane protein 1 (KREMEN1, p = 0.006). KREMEN1 is associated with Wnt signalling pathway which has been shown to play an important role for neurodegeneration in Parkinson's disease. CONCLUSIONS: Our method is time efficient and less sensitive to inclusion of neutral variants and direction effect of causal variants. It can narrow down a genomic region or a chromosome to a disease associated region. Using Parkinson's disease as a model, our method not only confirms association for a known gene but also identifies two genes previously found by other studies. In spite of many existing methods, we conclude that our method serves as an efficient alternative for exploring genomic data containing both rare and common variants.

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.

Body constitutions of traditional Chinese medicine caused a significant effect on irritable bowel syndrome.

BACKGROUND: According to the theory of traditional Chinese medicine (TCM), all types of body constitutions, except for the Gentleness (ie, the control group in our study), have disease susceptibility and affect the disease development process. This study attempted to investigate the relationship between TCM body constitutions and irritable bowel syndrome (IBS). METHODS: This cross-sectional study was based on Taiwan Biobank (TWB) and collected clinical data from 13 941 subjects aged 30 to 70. The results of the study showed that subjects with Yang-deficiency (N = 3161 subjects, odds ratio [OR] = 2.654, 95% CI = 1.740-3.910), Ying-deficiency (N = 3331 subjects, OR = 1.096, 95% CI = 0.627-1.782) or Stasis (N = 2335 subjects, OR = 1.680, 95% CI = 0.654-3.520) were more likely to have IBS. RESULTS: If the subjects with two or more TCM body constitutions: Yang-deficiency + Ying-deficiency (OR = 3.948, 95% CI = 2.742-5.560), Yang-deficiency + Stasis (OR = 2.312, 95% CI = 1.170-4.112), Ying-deficiency + Stasis (OR = 1.851, 95% CI = 0.828-3.567), or Yang-deficiency + Ying-deficiency + Stasis (OR = 3.826, 95% CI = 2.954-4.932) were also prone to IBS. CONCLUSION: These results confirmed the high correlation between TCM body constitutions and IBS. Because the current treatment for IBS is not entirely satisfactory, integrated traditional Chinese and Western medicine might provide patients with an alternative treatment option to alleviate IBS.

Incidence and Nature of Short-Term Adverse Events following COVID-19 Second Boosters: Insights from Taiwan's Universal Vaccination Strategy.

This study evaluates the incidence and characteristics of adverse events (AEs) following the second COVID-19 booster dose, leveraging Taiwan's distinctive approach of extending booster vaccinations to all citizens, unlike the targeted high-risk group strategies in other countries. Utilizing data from Taipei Veterans General Hospital's Vaccine Adverse Event Reporting System (VAERS) from 27 October 2022 to 19 January 2023, this research examines AEs in 441 out of 1711 booster recipients, considering factors like age, vaccine brands, and booster combinations. The findings revealed incidence rates (IRs) of 25.6% (95% CI: 21.1-30.8) after the first booster and 24.9% (95% CI: 20.5-30.0) after the second, mostly non-serious, with those having AEs post-first booster being five times more likely to report them again (incidence rate ratio, 5.02, p < 0.001). Significantly, switching from the mRNA1273 vaccine to another brand reduced AE risk by 18%. This study underscores that AEs are more repetitive than cumulative with additional booster doses, advocating for personalized vaccination strategies based on individual medical histories and previous vaccine reactions. These insights are valuable for healthcare providers in discussing potential AEs with patients, thereby improving vaccine compliance and public trust, and for policymakers in planning future booster vaccination strategies.

Decoding and reconstructing disease relations between dry eye and depression: a multimodal investigation comprising meta-analysis, genetic pathways and Mendelian randomization.

INTRODUCTION: The clinical presentations of dry eye disease (DED) and depression (DEP) often comanifest. However, the robustness and the mechanisms underlying this association were undetermined. OBJECTIVES: To this end, we set up a three-segment study that employed multimodality results (meta-analysis, genome-wide association study [GWAS] and Mendelian randomization [MR]) to elucidate the association, common pathways and causality between DED and DEP. METHODS: A meta-analysis comprising 26 case-control studies was first conducted to confirm the DED-DEP association. Next, we performed a linkage disequilibrium (LD)-adjusted GWAS and targeted phenotype association study (PheWAS) in East Asian TW Biobank (TWB) and European UK Biobank (UKB) populations. Single-nucleotide polymorphisms (SNPs) were further screened for molecular interactions and common pathways at the functional gene level. To further elucidate the activated pathways in DED and DEP, a systemic transcriptome review was conducted on RNA sequencing samples from the Gene Expression Omnibus. Finally, 48 MR experiments were implemented to examine the bidirectional causation between DED and DEP. RESULTS: Our meta-analysis showed that DED patients are associated with an increased DEP prevalence (OR = 1.83), while DEP patients have a concurrent higher risk of DED (OR = 2.34). Notably, cross-disease GWAS analysis revealed that similar genetic architecture (rG = 0.19) and pleiotropic functional genes contributed to phenotypes in both diseases. Through protein-protein interaction and ontology convergence, we summarized the pleiotropic functional genes under the ontology of immune activation, which was further validated by a transcriptome systemic review. Importantly, the inverse variance-weighted (IVW)-MR experiments in both TWB and UKB populations (p value <0.001) supported the bidirectional exposure-outcome causation for DED-to-DEP and DEP-to-DED. Despite stringent LD-corrected instrumental variable re-selection, the bidirectional causation between DED and DEP remained. CONCLUSION: With the multi-modal evidence combined, we consolidated the association and causation between DED and DEP.

Using maximal segmental score in genome-wide association studies.

Genome-wide association studies (GWAS) have become the method of choice for identifying disease susceptibility genes in common disease genetics research. Despite successes in these studies, much of the heritability remains unexplained due to lack of power and low resolution. High-density genotyping arrays can now screen more than 5 million genetic markers. As a result, multiple comparison has become an important issue especially in the era of next-generation sequencing. We propose to use a two-stage maximal segmental score procedure (MSS) which uses region-specific empirical P-values to identify genomic segments most likely harboring the disease gene. We develop scoring systems based on Fisher's P-value combining method to convert locus-specific significance levels into region-specific scores. Through simulations, our result indicated that MSS increased the power to detect genetic association as compared with conventional methods provided type I error was at 5%. We demonstrated the application of MSS on a publicly available case-control dataset of Parkinson's disease and replicated the findings in the literature. MSS provides an efficient exploratory tool for high-density association data in the current era of next-generation sequencing. R source codes to implement the MSS procedure are freely available at http://www.csjfann.ibms.sinica.edu.tw/EAG/program/programlist.htm.

Serum adipocyte fatty acid-binding protein levels in patients with critical illness are associated with insulin resistance and predict mortality.

INTRODUCTION: Hyperglycemia and insulin resistance are commonplace in critical illness, especially in patients with sepsis. Recently, several hormones secreted by adipose tissue have been determined to be involved in overall insulin sensitivity in metabolic syndrome-related conditions, including adipocyte fatty-acid binding protein (A-FABP). However, little is known about their roles in critical illness. On the other hand, there is evidence that several adipose tissue gene expressions change in critically ill patients. METHODS: A total of 120 patients (72 with sepsis, 48 without sepsis) were studied prospectively on admission to a medical ICU and compared with 45 healthy volunteers as controls. Various laboratory parameters and metabolic and inflammatory profiles were assessed within 48 hours after admission. Clinical data were collected from medical records. RESULTS: Compared with healthy controls, serum A-FABP concentrations were higher in all critically ill patients, and there was a trend of higher A-FABP in patients with sepsis. In multivariate correlation analysis in all critically ill patients, the serum A-FABP concentrations were independently related to serum creatinine, fasting plasma glucose, total cholesterol, TNF-alpha, albumin, and the Acute Physiology and Chronic Health Evaluation II scores. In survival analysis, higher A-FABP levels (> 40 ng/ml) were associated with an unfavorable overall survival outcome, especially in sepsis patients. CONCLUSIONS: Critically ill patients have higher serum A-FABP concentrations. Moreover, A-FABP may potentially serve as a prognostic biomarker in critically ill patients with sepsis.

TMEM132D and VIPR2 Polymorphisms as Genetic Risk Loci for Retinal Detachment: A Genome-Wide Association Study and Polygenic Risk Score Analysis.

Purpose: Retinal detachment (RD) is a sight-threatening ocular disease caused by separation of the neurosensory retina from the underlying retinal pigment epithelium layer. Its genetic basis is unclear because of a limited amount of data. In this study, we aimed to identify genetic risk loci associated with RD in participants without diabetes mellitus and to construct a polygenic risk score (PRS) to predict the risk of RD. Methods: A genome-wide association study was conducted using data from the Taiwan Biobank to identify RD risk loci. A total of 1533 RD cases and 106,270 controls were recruited, all of whom were Han Chinese. Replication studies were performed using data from the UK Biobank and Biobank Japan. To construct the PRS, a traditional clumping and thresholding method was performed and validated by fivefold cross-validation. Results: Two novel loci with significant associations were identified. These two genes were TMEM132D (lead single nucleotide polymorphism [SNP]: rs264498, adjusted-P = 7.18 × 10-9) and VIPR2 (lead SNP: rs3812305, adjusted-P = 8.38 × 10-9). The developed PRS was effective in discriminating individuals at high risk of RD with a dose-response relationship. The quartile with the highest risk had an odds ratio of 1244.748 compared to the lowest risk group (95% confidence interval, 175.174-8844.892). Conclusions: TMEM132D and VIPR2 polymorphisms are genetic candidates linked to RD in Han Chinese populations. Our proposed PRS was effective at discriminating high-risk from low-risk individuals.

On the use of multifactor dimensionality reduction (MDR) and classification and regression tree (CART) to identify haplotype-haplotype interactions in genetic studies.

Haplotype-based approaches may have greater power than single-locus analyses when the SNPs are in strong linkage disequilibrium with the risk locus. To overcome potential complexities owing to large numbers of haplotypes in genetic studies, we evaluated two data mining approaches, multifactor dimensionality reduction (MDR) and classification and regression tree (CART), with the concept of haplotypes considering their haplotype uncertainty to detect haplotype-haplotype (HH) interactions. In evaluation of performance for detecting HH interactions, MDR had higher power than CART, but MDR gave a slightly higher type I error. Additionally, we performed an HH interaction analysis with a publicly available dataset of Parkinson's disease and confirmed previous findings that the RET proto-oncogene is associated with the disease. In this study, we showed that using HH interaction analysis is possible to assist researchers in gaining more insight into identifying genetic risk factors for complex diseases.

Using Image Recognition to Process Unbalanced Data in Genetic Diseases From Biobanks.

With precision medicine as the goal, the human biobank of each country should be analyzed to determine the complete research results related to genetic diseases. In addition, with the increase in medical imaging data, automatic image processing with image recognition has been widely studied and applied in biomedicine. However, case-control data imbalance often occurs in human biobanks, which is usually solved by the statistical method SAIGE. Due to the huge amount of genetic data in human biobanks, the direct use of the SAIGE method often faces the problem of insufficient computer memory to support calculations and excessive calculation time. The other method is to use sampling to adjust the data to balance the case-control ratio, which is called Synthetic Minority Oversampling Technique (SMOTE). Our study employed the Manhattan plot and genetic disease information from the Taiwan Biobank to adjust the imbalance in the case-control ratio by SMOTE, called "TW-SMOTE." We further used a deep learning image recognition system to identify the TW-SMOTE. We found that TW-SMOTE can achieve the same results as that of SAIGE and the UK Biobank (UKB). The processing of the technical data can be equivalent to the use of data plots with a relatively large UKB sample size and achieve the same effect as that of SAIGE in addressing data imbalance.

Polygenic Risk Score Improves Cataract Prediction in East Asian Population.

Cataracts, characterized by crystalline lens opacities in human eyes, is the leading cause of blindness globally. Due to its multifactorial complexity, the molecular mechanisms remain poorly understood. Larger cohorts of genome-wide association studies (GWAS) are needed to investigate cataracts' genetic basis. In this study, a GWAS was performed on the largest Han population to date, analyzing a total of 7079 patients and 13,256 controls from the Taiwan Biobank (TWB) 2.0 cohort. Two cataract-associated SNPs with an adjustment of p < 1 × 10−7 in the older groups and nine SNPs with an adjustment of p < 1 × 10−6 in the younger group were identified. Except for the reported AGMO in animal models, most variations, including rs74774546 in GJA1 and rs237885 in OXTR, were not identified before this study. Furthermore, a polygenic risk score (PRS) was created for the young and old populations to identify high-risk cataract individuals, with areas under the receiver operating curve (AUROCs) of 0.829 and 0.785, respectively, after covariate adjustments. Younger individuals had 17.45 times the risk while older people had 10.97 times the risk when comparing individuals in the highest and lowest PRS quantiles. Validation analysis on an independent TWB1.0 cohort revealed AUROCs of 0.744 and 0.659.

Phenome-wide analysis of Taiwan Biobank reveals novel glycemia-related loci and genetic risks for diabetes.

To explore the complex genetic architecture of common diseases and traits, we conducted comprehensive PheWAS of ten diseases and 34 quantitative traits in the community-based Taiwan Biobank (TWB). We identified 995 significantly associated loci with 135 novel loci specific to Taiwanese population. Further analyses highlighted the genetic pleiotropy of loci related to complex disease and associated quantitative traits. Extensive analysis on glycaemic phenotypes (T2D, fasting glucose and HbA1c) was performed and identified 115 significant loci with four novel genetic variants (HACL1, RAD21, ASH1L and GAK). Transcriptomics data also strengthen the relevancy of the findings to metabolic disorders, thus contributing to better understanding of pathogenesis. In addition, genetic risk scores are constructed and validated for absolute risks prediction of T2D in Taiwanese population. In conclusion, our data-driven approach without a priori hypothesis is useful for novel gene discovery and validation on top of disease risk prediction for unique non-European population.

Hepatitis B virus persistent infection-related single nucleotide polymorphisms in HLA regions are associated with viral load in hepatoma families.

BACKGROUND: Genome-wide association studies from Asia indicate that HLA-DP and HLA-DQ loci are important in persistent hepatitis B virus (HBV) infections. One of the key elements for HBV-related carcinogenesis is persistent viral replication and inflammation. AIM: To examine genetic and nongenetic factors with persistent HBV infection and viral load in families with hepatocellular carcinoma (HCC). METHODS: The HCC families included 301 hepatitis B surface antigen (HBsAg) carriers and 424 noncarriers born before the nationwide vaccination program was initiated in 1984. Five HBV-related single nucleotide polymorphisms (SNPs) - rs477515, rs9272105, rs9276370, rs7756516, and rs9277535 - were genotyped. Factors associated with persistent HBV infection and viral load were analyzed by a generalized estimating equation. RESULTS: In the first-stage persistent HBV study, all SNPs except rs9272105 were associated with persistent infection. A significantly higher area under the reciprocal operating characteristic curve for nongenetic factors vs genetic factors (P < 0.001) suggests that the former play a major role in persistent HBV infection. In the second-stage viral load study, we added 8 HBsAg carriers born after 1984. The 309 HBsAg carriers were divided into low (n = 162) and high viral load (n = 147) groups with an HBV DNA cutoff of 105 cps/mL. Sex, relationship to the index case, rs477515, rs9272105, and rs7756516 were associated with viral load. Based on the receiver operating characteristic curve analysis, genetic and nongenetic factors affected viral load equally in the HCC family cohort (P = 0.3117). CONCLUSION: In these east Asian adults, the mechanism of persistent HBV infection-related SNPs was a prolonged viral replication phase.

Genetic factors of idiopathic central precocious puberty and their polygenic risk in early puberty.

OBJECTIVE,: To investigate the genetic characteristics of idiopathic central precocious puberty (ICPP) and validate its polygenic risk for early puberty. DESIGN AND METHODS: A bootstrap subsampling and genome-wide association study were performed on Taiwanese Han Chinese girls comprising 321 ICPP patients and 148 controls. Using previous GWAS data on pubertal timing, a replication study was performed. A validation group was also investigated for the weighted polygenic risk score (wPRS) of the risk of early puberty. RESULTS: A total of 105 SNPs for the risk of ICPP were identified, of which 22 yielded an area under the receiver operating characteristic curve of 0.713 for the risk of early puberty in the validation group. A replication study showed that 33 SNPs from previous GWAS data of pubertal timing were associated with the risk of ICPP (training group: P-value < 0.05). In the validation group, a cumulative effect was observed between the wPRS and the risk of early puberty in a dose-dependent manner (validation group: Cochran-Armitage trend test: P-value < 1.00E-04; wPRS quartile 2 (Q2) (odds ratio (OR) = 5.00, 95% CI: 1.55-16.16), and wPRS Q3 (OR = 11.67, 95% CI: 2.44-55.83)). CONCLUSIONS: This study reveals the ICPP genetic characteristics with 22 independent and 33 reported SNPs in the Han Chinese population from Taiwan. This study may contribute to understand the genetic features and underlying biological pathways that control pubertal timing and pathogenesis of ICPP and also to the identification of individuals with a potential genetic risk of early puberty.

Genome-Wide Polygenic Risk Score for Predicting High Risk Glaucoma Individuals of Han Chinese Ancestry.

Glaucoma is a progressive and irreversible blindness-causing disease. However, the underlying genetic factors and molecular mechanisms remain poorly understood. Previous genome-wide association studies (GWAS) have made tremendous progress on the SNP-based disease association and characterization. However, most of them were conducted for Europeans. Since differential genetic characteristics among ethnic groups were evident in glaucoma, it is worthwhile to complete its genetic landscape from the larger cohorts of Asian individuals. Here, we present a GWAS based on the Taiwan Biobank. Among 1013 glaucoma patients and 36,562 controls, we identified a total of 138 independent glaucoma-associated SNPs at the significance level of p < 1 × 10-5. After clumping genetically linked SNPs (LD clumping), 134 independent SNPs with p < 10-4 were recruited to construct a Polygenic Risk Score (PRS). The model achieved an area under the receiver operating characteristic curve (AUC) of 0.8387 (95% CI = [0.8269-0.8506]), and those within the top PRS quantile had a 45.48-fold increased risk of glaucoma compared with those within the lowest quantile. The PRS model was validated with an independent cohort that achieved an AUC of 0.7283, thereby showing the effectiveness of our polygenic risk score in predicting individuals in the Han Chinese population with higher glaucoma risks.