Genome-wide host-pathogen analyses reveal genetic interaction points in tuberculosis disease.

The genetics underlying tuberculosis (TB) pathophysiology are poorly understood. Human genome-wide association studies have failed so far to reveal reproducible susceptibility loci, attributed in part to the influence of the underlying Mycobacterium tuberculosis (Mtb) bacterial genotype on the outcome of the infection. Several studies have found associations of human genetic polymorphisms with Mtb phylo-lineages, but studies analysing genome-genome interactions are needed. By implementing a phylogenetic tree-based Mtb-to-human analysis for 714 TB patients from Thailand, we identify eight putative genetic interaction points (P < 5 × 10-8) including human loci DAP and RIMS3, both linked to the IFNγ cytokine and host immune system, as well as FSTL5, previously associated with susceptibility to TB. Many of the corresponding Mtb markers are lineage specific. The genome-to-genome analysis reveals a complex interactome picture, supports host-pathogen adaptation and co-evolution in TB, and has potential applications to large-scale studies across many TB endemic populations matched for host-pathogen genomic diversity.

TogoVar: A comprehensive Japanese genetic variation database.

TogoVar ( https://togovar.org ) is a database that integrates allele frequencies derived from Japanese populations and provides annotations for variant interpretation. First, a scheme to reanalyze individual-level genome sequence data deposited in the Japanese Genotype-phenotype Archive (JGA), a controlled-access database, was established to make allele frequencies publicly available. As more Japanese individual-level genome sequence data are deposited in JGA, the sample size employed in TogoVar is expected to increase, contributing to genetic study as reference data for Japanese populations. Second, public datasets of Japanese and non-Japanese populations were integrated into TogoVar to easily compare allele frequencies in Japanese and other populations. Each variant detected in Japanese populations was assigned a TogoVar ID as a permanent identifier. Third, these variants were annotated with molecular consequence, pathogenicity, and literature information for interpreting and prioritizing variants. Here, we introduce the newly developed TogoVar database that compares allele frequencies among Japanese and non-Japanese populations and describes the integrated annotations.

The NBDC-DDBJ imputation server facilitates the use of controlled access reference panel datasets in Japan.

Accurate genotype imputation requires large-scale reference panel datasets. When conducting genotype imputation on the Japanese population, researchers can use such datasets under collaborative studies or controlled access conditions in public databases. We developed the NBDC-DDBJ imputation server, which securely provides users with a web user interface to execute genotype imputation on the server. Our benchmarking analysis showed that the accuracy of genotype imputation was improved by leveraging controlled access datasets to increase the number of haplotypes available for analysis compared to using publicly available reference panels such as the 1000 Genomes Project. The NBDC-DDBJ imputation server facilitates the use of controlled access datasets for accurate genotype imputation.

Haplotype-specific PCR for NAT2 diplotyping.

N-acetyltransferase 2 (NAT2) is an enzyme that acetylates many kinds of drugs, including the antituberculosis drug isoniazid. The NAT2 gene is highly diverse across populations. An individual can be classified as having a slow acetylator (SA), an intermediate acetylator (IA), or a rapid acetylator (RA) phenotype based on its two haplotypes (diplotype) of NAT2. SA individuals are at a higher risk for isoniazid-induced hepatitis, while the RA phenotype contributes to failure in tuberculosis treatment. Being able to predict individual NAT2 phenotypes is important for dose adjustment of isoniazid. NAT2 haplotypes are commonly determined via an indirect method of statistical haplotype inference from SNP genotyping. Here, we report a direct NAT2 haplotyping method using haplotype-specific PCR (HS-PCR) for the 6 most commonly found NAT2 haplotypes: NAT2*4, NAT2*5B, NAT2*6A, NAT2*7B, NAT2*12A, and NAT2*13A. Validation of this HS-PCR method via comparison with a sequencing method in 650 Thai DNA samples (107 RA, 279 IA, and 264 SA samples) showed a concordance rate for diplotype calls of 99.23% (645/650 samples). The discordant results in 5 samples were due to 3 rare NAT2 haplotypes: NAT*5C (n = 3), NAT2*7C (n = 1), and NAT2*11A (n = 1). This novel HS-PCR method allows direct NAT2 diplotyping, enabling the implementation of NAT2 acetylator phenotypes in clinical pharmacogenetic testing.

An optimized genomic VCF workflow for precise identification of Mycobacterium tuberculosis cluster from cross-platform whole genome sequencing data.

Whole-genome sequencing (WGS) data allow for an inference of Mycobacterium tuberculosis (Mtb) clusters by using a pairwise genetic distance of ≤12 single nucleotide polymorphisms (SNPs) as a threshold. However, a problem of discrepancies in numbers of SNPs and genetic distance measurement is a great concern when combining WGS data from different next generation sequencing (NGS) platforms. We performed SNP variant calling on WGS data of 9 multidrug-resistant (MDR-TB), 3 extensively drug-resistant tuberculosis (XDR-TB) and a standard M. tuberculosis strain H37Rv from an Illumina/NextSeq500 and an Ion Torrent PGM. Variant calls were obtained using four different common variant calling tools, including Genome Analysis Toolkit (GATK) HaplotypeCaller (GATK-VCF workflow), GATK HaplotypeCaller and GenotypeGVCFs (GATK-GVCF workflow), SAMtools, and VarScan 2. Cross-platform pairwise SNP differences, minimum spanning networks and average nucleotide identity (ANI) were analysed to measure performance of the variant calling tools. Minimum pairwise SNP differences ranged from 2 to 14 SNPs when using GVCF workflow while maximum pairwise SNP differences ranged from 7 to 158 SNPs when using VarScan 2. ANI comparison between SNPs data from NextSeq500 and PGM of MDR-TB and XDR-TB showed maximum ANI of 99.7% and 99.0%, respectively, with GVCF workflow while the other SNP calling results showed lower ANI in a range of 98.6% to 95.1%. In this study, we suggest that the GVCF workflow showed the best performing variant caller to avoid cross-platform pairwise SNP differences.

Treatment of two mitochondrial disease patients with a combination of febuxostat and inosine that enhances cellular ATP.

Since mitochondria are energy-generating micro-organisms, most of the disorders in patients with mitochondrial diseases (mt-disease) are considered secondary to defects in ATP synthesis, although some other factors such as reactive oxygen species may be involved. A simultaneous oral administration of febuxostat and inosine was reported to elevate both hypoxanthine and ATP levels in peripheral blood. Based on those results, we attempted co-administration of febuxostat and inosine in two patients with mitochondrial disease: one patient with mitochondrial cardiomyopathy and the other patient with mitochondrial diabetes. In the former case, brain natriuretic peptide (BNP), which is a specific marker for heart failure, was decreased by 31%, and in the latter case, the insulinogenic index increased 3.1 times, suggesting the favorable action of the treatment. Considering that there is no effective treatment available for this disorder, the present therapy may be quite useful for the management of patients with mitochondrial diseases.

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.

NAT2 ultra-slow acetylator and risk of anti-tuberculosis drug-induced liver injury: a genotype-based meta-analysis.

BACKGROUND: NAT2 slow acetylator is a confirmed risk of anti-tuberculosis drug-induced liver injury (ATDILI). However, NAT2 ultra-slow acetylators, a new refinement among NAT2 slow acetylators, have been recently proposed. The patients with NAT2 genotypes of *6A/*6A, *6A/*7B and *7B/*7B are referred to in this group. OBJECTIVE: We aim to prove an association of the NAT2 ultra-slow acetylators with the risk of ATDILI. MATERIALS AND METHODS: Systematic review and meta-analysis were performed based on each NAT2 genotype and risk of ATDILI cases and also new classification of the ultra-slow acetylators up to 31 October 2016. Meta-analysis of 18 studies with 822 ATDILI cases and 4630 controls was carried out in the RevMan software, version 5.3 with fixed-effect (low heterogeneity) and random effect (moderate to high heterogeneity) methods. RESULTS: The strong associations between each NAT2 slow acetylator genotypes and ATDILI were confirmed in meta-analysis except for NAT2*5B/*5B [odds ratio (OR): 1.69; 95% confidence interval (CI): 0.96-2.95; P=0.0679]. The NAT2 ultra-slow acetylators contribute to higher risk of ATDILI (OR: 3.60; 95% CI: 2.30-5.63; P=1.76E-08) than all NAT2 slow acetylators (OR: 2.80; 95% CI: 2.20-3.57; P=5.73E-18) as well as fast acetylators. Additional in-vitro study using isoniazid as a substrate supports the existence of ultra-slow acetylator alleles (NAT2*6A and NAT2*7B). CONCLUSION: This is the first meta-analysis of NAT2 and the risk of ATDILI at the genotypic level. The result demonstrated that NAT2 ultra-slow acetylator genotypes will have the most effect on the increased risk of ATDILI.

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.

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.

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.

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.

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 (r(2) = 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.

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.