Search

While the advent of GWAS more than a decade ago has ushered in remarkable advances in our understanding of complex traits, the limitations of single-SNP analysis have also led to the development of several other approaches. Simulation studies have shown that the regional heritability mapping (RHM) method, which makes use of multiple adjacent SNPs jointly to estimate the genetic effect of a given region of the genome, generally has higher detection power than single-SNP GWAS. However, thus far its use has been mostly limited to agricultural settings, and its potential for the discovery of new genes in human diseases is yet to be fully exploited. In this study, by applying the RHM method to primary biliary cholangitis (PBC) in the Japanese population, we identified three novel loci (STAT4, ULK4, and KCNH5) at the genome-wide significance level, two of which (ULK4 and KCNH5) have not been found associated with PBC in any population previously. Notably, these genes could not be detected by using conventional single-SNP GWAS, highlighting the potential of the RHM method for the detection of new susceptibility loci in human diseases. These findings thereby provide strong empirical evidence that RHM is an effective and practical complementary approach to GWAS in this context. Also, liver tissue mRNA microarray analysis revealed higher gene expression levels in ULK4 in PBC patients (P, 原発性胆汁性胆管炎の新たな遺伝要因を同定 --ヒト全ゲノム領域へのRHM法による世界初の成果--. 京都大学プレスリリース. 2021-04-09.

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.

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.

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

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.

Article, Scientific reports.9 :102(2019)

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.

Neural stem cells (NSCs) are considered to be the cell of origin of glioblastoma multiforme (GBM). However, the genetic alterations that transform NSCs into glioma-initiating cells remain elusive. Using a unique transposon mutagenesis strategy that mutagenizes NSCs in culture, followed by additional rounds of mutagenesis to generate tumors in vivo, we have identified genes and signaling pathways that can transform NSCs into glioma-initiating cells. Mobilization of Sleeping Beauty transposons in NSCs induced the immortalization of astroglial-like cells, which were then able to generate tumors with characteristics of the mesenchymal subtype of GBM on transplantation, consistent with a potential astroglial origin for mesenchymal GBM. Sequence analysis of transposon insertion sites from tumors and immortalized cells identified more than 200 frequently mutated genes, including human GBM-associated genes, such as Met and Nf1 , and made it possible to discriminate between genes that function during astroglial immortalization vs. later stages of tumor development. We also functionally validated five GBM candidate genes using a previously undescribed high-throughput method. Finally, we show that even clonally related tumors derived from the same immortalized line have acquired distinct combinations of genetic alterations during tumor development, suggesting that tumor formation in this model system involves competition among genetically variant cells, which is similar to the Darwinian evolutionary processes now thought to generate many human cancers. This mutagenesis strategy is faster and simpler than conventional transposon screens and can potentially be applied to any tissue stem/progenitor cells that can be grown and differentiated in vitro.

Genetic materials are one of the most important and fundamental research resources for studying biological phenomena. Scientific need for genetic materials has been increasing and will never cease. Ever since it was established as RIKEN DNA Bank in 1987, the Gene Engineering Division of RIKEN BioResource Center (BRC) has been engaged in the collection, maintenance, storage, propagation, quality control, and distribution of genetic resources developed mainly by the Japanese research community. When RIKEN BRC was inaugurated in 2001, RIKEN DNA Bank was incorporated as one of its six Divisions, the Gene Engineering Division. The Gene Engineering Division was selected as a core facility for the genetic resources of mammalian and microbe origin by the National BioResource Project (NBRP) of the Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan in 2002. With support from the scientific community, the Division now holds over 3 million clones of genetic materials for distribution. The genetic resources include cloned DNAs, gene libraries (e.g., cDNA and genomic DNA cloned into phage, cosmid, BAC, phosmid, and YAC), vectors, hosts, recombinant viruses, and ordered library sets derived from animal cells, including human and mouse cells, microorganisms, and viruses. Recently genetic materials produced by a few MEXT national research projects were transferred to the Gene Engineering Division for further dissemination. The Gene Engineering Division performs rigorous quality control of reproducibility, restriction enzyme mapping and nucleotide sequences of clones to ensure the reproducibility of in vivo and in vitro experiments. Users can easily access our genetic materials through the internet and obtain the DNA resources for a minimal fee. Not only the materials, but also information of features and technology related to the materials are provided via the web site of RIKEN BRC. Training courses are also given to transfer the technology for handling viral vectors. RIKEN BRC supports scientists around the world in the use of valuable genetic materials.

Alternative splicing plays a prominent role in eukaryotic gene regulations that allow a single gene to generate the multiple mRNA products. The recent advent of GeneChipreg Human Exon 1.0 ST Array enables us to measure the exon expression profiles of human cells on a genome-wide scale. With this advent, analysis of functional gene regulation could be extended to detect not only differentially expressed genes, but also specific splicing events that occur in target cells, but not in normal controls. We address some statistical issues for the identification of biomarker splice variations with exon expression data. The proposed method involves the following steps: (1) Whole transcript analysis with the nonparametric analysis of variance (ANOVA) to identify potential biomarkers that present specific splice variations. (2) Meta-analysis for discriminating non-specific splice variations that are caused by clinical heterogeneity in the collected samples. In the analysis of human cells, controlling non-specific splicing factors is essential for success in the detection of biomarker splice variations because splice patterns are possibly affected by inter-individual differences in the collected samples. We demonstrate its utility and perform a whole transcript analysis of exon expression profiles of colorectal carcinoma.