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.

Effect of genetic variants in cell adhesion pathways on the biochemical recurrence in prostate cancer patients with radical prostatectomy.

The aberrant expression of cell adhesion molecules is a hallmark of epithelial-to-mesenchymal transition, resulting in the transformation of cancer cells to a more aggressive phenotype. This study investigated the association between genetic variants in cell adhesion pathways and the prognosis of patients with prostate cancer following radical prostatectomy (RP). A total of 18 haplotype-tagging single-nucleotide polymorphisms (SNPs) in eight cancer-related adhesion molecules were genotyped in 458 prostate cancer patients, followed by the replication of the top SNPs in an additional set of 185 patients. Log-rank test and multivariate Cox regression analysis adjusted for covariates were used to evaluate associations with the risk of biochemical recurrence (BCR) after RP. In the discovery set, four SNPs in CDH2 were marginally associated with BCR. Among these, CDH2 rs643555C > T was found to be associated with BCR in the replication set. Patients with rs643555TT genotype had a significantly shorter BCR-free survival compared with those with CC/CT genotypes in the combined analysis (adjusted hazard ratio 1.78, 95% confidence interval 1.19-2.67, P = 0.005). Additional analyses revealed that rs643555T was associated with higher expression of CDH2, and upregulated CDH2 was correlated with tumor aggressiveness and shortened BCR-free survival. In conclusion, rs643555C > T affects CDH2 expression, and thus influences BCR in localized prostate cancer patients treated with RP. CDH2 rs643555 may be a promising biomarker to identify patients at high risk of poor prostate cancer prognosis.

Determining population stratification and subgroup effects in association studies of rare genetic variants for nicotine dependence.

BACKGROUND: Rare variants (minor allele frequency < 1% or 5 %) can help researchers to deal with the confounding issue of 'missing heritability' and have a proven role in dissecting the etiology for human diseases and complex traits. METHODS: We extended the combined multivariate and collapsing (CMC) and weighted sum statistic (WSS) methods and accounted for the effects of population stratification and subgroup effects using stratified analyses by the principal component analysis, named here as 'str-CMC' and 'str-WSS'. To evaluate the validity of the extended methods, we analyzed the Genetic Architecture of Smoking and Smoking Cessation database, which includes African Americans and European Americans genotyped on Illumina Human Omni2.5, and we compared the results with those obtained with the sequence kernel association test (SKAT) and its modification, SKAT-O that included population stratification and subgroup effect as covariates. We utilized the Cochran-Mantel-Haenszel test to check for possible differences in single nucleotide polymorphism allele frequency between subgroups within a gene. We aimed to detect rare variants and considered population stratification and subgroup effects in the genomic region containing 39 acetylcholine receptor-related genes. RESULTS: The Cochran-Mantel-Haenszel test as applied to GABRG2 (P = 0.001) was significant. However, GABRG2 was detected both by str-CMC (P= 8.04E-06) and str-WSS (P= 0.046) in African Americans but not by SKAT or SKAT-O. CONCLUSIONS: Our results imply that if associated rare variants are only specific to a subgroup, a stratified analysis might be a better approach than a combined analysis.

Effect of adiponectin level and genetic variation of its receptors on diabetic retinopathy: A case-control study.

Adiponectin (APN) and its receptors have been reported to be associated with metabolic phenotypes. To better understand the effects of APN levels and its receptors on diabetic retinopathy (DR), we investigated the association of the plasma APN level and variations in APN-related genes with DR, individually and in combination.Patients with type 2 diabetes (T2D; N = 1604), above 20 years of age from the Taiwanese population participated in the study. Demographic information, blood pressure, and serological markers were recorded at enrollment. Genomic DNA was isolated and genotyped. The plasma APN levels were measured by enzyme-linked immunosorbent assay.T2D patients with DR (N = 632) had diabetes for a longer duration, and had higher HbA1c, and systolic and diastolic blood pressure compared to those without DR (N = 972) (P < .001, for all the parameters). Overall, 10 single nucleotide polymorphisms (SNPs) in ADIPOQ and CDH13 susceptibility loci were associated with DR. Gene risk score (GRS) was calculated based on 10 SNPs for each subject and the cumulative effect of genes was observed. Among the subjects with plasma APN level (N = 518), natural logarithm (LN) of APN (LN [APN]; odds ratio [OR] = 1.63, 95% confidence interval [CI] = 1.19-2.25) and GRS (OR = 1.90, 95% CI = 1.11-3.26 for middle range of GRS, and OR = 2.61, 95% CI = 1.48-4.59 for high range of GRS) were independent risk factors for DR after adjustment for other parameters.In conclusion, the plasma APN level and the genetic variations in adiponectin receptors were associated with DR.