Tatsuhiko Tsunoda

Functional haplotypes of PADI4, encoding citrullinating enzyme peptidylarginine deiminase 4, are associated with rheumatoid arthritis.

Individuals with rheumatoid arthritis frequently have autoantibodies to citrullinated peptides, suggesting the involvement of the peptidylarginine deiminases citrullinating enzymes (encoded by PADI genes) in rheumatoid arthritis. Previous linkage studies have shown that a susceptibility locus for rheumatoid arthritis includes four PADI genes but did not establish which PADI gene confers susceptibility to rheumatoid arthritis. We used a case-control linkage disequilibrium study to show that PADI type 4 is a susceptibility locus for rheumatoid arthritis (P = 0.000008). PADI4 was expressed in hematological and rheumatoid arthritis synovial tissues. We also identified a haplotype of PADI4 associated with susceptibility to rheumatoid arthritis that affected stability of transcripts and was associated with levels of antibody to citrullinated peptide in sera from individuals with rheumatoid arthritis. Our results imply that the PADI4 haplotype associated with susceptibility to rheumatoid arthritis increases production of citrullinated peptides acting as autoantigens, resulting in heightened risk of developing the disease.

Genome-wide association study identifies common variants at four loci as genetic risk factors for Parkinson's disease.

To identify susceptibility variants for Parkinsons disease (PD), we performed a genome-wide association study (GWAS) and two replication studies in a total of 2,011 cases and 18,381 controls from Japan. We identified a new susceptibility locus on 1q32 (P = 1.52 × 10−12) and designated this as PARK16, and we also identified BST1 on 4p15 as a second new risk locus (P = 3.94 × 10−9). We also detected strong associations at SNCA on 4q22 (P = 7.35 × 10−17) and LRRK2 on 12q12 (P = 2.72 × 10−8), both of which are implicated in autosomal dominant forms of parkinsonism. By comparing results of a GWAS performed on individuals of European ancestry, we identified PARK16, SNCA and LRRK2 as shared risk loci for PD and BST1 and MAPT as loci showing population differences. Our results identify two new PD susceptibility loci, show involvement of autosomal dominant parkinsonism loci in typical PD and suggest that population differences contribute to genetic heterogeneity in PD.

Functional SNPs in the lymphotoxin-α gene that are associated with susceptibility to myocardial infarction

By means of a large-scale, case-control association study using 92,788 gene-based single-nucleotide polymorphism (SNP) markers, we identified a candidate locus on chromosome 6p21 associated with susceptibility to myocardial infarction. Subsequent linkage-disequilibrium (LD) mapping and analyses of haplotype structure showed significant associations between myocardial infarction and a single 50 kb halpotype comprised of five SNPs in LTA (encoding lymphotoxin-α), NFKBIL1 (encoding nuclear factor of κ light polypeptide gene enhancer in B cells, inhibitor-like 1) and BAT1 (encoding HLA-B associated transcript 1). Homozygosity with respect to each of the two SNPs in LTA was significantly associated with increased risk for myocardial infarction (odds ratio = 1.78, χ2 = 21.6, P = 0.00000033; 1,133 affected individuals versus 1,006 controls). In vitro functional analyses indicated that one SNP in the coding region of LTA, which changed an amino-acid residue from threonine to asparagine (Thr26Asn), effected a twofold increase in induction of several cell-adhesion molecules, including VCAM1, in vascular smooth-muscle cells of human coronary artery. Moreover, the SNP, in intron 1 of LTA, enhanced the transcriptional level of LTA. These results indicate that variants in the LTA are risk factors for myocardial infraction and implicate LTA in the pathogenesis of the disorder.

SNPs in KCNQ1 are associated with susceptibility to type 2 diabetes in East Asian and European populations

We conducted a genome-wide association study using 207,097 SNP markers in Japanese individuals with type 2 diabetes and unrelated controls, and identified KCNQ1 (potassium voltage-gated channel, KQT-like subfamily, member 1) to be a strong candidate for conferring susceptibility to type 2 diabetes. We detected consistent association of a SNP in KCNQ1 (rs2283228) with the disease in several independent case-control studies (additive model P = 3.1 × 10−12; OR = 1.26, 95% CI = 1.18–1.34). Several other SNPs in the same linkage disequilibrium (LD) block were strongly associated with type 2 diabetes (additive model: rs2237895, P = 7.3 × 10−9; OR = 1.32, 95% CI = 1.20–1.45, rs2237897, P = 6.8 × 10−13; OR = 1.41, 95% CI = 1.29–1.55). The association of these SNPs with type 2 diabetes was replicated in samples from Singaporean (additive model: rs2237895, P = 8.5 × 10−3; OR = 1.14, rs2237897, P = 2.4 × 10−4; OR = 1.22) and Danish populations (additive model: rs2237895, P = 3.7 × 10−11; OR = 1.24, rs2237897, P = 1.2 × 10−4; OR = 1.36).

An intronic SNP in a RUNX1 binding site of SLC22A4 , encoding an organic cation transporter, is associated with rheumatoid arthritis

Rheumatoid arthritis is a common inflammatory disease with complex genetic components. We investigated the genetic contribution of the cytokine gene cluster in chromosome 5q31 to susceptibility to rheumatoid arthritis in the Japanese population by case-control linkage disequilibrium (LD) mapping using single nucleotide polymorphisms (SNPs). Here we report that there is significant association between rheumatoid arthritis and the organic cation transporter gene SLC22A4 (P = 0.000034). We show that expression of SLC22A4 is specific to hematological and immunological tissues and that SLC22A4 is also highly expressed in the inflammatory joints of mice with collagen-induced arthritis. A SNP affects the transcriptional efficiency of SLC22A4 in vitro, owing to an allelic difference in affinity to Runt-related transcription factor 1 (RUNX1), a transcriptional regulator in the hematopoietic system. A SNP in RUNX1 is also strongly associated with rheumatoid arthritis (P = 0.00035). Our data indicate that the regulation of SLC22A4 expression by RUNX1 is associated with susceptibility to rheumatoid arthritis, which may represent an example of an epistatic effect of two genes on this disorder.

Whole-genome sequencing of liver cancers identifies etiological influences on mutation patterns and recurrent mutations in chromatin regulators

Hepatocellular carcinoma (HCC) is the third leading cause of cancer-related death worldwide. We sequenced and analyzed the whole genomes of 27 HCCs, 25 of which were associated with hepatitis B or C virus infections, including two sets of multicentric tumors. Although no common somatic mutations were identified in the multicentric tumor pairs, their whole-genome substitution patterns were similar, suggesting that these tumors developed from independent mutations, although their shared etiological backgrounds may have strongly influenced their somatic mutation patterns. Statistical and functional analyses yielded a list of recurrently mutated genes. Multiple chromatin regulators, including ARID1A, ARID1B, ARID2, MLL and MLL3, were mutated in ∼50% of the tumors. Hepatitis B virus genome integration in the TERT locus was frequently observed in a high clonal proportion. Our whole-genome sequencing analysis of HCCs identified the influence of etiological background on somatic mutation patterns and subsequent carcinogenesis, as well as recurrent mutations in chromatin regulators in HCCs.

Integrated molecular analysis of clear-cell renal cell carcinoma

Clear-cell renal cell carcinoma (ccRCC) is the most prevalent kidney cancer and its molecular pathogenesis is incompletely understood. Here we report an integrated molecular study of ccRCC in which ≥100 ccRCC cases were fully analyzed by whole-genome and/or whole-exome and RNA sequencing as well as by array-based gene expression, copy number and/or methylation analyses. We identified a full spectrum of genetic lesions and analyzed gene expression and DNA methylation signatures and determined their impact on tumor behavior. Defective VHL-mediated proteolysis was a common feature of ccRCC, which was caused not only by VHL inactivation but also by new hotspot TCEB1 mutations, which abolished Elongin C–VHL binding, leading to HIF accumulation. Other newly identified pathways and components recurrently mutated in ccRCC included PI3K-AKT-mTOR signaling, the KEAP1-NRF2-CUL3 apparatus, DNA methylation, p53-related pathways and mRNA processing. This integrated molecular analysis unmasked new correlations between DNA methylation, gene mutation and/or gene expression and copy number profiles, enabling the stratification of clinical risks for patients with ccRCC.

A genome-wide association study identifies variants in the HLA-DP locus associated with chronic hepatitis B in Asians

Chronic hepatitis B is a serious infectious liver disease that often progresses to liver cirrhosis and hepatocellular carcinoma; however, clinical outcomes after viral exposure vary enormously among individuals. Through a two-stage genome-wide association study using 786 Japanese chronic hepatitis B cases and 2,201 controls, we identified a significant association of chronic hepatitis B with 11 SNPs in a region including HLA-DPA1 and HLA-DPB1. We validated these associations by genotyping two SNPs from the region in three additional Japanese and Thai cohorts consisting of 1,300 cases and 2,100 controls (combined P = 6.34 × 10−39 and 2.31 × 10−38, OR = 0.57 and 0.56, respectively). Subsequent analyses revealed risk haplotypes (HLA-DPA1*0202-DPB1*0501 and HLA-DPA1*0202-DPB1*0301, OR = 1.45 and 2.31, respectively) and protective haplotypes (HLA-DPA1*0103-DPB1*0402 and HLA-DPA1*0103-DPB1*0401, OR = 0.52 and 0.57, respectively). Our findings show that genetic variants in the HLA-DP locus are strongly associated with risk of persistent infection with hepatitis B virus.

Meta-analysis identifies six new susceptibility loci for atrial fibrillation

Atrial fibrillation is a highly prevalent arrhythmia and a major risk factor for stroke, heart failure and death. We conducted a genome-wide association study (GWAS) in individuals of European ancestry, including 6,707 with and 52,426 without atrial fibrillation. Six new atrial fibrillation susceptibility loci were identified and replicated in an additional sample of individuals of European ancestry, including 5,381 subjects with and 10,030 subjects without atrial fibrillation (P < 5 × 10−8). Four of the loci identified in Europeans were further replicated in silico in a GWAS of Japanese individuals, including 843 individuals with and 3,350 individuals without atrial fibrillation. The identified loci implicate candidate genes that encode transcription factors related to cardiopulmonary development, cardiac-expressed ion channels and cell signaling molecules.

A functional variant in FCRL3, encoding Fc receptor-like 3, is associated with rheumatoid arthritis and several autoimmunities

Rheumatoid arthritis is a common autoimmune disease with a complex genetic etiology. Here we identify a SNP in the promoter region of FCRL3, a member of the Fc receptor-like family, that is associated with susceptibility to rheumatoid arthritis (odds ratio = 2.15, P = 0.00000085). This polymorphism alters the binding affinity of nuclear factor-κB and regulates FCRL3 expression. We observed high FCRL3 expression on B cells and augmented autoantibody production in individuals with the disease-susceptible genotype. We also found associations between the SNP and susceptibility to autoimmune thyroid disease and systemic lupus erythematosus. FCRL3 may therefore have a pivotal role in autoimmunity.