jMorp: Japanese Multi-Omics Reference Panel update report 2023

Modern medicine is increasingly focused on personalized medicine, and multi-omics data is crucial in understanding biological phenomena and disease mechanisms. Each ethnic group has its unique genetic background with specific genomic variations influencing disease risk and drug response. Therefore, multi-omics data from specific ethnic populations are essential for the effective implementation of personalized medicine. Various prospective cohort studies, such as the UK Biobank, All of Us and Lifelines, have been conducted worldwide. The Tohoku Medical Megabank project was initiated after the Great East Japan Earthquake in 2011. It collects biological specimens and conducts genome and omics analyses to build a basis for personalized medicine. Summary statistical data from these analyses are available in the jMorp web database (https://jmorp.megabank.tohoku.ac.jp), which provides a multidimensional approach to the diversity of the Japanese population. jMorp was launched in 2015 as a public database for plasma metabolome and proteome analyses and has been continuously updated. The current update will significantly expand the scale of the data (metabolome, genome, transcriptome, and metagenome). In addition, the user interface and backend server implementations were rewritten to improve the connectivity between the items stored in jMorp. This paper provides an overview of the new version of the jMorp.

Tests of Light-Lepton Universality in Angular Asymmetries of B^{0}→D^{*-}ℓν Decays

We present the first comprehensive tests of the universality of the light leptons in the angular distributions of semileptonic B^{0}-meson decays to charged spin-1 charmed mesons. We measure five angular-asymmetry observables as functions of the decay recoil that are sensitive to lepton-universality-violating contributions. We use events where one neutral B is fully reconstructed in ϒ(4S)→BB[over ¯] decays in data corresponding to 189  fb^{-1} integrated luminosity from electron-positron collisions collected with the Belle II detector. We find no significant deviation from the standard model expectations.

Precise Measurement of the D_{s}^{+} Lifetime at Belle II

We measure the lifetime of the D_{s}^{+} meson using a data sample of 207  fb^{-1} collected by the Belle II experiment running at the SuperKEKB asymmetric-energy e^{+}e^{-} collider. The lifetime is determined by fitting the decay-time distribution of a sample of 116×10^{3} D_{s}^{+}→ϕπ^{+} decays. Our result is τ_{D_{s}^{+}}=(499.5±1.7±0.9)  fs, where the first uncertainty is statistical and the second is systematic. This result is significantly more precise than previous measurements.

Search for a τ^{+}τ^{-} Resonance in e^{+}e^{-}→μ^{+}μ^{-}τ^{+}τ^{-} Events with the Belle II Experiment

We report the first search for a nonstandard-model resonance decaying into τ pairs in e^{+}e^{-}→μ^{+}μ^{-}τ^{+}τ^{-} events in the 3.6-10  GeV/c^{2} mass range. We use a 62.8  fb^{-1} sample of e^{+}e^{-} collisions collected at a center-of-mass energy of 10.58 GeV by the Belle II experiment at the SuperKEKB collider. The analysis probes three different models predicting a spin-1 particle coupling only to the heavier lepton families, a Higgs-like spin-0 particle that couples preferentially to charged leptons (leptophilic scalar), and an axionlike particle, respectively. We observe no evidence for a signal and set exclusion limits at 90% confidence level on the product of cross section and branching fraction into τ pairs, ranging from 0.7 to 24 fb, and on the couplings of these processes. We obtain world-leading constraints on the couplings for the leptophilic scalar model for masses above 6.5  GeV/c^{2} and for the axionlike particle model over the entire mass range.

Measurement of CP Violation in B^{0}→K_{S}^{0}π^{0} Decays at Belle II

We report a measurement of the CP-violating parameters C and S in B^{0}→K_{S}^{0}π^{0} decays at Belle II using a sample of 387×10^{6}  BB[over ¯] events recorded in e^{+}e^{-} collisions at a center-of-mass energy corresponding to the ϒ(4S) resonance. These parameters are determined by fitting the proper decay-time distribution of a sample of 415 signal events. We obtain C=-0.04_{-0.15}^{+0.14}±0.05 and S=0.75_{-0.23}^{+0.20}±0.04, where the first uncertainties are statistical and the second are systematic.

Test of Light-Lepton Universality in the Rates of Inclusive Semileptonic B-Meson Decays at Belle II

We present the first measurement of the ratio of branching fractions of inclusive semileptonic B-meson decays, R(X_{e/μ})=B(B→Xeν)/B(B→Xμν), a precision test of electron-muon universality, using data corresponding to 189  fb^{-1} from electron-positron collisions collected with the Belle II detector. In events where the partner B meson is fully reconstructed, we use fits to the lepton momentum spectra above 1.3  GeV/c to obtain R(X_{e/μ})=1.007±0.009(stat)±0.019(syst), which is the most precise lepton-universality test of its kind and agrees with the standard-model expectation.

Search for an Invisible Z^{'} in a Final State with Two Muons and Missing Energy at Belle II

The L_{μ}-L_{τ} extension of the standard model predicts the existence of a lepton-flavor-universality-violating Z^{'} boson that couples only to the heavier lepton families. We search for such a Z^{'} through its invisible decay in the process e^{+}e^{-}→μ^{+}μ^{-}Z^{'}. We use a sample of electron-positron collisions at a center-of-mass energy of 10.58 GeV collected by the Belle II experiment in 2019-2020, corresponding to an integrated luminosity of 79.7  fb^{-1}. We find no excess over the expected standard-model background. We set 90%-confidence-level upper limits on the cross section for this process as well as on the coupling of the model, which ranges from 3×10^{-3} at low Z^{'} masses to 1 at Z^{'} masses of 8  GeV/c^{2}.

Prognostic value of increased intraepithelial lymphocytes and lymphocytic clonality in dogs with chronic enteropathy or small-cell lymphoma

The clinical significance of severe infiltration of small intraepithelial lymphocytes (IEL) and the results of polymerase chain reaction for antigen receptor rearrangement (PARR) in dogs with chronic enteropathy (CE) and small-cell lymphoma (SCL) are controversial. This cohort study aimed to evaluate the prognostic significance of the IEL and PARR results in dogs with CE or SCL. Although definitive diagnostic histopathological criteria for SCL in dogs have yet to be established, dogs with the histopathological findings of severe IEL infiltration were diagnosed with SCL in this study. One hundred and nineteen dogs were recruited, with 23 dogs classified as having SCL and 96 dogs as having CE. The positive rate of PARR was 59.6% (71/119) in the duodenum and 57.7% (64/111) in the ileum. Subsequently, three dogs with SCL and four dogs with CE developed large-cell lymphoma (LCL). The median overall survival (OS) of dogs with SCL was 700 days (range, 6-1410 days), and that of dogs with CE was not reached. In the log-rank test, shorter OS was observed in cases with histopathological SCL (P = 0.035), clonal TCRγ rearrangement in the duodenum (P < 0.012), and clonal IgH rearrangement in the ileum (P < 0.0001). The Cox proportional hazards model adjusted for sex and age showed that histopathological SCL (hazard ratio [HR] 1.74; 95% confidence interval [CI], 0.83-3.65), duodenal clonal TCRγ rearrangement (HR, 1.80; 95% CI, 0.86-3.75), and ileal clonal IgH rearrangement (HR, 2.28; 95% CI, 0.92-5.70) could shorten overall survival, although their 95% CIs included 1.0. These results indicate that severe IEL infiltration could be a useful histopathological feature for diagnosing SCL, and clonality-positive results could be a negative prognostic factor in dogs with CE. Furthermore, the development of LCL should be carefully monitored in dogs with CE and SCL.

Search for Lepton-Flavor-Violating τ Decays to a Lepton and an Invisible Boson at Belle II

We search for lepton-flavor-violating τ^{-}→e^{-}α and τ^{-}→μ^{-}α decays, where α is an invisible spin-0 boson. The search uses electron-positron collisions at 10.58 GeV center-of-mass energy with an integrated luminosity of 62.8  fb^{-1}, produced by the SuperKEKB collider and collected with the Belle II detector. We search for an excess in the lepton-energy spectrum of the known τ^{-}→e^{-}ν[over ¯]_{e}ν_{τ} and τ^{-}→μ^{-}ν[over ¯]_{μ}ν_{τ} decays. We report 95% confidence-level upper limits on the branching-fraction ratio B(τ^{-}→e^{-}α)/B(τ^{-}→e^{-}ν[over ¯]_{e}ν_{τ}) in the range (1.1-9.7)×10^{-3} and on B(τ^{-}→μ^{-}α)/B(τ^{-}→μ^{-}ν[over ¯]_{μ}ν_{τ}) in the range (0.7-12.2)×10^{-3} for α masses between 0 and 1.6  GeV/c^{2}. These results provide the most stringent bounds on invisible boson production from τ decays.

Observation of e^{+}e^{-}→ωχ_{bJ}(1P) and Search for X_{b}→ωϒ(1S) at sqrt[s] near 10.75 GeV

We study the processes e^{+}e^{-}→ωχ_{bJ}(1P) (J=0, 1, or 2) using samples at center-of-mass energies sqrt[s]=10.701, 10.745, and 10.805 GeV, corresponding to 1.6, 9.8, and 4.7  fb^{-1} of integrated luminosity, respectively. These data were collected with the Belle II detector during special operations of the SuperKEKB collider above the ϒ(4S) resonance. We report the first observation of ωχ_{bJ}(1P) signals at sqrt[s]=10.745  GeV. By combining Belle II data with Belle results at sqrt[s]=10.867  GeV, we find energy dependencies of the Born cross sections for e^{+}e^{-}→ωχ_{b1,b2}(1P) to be consistent with the shape of the ϒ(10753) state. These data indicate that the internal structures of the ϒ(10753) and ϒ(10860) states may differ. Including data at sqrt[s]=10.653  GeV, we also search for the bottomonium equivalent of the X(3872) state decaying into ωϒ(1S). No significant signal is observed for masses between 10.45 and 10.65  GeV/c^{2}.

Measurement of the Λ_{c}^{+} Lifetime

An absolute measurement of the Λ_{c}^{+} lifetime is reported using Λ_{c}^{+}→pK^{-}π^{+} decays in events reconstructed from data collected by the Belle II experiment at the SuperKEKB asymmetric-energy electron-positron collider. The total integrated luminosity of the data sample, which was collected at center-of-mass energies at or near the ϒ(4S) resonance, is 207.2  fb^{-1}. The result, τ(Λ_{c}^{+})=203.20±0.89±0.77  fs, where the first uncertainty is statistical and the second systematic, is the most precise measurement to date and is consistent with previous determinations.

Search for a Dark Photon and an Invisible Dark Higgs Boson in μ^{+}μ^{-} and Missing Energy Final States with the Belle II Experiment

The dark photon A^{'} and the dark Higgs boson h^{'} are hypothetical particles predicted in many dark sector models. We search for the simultaneous production of A^{'} and h^{'} in the dark Higgsstrahlung process e^{+}e^{-}→A^{'}h^{'} with A^{'}→μ^{+}μ^{-} and h^{'} invisible in electron-positron collisions at a center-of-mass energy of 10.58 GeV in data collected by the Belle II experiment in 2019. With an integrated luminosity of 8.34  fb^{-1}, we observe no evidence for signal. We obtain exclusion limits at 90% Bayesian credibility in the range of 1.7-5.0 fb on the cross section and in the range of 1.7×10^{-8}-200×10^{-8} on the effective coupling ϵ^{2}×α_{D} for the A^{'} mass in the range of 4.0  GeV/c^{2}<M_{A^{'}}<9.7  GeV/c^{2} and for the h^{'} mass M_{h^{'}}<M_{A^{'}}, where ϵ is the mixing strength between the standard model and the dark photon and α_{D} is the coupling of the dark photon to the dark Higgs boson. Our limits are the first in this mass range.

Study on ammonia generation from digested sludge by subcritical water treatment

Organic sludge has recently attracted attention as a renewable energy source. While the organic matter contained in sludge can be utilized as various renewable energy sources, its nitrogen component has limited use. In this study, subcritical water treatment was conducted to recover ammonia from digested sludge. While ammonia recovered via stripping is limited to soluble components, subcritical water treatment can convert solid components and dissolved organic nitrogen sludge into ammonia. Digested sludge was treated at several reaction temperatures, reaction pressures, treatment times, and oxygen ratios to determine the ammonia generation rate. Among the conditions tested in this study, an ammonium generation rate of 84.0% was obtained at 400 °C, 10 MPa, a treatment time of 5 min, and at an oxygen ratio of 1.2.

Sex-Specific Differences in the Transcriptome of the Human Dorsolateral Prefrontal Cortex in Schizophrenia

Schizophrenia presents clinical and biological differences between males and females. This study investigated transcriptional profiles in the dorsolateral prefrontal cortex (DLPFC) using postmortem data from the largest RNA-sequencing (RNA-seq) database on schizophrenic cases and controls. Data for 154 male and 113 female controls and 160 male and 93 female schizophrenic cases were obtained from the CommonMind Consortium. In the RNA-seq database, the principal component analysis showed that sex effects were small in schizophrenia. After we analyzed the impact of sex-specific differences on gene expression, the female group showed more significantly changed genes compared with the male group. Based on the gene ontology analysis, the female sex-specific genes that changed were overrepresented in the mitochondrion, ATP (phosphocreatine and adenosine triphosphate)-, and metal ion-binding relevant biological processes. An ingenuity pathway analysis revealed that the differentially expressed genes related to schizophrenia in the female group were involved in midbrain dopaminergic and γ-aminobutyric acid (GABA)-ergic neurons and microglia. We used methylated DNA-binding domain-sequencing analyses and microarray to investigate the DNA methylation that potentially impacts the sex differences in gene transcription using a maternal immune activation (MIA) murine model. Among the sex-specific positional genes related to schizophrenia in the PFC of female offspring from MIA, the changes in the methylation and transcriptional expression of loci ACSBG1 were validated in the females with schizophrenia in independent postmortem samples by real-time PCR and pyrosequencing. Our results reveal potential genetic risks in the DLPFC for the sex-dependent prevalence and symptomology of schizophrenia.

Convolutional neural network-based skin image segmentation model to improve classification of skin diseases in conventional and non-standardized picture images

BACKGROUND

For dermatological practices, non-standardized conventional photo images are taken and collected as a mixture of variable fields of the image view, including close-up images focusing on designated lesions and long-shot images including normal skin and background of the body surface. Computer-aided detection/diagnosis (CAD) models trained using non-standardized conventional photo images exhibit lower performance rates than CAD models that detect lesions in a localized small area, such as dermoscopic images.

OBJECTIVE

We aimed to develop a convolutional neural network (CNN) model for skin image segmentation to generate a skin disease image dataset suitable for CAD of multiple skin disease classification.

METHODS

We trained a DeepLabv3 + -based CNN segmentation model to detect skin and lesion areas and segmented out areas that satisfy the following conditions: more than 80% of the image will be the skin area, and more than 10% of the image will be the lesion area.

RESULTS

The generated CNN-segmented image database was examined using CAD of skin disease classification and achieved approximately 90% sensitivity and specificity to differentiate atopic dermatitis from malignant diseases and complications, such as mycosis fungoides, impetigo, and herpesvirus infection. The accuracy of skin disease classification in the CNN-segmented image dataset was almost equal to that of the manually cropped image dataset and higher than that of the original image dataset.

CONCLUSION

Our CNN segmentation model, which automatically extracts lesions and segmented images of the skin area regardless of image fields, will reduce the burden of physician annotation and improve CAD performance.

Genome-wide Association Study of Axial Length in Population-based Cohorts in Japan

Purpose

To elucidate the differences in ocular biometric parameters by generation and gender and to identify axial length (AL)-associated genetic variants in Japanese individuals, we analyzed Tohoku Medical Megabank Organization (ToMMo) Eye Study data.

Design

We designed the ToMMo Eye Study, examined AL variations, and conducted genome-wide association studies (GWASs).

Participants

In total, 33 483 participants aged > 18 years who were recruited into the community-based cohort (CommCohort) and the birth and three-generation cohort (BirThree Cohort) of the ToMMo Eye Study were examined.

Methods

Each participant was screened with an interview, ophthalmic examinations, and a microarray analysis. The GWASs were performed in 22 379 participants in the CommCohort (discovery stage) and 11 104 participants in the BirThree Cohort (replication stage). We evaluated the associations of single nucleotide polymorphisms (SNPs) with AL using a genome-wide significance threshold (5 × 10^-8) in each stage of the study and in the subsequent meta-analysis.

Main Outcome Measures

We identified the association of SNPs with AL and distributions of AL in right and left eyes and individuals of different sexes and ages.

Results

In the discovery stage, the mean AL of the right eye (23.99 mm) was significantly greater than that of the left eye (23.95 mm). This difference was reproducible across sexes and ages. The GWASs revealed 703 and 215 AL-associated SNPs with genome-wide significance in the discovery and validation stages, respectively, and many of the SNPs in the discovery stage were replicated in the validation stage. Validated SNPs and their associated loci were meta-analyzed for statistical significance (P < 5 × 10^-8). This study identified 1478 SNPs spread over 31 loci. Of the 31 loci, 5 are known AL loci, 15 are known refractive-error loci, 4 are known corneal-curvature loci, and 7 loci are newly identified loci that are not known to be associated with AL. Of note, some of them shared functional relationships with previously identified loci.

Conclusions

Our large-scale GWASs exploiting ToMMo Eye Study data identified 31 loci linked to variations in AL, 7 of which are newly reported in this article. The results revealed genetic heterogeneity and similarity in SNPs related to ethnic variations in AL.

Genetic loci for lung function in Japanese adults with adjustment for exhaled nitric oxide levels as airway inflammation indicator

Lung function reflects the ability of the respiratory system and is utilized for the assessment of respiratory diseases. Because type 2 airway inflammation influences lung function, genome wide association studies (GWAS) for lung function would be improved by adjustment with an indicator of the inflammation. Here, we performed a GWAS for lung function with adjustment for exhaled nitric oxide (FeNO) levels in two independent Japanese populations. Our GWAS with genotype imputations revealed that the RNF5/AGER locus including AGER rs2070600 SNP, which introduces a G82S substitution of AGER, was the most significantly associated with FEV1/FVC. Three other rare missense variants of AGER were further identified. We also found genetic loci with three candidate genes (NOS2, SPSB2 and RIPOR2) associated with FeNO levels. Analyses with the BioBank-Japan GWAS resource revealed genetic links of FeNO and asthma-related traits, and existence of common genetic background for allergic diseases and their biomarkers. Our study identified the genetic locus most strongly associated with airway obstruction in the Japanese population and three genetic loci associated with FeNO, an indicator of type 2 airway inflammation in adults.

Ehlers-Danlos syndrome type IV with a novel COL3A1 exon 14 skipping variation confirmed by Tohoku Medical Megabank Organization genomic database

A novel COL3A1 variant was identified in a Japanese case of Ehlers-Danlos syndrome type IV (EDS-IV) with a characteristic "Madonna" face, fragile uterus, and easy bruising in addition to a history of cavernous sinus fistula. We confirmed variable diameters of collagen fibrils in the dermis and decrease in type 3 collagen production from cultured fibroblasts. Genomic DNA sequencing of the COL3A1 region and COL3A1 cDNA sequence expressing in cultured fibroblasts identified that a nucleotide variation at c.951+2T>G on intron 14 leads to skipping of exon 14 in COL3A1 cDNA. The novel variation in the splice site of COL3A1 region g.IVS14+2T>G was not listed in the EDS-IV pathogenic genetic databases including Human Gene Mutation Database, ClinVar, and Leiden Open Variation Database. Using the whole genome sequence database of 8380 Japanese individuals reported by the Tohoku Medical Megabank Organization (ToMMo) cohort study, we also confirmed that COL3A1 g.IVS14+2T>G was not a common single nucleotide variation in the Japanese population, although 13 EDS-related COL3A1 variants were identified in the ToMMo database of 8380 Japanese individuals. These results demonstrated that our case of EDS-IV was a result of the novel variation of COL3A1 g.IVS14+2T>G. These statistical genetics approaches with the combination of the ToMMo database of 8380 Japanese individuals and pathogenic genetic databases are a useful method to confirm the uniqueness of novel variation in Japanese.

P–145 usefulness of morphokinetic data to predict pregnancy rates of day–6 blastocyst transfers

Study question

Can a scoring model based on morphokinetic data developed to predict pregnancy rates of day–5 blastocyst transfers (KIDSCORE™D5) predict pregnancy rates of day–6 blastocyst transfers?

Summary answer

KIDSCORE™D5 was able to predict the clinical pregnancy rates of embryo transfers done on day 6 with an area under the curve (AUC) of 0.72.

What is known already

KIDSCORE™D5 is a scoring model based on morphokinetic data developed to predict the pregnancy rates of day–5 blastocysts. In 2019, Regnier et al. reported that the AUC of KIDSCORE™D5 for predicting clinical pregnancy rates of day–5 blastocyst transfers was 0.6. However, as KIDSCORE™D5 is constructed based on morphological characteristics and developmental dynamics of day–5 blastocysts, it is unclear whether KIDSCORE™D5 can predict pregnancy rates of day–6 blastocyst transfers. Since there are many cases of day–6 blastocyst transfers, it is important to know if KIDSCORE™D5 can predict pregnancy rates of day–6 blastocyst transfers.

Study design, size, duration

This retrospective single-center study, which included 162 day–5 and 72 day–6 blastocyst transfers, respectively, was conducted at Takahashi Women’s clinic from January to December 2019. Blastocysts derived from 146 patients who underwent intracytoplasmic sperm injection. All blastocysts were cryopreserved and were transferred singly.

Participants/materials, setting, methods

We used EmbryoScope+™ (Vitrolife) for in-vitro culture and calculated KIDSCORE™D5 (ver.3) using Embryoviewer™ (Vitrolife). Blastocyst scoring was done from 1.0 to 9.9. Clinical pregnancy was defined as the presence of a gestational sac confirmed by transvaginal ultrasonography. Statistical analysis was performed with JMP Pro 15.00 (SAS). The relationship between KIDSCORE™D5 and clinical pregnancy was evaluated by the AUC using ROC curve analysis and multivariate analysis adjusted for patient age.

Main results and the role of chance

The mean KIDSCORE™D5 of day–5 and day–6 blastocysts was 7.1±1.7 and 3.7±1.5, respectively. KIDSCORE™D5 of day–6 blastocysts was significantly lower than that of day–5 blastocysts (p &lt; 0.0001, Wilcoxon test). ROC curve analysis showed that the KIDSCORE™D5 could predict clinical pregnancy rates with an AUC of 0.62 for day–5 blastocysts and 0.72 for day–6 blastocysts. The cut-off values for KIDSCORE™D5 were 5.7 and 4.9 for day–5 and day–6 blastocysts, respectively. Blastocysts above the cut off value on both day–5 and day–6 had a significantly higher pregnancy rate than those below the cut off value (day–5: 61.9% vs. 33.3%(p = 0.0023), day–6: 47.4% vs. 7.6%(p = 0.0003)). Multivariate analysis adjusted for patient age showed that KIDSCORE™D5 correlated with clinical pregnancy rates of days 5 and 6 of blastocyst transfer with AUCs of 0.66 and 0.73, respectively.

Limitations, reasons for caution

This study had a small sample size, and it was a retrospective single-center study. In addition, the relationship between KIDSCORE™D5 and clinical pregnancy rates may vary among facilities. Therefore, a prospective multicenter validation is necessary.

Wider implications of the findings: Our study results indicated that KIDSCORE™D5 predicted clinical pregnancy and that morphokinetic parameters related to clinical pregnancy were similar between day–5 and day–6 blastocysts. Hence, morphokinetic evaluation can serve as a criterion for deciding which of multiple day–6 blastocysts can be transferred.

Trial registration number

Not applicable

rs1944919 on chromosome 11q23.1 and its effector genes COLCA1/COLCA2 confer susceptibility to primary biliary cholangitis

Primary biliary cholangitis (PBC) is a chronic, progressive cholestatic liver disease in which intrahepatic bile ducts are destroyed by an autoimmune reaction. Our previous genome-wide association study (GWAS) identified chromosome 11q23.1 as a susceptibility gene locus for PBC in the Japanese population. Here, high-density association mapping based on single nucleotide polymorphism (SNP) imputation and in silico/in vitro functional analyses identified rs1944919 as the primary functional variant. Expression-quantitative trait loci analyses showed that the PBC susceptibility allele of rs1944919 was significantly associated with increased COLCA1/COLCA2 expression levels. Additionally, the effects of rs1944919 on COLCA1/COLCA2 expression levels were confirmed using genotype knock-in versions of cell lines constructed using the CRISPR/Cas9 system and differed between rs1944919-G/G clones and -T/T clones. To our knowledge, this is the first study to demonstrate the contribution of COLCA1/COLCA2 to PBC susceptibility.

A pro-diabetogenic mtDNA polymorphism in the mitochondrial-derived peptide, MOTS-c

Type 2 Diabetes (T2D) is an emerging public health problem in Asia. Although ethnic specific mtDNA polymorphisms have been shown to contribute to T2D risk, the functional effects of the mtDNA polymorphisms and the therapeutic potential of mitochondrial-derived peptides at the mtDNA polymorphisms are underexplored. Here, we showed an Asian-specific mitochondrial DNA variation m.1382A>C (rs111033358) leads to a K14Q amino acid replacement in MOTS-c, an insulin sensitizing mitochondrial-derived peptide. Meta-analysis of three cohorts (n = 27,527, J-MICC, MEC, and TMM) show that males but not females with the C-allele exhibit a higher prevalence of T2D. In J-MICC, only males with the C-allele in the lowest tertile of physical activity increased their prevalence of T2D, demonstrating a kinesio-genomic interaction. High-fat fed, male mice injected with MOTS-c showed reduced weight and improved glucose tolerance, but not K14Q-MOTS-c treated mice. Like the human data, female mice were unaffected. Mechanistically, K14Q-MOTS-c leads to diminished insulin-sensitization in vitro. Thus, the m.1382A>C polymorphism is associated with susceptibility to T2D in men, possibly interacting with exercise, and contributing to the risk of T2D in sedentary males by reducing the activity of MOTS-c.

Facial UV photo imaging for skin pigmentation assessment using conditional generative adversarial networks

Skin pigmentation is associated with skin damages and skin cancers, and ultraviolet (UV) photography is used as a minimally invasive mean for the assessment of pigmentation. Since UV photography equipment is not usually available in general practice, technologies emphasizing pigmentation in color photo images are desired for daily care. We propose a new method using conditional generative adversarial networks, named UV-photo Net, to generate synthetic UV images from color photo images. Evaluations using color and UV photo image pairs taken by a UV photography system demonstrated that pigment spots were well reproduced in synthetic UV images by UV-photo Net, and some of the reproduced pigment spots were difficult to be recognized in color photo images. In the pigment spot detection analysis, the rate of pigment spot areas in cheek regions for synthetic UV images was highly correlated with the rate for UV photo images (Pearson's correlation coefficient 0.92). We also demonstrated that UV-photo Net was effective for floating up pigment spots for photo images taken by a smartphone camera. UV-photo Net enables an easy assessment of pigmentation from color photo images and will promote self-care of skin damages and early signs of skin cancers for preventive medicine.

A genotype imputation method for de-identified haplotype reference information by using recurrent neural network

Genotype imputation estimates the genotypes of unobserved variants using the genotype data of other observed variants based on a collection of haplotypes for thousands of individuals, which is known as a haplotype reference panel. In general, more accurate imputation results were obtained using a larger size of haplotype reference panel. Most of the existing genotype imputation methods explicitly require the haplotype reference panel in precise form, but the accessibility of haplotype data is often limited, due to the requirement of agreements from the donors. Since de-identified information such as summary statistics or model parameters can be used publicly, imputation methods using de-identified haplotype reference information might be useful to enhance the quality of imputation results under the condition where the access of the haplotype data is limited. In this study, we proposed a novel imputation method that handles the reference panel as its model parameters by using bidirectional recurrent neural network (RNN). The model parameters are presented in the form of de-identified information from which the restoration of the genotype data at the individual-level is almost impossible. We demonstrated that the proposed method provides comparable imputation accuracy when compared with the existing imputation methods using haplotype datasets from the 1000 Genomes Project (1KGP) and the Haplotype Reference Consortium. We also considered a scenario where a subset of haplotypes is made available only in de-identified form for the haplotype reference panel. In the evaluation using the 1KGP dataset under the scenario, the imputation accuracy of the proposed method is much higher than that of the existing imputation methods. We therefore conclude that our RNN-based method is quite promising to further promote the data-sharing of sensitive genome data under the recent movement for the protection of individuals’ privacy.

Genotype imputation estimates the genotypes of unobserved variants using the genotype data of other observed variants based on a collection of genome data of a large number of individuals called a reference panel. In general, more accurate imputation results are obtained using a larger size of the reference panel. Although most of the existing imputation methods use the reference panel in an explicit form, the accessibility of genome data is often limited due to the requirement of agreements from the donors. We thus proposed a new imputation method that handles the reference panel as its model parameters by using bidirectional recurrent neural network. Since it is almost impossible to restore genome data at the individual-level from the model parameters, they can be shared publicly as the de-identified information even when the accessibility of the original reference panel is limited. We demonstrate that the proposed method provides comparable imputation accuracy with the existing methods. We also considered a scenario where a part of the genome data is made available only in de-identified form for the reference panel and have shown that the imputation accuracy of the proposed method is much higher than that of the existing methods under the scenario.

Integrated GWAS and mRNA Microarray Analysis Identified IFNG and CD40L as the Central Upstream Regulators in Primary Biliary Cholangitis

Genome‐wide association studies (GWASs) in European and East Asian populations have identified more than 40 disease‐susceptibility genes in primary biliary cholangitis (PBC). The aim of this study is to computationally identify disease pathways, upstream regulators, and therapeutic targets in PBC through integrated GWAS and messenger RNA (mRNA) microarray analysis. Disease pathways and upstream regulators were analyzed with ingenuity pathway analysis in data set 1 for GWASs (1,920 patients with PBC and 1,770 controls), which included 261 annotated genes derived from 6,760 single‐nucleotide polymorphisms (P < 0.00001), and data set 2 for mRNA microarray analysis of liver biopsy specimens (36 patients with PBC and 5 normal controls), which included 1,574 genes with fold change >2 versus controls (P < 0.05). Hierarchical cluster analysis and categorization of cell type–specific genes were performed for data set 2. There were 27 genes, 10 pathways, and 149 upstream regulators that overlapped between data sets 1 and 2. All 10 pathways were immune‐related. The most significant common upstream regulators associated with PBC disease susceptibility identified were interferon‐gamma (IFNG) and CD40 ligand (CD40L). Hierarchical cluster analysis of data set 2 revealed two distinct groups of patients with PBC by disease activity. The most significant upstream regulators associated with disease activity were IFNG and CD40L. Several molecules expressed in B cells, T cells, Kupffer cells, and natural killer–like cells were identified as potential therapeutic targets in PBC with reference to a recently reported list of cell type–specific gene expression in the liver. Conclusion: Our integrated analysis using GWAS and mRNA microarray data sets predicted that IFNG and CD40L are the central upstream regulators in both disease susceptibility and activity of PBC and identified potential downstream therapeutic targets.

Interlaboratory Evaluation of Two Enzyme-Linked Immunosorbent Assay Kits for the Determination of Crustacean Protein in Processed Foods

The labeling of foods containing material derived from crustaceans such as shrimp and crab is to become mandatory in Japan because of increases in the number of allergy patients. To ensure proper labeling, 2 novel sandwich enzyme-linked immunosorbent assay (ELISA) kits for the determination of crustacean protein in processed foods, the N kit (Nissui Pharmaceutical Co., Ltd, Ibaraki, Japan) and the M kit (Maruha Nichiro Holdings, Inc., Ibaraki, Japan), have been developed. Five types of model processed foods containing 10 and/or 11.9 g/g crustacean soluble protein were prepared for interlaboratory evaluation of the performance of these kits. The N kit displayed a relatively high level of reproducibility relative standard deviation (interlaboratory precision; 4.08.4 RSDR) and sufficient recovery (6586) for all the model processed foods. The M kit displayed sufficient reproducibility (17.620.5 RSDR) and a reasonably high level of recovery (82103). The repeatability relative standard deviation (RSDr) values regarding the detection of crustacean proteins in the 5 model foods were mostly &lt;5.1 RSDr for the N kit and 9.9 RSDr for the M kit. In conclusion, the results of this interlaboratory evaluation suggest that both these ELISA kits would be very useful for detecting crustacean protein in processed foods.

Determination of Walnut Protein in Processed Foods by Enzyme-Linked Immunosorbent Assay: Interlaboratory Study

Because food allergens from tree nuts, including walnuts, are a frequent cause of adverse food reactions for allergic patients, the labeling of foods containing ingredients derived from tree nuts is required in numerous countries. According to Japanese regulations, the labeling of food products containing walnuts is recommended. To ensure proper labeling, a novel sandwich ELISA kit for the determination of walnut protein in processed foods (Walnut Protein [2S-Albumin] Kit; Morinaga Institute of Biological Science, Inc.; walnut kit) has been developed. We prepared seven types of incurred samples (model processed foods: biscuits, bread, sponge cake, orange juice, jelly, chicken meatballs, and rice gruel) containing 10 g walnut soluble protein/g of food for use in interlaboratory evaluations of the walnut kit. The walnut kit displayed sufficient reproducibility relative standard deviations (interlaboratory precision: 5.89.9 RSDR) and a high level of recovery (81119) for all the incurred samples. All the repeatability relative standard deviation (RSDr) values for the incurred samples that were examined were less than 6.0. The results of this interlaboratory evaluation suggested that the walnut kit could be used as a precise and reliable tool for determination of walnut protein in processed foods.

Enzyme-Linked Immunosorbent Assay Kit for the Determination of Soybean Protein in Processed Foods: Interlaboratory Evaluation

The labeling of foods containing ingredients derived from soybean is recommended in Japan because of an increasing number of patients who are allergic to soybeans. To ensure proper labeling, a novel sandwich ELISA kit for the determination of soybean protein in processed foods (FASTKIT Ver. II, Soybean, Nippon Meat Packers, Inc.; soy kit) has been developed. Five types of incurred samples (model processed foods: rice gruel, sausage, sweet adzuki bean soup, sweet potato cake, and tomato sauce) containing 10 g soybean soluble protein/g food were prepared for use in interlaboratory evaluations of the soy kit. The soy kit displayed a sufficient RSDR value (interlaboratory precision: 9.313.4 RSDR) and a high level of recovery (97114) for all the incurred samples. The RSDr value for the incurred samples was mostly &lt;4.8. The results of this interlaboratory evaluation suggest that the soy kit can be used as a precise and reliable tool for the determination of soybean proteins in processed foods.

P1828Aortic calcification detected by computed tomography and aortic vulnerable plaques: aortic angioscopy study

Background

Aortic calcification is associated with atherosclerotic risk factors and an increased risk of death and cardiovascular disease. However, the relationships aortic calcification and aortic plaque instability are not yet elucidated. Recently, some reports showed non-obstructive aortic angioscopy seemed to visualize atherosclerotic changes of aortic wall more clearly compared with computed tomography (CT). The purpose of this study was to evaluate whether aortic calcification is associated with aortic vulnerable plaques in patients with cardiovascular disease.

Methods

We investigated 60 consecutive patients with confirmed or suspected coronary artery disease who underwent both aortic angioscopy and CT. The AC volume (ACV) was measured using the volume-rendering method by extracting the area >130 HU within the whole aorta. ACV index (ACVI) was defined as ACV divided by the body surface area. We evaluated the number of ruptured plaque (RP), ulceration and fissure by aortic angioscopy in the whole aorta. We excluded 4 hemodialysis patients. All patients were divided into the median value of ACVI.

Results

The mean age of patients was 68±10. The median of ACVI was 10.7 ml/m2 [3.9–22.7]. High ACVI patients had significantly greater number of RP, ulceration and atheromatous plaques detected by aortic angioscopy compared with those of low ACVI (2.2±2.7 vs 0.8±1.1, p=0.033, 1.6±1.2 vs 0.9±1.0, p=0.041, 4.0±3.1 vs 1.9±1.8, p=0.009, respectively). Furthermore, the patients without aortic calcification did not have RP at all. In a multivariate model, the number of the atheromatous plaques was independently associated with high ACVI (odds ratio 1.57, 95% confidence interval 1.07–2.69, p=0.018)

Conclusions

Aortic calcification detected by CT was related to aortic vulnerable plaques in patients with cardiovascular disease.

6125High wall shear stress predicts plaque rupture of the aortic arch: computational fluid dynamics model and non-obstructive general angioscopy study

Background

Wall shear stress (WSS) has been considered as a major determinant of aortic atherosclerosis. Recently, non-obstructive general angioscopy (NOGA) was developed to be able to visualize a variety of its atherosclerotic pathology, including in vivo ruptured plaque (RP) in the aorta. We, therefore, investigated the relationship between NOGA derived RP in the aortic arch and the stereographic distribution of WSS by using computational fluid dynamics modeling (CFD) on three-dimensional CT angiography (3D-CT).

Methods

We investigated 30 consecutive patients who underwent 3D-CT before and NOGA during coronary angiography. WSS in the aortic arch was measured with an application of CFD based on finite element method by using uniform inlet and outlet flow conditions. Aortic RP was detected by NOGA.

Results

The maximum and mean values of WSS were 67.2±29.2 Pa and 2.4±0.6 Pa. A total of 18 RPs was detected by NOGA. The patients with a distinct RP showed a significantly higher maximum WSS in the whole aortic arch, and the greater and lesser curvature of the aortic arch than those without it (73.3±29.0 Pa vs 50.4±15.2 Pa, p=0.035, 95.0±27.5 Pa vs 42.8±25.2 Pa, p=0.003, 70.8±29.3 Pa vs 46.1±11.9 Pa, p=0.013, respectively), whereas there was no significant difference in the mean WSS between those with and without it. In a multivariate analysis, the maximum value of WSS was an independent predictor of RP in the aortic arch (odds ratio 1.05, 95% confidence interval 1.01–1.13, p=0.019).

Representative picture of WSS and NOGA

Conclusions

Aortic RP detected by NOGA was strongly associated with the higher maximum WSS in the aortic arch derived by CFD using 3D-CT. Maximum WSS may explain the underlying mechanism of not only aortic atherosclerosis, but also aortic RP.

P1831Prediction of cardiovascular events by atheromatous plaques detected by non-obstructive general angioscopy: two-year results of EAST-NOGA Registry

Background

Non-obstructive general angioscopy (NOGA) has revealed the intimal damages or atheromatous plaques as well as its spontaneous rupture of the aorta. Recent study revealed that plaque debris or different size of cholesterol crystals were detected in the blood above the spontaneous ruptured aortic plaque observed by NOGA and these plaque materials might cause the peripheral organ damages as the embolic source. These various morphological changes may cause the acute aortic events or atheroembolic events on the peripheral organs, such as brain, kidney, peripheral artery and so on.

Purpose

EAST-NOGA (Evaluation of AtheroSclerotic and rupture events by Non-Obstructive General Angioscopy) is a multi-center prospective observational study to assess the relationship between the findings of NOGA and future cardiovascular events.

Methods

Five hundred and seventy-seven patients with atherosclerotic cardiovascular disease who underwent NOGA study. The major cerebrocardiovascular events including cardiovascular death, non-fatal myocardial infarction, non-fatal cerebral infarction, and acute aortic syndrome were accumulated during the 2-year follow-up after NOGA study.

Results

The median number of aortic atheromatous plaques was 6 [IQR: 3–12]. A total of 514 patients were followed up (89.1%). The mean follow-up duration was 757±120 days. Major adverse cardiovascular events developed in 23 (4.5%) during 2 years follow-up. Patients with MACE and cerebral infarction, had significantly greater number of aortic atheromatous plaques (11 [5–19] vs. 6 [3–11], p<0.001, 12 [4–20] vs. 6 [3–12], p=0.014, respectively). In a univariate analysis, the number of aortic atheromatous plaques and ruptured plaque were significant predictors of MACE (HR: 1.09 95% confidence interval 1.05–1.14, p<0.001) and (HR: 1.12, 95% confidence interval 1.02–1.23, p=0.02). In a multivariate logistic analysis, the number of aortic atheromatous plaques is one of the independent predictors of MACE (HR 1.05, 95% confidence interval 1.00–1.10, p=0.032).

Conclusion

The number of atheromatous plaques identified by NOGA has a significant relation to the onset of cerebral infarction, which suggest the atheromatous plaque were vulnerable and ruptured spontaneously, then cause the aortogenic cerebral infarction. The NOGA study would be useful for predicting the futured atheroembolic events.

Acknowledgement/Funding

None

A Polymorphism rs6726395 in Nrf2 Contributes to the Development of Emphysema-Associated Age in Smokers Without COPD

INTRODUCTION

Several studies have reported that single nucleotide polymorphisms (SNPs) in the gene encoding NF-E2-related factor 2 (Nrf2) contribute to airflow limitations in smokers without COPD. Although small airway lesions and emphysema contribute cooperatively to airflow limitation, the relationship between Nrf2 SNPs and the development of emphysema in smokers without COPD is not well understood.

METHODS

Healthy subjects who underwent an annual health checkup with computed tomography (CT) of the chest at Osaka City University Hospital were prospectively recruited. The percentage of low-attenuation area (%LAA) on chest CT was quantified, and correlations between %LAA, Nrf2 SNP [rs6726395 (G/A)] genotypes, and clinical characteristics were examined.

RESULTS

A total of 245 subjects without COPD [non-/light-smoker: 153 (62.4%) and smoker: 92 (37.6%)] were enrolled. The %LAA in the upper lung field was higher than that in the lower lung field (p < 0.001). The %LAA in smokers was significantly higher than that in non-/light-smokers (p = 0.021). The %LAA showed significant but weak correlation with age in all subjects (r = 0.141, p = 0.028). Divided by genotype, the %LAA of the upper lung field was significantly correlated with age in smokers with genotype GG (wild type) (r = 0.333, p = 0.022), but was not significantly correlated with age in smokers with genotype AG/AA. These correlations were not observed in non-/light smokers.

CONCLUSION

A polymorphism rs6726395 in Nrf2 can contribute to the development of emphysema-associated aging in smokers. The Nrf2 SNP may be a predictive factor for smoking-induced emphysema, and genotyping of Nrf2 SNP may serve as biomarker for emphysema prevention.

Construction of JRG (Japanese reference genome) with single-molecule real-time sequencing

In recent genome analyses, population-specific reference panels have indicated important. However, reference panels based on short-read sequencing data do not sufficiently cover long insertions. Therefore, the nature of long insertions has not been well documented. Here, we assembled a Japanese genome using single-molecule real-time sequencing data and characterized insertions found in the assembled genome. We identified 3691 insertions ranging from 100 bps to ~10,000 bps in the assembled genome relative to the international reference sequence (GRCh38). To validate and characterize these insertions, we mapped short-reads from 1070 Japanese individuals and 728 individuals from eight other populations to insertions integrated into GRCh38. With this result, we constructed JRGv1 (Japanese Reference Genome version 1) by integrating the 903 verified insertions, totaling 1,086,173 bases, shared by at least two Japanese individuals into GRCh38. We also constructed decoyJRGv1 by concatenating 3559 verified insertions, totaling 2,536,870 bases, shared by at least two Japanese individuals or by six other assemblies. This assembly improved the alignment ratio by 0.4% on average. These results demonstrate the importance of refining the reference assembly and creating a population-specific reference genome. JRGv1 and decoyJRGv1 are available at the JRG website.

Researchers in Japan have assembled a Japanese reference genome, which includes sequences missing from the international reference genome, as well as others specific to East Asian populations. A team led by Masao Nagasaki and Masayuki Yamamoto sequenced a Japanese individual using a method, which produces longer sequences than previous technologies. Using this approach, they identified thousands of sequences spanning 2.5 million bases, which were absent in the international reference genome. Many of these were sequences able to move within the genome. They showed that the majority of these sequences are also present in early humans and chimpanzees, demonstrating that their absence from the current reference is due to deletions or limitations of earlier sequencing methodologies. In addition to providing a population-specific reference, these findings demonstrate the importance of continually improving the international reference genome.

Genome analyses for the Tohoku Medical Megabank Project towards establishment of personalized healthcare

Personalized healthcare (PHC) based on an individual's genetic make-up is one of the most advanced, yet feasible, forms of medical care. The Tohoku Medical Megabank (TMM) Project aims to combine population genomics, medical genetics and prospective cohort studies to develop a critical infrastructure for the establishment of PHC. To date, a TMM CommCohort (adult general population) and a TMM BirThree Cohort (birth+three-generation families) have conducted recruitments and baseline surveys. Genome analyses as part of the TMM Project will aid in the development of a high-fidelity whole-genome Japanese reference panel, in designing custom single-nucleotide polymorphism (SNP) arrays specific to Japanese, and in estimation of the biological significance of genetic variations through linked investigations of the cohorts. Whole-genome sequencing from >3,500 unrelated Japanese and establishment of a Japanese reference genome sequence from long-read data have been done. We next aim to obtain genotype data for all TMM cohort participants (>150,000) using our custom SNP arrays. These data will help identify disease-associated genomic signatures in the Japanese population, while genomic data from TMM BirThree Cohort participants will be used to improve the reference genome panel. Follow-up of the cohort participants will allow us to test the genetic markers and, consequently, contribute to the realization of PHC.

POGLUT1, the putative effector gene driven by rs2293370 in primary biliary cholangitis susceptibility locus chromosome 3q13.33

Primary biliary cholangitis (PBC) is a chronic and cholestatic autoimmune liver disease caused by the destruction of intrahepatic small bile ducts. Our previous genome-wide association study (GWAS) identified six susceptibility loci for PBC. Here, in order to further elucidate the genetic architecture of PBC, a GWAS was performed on an additional independent sample set, then a genome-wide meta-analysis with our previous GWAS was performed based on a whole-genome single nucleotide polymorphism (SNP) imputation analysis of a total of 4,045 Japanese individuals (2,060 cases and 1,985 healthy controls). A susceptibility locus on chromosome 3q13.33 (including ARHGAP31, TMEM39A, POGLUT1, TIMMDC1, and CD80) was previously identified both in the European and Chinese populations and was replicated in the Japanese population (OR = 0.7241, P = 3.5 × 10^-9). Subsequent in silico and in vitro functional analyses identified rs2293370, previously reported as the top-hit SNP in this locus in the European population, as the primary functional SNP. Moreover, e-QTL analysis indicated that the effector gene of rs2293370 was Protein O-Glucosyltransferase 1 (POGLUT1) (P = 3.4 × 10^-8). This is the first study to demonstrate that POGLUT1 and not CD80 is the effector gene regulated by the primary functional SNP rs2293370, and that increased expression of POGLUT1 might be involved in the pathogenesis of PBC.

NFKB1 and MANBA Confer Disease Susceptibility to Primary Biliary Cholangitis via Independent Putative Primary Functional Variants

BACKGROUND & AIMS

Primary biliary cholangitis (PBC) is a chronic and cholestatic liver disease that eventually leads to cirrhosis and hepatic failure. We recently identified several susceptibility genes included NFKB1 and MANBA for PBC in the Japanese population by genome-wide association study. However, the primary functional variants in the NFKB1/MANBA region and the molecular mechanism for conferring disease susceptibility to PBC have not yet been clarified.

METHODS

We performed high-density association mapping based on a single-nucleotide polymorphism (SNP) imputation analysis, using data from a whole-genome sequence reference panel of 1070 Japanese individuals and the previous genome-wide association study (1389 PBC patients, 1508 healthy controls). Among SNPs (P < 5.0 × 10^-7) in the NFKB1/MANBA region, putative primary functional variants and the molecular mechanism for conferring disease susceptibility to PBC were identified by in silico/in vitro functional analysis.

RESULTS

Among the SNPs in the NFKB1/MANBA region, rs17032850 and rs227361, which changed the binding of transcription factors lymphoid enhancer-binding factor 1 (LEF-1) and retinoid X receptor α (RXRα), respectively, were identified as putative primary functional variants that regulate gene expression. In addition, expression-quantitative trait locus data and gene editing using a clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9 system supported the potential role of rs17032850 and rs227361 in regulating NFKB1 and MANBA expression, respectively.

CONCLUSIONS

We identified independent putative primary functional variants in NFKB1/MANBA and showed the distinct molecular mechanism by which each putative primary functional variant conferred susceptibility to PBC. Our approach was useful to dissect the pathogenesis not only of PBC, but also other digestive diseases in which NFKB1/MANBA has been reported as a susceptibility locus.

Regional genetic differences among Japanese populations and performance of genotype imputation using whole-genome reference panel of the Tohoku Medical Megabank Project

Background

Genotype imputation from single-nucleotide polymorphism (SNP) genotype data using a haplotype reference panel consisting of thousands of unrelated individuals from populations of interest can help to identify strongly associated variants in genome-wide association studies. The Tohoku Medical Megabank (TMM) project was established to support the development of precision medicine, together with the whole-genome sequencing of 1070 human genomes from individuals in the Miyagi region (Northeast Japan) and the construction of the 1070 Japanese genome reference panel (1KJPN). Here, we investigated the performance of 1KJPN for genotype imputation of Japanese samples not included in the TMM project and compared it with other population reference panels.

Results

We found that the 1KJPN population was more similar to other Japanese populations, Nagahama (south-central Japan) and Aki (Shikoku Island), than to East Asian populations in the 1000 Genomes Project other than JPT, suggesting that the large-scale collection (more than 1000) of Japanese genomes from the Miyagi region covered many of the genetic variations of Japanese in mainland Japan. Moreover, 1KJPN outperformed the phase 3 reference panel of the 1000 Genomes Project (1KGPp3) for Japanese samples, and IKJPN showed similar imputation rates for the TMM and other Japanese samples for SNPs with minor allele frequencies (MAFs) higher than 1%.

Conclusions

1KJPN covered most of the variants found in the samples from areas of the Japanese mainland outside the Miyagi region, implying 1KJPN is representative of the Japanese population’s genomes. 1KJPN and successive reference panels are useful genome reference panels for the mainland Japanese population. Importantly, the addition of whole genome sequences not included in the 1KJPN panel improved imputation efficiencies for SNPs with MAFs under 1% for samples from most regions of the Japanese archipelago.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4942-0) contains supplementary material, which is available to authorized users.

Strong Association of the HLA-DR/DQ Locus with Childhood Steroid-Sensitive Nephrotic Syndrome in the Japanese Population

Background Nephrotic syndrome is the most common cause of chronic glomerular disease in children. Most of these patients develop steroid-sensitive nephrotic syndrome (SSNS), but the loci conferring susceptibility to childhood SSNS are mainly unknown.Methods We conducted a genome-wide association study (GWAS) in the Japanese population; 224 patients with childhood SSNS and 419 adult healthy controls were genotyped using the Affymetrix Japonica Array in the discovery stage. Imputation for six HLA genes (HLA-A, -C, -B, -DRB1, -DQB1, and -DPB1) was conducted on the basis of Japanese-specific references. We performed genotyping for HLA-DRB1/-DQB1 using a sequence-specific oligonucleotide-probing method on a Luminex platform. Whole-genome imputation was conducted using a phased reference panel of 2049 healthy Japanese individuals. Replication was performed in an independent Japanese sample set including 216 patients and 719 healthy controls. We genotyped candidate single-nucleotide polymorphisms using the DigiTag2 assay.Results The most significant association was detected in the HLA-DR/DQ region and replicated (rs4642516 [minor allele G], combined P_allelic=7.84×10^-23; odds ratio [OR], 0.33; 95% confidence interval [95% CI], 0.26 to 0.41; rs3134996 [minor allele A], combined P_allelic=1.72×10^-25; OR, 0.29; 95% CI, 0.23 to 0.37). HLA-DRB1*08:02 (Pc=1.82×10^-9; OR, 2.62; 95% CI, 1.94 to 3.54) and HLA-DQB1*06:04 (Pc=2.09×10^-12; OR, 0.10; 95% CI, 0.05 to 0.21) were considered primary HLA alleles associated with childhood SSNS. HLA-DRB1*08:02-DQB1*03:02 (Pc=7.01×10^-11; OR, 3.60; 95% CI, 2.46 to 5.29) was identified as the most significant genetic susceptibility factor.Conclusions The most significant association with childhood SSNS was detected in the HLA-DR/DQ region. Further HLA allele/haplotype analyses should enhance our understanding of molecular mechanisms underlying SSNS.