Kyuto Sonehara, MD, PhD
曽根原 究人

†: equal contribution

1. 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

   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.

2. 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.

3. 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.

4. 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 8:1079–1094, 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.

5. 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.

6. 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.

7. 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.

8. 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.