Comprehensive human leukocyte antigen genotyping of patients with type 1 diabetes mellitus in Taiwan.

BACKGROUND: Type 1 diabetes (T1D) mellitus is an autoimmune disorder involving both complex genetic and environmental factors. The incidence rates are low in Asian countries, and the specific, explanatory genetic factors underlying this have been investigated. The aim of this study was to elucidate the association of human leukocyte antigen (HLA) alleles/haplotypes with T1D in Taiwan. METHODS: We performed direct comprehensive genotyping of 6 classical HLA loci (HLA-A, -B, -C, -DPB1, -DQB1, and -DRB1) to 4-digit resolution in 104 unrelated T1D patients and 504 controls. Twenty-four of the 104 patients also exhibited thyroid autoimmunity. RESULTS: Three major susceptibility haplotypes were identified: DRB1*03:01-DQB1*02:01 (odds ratio [OR] = 5.39 under the dominant model, P = 2.3 × 10-13 ), DRB1*04:05-DQB1*04:01 (OR = 2.44, P = 5.0 × 10-4 ), and DRB1*09:01-DQB1*03:03 (OR = 2.02, P = 1.4 × 10-3 ); one protective haplotype was identified: DRB1*08:03-DQB1*06:01 (OR = 0.10, P = 1.6 × 10-3 ). DRB1*03:01-DQB1*02:01, the major T1D susceptibility haplotype, was found at a lower frequency in T1D patients with thyroid autoimmunity. The T1D protective allele DRB1*12:02 was shown to be protective against Graves' disease in our previous report. CONCLUSION: In addition to clarifying the roles of several known T1D HLA alleles and haplotypes, we discovered that the DRB1*08:03-DQB1*06:01 haplotype is protective against T1D. The DRB1*12:02 allele protected against both T1D and Graves' disease.

Susceptible genes of restless legs syndrome in migraine.

Objective Several genetic variants have been found to increase the risk of restless legs syndrome (RLS). The aim of the present study was to determine if these genetic variants were also associated with the comorbidity of RLS and migraine in patients. Methods Thirteen single-nucleotide polymorphisms (SNPs) at six RLS risk loci ( MEIS1, BTBD9, MAP2K5, PTPRD, TOX3, and an intergenic region on chromosome 2p14) were genotyped in 211 migraine patients with RLS and 781 migraine patients without RLS. Association analyses were performed for the overall cohort, as well as for the subgroups of patients who experienced migraines with and without aura and episodic migraines (EMs) vs. chronic migraines (CMs). In order to verify which genetic markers were potentially related to the incidence of RLS in migraine patients, multivariate regression analyses were also performed. Results Among the six tested loci, only MEIS1 was significantly associated with RLS. The most significant SNP of MEIS1, rs2300478, increased the risk of RLS by 1.42-fold in the overall cohort ( p = 0.0047). In the subgroup analyses, MEIS1 augmented the risk of RLS only in the patients who experienced EMs (odds ratio (OR) = 1.99, p = 0.0004) and not those experiencing CMs. Multivariate regression analyses further showed that rs2300478 in MEIS1 (OR = 1.39, p = 0.018), a CM diagnosis (OR = 1.52, p = 0.022), and depression (OR = 1.86, p = 0.005) were independent predictors of RLS in migraine. Conclusions MEIS1 variants were associated with an increased risk of RLS in migraine patients. It is possible that an imbalance in iron homeostasis and the dopaminergic system may represent a link between RLS incidence and migraines.

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

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

DRD2 haplotype associated with negative symptoms and sustained attention deficits in Han Chinese with schizophrenia in Taiwan.

Previous studies have reported significant associations between schizophrenia and the dopamine receptor D2 gene (DRD2) variants. The relationships between DRD2 and clinical phenotypes are of particular interest because DRD2 has been shown to associate with treatment response and prefrontal dopamine transmission. Glatt et al. reported significant associations between schizophrenia and DRD2 variants (two single-nucleotide polymorphisms (SNPs) rs1079727 and rs2283265, and two haplotypes, block 3 (rs1079727(A)-rs2440390(C)-rs2283265(G)) and block 4 (rs1801028(G)-rs1110977(A)-rs1124492(C)-rs2734841 (T))) in 2408 Han Chinese individuals in Taiwan. To further investigate the relationships between the SNPs/haplotypes of DRD2 and clinical symptoms and neuropsychological function, we compared the quantitative phenotypes in patients with risk alleles/haplotypes and those without. The results showed that the A allele of rs1079727, G allele of rs2283265, A allele of rs1124492 and the risk haplotype (A-C-G) of block 3 were associated with more severe negative symptoms. With regard to neuropsychological performance, the risk haplotype (G-A-C-T) of block 4 was associated with poorer performance in the sustained attention task. Our results imply that the genetic variants of DRD2 might not only have a role in susceptibility to schizophrenia, but also influence the phenotypes of negative symptoms and sustained attention in schizophrenia. This association warrants further validation.

Accelerated cardiovascular risk after viral clearance in hepatitis C patients with the NAMPT-rs61330082 TT genotype: An 8-year prospective cohort study.

Involvement of extracellular nicotinamide phosphoribosyltransferase (eNAMPT, i.e., visfatin or pre-B-cell colony-enhancing factor), a cancer metabokine, in chronically hepatitis C virus (HCV)-infected (CHC) patients with sustained virological responses (SVRs) remains elusive. This 8-year prospective cohort study evaluated eNAMPT profiles of 842 consecutive CHC patients, including 519 who had completed an anti-HCV therapy course and pre-therapy and 24-week post-therapy surveys. For 842 patients, pre-therapy associations were HCV RNA, homeostatic model assessment for insulin resistance (HOMA-IR) index, and body mass index with eNAMPT levels, and NAMPT-rs61330082 T allele with total cholesterol levels. NAMPT-rs10953502, NAMPT-rs2058539, and NAMPT-rs61330082 were in a linkage disequilibrium block, which was associated with total cholesterol levels. Compared to pre-therapy levels, at 24 weeks post-therapy, decreased eNAMPT and increased lipid levels were observed in SVR patients (n = 427). Among SVR patients, higher cumulative incidences of cardiovascular events occurred in those with a NAMPT-rs61330082 TT genotype than those with non-TT genotypes (28.2% vs. 8.4%, p < 0.001). NAMPT-rs61330082 TT genotype was independently associated with incident cardiovascular events (95% CI hazard ratio (HR): 1.88-10.37; HR: 4.415); no eNAMPT profiles were associated with incident malignancies. Of CHC patients, hepatic vascular endothelial cells and baseline peripheral leukocytes expressed higher eNAMPT levels than controls, and peripheral eNAMPT-positive leukocyte proportions decreased after SVR. During HCV infection, eNAMPT involvement in glucose metabolism was modulated by HCV RNA linked to lipid metabolism and NAMPT-associated SNPs. Hepatic endothelial cells and peripheral leukocytes potentially secrete eNAMPT. Caution is required for incident cardiovascular events in SVR patients with NAMPT-rs61330082 TT genotype.

Modeling expression quantitative trait loci in data combining ethnic populations.

BACKGROUND: Combining data from different ethnic populations in a study can increase efficacy of methods designed to identify expression quantitative trait loci (eQTL) compared to analyzing each population independently. In such studies, however, the genetic diversity of minor allele frequencies among populations has rarely been taken into account. Due to the fact that allele frequency diversity and population-level expression differences are present in populations, a consensus regarding the optimal statistical approach for analysis of eQTL in data combining different populations remains inconclusive. RESULTS: In this report, we explored the applicability of a constrained two-way model to identify eQTL for combined ethnic data that might contain genetic diversity among ethnic populations. In addition, gene expression differences resulted from ethnic allele frequency diversity between populations were directly estimated and analyzed by the constrained two-way model. Through simulation, we investigated effects of genetic diversity on eQTL identification by examining gene expression data pooled from normal quantile transformation of each population. Using the constrained two-way model to reanalyze data from Caucasians and Asian individuals available from HapMap, a large number of eQTL were identified with similar genetic effects on the gene expression levels in these two populations. Furthermore, 19 single nucleotide polymorphisms with inter-population differences with respect to both genotype frequency and gene expression levels directed by genotypes were identified and reflected a clear distinction between Caucasians and Asian individuals. CONCLUSIONS: This study illustrates the influence of minor allele frequencies on common eQTL identification using either separate or combined population data. Our findings are important for future eQTL studies in which different datasets are combined to increase the power of eQTL identification.

Effects of sex and generation on hepatitis B viral load in families with hepatocellular carcinoma.

AIM: To explore factors associated with persistent hepatitis B virus (HBV) infection in a cohort of hepatocellular carcinoma (HCC)-affected families and then investigate factors that correlate with individual viral load among hepatitis B surface antigen (HBsAg)-positive relatives. METHODS: We evaluated non-genetic factors associated with HBV replication in relatives of patients with HCC. Relatives of 355 HCC cases were interviewed using a structured questionnaire. Demographics, relationship to index case, HBsAg status of mothers and index cases were evaluated for association with the HBV persistent infection or viral load by generalized estimating equation analysis. RESULTS: Among 729 relatives enrolled, parent generation (P = 0.0076), index generation (P = 0.0044), mothers positive for HBsAg (P = 0.0007), and HBsAg-positive index cases (P = 5.98 × 10-8) were associated with persistent HBV infection. Factors associated with HBV viral load were evaluated among 303 HBsAg-positive relatives. Parent generation (P = 0.0359) and sex (P = 0.0007) were independent factors associated with HBV viral load. The intra-family HBV viral load was evaluated in families clustered with HBsAg-positive siblings. An intra-family trend of similar HBV viral load was found for 27 of 46 (58.7%) families. Male offspring of HBsAg-positive mothers (P = 0.024) and older siblings were associated with high viral load. CONCLUSION: Sex and generation play important roles on HBV viral load. Maternal birth age and nutritional changes could be the reasons of viral load difference between generations.

Association of circadian genes with diurnal blood pressure changes and non-dipper essential hypertension: a genetic association with young-onset hypertension.

Recent studies have suggested that circadian genes have important roles in maintaining the circadian rhythm of the cardiovascular system. However, the associations between diurnal BP changes and circadian genes remain undetermined. We conducted a genetic association study of young-onset hypertension, in which 24-h ambulatory blood pressure (BP) monitoring was performed. A total of 23 tag single-nucleotide polymorphisms (SNPs) on 11 genes involved in circadian rhythms were genotyped for correlations with diurnal BP variation phenotypes. A permutation test was used to correct for multiple testing. Five tag SNPs within five loci, including rs3888170 in NPAS2, rs6431590 in PER2, rs1410225 in RORßß, rs3816358 in BMAL1 and rs10519096 in RORα, were significantly associated with the non-dipper phenotype in 372 young hypertensive patients. A genetic risk score was generated by counting the risk alleles and effects for each individual. Genotyping was performed in an additional independent set of 619 young-onset hypertensive subjects. Altogether, non-dippers had a higher weighted genetic risk score than dippers (1.67±0.56 vs. 1.54±0.55, P<0.001), and the additive genetic risk score also indicated a graded association with decreased diurnal BP changes (P=0.006), as well as a non-dipper phenotype (P=0.031). After multivariable logistic analysis, only the circadian genetic risk score (odds ratio (OR), 1550; 95% confidence interval (CI), 1.225-1.961, P<0.001) and the use of ß-blockers (OR, 1.519; 95% CI, 1.164-1.982, P=0.003) were independently associated with the presence of non-dippers among subjects with young-onset hypertension. Genetic variants in circadian genes were associated with the diurnal phenotype of hypertension, suggesting a genetic association with diurnal BP changes in essential hypertension.

ANXA7, PPP3CB, DNAJC9, and ZMYND17 genes at chromosome 10q22 associated with the subgroup of schizophrenia with deficits in attention and executive function.

BACKGROUND: A genome scan of Taiwanese schizophrenia families suggested linkage to chromosome 10q22.3. We aimed to find the candidate genes in this region. METHODS: A total of 476 schizophrenia families were included. Hierarchical clustering method was used for clustering families to homogeneous subgroups according to their performances of sustained attention and executive function. Association analysis was performed using family-based association testing and TRANSMIT. Candidate associated regions were identified using the longest significance run method. The relative messenger RNA expression level was determined using real-time reverse transcriptase polymerase chain reaction. RESULTS: First, we genotyped 18 microsatellite markers between D10S1432 and D10S1239. The maximum nonparametric linkage score was 2.79 on D10S195. Through family clustering, we found the maximum nonparametric linkage score was 3.70 on D10S195 in the family cluster with deficits in attention and executive function. Second, we genotyped 79 single nucleotide polymorphisms between D10S1432 and D10S580 in 90 attention deficit and execution deficit families. Association analysis indicated significant transmission distortion for nine single nucleotide polymorphisms. Using the longest significance run method, we identified a 427-kilobase region as a significant candidate region, which encompasses nine genes. Third, we studied messenger RNA expression of these nine genes in Epstein-Barr virus-transformed lymphoblastic cells. In schizophrenic patients, there was significantly lower expression of ANXA7, PPP3CB, and DNAJC9 and significantly higher expression of ZMYND17. CONCLUSIONS: ANXA7, PPP3CB, DNAJC9, and ZMYND17 genes are potential candidate genes for schizophrenia, especially in patients with deficits in sustained attention and executive function. The responsible functional variants remained to be clarified.

Ectopic and high CXCL13 chemokine expression in myasthenia gravis with thymic lymphoid hyperplasia.

Myasthenia gravis (MG) is an antibody and complement mediated autoimmune disease. Serum CXC chemokine ligand 13 (CXCL13) was found to be elevated in MG patients and high CXCL13 level was associated with severe clinical stages, especially in females with thymic lymphoid hyperplasia. Both protein and mRNA of CXCL13 and CXC chemokine receptor 5 (CXCR5) in the thymic tissues were significantly higher in MG patients with lymphoid hyperplasia than those with thymoma. Our data indicated that serum CXCL13 can be used as an index of disease severity and ectopic thymic expression of CXCL13 might be associated with aberrant cell trafficking to the thymus of MG.

Clinical implications of chromosomal abnormalities in gastric adenocarcinomas.

Gastric carcinoma (GC) is one of the most common malignancies worldwide and has a very poor prognosis. Genetic imbalances in 62 primary gastric adenocarcinomas of various histopathologic types and pathologic stages and six gastric cancer-derived cell lines were analyzed by comparative genomic hybridization, and the relationship of genomic abnormalities to clinical features in primary GC was evaluated at a genome-wide level. Eighty-four percent of the tumors and all six cell lines showed DNA copy number changes. The recurrent chromosomal abnormalities including gains at 15 regions and losses at 8 regions were identified. Statistical analyses revealed that gains at 17q24-qter (53%), 20q13-qter (48%), 1p32-p36 (42%), 22q12-qter (27%), 17p13-pter (24%), 16p13-pter (21%), 6p21-pter (19%), 20p12-pter (19%), 7p21-pter (18%), 3q28-qter (8%), and 13q13-q14 (8%), and losses at 18q12-qter (11%), 3p12 (8%), 3p25-pter (8%), 5q14-q23 (8%), and 9p21-p23 (5%), are associated with unique patient or tumor-related features. GCs of differing histopathologic features were shown to be associated with distinct patterns of genetic alterations, supporting the notion that they evolve through distinct genetic pathways. Metastatic tumors were also associated with specific genetic changes. These regions may harbor candidate genes involved in the pathogenesis of this malignancy.