Genomic insights into mRNA COVID-19 vaccines efficacy: Linking genetic polymorphisms to waning immunity.

Genetic polymorphisms have been linked to the differential waning of vaccine-induced immunity against COVID-19 following vaccination. Despite this, evidence on the mechanisms behind this waning and its implications for vaccination policy remains limited. We hypothesize that specific gene variants may modulate the development of vaccine-initiated immunity, leading to impaired immune function. This study investigates genetic determinants influencing the sustainability of immunity post-mRNA vaccination through a genome-wide association study (GWAS). Utilizing a hospital-based, test negative case-control design, we enrolled 1,119 participants from the Taiwan Precision Medicine Initiative (TPMI) cohort, all of whom completed a full mRNA COVID-19 vaccination regimen and underwent PCR testing during the Omicron outbreak. Participants were classified into breakthrough and protected groups based on PCR results. Genetic samples were analyzed using SNP arrays with rigorous quality control. Cox regression identified significant single nucleotide polymorphisms (SNPs) associated with breakthrough infections, affecting 743 genes involved in processes such as antigenic protein translation, B cell activation, and T cell function. Key genes identified include CD247, TRPV1, MYH9, CCL16, and RPTOR, which are vital for immune responses. Polygenic risk score (PRS) analysis revealed that individuals with higher PRS are at greater risk of breakthrough infections post-vaccination, demonstrating a high predictability (AUC = 0.787) in validating population. This finding confirms the significant influence of genetic variations on the durability of immune responses and vaccine effectiveness. This study highlights the importance of considering genetic polymorphisms in evaluating vaccine-induced immunity and proposes potential personalized vaccination strategies by tailoring regimens to individual genetic profiles.

Biomedical literature mining: graph kernel-based learning for gene-gene interaction extraction.

The supervised machine learning method is often used for biomedical relationship extraction. The disadvantage is that it requires much time and money to manually establish an annotated dataset. Based on distant supervision, the knowledge base is combined with the corpus, thus, the training corpus can be automatically annotated. As many biomedical databases provide knowledge bases for study with a limited number of annotated corpora, this method is practical in biomedicine. The clinical significance of each patient's genetic makeup can be understood based on the healthcare provider's genetic database. Unfortunately, the lack of previous biomedical relationship extraction studies focuses on gene-gene interaction. The main purpose of this study is to develop extraction methods for gene-gene interactions that can help explain the heritability of human complex diseases. This study referred to the information on gene-gene interactions in the KEGG PATHWAY database, the abstracts in PubMed were adopted to generate the training sample set, and the graph kernel method was adopted to extract gene-gene interactions. The best assessment result was an F1-score of 0.79. Our developed distant supervision method automatically finds sentences through the corpus without manual labeling for extracting gene-gene interactions, which can effectively reduce the time cost for manual annotation data; moreover, the relationship extraction method based on a graph kernel can be successfully applied to extract gene-gene interactions. In this way, the results of this study are expected to help achieve precision medicine.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Genome-Wide and Candidate Gene Association Analyses Identify a 14-SNP Combination for Hypertension in Patients With Type 2 Diabetes.

BACKGROUND: High blood pressure is common and comorbid with type 2 diabetes (T2D). Almost 50% of patients with T2D have high blood pressure. Patients with both conditions of hypertension (HTN) and T2D are at risk for cardiovascular diseases and mortality. The study aim was to investigate genetic risk factors for HTN in T2D patients. METHODS: This study included 999 T2D (cohort 1) patients for the first genome scan stage and 922 T2D (cohort 2) patients for the replication stage. Here, we investigated the genetic susceptibility and cumulative weighted genetic risk score for HTN in T2D patients of Han Chinese descent in Taiwan. RESULTS: Thirty novel genetic single nucleotide polymorphisms (SNPs) were associated with HTN in T2D after adjusting for age and body mass index (P value <1 × 10-4). Eight blood pressure-related and/or HTN-related genetic SNPs were associated with HTN in T2D after adjusting for age and body mass index (P value <0.05). Linkage disequilibrium and cumulative weighted genetic risk score analyses showed that 14 of the 38 SNPs were associated with risk of HTN in a dose-dependent manner in T2D (Cochran-Armitage trend test: P value <0.0001). The 14-SNP cumulative weighted genetic risk score was also associated with increased regression tendency of systolic blood pressure in T2D (SBP = 122.05 + 0.8 × weighted genetic risk score; P value = 0.0001). CONCLUSIONS: A cumulative weighted genetic risk score composed of 14 SNPs is important for HTN, increased tendency of systolic blood pressure, and may contribute to HTN risk in T2D in Taiwan.

Maximal Segmental Score Method for Localizing Recessive Disease Variants Based on Sequence Data.

BACKGROUND: Due to the affordability of whole-genome sequencing, the genetic association design can now address rare diseases. However, some common statistical association methods only consider homozygosity mapping and need several criteria, such as sliding windows of a given size and statistical significance threshold setting, such as P-value < 0.05 to achieve good power in rare disease association detection. METHODS: Our region-specific method, called expanded maximal segmental score (eMSS), converts p-values into continuous scores based on the maximal segmental score (MSS) (Lin et al., 2014) for detecting disease-associated segments. Our eMSS considers the whole genome sequence data, not only regions of homozygosity in candidate genes. Unlike sliding window methods of a given size, eMSS does not need predetermined parameters, such as window size or minimum or maximum number of SNPs in a segment. The performance of eMSS was evaluated by simulations and real data analysis for autosomal recessive diseases multiple intestinal atresia (MIA) and osteogenesis imperfecta (OI), where the number of cases is extremely small. For the real data, the results by eMSS were compared with a state-of-the-art method, HDR-del (Imai et al., 2016). RESULTS: Our simulation results show that eMSS had higher power as the number of non-causal haplotype blocks decreased. The type I error for eMSS under different scenarios was well controlled, p < 0.05. For our observed data, the bone morphogenetic protein 1 (BMP1) gene on chromosome 8, the Violaxanthin de-epoxidase-related chloroplast (VDR) gene on chromosome 12 associated with OI, and the tetratricopeptide repeat domain 7A (TTC7A) gene on chromosome 2 associated with MIA have previously been identified as harboring the relevant pathogenic mutations. CONCLUSIONS: When compared to HDR-del, our eMSS is powerful in analyzing even small numbers of recessive cases, and the results show that the method can further reduce numbers of candidate variants to a very small set of susceptibility pathogenic variants underlying OI and MIA. When we conduct whole-genome sequence analysis, eMSS used 3/5 the computation time of HDR-del. Without additional parameters needing to be set in the segment detection, the computational burden for eMSS is lower compared with that in other region-specific approaches.

Genetic Architecture Associated With Familial Short Stature.

CONTEXT: Human height is an inheritable, polygenic trait under complex and multilocus genetic regulation. Familial short stature (FSS; also called genetic short stature) is the most common type of short stature and is insufficiently known. OBJECTIVE: To investigate the FSS genetic profile and develop a polygenic risk predisposition score for FSS risk prediction. DESIGN AND SETTING: The FSS participant group of Han Chinese ancestry was diagnosed by pediatric endocrinologists in Taiwan. PATIENTS AND INTERVENTIONS: The genetic profiles of 1163 participants with FSS were identified by using a bootstrapping subsampling and genome-wide association studies (GWAS) method. MAIN OUTCOME MEASURES: Genetic profile, polygenic risk predisposition score for risk prediction. RESULTS: Ten novel genetic single nucleotide polymorphisms (SNPs) and 9 reported GWAS human height-related SNPs were identified for FSS risk. These 10 novel SNPs served as a polygenic risk predisposition score for FSS risk prediction (area under the curve: 0.940 in the testing group). This FSS polygenic risk predisposition score was also associated with the height reduction regression tendency in the general population. CONCLUSION: A polygenic risk predisposition score composed of 10 genetic SNPs is useful for FSS risk prediction and the height reduction tendency. Thus, it might contribute to FSS risk in the Han Chinese population from Taiwan.

Lack of association of genetic variants for diabetic retinopathy in Taiwanese patients with diabetic nephropathy.

OBJECTIVE: Diabetic nephropathy (DN) and diabetic retinopathy (DR) comprise major microvascular complications of diabetes that occur with a high concordance rate in patients and are considered to potentially share pathogeneses. In this case-control study, we sought to investigate whether DR-related single nucleotide polymorphisms (SNPs) exert pleiotropic effects on renal function outcomes among patients with diabetes. RESEARCH DESIGN AND METHODS: A total of 33 DR-related SNPs were identified by replicating published SNPs and via a genome-wide association study. Furthermore, we assessed the cumulative effects by creating a weighted genetic risk score and evaluated the discriminatory and prediction ability of these genetic variants using DN cases according to estimated glomerular filtration rate (eGFR) status along with a cohort with early renal functional decline (ERFD). RESULTS: Multivariate logistic regression models revealed that the DR-related SNPs afforded no individual or cumulative genetic effect on the nephropathy risk, eGFR status or ERFD outcome among patients with type two diabetes in Taiwan. CONCLUSION: Our findings indicate that larger studies would be necessary to clearly ascertain the effects of individual genetic variants and further investigation is also required to identify other genetic pathways underlying DN.