The Dynamics of Transcriptional Activation by Hepatic Reprogramming Factors.

Specific combinations of two transcription factors (Hnf4α plus Foxa1, Foxa2, or Foxa3) can induce direct conversion of mouse fibroblasts into hepatocyte-like cells. However, the molecular mechanisms underlying hepatic reprogramming are largely unknown. Here, we show that the Foxa protein family members and Hnf4α sequentially and cooperatively bind to chromatin to activate liver-specific gene expression. Although all Foxa proteins bind to and open regions of closed chromatin as pioneer factors, Foxa3 has the unique potential of transferring from the distal to proximal regions of the transcription start site of target genes, binding RNA polymerase II, and co-traversing target genes. These distinctive characteristics of Foxa3 are essential for inducing the hepatic fate in fibroblasts. Similar functional coupling of transcription factors to RNA polymerase II may occur in other contexts whereby transcriptional activation can induce cell differentiation.

rs10924104 in the expression enhancer motif of CD58 confers susceptibility to human autoimmune diseases.

CD58 plays roles in cell adhesion and co-stimulation with antigen presentation from major histocompatibility complex class II on antigen-presenting cells to T-cell antigen receptors on naïve T cells. CD58 reportedly contributes to the development of various human autoimmune diseases. Recently, genome-wide association studies (GWASs) identified CD58 as a susceptibility locus for autoimmune diseases such as systemic lupus erythematosus (SLE), multiple sclerosis (MS), and primary biliary cholangitis (PBC). However, the primary functional variant and molecular mechanisms of susceptibility to autoimmune diseases in the CD58 locus were not clarified. Here, rs10924104, located in the ZNF35-binding motif within the gene expression regulatory motif, was identified as the primary functional variant for SLE, MS, and PBC among genetic variants showing stronger linkage disequilibrium (LD) with GWAS-lead variants in the CD58 locus. Expression-quantitative trait locus (e-QTL) data for each distinct blood cell type and in vitro functional analysis using the CRISPR/Cas9 system corroborated the functional role of rs10924104 in the upregulation of CD58 transcription by the disease-risk allele. Additionally, the strength of disease susceptibility observed in the CD58 locus could be accounted for by the strength of LD between rs10924104 and each GWAS-lead variant. In conclusion, the present study demonstrated for the first time the existence of a shared autoimmune disease-related primary functional variant (i.e., rs10924104) that regulates the expression of CD58. Clarifying the molecular mechanism of disease susceptibility derived from such a shared genetic background is important for understanding human autoimmune diseases and human immunology.

HLA-DQA1*05 and upstream variants of PPARGC1B are associated with infliximab persistence in Japanese Crohn's disease patients.

Recently, the HLA-DQA1*05 (rs2097432) genetic variation has been reported to be linked to early infliximab (IFX) treatment failure in the Caucasian Crohn's disease (CD) population, but that evidence is scarce in the Asian population. This study aimed to investigate the relationship between rs2097432 and the cumulative discontinuation-free time of IFX (IFX persistence) in 189 Japanese biologics-naive CD patients. We also performed a genome-wide association study (GWAS) to discover novel genetic predictors for IFX persistence. The C allele of rs2097432 significantly increased the risk of early discontinuation of IFX [Hazard ratio (HR) = 2.23 and P-value = 0.026]. In GWAS, one locus tagged by rs73277969, located upstream of PPARGC1B which attenuates macrophage-mediated inflammation, reached genome-wide significance (HR = 6.04 and P-value = 7.93E-9). Pathway analysis suggested association of signaling by PDGF and FCGR activation signaling with IFX persistence (P-value = 8.56E-5 and 5.80E-4, respectively).

CFH-CFHR1 hybrid genes in two cases of atypical hemolytic uremic syndrome.

Atypical hemolytic uremic syndrome (aHUS) is a rare complement-mediated disease that manifests as the triad of thrombotic microangiopathy. We identified two aHUS patients with neither anti-complement factor H (CFH) antibodies nor causative variants of seven aHUS-related genes (CFH, CFI, CFB, C3, MCP, THBD, and DGKE); however, their plasma showed increased levels of hemolysis by hemolytic assay, which strongly suggests CFH-related abnormalities. Using a copy number variation (CNV) analysis of the CFH/CFHR gene cluster, we identified CFH-CFHR1 hybrid genes in these patients. We verified the absence of aHUS-related abnormal CNVs of the CFH gene in control genomes of 2036 individuals in the general population, which suggests that pathogenicity is related to these hybrid genes. Our study emphasizes that, for patients suspected of having aHUS, it is important to perform an integrated analysis based on a clinical examination, functional analysis, and detailed genetic investigation.

rs2013278 in the multiple immunological-trait susceptibility locus CD28 regulates the production of non-functional splicing isoforms.

BACKGROUND: Ligation of CD28 with ligands such as CD80 or CD86 provides a critical second signal alongside antigen presentation by class II major histocompatibility complex expressed on antigen-presenting cells through the T cell antigen receptor for naïve T cell activation. A number of studies suggested that CD28 plays an important role in the pathogenesis of various human diseases. Recent genome-wide association studies (GWASs) identified CD28 as a susceptibility locus for lymphocyte and eosinophil counts, multiple sclerosis, ulcerative colitis, celiac disease, rheumatoid arthritis, asthma, and primary biliary cholangitis. However, the primary functional variant and molecular mechanisms of disease susceptibility in this locus remain to be elucidated. This study aimed to identify the primary functional variant from thousands of genetic variants in the CD28 locus and elucidate its functional effect on the CD28 molecule. RESULTS: Among the genetic variants exhibiting stronger linkage disequilibrium (LD) with all GWAS-lead variants in the CD28 locus, rs2013278, located in the Rbfox binding motif related to splicing regulation, was identified as a primary functional variant related to multiple immunological traits. Relative endogenous expression levels of CD28 splicing isoforms (CD28i and CD28Δex2) compared with full-length CD28 in allele knock-in cell lines generated using CRISPR/Cas9 were directly regulated by rs2013278 (P < 0.05). Although full-length CD28 protein expressed on Jurkat T cells showed higher binding affinity for CD80/CD86, both CD28i and CD28Δex2 encoded loss-of-function isoforms. CONCLUSION: The present study demonstrated for the first time that CD28 has a shared disease-related primary functional variant (i.e., rs2013278) that regulates the CD28 alternative splicing that generates loss-of-function isoforms. They reduce disease risk by inducing anergy of effector T cells that over-react to autoantigens and allergens.

Genome-Wide Association Study and Transcriptome of Japanese Patients with Developmental Dysplasia of the Hip Demonstrates an Association with the Ferroptosis Signaling Pathway.

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

Prediction Model with HLA-A*33:03 Reveals Number of Days to Develop Liver Cancer from Blood Test.

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

X Chromosome Contribution to the Genetic Architecture of Primary Biliary Cholangitis.

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.

rs9459874 and rs1012656 in CCR6/FGFR1OP confer susceptibility to primary biliary cholangitis.

Primary biliary cholangitis (PBC) is a chronic cholestatic autoimmune liver disease that appears to be strongly influenced by genetic factors. Recently, an international meta-analysis of genome-wide association studies (GWAS) identified CC-Motif Chemokine Receptor-6 (CCR6) and FGFR1 Oncogene-Partner (FGFR1OP) as PBC-susceptibility genes. However, the lead single nucleotide polymorphisms (SNPs) of CCR6/FGFR1OP showed low linkage disequilibrium with each other in East Asian and European populations. Additionally, the primary functional variants and the molecular mechanisms responsible for PBC-susceptibility remain unclear. Here, among the PBC-susceptibility SNPs identified by high-density association mapping in our previous meta-GWAS (Patients: n = 10,516; healthy controls: n = 20,772) within the CCR6/FGFR1OP locus, rs9459874 and rs1012656 were identified as primary functional variants. These functional variants accounted for the effects of GWAS-identified lead SNPs in CCR6/FGFR1OP. Additionally, the roles of rs9459874 and rs1012656 in regulating FGFR1OP transcription and CCR6 translation, respectively, were supported by expression quantitative trait loci (eQTL) analysis and gene editing technology using the CRISPR/Cas9 system. Immunohistochemistry showed higher expression of CCR6 protein in the livers of patients with PBC than in those of a non-diseased control. In conclusion, we identified primary functional variants in CCR6/FGFR1OP and revealed the molecular mechanisms by which these variants confer PBC-susceptibility in an eQTL-dependent or -independent manner. The approach in this study is applicable for the elucidation of the pathogenesis of other autoimmune disorders in which CCR6/FGFR1OP is known as a susceptibility locus, as well as PBC.

Genome-wide Survival Analysis for Macular Neovascularization Development in Central Serous Chorioretinopathy Revealed Shared Genetic Susceptibility with Polypoidal Choroidal Vasculopathy.

PURPOSE: To identify susceptibility genes for macular neovascularization (MNV) development in central serous chorioretinopathy (CSC). DESIGN: Genome-wide survival analysis using a longitudinal cohort study. PARTICIPANTS: We included 402 and 137 patients with CSC but without MNV at their first visit from the Kyoto CSC Cohort and Kobe CSC dataset, respectively. All patients underwent detailed ophthalmic examinations, including multimodal imaging, such as fundus autofluorescence, spectral-domain OCT, and fluorescein angiography/indocyanine green angiography or OCT angiography. METHODS: We conducted a genome-wide survival analysis using the Kyoto CSC Cohort. We applied the Cox proportional hazard model to adjust for age, sex, and the first principal component. Single nucleotide polymorphisms (SNPs) with P values < 1.0 × 10-5 were carried forward to the replication in the Kobe CSC dataset. Moreover, we evaluated the contribution of previously reported age-related macular degeneration (AMD) susceptibility loci. We used FUMA and ToppFun for the functional enrichment analysis. MAIN OUTCOME MEASURES: The association between SNPs and MNV development in patients with CSC. RESULTS: Rs370974631 near ARMS2 displayed a genome-wide significant association in the meta-analysis of discovery and replication result (hazard ratio [HR]meta, 3.63; Pmeta = 5.76 × 10-9). Among previously reported AMD susceptibility loci, we additionally identified CFH rs800292 (HR, 0.39, P = 2.55 × 10-4), COL4A3 rs4276018 (HR, 0.26, P = 1.56 × 10-3), and B3GALTL rs9564692 (HR, 0.56, P = 8.30 × 10-3) as susceptibility loci for MNV development in CSC. The functional enrichment analysis revealed significant enrichment of 8 pathways (GO:0051561, GO:0036444, GO:0008282, GO:1990246, GO:0015272, GO:0030955, GO:0031420, and GO:0005242) related to ion transport. CONCLUSIONS: ARMS2, CFH, COL4A3, and B3GALTL were identified as susceptibility genes for MNV development in CSC. These 4 genes are known as susceptibility genes for AMD, whereas COL4A3 and B3GALTL were previously reported to be polypoidal choroidal vasculopathy (PCV)-specific susceptibility genes. Our findings revealed the shared genetic susceptibility between PCV and MNV secondary to CSC.

An international genome-wide meta-analysis of primary biliary cholangitis: Novel risk loci and candidate drugs.

BACKGROUNDS & AIMS: Primary biliary cholangitis (PBC) is a chronic liver disease in which autoimmune destruction of the small intrahepatic 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 5 European and 2 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, which each play key roles in immunity. Pathway analysis reiterated the likely importance of pattern recognition receptor and TNF signalling, JAK-STAT signalling, and differentiation of T helper (TH)1 and TH17 cells in the pathogenesis of this disease. Drug efficacy screening identified several medications predicted to be therapeutic in PBC, some of which are 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, the UK, or the 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.

Variants That Affect Function of Calcium Channel TRPV6 Are Associated With Early-Onset Chronic Pancreatitis.

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.

Practical guide for managing large-scale human genome data in research.

Studies in human genetics deal with a plethora of human genome sequencing data that are generated from specimens as well as available on public domains. With the development of various bioinformatics applications, maintaining the productivity of research, managing human genome data, and analyzing downstream data is essential. This review aims to guide struggling researchers to process and analyze these large-scale genomic data to extract relevant information for improved downstream analyses. Here, we discuss worldwide human genome projects that could be integrated into any data for improved analysis. Obtaining human whole-genome sequencing data from both data stores and processes is costly; therefore, we focus on the development of data format and software that manipulate whole-genome sequencing. Once the sequencing is complete and its format and data processing tools are selected, a computational platform is required. For the platform, we describe a multi-cloud strategy that balances between cost, performance, and customizability. A good quality published research relies on data reproducibility to ensure quality results, reusability for applications to other datasets, as well as scalability for the future increase of datasets. To solve these, we describe several key technologies developed in computer science, including workflow engine. We also discuss the ethical guidelines inevitable for human genomic data analysis that differ from model organisms. Finally, the future ideal perspective of data processing and analysis is summarized.

Long non-coding RNA MIR4300HG polymorphisms are associated with postoperative nausea and vomiting: a genome-wide association study.

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.

Study Profile of the Tohoku Medical Megabank Community-Based Cohort Study.

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

Common risk variants in NPHS1 and TNFSF15 are associated with childhood steroid-sensitive nephrotic syndrome.

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 (Pmeta=6.71E-28, OR=1.88) and TNFSF15 (Pmeta=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).

Estimating carrier frequencies of newborn screening disorders using a whole-genome reference panel of 3552 Japanese individuals.

Incidence rates of Mendelian diseases vary among ethnic groups, and frequencies of variant types of causative genes also vary among human populations. In this study, we examined to what extent we can predict population frequencies of recessive disorders from genomic data, and explored better strategies for variant interpretation and classification. We used a whole-genome reference panel from 3552 general Japanese individuals constructed by the Tohoku Medical Megabank Organization (ToMMo). Focusing on 32 genes for 17 congenital metabolic disorders included in newborn screening (NBS) in Japan, we identified reported and predicted pathogenic variants through variant annotation, interpretation, and multiple ways of classifications. The estimated carrier frequencies were compared with those from the Japanese NBS data based on 1,949,987 newborns from a previous study. The estimated carrier frequency based on genomic data with a recent guideline of variant interpretation for the PAH gene, in which defects cause hyperphenylalaninemia (HPA) and phenylketonuria (PKU), provided a closer estimate to that by the observed incidence than the other methods. In contrast, the estimated carrier frequencies for SLC25A13, which causes citrin deficiency, were much higher compared with the incidence rate. The results varied greatly among the 11 NBS diseases with single responsible genes; the possible reasons for departures from the carrier frequencies by reported incidence rates were discussed. Of note, (1) the number of pathogenic variants increases by including additional lines of evidence, (2) common variants with mild effects also contribute to the actual frequency of patients, and (3) penetrance of each variant remains unclear.

Omics research project on prospective cohort studies from the Tohoku Medical Megabank Project.

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.

Construction of full-length Japanese reference panel of class I HLA genes with single-molecule, real-time sequencing.

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.

A 12-kb structural variation in progressive myoclonic epilepsy was newly identified by long-read whole-genome sequencing.

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.