Licht Toyo-oka

Pathogen lineage-based genome-wide association study identified CD53 as susceptible locus in tuberculosis.

Tuberculosis (TB) is known to be affected by host genetic factors. We reported a specific genetic risk factor through a genome-wide association study (GWAS) that focused on young age onset TB. In this study, we further focused on the heterogeneity of Mycobacterium tuberculosis (M. tb) lineages and assessed its possible interaction with age at onset on host genetic factors. We identified the pathogen lineage in 686 Thai TB cases and GWAS stratified by both infected pathogen lineage information and age at onset revealed a genome-wide significant association of one single-nucleotide polymorphism (SNP) on chromosome 1p13, which was specifically associated with non-Beijing lineage-infected old age onset cases (P=2.54E-08, OR=1.74 (95% CI=1.43–2.12)), when we compared them to the population-matched healthy controls. This SNP locates near the CD53 gene, which encodes a leukocyte surface glycoprotein. Interestingly, the expression of CD53 was also correlated with the patients’ active TB status. This is the first report of a pathogen lineage-based genome-wide association study. The results suggested that host genetic risk in TB is depended upon the pathogen genetic background and demonstrate the importance of analyzing the interaction between host and pathogen genomes in TB.

Strain-based HLA association analysis identified HLA-DRB1*09:01 associated with modern strain tuberculosis

Tuberculosis (TB) occurs as a result of complex interactions between the host immune system and pathogen virulence factors. Human leukocyte antigen (HLA) class II molecules play an important role in the host immune system. However, no study has assessed the association between HLA class II genes and susceptibility to TB caused by specific strains. This study investigated the possible association of HLA class II genes with TB caused by modern and ancient Mycobacterium tuberculosis (MTB). The study included 682 patients with TB and 836 control subjects who were typed for HLA‐DRB1 and HLA‐DQB1 alleles. MTB strains were classified using a large sequence polymorphism typing method. Association analysis was performed using common HLA alleles and haplotypes in different MTB strains. HLA association analysis of patients infected with modern MTB strains showed significant association for HLA‐DRB1*09:01 (odds ratio [OR] = 1.82; P‐value = 9.88 × 10−4) and HLA‐DQB1*03:03 alleles (OR = 1.76; P‐value = 1.31 × 10−3) with susceptibility to TB. Haplotype analysis confirmed that these alleles were in strong linkage disequilibrium and did not exert an interactive effect. Thus, the results of this study showed an association between HLA class II genes and susceptibility to TB caused by modern MTB strains, suggesting the importance of strain‐specific analysis to determine susceptibility genes associated with TB.

Genome-Wide Association Study Identifies Risk Variants for Lichen Planus in Patients With Hepatitis C Virus Infection

Background & Aims There is a close relationship between hepatitis C virus (HCV) infection and lichen planus, a chronic inflammatory mucocutaneous disease. We performed a genome‐wide association study (GWAS) to identify genetic variants associated with HCV‐related lichen planus. Methods We conducted a GWAS of 261 patients with HCV infection treated at a tertiary medical center in Japan from October 2007 through January 2013; a total of 71 had lichen planus and 190 had normal oral mucosa. We validated our findings in a GWAS of 38 patients with HCV‐associated lichen planus and 7 HCV‐infected patients with normal oral mucosa treated at a medical center in Italy. Results Single‐nucleotide polymorphisms in NRP2 (rs884000) and IGFBP4 (rs538399) were associated with risk of HCV‐associated lichen planus (P < 1 × 10−4). We also found an association between a single‐nucleotide polymorphism in the HLA‐DR/DQ genes (rs9461799) and susceptibility to HCV‐associated lichen planus. The odds ratios for the minor alleles of rs884000, rs538399, and rs9461799 were 3.25 (95% confidence interval, 1.95–5.41), 0.40 (95% confidence interval, 0.25–0.63), and 2.15 (95% confidence interval, 1.41–3.28), respectively. Conclusions In a GWAS of Japanese patients with HCV infection, we replicated associations between previously reported polymorphisms in HLA class II genes and risk for lichen planus. We also identified single‐nucleotide polymorphisms in NRP2 and IGFBP4 loci that increase and reduce risk of lichen planus, respectively. These genetic variants might be used to identify patients with HCV infection who are at risk for lichen planus.

The first genome-wide association study identifying new susceptibility loci for obstetric antiphospholipid syndrome

Antiphospholipid syndrome (APS) is the most important treatable cause of recurrent pregnancy loss. The live birth rate is limited to only 70–80% in patients with APS undergoing established anticoagulant therapy. Lupus anticoagulant (LA), but not anticardiolipin antibody (aCL), was found to predict adverse pregnancy outcome. Recent genome-wide association studies (GWAS) of APS focusing on aCL have shown that several molecules may be involved. This is the first GWAS for obstetric APS focusing on LA. A GWAS was performed to compare 115 Japanese patients with obstetric APS, diagnosed according to criteria of the International Congress on APS, and 419 healthy individuals. Allele or genotype frequencies were compared in a total of 426 344 single-nucleotide polymorphisms (SNPs). Imputation analyses were also performed for the candidate regions detected by the GWAS. One SNP (rs2288493) located on the 3′-UTR of TSHR showed an experiment-wide significant APS association (P=7.85E-08, OR=6.18) under a recessive model after Bonferroni correction considering the number of analyzed SNPs. Another SNP (rs79154414) located around the C1D showed a genome-wide significant APS association (P=4.84E-08, OR=6.20) under an allelic model after applying the SNP imputation. Our findings demonstrate that a specific genotype of TSHR and C1D genes can be a risk factor for obstetric APS.

Fine-mapping analysis revealed complex pleiotropic effect and tissue-specific regulatory mechanism of TNFSF15 in primary biliary cholangitis, Crohn’s disease and leprosy

Genetic polymorphism within the 9q32 locus is linked with increased risk of several diseases, including Crohn’s disease (CD), primary biliary cholangitis (PBC) and leprosy. The most likely disease-causing gene within 9q32 is TNFSF15, which encodes the pro-inflammatory cytokine TNF super-family member 15, but it was unknown whether these disparate diseases were associated with the same genetic variance in 9q32, and how variance within this locus might contribute to pathology. Using genetic data from published studies on CD, PBC and leprosy we revealed that bearing a T allele at rs6478108/rs6478109 (r2 = 1) or rs4979462 was significantly associated with increased risk of CD and decreased risk of leprosy, while the T allele at rs4979462 was associated with significantly increased risk of PBC. In vitro analyses showed that the rs6478109 genotype significantly affected TNFSF15 expression in cells from whole blood of controls, while functional annotation using publicly-available data revealed the broad cell type/tissue-specific regulatory potential of variance at rs6478109 or rs4979462. In summary, we provide evidence that variance within TNFSF15 has the potential to affect cytokine expression across a range of tissues and thereby contribute to protection from infectious diseases such as leprosy, while increasing the risk of immune-mediated diseases including CD and PBC.

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.

Polymorphisms in the TMEM132D region are associated with panic disorder in HLA-DRB1*13:02-negative individuals of a Japanese population

We herein report an association between TMEM132D and panic disorder (PD) in a Japanese population, evaluating the effects of HLA-DRB1*13:02, which we previously reported as a susceptibility genetic factor for PD. SNPs in TMEM132D showed significant associations with PD in subjects without HLA-DRB1*13:02 (rs4759997; P=5.02×10−6, odds ratio=1.50) but not in those with the HLA allele. TMEM132D might have a role in the development of PD in subjects without HLA-DRB1*13:02.

Evidence for Host-Bacterial Co-evolution via Genome Sequence Analysis of 480 Thai Mycobacterium tuberculosis Lineage 1 Isolates

Tuberculosis presents a global health challenge. Mycobacterium tuberculosis is divided into several lineages, each with a different geographical distribution. M. tuberculosis lineage 1 (L1) is common in the high-burden areas in East Africa and Southeast Asia. Although the founder effect contributes significantly to the phylogeographic profile, co-evolution between the host and M. tuberculosis may also play a role. Here, we reported the genomic analysis of 480 L1 isolates from patients in northern Thailand. The studied bacterial population was genetically diverse, allowing the identification of a total of 18 sublineages distributed into three major clades. The majority of isolates belonged to L1.1 followed by L1.2.1 and L1.2.2. Comparison of the single nucleotide variant (SNV) phylogenetic tree and the clades defined by spoligotyping revealed some monophyletic clades representing EAI2_MNL, EAI2_NTM and EAI6_BGD1 spoligotypes. Our work demonstrates that ambiguity in spoligotype assignment could be partially resolved if the entire DR region is investigated. Using the information to map L1 diversity across Southeast Asia highlighted differences in the dominant strain-types in each individual country, despite extensive interactions between populations over time. This finding supported the hypothesis that there is co-evolution between the bacteria and the host, and have implications for tuberculosis disease control.