Kyuto Sonehara, MD, PhD
曽根原 究人

†: equal contribution

1. Germline variants and mosaic chromosomal alterations affect COVID-19 vaccine immunogenicity

   †Sonehara K., †Uwamino Y., Saiki R., Takeshita M., Namba S., Uno S., and 56 more authors

   Cell Genomics, Mar 2025

   Abstract PubMed

   Vaccine immunogenicity is influenced by the vaccinee’s genetic background. Here, we perform a genome-wide association study of vaccine-induced SARS-CoV-2-specific immunoglobulin G (IgG) antibody titers and T cell immune responses in 1,559 mRNA-1273 and 537 BNT162b2 vaccinees of Japanese ancestry. SARS-CoV-2-specific antibody titers are associated with the immunoglobulin heavy chain (IGH) and major histocompatibility complex (MHC) locus, and T cell responses are associated with MHC. The lead variants at IGH contain a population-specific missense variant (rs1043109-C; p.Leu192Val) in the immunoglobulin heavy constant gamma 1 gene (IGHG1), with a strong decreasing effect (β = −0.54). Antibody-titer-associated variants modulate circulating immune regulatory proteins (e.g., LILRB4 and FCRL6). Age-related hematopoietic expanded mosaic chromosomal alterations (mCAs) affecting MHC and IGH also impair antibody production. MHC-/IGH-affecting mCAs confer infectious and immune disease risk, including sepsis and Graves’ disease. Impacts of expanded mosaic loss of chromosomes X/Y on these phenotypes were examined. Altogether, both germline and somatic mutations contribute to adaptive immunity functions.

2. Common and rare genetic variants predisposing females to unexplained recurrent pregnancy loss

   †Sonehara K., †Yano Y., Naito T., Goto S., Yoshihara H., Otani T., and 6 more authors

   Nature Communications 15:5744, Jul 2024

   Abstract PubMed

   Recurrent pregnancy loss (RPL) is a major reproductive health issue with multifactorial causes, affecting 2.6% of all pregnancies worldwide. Nearly half of the RPL cases lack clinically identifiable causes (e.g., antiphospholipid syndrome, uterine anomalies, and parental chromosomal abnormalities), referred to as unexplained RPL (uRPL). Here, we perform a genome-wide association study focusing on uRPL in 1,728 cases and 24,315 female controls of Japanese ancestry. We detect significant associations in the major histocompatibility complex (MHC) region at 6p21 (lead variant=rs9263738; P = 1.4 × 10−10; odds ratio [OR] = 1.51 [95% CI: 1.33–1.72]; risk allele frequency = 0.871). The MHC associations are fine-mapped to the classical HLA alleles, HLA-C*12:02, HLA-B*52:01, and HLA-DRB1*15:02 (P = 1.1 × 10−10, 1.5 × 10−10, and 1.2 × 10−9, respectively), which constitute a population-specific common long-range haplotype with a protective effect (P = 2.8 × 10−10; OR = 0.65 [95% CI: 0.57–0.75]; haplotype frequency=0.108). Genome-wide copy-number variation (CNV) calling demonstrates rare predicted loss-of-function (pLoF) variants of the cadherin-11 gene (CDH11) conferring the risk of uRPL (P = 1.3 × 10−4; OR = 3.29 [95% CI: 1.78–5.76]). Our study highlights the importance of reproductive immunology and rare variants in the uRPL etiology.

3. A second update on mapping the human genetic architecture of COVID-19

   COVID-19 Host Genetics Initiative .

   Nature 621:E7–E26, Sep 2023

   PubMed
4. Genome-wide association study identifies risk loci within the major histocompatibility complex region for Hunner-type interstitial cystitis

   †Akiyama Y., †Sonehara K., Maeda D., Katoh H., Naito T., Yamamoto K., and 7 more authors

   Cell Reports Medicine 4:101114, Jul 2023

   Abstract PubMed

   Hunner-type interstitial cystitis (HIC) is a rare, chronic inflammatory disease of the urinary bladder with unknown etiology and genetic background. Here, we conduct a genome-wide association study of 144 patients with HIC and 41,516 controls of Japanese ancestry. The genetic variant, rs1794275, in the major histocompatibility complex (MHC) region (chromosome 6p21.3) is associated with HIC risk (odds ratio [OR] = 2.32; p = 3.4 × 10⁻⁹). The association is confirmed in a replication set of 26 cases and 1,026 controls (p = 0.014). Fine mapping demonstrates the contribution to the disease risk of a completely linked haplotype of three human leukocyte antigen HLA-DQβ1 amino acid positions, 71, 74, and 75 (OR = 1.94; p = 5 × 10⁻⁸) and of HLA-DPβ1 amino acid position 178, which tags HLA-DPB1∗04:02 (OR = 2.35; p = 7.5 × 10⁻⁸). The three HLA-DQβ1 amino acid positions are located together at the peptide binding groove, suggesting their functional importance in antigen presentation. Our study reveals genetic contributions to HIC risk that may be associated with class II MHC molecule antigen presentation.

5. Reconstruction of the personal information from human genome reads in gut metagenome sequencing data

   Tomofuji Y., Sonehara K., Kishikawa T., Maeda Y., Ogawa K., Kawabata S., and 21 more authors

   Nature Microbiology 8:1079–1094, Jun 2023

   Abstract PubMed

   Human DNA present in faecal samples can result in a small number of human reads in gut shotgun metagenomic sequencing data. However, it is presently unclear how much personal information can be reconstructed from such reads, and this has not been quantitatively evaluated. Such a quantitative evaluation is necessary to clarify the ethical concerns related to data sharing and to enable efficient use of human genetic information in stool samples, such as for research and forensics. Here we used genomic approaches to reconstruct personal information from the faecal metagenomes of 343 Japanese individuals with associated human genotype data. Genetic sex could be accurately predicted based on the sequencing depth of sex chromosomes for 97.3% of the samples. Individuals could be re-identified from the matched genotype data based on human reads recovered from the faecal metagenomic data with 93.3% sensitivity using a likelihood score-based method. This method also enabled us to predict the ancestries of 98.3% of the samples. Finally, we performed ultra-deep shotgun metagenomic sequencing of five faecal samples as well as whole-genome sequencing of blood samples. Using genotype-calling approaches, we demonstrated that the genotypes of both common and rare variants could be reconstructed from faecal samples. This included clinically relevant variants. Our approach can be used to quantify personal information contained within gut metagenome data.

6. DOCK2 is involved in the host genetics and biology of severe COVID-19

   Namkoong H., Edahiro R., Takano T., Nishihara H., Shirai Y., Sonehara K., and 481 more authors

   Nature 609:754–760, Sep 2022

   Abstract PubMed

   Identifying the host genetic factors underlying severe COVID-19 is an emerging challenge1–5. Here we conducted a genome-wide association study (GWAS) involving 2,393 cases of COVID-19 in a cohort of Japanese individuals collected during the initial waves of the pandemic, with 3,289 unaffected controls. We identified a variant on chromosome 5 at 5q35 (rs60200309-A), close to the dedicator of cytokinesis 2 gene (DOCK2), which was associated with severe COVID-19 in patients less than 65 years of age. This risk allele was prevalent in East Asian individuals but rare in Europeans, highlighting the value of genome-wide association studies in non-European populations. RNA-sequencing analysis of 473 bulk peripheral blood samples identified decreased expression of DOCK2 associated with the risk allele in these younger patients. DOCK2 expression was suppressed in patients with severe cases of COVID-19. Single-cell RNA-sequencing analysis (n = 61 individuals) identified cell-type-specific downregulation of DOCK2 and a COVID-19-specific decreasing effect of the risk allele on DOCK2 expression in non-classical monocytes. Immunohistochemistry of lung specimens from patients with severe COVID-19 pneumonia showed suppressed DOCK2 expression. Moreover, inhibition of DOCK2 function with CPYPP increased the severity of pneumonia in a Syrian hamster model of SARS-CoV-2 infection, characterized by weight loss, lung oedema, enhanced viral loads, impaired macrophage recruitment and dysregulated type I interferon responses. We conclude that DOCK2 has an important role in the host immune response to SARS-CoV-2 infection and the development of severe COVID-19, and could be further explored as a potential biomarker and/or therapeutic target.

7. Genetic footprints of assortative mating in the Japanese population

   Yamamoto K., Sonehara K., Namba S., Konuma T., Masuko H., Miyawaki S., and 5 more authors

   Nature Human Behaviour 609:754–760, Sep 2022

   Abstract PubMed

   Assortative mating (AM) is a pattern characterized by phenotypic similarities between mating partners. Detecting the evidence of AM has been challenging due to the lack of large-scale datasets that include phenotypic data on both partners, especially in populations of non-European ancestries. Gametic phase disequilibrium between trait-associated alleles is a signature of parental AM on a polygenic trait, which can be detected even without partner data. Here, using polygenic scores for 81 traits in the Japanese population using BioBank Japan Project genome-wide association studies data (n = 172,270), we found evidence of AM on the liability to type 2 diabetes and coronary artery disease, as well as on dietary habits. In cross-population comparison using United Kingdom Biobank data (n = 337,139) we found shared but heterogeneous impacts of AM between populations.

8. A common deletion at BAK1 reduces enhancer activity and confers risk of intracranial germ cell tumors

   Sonehara K., Kimura Y., Nakano Y., Ozawa T., Takahashi M., Suzuki K., and 22 more authors

   Nature Communications 13:4478, Aug 2022

   Abstract PubMed

   Intracranial germ cell tumors (IGCTs) are rare brain neoplasms that mainly occur in children and adolescents with a particularly high incidence in East Asian populations. Here, we conduct a genome-wide association study (GWAS) of 133 patients with IGCTs and 762 controls of Japanese ancestry. A common 4-bp deletion polymorphism in an enhancer adjacent to BAK1 is significantly associated with the disease risk (rs3831846; P = 2.4 × 10⁻⁹, odds ratio = 2.46 [95% CI: 1.83–3.31], minor allele frequency = 0.43). Rs3831846 is in strong linkage disequilibrium with a testicular GCTs susceptibility variant rs210138. In-vitro reporter assays reveal rs3831846 to be a functional variant attenuating the enhancer activity, suggesting its contribution to IGCTs predisposition through altering BAK1 expression. Risk alleles of testicular GCTs derived from the European GWAS show significant positive correlations in the effect sizes with the Japanese IGCTs GWAS (P = 1.3 × 10⁻⁴, Spearman’s ρ = 0.48). These results suggest the shared genetic susceptibility of GCTs beyond ethnicity and primary sites.

9. Genetic architecture of microRNA expression and its link to complex diseases in the Japanese population

   †Sonehara K., †Sakaue S., Maeda Y., Hirata J., Kishikawa T., Yamamoto K., and 6 more authors

   Human Molecular Genetics 31:1806–1820, Jun 2022

   Abstract PubMed

   Understanding the genetic effects on non-coding RNA (ncRNA) expression facilitates functional characterization of disease-associated genetic loci. Among several classes of ncRNAs, microRNAs (miRNAs) are key post-transcriptional gene regulators. Despite its biological importance, previous studies on the genetic architecture of miRNA expression focused mostly on the European individuals, underrepresented in other populations. Here, we mapped miRNA expression quantitative trait loci (miRNA-eQTL) for 343 miRNAs in 141 Japanese using small RNA sequencing and whole-genome sequencing, identifying 1275 cis-miRNA-eQTL variants for 40 miRNAs (false discovery rate < 0.2). Of these, 25 miRNAs having eQTL were unreported in the European studies, including 5 miRNAs with their lead variant monomorphic in the European populations, which demonstrates the value of miRNA-eQTL analysis in diverse ancestral populations. MiRNAs with eQTL effect showed allele-specific expression (ASE; e.g. miR-146a-3p), and ASE analysis further detected cis-regulatory variants not captured by the conventional miRNA-eQTL mapping (e.g. miR-933). We identified a copy number variation associated with miRNA expression (e.g. miR-570-3p, P = 7.2 ×10^-6), which contributes to a more comprehensive landscape of miRNA-eQTLs. To elucidate a post-transcriptional modification in miRNAs, we created a catalog of miRNA-editing sites, including 10 canonical and 6 non-canonical sites. Finally, by integrating the miRNA-eQTLs and Japanese genome-wide association studies of 25 complex traits (mean n = 192 833), we conducted a transcriptome-wide association study, identifying miR-1908-5p as a potential mediator for adult height, colorectal cancer and type 2 diabetes (P < 9.1 ×10^-5). Our study broadens the population diversity in ncRNA-eQTL studies and contributes to functional annotation of disease-associated loci found in non-European populations.

10. Genome-Wide Association Study of Intracranial Artery Stenosis Followed by Phenome-Wide Association Study

    †Dofuku S., †Sonehara K., Miyawaki S., Sakaue S., Imai H., Shimizu M., and 20 more authors

    Translational Stroke Research 31:1806–1820, Jun 2022

    Abstract PubMed

    The genetic background of intracranial artery stenosis (ICAS), a major cause of ischemic stroke, remains elusive. We performed the world’s first genome-wide association study (GWAS) of ICAS using DNA samples from Japanese subjects, to identify the genetic factors associated with ICAS and their correlation with clinical features. We also conducted a phenome-wide association study (PheWAS) of the top variant identified via GWAS to determine its association with systemic disease. The GWAS involved 408 patients with ICAS and 349 healthy controls and utilized an Asian Screening Array of venous blood samples. The PheWAS was performed using genotypic and phenotypic data of the Biobank Japan Project, which contained information on 46 diseases and 60 quantitative trait data from > 150,000 Japanese individuals. The GWAS revealed that the East Asian-specific functional variant of RNF213, rs112735431 (c.14429G > A, p.Arg4810Lys), was associated with ICAS (odds ratio, 12.3; 95% CI 5.5 to 27.5; P = 7.8 × 10⁻¹⁰). Stratified analysis within ICAS cases demonstrated that clinical features of those with and without the risk allele were different. PheWAS indicated that high blood pressure and angina were significantly associated with RNF213 rs112735431. The first GWAS of ICAS, which stratifies subpopulations within the ICAS cases with distinct clinical features, revealed that RNF213 rs112735431 was the most significant variant associated with ICAS. Thus, RNF213 rs112735431 shows potential as an important clinical biomarker that characterizes pleiotropic risk in various vascular diseases, such as blood pressure and angina, thereby facilitating personalized medicine for systemic vascular diseases in East Asian populations.

11. Mapping the human genetic architecture of COVID-19

    COVID-19 Host Genetics Initiative .

    Nature 600:472–477, Dec 2021

    Abstract PubMed

    The genetic make-up of an individual contributes to the susceptibility and response to viral infection. Although environmental, clinical and social factors have a role in the chance of exposure to SARS-CoV-2 and the severity of COVID-191,2, host genetics may also be important. Identifying host-specific genetic factors may reveal biological mechanisms of therapeutic relevance and clarify causal relationships of modifiable environmental risk factors for SARS-CoV-2 infection and outcomes. We formed a global network of researchers to investigate the role of human genetics in SARS-CoV-2 infection and COVID-19 severity. Here we describe the results of three genome-wide association meta-analyses that consist of up to 49,562 patients with COVID-19 from 46 studies across 19 countries. We report 13 genome-wide significant loci that are associated with SARS-CoV-2 infection or severe manifestations of COVID-19. Several of these loci correspond to previously documented associations to lung or autoimmune and inflammatory diseases3-7. They also represent potentially actionable mechanisms in response to infection. Mendelian randomization analyses support a causal role for smoking and body-mass index for severe COVID-19 although not for type II diabetes. The identification of novel host genetic factors associated with COVID-19 was made possible by the community of human genetics researchers coming together to prioritize the sharing of data, results, resources and analytical frameworks. This working model of international collaboration underscores what is possible for future genetic discoveries in emerging pandemics, or indeed for any complex human disease.

12. Obelisc: an identical-by-descent mapping tool based on SNP streak

    Sonehara K., and Okada Y.

    Bioinformatics (Oxford, England) 36:5567–5570, Apr 2021

    Abstract PubMed

    MOTIVATION: Genetic linkage analysis has made a huge contribution to the genetic mapping of Mendelian diseases. However, most previously available linkage analysis methods have limited applicability. Since parametric linkage analysis requires predefined model of inheritance with a fixed set of parameters, it is inapplicable without fully structured pedigree information. Furthermore, the analytical results are dependent on the specification of model parameters. While non-parametric linkage analysis can avoid these problems, the runs of homozygosity (ROH) mapping, a widely used non-parametric linkage analysis method, can only deal with recessive inheritance. The implementation of non-parametric linkage analyses capable of dealing with both dominant and recessive inheritance has been required.
    RESULTS: We have developed the Obelisc (Observational linkage scan), a flexibly applicable user-friendly non-parametric linkage analysis tool, which also provides an intuitive visualization of the analytical results. Obelisc is based on the SNP streak approach, which does not require any predefined inheritance model with parameters. In contrast to the ROH mapping, the SNP streak approach is applicable to both dominant and recessive traits. To illustrate the performance of Obelisc, we generated a pseudo-pedigree from the publicly available BioBank Japan Project genome-wide genotype dataset (n > 180 000). By applying Obelisc to this pseudo-pedigree, we successfully identified the regions with inherited identical-by-descent haplotypes shared among the members of the pseudo-pedigree, which was validated by the population-based haplotype phasing approach.
    AVAILABILITY AND IMPLEMENTATION: Obelisc is feely available at https://github.com/qsonehara/Obelisc as a python package with example datasets. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.