Search

Recent advances have enabled the direct induction of human tissue-specific stem and progenitor cells from differentiated somatic cells. However, it is not known whether human hepatic progenitor cells (hHepPCs) can be generated from other cell types by direct lineage reprogramming with defined transcription factors. Here, we show that a set of three transcription factors, FOXA3, HNF1A, and HNF6, can induce human umbilical vein endothelial cells to directly acquire the properties of hHepPCs. These induced hHepPCs (hiHepPCs) propagate in long-term monolayer culture and differentiate into functional hepatocytes and cholangiocytes by forming cell aggregates and cystic epithelial spheroids, respectively, under three-dimensional culture conditions. After transplantation, hiHepPC-derived hepatocytes and cholangiocytes reconstitute damaged liver tissues and support hepatic function. The defined transcription factors also induce hiHepPCs from endothelial cells circulating in adult human peripheral blood. These expandable and bipotential hiHepPCs may be useful in the study and treatment of human liver diseases., 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.

Purpose: The epithelial-to-mesenchymal transition, the major process by which some cancer cells convert from an epithelial phenotype to a mesenchymal one, has been suggested to drive chemo-resistance and/or metastasis in patients with cancer. However, only a few studies have demonstrated the presence of CD45/CD326 doubly-positive cells (CD45/CD326 DPC) in cancer. We deployed a combination of cell surface markers to elucidate the phenotypic heterogeneity in non–small cell lung cancer (NSCLC) cells and identified a new subpopulation that is doubly-positive for epithelial and non–epithelial cell-surface markers in both NSCLC cells and patients' malignant pleural effusions. Experimental Design: We procured a total of 39 patients' samples, solid fresh lung cancer tissues from 21 patients and malignant pleural effusion samples from 18 others, and used FACS and fluorescence microscopy to check their surface markers. We also examined the EGFR mutations in patients with known acquired EGFR mutations. Results: Our data revealed that 0.4% to 17.9% of the solid tumor tissue cells and a higher percentage of malignant pleural effusion cells harbored CD45/CD326 DPC expressing both epithelial and nonepithelial surface markers. We selected 3 EGFR mutation patients and genetically confirmed that the newly identified cell population really originated from cancer cells. We also found that higher proportions of CD45/CD326 DPC are significantly associated with poor prognosis. Conclusions: In conclusion, varying percentages of CD45/CD326 DPC exist in both solid cancer tissue and malignant pleural effusion in patients with NSCLC. This CD45/CD326 doubly-positive subpopulation can be an important key to clinical management of patients with NSCLC.

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., Graphical abstract

In time-series analysis, state-space models (SSMs) have been widely used to estimate the conditional probability distributions of hidden variables and parameter values, as well as to understand structures that can generate the data. For example, the Kalman filter is used to analytically calculate the conditional probability distribution on linear SSMs in terms of minimizing the variance, and several extensions, such as the unscented Kalman filter and particle filter, have been applied to calculate the approximate distribution on nonlinear SSMs. Recently, the approximate Bayesian computation (ABC) has been applied to such time-series filtering to handle intractable likelihoods; however, it remains problematic with respect to consistently achieving a reduction of the estimation bias, evaluating the validity of the models, and dealing with replicated observations. To address these problems, in this paper, we propose a novel method combined with the kernel ABC to perform filtering, parameter estimation, and the model evaluation in SSMs. Simulation studies show that the proposed method produces inference that is comparable to other ABC methods, with the advantage of not requiring a careful calibration of the ABC threshold. In addition, we evaluate the performance of the model-selection capability using true and competitive models on synthetic data from nonlinear SSMs. Finally, we apply the proposed method to real data in rat circadian oscillations, and demonstrated the usefulness in practical situations.

Genetic polymorphisms contribute to inter-individual variability in the metabolism of multiple clinical drugs, including warfarin, thiopurines, primaquine, and aminoglycosides. A rapid and sensitive clinical assessment of various genome biomarkers is, therefore, required to predict the individual responsiveness of each patient to these drugs. In this study, we developed a novel genotyping method for the detection of nine pharmacogene variants that are important in the prediction of drug efficiency and toxicity. This genotyping method uses competitive allele-specific PCR and a single-stranded tag hybridization chromatographic printed-array strip (STH-PAS) that can unambiguously determine the presence or absence of the gene variant by displaying visible blue lines on the chromatographic printed-array strip. Notably, the results of our STH-PAS method were in 100% agreement with those obtained using standard Sanger sequencing and KASP assay genotyping methods for CYP4F2 gene deletion. Moreover, the results were obtained within 90 min, including the PCR amplification and signal detection processes. The sensitive and rapid nature of this novel method make it ideal for clinical genetic testing to predict drug efficacy and toxicity, and in doing so will aid in the development of individualized medicine and better patient care.

Background & aims: Genome-wide association studies in primary biliary cholangitis (PBC) have failed to find X chromosome (chrX) variants associated with the disease. Here, we specifically explore the chrX contribution to PBC, a sexually dimorphic complex autoimmune disease. Methods: We performed a chrX-wide association study, including genotype data from 5 genome-wide association studies (from Italy, United Kingdom, Canada, China, and Japan; 5244 case patients and 11,875 control individuals). Results: Single-marker association analyses found approximately 100 loci displaying P < 5 × 10-4, with the most significant being a signal within the OTUD5 gene (rs3027490; P = 4.80 × 10-6; odds ratio [OR], 1.39; 95% confidence interval [CI], 1.028-1.88; Japanese cohort). Although the transethnic meta-analysis evidenced only a suggestive signal (rs2239452, mapping within the PIM2 gene; OR, 1.17; 95% CI, 1.09-1.26; P = 9.93 × 10-8), the population-specific meta-analysis showed a genome-wide significant locus in East Asian individuals pointing to the same region (rs7059064, mapping within the GRIPAP1 gene; P = 6.2 × 10-9; OR, 1.33; 95% CI, 1.21-1.46). Indeed, rs7059064 tags a unique linkage disequilibrium block including 7 genes: TIMM17B, PQBP1, PIM2, SLC35A2, OTUD5, KCND1, and GRIPAP1, as well as a superenhancer (GH0XJ048933 within OTUD5) targeting all these genes. GH0XJ048933 is also predicted to target FOXP3, the main T-regulatory cell lineage specification factor. Consistently, OTUD5 and FOXP3 RNA levels were up-regulated in PBC case patients (1.75- and 1.64-fold, respectively). Conclusions: This work represents the first comprehensive study, to our knowledge, of the chrX contribution to the genetics of an autoimmune liver disease and shows a novel PBC-related genome-wide significant locus.

[BACKGROUND & AIMS] Primary biliary cholangitis (PBC) is a chronic liver disease in which autoimmune destruction of the small intra-hepatic bile ducts eventually leads to cirrhosis. Many patients have inadequate response to licensed medications, motivating the search for novel therapies. Previous genome-wide association studies (GWAS) and meta-analyses (GWMA) of PBC have identified numerous risk loci for this condition, providing insight into its aetiology. We undertook the largest GWMA of PBC to date, aiming to identify additional risk loci and prioritise candidate genes for in silico drug efficacy screening. [METHODS] We combined new and existing genotype data for 10, 516 cases and 20, 772 controls from five European and two East Asian cohorts. [RESULTS] We identified 56 genome-wide significant loci (20 novel) including 46 in European, 13 in Asian, and 41 in combined cohorts; and a 57th genome-wide significant locus (also novel) in conditional analysis of the European cohorts. Candidate genes at newly identified loci include FCRL3, INAVA, PRDM1, IRF7, CCR6, CD226, and IL12RB1, each having key roles in immunity. Pathway analysis reiterated the likely importance of pattern recognition receptor and TNF signalling, Jak-STAT signalling, and differentiation of TH1 and TH17 cells in the pathogenesis of this disease. Drug efficacy screening identified several medications predicted to be therapeutic in PBC, some well-established in the treatment of other autoimmune disorders. [CONCLUSIONS] This study has identified additional risk loci for PBC, provided a hierarchy of agents that could be trialled in this condition, and emphasised the value of genetic and genomic approaches to drug discovery in complex disorders. [Lay summary] Primary biliary cholangitis (PBC) is a chronic liver disease that eventually leads to cirrhosis. In this study, we analysed genetic information from 10, 516 people with PBC and 20, 772 healthy individuals recruited in Canada, China, Italy, Japan, UK, or USA. We identified several genetic regions associated with PBC. Each of these regions contains several genes. For each region, we used diverse sources of evidence to help us choose the gene most likely to be involved in causing PBC. We used these ‘candidate genes’ to help us identify medications that are currently used for treatment of other conditions, which might also be useful for treatment of PBC., 原発性胆汁性胆管炎のゲノムワイド関連解析 --国際メタ解析による新規疾患感受性遺伝子と治療薬候補の同定--. 京都大学プレスリリース. 2021-06-28.

Rs671 in the aldehyde dehydrogenase 2 gene (ALDH2) is the cause of Asian alcohol flushing response after drinking. ALDH2 detoxifies endogenous aldehydes, which are the major source of DNA damage re...

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.

Amyotrophic lateral sclerosis (ALS) is a devastating progressive motor neuron disease that affects people of all ethnicities. Approximately 90% of ALS cases are sporadic and thought to have multifactorial pathogenesis. To understand the genetics of sporadic ALS, we conducted a genome-wide association study using 1,173 sporadic ALS cases and 8,925 controls in a Japanese population. A combined meta-analysis of our Japanese cohort with individuals of European ancestry revealed a significant association at the ACSL5 locus (top SNP p = 2.97 × 10−8). We validated the association with ACSL5 in a replication study with a Chinese population and an independent Japanese population (1941 ALS cases, 3821 controls; top SNP p = 1.82 × 10−4). In the combined meta-analysis, the intronic ACSL5 SNP rs3736947 showed the strongest association (p = 7.81 × 10−11). Using a gene-based analysis of the full multi-ethnic dataset, we uncovered additional genes significantly associated with ALS: ERGIC1, RAPGEF5, FNBP1, and ATXN3. These results advance our understanding of the genetic basis of sporadic ALS., Gen Sobue, Masao Nagasaki and colleagues report a genome-wide association study for amyotrophic lateral sclerosis (ALS) in a large, multi-ethnic cohort comprising Japanese, Chinese, and European ancestry populations. They find a significant association to variants within the ACSL5 gene and identify novel associations with 4 additional genes using a gene-based approach.

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.

Panel sequencing of susceptibility genes for hereditary breast and ovarian cancer (HBOC) syndrome has uncovered numerous germline variants; however, their pathogenic relevance and ethnic diversity remain unclear. Here, we examined the prevalence of germline variants among 568 Japanese patients with BRCA1/2-wildtype HBOC syndrome and a strong family history. Pathogenic or likely pathogenic variants were identified on 12 causal genes for 37 cases (6.5%), with recurrence for 4 SNVs/indels and 1 CNV. Comparisons with non-cancer east-Asian populations and European familial breast cancer cohorts revealed significant enrichment of PALB2, BARD1, and BLM mutations. Younger onset was associated with but not predictive of these mutations. Significant somatic loss-of-function alterations were confirmed on the wildtype alleles of genes with germline mutations, including PALB2 additional somatic truncations. This study highlights Japanese-associated germline mutations among patients with BRCA1/2 wildtype HBOC syndrome and a strong family history, and provides evidence for the medical care of this high-risk population.

To understand the genetics of steroid-sensitive nephrotic syndrome (SSNS), we conducted a genome-wide association study in 987 childhood SSNS patients and 3,206 healthy controls with Japanese ancestry. Beyond known associations in the HLA-DR/DQ region, common variants in NPHS1-KIRREL2 (rs56117924, P=4.94E-20, odds ratio (OR) =1.90) and TNFSF15 (rs6478109, P=2.54E-8, OR=0.72) regions achieved genome-wide significance and were replicated in Korean, South Asian and African populations. Trans-ethnic meta-analyses including Japanese, Korean, South Asian, African, European, Hispanic and Maghrebian populations confirmed the significant associations of variants in NPHS1-KIRREL2 (P(meta)=6.71E-28, OR=1.88) and TNFSF15 (P(meta)=5.40E-11, OR=1.33) loci. Analysis of the NPHS1 risk alleles with glomerular NPHS1 mRNA expression from the same person revealed allele specific expression with significantly lower expression of the transcript derived from the risk haplotype (Wilcox test p=9.3E-4). Because rare pathogenic variants in NPHS1 cause congenital nephrotic syndrome of the Finnish type (CNSF), the present study provides further evidence that variation along the allele frequency spectrum in the same gene can cause or contribute to both a rare monogenic disease (CNSF) and a more complex, polygenic disease (SSNS).

Primary open-angle glaucoma (POAG) can develop even within normal ranges of intraocular pressure, and this type of glaucoma (so-called ‘normal-tension glaucoma [NTG]’) is highly prevalent in East Asia including Korea and Japan. We conducted exome chip analysis to identify low-frequency and rare variants associated with POAG from the primary cohort (309 POAG patients and 5,400 control, all Koreans). For replication, Korean (310 POAG patients and 5,612 controls) and Japanese (565 POAG patients and 1,104 controls) cohorts were further investigated by targeted genotyping. SNP rs116121322 in LRRC27 showed nominally significant association with POAG in the discovery cohort (OR = 29.85, P = 2E–06). This SNP was validated in the Korean replication cohort but only in the NTG subgroups (OR = 9.86, P = 0.007). Japanese replication cohort did not show significant association with POAG (P .00.44). However, the meta-analysis in the entire cohort revealed significant association of rs116121322 with POAG (ORcombined = 10.28, Pcombined = 1.4E–07). The LRRC27 protein expression was confirmed from human trabecular meshwork cells. For gene-based testing, METTL20 showed a significant association in POAG (Pcombined = 0.002) and in the subgroup of NTG (Pcombined = 0.02), whereas ZNF677 were significantly associated with only in the subgroup of high-tension glaucoma (Pcombined = 1.5E–06). Our findings may provide further genetic backgrounds into the pathogenesis of POAG, especially for the patients who have lower baseline intraocular pressures.

Genetic variations within cytochrome P450 2B6 (CYP2B6) contribute to inter-individual variation in the metabolism of clinically important drugs, including cyclophosphamide, bupropion, methadone and efavirenz (EFZ). In this study, we performed an in vitro analysis of 40 CYP2B6 allelic variant proteins including seven novel variants identified in 1070 Japanese individuals. Wild-type and 39 variant proteins were heterologously expressed in 293FT cells to estimate the kinetic parameters (Km, Vmax, and CLint) of EFZ 8-hydroxylation and 7-ethoxy-4-trifluoromethylcoumarin (7-ETC) O-deethylation activities. The concentrations of CYP2B6 variant holo-enzymes were measured by using carbon monoxide (CO)-reduced difference spectroscopy, and the wild-type and 28 variants showed a peak at 450 nm. The kinetic parameters were measured for the wild-type and 24 variant proteins. The values for the remaining 15 variants could not be determined because the enzymatic activity was not detected at the highest substrate concentration used. Compared to wild-type, six variants showed significantly decreased EFZ 8-hydroxylation CLint values, while these values were significantly increased in another six variants, including CYP2B6.6. Although 7-ETC O-deethylation CLint values of CYP2B6 variants did not differ significantly from that of CYP2B6.1, the CLint ratios obtained for 7-ETC O-deethylation were highly correlated with EFZ 8-hydroxylation. Furthermore, three-dimensional structural modeling analysis was performed to elucidate the mechanism of changes in the kinetics of CYP2B6 variants. Our findings could provide evidence of the specific metabolic activities of the CYP2B6 proteins encoded by these variant alleles.

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 Pallelic=7.84×10-23; odds ratio [OR], 0.33; 95% confidence interval [95% CI], 0.26 to 0.41; rs3134996 [minor allele A], combined Pallelic=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.

Background Despite NUDT15 variants showing significant association with thiopurine-induced adverse events (AEs) in Asians, it remains unclear which variants of NUDT15 or whether additional genetic variants should be tested to predict AEs. To clarify the best pharmacogenetic test to be used clinically, we performed association studies of NUDT15 variants and haplotypes with AEs, genome-wide association study (GWAS) to discover additional variants, and ROC analysis to select the model to predict severe AEs. Methods Overall, 2630 patients with inflammatory bowel disease (IBD) were enrolled and genotyped for NUDT15 codon 139; 1291 patients were treated with thiopurines. diplotypes were analyzed in 970 patients, and GWASs of AEs were performed with 1221 patients using population-optimized genotyping array and imputation. Results We confirmed the association of NUDT15 p.Arg139Cys with leukopenia and alopecia (p = 2.20E−63, 1.32E−69, OR = 6.59, 12.1, respectively), and found a novel association with digestive symptoms (p = 6.39E−04, OR = 1.89). Time to leukopenia was significantly shorter, and when leukopenia was diagnosed, thiopurine doses were significantly lower in Arg/Cys and Cys/Cys than in Arg/Arg. In GWASs, no additional variants were found to be associated with thiopurine-induced AEs. Despite strong correlation of leukopenia frequency with estimated enzyme activities based on the diplotypes (r2 = 0.926, p = 0.0087), there were no significant differences in the AUCs of diplotypes from those of codon 139 to predict severe AEs (AUC = 0.916, 0.921, for acute severe leukopenia, AUC = 0.990, 0.991, for severe alopecia, respectively). Conclusions Genotyping of NUDT15 codon 139 was sufficient to predict acute severe leukopenia and alopecia in Japanese patients with IBD. Electronic supplementary material The online version of this article (10.1007/s00535-018-1486-7) contains supplementary material, which is available to authorized users.

Dihydropyrimidine dehydrogenase (DPD, EC 1.3.1.2), encoded by the DPYD gene, is the rate-limiting enzyme in the degradation pathway of endogenous pyrimidine and fluoropyrimidine drugs such as 5-fluorouracil (5-FU). DPD catalyzes the reduction of uracil, thymine, and 5-FU. In Caucasians, DPYD mutations, including DPYD*2A, DPYD*13, c.2846A>T, and c.1129-5923C>G/hapB3, are known to contribute to interindividual variations in the toxicity of 5-FU; however, none of these DPYD polymorphisms has been identified in the Asian population. Recently, 21 DPYD allelic variants, including some novel single-nucleotide variants (SNVs), were identified in 1070 healthy Japanese individuals by analyzing their whole-genome sequences (WGSs), but the functional alterations caused by these variants remain unknown. In this study, in vitro analysis was performed on 22 DPD allelic variants by transiently expressing wild-type DPD and 21 DPD variants in 293FT cells and characterizing their enzymatic activities using 5-FU as a substrate. DPD expression levels and dimeric forms were determined using immunoblotting and blue-native PAGE, respectively. Additionally, the values of three kinetic parameters-the Michaelis constant (Km ), maximum velocity (Vmax ), and intrinsic clearance (CLint = Vmax/Km )-were determined for the reduction of 5-FU. Eleven variants exhibited significantly decreased intrinsic clearance compared with wild-type DPD. Moreover, the band patterns observed in the immunoblots of blue-native gels indicated that DPD dimerization is required for enzymatic activity in DPD. Thus, the detection of rare DPYD variants might facilitate severe adverse effect prediction of 5-FU-based chemotherapy in the Japanese population.

Background and aim Patients with inflammatory bowel disease (IBD) are at a high risk of low bone mineral density (BMD). Reportedly, clinical and genetic factors cause low BMD in Caucasians; however, studies in non-Caucasian populations remain scarce. Methods Clinical risk factors for low BMD were investigated in 266 Japanese patients with IBD, and a genome-wide association analysis (GWAS) was performed using linear regression with associated clinical factors as covariates. Genotyping was performed using a population-optimized genotyping array (Japonica array® ). After quality control, the genotype data of 4 384 682 single-nucleotide polymorphisms (SNPs) from 254 patients with IBD were used for GWAS. Results Body mass index, age, and disease duration were independently associated with the BMD of the femoral neck (P = 1.41E - 13, 1.04E - 5, and 1.58E - 3, respectively), and body mass index and sex were associated with the BMD of the lumbar spine (P = 6.90E - 10 and 6.84E - 3, respectively). In GWAS, 118 and 42 candidate SNPs of the femoral neck and lumbar spine, respectively, were identified. Among 118, 111 candidate SNPs of the femoral neck were located within the SLC22A23 gene, which is a known IBD susceptibility gene (minimum P = 1.42E - 07). Among 42, 18 candidate SNPs of the lumbar spine were located within the MECOM gene, which is associated with osteopenia (minimum P = 5.86E - 07). Interestingly, none of the known loci showed a significant association with BMD. Conclusions Although clinical risk factors for low BMD in IBD were similar to those in the general population, genetic risk factors were rather different.

Approximately 10–20% of patients with primary biliary cholangitis (PBC) progress to jaundice stage regardless of treatment with ursodeoxycholic acid and bezafibrate. In this study, we performed a GWAS and a replication study to identify genetic variants associated with jaundice-stage progression in PBC using a total of 1,375 patients (1,202 early-stage and 173 jaundice-stage) in a Japanese population. SNP rs13720, which is located in the 3′UTR of cathepsin Z (CTSZ), showed the strongest association (odds ratio [OR] = 2.15, P = 7.62 × 10−7) with progression to jaundice stage in GWAS. High-density association mapping at the CTSZ and negative elongation factor complex member C/D (NELFCD) loci, which are located within a strong linkage disequilibrium (LD) block, revealed that an intronic SNP of CTSZ, rs163800, was significantly associated with jaundice-stage progression (OR = 2.16, P = 8.57 × 10−8). In addition, eQTL analysis and in silico functional analysis indicated that genotypes of rs163800 or variants in strong LD with rs163800 influence expression levels of both NELFCD and CTSZ mRNA. The present novel findings will contribute to dissect the mechanism of PBC progression and also to facilitate the development of therapies for PBC patients who are resistant to current therapies.

Aim Somatic mutations in the human brain are hypothesized to contribute to the functional diversity of brain cells as well as the pathophysiology of neuropsychiatric diseases. However, there are still few reports on somatic mutations in non-neoplastic human brain tissues. This study attempted to unveil the landscape of somatic mutations in the human brain. Methods We explored the landscape of somatic mutations in human brain tissues derived from three individuals with no neuropsychiatric diseases by whole-genome deep sequencing at a depth of around 100. The candidate mutations underwent multi-layered filtering, and were validated by ultra-deep target amplicon sequencing at a depth of around 200 000. Results Thirty-one somatic mutations were identified in the human brain, demonstrating the utility of whole-genome sequencing of bulk brain tissue. The mutations were enriched in neuron-expressed genes, and two-thirds of the identified somatic single nucleotide variants in the brain tissues were cytosine-to-thymine transitions, half of which were in CpG dinucleotides. Conclusion Our developed filtering and validation approaches will be useful to identify somatic mutations in the human brain. The vulnerability of neuron-expressed genes to mutational events suggests their potential relevance to neuropsychiatric diseases.

Clarifying allele frequencies of disease-related genetic variants in a population is important in genomic medicine; however, such data is not yet available for the Japanese population. To estimate frequencies of actionable pathogenic variants in the Japanese population, we examined the reported pathological variants in genes recommended by the American College of Medical Genetics and Genomics (ACMG) in our reference panel of genomic variations, 2KJPN, which was created by whole-genome sequencing of 2049 individuals of the resident cohort of the Tohoku Medical Megabank Project. We searched for pathogenic variants in 2KJPN for 57 autosomal ACMG-recommended genes responsible for 26 diseases and then examined their frequencies. By referring to public databases of pathogenic variations, we identified 143 reported pathogenic variants in 2KJPN for the 57 ACMG recommended genes based on a classification system. At the individual level, 21% of the individuals were found to have at least one reported pathogenic allele. We then conducted a literature survey to review the variants and to check for evidence of pathogenicity. Our results suggest that a substantial number of people have reported pathogenic alleles for the ACMG genes, and reviewing variants is indispensable for constructing the information infrastructure of genomic medicine for the Japanese population.

Ovarian clear cell carcinoma (OCCC) is the most refractory subtype of ovarian cancer and more prevalent in Japanese than Caucasians (25% and 5% of all ovarian cancer, respectively). The aim of this study is to discover the genomic alterations that may cause OCCC and effective molecular targets for chemotherapy. Paired genomic DNAs of 48 OCCC tissues and corresponding noncancerous tissues were extracted from formalin-fixed, paraffin embedded specimens collected between 2007 and 2015 at Tohoku University Hospital. All specimens underwent exome sequencing and the somatic genetic alterations were identified. We divided the cases into three clusters based on the mutation spectra. Clinical characteristics such as age of onset and endometriosis are similar among the clusters but one cluster shows mutations related to APOBEC activation, indicating its contribution to subset of OCCC cases. There are three hypermutated cases (showing 12-fold or higher somatic mutations than the other 45 cases) and they have germline and somatic mismatch repair gene alterations. The frequently mutated genes are ARID1A (66.7%), PIK3CA (50%), PPP2R1A (18.8%), and KRAS (16.7%). Somatic mutations important for selection of chemotherapeutic agents, such as BRAF, ERBB2, PDGFRB, PGR, and KRAS are found in 27.1% of OCCC cases, indicating clinical importance of exome analysis for OCCC. Our study suggests that the genetic instability caused by either mismatch repair defect or activation of APOBEC play critical roles in OCCC carcinogenesis.

Background With the goal of realizing genome-based personalized healthcare, we have developed a biobank that integrates personal health, genome, and omics data along with biospecimens donated by volunteers of 150,000. Such a large-scale of data integration involves obvious risks of privacy violation. The research use of personal genome and health information is a topic of global discussion with regard to the protection of privacy while promoting scientific advancement. The present paper reports on our plans, current attempts, and accomplishments in addressing security problems involved in data sharing to ensure donor privacy while promoting scientific advancement. Methods Biospecimens and data have been collected in prospective cohort studies with the comprehensive agreement. The sample size of 150,000 participants was required for multiple researches including genome-wide screening of gene by environment interactions, haplotype phasing, and parametric linkage analysis. Results We established the T ohoku M edical M egabank (TMM) data sharing policy: a privacy protection rule that requires physical, personnel, and technological safeguards against privacy violation regarding the use and sharing of data. The proposed policy refers to that of NCBI and that of the Sanger Institute. The proposed policy classifies shared data according to the strength of re-identification risks. Local committees organized by TMM evaluate re-identification risk and assign a security category to a dataset. Every dataset is stored in an assigned segment of a supercomputer in accordance with its security category. A security manager should be designated to handle all security problems at individual data use locations. The proposed policy requires closed networks and IP-VPN remote connections. Conclusion The mission of the biobank is to distribute biological resources most productively. This mission motivated us to collect biospecimens and health data and simultaneously analyze genome/omics data in-house. The biobank also has the mission of improving the quality and quantity of the contents of the biobank. This motivated us to request users to share the results of their research as feedback to the biobank. The TMM data sharing policy has tackled every security problem originating with the missions. We believe our current implementation to be the best way to protect privacy in data sharing. Electronic supplementary material The online version of this article (doi:10.1186/s12911-017-0494-5) contains supplementary material, which is available to authorized users.

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.

Background & Aims Changes in pancreatic calcium levels affect secretion and might be involved in development of chronic pancreatitis (CP). We investigated the association of CP with the transient receptor potential cation channel subfamily V member 6 gene (TRPV6), which encodes a Ca2+-selective ion channel, in an international cohort of patients and in mice. Methods We performed whole-exome DNA sequencing from a patient with idiopathic CP and from his parents, who did not have CP. We validated our findings by sequencing DNA from 300 patients with CP (not associated with alcohol consumption) and 1070 persons from the general population in Japan (control individuals). In replication studies, we sequenced DNA from patients with early-onset CP (20 years or younger) not associated with alcohol consumption from France (n = 470) and Germany (n = 410). We expressed TRPV6 variants in HEK293 cells and measured their activity using Ca2+ imaging assays. CP was induced by repeated injections of cerulein in TRPV6mut/mut mice. Results We identified the variants c.629C>T (p.A210V) and c.970G>A (p.D324N) in TRPV6 in the index patient. Variants that affected function of the TRPV6 product were found in 13 of 300 patients (4.3%) and 1 of 1070 control individuals (0.1%) from Japan (odds ratio [OR], 48.4; 95% confidence interval [CI], 6.3–371.7; P = 2.4 × 10–8). Twelve of 124 patients (9.7%) with early-onset CP had such variants. In the replication set from Europe, 18 patients with CP (2.0%) carried variants that affected the function of the TRPV6 product compared with 0 control individuals (P = 6.2 × 10–8). Variants that did not affect the function of the TRPV6 product (p.I223T and p.D324N) were overrepresented in Japanese patients vs control individuals (OR, 10.9; 95% CI, 4.5–25.9; P = 7.4 × 10–9 for p.I223T and P = .01 for p.D324N), whereas the p.L299Q was overrepresented in European patients vs control individuals (OR, 3.0; 95% CI, 1.9–4.8; P = 1.2 × 10–5). TRPV6mut/mut mice given cerulein developed more severe pancreatitis than control mice, as shown by increased levels of pancreatic enzymes, histologic alterations, and pancreatic fibrosis. Conclusions We found that patients with early-onset CP not associated with alcohol consumption carry variants in TRPV6 that affect the function of its product, perhaps by altering Ca2+ balance in pancreatic cells. TRPV6 regulates Ca2+ homeostasis and pancreatic inflammation.

Article, Scientific reports.9 :102(2019)

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.

We report a family with progressive myoclonic epilepsy who underwent whole-exome sequencing but was negative for pathogenic variants. Similar clinical courses of a devastating neurodegenerative phenotype of two affected siblings were highly suggestive of a genetic etiology, which indicates that the survey of genetic variation by whole-exome sequencing was not comprehensive. To investigate the presence of a variant that remained unrecognized by standard genetic testing, PacBio long-read sequencing was performed. Structural variant (SV) detection using low-coverage (6×) whole-genome sequencing called 17,165 SVs (7,216 deletions and 9,949 insertions). Our SV selection narrowed down potential candidates to only five SVs (two deletions and three insertions) on the genes tagged with autosomal recessive phenotypes. Among them, a 12.4-kb deletion involving the CLN6 gene was the top candidate because its homozygous abnormalities cause neuronal ceroid lipofuscinosis. This deletion included the initiation codon and was found in a GC-rich region containing multiple repetitive elements. These results indicate the presence of a causal variant in a difficult-to-sequence region and suggest that such variants that remain enigmatic after the application of current whole-exome sequencing technology could be uncovered by unbiased application of long-read whole-genome sequencing.

Background and Aims Genome-wide association studies [GWASs] of European populations have identified numerous susceptibility loci for Crohn’s disease [CD]. Susceptibility genes differ by ethnicity, however, so GWASs specific for Asian populations are required. This study aimed to clarify the Japanese-specific genetic background for CD by a GWAS using the Japonica array [JPA] and subsequent imputation with the 1KJPN reference panel. Methods Two independent Japanese case/control sets (Tohoku region [379 CD patients, 1621 controls] and Kyushu region [334 CD patients, 462 controls]) were included. GWASs were performed separately for each population, followed by a meta-analysis. Two additional replication sets [254 + 516 CD patients and 287 + 565 controls] were analysed for top hit single nucleotide polymorphisms [SNPs] from novel genomic regions. Results Genotype data of 4 335 144 SNPs from 713 Japanese CD patients and 2083 controls were analysed. SNPs located in TNFSF15 (rs78898421, Pmeta = 2.59 × 10−26, odds ratio [OR] = 2.10), HLA-DQB1 [rs184950714, pmeta = 3.56 × 10−19, OR = 2.05], ZNF365, and 4p14 loci were significantly associated with CD in Japanese individuals. Replication analyses were performed for four novel candidate loci [p Conclusions RAP1A is a novel susceptibility locus for CD in the Japanese population.

Background/objectives The UBR1 gene encodes the enzyme ubiquitin-protein ligase E3 component n-recognin 1. Loss-of-function mutations in the UBR1 gene cause Johanson-Blizzard syndrome, which involves pancreatic exocrine insufficiency. No previous studies have examined an association of UBR1 variants with pancreatitis, in part due to the large size of the gene. This study aimed to clarify whether UBR1 variants are associated with chronic pancreatitis (CP) by the application of targeted next generation sequencing. Methods Exon sequences of the UBR1 gene from 389 Japanese patients with CP (188 idiopathic, 172 alcoholic, 20 hereditary, 9 familial) were captured by the HaloPlex target enrichment technology and subjected to next generation sequencing. Results Ninety nine point two % of the coding regions of the UBR1 gene could be sequenced by ≥ 20 reads with a mean read depth of 595 and a median depth of 399. Fifteen non-synonymous variants including three novel ones [c.4514T > C (p.I1505T), c.4828C > G (p.H1610D) and c.4856A > T (p.D1619V)] and two synonymous variants were identified in the exonic regions. The frequency of any non-synonymous or synonymous variants was not different between the patients with CP and controls. Conclusions Variants in the UBR1 gene were not associated with CP in Japan.

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.

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.

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.

Population-based prospective cohort studies are indispensable for modern medical research as they provide important knowledge on the influences of many kinds of genetic and environmental factors on the cause of disease. Although traditional cohort studies are mainly conducted using questionnaires and physical examinations, modern cohort studies incorporate omics and genomic approaches to obtain comprehensive physical information, including genetic information. Here, we report the design and midterm results of multi-omics analysis on population-based prospective cohort studies from the Tohoku Medical Megabank (TMM) Project. We have incorporated genomic and metabolomic studies in the TMM cohort study as both metabolome and genome analyses are suitable for high-throughput analysis of large-scale cohort samples. Moreover, an association study between the metabolome and genome show that metabolites are an important intermediate phenotype connecting genetic and lifestyle factors to physical and pathologic phenotypes. We apply our metabolome and genome analyses to large-scale cohort samples in the following studies.

Human leukocyte antigen (HLA) is a gene complex known for its exceptional diversity across populations, importance in organ and blood stem cell transplantation, and associations of specific alleles with various diseases. We constructed a Japanese reference panel of class I HLA genes (ToMMo HLA panel), comprising a distinct set of HLA-A, HLA-B, HLA-C, and HLA-H alleles, by single-molecule, real-time (SMRT) sequencing of 208 individuals included in the 1070 whole-genome Japanese reference panel (1KJPN). For high-quality allele reconstruction, we developed a novel pipeline, Primer-Separation Assembly and Refinement Pipeline (PSARP), in which the SMRT sequencing and additional short-read data were used. The panel consisted of 139 alleles, which were all extended from known IPD-IMGT/HLA sequences, contained 40 with novel variants, and captured more than 96.5% of allelic diversity in 1KJPN. These newly available sequences would be important resources for research and clinical applications including high-resolution HLA typing, genetic association studies, and analyzes of cis-regulatory elements.

Purpose This study evaluated specific relationships between pathogenic mechanisms and genetic polymorphisms in primary open-angle glaucoma (POAG). We analyzed the morphologies of trabeculectomy specimens obtained from patients with familial POAG. Methods We used light microscopy and transmission electron microscopy to examine specimens obtained from 17 eyes of 14 patients with familial POAG. We also conducted exome analyses of two families and used targeted Sanger sequencing to analyze samples obtained from the remaining patients. Results The POAG cases examined in this study were divided into two groups based on morphologic characteristics. Group A eyes (7 eyes from 5 patients) had an abnormally thick trabecular meshwork (TM), whereas group B eyes (10 eyes from 9 patients) had a TM of normal thickness. The characteristics of the outflow routes in group A eyes were remarkable and included apoptotic TM cells, abnormally thickened TM basement membranes, fused TM beams, and occluded Schlemm's canals. All group A patients harbored mutations (F369L, P370L, T377M, and T448P) in the myocilin (MYOC) gene that were not found in group B patients. Conclusions Although age matching of morphologic changes in the outflow routes was impossible due to the small sample size, this study suggests that abnormal TM cells may cause sequential damage in abnormally thickened TM basement membranes, TM cell apoptosis, TM beam fusion, and the occlusion of Schlemm's canals. The four detected MYOC mutations appeared to be associated with morphologic changes in the TM and the underlying pathogenesis of a subtype of familial POAG.

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 −8, rs12946510 was identified as the functional variant that influences gene expression via alteration of Forkhead box protein O1 (FOXO1) binding affinity in vitro. Moreover, expression-quantitative trait locus (e-QTL) analyses showed that the PBC susceptibility allele of rs12946510 was significantly associated with lower endogenous expression of ORMDL3 and GSDMB in whole blood and spleen. This study not only identified the functional variant in chr.17q12-21 and its molecular mechanism through which it conferred susceptibility to PBC, but it also illustrated an efficient systematic approach for post-GWAS analysis that is applicable to other complex diseases.

BACKGROUND. Morphological patterns of Paneth cells are a prognostic biomarker in Western Crohn’s disease (CD) patients, and are associated with autophagy-associated ATG16L1 and NOD2 variants. We hypothesized that genetic determinants of Paneth cell phenotype in other ethnic CD cohorts are distinct but also involved in autophagy. METHODS. We performed a hypothesis-driven analysis of 56 single nucleotide polymorphisms (SNPs) associated with CD susceptibility or known to affect Paneth cell function in 110 Japanese CD patients who underwent ileal resection. We subsequently performed a genome-wide association analysis. Paneth cell phenotype was determined by defensin-5 immunofluorescence. Selected genotype–Paneth cell defect correlations were compared to a Western CD cohort (n = 164). RESULTS. The average percentage of abnormal Paneth cells in Japanese CD was similar to Western CD (P = 0.87), and abnormal Paneth cell phenotype was also associated with early recurrence (P = 0.013). In contrast to Western CD, ATG16L1 T300A was not associated with Paneth cell defect in Japanese CD (P = 0.20). Among the 56 selected SNPs, only LRRK2 M2397T showed significant association with Paneth cell defect (P = 3.62 × 10–4), whereas in the Western CD cohort it was not (P = 0.76). Pathway analysis of LRRK2 and other candidate genes with P less than 5 × 10–4 showed connections with known CD susceptibility genes and links to autophagy and TNF-α networks. CONCLUSIONS. We found dichotomous effects of ATG16L1 and LRRK2 on Paneth cell defect between Japanese and Western CD. Genes affecting Paneth cell phenotype in Japanese CD were also associated with autophagy. Paneth cell phenotype also predicted prognosis in Japanese CD. FUNDING. Helmsley Charitable Trust, Doris Duke Foundation (grant 2014103), Japan Society for the Promotion of Science (KAKENHI grants JP15H04805 and JP15K15284), Crohn’s and Colitis Foundation grant 274415, NIH (grants 1R56DK095820, K01DK109081, and UL1 TR000448).

PURPOSE 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. METHODS 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. RESULTS 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). CONCLUSION 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.

Background In the estimation of repeat numbers in a short tandem repeat (STR) region from high-throughput sequencing data, two types of strategies are mainly taken: a strategy based on counting repeat patterns included in sequence reads spanning the region and a strategy based on estimating the difference between the actual insert size and the insert size inferred from paired-end reads. The quality of sequence alignment is crucial, especially in the former approaches although usual alignment methods have difficulty in STR regions due to insertions and deletions caused by the variations of repeat numbers. Results We proposed a new dynamic programming based realignment method named STR-realigner that considers repeat patterns in STR regions as prior knowledge. By allowing the size change of repeat patterns with low penalty in STR regions, accurate realignment is expected. For the performance evaluation, publicly available STR variant calling tools were applied to three types of aligned reads: synthetically generated sequencing reads aligned with BWA-MEM, those realigned with STR-realigner, those realigned with ReviSTER, and those realigned with GATK IndelRealigner. From the comparison of root mean squared errors between estimated and true STR region size, the results for the dataset realigned with STR-realigner are better than those for other cases. For real data analysis, we used a real sequencing dataset from Illumina HiSeq 2000 for a parent-offspring trio. RepeatSeq and lobSTR were applied to the sequence reads for these individuals aligned with BWA-MEM, those realigned with STR-realigner, ReviSTER, and GATK IndelRealigner. STR-realigner shows the best performance in terms of consistency of the size of estimated STR regions in Mendelian inheritance. Root mean squared error values were also calculated from the comparison of these estimated results with STR region sizes obtained from high coverage PacBio sequencing data, and the results from the realigned sequencing data with STR-realigner showed the least (the best) root mean squared error value. Conclusions The effectiveness of the proposed realignment method for STR regions was verified from the comparison with an existing method on both simulation datasets and real whole genome sequencing dataset.

Epigenome-wide association studies, which searches for blood-based DNA methylation signatures associated with environmental exposures and/or disease susceptibilities, is a promising approach to a better understanding of the molecular aetiology of common diseases. To carry out large-scale epigenome-wide association studies while avoiding false negative detection, an efficient strategy to determine target CpG sites for microarray-based or sequencing-based DNA methylation profiling is essentially needed. Here, we propose and validate a hypothesis that a strategy focusing on CpG sites with high DNA methylation level variability may attain an improved efficacy. Through whole-genome bisulfite sequencing of purified blood cells collected from > 100 apparently healthy subjects, we identified ~2.0 million inter-individually variable CpG sites as potential targets. The efficacy of our strategy was estimated to be 3.7-fold higher than that of the most frequently used strategy. Our catalogue of inter-individually variable CpG sites will accelerate the discovery of clinically relevant DNA methylation biomarkers in future epigenome-wide association studies., 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.

Background Genome-wide association studies have revealed associations between single-nucleotide polymorphisms (SNPs) and phenotypes such as disease symptoms and drug tolerance. To address the small sample size for rare variants, association studies tend to group gene or pathway level variants and evaluate the effect on the set of variants. One of such strategies, known as the sequential kernel association test (SKAT), is a widely used collapsing method. However, the reported p-values from SKAT tend to be biased because the asymptotic property of the statistic is used to calculate the p-value. Although this bias can be corrected by applying permutation procedures for the test statistics, the computational cost of obtaining p-values with high resolution is prohibitive. Results To address this problem, we devise an adaptive SKAT procedure termed AP-SKAT that efficiently classifies significant SNP sets and ranks them according to the permuted p-values. Our procedure adaptively stops the permutation test when the significance level is outside some confidence interval of the estimated p-value for a binomial distribution. To evaluate the performance, we first compare the power and sample size calculation and the type I error rates estimate of SKAT, SKAT-O, and the proposed procedure using genotype data in the SKAT R package and from 1000 Genome Project. Through computational experiments using whole genome sequencing and SNP array data, we show that our proposed procedure is highly efficient and has comparable accuracy to the standard procedure. Conclusions For several types of genetic data, the developed procedure could achieve competitive power and sample size under small and large sample size conditions with controlling considerable type I error rates, and estimate p-values of significant SNP sets that are consistent with those estimated by the standard permutation test within a realistic time. This demonstrates that the procedure is sufficiently powerful for recent whole genome sequencing and SNP array data with increasing numbers of phenotypes. Additionally, this procedure can be used in other association tests by employing alternative methods to calculate the statistics.

Background Two types of approaches are mainly considered for the repeat number estimation in short tandem repeat (STR) regions from high-throughput sequencing data: approaches directly counting repeat patterns included in sequence reads spanning the region and approaches based on detecting the difference between the insert size inferred from aligned paired-end reads and the actual insert size. Although the accuracy of repeat numbers estimated with the former approaches is high, the size of target STR regions is limited to the length of sequence reads. On the other hand, the latter approaches can handle STR regions longer than the length of sequence reads. However, repeat numbers estimated with the latter approaches is less accurate than those with the former approaches. Results We proposed a new statistical model named coalescentSTR that estimates repeat numbers from paired-end read distances for multiple individuals simultaneously by connecting the read generative model for each individual with their genealogy. In the model, the genealogy is represented by handling coalescent trees as hidden variables, and the summation of the hidden variables is taken on coalescent trees sampled based on phased genotypes located around a target STR region with Markov chain Monte Carlo. In the sampled coalescent trees, repeat number information from insert size data is propagated, and more accurate estimation of repeat numbers is expected for STR regions longer than the length of sequence reads. For finding the repeat numbers maximizing the likelihood of the model on the estimation of repeat numbers, we proposed a state-of-the-art belief propagation algorithm on sampled coalescent trees. Conclusions We verified the effectiveness of the proposed approach from the comparison with existing methods by using simulation datasets and real whole genome and whole exome data for HapMap individuals analyzed in the 1000 Genomes Project.

Relationship between structural variants of enzymes and metabolic phenotypes in human population was investigated based on the association study of metabolite quantitative traits with whole genome sequence data for 512 individuals from a population cohort. We identified five significant associations between metabolites and non-synonymous variants. Four of these non-synonymous variants are located in enzymes involved in metabolic disorders, and structural analyses of these moderate non-synonymous variants demonstrate that they are located in peripheral regions of the catalytic sites or related regulatory domains. In contrast, two individuals with larger changes of metabolite levels were also identified, and these individuals retained rare variants, which caused non-synonymous variants located near the catalytic site. These results are the first demonstrations that variant frequency, structural location, and effect for phenotype correlate with each other in human population, and imply that metabolic individuality and susceptibility for diseases may be elicited from the moderate variants and much more deleterious but rare variants.

Genetic polymorphism within the 9q32 locus is linked with increased risk of several diseases, including Crohn’s disease (CD), primary biliary cholangitis (PBC) and leprosy. The most likely disease-causing gene within 9q32 is TNFSF15, which encodes the pro-inflammatory cytokine TNF super-family member 15, but it was unknown whether these disparate diseases were associated with the same genetic variance in 9q32, and how variance within this locus might contribute to pathology. Using genetic data from published studies on CD, PBC and leprosy we revealed that bearing a T allele at rs6478108/rs6478109 (r2 = 1) or rs4979462 was significantly associated with increased risk of CD and decreased risk of leprosy, while the T allele at rs4979462 was associated with significantly increased risk of PBC. In vitro analyses showed that the rs6478109 genotype significantly affected TNFSF15 expression in cells from whole blood of controls, while functional annotation using publicly-available data revealed the broad cell type/tissue-specific regulatory potential of variance at rs6478109 or rs4979462. In summary, we provide evidence that variance within TNFSF15 has the potential to affect cytokine expression across a range of tissues and thereby contribute to protection from infectious diseases such as leprosy, while increasing the risk of immune-mediated diseases including CD and PBC.

Background & Aims There is still a risk for hepatocellular carcinoma (HCC) development after eradication of hepatitis C virus (HCV) infection with antiviral agents. We investigated genetic factors associated with the development of HCC in patients with a sustained virologic response (SVR) to treatment for chronic HCV infection. Methods We obtained genomic DNA from 457 patients in Japan with a SVR to interferon-based treatment for chronic HCV infection from 2007 through 2015. We conducted a genome-wide association study (GWAS), followed by a replication analysis of 79 candidate single nucleotide polymorphisms (SNPs) in an independent set of 486 patients in Japan. The study end point was HCC diagnosis or confirmation of lack of HCC (at follow-up examinations until December 2014 in the GWAS cohort, and until January 2016 in the replication cohort). We collected clinical and laboratory data from all patients. We analyzed expression levels of candidate gene variants in human hepatic stellate cells, rats with steatohepatitis caused by a choline-deficient L-amino acid-defined diet, and a mouse model of liver injury caused by administration of carbon tetrachloride. We also analyzed expression levels in liver tissues of patients with chronic HCV infection with different stages of fibrosis or tumors vs patients without HCV infection (controls). Results We found a strong association between the SNP rs17047200, located within the intron of the tolloid like 1 gene ( TLL1 ) on chromosome 4, and development of HCC; there was a genome-wide level of significance when the results of the GWAS and replication study were combined (odds ratio, 2.37; P = 2.66 × 10 −8 ). Multivariate analysis showed rs17047200 AT/TT to be an independent risk factor for HCC (hazard ratio, 1.78; P = .008), along with male sex, older age, lower level of albumin, advanced stage of hepatic fibrosis, presence of diabetes, and higher post-treatment level of α-fetoprotein. Combining the rs17047200 genotype with other factors, we developed prediction models for HCC development in patients with mild or advanced hepatic fibrosis. Levels of TLL1 messenger RNA (mRNA) in human hepatic stellate cells increased with activation. Levels of Tll1 mRNA increased in liver tissues of rodents with hepatic fibrogenesis compared with controls. Levels of TLL1 mRNA increased in liver tissues of patients with progression of fibrosis. Gene expression levels of TLL1 short variants, including isoform 2, were higher in patients with rs17047200 AT/TT. Conclusions In a GWAS, we identified the association between the SNP rs17047200, within the intron of TLL1 , and development of HCC in patients who achieved an SVR to treatment for chronic HCV infection. We found levels of Tll1/TLL1 mRNA to be increased in rodent models of liver injury and liver tissues of patients with fibrosis, compared with controls. We propose that this SNP might affect splicing of TLL1 mRNA, yielding short variants with high catalytic activity that accelerates hepatic fibrogenesis and carcinogenesis. Further studies are needed to determine how rs17047200 affects TLL1 mRNA levels, splicing, and translation, as well as the prevalence of this variant among other patients with HCC. Tests for the TLL1 SNP might be used to identify patients at risk for HCC after an SVR to treatment of HCV infection.

To analyze cell lineage-specific transitions in global transcriptional and epigenetic changes during retinogenesis, we purified retinal cells from normal mice during postnatal development into two fractions, namely, photoreceptors and other retinal cells, based on Cd73 expression, and performed RNA sequencing and ChIP sequencing of H3K27me3 and H3K4me3. Genes expressed in the photoreceptor lineage were marked with H3K4me3 in the Cd73-positive cell fraction; however, the level of H3K27me3 was very low in both Cd73-positive and -negative populations. H3K27me3 may be involved in spatio-temporal onset of a subset of bipolar-related genes. Subsets of genes expressed in amacrine and retinal ganglion cells, which are early-born retinal cell types, were suggested to be maintained in a silent state by H3K27me3 during late-stage retinogenesis. In the outer nuclear layer, upregulation of Rho and rod-related genes were observed in Ezh2-ablated retina, suggesting a role for H3K27me3 in the maintenance of proper expression levels. Taken together, our data on the transition of lineage-specific molecular signatures during development suggest that histone methylation is involved in retinal differentiation and maintenance through cell lineage-specific mechanisms.

A genome-wide association study (GWAS) for cold medicine-related Stevens–Johnson syndrome (CM-SJS) with severe ocular complications (SOC) was performed in a Japanese population. A recently developed ethnicity-specific array with genome-wide imputation that was based on the whole-genome sequences of 1070 unrelated Japanese individuals was used. Validation analysis with additional samples from Japanese individuals and replication analysis using samples from Korean individuals identified two new susceptibility loci on chromosomes 15 and 16. This study might suggest the usefulness of GWAS using the ethnicity-specific array and genome-wide imputation based on large-scale whole-genome sequences. Our findings contribute to the understanding of genetic predisposition to CM-SJS with SOC.