Whole-genome and -exome resequencing using next-generation sequencers is a powerful approach for identifying genomic variations that are associated with diseases. However, systematic strategies for prioritizing causative variants from many candidates to explain the disease phenotype are still far from being established, because the population-specific frequency spectrum of genetic variation has not been characterized. Here, we have collected exomic genetic variation from 1208 Japanese individuals through a collaborative effort, and aggregated the data into a prevailing catalog. In total, we identified 156 622 previously unreported variants. The allele frequencies for the majority (88.8%) were lower than 0.5% in allele frequency and predicted to be functionally deleterious. In addition, we have constructed a Japanese-specific major allele reference genome by which the number of unique mapping of the short reads in our data has increased 0.045% on average. Our results illustrate the importance of constructing an ethnicity-specific reference genome for identifying rare variants. All the collected data were centralized to a newly developed database to serve as useful resources for exploring pathogenic variations. Public access to the database is available at http://www.genome.med.kyoto-u.ac.jp/SnpDB/.

The accurate typing of human leukocyte antigen (HLA) alleles is critical for a variety of medical applications, such as genomic studies of multifactorial diseases, including immune system and inflammation‐related disorders, and donor selection in organ transplantation and regenerative medicine. Here, we developed a new algorithm for determining HLA alleles using next‐generation sequencing (NGS) results. The method consists of constructing an extensive dictionary of HLA alleles, precise mapping of the NGS reads, and calculating a score based on weighted read counts to select the most suitable pair of alleles. The developed algorithm compares the score of all allele pairs, taking into account variation not only in the domain for antigen presentation (G‐DOMAIN), but also outside this domain. Using this method, HLA alleles could be determined with 6‐digit precision. We showed that our method was more accurate than other NGS‐based methods and revealed limitations of the conventional HLA typing technologies. Furthermore, we determined the complete genomic sequence of an HLA‐A‐like‐pseudogene when we assembled NGS reads that had caused arguable typing, and found its identity with HLA‐Y*02:01. The accuracy of the HLA‐A allele typing was improved after the HLA‐Y*02:01 sequence was included in the HLA allele dictionary.

Spermidine (SPD) delays age-related pathologies in various organisms. SPD supplementation overcame the impaired immunotherapy against tumors in aged mice by increasing mitochondrial function and activating CD8 ⁺ T cells. Treatment of naïve CD8 ⁺ T cells with SPD acutely enhanced fatty acid oxidation. SPD conjugated to beads bound to the mitochondrial trifunctional protein (MTP). In the MTP complex, synthesized and purified from Escherichia coli , SPD bound to the α and β subunits of MTP with strong affinity and allosterically enhanced their enzymatic activities. T cell–specific deletion of the MTP α subunit abolished enhancement of programmed cell death protein 1 (PD-1) blockade immunotherapy by SPD, indicating that MTP is required for SPD-dependent T cell activation.

Nonalcoholic fatty liver disease (NAFLD) includes a broad range of liver pathologies from simple steatosis to cirrhosis and fibrosis, in which a subtype accompanying hepatocyte degeneration and fibrosis is classified as nonalcoholic steatohepatitis (NASH). NASH accounts for approximately 10–30% of NAFLD and causes a higher frequency of liver-related death, and its progression of NASH has been considered to be complex involving multiple genetic factors interacting with the environment and lifestyle. To identify genetic factors related to NAFLD in the Japanese, we performed a genome-wide association study recruiting 529 histologically diagnosed NAFLD patients and 932 population controls. A significant association was observed for a cluster of SNPs in PNPLA3 on chromosome 22q13 with the strongest p-value of 1.4×10⁻¹⁰ (OR = 1.66, 95%CI: 1.43–1.94) for rs738409. Rs738409 also showed the strongest association (p = 3.6×10⁻⁶) with the histological classifications proposed by Matteoni and colleagues based on the degree of inflammation, ballooning degeneration, fibrosis and Mallory-Denk body. In addition, there were marked differences in rs738409 genotype distributions between type4 subgroup corresponding to NASH and the other three subgroups (p = 4.8×10⁻⁶, OR = 1.96, 95%CI: 1.47–2.62). Moreover, a subgroup analysis of NAFLD patients against controls showed a significant association of rs738409 with type4 (p = 1.7×10⁻¹⁶, OR = 2.18, 95%CI: 1.81–2.63) whereas no association was obtained for type1 to type3 (p = 0.41). Rs738409 also showed strong associations with three clinical traits related to the prognosis of NAFLD, namely, levels of hyaluronic acid (p = 4.6×10⁻⁴), HbA1c (p = 0.0011) and iron deposition in the liver (p = 5.6×10⁻⁴). With these results we clearly demonstrated that Matteoni type4 NAFLD is both a genetically and clinically different subset from the other spectrums of the disease and that the PNPLA3 gene is strongly associated with the progression of NASH in Japanese population.

Germ line DDX41 variants have been implicated in late-onset myeloid neoplasms (MNs). Despite an increasing number of publications, many important features of DDX41-mutated MNs remain to be elucidated. Here we performed a comprehensive characterization of DDX41-mutated MNs, enrolling a total of 346 patients with DDX41 pathogenic/likely-pathogenic (P/LP) germ line variants and/or somatic mutations from 9082 MN patients, together with 525 first-degree relatives of DDX41-mutated and wild-type (WT) patients. P/LP DDX41 germ line variants explained ∼80% of known germ line predisposition to MNs in adults. These risk variants were 10-fold more enriched in Japanese MN cases (n = 4461) compared with the general population of Japan (n = 20 238). This enrichment of DDX41 risk alleles was much more prominent in male than female (20.7 vs 5.0). P/LP DDX41 variants conferred a large risk of developing MNs, which was negligible until 40 years of age but rapidly increased to 49% by 90 years of age. Patients with myelodysplastic syndromes (MDS) along with a DDX41-mutation rapidly progressed to acute myeloid leukemia (AML), which was however, confined to those having truncating variants. Comutation patterns at diagnosis and at progression to AML were substantially different between DDX41-mutated and WT cases, in which none of the comutations affected clinical outcomes. Even TP53 mutations made no exceptions and their dismal effect, including multihit allelic status, on survival was almost completely mitigated by the presence of DDX41 mutations. Finally, outcomes were not affected by the conventional risk stratifications including the revised/molecular International Prognostic Scoring System. Our findings establish that MDS with DDX41-mutation defines a unique subtype of MNs that is distinct from other MNs.

Papillary thyroid cancer (PTC) among individuals exposed to radioactive iodine in their childhood or adolescence is a major internationally recognized health consequence of the Chernobyl accident. To identify genetic determinants affecting individual susceptibility to radiation-related PTC, we conducted a genome-wide association study employing Belarusian patients with PTC aged 0–18 years at the time of accident and age-matched Belarusian control subjects. Two series of genome scans were performed using independent sample sets, and association with radiation-related PTC was evaluated. Meta-analysis by the Mantel–Haenszel method combining the two studies identified four SNPs at chromosome 9q22.33 showing significant associations with the disease (Mantel–Haenszel P : mhp = 1.7 × 10 ⁻⁹ to 4.9 × 10 ⁻⁹ ). The association was further reinforced by a validation analysis using one of these SNP markers, rs965513, with a new set of samples (overall mhp = 4.8 × 10 ⁻¹² , OR = 1.65, 95% CI: 1.43–1.91). Rs965513 is located 57-kb upstream to FOXE1 , a thyroid-specific transcription factor with pivotal roles in thyroid morphogenesis and was recently reported as the strongest genetic risk marker of sporadic PTC in European populations. Of interest, no association was obtained between radiation-related PTC and rs944289 (mhp = 0.17) at 14p13.3 which showed the second strongest association with sporadic PTC in Europeans. These results show that the complex pathway underlying the pathogenesis may be partly shared by the two etiological forms of PTC, but their genetic components do not completely overlap each other, suggesting the presence of other unknown etiology-specific genetic determinants in radiation-related PTC.