Further development of multiplex single nucleotide polymorphism typing method, the DigiTag2 assay

A common variant rs2054564 in ADAMST17 is associated with susceptibility to lumbar spondylosis

The molecular pathophysiology underlying lumbar spondylosis development remains unclear. To identify genetic factors associated with lumbar spondylosis, we conducted a genome-wide association study using 83 severe lumbar spondylosis cases and 182 healthy controls and identified 65 candidate disease-associated single nucleotide polymorphisms (SNPs). Replication analysis in 510 case and 911 control subjects from five independent Japanese cohorts identified rs2054564, located in intron 7 of ADAMTS17, as a disease-associated SNP with a genome-wide significance threshold (P = 1.17 × 10-11, odds ratio = 1.92). This association was significant even after adjustment of age, sex, and body mass index (P = 7.52 × 10-11). A replication study in a Korean cohort, including 123 case and 319 control subjects, also verified the significant association of this SNP with severe lumbar spondylosis. Immunohistochemistry revealed that fibrillin-1 (FBN1) and ADAMTS17 were co-expressed in the annulus fibrosus of intervertebral discs (IVDs). ADAMTS17 overexpression in MG63 cells promoted extracellular microfibrils biogenesis, suggesting the potential role of ADAMTS17 in IVD function through interaction with fibrillin fibers. Finally, we provided evidence of FBN1 involvement in IVD function by showing that lumbar IVDs in patients with Marfan syndrome, caused by heterozygous FBN1 gene mutation, were significantly more degenerated. We identified a common SNP variant, located in ADAMTS17, associated with susceptibility to lumbar spondylosis and demonstrated the potential role of the ADAMTS17-fibrillin network in IVDs in lumbar spondylosis development.

Genome-wide association study of idiopathic hypersomnia in a Japanese population

Idiopathic hypersomnia (IH) is a rare sleep disorder characterized by excessive daytime sleepiness, great difficulty upon awakening, and prolonged sleep time. In contrast to narcolepsy type 1, which is a well-recognized hypersomnia, the etiology of IH remains poorly understood. No susceptibility loci for IH have been identified, although familial aggregations have been observed among patients with IH. Narcolepsy type 1 is strongly associated with human leukocyte antigen (HLA)-DQB1*06:02; however, no significant associations between IH and HLA alleles have been reported. To identify genetic variants that affect susceptibility to IH, we performed a genome-wide association study (GWAS) and two replication studies involving a total of 414 Japanese patients with IH and 6587 healthy Japanese individuals. A meta-analysis of the three studies found no single-nucleotide polymorphisms (SNPs) that reached the genome-wide significance level. However, we identified several candidate SNPs for IH. For instance, a common genetic variant (rs2250870) within an intron of PDE9A was suggestively associated with IH. rs2250870 was significantly associated with expression levels of PDE9A in not only whole blood but also brain tissues. The leading SNP in the PDE9A region was the same in associations with both IH and PDE9A expression. PDE9A is a potential target in the treatment of several brain diseases, such as depression, schizophrenia, and Alzheimer's disease. It will be necessary to examine whether PDE9A inhibitors that have demonstrated effects on neurophysiologic and cognitive function can contribute to the development of new treatments for IH, as higher expression levels of PDE9A were observed with regard to the risk allele of rs2250870. The present study constitutes the first GWAS of genetic variants associated with IH. A larger replication study will be required to confirm these associations.

Supplementary Information

The online version contains supplementary material available at 10.1007/s41105-021-00349-2.

Associations between peroxisome proliferator-activated receptor γ (PPAR-γ) polymorphisms and serum lipids: Two cross-sectional studies of community-dwelling adults

Dyslipidemia is a well-established risk factor for cardiovascular disease. Experimental studies have reported that peroxisome proliferator-activated receptor γ (PPAR-γ) regulates adipocyte differentiation, lipid storage, and glucose metabolism. Therefore, we examined the associations between PPAR-γ polymorphisms (rs1801282, rs3856806, rs12497191, rs1151999, and rs1152003) and serum lipids in two cross-sectional studies. In the Shizuoka area of the Japan Multi-Institutional Collaborative Cohort Study, we examined 4,952 participants (3,356 men and 1,596 women) in a baseline survey and 2,245 participants (1,550 men and 695 women) in a second survey 5 years later. Outcome measures were the prevalence of dyslipidemia (low-density lipoprotein-cholesterol [LDL-C] ≥ 140 mg/dl, high-density lipoprotein-cholesterol < 40 mg/dl, triglycerides ≥ 150 mg/dl, and/or use of cholesterol-lowering drugs) and the prevalence of high LDL-C (LDL-C ≥ 140 mg/dl and/or use of cholesterol-lowering drugs). Multivariate odds ratios (ORs) were estimated by using unconditional logistic regression models. A total of 2,114 and 1,431 individuals (42.7% and 28.9%) had dyslipidemia and high LDL-C in the baseline survey, respectively, as did 933 and 716 (41.6% and 31.9%), respectively, in the second survey. In the baseline study, compared with major allele homozygotes, minor allele homozygotes of rs3856806 and rs12497191 had a 42% (OR, 0.58; 95% confidence interval (CI), 0.39-0.85) and 23% (OR, 0.77; 95% CI, 0.60-0.99) lower risk of dyslipidemia, respectively, after adjustment for potential confounding factors. In addition, minor allele homozygotes of rs3856806 had a 45% (OR, 0.55; 95% CI, 0.35-0.86) lower risk of high LDL-C. Similar risk reductions were found in the second survey. In conclusion, rs3856806 and rs12497191 polymorphisms may be related to a lower risk of dyslipidemia and high LDL-C.

One Single Nucleotide Polymorphism of the TRPM2 Channel Gene Identified as a Risk Factor in Bipolar Disorder Associates with Autism Spectrum Disorder in a Japanese Population

The transient receptor potential melastatin 2 (TRPM2) is a non-specific cation channel, resulting in Ca2+ influx at warm temperatures from 34 °C to 47 °C, thus including the body temperature range in mammals. TRPM2 channels are activated by β-NAD+, ADP-ribose (ADPR), cyclic ADPR, and 2'-deoxyadenosine 5'-diphosphoribose. It has been shown that TRPM2 cation channels and CD38, a type II or type III transmembrane protein with ADP-ribosyl cyclase activity, simultaneously play a role in heat-sensitive and NAD+ metabolite-dependent intracellular free Ca2+ concentration increases in hypothalamic oxytocinergic neurons. Subsequently, oxytocin (OT) is released to the brain. Impairment of OT release may induce social amnesia, one of the core symptoms of autism spectrum disorder (ASD). The risk of single nucleotide polymorphisms (SNPs) and variants of TRPM2 have been reported in bipolar disorder, but not in ASD. Therefore, it is reasonable to examine whether SNPs or haplotypes in TRPM2 are associated with ASD. Here, we report a case-control study with 147 ASD patients and 150 unselected volunteers at Kanazawa University Hospital in Japan. The sequence-specific primer-polymerase chain reaction method together with fluorescence correlation spectroscopy was applied. Of 14 SNPs examined, one SNP (rs933151) displayed a significant p-value (OR = 0.1798, 95% CI = 0.039, 0.83; Fisher's exact test; p = 0.0196). The present research data suggest that rs93315, identified as a risk factor for bipolar disorder, is a possible association factor for ASD.

Strong Association of the HLA-DR/DQ Locus with Childhood Steroid-Sensitive Nephrotic Syndrome in the Japanese Population

Background Nephrotic syndrome is the most common cause of chronic glomerular disease in children. Most of these patients develop steroid-sensitive nephrotic syndrome (SSNS), but the loci conferring susceptibility to childhood SSNS are mainly unknown.Methods We conducted a genome-wide association study (GWAS) in the Japanese population; 224 patients with childhood SSNS and 419 adult healthy controls were genotyped using the Affymetrix Japonica Array in the discovery stage. Imputation for six HLA genes (HLA-A, -C, -B, -DRB1, -DQB1, and -DPB1) was conducted on the basis of Japanese-specific references. We performed genotyping for HLA-DRB1/-DQB1 using a sequence-specific oligonucleotide-probing method on a Luminex platform. Whole-genome imputation was conducted using a phased reference panel of 2049 healthy Japanese individuals. Replication was performed in an independent Japanese sample set including 216 patients and 719 healthy controls. We genotyped candidate single-nucleotide polymorphisms using the DigiTag2 assay.Results The most significant association was detected in the HLA-DR/DQ region and replicated (rs4642516 [minor allele G], combined P_allelic=7.84×10^-23; odds ratio [OR], 0.33; 95% confidence interval [95% CI], 0.26 to 0.41; rs3134996 [minor allele A], combined P_allelic=1.72×10^-25; OR, 0.29; 95% CI, 0.23 to 0.37). HLA-DRB1*08:02 (Pc=1.82×10^-9; OR, 2.62; 95% CI, 1.94 to 3.54) and HLA-DQB1*06:04 (Pc=2.09×10^-12; OR, 0.10; 95% CI, 0.05 to 0.21) were considered primary HLA alleles associated with childhood SSNS. HLA-DRB1*08:02-DQB1*03:02 (Pc=7.01×10^-11; OR, 3.60; 95% CI, 2.46 to 5.29) was identified as the most significant genetic susceptibility factor.Conclusions The most significant association with childhood SSNS was detected in the HLA-DR/DQ region. Further HLA allele/haplotype analyses should enhance our understanding of molecular mechanisms underlying SSNS.

A multiplex RNA quantification method to determine the absolute amounts of mRNA without reverse transcription

We have developed a highly sensitive microarray-based method that determines the absolute amounts of mRNA in a total RNA sample in a multiplex manner without reverse transcription. This direct mRNA measurement promotes high-throughput testing and reduces bias in transcriptome analyses. Furthermore, quantification of the absolute amount of mRNA allows transcriptome analysis without common controls or additional, complicated normalization. The method, called Photo-DEAN, was validated using chemically synthesized RNAs of known quantities and mouse liver total RNA samples. We found that the absolute amounts of mRNA were successfully measured without the cDNA synthesis step, with a sensitivity of 15 zmol achieved in 7 h.

Modification of the Associations of Alcohol Intake With Serum Low-Density Lipoprotein Cholesterol and Triglycerides by ALDH2 and ADH1B Polymorphisms in Japanese Men

Background

Although beneficial associations have been reported between moderate alcohol intake and the serum lipid profile, it is unclear whether polymorphisms in alcohol-metabolizing enzymes can modify these associations. Here, we assessed the effects of ADH1B His48Arg (rs1229984), ALDH2 Glu504Lys (rs671), and their combination on these associations. Furthermore, we examined if the findings for ALDH2 could be replicated.

Methods

We categorized 889 male participants in the Japan Multi-Institutional Collaborative Cohort (J-MICC) Study into two groups based on presence or absence of minor allele(s) or four groups based on genotype combinations. We performed regression analyses of serum lipid concentrations on alcohol intake, with multivariable adjustment. The replication study was conducted among 2,562 men in the Shizuoka part of the J-MICC Study.

Results

The ALDH2 Glu/Lys or Lys/Lys groups showed significant decreases in serum low-density lipoprotein (LDL) cholesterol with increasing alcohol consumption; the coefficient per intake increase of 10 g/day was −2.49 mg/dL (95% confidence interval [CI], −3.85 to −1.13), and a significant interaction with the polymorphism was confirmed (P for interaction = 0.006). This inverse correlation was more evident among the ADH1B His/His + ALDH2 Glu/Lys or Lys/Lys groups (−3.24 mg/dL, 95% CI, −5.03 to −1.45). Serum triglycerides were positively associated with alcohol consumption in the ADH1B His/His group (P for interaction = 0.020). The stronger association between serum LDL cholesterol and alcohol consumption in the ALDH2 Glu/Lys or Lys/Lys groups was replicated.

Conclusions

The ALDH2 Glu504Lys polymorphism can modify the association between alcohol intake and serum LDL cholesterol in Japanese men.

Association of genetic polymorphisms with erythrocyte traits: Verification of SNPs reported in a previous GWAS in a Japanese population

Erythrocyte count and volume are the commonly used hematological indices for anemia that change in various diseases. To date, however, only one study ever exists that addressed erythrocyte trait-associated single nucleotide polymorphisms (SNPs) in a Japanese population. Because that study was performed in patients with various diseases, we confirmed the reported associations in a general population. Participants in the current study were from the Shizuoka component of the Japan Multi-Institutional Collaborative Cohort Study, which included 4971 men and women aged 35 to 69years who were recruited between 2006 and 2007. We analyzed the association of seven selected SNPs with the following erythrocyte traits: red blood cell count, hemoglobin (Hb) and hematocrit (Ht) levels, mean corpuscular volume, mean corpuscular hemoglobin, and mean corpuscular hemoglobin concentration. The erythrocyte traits were regressed on a number of minor alleles of selected SNPs. Then we compared our findings with those from a genome-wide association study performed in a Japanese population. We replicated the association of ABO rs495828, PDGFRA-KIT rs218237, USP49-MED20-BSYL-CCND3 rs3218097, C6orf182-CD164 rs11966072, TERT rs2736100, and TMPRSS6 rs5756504 with erythrocyte traits in our independent Japanese population. In addition, we found a significant interaction between TERT rs2736100 and smoking habit that affected Hb and Ht levels.

Association between kidney function and genetic polymorphisms in atherosclerotic and chronic kidney diseases: A cross-sectional study in Japanese male workers

Background

Several single nucleotide polymorphisms (SNPs) have been implicated in the predisposition to chronic kidney disease (CKD). Atherosclerotic disease is deeply involved in the incidence of CKD; however, whether SNPs related to arteriosclerosis are involved in CKD remains unclear. This study aimed to identify SNPs associated with CKD and to examine whether risk allele accumulation is associated with CKD.

Methods

We conducted a cross-sectional study using data of 4814 male workers to examine the association between estimated glomerular filtration rate (eGFR) and 59 candidate polymorphisms (17 CKD, 42 atherosclerotic diseases). We defined the genetic risk score (GRS) as the total number of risk alleles that showed a significant association in this analysis and examined the relationship with CKD (eGFR < 60 ml/min/1.73m²). Multivariate logistic regression, discrimination by area under the receiver operating characteristic curve, integrated discrimination improvement (IDI), and category-free net reclassification improvement (cNRI) were evaluated.

Results

In total, 432 participants were categorized as having CKD. We found eight candidate SNPs with P value < 0.05 (CX3CR1 rs3732379, SHROOM3 rs17319721, MTP rs1800591, PIP5K1B rs4744712, APOA5 rs662799, BRAP rs3782886, SPATA5L1 rs2467853, and MCP1 rs1024611) in the multivariate linear regression adjusted for age, body mass index, systolic blood pressure, and fasting blood glucose. Among these eight SNPs, BRAP rs3782886 and SPATA5L1 rs2467853 were significantly associated with eGFR (false discovery rate < 0.05). GRS was significantly associated with CKD (odds ratio, 1.17; 95% confidence interval, 1.09–1.26). C-statisics improved from 0.775 to 0.780 but showed no statistical significance. However, adding GRS significantly improved IDI and cNRI (0.0057, P = 0.0028, and 0.212, P < 0.001, respectively).

Conclusions

After adjustment for clinical factors, kidney function was associated with BRAP rs3782886 and SPATA5L1 rs2467853 and the GRS for CKD that we developed was associated CKD.

Interferon regulatory factor 6 variants affect nasolabial morphology in East Asian populations

OBJECTIVE

The interferon regulatory factor 6 gene (IRF6) is one of the most conspicuous genes among a large number of candidate risk genes for non-syndromic cleft lip with or without cleft palate, which is considered to be a multifactorial defect. Variants of IRF6 are also suggested to affect normal craniofacial variations, especially in the area of the nose and the upper lip. In the present study, we used lateral cephalograms to establish the relationship between IRF6 and sagittal nasolabial morphology in healthy East Asian subjects.

DESIGN

Genomic DNA was extracted from 215 Japanese and 226 Korean individuals, and genotyped for five IRF6 single nucleotide polymorphisms (SNPs): rs17389541, rs642961, rs2013162, rs2235371, and rs7802. These SNPs were tested by multiple regression analyses for their association with craniofacial measurements obtained from lateral cephalometrics.

RESULTS

We detected a significant association between the derived variants, rs2013162 and rs2235371 and the distances between a facial bone plane indicated by distance from Nasion and Point A (NA plane) to soft tissue landmarks; the Subalare (NA-Sbal) and the Subnasale (NA-Sn) in the sagittal plane.

CONCLUSION

Our results indicate that IRF6 variants play an important role in the normal range of variation in nasolabial soft-tissue morphology.

Contribution of FGFR1 Variants to Craniofacial Variations in East Asians

FGFR1 plays an important role in the development of the nervous system as well as the regulation of the skeletal development and bone homeostasis. Mutations in FGFR1 genes affect skull development, specifically suture and synchondrosis, resulting in craniosynostosis and facial abnormalities. We examined subjects with normal skull morphology for genetic polymorphisms that might be associated with normal craniofacial variations. Genomic DNA was obtained from 216 Japanese and 227 Korean subjects. Four FGFR1 SNPs, namely, rs881301, rs6996321, rs4647905, and rs13317, were genotyped. These SNPs were tested for association with craniofacial measurements obtained from lateral and posteroanterior cephalometries, in which principle component analysis was performed to compress the data of the craniofacial measurements. We observed that SNPs rs13317 and rs6996321 were correlated with the overall head size and midfacial development, indicating that FGFR1 SNPs played crucial roles in the normal variation of human craniofacial morphology. Subjects with the derived alleles of SNPs rs13317 and rs6996321 had a small face and a facial pattern associated with a retruded midface and relatively wide-set eyes. These facial features were similar to but were milder than those of individuals with Pfeiffer syndrome, which is caused by a dysfunctional mutation in FGFR1.

Genome-Wide Association Study Identifies ZNF354C Variants Associated with Depression from Interferon-Based Therapy for Chronic Hepatitis C

The therapeutic use of interferon (IFN) is known to cause depression that frequently interrupts treatment. To identify genetic variants associated with IFN-induced depression, we conducted a genome-wide association study (GWAS) of 224 Japanese chronic hepatitis C patients receiving IFN-based therapy in a multicenter prospective study and stratified them into two groups according to the Beck Depression Inventory, Second Edition (BDI-II) score. In the GWAS stage, we selected 42 candidate single nucleotide polymorphisms (SNPs) to perform replication analysis in an independent set of 160 subjects. The SNP rs1863918 in strong linkage disequilibrium with SNPs located around the Zinc finger 354C (ZNF354C) gene on chromosome 5 showed a significant association when the results of GWAS and replication were combined (odds ratio = 2.55, P = 7.89×10-8 in the allele frequency model), suggesting that the rs1863918 T allele was associated with IFN-induced depression. Furthermore, logistic regression analysis showed that rs1863918 T allele, a history of depression, and younger age were independent predictive factors for IFN-induced depression. Interestingly, western blotting and immunofluorescence showed that ZNF354C was highly expressed in the hippocampus in mice, a region implicated in the pathology of psychiatric symptoms. In conclusion, we identified rs1863918 as significantly associated with IFN-induced depression, and revealed that the candidate gene ZNF354C is highly expressed in the hippocampus of mice. Our data might be useful for elucidating the pathogenic mechanisms of depression induced by drugs including IFN.

Cold medicine-related Stevens-Johnson syndrome/toxic epidermal necrolysis with severe ocular complications-phenotypes and genetic predispositions

Stevens–Johnson syndrome (SJS) is an acute inflammatory vesiculobullous reaction of the skin and mucosa, such as the ocular surface, oral cavity, and genitals. In patients with extensive skin detachment and a poor prognosis, the condition is called toxic epidermal necrolysis (TEN). Severe ocular complications (SOCs) appear in some–but not all–SJS/TEN patients who are diagnosed by dermatologists, and cold medicines including multi-ingredient cold medications and nonsteroidal anti-inflammatory drugs are the main causative drugs particularly for SJS/TEN with SOCs and all SJS and TEN. In this review, we focus on the genetic predisposition of cold medicine-related SJS/TEN (CM-SJS/TEN) with SOCs. CM-SJS/TEN with SOCs was strongly associated with HLA-A*02:06 and significantly associated with HLA-B*44:03 in Japanese individuals, significantly associated with HLA-B*44:03 in Indian and Brazilian individuals, and associated with HLA-A*02:06 in Korean individuals. In the first genome-wide association study (GWAS), we found an association between the prostaglandin E receptor 3 (PTGER3) gene and SJS/TEN with SOCs. In this study, we focused on CM-SJS/TEN with SOCs and found that the association of CM-SJS/TEN with SOCs became stronger than all SJS/TEN with SOCs. In the second GWAS, we found an association between the IKZF1 gene and CM-SJS/TEN with SOCs not only in Japanese, but also in Korean and Indian populations. Moreover, we found that TSHZ2 gene single nucleotide polymorphisms (SNPs) also showed especially low p values in the Japanese population; however, this association was not found in the Korean population. Furthermore, we investigated the interaction between susceptibility genes, and found multiplicative interactions of HLA-A*02:06 and TLR3 SNPs and additive interactions of HLA-A*02:06 and PTGER3 SNPs.

Association of interactions between dietary salt consumption and hypertension-susceptibility genetic polymorphisms with blood pressure among Japanese male workers

BACKGROUND

Blood pressure is influenced by hereditary factors and dietary habits. The objective of this study was to examine the effect of dietary salt consumption and single-nucleotide polymorphisms (SNPs) on blood pressure (BP).

METHODS

This was a cross-sectional analysis of 2728 male participants who participated in a health examination in 2009. Average dietary salt consumption was estimated using electronically collected meal purchase data from cafeteria. A multivariate analysis, adjusting for clinically relevant factors, was conducted to examine whether the effect on BP of salt consumption, SNPs, and interaction between salt consumption and each SNP. This study examined the SNPs AGT rs699 (Met235Thr), ADD1 rs4961 (Gly460Trp), NPPA rs5063 (Val32Met), GPX1 rs1050450 (Pro198Leu), and AGTR1 rs5186 (A1166C) in relation to hypertension and salt sensitivity.

RESULTS

BP was not significantly associated with SNPs or salt consumption. The interaction between salt consumption and SNPs with systolic BP showed a significant association in NPPA rs5063 (Val32Met) (P = 0.023) and a marginal trend toward significance in rs4961 and rs1050450 (P = 0.060 and 0.067, respectively).

CONCLUSION

The effect of salt consumption on BP differed by genotype. Dietary salt consumption and genetic variation can predict a high risk of hypertension.

A Study of Single Nucleotide Polymorphisms of the SLC19A1/RFC1 Gene in Subjects with Autism Spectrum Disorder

Autism Spectrum Disorder (ASD) is a group of neurodevelopmental disorders with complex genetic etiology. Recent studies have indicated that children with ASD may have altered folate or methionine metabolism, suggesting that the folate-methionine cycle may play a key role in the etiology of ASD. SLC19A1, also referred to as reduced folate carrier 1 (RFC1), is a member of the solute carrier group of transporters and is one of the key enzymes in the folate metabolism pathway. Findings from multiple genomic screens suggest the presence of an autism susceptibility locus on chromosome 21q22.3, which includes SLC19A1. Therefore, we performed a case-control study in a Japanese population. In this study, DNA samples obtained from 147 ASD patients at the Kanazawa University Hospital in Japan and 150 unrelated healthy Japanese volunteers were examined by the sequence-specific primer-polymerase chain reaction method pooled with fluorescence correlation spectroscopy. p < 0.05 was considered to represent a statistically significant outcome. Of 13 single nucleotide polymorphisms (SNPs) examined, a significant p-value was obtained for AA genotype of one SNP (rs1023159, OR = 0.39, 95% CI = 0.16-0.91, p = 0.0394; Fisher's exact test). Despite some conflicting results, our findings supported a role for the polymorphism rs1023159 of the SLC19A1 gene, alone or in combination, as a risk factor for ASD. However, the findings were not consistent after multiple testing corrections. In conclusion, although our results supported a role of the SLC19A1 gene in the etiology of ASD, it was not a significant risk factor for the ASD samples analyzed in this study.

The impact of PNPLA3 and JAZF1 on hepatocellular carcinoma in non-viral hepatitis patients with type 2 diabetes mellitus

BACKGROUND

Type 2 diabetes mellitus (T2DM) is an established independent risk factor for hepatocellular carcinoma (HCC). T2DM is associated with non-alcoholic steatohepatitis (NASH), which is a major cause of non-HBV and non-HCV-related HCC; nevertheless, it has been difficult to identify those patients with T2DM who have a high risk of developing HCC. The aim of this study was to identify genetic determinants that predispose T2DM patients to HCC by genotyping T2DM susceptibility loci and PNPLA3.

METHODS

We recruited 389 patients with T2DM who satisfied the following three criteria: negative for HBs-Ag and anti-HCV Ab, alcohol intake <60 g/day, and history of T2DM >10 years. These patients were divided into two groups: T2DM patients with HCC (DM-HCC, n = 59) or those without HCC (DM-non-HCC, n = 330). We genotyped 51 single-nucleotide polymorphisms (SNPs) previously reported as T2DM or NASH susceptibility loci (PNPLA3) compared between the DM-HCC and DM-non-HCC groups with regard to allele frequencies at each SNP.

RESULTS

The SNP rs738409 located in PNPLA3 was the greatest risk factor associated with HCC. The frequency of the PNPLA3 G allele was significantly higher among DM-HCC individuals than DM-non-HCC individuals (OR 2.53, p = 1.05 × 10(-5)). Among individuals homozygous for the PNPLA3 G allele (n = 115), the frequency of the JAZF1 rs864745 G allele was significantly higher among DM-HCC individuals than DM-non-HCC individuals (OR 3.44, p = 0.0002).

CONCLUSIONS

PNPLA3 and JAZF1 were associated with non-HBV and non-HCV-related HCC development among Japanese patients with T2DM.

A polymorphism in CCR1/CCR3 is associated with narcolepsy

Etiology of narcolepsy-cataplexy involves multiple genetic and environmental factors. While the human leukocyte antigen (HLA)-DRB1*15:01-DQB1*06:02 haplotype is strongly associated with narcolepsy, it is not sufficient for disease development. To identify additional, non-HLA susceptibility genes, we conducted a genome-wide association study (GWAS) using Japanese samples. An initial sample set comprising 409 cases and 1562 controls was used for the GWAS of 525,196 single nucleotide polymorphisms (SNPs) located outside the HLA region. An independent sample set comprising 240 cases and 869 controls was then genotyped at 37 SNPs identified in the GWAS. We found that narcolepsy was associated with a SNP in the promoter region of chemokine (C-C motif) receptor 1 (CCR1) (rs3181077, P=1.6×10(-5), odds ratio [OR]=1.86). This rs3181077 association was replicated with the independent sample set (P=0.032, OR=1.36). We measured mRNA levels of candidate genes in peripheral blood samples of 38 cases and 37 controls. CCR1 and CCR3 mRNA levels were significantly lower in patients than in healthy controls, and CCR1 mRNA levels were associated with rs3181077 genotypes. In vitro chemotaxis assays were also performed to measure monocyte migration. We observed that monocytes from carriers of the rs3181077 risk allele had lower migration indices with a CCR1 ligand. CCR1 and CCR3 are newly discovered susceptibility genes for narcolepsy. These results highlight the potential role of CCR genes in narcolepsy and support the hypothesis that patients with narcolepsy have impaired immune function.

Genetic structure and relationships of 16 Asian and European cattle populations using DigiTag2 assay

In this study, we genotyped 117 autosomal single nucleotide polymorphisms using a DigiTag2 assay to assess the genetic diversity, structure and relationships of 16 Eurasian cattle populations, including nine cattle breeds and seven native cattle. Phylogenetic and principal component analyses showed that Bos taurus and Bos indicus populations were clearly distinguished, whereas Japanese Shorthorn and Japanese Polled clustered with European populations. Furthermore, STRUCTURE analysis demonstrated the distinct separation between Bos taurus and Bos indicus (K=2), and between European and Asian populations (K=3). In addition, Japanese Holstein exhibited an admixture pattern with Asian and European cattle (K=3‐5). Mongolian (K=13‐16) and Japanese Black (K=14‐16) populations exhibited admixture patterns with different ancestries. Bos indicus populations exhibited a uniform genetic structure at K=2‐11, thereby suggesting that there are close genetic relationships among Bos indicus populations. However, the Bhutan and Bangladesh populations formed a cluster distinct from the other Bos indicus populations at K=12‐16. In conclusion, our study could sufficiently explain the genetic construction of Asian cattle populations, including: (i) the close genetic relationships among Bos indicus populations; (ii) the genetic influences of European breeds on Japanese breeds; (iii) the genetic admixture in Japanese Holstein, Mongolian and Japanese Black cattle; and (iv) the genetic subpopulations in Southeast Asia.

HLA-A*02:06 and PTGER3 polymorphism exert additive effects in cold medicine-related Stevens-Johnson syndrome with severe ocular complications

We previously reported that PTGER3 (prostaglandin E receptor 3 (subtype EP3)) single-nucleotide polymorphisms (SNPs) were associated with Stevens–Johnson syndrome (SJS)/toxic epidermal necrolysis (TEN) with severe ocular complications (SOC). We also documented that approximately 80% of our SJS/TEN patients had taken cold medicines within several days before disease onset, and we thus designated them cold medicine-related SJS/TEN (CM-SJS/TEN) patients. Moreover, we reported that HLA-A*02:06 with TLR3 polymorphisms exerted more than additive effects in SJS/TEN with SOC. In this study, we focused on CM-SJS/TEN with SOC and analyzed the association with PTGER3 SNPs and an interactive effect between PTGER3 SNPs and HLA-A*02:06 in not only the Japanese but also the Korean population. In the Japanese population, PTGER3 SNP rs1327464 was most significantly associated with CM-SJS/TEN with SOC (G versus A; odds ratio (OR)=0.232, P=7.92×10⁻¹⁰), and we found an interaction with additive effects between HLA-A*02:06 and the high-risk genotypes PTGER3 rs1327464 GA or AA (OR=10.8, P=2.56×10⁻⁷). We also found a significant association between Korean CM-SJS/TEN with SOC and PTGER3 SNP rs1327464 (GG versus GA+AA, OR=0.246, P=0.00101), and we detected an additive effect between HLA-A*02:06 and the high-risk genotypes PTGER3 rs1327464 GA or AA (OR=14.2, P=5.58×10⁻⁶).

Association Study between the CD157/BST1 Gene and Autism Spectrum Disorders in a Japanese Population

CD157, also referred to as bone marrow stromal cell antigen-1 (BST-1), is a glycosylphosphatidylinositol-anchored molecule that promotes pre-B-cell growth. Previous studies have reported associations between single-nucleotide polymorphisms (SNPs) of the CD157/BST1 gene with Parkinson's disease. In an attempt to determine whether SNPs or haplotypes in the CD157/BST1 are associated with other brain disorders, we performed a case-control study including 147 autism spectrum disorder (ASD) patients at Kanazawa University Hospital in Japan and 150 unselected Japanese volunteers by the sequence-specific primer-polymerase chain reaction method combined with fluorescence correlation spectroscopy. Of 93 SNPs examined, two SNPs showed significantly higher allele frequencies in cases with ASDs than in unaffected controls (rs4301112, OR = 6.4, 95% CI = 1.9 to 22, p = 0.0007; and rs28532698, OR = 6.2, 95% CI = 1.8 to 21, p = 0.0012; Fisher's exact test; p < 0.002 was considered significant after multiple testing correction). In addition, CT genotype in rs10001565 was more frequently observed in the ASD group than in the control group (OR = 15, 95% CI = 2.0 to 117, p = 0.0007; Fisher's exact test). The present data indicate that genetic variation of the CD157/BST1 gene might confer susceptibility to ASDs.

Genome-wide association study identifies a PSMD3 variant associated with neutropenia in interferon-based therapy for chronic hepatitis C

Cytopenia during interferon-based (IFN-based) therapy for chronic hepatitis C (CHC) often necessitates reduction of doses of drugs and premature withdrawal from therapy resulting in poor response to treatment. To identify genetic variants associated with IFN-induced neutropenia, we conducted a genome-wide association study (GWAS) in 416 Japanese CHC patients receiving IFN-based therapy. Based on the results, we selected 192 candidate single nucleotide polymorphisms (SNPs) to carry out a replication analysis in an independent set of 404 subjects. The SNP rs2305482, located in the intron region of the PSMD3 gene on chromosome 17, showed a strong association when the results of GWAS and the replication stage were combined (OR = 2.18, P = 3.05 × 10(-7) in the allele frequency model). Logistic regression analysis showed that rs2305482 CC and neutrophil count at baseline were independent predictive factors for IFN-induced neutropenia (OR = 2.497, P = 0.0072 and OR = 0.998, P < 0.0001, respectively). Furthermore, rs2305482 genotype was associated with the doses of pegylated interferon (PEG-IFN) that could be tolerated in hepatitis C virus genotype 1-infected patients treated with PEG-IFN plus ribavirin, but not with treatment efficacy. Our results suggest that genetic testing for this variant might be useful for establishing personalized drug dosing in order to minimize drug-induced adverse events.

IKZF1, a new susceptibility gene for cold medicine-related Stevens-Johnson syndrome/toxic epidermal necrolysis with severe mucosal involvement

BACKGROUND

Stevens-Johnson syndrome (SJS) and its severe form, toxic epidermal necrolysis (TEN), are acute inflammatory vesiculobullous reactions of the skin and mucous membranes, including the ocular surface, oral cavity, and genitals. These reactions are very rare but are often associated with inciting drugs, infectious agents, or both.

OBJECTIVE

We sought to identify susceptibility loci for cold medicine-related SJS/TEN (CM-SJS/TEN) with severe mucosal involvement (SMI).

METHODS

A genome-wide association study was performed in 808 Japanese subjects (117 patients with CM-SJS/TEN with SMI and 691 healthy control subjects), and subsequent replication studies were performed in 204 other Japanese subjects (16 cases and 188 control subjects), 117 Korean subjects (27 cases and 90 control subjects), 76 Indian subjects (20 cases and 56 control subjects), and 174 Brazilian subjects (39 cases and 135 control subjects).

RESULTS

In addition to the most significant susceptibility region, HLA-A, we identified IKZF1, which encodes Ikaros, as a novel susceptibility gene (meta-analysis, rs4917014 [G vs. T]; odds ratio, 0.5; P = 8.5 × 10(-11)). Furthermore, quantitative ratios of the IKZF1 alternative splicing isoforms Ik1 and Ik2 were significantly associated with rs4917014 genotypes.

CONCLUSION

We identified IKZF1 as a susceptibility gene for CM-SJS/TEN with SMI not only in Japanese subjects but also in Korean and Indian subjects and showed that the Ik2/Ik1 ratio might be influenced by IKZF1 single nucleotide polymorphisms, which were significantly associated with susceptibility to CM-SJS/TEN with SMI.

Identification of an interaction between VWF rs7965413 and platelet count as a novel risk marker for metabolic syndrome: an extensive search of candidate polymorphisms in a case-control study

Although many single nucleotide polymorphisms (SNPs) have been identified to be associated with metabolic syndrome (MetS), there was only a slight improvement in the ability to predict future MetS by the simply addition of SNPs to clinical risk markers. To improve the ability to predict future MetS, combinational effects, such as SNP—SNP interaction, SNP—environment interaction, and SNP—clinical parameter (SNP × CP) interaction should be also considered. We performed a case-control study to explore novel SNP × CP interactions as risk markers for MetS based on health check-up data of Japanese male employees. We selected 99 SNPs that were previously reported to be associated with MetS and components of MetS; subsequently, we genotyped these SNPs from 360 cases and 1983 control subjects. First, we performed logistic regression analyses to assess the association of each SNP with MetS. Of these SNPs, five SNPs were significantly associated with MetS (P < 0.05): LRP2 rs2544390, rs1800592 between UCP1 and TBC1D9, APOA5 rs662799, VWF rs7965413, and rs1411766 between MYO16 and IRS2. Furthermore, we performed multiple logistic regression analyses, including an SNP term, a CP term, and an SNP × CP interaction term for each CP and SNP that was significantly associated with MetS. We identified a novel SNP × CP interaction between rs7965413 and platelet count that was significantly associated with MetS [SNP term: odds ratio (OR) = 0.78, P = 0.004; SNP × CP interaction term: OR = 1.33, P = 0.001]. This association of the SNP × CP interaction with MetS remained nominally significant in multiple logistic regression analysis after adjustment for either the number of MetS components or MetS components excluding obesity. Our results reveal new insight into platelet count as a risk marker for MetS.

A genome-wide association study for diabetic retinopathy in a Japanese population: potential association with a long intergenic non-coding RNA

Elucidation of the genetic susceptibility factors for diabetic retinopathy (DR) is important to gain insight into the pathogenesis of DR, and may help to define genetic risk factors for this condition. In the present study, we conducted a three-stage genome-wide association study (GWAS) to identify DR susceptibility loci in Japanese patients, which comprised a total of 837 type 2 diabetes patients with DR (cases) and 1,149 without DR (controls). From the stage 1 genome-wide scan of 446 subjects (205 cases and 241 controls) on 614,216 SNPs, 249 SNPs were selected for the stage 2 replication in 623 subjects (335 cases and 288 controls). Eight SNPs were further followed up in a stage 3 study of 297 cases and 620 controls. The top signal from the present association analysis was rs9362054 in an intron of RP1-90L14.1 showing borderline genome-wide significance (P_met = 1.4×10⁻⁷, meta-analysis of stage 1 and stage 2, allele model). RP1-90L14.1 is a long intergenic non-coding RNA (lincRNA) adjacent to KIAA1009/QN1/CEP162 gene; CEP162 plays a critical role in ciliary transition zone formation before ciliogenesis. The present study raises the possibility that the dysregulation of ciliary-associated genes plays a role in susceptibility to DR.

Serum interferon-gamma-inducible protein-10 concentrations and IL28B genotype associated with responses to pegylated interferon plus ribavirin with and without telaprevir for chronic hepatitis C

AIM

Several studies have shown that high pretreatment concentrations of serum interferon-γ-inducible protein-10 (IP-10) are correlated with non-response to pegylated interferon (PEG-IFN) plus ribavirin (RBV) for chronic hepatitis C (CHC). However, there are few reports on their effect on the Asian population.

METHODS

We enrolled 104 Japanese genotype 1 CHC individuals treated with PEG-IFN/RBV and 45 with PEG-IFN/RBV/telaprevir, and evaluated the impact of pretreatment serum IP-10 concentrations on their virological responses.

RESULTS

The pretreatment serum IP-10 concentrations were not correlated with IL28B genotype. The receiver-operator curve analysis determined the cut-off value of IP-10 for predicting a sustained virological response (SVR) as 300 pg/mL. In multivariate analysis, the IL28B favorable genotype and IP-10 concentration of less than 300 pg/mL were independent factors for predicting SVR. In a subgroup of patients with the IL28B favorable genotype, the SVR rate was higher in the patients with IP-10 of less than 300 than in those with 300 pg/mL or more, whereas no patient with the IL28B unfavorable genotype and IP-10 of 300 pg/mL or more achieved SVR. Among the patients treated with PEG-IFN/RBV/telaprevir, low pretreatment concentrations of serum IP-10 were associated with a very rapid virological response, defined as undetectable HCV RNA at week 2 after the start of therapy.

CONCLUSION

Pretreatment serum IP-10 concentrations are associated with treatment efficacy in PEG-IFN/RBV and with early viral kinetics of hepatitis C virus in PEG-IFN/RBV/telaprevir therapy.

Association of functional polymorphisms in interferon regulatory factor 2 (IRF2) with susceptibility to systemic lupus erythematosus: a case-control association study

Interferon regulatory factor 2 (IRF2) negatively regulates type I interferon (IFN) responses, while it plays a role in induction of Th1 differentiation. Previous linkage and association studies in European-American populations suggested genetic role of IRF2 in systemic lupus erythematosus (SLE); however, this observation has not yet been confirmed. No studies have been reported in the Asian populations. Here we investigated whether IRF2 polymorphisms contribute to susceptibility to SLE in a Japanese population. Association study of 46 IRF2 tag single nucleotide polymorphisms (SNPs) detected association of an intronic SNP, rs13146124, with SLE. When the association was analyzed in 834 Japanese patients with SLE and 817 healthy controls, rs13146124 T was significantly increased in SLE compared with healthy controls (dominant model, P = 5.4×10⁻⁴, Bonferroni-corrected P [Pc] = 0.026, odds ratio [OR] 1.48, 95% confidence interval [CI] 1.18–1.85). To find causal SNPs, resequencing was performed by next-generation sequencing. Twelve polymorphisms in linkage disequilibrium with rs13146124 (r²: 0.30–1.00) were identified, among which significant association was observed for rs66801661 (allele model, P = 7.7×10⁻⁴, Pc = 0.037, OR 1.53, 95%CI 1.19–1.96) and rs62339994 (dominant model, P = 9.0×10⁻⁴, Pc = 0.043, OR 1.46, 95%CI 1.17–1.82). The haplotype carrying both of the risk alleles (rs66801661A–rs62339994A) was significantly increased in SLE (P = 9.9×10⁻⁴), while the haplotype constituted by both of the non-risk alleles (rs66801661G–rs62339994G) was decreased (P = 0.0020). A reporter assay was carried out to examine the effect of the IRF2 haplotypes on the transcriptional activity, and association of the IRF2 risk haplotype with higher transcriptional activity was detected in Jurkat T cells under IFNγ stimulation (Tukey's test, P = 1.2×10⁻⁴). In conclusion, our observations supported the association of IRF2 with susceptibility to SLE, and the risk haplotype was suggested to be associated with transcriptional activation of IRF2.

A PCR method for VKORC1 G-1639A and CYP2C9 A1075C genotyping useful to warfarin therapy among Japanese

Warfarin is widely prescribed for patients with the risk of thromboembolism around the world. The inter-individual and inter-racial differences in appropriate dosage depend highly on age, body weight, and genetic factors. A lot of studies including genome-wide association studies revealed that vitamin K epoxide reductase complex, subunit 1 (VKORC1) G-1639A and Cytochrome P450 (CYP) 2C9 A1075C are the most strong genetic factors for determining warfarin effects in Asians and Africans. Since we developed a quick and inexpensive genotyping method, polymerase chain reaction with confronting two-pair primers (PCR-CTPP), the method was applied for these genotypes to examine the possibility to clinical use. Subjects were 436 examinees (117 males and 319 females, aged 32 to 85 years) who attended a health checkup program in Japan. The PCR-CTPP for VKORC1 G-1639A and CYP2C9 A1075C was conducted for the subjects, as well as the samples genotyped by DigiTag2 method. The allele frequencies of VKORC1 G-1639A were 0.085 for G and 0.915 for A, and those of CYP2C9 A1075C were 0.979 for A and 0.021 for C, being in Hardy-Weinberg equilibrium (p = 0.658 and p = 0.514, respectively). These frequencies were similar to those reported in the HapMap project. Genotyping for both SNPs by PCR-CTPP was replicated by DigiTag2 method. Our results indicated that the PCR-CTPP could be one of the alternative methods for genotyping VKORC1 G-1639A and CYP2C9 A1075C for Asians and Africans with similar allele frequencies to Japanese.

Tag/hybridization-based sensitive detection of polymerase chain reaction products

The polymerase chain reaction (PCR) is an important technology to amplify a single copy or a few copies of DNA segment in genomic DNAs, visualizing the segment as DNA fragment. Thus, PCR is frequently used in various examinations such as detection of bacteria and fungi in the food industry. Here, we report a simple and sensitive method for detection of PCR products using single-strand tag sequence and hybridization of the tag sequence to the complementary tag sequence immobilized on solid material (STH). The detection sensitivity was found to be at least 50 times higher than electrophoresis/ethidium bromide (EtBr) visualization for approximately a 500-bp fragment and higher than the ordinary hybridization, that is, hybridization of denatured PCR product to probe sequence immobilized on solid material.

ITPA genetic variants influence efficacy of PEG-IFN/RBV therapy in older patients infected with HCV genotype 1 and favourable IL28B type

Inosine triphosphatase (ITPA) genetic variants are strongly associated with ribavirin (RBV)-induced anaemia during pegylated interferon (PEG-IFN) plus RBV therapy. However, the treatment efficacy of ITPA genetic variants has not been fully explored. We enrolled 309 individuals infected with hepatitis C virus genotype 1, who were treated with PEG-IFN plus RBV for 48 weeks. The ITPA SNP: rs1127354 and IL28B SNP: rs8099917 were genotyped. We examined the risk factors for severe anaemia up to week 12 after the start of treatment and treatment efficacy. The incidence of severe anaemia, ≥ 3 g/dL reduction or <10 g/dL of haemoglobin (Hb) up to week 12, was more frequent in patients with CC at rs1127354 [65% (145/224), 33% (73/224)] than in those with CA/AA [25% (21/85), 6% (8/85)] (P < 0.0001). ITPA genotype, pretreatment Hb level and age were independent predictive factors for severe anaemia: Hb < 10 g/dL. In IL28B favourable type, the sustained virologic response rate was higher in ≥ 60-year-old patients with CA/AA than in those with CC [71% (22/31) vs 40% (26/65), P = 0.005], although there was no significant difference in treatment efficacy according to ITPA genetic variants in the <60-year-old patients. The proportion of patients administered ≥ 80% of the dosage of RBV was significantly higher in the patients with CA/AA than in those with CC (P = 0.025), resulting in a lower relapse rate. In conclusion, ITPA genetic variants were associated with severe RBV-induced anaemia and could influence the efficacy of PEG-IFN plus RBV treatment among elderly patients with IL28B favourable type.

Novel DNA chip based on a modified DigiTag2 assay for high-throughput species identification and genotyping of Mycobacterium tuberculosis complex isolates

A multipurpose high-throughput genotyping tool for the assessment of recent epidemiological data and evolutional pattern in Mycobacterium tuberculosis complex (MTBC) clinical isolates was developed in this study. To facilitate processing, 51 highly informative single nucleotide polymorphisms (SNPs) were selected for discriminating the clinically most relevant MTBC species and genotyping M. tuberculosis into its principle genetic groups (PGGs) and SNP cluster groups (SCGs). Because of the high flexibility of the DigiTag2 assay, the identical protocol and DNA array containing the identical set of probes were applied to the highly GC-rich mycobacterial genome. The specific primers with multiplex amplification and hybridization conditions based on the DigiTag2 principle were optimized and evaluated with 14 MTBC reference strains, 4 nontuberculous mycobacteria (NTM) isolates, and 322 characterized M. tuberculosis clinical isolates. The DNA chip that was developed revealed a 99.85% call rate, a 100% conversion rate, and 99.75% reproducibility. For the accuracy rate, 98.94% of positive calls were consistent with previous molecular characterizations. Our cost-effective technology was capable of simultaneously identifying the MTBC species and the genotypes of 96 M. tuberculosis clinical isolates within 6 h using only simple instruments, such as a thermal cycler, a hybridization oven, and a DNA chip scanner, and less technician skill was required than for other techniques. We demonstrate this approach's potential as a simple, flexible, and rapid tool for providing clearer information regarding circulating MTBC isolates.

DNA polymorphism in the FKBP5 gene affects impulsivity in intertemporal choice

INTRODUCTION

Impulsivity in intertemporal choice has been operationalized as "delay discounting", referring to the preference for a sooner, smaller reward. FK506 binding protein 5 (FKBP5) is a co-chaperone of the glucocorticoid receptor (GR). FKBP5 overexpression causes GR resistance, resulting in increased plasma cortisol levels. High cortisol levels are associated with low impulsivity in intertemporal choice. The aim of this study was to explore the effect of single nucleotide polymorphisms (SNPs) in FKBP5 on delay discounting.

METHODS

The participants consisted of 91 healthy Japanese people (66 males and 25 females with a mean age of 40.9 ± 6.9 years). Each participant completed Kirby's monetary choice questionnaire (MCQ) and donated a whole blood sample. Five SNPs in FKBP5 were genotyped using the DigiTag2. SNP linear regression analyses with 100,000 permutations were conducted for the hyperbolic time-discount rate (k).

RESULTS

Two SNPs were excluded from analysis because of their low minor allelic frequencies. The SNP rs1360780 showed a significant association; participants with more minor alleles (T) were less impulsive in intertemporal choice for delayed gain (multiplicity-corrected P = 0.047).

DISCUSSION

The significant SNP rs1360780 is located in the region adjacent to the hormone response element (HRE)-binding sequence where transcription factors bind and alter the transcription of FKBP5. A minor allele (T) of rs1360780, which causes FKBP5 overexpression, may reduce impulsivity in intertemporal choice (i.e. delay discounting) via GR resistance and the subsequent high cortisol levels. This is the first study to demonstrate an association between FKBP5 and impulsivity in intertemporal choice.

A functional SNP upstream of the beta-2 adrenergic receptor gene (ADRB2) is associated with obesity in Oceanic populations

OBJECTIVE:

Obesity is a growing health concern in the Oceanic populations. To investigate the genetic factors associated with adult obesity in the Oceanic populations, the association of single nucleotide polymorphisms (SNPs) of the beta-2 adrenergic receptor (ADRB2) gene with obesity was examined in 694 adults living in Tonga and Solomon Islands.

RESULTS:

A screening for variation in 16 Oceanic subjects detected 17 SNPs in the entire region of ADRB2, of which nine SNPs including two non-synonymous ones, rs1042713 (Arg16Gly) and rs1042714 (Gln27Glu), were further genotyped for all subjects. The rs34623097-A allele, at a SNP located upstream of ADRB2, showed the strongest association with risk for obesity in a logistic regression analysis adjusted for age, sex, and population (P=5.6 × 10⁻⁴, odds ratio [OR]=2.5, 95% confidence interval [CI]=1.5–4.2). The 27Glu was also significantly associated with obesity in the single-point association analysis (P=0.013, OR=2.0, 95%CI=1.2–3.4); however, this association was no longer significant after adjustment for rs34623097 since these SNPs were in linkage disequilibrium with each other. A copy of the obesity-risk allele, rs34623097-A, led to a 1.6 kg/m² increase in body mass index (BMI; defined as weight in kilograms divided by height in meters squared) (P=0.0019). A luciferase reporter assay indicated that rs34623097-A reduced the transcriptional activity of the luciferase reporter gene by approximately 10% compared with rs34623097-G. An electrophoretic mobility shift assay demonstrated that rs34623097 modulated the binding affinity with nuclear factors. An evolutionary analysis implies that a G>A mutation at rs34623097 occurred in the Neandertal genome and then the rs34623097-A allele flowed into the ancestors of present-day humans.

CONCLUSION:

The present results suggest that rs34623097-A, which would lead to lower expression of ADRB2, contributes to the onset of obesity in the Oceanic populations.

Changes in sequences of core region, interferon sensitivity-determining region and interferon and ribavirin resistance-determining region of hepatitis C virus genotype 1 during interferon-alpha and ribavirin therapy, and efficacy of retreatment

AIM

  Some regions associated with sensitivity to interferon-α and ribavirin have been identified in the hepatitis C virus (HCV) genome, including amino acid 70 in the core region (core a.a. 70), a.a. 2209-2248 (interferon sensitivity-determining region, ISDR) and a.a. 2334-2379 (interferon and ribavirin resistance-determining region, IRRDR).

METHODS

  We examined changes in the sequences of these regions in 25 patients with chronic HCV genotype 1 infection who had not had sustained virological response (SVR) to interferon-α and ribavirin for 24-48 weeks and subsequently received retreatment for 48-72 weeks.

RESULTS

  At baseline, the core a.a. 70 was mutant (resistant) type in seven patients. At the start of retreatment, the core a.a. 70 had changed from sensitive to resistant type in 2 patients, and SVR was not achieved by retreatment. The ISDR variations were resistant type (0-1 mutations) in 17 patients at baseline. After 2 weeks of treatment, amino acid change was found in two patients; in one, the substitutions returned to baseline status after treatment, and in the other, the substitution persisted. At the start of retreatment, ISDR sequences had changed from resistant to sensitive type in two patients and SVR was achieved and from sensitive to resistant type in three patients and SVR was not achieved. The IRRDR variations were resistant type (<6 mutations) in 19 patients at baseline and at the start of retreatment.

CONCLUSION

  Sequences of the core region and ISDR sometimes change during anti-HCV therapy, potentially affecting the outcomes of retreatment.

Common variants on chromosome 9p21 are associated with normal tension glaucoma

Although intraocular pressure (IOP) is the most definitive cause of glaucoma, a subtype of open angle glaucoma (OAG) termed normal tension glaucoma (NTG), which occurs in spite of normal IOP, accounts for a large part of glaucoma cases, especially in Japan. To find common genetic variants contributing to NTG in Japanese patients, we conducted a genome-wide association study (GWAS). We performed the first screening for 531,009 autosomal SNPs with a discovery cohort of 286 cases and 557 controls, and then a second screening for the top 30 suggestive loci in an independent cohort of 183 cases and 514 controls. Our findings identified a significantly associated SNP; rs523096 [combined p-value = 7.40× 10⁻⁸, odds ratio (OR)  = 2.00 with 95% confidence interval (CI) 1.55–2.58] located 10 kbp upstream of CDKN2B on chromosome 9p21. Moreover, analysis of another independent case-control set successfully replicated the results of the screening studies (combined values of all 3 stages p = 4.96 × 10⁻¹¹, OR  = 2.13 with 95% CI 1.69–2.68). The SNPs near rs523096 were recently reported to be associated with OAG associated with elevated IOP in primary open-angle glaucoma (POAG), the predominant subtype of glaucoma in Caucasian populations. Our results revealed that the 9p21 locus is also associated with NTG in Japanese. In addition, we identified SNPs more strongly associated with NTG.

Genome-wide association study confirming association of HLA-DP with protection against chronic hepatitis B and viral clearance in Japanese and Korean

Hepatitis B virus (HBV) infection can lead to serious liver diseases, including liver cirrhosis (LC) and hepatocellular carcinoma (HCC); however, about 85-90% of infected individuals become inactive carriers with sustained biochemical remission and very low risk of LC or HCC. To identify host genetic factors contributing to HBV clearance, we conducted genome-wide association studies (GWAS) and replication analysis using samples from HBV carriers and spontaneously HBV-resolved Japanese and Korean individuals. Association analysis in the Japanese and Korean data identified the HLA-DPA1 and HLA-DPB1 genes with P(meta) = 1.89×10⁻¹² for rs3077 and P(meta) = 9.69×10⁻¹⁰ for rs9277542. We also found that the HLA-DPA1 and HLA-DPB1 genes were significantly associated with protective effects against chronic hepatitis B (CHB) in Japanese, Korean and other Asian populations, including Chinese and Thai individuals (P(meta) = 4.40×10⁻¹⁹ for rs3077 and P(meta) = 1.28×10⁻¹⁵ for rs9277542). These results suggest that the associations between the HLA-DP locus and the protective effects against persistent HBV infection and with clearance of HBV were replicated widely in East Asian populations; however, there are no reports of GWAS in Caucasian or African populations. Based on the GWAS in this study, there were no significant SNPs associated with HCC development. To clarify the pathogenesis of CHB and the mechanisms of HBV clearance, further studies are necessary, including functional analyses of the HLA-DP molecule.

Response-guided therapy for patients with chronic hepatitis who have high viral loads of hepatitis C virus genotype 2

AIM

  We evaluated the efficacy of response-guided therapy in patients with hepatitis C virus (HCV) genotype 2.

METHODS

  We studied 105 patients with an HCV genotype 2 load of higher than 5.0 Log IU/mL who received more than 75% of the target dose of pegylated interferon plus ribavirin. Among patients with rapid viral response (RVR; no HCV RNA detected at week 4), 14 selected 16 weeks of therapy (group A), and 28 selected 24 weeks of therapy (group B). Among non-RVR patients, 40 selected 24 weeks of therapy (group C), and 19 selected 48 weeks of therapy (group D).

RESULTS

  All patients in group A and B achieved a sustained viral response (SVR). Clinical characteristics did not differ significantly between groups C and D. However, the proportion of patients in whom HCV RNA disappeared at a later week after starting treatment was higher in group D (P = 0.0578). SVR rate was 73% in C, and 79% in D. Among patients in whom HCV RNA disappeared between weeks 5 and 8, SVR was achieved in 28 (82%) of 34 patients in C and 10 (91%) of 11 patients in D. Among patients whose HCV RNA disappeared between weeks 9 and 12, SVR was achieved in one (20%) of five patients in C and five (63%) of eight patients in D (not statistically significant).

CONCLUSIONS

  16 weeks of combination therapy could achieve an adequate antiviral effect for RVR patients. Extending therapy could not significantly improve SVR rate in non-RVR patients.

Association of IFNGR2 gene polymorphisms with pulmonary tuberculosis among the Vietnamese

Interferon-γ (IFN-γ) is a key molecule of T helper 1 (Th1)-immune response against tuberculosis (TB), and rare genetic defects of IFN-γ receptors cause disseminated mycobacterial infection. The aim of the present study was to investigate whether genetic polymorphisms found in the Th1-immune response genes play a role in TB. In our study, DNA samples were collected from two series of cases including 832 patients with new smear-positive TB and 506 unrelated individuals with no history of TB in the general population of Hanoi, Vietnam. Alleles of eight microsatellite markers located around Th1-immune response-related genes and single nucleotide polymorphisms near the promising microsatellites were genotyped. A set of polymorphisms within the interferon gamma receptor 2 gene (IFNGR2) showed a significant association with protection against TB (P = 0.00054). Resistant alleles tend to be less frequently found in younger age at diagnosis (P = 0.011). Luciferase assays revealed high transcriptional activity of the promoter segment in linkage disequilibrium with resistant alleles. We conclude that the polymorphisms of IFNGR2 may confer resistance to the TB development of newly infected individuals. Contribution of the genetic factors to TB appeared to be different depending on age at diagnosis.

Highly parallel and short-acting amplification with locus-specific primers to detect single nucleotide polymorphisms by the DigiTag2 assay

The DigiTag2 assay enables analysis of a set of 96 SNPs using Kapa 2GFast HotStart DNA polymerase with a new protocol that has a total running time of about 7 hours, which is 6 hours shorter than the previous protocol. Quality parameters (conversion rate, call rate, reproducibility and concordance) were at the same levels as when genotype calls were acquired using the previous protocol. Multiplex PCR with 192 pairs of locus-specific primers was available for target preparation in the DigiTag2 assay without the optimization of reaction conditions, and quality parameters had the same levels as those acquired with 96-plex PCR. The locus-specific primers were able to achieve sufficient (concentration of target amplicon ≥5 nM) and specific (concentration of unexpected amplicons <2 nM) amplification within 2 hours, were also able to achieve detectable amplifications even when working in a 96-plex or 192-plex form. The improved DigiTag2 assay will be an efficient platform for screening an intermediate number of SNPs (tens to hundreds of sites) in the replication analysis after genome-wide association study. Moreover, highly parallel and short-acting amplification with locus-specific primers may thus facilitate widespread application to other PCR-based assays.

Association of TLR polymorphisms with development of tuberculosis in Indonesian females

Tuberculosis (TB) is caused by Mycobacterium tuberculosis and is a major cause of morbidity and mortality worldwide. Many candidate genes have been investigated for a possible association with TB. Toll-like receptors (TLRs) are known to play important roles in human innate immune systems. Polymorphisms in and functions of TLRs have been investigated to identify associations with specific infectious diseases, including TB. Here, we examined whether single-nucleotide polymorphisms (SNPs) in TLRs and genes in TLR signaling were associated with TB susceptibility in Indonesian and Vietnamese populations. A statistically significant association was observed between TB susceptibility in a classified Indonesian female group and rs352139, an SNP located in the intron of TLR9, using the genotype (P = 2.76E-04) and recessive (AA vs AG+GG, P = 2.48E-04, odds ratio = 1.827, 95% confidence interval = 1.321-2.526) models. Meta-analysis of the Indonesian and Vietnamese populations showed that rs352139 was significantly associated with TB in the recessive model. This finding indicated that a TLR9 polymorphism might have an important role in the susceptibility to M. tuberculosis in Asian populations.

Favorable factors for re-treatment with pegylated interferon α2a plus ribavirin in patients with high viral loads of genotype 1 hepatitis C virus

AIM

  Effect of re-treatment for pegylated interferon (PEG-IFN) plus ribavirin was not fully evaluated. We examined the effects of re-treatment with PEG-IFN plus ribavirin in patients with high viral loads of genotype 1 hepatitis C virus who failed to achieve a sustained virological response (SVR) with combination therapy.

METHODS

  We examined 38 patients who were re-treated with PEG-IFN α2a plus ribavirin for more than 60 weeks, among whom 14 were non-responders and 24 were relapsers after previous treatment with PEG-IFN α2b plus ribavirin. IL28B genotyping was done in 21 patients.

RESULTS

  The overall SVR rate was 34%. Analysis of baseline characteristics showed that the relapsers had a significantly higher SVR rate than the non-responders (50.0%, 12/24 vs. 7.1%, 1/14, respectively, P = 0.012) The SVR rates of re-treated patients who had turned hepatitis C virus (HCV) RNA-negative at weeks 8, 12, 24, and 48 of the previous therapy were 67% (4/6), 67% (4/6), 29% (2/7), and 25% (1/4), respectively. Re-treatment achieved an SVR in five of 12 patients with IL28B major alleles and three of nine patients with IL28B minor alleles. During the re-treatment, patients with complete viral suppression at week-12 achieved a significantly higher SVR rate (P = 0.001).

CONCLUSIONS

  Re-treatment with PEG-IFN α2a plus ribavirin therapy is effective in patients who relapse after a course of PEG-IFN α2b plus ribavirin therapy. Re-treatment is a particularly useful option for patients who achieve early viral clearance during previous therapy.

Genome-wide association study identified ITPA/DDRGK1 variants reflecting thrombocytopenia in pegylated interferon and ribavirin therapy for chronic hepatitis C

Hematologic abnormalities during current therapy with pegylated interferon and ribavirin (PEG-IFN/RBV) for chronic hepatitis C (CHC) often necessitate dose reduction and premature withdrawal from therapy. The aim of this study was to identify host factors associated with IFN-induced thrombocytopenia by genome-wide association study (GWAS). In the GWAS stage using 900K single-nucleotide polymorphism (SNP) microarrays, 303 Japanese CHC patients treated with PEG-IFN/RBV therapy were genotyped. One SNP (rs11697186) located on DDRGK1 gene on chromosome 20 showed strong associations in the minor-allele-dominant model with the decrease of platelet counts in response to PEG-IFN/RBV therapy [P = 8.17 × 10(-9); odds ratio (OR) = 4.6]. These associations were replicated in another sample set (n = 391) and the combined P-values reached 5.29 × 10(-17) (OR = 4.5). Fine mapping with 22 SNPs around DDRGK1 and ITPA genes showed that rs11697186 at the GWAS stage had a strong linkage disequilibrium with rs1127354, known as a functional variant in the ITPA gene. The ITPA-AA/CA genotype was independently associated with a higher degree of reduction in platelet counts at week 4 (P < 0.0001), as well as protection against the reduction in hemoglobin, whereas the CC genotype had significantly less reduction in the mean platelet counts compared with the AA/CA genotype (P < 0.0001 for weeks 2, 4, 8, 12), due to a reactive increase of the platelet count through weeks 1-4. Our present results may provide a valuable pharmacogenetic diagnostic tool for tailoring PEG-IFN/RBV dosing to minimize drug-induced adverse events.

Genetic relationships between Japanese native and commercial breeds using 70 chicken autosomal SNP genotypes by the DigiTag2 assay

Recently, single nucleotide polymorphisms (SNPs) have been used to identify genes or genomic regions responsible for economic traits, including genetic diseases in domestic animals, and to examine genetic diversity of populations. In this study, we genotyped 70 chicken autosomal SNPs using DigiTag2 assay to understand the genetic structure of the Japanese native chicken breeds Satsumadori and Ingie, and the relationship of these breeds with other established breeds, Rhode Island Red (RIR), commercial broiler and layer. Five breeds, each consisting of approximately 20 chickens, were subjected to the assay, revealing the following: Average expected heterozygosities of broiler, Satsumadori, RIR, layer and Ingie were 0.265, 0.254, 0.244, 0.179 and 0.176, respectively. Phylogenetic analysis using the concatenated 70 autosomal SNP genotypes distinguished all chickens and formed clusters of chickens belonging to the respective breeds. In addition, the 2-D scatter plot of the first two principal components was consistent with the phylogenic tree. Taken together with the pairwise F(st) distances, broiler and RIR were closely positioned near each other, while Ingie was positioned far from the other breeds. Structure analysis revealed that the probable number of genetic clusters (K) was six and four with maximum likelihood and ΔK values, respectively. The clustering with maximum likelihood revealed that, in addition to the clustering of the other five breeds, the Satsumadori was subdivided into two genetic clusters. The clustering with ΔK value indicated that the broiler and Rhode Island Red were assigned to the same genetic cluster.

Common variation in GPC5 is associated with acquired nephrotic syndrome

Severe proteinuria is a defining factor of nephrotic syndrome irrespective of the etiology. Investigation of congenital nephrotic syndrome has shown that dysfunction of glomerular epithelial cells (podocytes) plays a crucial role in this disease. Acquired nephrotic syndrome is also assumed to be associated with podocyte injury. Here we identify an association between variants in GPC5, encoding glypican-5, and acquired nephrotic syndrome through a genome-wide association study and replication analysis (P value under a recessive model (P(rec)) = 6.0 × 10(-11), odds ratio = 2.54). We show that GPC5 is expressed in podocytes and that the risk genotype is associated with higher expression. We further show that podocyte-specific knockdown and systemic short interfering RNA injection confers resistance to podocyte injury in mouse models of nephrosis. This study identifies GPC5 as a new susceptibility gene for nephrotic syndrome and implicates GPC5 as a promising therapeutic target for reducing podocyte vulnerability in glomerular disease.

Multiplex cDNA quantification method that facilitates the standardization of gene expression data

Microarray-based gene expression measurement is one of the major methods for transcriptome analysis. However, current microarray data are substantially affected by microarray platforms and RNA references because of the microarray method can provide merely the relative amounts of gene expression levels. Therefore, valid comparisons of the microarray data require standardized platforms, internal and/or external controls and complicated normalizations. These requirements impose limitations on the extensive comparison of gene expression data. Here, we report an effective approach to removing the unfavorable limitations by measuring the absolute amounts of gene expression levels on common DNA microarrays. We have developed a multiplex cDNA quantification method called GEP-DEAN (Gene expression profiling by DCN-encoding-based analysis). The method was validated by using chemically synthesized DNA strands of known quantities and cDNA samples prepared from mouse liver, demonstrating that the absolute amounts of cDNA strands were successfully measured with a sensitivity of 18 zmol in a highly multiplexed manner in 7 h.

Pre-treatment prediction of response to pegylated-interferon plus ribavirin for chronic hepatitis C using genetic polymorphism in IL28B and viral factors

BACKGROUND & AIMS

Pegylated interferon and ribavirin (PEG-IFN/RBV) therapy for chronic hepatitis C virus (HCV) genotype 1 infection is effective in 50% of patients. Recent studies revealed an association between the IL28B genotype and treatment response. We aimed to develop a model for the pre-treatment prediction of response using host and viral factors.

METHODS

Data were collected from 496 patients with HCV genotype 1 treated with PEG-IFN/RBV at five hospitals and universities in Japan. IL28B genotype and mutations in the core and IFN sensitivity determining region (ISDR) of HCV were analyzed to predict response to therapy. The decision model was generated by data mining analysis.

RESULTS

The IL28B polymorphism correlated with early virological response and predicted null virological response (NVR) (odds ratio=20.83, p<0.0001) and sustained virological response (SVR) (odds ratio=7.41, p<0.0001) independent of other covariates. Mutations in the ISDR predicted relapse and SVR independent of IL28B. The decision model revealed that patients with the minor IL28B allele and low platelet counts had the highest NVR (84%) and lowest SVR (7%), whereas those with the major IL28B allele and mutations in the ISDR or high platelet counts had the lowest NVR (0-17%) and highest SVR (61-90%). The model had high reproducibility and predicted SVR with 78% specificity and 70% sensitivity.

CONCLUSIONS

The IL28B polymorphism and mutations in the ISDR of HCV were significant pre-treatment predictors of response to PEG-IFN/RBV. The decision model, including these host and viral factors may support selection of optimum treatment strategy for individual patients.