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

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

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

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

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

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

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

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

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

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

Genetic insights into ossification of the posterior longitudinal ligament of the spine.

Ossification of the posterior longitudinal ligament of the spine (OPLL) is an intractable disease leading to severe neurological deficits. Its etiology and pathogenesis are primarily unknown. The relationship between OPLL and comorbidities, especially type 2 diabetes (T2D) and high body mass index (BMI), has been the focus of attention; however, no trait has been proven to have a causal relationship. We conducted a meta-analysis of genome-wide association studies (GWASs) using 22,016 Japanese individuals and identified 14 significant loci, 8 of which were previously unreported. We then conducted a gene-based association analysis and a transcriptome-wide Mendelian randomization approach and identified three candidate genes for each. Partitioning heritability enrichment analyses observed significant enrichment of the polygenic signals in the active enhancers of the connective/bone cell group, especially H3K27ac in chondrogenic differentiation cells, as well as the immune/hematopoietic cell group. Single-cell RNA sequencing of Achilles tendon cells from a mouse Achilles tendon ossification model confirmed the expression of genes in GWAS and post-GWAS analyses in mesenchymal and immune cells. Genetic correlations with 96 complex traits showed positive correlations with T2D and BMI and a negative correlation with cerebral aneurysm. Mendelian randomization analysis demonstrated a significant causal effect of increased BMI and high bone mineral density on OPLL. We evaluated the clinical images in detail and classified OPLL into cervical, thoracic, and the other types. GWAS subanalyses identified subtype-specific signals. A polygenic risk score for BMI demonstrated that the effect of BMI was particularly strong in thoracic OPLL. Our study provides genetic insight into the etiology and pathogenesis of OPLL and is expected to serve as a basis for future treatment development.

Combining polygenic risk scores and human leukocyte antigen variants for personalized risk assessment of type 1 diabetes in the Taiwanese population.

AIMS: To analyse the genome-wide association study (GWAS) data of patients with type 1 diabetes mellitus (T1D) in order to develop a risk score for the genetic effects on T1D risk and age at diagnosis in the Taiwanese population. MATERIALS AND METHODS: We selected 610 patients with T1D and 2511 healthy individuals from an electronic medical record database of more than 300 000 individuals with genetic information, analysed their GWAS data, and developed a polygenic risk score (PRS). RESULTS: The PRS, based on 149 selected single-nucleotide polymorphisms, could effectively predict T1D risk. A PRS increase was associated with increased T1D risk (odds ratio [OR] 2.09, 95% confidence interval [CI] 1.72-2.55). Moreover, a 1-unit increase in standardized T1D PRS decreased the age at diagnosis by 0.74 years. Combined PRS and human leukocyte antigen (HLA) DQA1*03:02-DQA1*05:01 genotypes could accurately predict T1D risk. In multivariable models, HLA variants and PRS were independent risk factors for T1D risk (OR 3.76 [95% CI 1.54-9.16] and 1.71 [95% CI 1.37-2.13] for HLA DQA1*03:02-DQA1*05:01 and PRS, respectively). In a limited study population of those aged ≤18 years, PRS remained significantly associated with T1D risk. The association between T1D PRS and age at diagnosis was more obvious among males and patients aged ≤18 years. CONCLUSIONS: Polygenic risk score and HLA variations enable personalized risk estimates, enhance newborn screening efficiency for ketoacidosis prevention, and addresses the gap in data on T1D prediction in isolated Asian populations.

Association of combination antiretroviral therapy with risk of neurological diseases in patients with HIV/AIDS in Taiwan: a nested case-control study.

Heterogeneous neurocognitive impairment remains an important issue, even in the era of combination antiretroviral therapy (cART), with an incidence ranging from 15% to 65%. Although ART drugs with higher penetration scores to the central nervous system (CNS) show better HIV replication control in the CNS, the association between CNS penetration effectiveness (CPE) scores and neurocognitive impairment remains inconclusive. To explore whether ART exposure is associated with the risk of neurological diseases among patients with HIV/AIDS, this study in Taiwan involved 2,571 patients with neurological diseases and 10,284 matched, randomly selected patients without neurological diseases between 2010 and 2017. A conditional logistic regression model was used in this study. The parameters for ART exposure included ART usage, timing of exposure, cumulative defined daily dose (DDD), adherence, and cumulative CPE score. Incident cases of neurological diseases, including CNS infections, cognitive disorders, vasculopathy, and peripheral neuropathy, were obtained from the National Health Insurance Research Database in Taiwan. Odds ratios (ORs) for the risk of neurological diseases were conducted using a multivariate conditional logistic regression model. Patients with a history of past exposure (OR: 1.68, 95% confidence interval [CI]:1.22-2.32), low cumulative DDDs (< 2,500) (OR: 1.28, 95% CI: 1.15-1.42), low adherence (0 < adherence (ADH) ≤ 0.8) (OR: 1.46, 95% CI: 1.30-1.64), or high cumulative CPE scores (>14) (OR: 1.34, 95% CI: 1.14-1.57) had a high risk of neurological diseases. When stratified by classes of ART drugs, patients with low cumulative DDDs or low adherence had a high risk of neurological diseases, including NRTIs, PIs, NNRTIs, INSTIs, and multi-drug tablets. Subgroup analyses also suggested that patients with low cumulative DDDs or low adherence had a high risk of neurological diseases when they had high cumulative CPE scores. Patients with high cumulative DDDs or medication adherence were protected against neurological diseases only when they had low cumulative CPE scores (≤ 14). Patients may be at risk for neurological diseases when they have low cumulative DDDs, low adherence, or usage with high cumulative CPE scores. Continuous usage and low cumulative CPE scores of ART drugs may benefit neurocognitive health in patients with HIV/AIDS.

Effect of Chinese herbal medicine therapy on risks of all-cause mortality, infections, parasites, and circulatory-related mortality in HIV/AIDS patients with neurological diseases.

Introduction: Long-term living with human immunodeficiency virus (HIV) and/or antiretroviral therapy (ART) is associated with various adverse effects, including neurocognitive impairment. Heterogeneous neurocognitive impairment remains an important issue, affecting between 15-65% of human immunodeficiency virus infection and acquired immunodeficiency syndrome (HIV/AIDS) patients and resulting in work performance, safety, and health-related outcomes that have a heavy economic burden. Methods: We identified 1,209 HIV/AIDS patients with neurological diseases during 2010-2017. The Kaplan-Meier method, log-rank test, and Cox proportional hazards model were used to analyze 308 CHM users and 901 non-CHM users within this population. Major CHM clusters were determined using association rule mining and network analysis. Results and Discussion: Results showed that CHM users had a 70% lower risk of all-cause mortality (adjusted hazard ratio (aHR) = 0.30, 95% confidence interval (CI):0.16-0.58, p < 0.001) (p = 0.0007, log-rank test). Furthermore, CHM users had an 86% lower risk of infections, parasites, and circulatory-related mortality (aHR = 0.14, 95% confidence interval (CI):0.04-0.46, p = 0.001) (p = 0.0010, log-rank test). Association rule mining and network analysis showed that two CHM clusters were important for patients with neurological diseases. In the first CHM cluster, Huang Qin (HQ; root of Scutellaria baicalensis Georgi), Gan Cao (GC; root of Glycyrrhiza uralensis Fisch.), Huang Lian (HL; root of Coptis chinensis Franch.), Jie Geng (JG; root of Platycodon grandiflorus (Jacq.) A.DC.), and Huang Bai (HB; bark of Phellodendron amurense Rupr.) were identified as important CHMs. Among them, the strongest connection strength was identified between the HL and HQ. In the second CHM cluster, Suan-Zao-Ren-Tang (SZRT) and Ye Jiao Teng (YJT; stem of Polygonum multiflorum Thunb.) were identified as important CHMs with the strongest connection strength. CHMs may thus be effective in treating HIV/AIDS patients with neurological diseases, and future clinical trials are essential for the prevention of neurological dysfunction in the population.

Efficacy of HMJ-38, a new quinazolinone analogue, against the gemcitabine-resistant MIA-PaCa-2 pancreatic cancer cells.

Gemcitabine is frequently utilized to treat pancreatic cancer. The purpose of our study was to create a gemcitabine-resistant MIA-PaCa-2 pancreatic cancer cell line (MIA-GR100) and to evaluate the anti-pancreatic cancer efficacy of HMJ-38, a new quinazolinone analogue. Compared to their parental counterparts, MIA-PaCa-2, established MIA-GR100 cells were less sensitive to gemcitabine. MIA-GR100 cell viability was not affected by 10, 50 and 100 nM gemcitabine concentrations. HMJ-38 reduced MIA-GR100 cell growth and induced autophagy and apoptosis. When stained with monodansylcadaverine (MDC), acridine orange (AO), and terminal deoxynucleotide transferase dUTP nick end labeling (TUNEL), MIA-GR100 cells shrunk, punctured their membranes, and produced autophagy vacuoles and apoptotic bodies. Combining chloroquine (CQ) and 3-methyladenine (3-MA) with HMJ-38 dramatically reduced cell viability, indicating that autophagy function as a cytoprotective mechanism. MIA-GR100 cells treated with both z-VAD-FMK and HMJ-38 were much more viable than those treated with HMJ-38 alone. HMJ-38 promotes apoptosis in MIA-GR100 cells by activating caspases. Epidermal growth factor receptor (EGFR) is one of HMJ-38's principal targets, as determined via in silico target screening with network prediction. HMJ-38 also inhibited EGFR kinase activity and EGFR-associated signaling in MIA-GR100 cells. HMJ-38 may be an effective chemotherapeutic adjuvant for gemcitabine-resistant pancreatic cancer cells, in which it induces an antitumor response.

Comprehensive characterization of pharmacogenes in a Taiwanese Han population.

Pharmacogenetic (PGx) testing has not been well adopted in current clinical practice. The phenotypic distribution of clinically relevant pharmacogenes remains to be fully characterized in large population cohorts. In addition, no study has explored actionable PGx alleles in the East Asian population at a large scale. This study comprehensively analyzed 14 actionable pharmacogene diplotypes and phenotypes in 172,854 Taiwanese Han individuals by using their genotype data. Furthermore, we analyzed data from electronic medical records to investigate the effect of the actionable phenotypes on the individuals. The PGx phenotype frequencies were comparable between our cohort and the East Asian population. Overall, 99.9% of the individuals harbored at least one actionable PGx phenotype, and 29% of them have been prescribed a drug to which they may exhibit an atypical response. Our findings can facilitate the clinical application of PGx testing and the optimization of treatment and dosage individually.

Genome-wide association study of hyperthyroidism based on electronic medical record from Taiwan.

Excess thyroid hormones have complex metabolic effects, particularly hyperthyroidism, and are associated with various cardiovascular risk factors. Previous candidate gene studies have indicated that genetic variants may contribute to this variable response. Electronic medical record (EMR) biobanks containing clinical and genomic data on large numbers of individuals have great potential to inform the disease comorbidity development. In this study, we combined electronic medical record (EMR) -derived phenotypes and genotype information to conduct a genome-wide analysis of hyperthyroidism in a 35,009-patient cohort in Taiwan. Diagnostic codes were used to identify 2,767 patients with hyperthyroidism. Our genome-wide association study (GWAS) identified 44 novel genomic risk markers in 10 loci on chromosomes 2, 6, and 14 (P < 5 × 10-14), including CTLA4, HCP5, HLA-B, POU5F1, CCHCR1, HLA-DRA, HLA-DRB9, TSHR, RPL17P3, and CEP128. We further conducted a comorbidity analysis of our results, and the data revealed a strong correlation between hyperthyroidism patients with thyroid storm and stroke. In this study, we demonstrated application of the PheWAS using large EMR biobanks to inform the comorbidity development in hyperthyroidism patients. Our data suggest significant common genetic risk factors in patients with hyperthyroidism. Additionally, our results show that sex, body mass index (BMI), and thyroid storm are associated with an increased risk of stroke in subjects with hyperthyroidism.

Polygenic Risk Score of Adolescent Idiopathic Scoliosis for Potential Clinical Use.

Adolescent idiopathic scoliosis (AIS) is a common disease causing three-dimensional spinal deformity in as many as 3% of adolescents. Development of a method that can accurately predict the onset and progression of AIS is an immediate need for clinical practice. Because the heritability of AIS is estimated as high as 87.5% in twin studies, prediction of its onset and progression based on genetic data is a promising option. We show the usefulness of polygenic risk score (PRS) for the prediction of onset and progression of AIS. We used AIS genomewide association study (GWAS) data comprising 79,211 subjects in three cohorts and constructed a PRS based on association statistics in a discovery set including 31,999 female subjects. After calibration using a validation data set, we applied the PRS to a test data set. By integrating functional annotations showing heritability enrichment in the selection of variants, the PRS demonstrated an association with AIS susceptibility (p = 3.5 × 10-40 with area under the receiver-operating characteristic [AUROC] = 0.674, sensitivity = 0.644, and specificity = 0.622). The decile with the highest PRS showed an odds ratio of as high as 3.36 (p = 1.4 × 10-10 ) to develop AIS compared with the fifth in decile. The addition of a predictive model with only a single clinical parameter (body mass index) improved predictive ability for development of AIS (AUROC = 0.722, net reclassification improvement [NRI] 0.505 ± 0.054, p = 1.6 × 10-8 ), potentiating clinical use of the prediction model. Furthermore, we found the Cobb angle (CA), the severity measurement of AIS, to be a polygenic trait that showed a significant genetic correlation with AIS susceptibility (rg = 0.6, p = 3.0 × 10-4 ). The AIS PRS demonstrated a significant association with CA. These results indicate a shared polygenic architecture between onset and progression of AIS and the potential usefulness of PRS in clinical settings as a predictor to promote early intervention of AIS and avoid invasive surgery. © 2021 American Society for Bone and Mineral Research (ASBMR).

Comparison of multiple imputation algorithms and verification using whole-genome sequencing in the CMUH genetic biobank.

A genome-wide association study (GWAS) can be conducted to systematically analyze the contributions of genetic factors to a wide variety of complex diseases. Nevertheless, existing GWASs have provided highly ethnic specific data. Accordingly, to provide data specific to Taiwan, we established a large-scale genetic database in a single medical institution at the China Medical University Hospital. With current technological limitations, microarray analysis can detect only a limited number of single-nucleotide polymorphisms (SNPs) with a minor allele frequency of >1%. Nevertheless, imputation represents a useful alternative means of expanding data. In this study, we compared four imputation algorithms in terms of various metrics. We observed that among the compared algorithms, Beagle5.2 achieved the fastest calculation speed, smallest storage space, highest specificity, and highest number of high-quality variants. We obtained 15,277,414 high-quality variants in 175,871 people by using Beagle5.2. In our internal verification process, Beagle5.2 exhibited an accuracy rate of up to 98.75%. We also conducted external verification. Our imputed variants had a 79.91% mapping rate and 90.41% accuracy. These results will be combined with clinical data in future research. We have made the results available for researchers to use in formulating imputation algorithms, in addition to establishing a complete SNP database for GWAS and PRS researchers. We believe that these data can help improve overall medical capabilities, particularly precision medicine, in Taiwan.

Functional correlations between CXCL10/IP10 gene polymorphisms and risk of Kawasaki disease.

BACKGROUND: Kawasaki disease (KD) is an acute systemic vasculitis syndrome with unknown pathogen. The immune system has been suggested to involve in the pathogenesis in KD. IP10 is a chemoattractant for initiating T-cell activation. The aim of this study was to investigate the association between genetic polymorphisms of IP10 and KD. METHODS: A total of 354 KD patients and 1,709 control subjects (709 subjects in cohort 1 and 1,000 subjects in cohort 2) were enrolled in this study. Four tagging single nucleotide polymorphisms (rs3921, rs4256246, rs4508917, and rs4386624) were chosen for genotyping. RESULTS: Our results indicated that CC genotype of rs3921 and GG genotype of rs4386624 had higher frequency in KD patients compared to control. In addition, higher plasma IP10 level was observed in CC genotype of rs3921 than CG genotype and GG genotype. C/G haplotype carriers of rs3921/rs4386624 had 5.48-fold risk for KD compared to G/C haplotype carriers. Two-locus analysis further showed the combinatorial effects of rs3921 and rs4386624 in KD susceptibility. CONCLUSIONS: This study indicated the close correlation between IP10 and the risk of Kawasaki disease.

Pharmacogenomics Study for Raloxifene in Postmenopausal Female with Osteoporosis.

Osteoporosis is characterized by decreased bone mineral density and increased risk of fracture. Raloxifene is one of the treatments of osteoporosis. However, the responses were variable among patients. Previous studies revealed that the genetic variants are involved in the regulation of treatment outcomes. To date, studies that evaluate the influence of genes across all genome on the raloxifene treatment response are still limited. In this study, a total of 41 postmenopausal osteoporosis patients under regular raloxifene treatment were included. Gene-based analysis using MAGMA was applied to investigate the genetic association with the bone mineral density response to raloxifene at the lumbar spine or femoral neck site. Results from gene-based analysis indicated several genes (GHRHR, ABHD8, and TMPRSS6) related to the responses of raloxifene. Besides, the pathways of iron ion homeostasis, osteoblast differentiation, and platelet morphogenesis were enriched which implies that these pathways might be relatively susceptible to raloxifene treatment outcome. Our study provided a novel insight into the response to raloxifene.

Genomic interrogation of familial short stature contributes to the discovery of the pathophysiological mechanisms and pharmaceutical drug repositioning.

BACKGROUND: Genetic factors, dysregulation in the endocrine system, cytokine and paracrine factors are implicated in the pathogenesis of familial short stature (FSS). Nowadays, the treatment choice for FSS is limited, with only recombinant human growth hormone (rhGH) being available. METHODS: Herein, starting from the identification of 122 genetic loci related to FSS, we adopted a genetic-driven drug discovery bioinformatics pipeline based on functional annotation to prioritize crucial biological FSS-related genes. These genes were suggested to be potential targets for therapeutics. RESULTS: We discovered five druggable subnetworks, which contained seven FSS-related genes and 17 druggable targerts. CONCLUSIONS: This study provides a valuable drug repositioning accompanied by corresponding targetable gene clusters for FSS therapy.

Genome-wide association study identifies 14 previously unreported susceptibility loci for adolescent idiopathic scoliosis in Japanese.

Adolescent idiopathic scoliosis (AIS) is the most common pediatric spinal deformity. Several AIS susceptibility loci have been identified; however, they could explain only a small proportion of AIS heritability. To identify additional AIS susceptibility loci, we conduct a meta-analysis of the three genome-wide association studies consisting of 79,211 Japanese individuals. We identify 20 loci significantly associated with AIS, including 14 previously not reported loci. These loci explain 4.6% of the phenotypic variance of AIS. We find 21 cis-expression quantitative trait loci-associated genes in seven of the fourteen loci. By a female meta-analysis, we identify additional three significant loci. We also find significant genetic correlations of AIS with body mass index and uric acid. The cell-type specificity analyses show the significant heritability enrichment for AIS in multiple cell-type groups, suggesting the heterogeneity of etiology and pathogenesis of AIS. Our findings provide insights into etiology and pathogenesis of AIS.

Association of endothelin genetic variants and hospitalized infection complications in end-stage renal disease (ESRD) patients.

BACKGROUND: Infection is the second most common cause of mortality for patients with end-stage renal disease (ESRD), accompanying with immune dysfunction. Endothelin (EDN) is known to be related to inflammation; however, it is unknown whether genetic variants of the EDN gene family are associated with increased risk of hospitalized infection events. METHODS: Nineteen tagging single-nucleotide polymorphisms (tSNPs) of the EDN gene family were selected for genotyping a cohort of 190 ESRD patients. Patient demographics were recorded, the subtypes of infection events were identified, and association analysis between the EDN genetic variants and hospitalized infection events was performed. RESULTS: In this study, 106 patients were hospitalized for infection events. The leading events were pneumonia, bacteremia, and cellulitis. The minor allele of rs260741, rs197173, and rs926632 SNPs of EDN3 were found to be associated with reduced risk of hospitalized bacteremia events. CONCLUSIONS: The minor allele of rs260741, rs197173, and rs926632 in EDN3 were associated with reduced risk of hospitalized bacteremia events in ESRD patients.

Meta-Analysis of Genome-Wide Association Studies Identifies Three Loci Associated With Stiffness Index of the Calcaneus.

The stiffness index (SI) from quantitative ultrasound measurements is a good indicator of BMD and may be used to predict the risk of osteoporotic fracture. We conducted a genomewide association study (GWAS) for SI using 7742 individuals from the Taiwan Biobank, followed by a replication study in a Korean population (n = 2955). Approximately 6.1 million SNPs were subjected to association analysis, and SI-associated variants were identified. We further conducted a meta-analysis of Taiwan Biobank significant SNPs with a Korean population-based cohort. Candidate genes were prioritized according to epigenetic annotations, gene ontology, protein-protein interaction, GWAS catalog, and expression quantitative trait loci analyses. Our results revealed seven significant single-nucleotide polymorphisms (SNPs) within three loci: 7q31.31, 17p13.3, and 11q14.2. Conditional analysis showed that three SNPs, rs2536195 (CPED1/WNT16), rs1231207 (SMG6), and rs4944661 (LOC10050636/TMEM135), were the most important signals within these regions. The associations for the three SNPs were confirmed in a UK Biobank estimated BMD GWAS; these three cytobands were replicated successfully after a meta-analysis with a Korean population cohort as well. However, two SNPs were not replicated. After prioritization, we identified two novel genes, RAB15 and FNTB, as strong candidates for association with SI. Our study identified three SI-associated SNPs and two novel SI-related genes. Overall, these results provide further insight into the genetic architecture of osteoporosis. Further studies in larger East Asian populations are needed. © 2019 American Society for Bone and Mineral Research.