Search

Abstract Cell states are regulated by the response of signaling pathways to receptor ligand-binding and intercellular interactions. High-resolution imaging has been attempted to explore the dynamics of these processes and, recently, multiplexed imaging has profiled cell states by achieving a comprehensive acquisition of spatial protein information from cells. However, the specificity of antibodies is still compromised when visualizing activated signals. Here, we develop Precise Emission Canceling Antibodies (PECAbs) that have cleavable fluorescent labeling. PECAbs enable high-specificity sequential imaging using hundreds of antibodies, allowing for reconstruction of the spatiotemporal dynamics of signaling pathways. Additionally, combining this approach with seq-smFISH can effectively classify cells and identify their signal activation states in human tissue. Overall, the PECAb system can serve as a comprehensive platform for analyzing complex cell processes.

Leukocyte immunoglobulin (Ig)-like receptors (LILRs) on human chromosome 19q13.4 encode 11 immunoglobulin superfamily receptors, exhibiting genetic diversity within and between human populations. Among the LILR genes, the genomic region surrounding LILRB3 and LILRA6 has yet to be fully characterized due to their significant sequence homology, which makes it difficult to differentiate between them. To examine the LILRB3 and LILRA6 genomic region, a tool named JoGo-LILR CN Caller, which can call copy number from short-read whole genome sequencing (srWGS) data, was applied to an extensive international srWGS dataset comprising 2,504 samples. During this process, a previously unreported loss of both LILRB3 and LILRA6 was detected in three samples. Using long-read sequencing of these samples, we have discovered a novel large deletion (33,692 bp) in the LILRB3 and LILRA6 genomic regions in the Japanese population. This deletion spanned three genes, LILRB3, LILRA6, and LILRB5, resulting in LILRB3 exons 12-13 being located immediately downstream of LILRB5 exons 1-12 with the loss of LILRA6, suggesting the potential expression of a hybrid gene between LILRB5 and LILRB3 (LILRB5-3). Transcription and subsequent translation of the LILRB5-3 hybrid gene were also verified. The hybrid junction was located within the intracellular domain, resulting in an LILRB5 extracellular domain fused to a partial LILRB3 intracellular domain with three immunoreceptor tyrosine-based inhibitory motifs (ITIMs), suggesting that LILRB5-3 acquired a novel signaling function. Further application of the JoGo-LILR tool to srWGS samples suggested the presence of the LILRB5-3 hybrid gene in the CEU population. Our findings provide insight into the genetic and functional diversity of the LILR family.

Abstract In the field of genomic medical research, the amount of large-scale information continues to increase due to advances in measurement technologies, such as high-performance sequencing and spatial omics, as well as the progress made in genomic cohort studies involving more than one million individuals. Therefore, researchers require more computational resources to analyze this information. Here, we introduce a hybrid cloud system consisting of an on-premise supercomputer, science cloud, and public cloud at the Kyoto University Center for Genomic Medicine in Japan as a solution. This system can flexibly handle various heterogeneous computational resource-demanding bioinformatics tools while scaling the computational capacity. In the hybrid cloud system, we demonstrate the way to properly perform joint genotyping of whole-genome sequencing data for a large population of 11,238, which can be a bottleneck in sequencing data analysis. This system can be one of the reference implementations when dealing with large amounts of genomic medical data in research centers and organizations.

Abstract Secure storage and secondary use of individual human genome data is increasingly important for genome research and personalized medicine. Currently, it is necessary to store the whole genome sequencing information (FASTQ data), which enables detections of de novo mutations and structural variations in the analysis of hereditary diseases and cancer. Furthermore, bioinformatics tools to analyze FASTQ data are frequently updated to improve the precision and recall of detected variants. However, existing secure secondary use of data, such as multi-party computation or homomorphic encryption, can handle only a limited algorithms and usually requires huge computational resources. Here, we developed a high-performance one-stop system for large-scale genome data analysis with secure secondary use of the data by the data owner and multiple users with different levels of data access control. Our quantum secure cloud system is a distributed secure genomic data analysis system (DSGD) with a “trusted server” built on a quantum secure cloud, the information-theoretically secure Tokyo QKD Network. The trusted server will be capable of deploying and running a variety of sequencing analysis hardware, such as GPUs and FPGAs, as well as CPU-based software. We demonstrated that DSGD achieved comparable throughput with and without encryption on the trusted server Therefore, our system is ready to be installed at research institutes and hospitals that make diagnoses based on whole genome sequencing on a daily basis.

Abstract Background Ligation of CD28 with ligands such as CD80 or CD86 provides a critical second signal alongside antigen presentation by class II major histocompatibility complex expressed on antigen-presenting cells through the T cell antigen receptor for naïve T cell activation. A number of studies suggested that CD28 plays an important role in the pathogenesis of various human diseases. Recent genome-wide association studies (GWASs) identified CD28 as a susceptibility locus for lymphocyte and eosinophil counts, multiple sclerosis, ulcerative colitis, celiac disease, rheumatoid arthritis, asthma, and primary biliary cholangitis. However, the primary functional variant and molecular mechanisms of disease susceptibility in this locus remain to be elucidated. This study aimed to identify the primary functional variant from thousands of genetic variants in the CD28 locus and elucidate its functional effect on the CD28 molecule. Results Among the genetic variants exhibiting stronger linkage disequilibrium (LD) with all GWAS-lead variants in the CD28 locus, rs2013278, located in the Rbfox binding motif related to splicing regulation, was identified as a primary functional variant related to multiple immunological traits. Relative endogenous expression levels of CD28 splicing isoforms (CD28i and CD28Δex2) compared with full-length CD28 in allele knock-in cell lines generated using CRISPR/Cas9 were directly regulated by rs2013278 (P

Primary biliary cholangitis (PBC) is a rare autoimmune disease with a clear predisposition for human leukocyte antigen (HLA)-DR/DQ-associated loss of immune tolerance for the E2 component of the pyruvate dehydrogenase complex. Three-field-resolution HLA imputation of 1,670 Japanese PBC patients and 2,328 healthy controls was conducted using Japanese population-specific HLA reference panels. Eighteen previously reported Japanese PBC-associated HLA alleles were confirmed and extended to 3-field-resolution, including HLA-DRB1*08:03 to HLA-DRB1*08:03:02, HLA-DQB1*03:01 to HLA-DQB1*03:01:01, HLA-DQB1*04:01 to HLA-DQB1*04:01:01 and HLA-DQB1*06:04 to HLA-DQB1*06:04:01. In addition, additional significant novel HLA alleles were identified, including 3 novel susceptible HLA-DQA1 alleles: HLA-DQA1*03:03:01, HLA-DQA1*04:01:01, HLA-DQA1*01:04:01 and 1 novel protective HLA-DQA1 allele, HLA-DQA1*05:05:01. In addition, PBC patients carrying HLA-DRB1*15:01:01 and HLA-DQA1*03:03:01 would have a higher predisposition toward developing concomitant autoimmune hepatitis (AIH). Further, late-stage and symptomatic PBC shared the same susceptible HLA alleles of HLA-A*26:01:01, HLA-DRB1*09:01:02 and HLA-DQB1*03:03:02. Lastly, HLA-DPB1*05:01:01 was identified as a potential risk HLA allele for development of hepatocellular carcinoma (HCC) in PBC patients. In conclusion, we have extended the current knowledge of HLA allele associations to 3-field resolution and identified novel HLA allele associations with predisposition risk, staging, symptomatic state, and AIH and HCC events for Japanese PBC patients.

Abstract The elucidation of dynamic metabolomic changes during gestation is particularly important for the development of methods to evaluate pregnancy status or achieve earlier detection of pregnancy-related complications. Some studies have constructed models to evaluate pregnancy status and predict gestational age using omics data from blood biospecimens; however, less invasive methods are desired. Here we propose a model to predict gestational age, using urinary metabolite information. In our prospective cohort study, we collected 2741 urine samples from 187 healthy pregnant women, 23 patients with hypertensive disorders of pregnancy, and 14 patients with spontaneous preterm birth. Using gas chromatography-tandem mass spectrometry, we identified 184 urinary metabolites that showed dynamic systematic changes in healthy pregnant women according to gestational age. A model to predict gestational age during normal pregnancy progression was constructed; the correlation coefficient between actual and predicted weeks of gestation was 0.86. The predicted gestational ages of cases with hypertensive disorders of pregnancy exhibited significant progression, compared with actual gestational ages. This is the first study to predict gestational age in normal and complicated pregnancies by using urinary metabolite information. Minimally invasive urinary metabolomics might facilitate changes in the prediction of gestational age in various clinical settings.

Abstract Stevens–Johnson syndrome (SJS) and its severe condition with extensive skin detachment and a poor prognosis, toxic epidermal necrolysis (TEN), are immunologically mediated severe cutaneous reactions of the skin and mucous membranes such as the ocular surface. Genetic variations on the HLA-A and other autosomal genes have been identified as risk factors for cold medicine-related SJS/TEN with severe ocular complications (CM-SJS/TEN with SOC). Using a whole-genome sequencing (WGS) approach, we explored other susceptible variants of CM-SJS/TEN with SOC, especially among rare variants and structural variants (SVs). WGS was performed on samples from 133 patients with CM-SJS/TEN with SOC and 418 healthy controls to obtain single nucleotide polymorphisms (SNPs) and SVs. Genome-wide association tests were performed with these variants. Our genome-wide association test reproduced the associations of the common variants of HLA-A and loci on chromosome 16q12.1. We also identified novel associations of SVs on these loci and an aggregation of rare coding variants on the TPRM8 gene. In silico gene expression analysis on the HLA-A locus revealed that the SNP (rs12202296), which was significantly associated with susceptibility to CM-SJS/TEN with SOC, was correlated to an increase in HLA-A expression levels in the whole blood (P = 2.9 × 10−17), from the GTEx database. The majority of variants that were significantly associated with CM-SJS/TEN with SOC were found in non-coding regions, indicating the regulatory role of genetic variations in the pathogenesis of CM-SJS/TEN with SOC.

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

The development of liver cancer in patients with hepatitis B is a major problem, and several models have been reported to predict the development of liver cancer. However, no predictive model involving human genetic factors has been reported to date. For the items incorporated in the prediction model reported so far, we selected items that were significant in predicting liver carcinogenesis in Japanese patients with hepatitis B and constructed a prediction model of liver carcinogenesis by the Cox proportional hazard model with the addition of Human Leukocyte Antigen (HLA) genotypes. The model, which included four items—sex, age at the time of examination, alpha-fetoprotein level (log10AFP) and presence or absence of HLA-A*33:03—revealed an area under the receiver operating characteristic curve (AUROC) of 0.862 for HCC prediction within 1 year and an AUROC of 0.863 within 3 years. A 1000 repeated validation test resulted in a C-index of 0.75 or higher, or sensitivity of 0.70 or higher, indicating that this predictive model can distinguish those at high risk of developing liver cancer within a few years with high accuracy. The prediction model constructed in this study, which can distinguish between chronic hepatitis B patients who develop hepatocellular carcinoma (HCC) early and those who develop HCC late or not, is clinically meaningful.

This study examined the association between developmental dysplasia of the hip (DDH) and disease-associated loci in a Japanese cohort. A genome-wide association study (GWAS) of 238 Japanese patients with DDH and 2044 healthy individuals was performed. As a replicate, GWAS was also conducted on the UK Biobank data with 3315 cases and matched 74,038 controls. Gene set enrichment analyses (GSEAs) of both the genetics and transcriptome of DDH were performed. Transcriptome analysis of cartilage specimens from DDH-associated osteoarthritis and femoral neck fractures was performed as a control. Most of the lead variants were very low-frequency ones in the UK, and variants in the Japanese GWAS could not be replicated with the UK GWAS. We assigned DDH-related candidate variants to 42 and 81 genes from the Japanese and UK GWASs, respectively, using functional mapping and annotation. GSEA of gene ontology, disease ontology, and canonical pathways identified the most enriched pathway to be the ferroptosis signaling pathway, both in the Japanese gene set as well as the Japanese and UK merged set. Transcriptome GSEA also identified significant downregulation of genes in the ferroptosis signaling pathway. Thus, the ferroptosis signaling pathway may be associated with the pathogenic mechanism of DDH.

Background: We established a community-based cohort study to assess the long-term impact of the Great East Japan Earthquake on disaster victims and gene-environment interactions on the incidence of major diseases, such as cancer and cardiovascular diseases. Methods: We asked participants to join our cohort in the health check-up settings and assessment center based settings. Inclusion criteria were aged 20 years or over and living in Miyagi or Iwate Prefecture. We obtained information on lifestyle, effect of disaster, blood, and urine information (Type 1 survey), and some detailed measurements (Type 2 survey), such as carotid echography and calcaneal ultrasound bone mineral density. All participants agreed to measure genome information and to distribute their information widely. Results: As a result, 87,865 gave their informed consent to join our study. Participation rate at health check-up site was about 70%. The participants in the Type 1 survey were more likely to have psychological distress than those in the Type 2 survey, and women were more likely to have psychological distress than men. Additionally, coastal residents were more likely to have higher degrees of psychological distress than inland residents, regardless of sex. Conclusion: This cohort comprised a large sample size and it contains information on the natural disaster, genome information, and metabolome information. This cohort also had several detailed measurements. Using this cohort enabled us to clarify the long-term effect of the disaster and also to establish personalized prevention based on genome, metabolome, and other omics information.

Background: Platinum derivatives are important treatment options for patients with esophageal carcinoma (EC), and a predictive marker for platinum-based therapy is needed for precision medicine. Patients and methods: This study contained two cohorts consisting of EC patients treated using platinum-based chemoradiation therapy (CRT) as the first-line and another external cohort of nationwide clinicogenomic data from the BioBank Japan (BBJ). Results: Genome-wide association study (GWAS) of therapeutic outcomes, refractory disease or not, following platinum-based CRT as first-line in 94 patients in the first cohort suggested the association of 89 SNPs using p

The conditions to induce human hepatic progenitor cells from other cell types are unclear. Here, the authors reprogram human endothelial cells to hepatic progenitor cells by expressing FOXA3, HNF1A and HNF6, capable of giving rise to hepatocytes and cholangiocytes that reconstitute damaged liver tissues on transplantation.

Abstract Background Genetic factors such as single-nucleotide polymorphisms (SNPs) play a key role in the development of postoperative nausea and vomiting (PONV). However, previous findings are not widely applicable to different populations because of population-specific genetic variation. We developed a Japanese-specific DNA microarray for high-throughput genotyping. The aim of the current study was to identify SNPs associated with PONV on a genome-wide scale using this microarray in a sample of Japanese surgical patients. Methods Associations between 659,636 SNPs and the incidence of PONV 24 h after surgery in a limited sample of 24 female patients were assessed using the microarray. After imputation of genotypes at 24,330,529 SNPs, 78 SNPs were found to be associated with the incidence of PONV. We chose 4 of the 78 SNPs to focus on by in silico functional annotation. Finally, we genotyped these 4 candidate SNPs in 255 patients using real-time PCR to verify association with the incidence of PONV. Results The T > C variant of rs11232965 in the long non-coding RNA MIR4300HG was significantly associated with reduced incidence of PONV among genotypes and between alleles (p = 0.01 and 0.007). Conclusions We identified a novel SNP (rs11232965) in the long non-coding RNA MIR4300HG that is associated with PONV. The rs11232965-SNP variant (T > C) is protective against the incidence of PONV. Trial registration This study was registered at the UMIN Clinical Trials Registry (Identifier: UMIN000022903 , date of registration: June 27, 2016, retrospectively registered.

Nonalcoholic fatty liver disease (NAFLD) is supposed to manifest its metabolic phenotype in the liver, but it is common to have lean individuals diagnosed with NAFLD, known as lean NAFLD. We conducted a two‐stage analysis to identify NAFLD‐associated loci in Japanese patients. In stage I, 275 metabolically healthy normal‐weight patients with NAFLD were compared with 1,411 non‐NAFLD controls adjusted for age, sex, and alcohol consumption by a genome‐wide association study (GWAS). In stage II, human leukocyte antigen (HLA) in chromosome 6 (chr6) (P = 6.73E‐08), microRNA (MIR) MIR548F3 in chr7 (P = 4.25E‐07), myosin light chain 2 (MYL2) in chr12 (P = 4.39E‐07), and glycoprotein precursor (GPC)6 in chr13 (P = 5.43E‐07), as suggested by the GWAS, were assessed by single nucleotide polymorphism (SNP) association analysis of whole NAFLD against non‐NAFLD in 9,726 members of the general population. A minor allele of the secondary lead SNP in chr6, rs2076529, was significantly associated (odds ratio [OR], 1.19; 95% confidence interval [CI], 1.11‐1.28; P = 2.10E‐06) and the lead SNP in chr7 was weakly associated (OR 1.15; 95% CI, 1.04‐1.27; P = 6.19E‐03) with increased NAFLD risk. Imputation‐based typing of HLA showed a significant difference in the distribution of HLA‐B, HLA‐DR‐beta chain 1 (DRB1), and HLA‐DQ‐beta chain 1 (DQB1) alleles in lean NAFLD GWAS. Next‐generation sequence‐based typing of HLA in 5,649 members of the general population replicated the significant difference of HLA‐B allele distribution and the significant increase of the HLA‐B*54:01 allele in whole NAFLD. Fecal metagenomic analysis of 3,420 members of the general population showed significant dissimilarity in beta‐diversity analysis of rs2076529 and HLA‐B*54:01 allele carriers from noncarriers. Veillonellaceae was increased but Verrucomicrobia was decreased in rs2076529 minor allele and HLA‐B*54:01 allele carriers as in NAFLD. Conclusion: HLA was identified as a novel locus associated with NAFLD susceptibility, which might be affected by the alteration of gut microbiota.

Over the past few decades, obesity has become a public health issue of global concern. Even though disparities exist between human populations, e.g., the higher liver fat content of the Japanese despite a lower body mass index (BMI), studies on the genetics of obesity still largely focus on populations of European descent, leading to a dearth of genetic data on non-European populations. In this context, this study aimed to establish a broad picture of the genetic attributes of the Japanese population, by examining a representative sample of 18,889 individuals participating in the Tohoku Medical Megabank Project cohort. We applied linear mixed model methods to 17 traits related to obesity and associated diseases to estimate the heritabilities explained by common genetic variants and the genetic correlations between each pair of traits. These analyses allowed us to quantify the SNP heritability of health indicators such as BMI (0.248 ± 0.032) and HDL cholesterol (0.324 ± 0.031), and to provide one of the few estimates of the SNP heritability of cystatin C in unrelated individuals (0.260 ± 0.025). We discuss potential differences between the Japanese and people of European ancestry with respect to the genetic correlations between urinary biomarkers and adiposity traits, for which large estimates were obtained. For instance, the genetic correlations between urine potassium level and the values for weight, BMI, waist circumference, and waist-to-height ratio ranged from 0.290 to 0.559, much higher than the corresponding estimates in the UK Biobank.

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 2 versus controls (P

Abstract Background Kawasaki disease (KD) is a systemic vasculitis which may be associated with coronary artery aneurysms. A notable risk factor for the development of coronary artery aneurysms is resistance to intravenous immunoglobulin (IVIG) therapy, which comprises standard treatment for the acute phase of KD. The cause of IVIG resistance in KD is largely unknown; however, the contribution of genetic factors, especially variants in immune-related genes, has been suspected. Methods To explore genetic variants related to IVIG-unresponsiveness, we designated KD patients who did not respond to both first and second courses of IVIG therapy as IVIG-unresponsive patients. Using genomic DNA from 30 IVIG-unresponsive KD patients, we performed pooled genome sequencing targeting 39 immune-related cytokine receptor genes. Results The single nucleotide variant (SNV), rs563535954 (located in the IL4R locus), was concentrated in IVIG-unresponsive KD patients. Individual genotyping showed that the minor allele of rs563535954 was present in 4/33 patients with IVIG-unresponsive KD, compared with 20/1063 individuals in the Japanese genome variation database (odds ratio = 7.19, 95% confidence interval 2.43–21.47). Furthermore, the minor allele of rs563535954 was absent in 42 KD patients who responded to IVIG treatment (P = 0.0337), indicating that a low-frequency variant, rs563535954, is associated with IVIG-unresponsiveness in KD patients. Although rs563535954 is located in the 3′-untranslated region of IL4R, there was no alternation in IL4R expression associated with the mior allele of rs563535954. However, IVIG-unresponsive patients that exhibited the minor allele of rs563535954 tended to be classified into the low-risk group (based on previously reported risk scores) for prediction of IVIG-resistance. Therefore, IVIG-unresponsiveness associated with the minor allele of rs563535954 might differ from IVIG-unresponsiveness associated with previous risk factors used to evaluate IVIG-unresponsiveness in KD. Conclusion These findings suggest that the SNV rs563535954 could serve as a predictive indicator of IVIG-unresponsiveness, thereby improving the sensitivity of risk scoring systems, and may aid in prevention of coronary artery lesions in KD patients.

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. Keywords: PBC, GWAS, Disease Susceptibility Gene, e-QTL, CRISPR/Cas9

Background Although chronic thromboembolic pulmonary hypertension (CTEPH) and acute pulmonary embolism (APE) share some clinical manifestations, a limited proportion of patients with CTEPH have a history of APE. Moreover, in histopathologic studies, it has been revealed that pulmonary vasculature lesions similar to pulmonary arterial hypertension existed in patients with CTEPH. Thus, it remains unknown whether these 3 disorders also share genetic backgrounds. Methods and Results Whole exome screening was performed with DNA isolated from 51 unrelated patients with CTEPH of Japanese ancestry. The frequency of genetic variants associated with pulmonary arterial hypertension or APE in patients with CTEPH was compared with those in the integrative Japanese Genome Variation Database 3.5KJPN. Whole exome screening analysis showed 17 049 nonsynonymous variants in patients with CTEPH. Although we found 6 nonsynonymous variants that are associated with APE in patients with CTEPH, there was no nonsynonymous variant associated with pulmonary arterial hypertension. Patients with CTEPH with a history of APE had nonsynonymous variants of F5, which encodes factor V. In contrast, patients with CTEPH without a history of APE had a nonsynonymous variant of THBD, which encodes thrombomodulin. Moreover, thrombin‐activatable fibrinolysis inhibitor, which is one of the pathogenic proteins in CTEPH, was significantly more activated in those who had the variants of THBD compared with those without it. Conclusions These results provide the first evidence that patients with CTEPH have some variants associated with APE, regardless of the presence or absence of a history of APE. Furthermore, the variants might be different between patients with CTEPH with and without a history of APE.

Abstract The associations of single nucleotide polymorphisms (SNPs) in PLA2R1 and HLA-DQA1, as well as HLA-DRB1*15:01-DQB1*06:02 haplotype with idiopathic membranous nephropathy (IMN) is well known. However, the primary associations of these loci still need to be determined. We used Japanese-specific SNP genotyping array and imputation using 2,048 sequenced Japanese samples to fine-map PLA2R1 region in 98 patients and 413 controls. The most significant SNPs were replicated in a separate sample set of 130 patients and 288 controls. A two-SNP haplotype of intronic and missense SNPs showed the strongest association. The intronic SNP is strongly associated with PLA2R1 expression in the Genotype-Tissue Expression (GTEx) database, and the missense SNP is predicted to alter peptide binding with HLA-DRB1*15:01 by the Immune Epitope Database (IEDB). In HLA region, we performed relative predispositional effect (RPE) tests and identified additional risk alleles in both HLA-DRB1 and HLA-DQB1. We collapsed the risk alleles in each of HLA-DRB1 and HLA-DQB1 into single risk alleles. Reciprocal conditioning of these collapsed risk alleles showed more residual significance for HLA-DRB1 collapsed risk than HLA-DQB1 collapsed risk. These results indicate that changes in the expression levels of structurally different PLA2R protein confer risk for IMN in the presence of risk HLA-DRB1 alleles.

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

Abstract Numerous genome-wide association studies (GWAS) have been performed to identify susceptibility genes to various human complex diseases. However, in many cases, neither a functional variant nor a disease susceptibility gene have been clarified. Here, we show an efficient approach for identification of a functional variant in a primary biliary cholangitis (PBC)-susceptible region, chromosome 17q12-21 (ORMDL3-GSDMB-ZPBP2-IKZF3). High-density association mapping was carried out based on SNP imputation analysis by using the whole-genome sequence data from a reference panel of 1,070 Japanese individuals (1KJPN), together with genotype data from our previous GWAS (PBC patients: n = 1,389; healthy controls: n = 1,508). Among 23 single nucleotide polymorphisms (SNPs) with P

Abstract Retinal neurons and Müller glia are generated from a common population of multipotent retinal progenitor cells. We purposed to identify Müller glia-specific molecular signatures during retinal development. Using transgenic mice carrying the Hes1 promoter (pHes1) followed by EGFP, we purified EGFP-positive Müller glia and other EGFP-negative retinal cells from developing retinas and subjected them to RNA sequencing analysis. Gene expression pattern of EGFP-positive cell was similar to genes expressed in retinal progenitors, and they were downregulated in other cell lineages. Then, we examined the modification profiles of H3K27me3 and H3K4me3 by referring to chromatin immunoprecipitation-sequencing data of rods and other cells. Clustering of the H3K4me3 and H3K27me3 values followed by ontology analysis revealed a high incidence of transcription factors including Hes1 in clusters with high H3K27me3 levels. Hes1 expression level decreased dramatically, and the H3K27me3 level at the Hes1-locus was upregulated strongly during retinal development. Furthermore, the Hes1 expression level was upregulated in an Ezh2-knockout retina. These results suggest that downregulation of Müller glia-related genes in other lineage rather than upregulation of them in Müller glia contributed Müller-specific molecular features, and a role for modified H3K27me3 in suppressing Müller glia-related genes in other retinal cell lineages to avoid unfavorable expression.

Epigenomics: Finding DNA modifications that predict disease Patterns of chemical modifications on DNA can predict the risk of certain diseases, but the challenge is knowing where to look for them. Atsushi Shimizu from Iwate Medical University in Japan and his colleagues determined the location of these modifications, methyl groups added to DNA’s cytosine bases, in over hundred healthy people and found two million regions that vary widely between individuals. They used previous data to show that these regions are far more likely to change during the course of disease than genomic loci that are unchanged in all subjects. Current methods to look for disease-specific methylation changes can only profile at most 20% of the genome and Shimizu’s work will help focus on the 20% that matter for identifying disease risk and lead to better diagnosis and prognosis.

[This corrects the article DOI: 10.1371/journal.pone.0194036.].

Forensic DNA typing is widely used to identify missing persons and plays a central role in forensic profiling. DNA typing usually uses capillary electrophoresis fragment analysis of PCR amplification products to detect the length of short tandem repeat (STR) markers. Here, we analyzed whole genome data from 1,070 Japanese individuals generated using massively parallel short-read sequencing of 162 paired-end bases. We have analyzed 843,473 STR loci with two to six basepair repeat units and cataloged highly polymorphic STR loci in the Japanese population. To evaluate the performance of the cataloged STR loci, we compared 23 STR loci, widely used in forensic DNA typing, with capillary electrophoresis based STR genotyping results in the Japanese population. Seventeen loci had high correlations and high call rates. The other six loci had low call rates or low correlations due to either the limitations of short-read sequencing technology, the bioinformatics tool used, or the complexity of repeat patterns. With these analyses, we have also purified the suitable 218 STR loci with four basepair repeat units and 53 loci with five basepair repeat units both for short read sequencing and PCR based technologies, which would be candidates to the actual forensic DNA typing in Japanese population.

BACKGROUND:Inflammatory bowel disease (IBD) has an unknown etiology; however, accumulating evidence suggests that IBD is a multifactorial disease influenced by a combination of genetic and environmental factors. The influence of genetic variants on DNA methylation in cis and cis effects on expression have been demonstrated. We hypothesized that IBD susceptibility single-nucleotide polymorphisms (SNPs) regulate susceptibility gene expressions in cis by regulating DNA methylation around SNPs. For this, we determined cis-regulated allele-specific DNA methylation (ASM) around IBD susceptibility genes in CD4+ effector/memory T cells (Tem) in lamina propria mononuclear cells (LPMCs) in patients with IBD and examined the association between the ASM SNP genotype and neighboring susceptibility gene expressions. METHODS:CD4+ effector/memory T cells (Tem) were isolated from LPMCs in 15 Japanese IBD patients (ten Crohn's disease [CD] and five ulcerative colitis [UC] patients). ASM analysis was performed by methylation-sensitive SNP array analysis. We defined ASM as a changing average relative allele score ([Formula: see text]) >0.1 after digestion by methylation-sensitive restriction enzymes. Among SNPs showing [Formula: see text] >0.1, we extracted the probes located on tag-SNPs of 200 IBD susceptibility loci and around IBD susceptibility genes as candidate ASM SNPs. To validate ASM, bisulfite-pyrosequencing was performed. Transcriptome analysis was examined in 11 IBD patients (seven CD and four UC patients). The relation between rs36221701 genotype and neighboring gene expressions were analyzed. RESULTS:We extracted six candidate ASM SNPs around IBD susceptibility genes. The top of [Formula: see text] (0.23) was rs1130368 located on HLA-DQB1. ASM around rs36221701 ([Formula: see text] = 0.14) located near SMAD3 was validated using bisulfite pyrosequencing. The SMAD3 expression was significantly associated with the rs36221701 genotype (p = 0.016). CONCLUSIONS:We confirmed the existence of cis-regulated ASM around IBD susceptibility genes and the association between ASM SNP (rs36221701) genotype and SMAD3 expression, a susceptibility gene for IBD. These results give us supporting evidence that DNA methylation mediates genetic effects on disease susceptibility.

To test the genetic association between Japanese patients with primary open-angle glaucoma (POAG) and the previously reported POAG susceptibility loci and to perform genotype-phenotype analysis.Genetic associations for 27 SNPs from 16 loci previously linked to POAG were assessed using genome-wide SNP data of the primary cohort (565 Japanese POAG patients and 1,104 controls). Reproducibility of the assessment was tested in 607 POAG cases and 455 controls (second cohort) with a targeted genotyping approach. For POAG-associated variants, a genotype-phenotype correlation study (additive, dominant, recessive model) was performed using the objective clinical data derived from 598 eyes of 598 POAG patients.Among 27 SNPs from 16 loci previously linked to POAG, genotypes for total of 20 SNPs in 13 loci were available for targeted association study. Among 8 SNPs in 3 loci that showed at least nominal association (P < 5.00E-02) in the primary cohort, a representative SNP for each loci (rs2157719 for CDKN2B-AS1, rs33912345 for SIX6, and rs9913911 for GAS7) were selected. For these SNPs the association was found significant in both the second cohort analysis and meta-analysis. The genotype-phenotype analysis revealed significant correlations between CDKN2B-AS1 (rs2157719) and decreased intraocular pressure (β = -6.89 mmHg, P = 1.70E-04; dominant model) after multiple corrections. In addition, nominal correlation was observed between CDKN2B-AS1 (rs2157719) and optic nerve head blood flow (β = -0.54 and -0.67 arbitrary units (AU), P = 2.00E-02 and 1.39E-02), between SIX6 (rs33912345) and decreased total peripapillary retinal nerve fiber layer thickness (β = -2.16 and -2.82 μm, P = 4.68E-02 and 2.40E-02, additive and recessive model, respectively) and increased optic nerve head blood flow (β = 0.44 AU, P = 2.20E-02; additive model) and between GAS7 (rs9913911) and increased cup volume (β = 0.03 mm3, P = 4.60E-02) and mean cup depth (β = 0.03 mm3, P = 4.11E-02; additive model) and decreased pattern standard deviation (β = -0.87 dB, P = 2.44E-02; dominant model).The association between SNPs near GAS7 and POAG was found in Japanese patients for the first time. Clinical characterization of the risk variants is an important step toward understanding the pathology of the disease and optimizing treatment of patients with POAG.