Genotype imputation methods for whole and complex genomic regions utilizing deep learning technology

Naito T., Okada Y.. Genotype imputation methods for whole and complex genomic regions utilizing deep learning technology. Journal of Human Genetics , (2024); https://doi.org/10.1038/s10038-023-01213-6.The imputation of unmeasured genotypes is essential in human genetic research, particularly in enhancing the power of genome-wide association studies and conducting subsequent fine-mapping. Recently, several deep learning-based genotype imputation methods for genome-wide variants with the capability of learning complex linkage disequilibrium patterns have been developed. Additionally, deep learning-based imputation has been applied to a distinct genomic region known as the major histocompatibility complex, referred to as HLA imputation. Despite their various advantages, the current deep learning-based genotype imputation methods do have certain limitations and have not yet become standard. These limitations include the modest accuracy improvement over statistical and conventional machine learning-based methods. However, their benefits include other aspects, such as their “reference-free” nature, which ensures complete privacy protection, and their higher computational efficiency. Furthermore, the continuing evolution of deep learning technologies is expected to contribute to further improvements in prediction accuracy and usability in the future

Genome Wide Association Study of Age at Menarche in the Japanese Population

Age at menarche (AAM) is a complex trait involving both genetic and environmental factors. To identify the genetic factors associated with AAM, we conducted a large-scale meta-analysis of genome-wide association studies using more than 15,000 Japanese female samples. Here, we identified an association between SNP (single nucleotide polymorphism) rs364663 at the LIN28B locus and AAM, with a P-value of 5.49×10−7 and an effect size of 0.089 (year). We also evaluated 33 SNPs that were previously reported to be associated with AAM in women of European ancestry. Among them, two SNPs rs4452860 and rs7028916 in TMEM38B indicated significant association with AAM in the same directions as reported in previous studies (P = 0.0013 with an effect size of 0.051) even after Bonferroni correction for the 33 SNPs. In addition, six loci in or near CCDC85A, LOC100421670, CA10, ZNF483, ARNTL, and RXRG exhibited suggestive association with AAM (P<0.05). Our findings elucidated the impact of genetic variations on AAM in the Japanese population

Genome Wide Association Study of Age at Menarche in the Japanese Population

Age at menarche (AAM) is a complex trait involving both genetic and environmental factors. To identify the genetic factors associated with AAM, we conducted a large-scale meta-analysis of genome-wide association studies using more than 15,000 Japanese female samples. Here, we identified an association between SNP (single nucleotide polymorphism) rs364663 at the LIN28B locus and AAM, with a P-value of 5.49×10−7 and an effect size of 0.089 (year). We also evaluated 33 SNPs that were previously reported to be associated with AAM in women of European ancestry. Among them, two SNPs rs4452860 and rs7028916 in TMEM38B indicated significant association with AAM in the same directions as reported in previous studies (P = 0.0013 with an effect size of 0.051) even after Bonferroni correction for the 33 SNPs. In addition, six loci in or near CCDC85A, LOC100421670, CA10, ZNF483, ARNTL, and RXRG exhibited suggestive association with AAM (P<0.05). Our findings elucidated the impact of genetic variations on AAM in the Japanese population.</p

Loci associated with N-glycosylation of human IgG are not associated with rheumatoid arthritis: a Mendelian randomisation study

Objectives: A recent study identified 16 genetic variants associated with N-glycosylation of human IgG. Several of the genomic regions where these single nucleotide polymorphisms (SNPs) reside have also been associated with autoimmune disease (AID) susceptibility, suggesting there may be pleiotropy (genetic sharing) between loci controlling both N-glycosylation and AIDs. We investigated this by testing variants associated with levels of IgG N-glycosylation for association with rheumatoid arthritis (RA) susceptibility using a Mendelian randomisation study, and testing a subset of these variants in a less well-powered study of treatment response and severity. Methods: SNPs showing association with IgG N-glycosylation were analysed for association with RA susceptibility in 14 361 RA cases and 43 923 controls. Five SNPs were tested for association with response to anti-tumour necrosis factor (TNF) therapy in 1081 RA patient samples and for association with radiological disease severity in 342 patients. Results: Only one SNP (rs9296009) associated with N-glycosylation showed an association (p=6.92×10–266) with RA susceptibility, although this was due to linkage disequilibrium with causal human leukocyte antigen (HLA) variants. Four regions of the genome harboured SNPs associated with both traits (shared loci); although statistical analysis indicated that the associations observed for the two traits are independent. No SNPs showed association with response to anti-TNF therapy. One SNP rs12342831 was modestly associated with Larsen score (p=0.05). Conclusions: In a large, well-powered cohort of RA patients, we show SNPs driving levels of N-glycosylation have no association with RA susceptibility, indicating colocalisation of associated SNPs are not necessarily indicative of a shared genetic background or a role for glycosylation in disease susceptibility

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

Sonehara K., Yano Y., Naito T., et al. Common and rare genetic variants predisposing females to unexplained recurrent pregnancy loss. Nature communications 15, 5744 (2024); https://doi.org/10.1038/s41467-024-49993-5.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

Variants at HLA-A , HLA-C , and HLA-DQB1 Confer Risk of Psoriasis Vulgaris in Japanese

Psoriasis vulgaris (PsV) is an autoimmune disease of skin and joints with heterogeneity in epidemiologic and genetic landscapes of global populations. We conducted an initial genome-wide association study and a replication study of PsV in the Japanese population (606 PsV cases and 2,052 controls). We identified significant associations of the single nucleotide polymorphisms with PsV risk at TNFAIP3-interacting protein 1and the major histocompatibility complex region (P = 3.7 × 10−10 and 6.6 × 10−15, respectively). By updating the HLA imputation reference panel of Japanese (n = 908) to expand HLA gene coverage, we fine-mapped the HLA variants associated with PsV risk. Although we confirmed the PsV risk of HLA-C*06:02 (odds ratio = 6.36, P = 0.0015), its impact was relatively small compared with those in other populations due to rare allele frequency in Japanese (0.4% in controls). Alternatively, HLA-A*02:07, which corresponds to the cysteine residue at HLA-A amino acid position 99 (HLA-A Cys99), demonstrated the most significant association with PsV (odds ratio = 4.61, P = 1.2 × 10–10). In addition to HLA-A*02:07 and HLA-C*06:02, stepwise conditional analysis identified an independent PsV risk of HLA-DQβ1 Asp57 (odds ratio = 2.19, P = 1.9 × 10–6). Our PsV genome-wide association study in Japanese highlighted the genetic architecture of PsV, including the identification of HLA risk variants

Multi-ancestry meta-analysis of asthma identifies novel associations and highlights the value of increased power and diversity

Publisher Copyright: © 2022 The Author(s)Asthma is a complex disease that varies widely in prevalence across populations. The extent to which genetic variation contributes to these disparities is unclear, as the genetics underlying asthma have been investigated primarily in populations of European descent. As part of the Global Biobank Meta-analysis Initiative, we conducted a large-scale genome-wide association study of asthma (153,763 cases and 1,647,022 controls) via meta-analysis across 22 biobanks spanning multiple ancestries. We discovered 179 asthma-associated loci, 49 of which were not previously reported. Despite the wide range in asthma prevalence among biobanks, we found largely consistent genetic effects across biobanks and ancestries. The meta-analysis also improved polygenic risk prediction in non-European populations compared with previous studies. Additionally, we found considerable genetic overlap between age-of-onset subtypes and between asthma and comorbid diseases. Our work underscores the multi-factorial nature of asthma development and offers insight into its shared genetic architecture.Peer reviewe

Genome-wide association study on meningioma risk in Japan: a multicenter prospective study

The version of record of this article, first published in Journal of Neuro-Oncology, is available online at Publisher’s website: https://doi.org/10.1007/s11060-024-04727-x.Purpose: Although meningiomas are the most common primary intracranial tumors, their genetic etiologies have not been fully elucidated. To date, only two genome-wide association studies (GWASs) have focused on European ancestries, despite ethnic differences in the incidence of meningiomas. The aim of this study was to conduct the first GWAS of Japanese patients with meningiomas to identify the SNPs associated with meningioma susceptibility. Methods: In this multicenter prospective case-control study, we studied 401 Japanese patients with meningioma admitted in five institutions in Japan, and 50,876 control participants of Japanese ancestry enrolled in Biobank Japan. Results: The quality control process yielded 536,319 variants and imputation resulted in 8,224,735 variants on the autosomes and 224,820 variants on the X chromosomes. This GWAS eventually revealed no genetic variants with genome-wide significance (P < 5 × 10 − 8) and observed no significant association in the previously reported risk variants rs11012732 and rs2686876 due to low minor allele frequency in the Japanese population. Conclusion: This is the first GWAS of meningiomas in East Asian populations and is expected to contribute to the development of GWAS research for meningiomas

A single centre retrospective analysis of AECG classification criteria for primary Sjogren\u27s syndrome based on 112 minor salivary gland biopsies in a Japanese population.

OBJECTIVE: To assess the usefulness and performance of the American European Consensus Group (AECG) criteria based on minor salivary gland biopsy (MSGB) in Japanese patients with primary SS. METHODS: Among 208 MSGB cases, we retrospectively selected 112 subjects who satisfied the complete set of AECG classification criteria. Of the 112 subjects studied, 63 primary SS patients and 49 non-SS group subjects were classified according to the AECG criteria. The contribution of subjective and objective components was statistically analysed. RESULTS: Sex, dry eye, Saxon test, Schirmer\u27s test, anti-SSA/Ro antibody, MSGB grading and sialography statistically contributed to the diagnosis. Multiple logistic regression analysis showed that positive MSGB [odds ratio (OR) 105; 95% CI 13, 849), positive anti-SSA/Ro antibody (OR 96; 95% CI 10, 923), a positive Saxon test (OR 46; 95% CI, 6, 340) and the existence of dry eye (OR 8, 95% CI 2, 43) were associated with the diagnosis of primary SS. Among the components of the AECG criteria, MSGB and anti-SSA/Ro antibody were very strong contributors. Furthermore, the abnormal-finding positive rate in sialography significantly correlated with MSGB grading (P-value for trend = 0.0006), although other subjective and objective components were not associated with MSGB grading. CONCLUSION: The usefulness of the AECG criteria for Japanese primary SS patients was confirmed

Association of Common Variants in TNFRSF13B, TNFSF13, and ANXA3 with Serum Levels of Non-Albumin Protein and Immunoglobulin Isotypes in Japanese

We performed a genome-wide association study (GWAS) on levels of serum total protein (TP), albumin (ALB), and non-albumin protein (NAP). We analyzed SNPs on autosomal chromosomes using data from 9,103 Japanese individuals, followed by a replication study of 1,600 additional individuals. We confirmed the previously- reported association of GCKR on chromosome 2p23.3 with serum ALB (rs1260326, Pmeta = 3.1×10−9), and additionally identified the significant genome-wide association of rs4985726 in TNFRSF13B on 17p11.2 with both TP and NAP (Pmeta = 1.2×10−14 and 7.1×10−24, respectively). For NAP, rs3803800 and rs11552708 in TNFSF13 on 17p13.1 (Pmeta = 7.2×10−15 and 7.5×10−10, respectively) as well as rs10007186 on 4q21.2 near ANXA3 (Pmeta = 1.3×10−9) also indicated significant associations. Interestingly, TNFRSF13B and TNFSF13 encode a tumor necrosis factor (TNF) receptor and its ligand, which together constitute an important receptor-ligand axis for B-cell homeostasis and immunoglobulin production. Furthermore, three SNPs, rs4985726, rs3803800, and rs11552708 in TNFRSF13B and TNFSF13, were indicated to be associated with serum levels of IgG (P<2.3×10−3) and IgM (P<0.018), while rs3803800 and rs11552708 were associated with IgA (P<0.013). Rs10007186 in 4q21.2 was associated with serum levels of IgA (P = 0.036), IgM (P = 0.019), and IgE (P = 4.9×10−4). Our results should add interesting knowledge about the regulation of major serum components