jMorp: Japanese Multi-Omics Reference Panel update report 2023

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

Functional Characterization of 29 Cytochrome P450 4F2 Variants Identified in a Population of 8380 Japanese Subjects and Assessment of Arachidonic Acid ω-Hydroxylation

Cytochrome P450 4F2 (CYP4F2) is an enzyme that is involved in the metabolism of arachidonic acid (AA), vitamin E and K, and xenobiotics including drugs. CYP4F2*3 polymorphism (rs2108622; c.1297G>A; p.Val433Met) has been associated with hypertension, ischemic stroke, and variation in the effectiveness of the anticoagulant drug warfarin. In this study, we characterized wild-type CYP4F2 and 28 CYP4F2 variants, including a Val433Met substitution, detected in 8380 Japanese subjects. The CYP4F2 variants were heterologously expressed in 293FT cells to measure the concentrations of CYP4F2 variant holoenzymes using carbon monoxide-reduced difference spectroscopy, where the wild type and 18 holoenzyme variants showed a peak at 450 nm. Kinetic parameters [Vmax , substrate concentration producing half of Vmax (S50 ), and intrinsic clearance (CL int ) as Vmax /S50 ] of AA ω-hydroxylation were determined for the wild type and 21 variants with enzyme activity. Compared with the wild type, two variants showed significantly decreased CL int values for AA ω-hydroxylation. The values for seven variants could not be determined because no enzymatic activity was detected at the highest substrate concentration used. Three-dimensional structural modeling was performed to determine the reason for reduced enzymatic activity of the CYP4F2 variants. Our findings contribute to a better understanding of CYP4F2 variant-associated diseases and possible future therapeutic strategies. SIGNIFICANCE STATEMENT: CYP4F2 is involved in the metabolism of arachidonic acid and vitamin K, and CYP4F2*3 polymorphisms have been associated with hypertension and variation in the effectiveness of the anticoagulant drug warfarin. This study presents a functional analysis of 28 CYP4F2 variants identified in Japanese subjects, demonstrating that seven gene polymorphisms cause loss of CYP4F2 function, and proposes structural changes that lead to altered function.

Functional Characterization of 12 Dihydropyrimidinase Allelic Variants in Japanese Individuals for the Prediction of 5-Fluorouracil Treatment-Related Toxicity

The drug 5-fluorouracil (5-FU) is the first-choice chemotherapeutic agent against advanced-stage cancers. However, 10% to 30% of treated patients experience grade 3 to 4 toxicity. The deficiency of dihydropyrimidinase (DHPase), which catalyzes the second step of the 5-FU degradation pathway, is correlated with the risk of developing toxicity. Thus, genetic polymorphisms within DPYS, the DHPase-encoding gene, could potentially serve as predictors of severe 5-FU-related toxicity. We identified 12 novel DPYS variants in 3554 Japanese individuals, but the effects of these mutations on function remain unknown. In the current study, we performed in vitro enzymatic analyses of the 12 newly identified DHPase variants. Dihydrouracil or dihydro-5-FU hydrolytic ring-opening kinetic parameters, K_m and V_max , and intrinsic clearance (CL_int = V_max /K_m ) of the wild-type DHPase and eight variants were measured. Five of these variants (R118Q, H295R, T418I, Y448H, and T513A) showed significantly reduced CL_int compared with that in the wild-type. The parameters for the remaining four variants (V59F, D81H, T136M, and R490H) could not be determined as dihydrouracil and dihydro-5-FU hydrolytic ring-opening activity was undetectable. We also determined DHPase variant protein stability using cycloheximide and bortezomib. The mechanism underlying the observed changes in the kinetic parameters was clarified using blue-native polyacrylamide gel electrophoresis and three-dimensional structural modeling. The results suggested that the decrease or loss of DHPase enzymatic activity was due to reduced stability and oligomerization of DHPase variant proteins. Our findings support the use of DPYS polymorphisms as novel pharmacogenomic markers for predicting severe 5-FU-related toxicity in the Japanese population. SIGNIFICANCE STATEMENT: DHPase contributes to the degradation of 5-fluorouracil, and genetic polymorphisms that cause decreased activity of DHPase can cause severe toxicity. In this study, we performed functional analysis of 12 DHPase variants in the Japanese population and identified 9 genetic polymorphisms that cause reduced DHPase function. In addition, we found that the ability to oligomerize and the conformation of the active site are important for the enzymatic activity of DHPase.

Construction of a trio-based structural variation panel utilizing activated T lymphocytes and long-read sequencing technology

Long-read sequencing technology enable better characterization of structural variants (SVs). To adapt the technology to population-scale analyses, one critical issue is to obtain sufficient amount of high-molecular-weight genomic DNA. Here, we propose utilizing activated T lymphocytes, which can be established efficiently in a biobank to stably supply high-grade genomic DNA sufficiently. We conducted nanopore sequencing of 333 individuals constituting 111 trios with high-coverage long-read sequencing data (depth 22.2x, N50 of 25.8 kb) and identified 74,201 SVs. Our trio-based analysis revealed that more than 95% of the SVs were concordant with Mendelian inheritance. We also identified SVs associated with clinical phenotypes, all of which appear to be stably transmitted from parents to offspring. Our data provide a catalog of SVs in the general Japanese population, and the applied approach using the activated T-lymphocyte resource will contribute to biobank-based human genetic studies focusing on SVs at the population scale.

Long-read sequencing on activated T-cells from a sample of 333 Japanese individuals (representing 111 parent-offspring trios) provides a useful reference dataset of structural variation in the Japanese population.

Importance of Rare DPYD Genetic Polymorphisms for 5-Fluorouracil Therapy in the Japanese Population

Dihydropyrimidine dehydrogenase (DPD), encoded by the DPYD gene, is the rate-limiting enzyme in 5-fluorouracil (5-FU) degradation. In Caucasians, four DPYD risk variants are recognized to be responsible for interindividual variations in the development of 5-FU toxicity. However, these risk variants have not been identified in Asian populations. Recently, 41 DPYD allelic variants, including 15 novel single nucleotide variants, were identified in 3,554 Japanese individuals by analyzing their whole-genome sequences; however, the effects of these variants on DPD enzymatic activity remain unknown. In the present study, an in vitro analysis was performed on 41 DPD allelic variants and three DPD risk variants to elucidate the changes in enzymatic activity. Wild-type and 44 DPD-variant proteins were heterologously expressed in 293FT cells. DPD expression levels and dimerization of DPD were determined by immunoblotting after SDS-PAGE and blue native PAGE, respectively. The enzymatic activity of DPD was evaluated by quantification of dihydro-5-FU, a metabolite of 5-FU, using high-performance liquid chromatography-tandem mass spectrometry. Moreover, we used 3D simulation modeling to analyze the effect of amino acid substitutions on the conformation of DPD. Among the 41 DPD variants, seven exhibited drastically decreased intrinsic clearance (CL_int) compared to the wild-type protein. Moreover, R353C and G926V exhibited no enzymatic activity, and the band patterns observed in the immunoblots after blue native PAGE indicated that DPD dimerization is required for its enzymatic activity. Our data suggest that these variants may contribute to the significant inter-individual variability observed in the pharmacokinetics and pharmacodynamics of 5-FU. In our study, nine DPD variants exhibited drastically decreased or no enzymatic activity due to dimerization inhibition or conformational changes in each domain. Especially, the rare DPYD variants, although at very low frequencies, may serve as important pharmacogenomic markers associated with the severe 5-FU toxicity in Japanese population.

A Pilot Study for Return of Individual Pharmacogenomic Results to Population-Based Cohort Study Participants

Introduction:

Pharmacogenomic (PGx) testing results provide valuable information on drug selection and appropriate dosing, maximization of efficacy, and minimization of adverse effects. Although the number of large-scale, next-generation-sequencing-based PGx studies has recently increased, little is known about the risks and benefits of returning PGx results to ostensibly healthy individuals in research settings.

Methods:

Single-nucleotide variants of three actionable PGx genes, namely, MT-RNR1, CYP2C19, and NUDT15, were returned to 161 participants in a population-based Tohoku Medical Megabank project. Informed consent was obtained from the participants after a seminar on the outline of this study. The results were sent by mail alongside sealed information letter intended for clinicians. As an exception, genetic counseling was performed for the MT-RNR1 m.1555A > G variant carriers by a medical geneticist, and consultation with an otolaryngologist was encouraged. Questionnaire surveys (QSs) were conducted five times to evaluate the participants’ understanding of the topic, psychological impact, and attitude toward the study.

Results:

Whereas the majority of participants were unfamiliar with the term PGx, and none had undergone PGx testing before the study, more than 80% of the participants felt that they could acquire basic PGx knowledge sufficient to understand their genomic results and were satisfied with their potential benefit and use in future prescriptions. On the other hand, some felt that the PGx concepts or terminology was difficult to fully understand and suggested that in-person return of the results was desirable.

Conclusions:

These results collectively suggest possible benefits of returning preemptive PGx information to ostensibly healthy cohort participants in a research setting.

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

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

Comparison of Kit-Based Metabolomics with Other Methodologies in a Large Cohort, towards Establishing Reference Values

Metabolic profiling is an omics approach that can be used to observe phenotypic changes, making it particularly attractive for biomarker discovery. Although several candidate metabolites biomarkers for disease expression have been identified in recent clinical studies, the reference values of healthy subjects have not been established. In particular, the accuracy of concentrations measured by mass spectrometry (MS) is unclear. Therefore, comprehensive metabolic profiling in large-scale cohorts by MS to create a database with reference ranges is essential for evaluating the quality of the discovered biomarkers. In this study, we tested 8700 plasma samples by commercial kit-based metabolomics and separated them into two groups of 6159 and 2541 analyses based on the different ultra-high-performance tandem mass spectrometry (UHPLC-MS/MS) systems. We evaluated the quality of the quantified values of the detected metabolites from the reference materials in the group of 2541 compared with the quantified values from other platforms, such as nuclear magnetic resonance (NMR), supercritical fluid chromatography tandem mass spectrometry (SFC-MS/MS) and UHPLC-Fourier transform mass spectrometry (FTMS). The values of the amino acids were highly correlated with the NMR results, and lipid species such as phosphatidylcholines and ceramides showed good correlation, while the values of triglycerides and cholesterol esters correlated less to the lipidomics analyses performed using SFC-MS/MS and UHPLC-FTMS. The evaluation of the quantified values by MS-based techniques is essential for metabolic profiling in a large-scale cohort.

Japonica Array NEO with increased genome-wide coverage and abundant disease risk SNPs

Ethnic-specific SNP arrays are becoming more important to increase the power of genome-wide association studies in diverse population. In the Tohoku Medical Megabank Project, we have been developing a series of Japonica Arrays (JPA) for genotyping participants based on reference panels constructed from whole-genome sequence data of the Japanese population. Here, we designed a novel version of the SNP array for the Japanese population, called Japonica Array NEO (JPA NEO), comprising a total of 666,883 markers. Among them, 654,246 tag SNPs of autosomes and X chromosome were selected from an expanded reference panel of 3,552 Japanese, 3.5KJPNv2, using pairwise r² of linkage disequilibrium measures. Additionally, 28,298 markers were included for the evaluation of previously identified disease risk markers from the literature and databases, and those present in the Japanese population were extracted using the reference panel. Through genotyping 286 Japanese samples, we found that the imputation quality r² and INFO score in the minor allele frequency bin >2.5–5% were >0.9 and >0.8, respectively, and >12 million markers were imputed with an INFO score >0.8. From these results, JPA NEO is a promising tool for genotyping the Japanese population with genome-wide coverage, contributing to the development of genetic risk scores.

Construction and integration of three de novo Japanese human genome assemblies toward a population-specific reference

The complete human genome sequence is used as a reference for next-generation sequencing analyses. However, some ethnic ancestries are under-represented in the reference genome (e.g., GRCh37) due to its bias toward European and African ancestries. Here, we perform de novo assembly of three Japanese male genomes using > 100× Pacific Biosciences long reads and Bionano Genomics optical maps per sample. We integrate the genomes using the major allele for consensus and anchor the scaffolds using genetic and radiation hybrid maps to reconstruct each chromosome. The resulting genome sequence, JG1, is contiguous, accurate, and carries the Japanese major allele at most loci. We adopt JG1 as the reference for confirmatory exome re-analyses of seven rare-disease Japanese families and find that re-analysis using JG1 reduces total candidate variant calls versus GRCh37 while retaining disease-causing variants. These results suggest that integrating multiple genomes from a single population can aid genome analyses of that population.

jMorp updates in 2020: large enhancement of multi-omics data resources on the general Japanese population

In the Tohoku Medical Megabank project, genome and omics analyses of participants in two cohort studies were performed. A part of the data is available at the Japanese Multi Omics Reference Panel (jMorp; https://jmorp.megabank.tohoku.ac.jp) as a web-based database, as reported in our previous manuscript published in Nucleic Acid Research in 2018. At that time, jMorp mainly consisted of metabolome data; however, now genome, methylome, and transcriptome data have been integrated in addition to the enhancement of the number of samples for the metabolome data. For genomic data, jMorp provides a Japanese reference sequence obtained using de novo assembly of sequences from three Japanese individuals and allele frequencies obtained using whole-genome sequencing of 8,380 Japanese individuals. In addition, the omics data include methylome and transcriptome data from ∼300 samples and distribution of concentrations of more than 755 metabolites obtained using high-throughput nuclear magnetic resonance and high-sensitivity mass spectrometry. In summary, jMorp now provides four different kinds of omics data (genome, methylome, transcriptome, and metabolome), with a user-friendly web interface. This will be a useful scientific data resource on the general population for the discovery of disease biomarkers and personalized disease prevention and early diagnosis.

Functional Characterization of 40 CYP3A4 Variants by Assessing Midazolam 1'-Hydroxylation and Testosterone 6β-Hydroxylation

CYP3A4 is among the most abundant liver and intestinal drug-metabolizing cytochrome P450 enzymes, contributing to the metabolism of more than 30% of clinically used drugs. Therefore, interindividual variability in CYP3A4 activity is a frequent cause of reduced drug efficacy and adverse effects. In this study, we characterized wild-type CYP3A4 and 40 CYP3A4 variants, including 11 new variants, detected among 4773 Japanese individuals by assessing CYP3A4 enzymatic activities for two representative substrates (midazolam and testosterone). The reduced carbon monoxide-difference spectra of wild-type CYP3A4 and 31 CYP3A4 variants produced with our established mammalian cell expression system were determined by measuring the increase in maximum absorption at 450 nm after carbon monoxide treatment. The kinetic parameters of midazolam and testosterone hydroxylation by wild-type CYP3A4 and 29 CYP3A4 variants (K _m , k _cat , and catalytic efficiency) were determined, and the causes of their kinetic differences were evaluated by three-dimensional structural modeling. Our findings offer insight into the mechanism underlying interindividual differences in CYP3A4-dependent drug metabolism. Moreover, our results provide guidance for improving drug administration protocols by considering the information on CYP3A4 genetic polymorphisms. SIGNIFICANCE STATEMENT: CYP3A4 metabolizes more than 30% of clinically used drugs. Interindividual differences in drug efficacy and adverse-effect rates have been linked to ethnicity-specific differences in CYP3A4 gene variants in Asian populations, including Japanese individuals, indicating the presence of CYP3A4 polymorphisms resulting in the increased expression of loss-of-function variants. This study detected alterations in CYP3A4 activity due to amino acid substitutions by assessing the enzymatic activities of coding variants for two representative CYP3A4 substrates.

Enhancer remodeling promotes tumor-initiating activity in NRF2-activated non-small cell lung cancers

Transcriptional dysregulation, which can be caused by genetic and epigenetic alterations, is a fundamental feature of many cancers. A key cytoprotective transcriptional activator, NRF2, is often aberrantly activated in non-small cell lung cancers (NSCLCs) and supports both aggressive tumorigenesis and therapeutic resistance. Herein, we find that persistently activated NRF2 in NSCLCs generates enhancers at gene loci that are not normally regulated by transiently activated NRF2 under physiological conditions. Elevated accumulation of CEBPB in NRF2-activated NSCLCs is found to be one of the prerequisites for establishment of the unique NRF2-dependent enhancers, among which the NOTCH3 enhancer is shown to be critical for promotion of tumor-initiating activity. Enhancer remodeling mediated by NRF2-CEBPB cooperativity promotes tumor-initiating activity and drives malignancy of NRF2-activated NSCLCs via establishment of the NRF2-NOTCH3 regulatory axis.

Aberrant activation of NRF2 in cancer cells contributes to tumorigenicity and therapeutic resistance. Here, the authors show that NRF2 cooperates with CEBPB and remodels enhancers to confer tumor-initiating activity on NRF2- activated non-small cell lung cancers.

Identification of critical genetic variants associated with metabolic phenotypes of the Japanese population

We performed a metabolome genome-wide association study for the Japanese population in the prospective cohort study of Tohoku Medical Megabank. By combining whole-genome sequencing and nontarget metabolome analyses, we identified a large number of novel associations between genetic variants and plasma metabolites. Of the identified metabolite-associated genes, approximately half have already been shown to be involved in various diseases. We identified metabolite-associated genes involved in the metabolism of xenobiotics, some of which are from intestinal microorganisms, indicating that the identified genetic variants also markedly influence the interaction between the host and symbiotic bacteria. We also identified five associations that appeared to be female-specific. A number of rare variants that influence metabolite levels were also found, and combinations of common and rare variants influenced the metabolite levels more profoundly. These results support our contention that metabolic phenotyping provides important insights into how genetic and environmental factors provoke human diseases.

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

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

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

COXPRESdb v7: a gene coexpression database for 11 animal species supported by 23 coexpression platforms for technical evaluation and evolutionary inference

The advent of RNA-sequencing and microarray technologies has led to rapid growth of transcriptome data generated for a wide range of organisms, under various cellular, organ and individual conditions. Since the number of possible combinations of intercellular and extracellular conditions is almost unlimited, cataloging all transcriptome conditions would be an immeasurable challenge. Gene coexpression refers to the similarity of gene expression patterns under various conditions, such as disease states, tissue types, and developmental stages. Since the quality of gene coexpression data depends on the quality and quantity of transcriptome data, timely usage of the growing data is key to promoting individual research in molecular biology. COXPRESdb (http://coxpresdb.jp) is a database providing coexpression information for 11 animal species. One characteristic feature of COXPRESdb is its ability to compare multiple coexpression data derived from different transcriptomics technologies and different species, which strongly reduces false positive relationships in individual gene coexpression data. Here, we summarized the current version of this database, including 23 coexpression platforms with the highest-level quality till date. Using various functionalities in COXPRESdb, the new coexpression data would support a broader area of research from molecular biology to medical sciences.

jMorp: Japanese Multi Omics Reference Panel

We developed jMorp, a new database containing metabolome and proteome data for plasma obtained from >5000 healthy Japanese volunteers from the Tohoku Medical Megabank Cohort Study, which is available at https://jmorp.megabank.tohoku.ac.jp. Metabolome data were measured by proton nuclear magnetic resonance (NMR) and liquid chromatography–mass spectrometry (LC–MS), while proteome data were obtained by nanoLC–MS. We released the concentration distributions of 37 metabolites identified by NMR, distributions of peak intensities of 257 characterized metabolites by LC–MS, and observed frequencies of 256 abundant proteins. Additionally, correlation networks for the metabolites can be observed using an interactive network viewer. Compared with some existing databases, jMorp has some unique features: (i) Metabolome data were obtained using a single protocol in a single institute, ensuring that measurement biases were significantly minimized; (ii) The database contains large-scale data for healthy volunteers with various health records and genome data and (iii) Correlations between metabolites can be easily observed using the graphical viewer. Metabolites data are becoming important intermediate markers for evaluating the health states of humans, and thus jMorp is an outstanding resource for a wide range of researchers, particularly those in the fields of medical science, applied molecular biology, and biochemistry.

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

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

Pathogenic mutations identified by a multimodality approach in 117 Japanese Fanconi anemia patients

Fanconi anemia is a rare recessive disease characterized by multiple congenital abnormalities, progressive bone marrow failure, and a predisposition to malignancies. It results from mutations in one of the 22 known FANC genes. The number of Japanese Fanconi anemia patients with a defined genetic diagnosis was relatively limited. In this study, we reveal the genetic subtyping and the characteristics of mutated FANC genes in Japan and clarify the genotype-phenotype correlations. We studied 117 Japanese patients and successfully subtyped 97% of the cases. FANCA and FANCG pathogenic variants accounted for the disease in 58% and 25% of Fanconi anemia patients, respectively. We identified one FANCA and two FANCG hot spot mutations, which are found at low percentages (0.04-0.1%) in the whole-genome reference panel of 3,554 Japanese individuals (Tohoku Medical Megabank). FANCB was the third most common complementation group and only one FANCC case was identified in our series. Based on the data from the Tohoku Medical Megabank, we estimate that approximately 2.6% of Japanese are carriers of disease-causing FANC gene variants, excluding missense mutations. This is the largest series of subtyped Japanese Fanconi anemia patients to date and the results will be useful for future clinical management.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

ATTED-II in 2018: A Plant Coexpression Database Based on Investigation of the Statistical Property of the Mutual Rank Index

ATTED-II (http://atted.jp) is a coexpression database for plant species to aid in the discovery of relationships of unknown genes within a species. As an advanced coexpression analysis method, multispecies comparisons have the potential to detect alterations in gene relationships within an evolutionary context. However, determining the validity of comparative coexpression studies is difficult without quantitative assessments of the quality of coexpression data. ATTED-II (version 9) provides 16 coexpression platforms for nine plant species, including seven species supported by both microarray- and RNA sequencing (RNAseq)-based coexpression data. Two independent sources of coexpression data enable the assessment of the reproducibility of coexpression. The latest coexpression data for Arabidopsis (Ath-m.c7-1 and Ath-r.c3-0) showed the highest reproducibility (Jaccard coefficient = 0.13) among previous coexpression data in ATTED-II. We also investigated the statistical basis of the mutual rank (MR) index as a coexpression measure by bootstrap sampling of experimental units. We found that the error distribution of the logit-transformed MR index showed normality with equal variances for each coexpression platform. Because the MR error was strongly correlated with the number of samples for the coexpression data, typical confidence intervals for the MR index can be estimated for any coexpression platform. These new, high-quality coexpression data can be analyzed with any tool in ATTED-II and combined with external resources to obtain insight into plant biology.

NCMine: Core-peripheral based functional module detection using near-clique mining

Motivation: The identification of functional modules from protein–protein interaction (PPI) networks is an important step toward understanding the biological features of PPI networks. The detection of functional modules in PPI networks is often performed by identifying internally densely connected subnetworks, and often produces modules with “core” and “peripheral” proteins. The core proteins are the ones having dense connections to each other in a module. The difference between core and peripheral proteins is important to understand the functional roles of proteins in modules, but there are few methods to explicitly elucidate the internal structure of functional modules at gene level.

Results: We propose NCMine, which is a novel network clustering method and visualization tool for the core-peripheral structure of functional modules. It extracts near-complete subgraphs from networks based on a node-weighting scheme using degree centrality, and reports subgroups as functional modules. We implemented this method as a plugin of Cytoscape, which is widely used to visualize and analyze biological networks. The plugin allows users to extract functional modules from PPI networks and interactively filter modules of interest. We applied the method to human PPI networks, and found several examples with the core-peripheral structure of modules that may be related to cancer development.

Availability and Implementation: The Cytoscape plugin and tutorial are available at Cytoscape AppStore. (http://apps.cytoscape.org/apps/ncmine).

Contact:kengo@ecei.tohoku.ac.jp

Supplementary information:Supplementary data are available at Bioinformatics online.

ATTED-II in 2016: A Plant Coexpression Database Towards Lineage-Specific Coexpression

ATTED-II (http://atted.jp) is a coexpression database for plant species with parallel views of multiple coexpression data sets and network analysis tools. The user can efficiently find functional gene relationships and design experiments to identify gene functions by reverse genetics and general molecular biology techniques. Here, we report updates to ATTED-II (version 8.0), including new and updated coexpression data and analysis tools. ATTED-II now includes eight microarray- and six RNA sequencing-based coexpression data sets for seven dicot species (Arabidopsis, field mustard, soybean, barrel medick, poplar, tomato and grape) and two monocot species (rice and maize). Stand-alone coexpression analyses tend to have low reliability. Therefore, examining evolutionarily conserved coexpression is a more effective approach from the viewpoints of reliability and evolutionary importance. In contrast, the reliability of species-specific coexpression data remains poor. Our assessment scores for individual coexpression data sets indicated that the quality of the new coexpression data sets in ATTED-II is higher than for any previous coexpression data set. In addition, five species (Arabidopsis, soybean, tomato, rice and maize) in ATTED-II are now supported by both microarray- and RNA sequencing-based coexpression data, which has increased the reliability. Consequently, ATTED-II can now provide lineage-specific coexpression information. As an example of the use of ATTED-II to explore lineage-specific coexpression, we demonstrate monocot- and dicot-specific coexpression of cell wall genes. With the expanded coexpression data for multilevel evaluation, ATTED-II provides new opportunities to investigate lineage-specific evolution in plants.

COXPRESdb in 2015: coexpression database for animal species by DNA-microarray and RNAseq-based expression data with multiple quality assessment systems

The COXPRESdb (http://coxpresdb.jp) provides gene coexpression relationships for animal species. Here, we report the updates of the database, mainly focusing on the following two points. For the first point, we added RNAseq-based gene coexpression data for three species (human, mouse and fly), and largely increased the number of microarray experiments to nine species. The increase of the number of expression data with multiple platforms could enhance the reliability of coexpression data. For the second point, we refined the data assessment procedures, for each coexpressed gene list and for the total performance of a platform. The assessment of coexpressed gene list now uses more reasonable P-values derived from platform-specific null distribution. These developments greatly reduced pseudo-predictions for directly associated genes, thus expanding the reliability of coexpression data to design new experiments and to discuss experimental results.

BioHackathon series in 2011 and 2012: penetration of ontology and linked data in life science domains

The application of semantic technologies to the integration of biological data and the interoperability of bioinformatics analysis and visualization tools has been the common theme of a series of annual BioHackathons hosted in Japan for the past five years. Here we provide a review of the activities and outcomes from the BioHackathons held in 2011 in Kyoto and 2012 in Toyama. In order to efficiently implement semantic technologies in the life sciences, participants formed various sub-groups and worked on the following topics: Resource Description Framework (RDF) models for specific domains, text mining of the literature, ontology development, essential metadata for biological databases, platforms to enable efficient Semantic Web technology development and interoperability, and the development of applications for Semantic Web data. In this review, we briefly introduce the themes covered by these sub-groups. The observations made, conclusions drawn, and software development projects that emerged from these activities are discussed.

ATTED-II in 2014: evaluation of gene coexpression in agriculturally important plants

ATTED-II (http://atted.jp) is a database of coexpressed genes that was originally developed to identify functionally related genes in Arabidopsis and rice. Herein, we describe an updated version of ATTED-II, which expands this resource to include additional agriculturally important plants. To improve the quality of the coexpression data for Arabidopsis and rice, we included more gene expression data from microarray and RNA sequencing studies. The RNA sequencing-based coexpression data now cover 94% of the Arabidopsis protein-encoding genes, representing a substantial increase from previously available microarray-based coexpression data (76% coverage). We also generated coexpression data for four dicots (soybean, poplar, grape and alfalfa) and one monocot (maize). As both the quantity and quality of expression data for the non-model species are generally poorer than for the model species, we verified coexpression data associated with these new species using multiple methods. First, the overall performance of the coexpression data was evaluated using gene ontology annotations and the coincidence of a genomic feature. Secondly, the reliability of each guide gene was determined by comparing coexpressed gene lists between platforms. With the expanded and newly evaluated coexpression data, ATTED-II represents an important resource for identifying functionally related genes in agriculturally important plants.

COXPRESdb: a database of comparative gene coexpression networks of eleven species for mammals

Coexpressed gene databases are valuable resources for identifying new gene functions or functional modules in metabolic pathways and signaling pathways. Although coexpressed gene databases are a fundamental platform in the field of plant biology, their use in animal studies is relatively limited. The COXPRESdb (http://coxpresdb.jp) provides coexpression relationships for multiple animal species, as comparisons of coexpressed gene lists can enhance the reliability of gene coexpression determinations. Here, we report the updates of the database, mainly focusing on the following two points. First, we updated our coexpression data by including recent microarray data for the previous seven species (human, mouse, rat, chicken, fly, zebrafish and nematode) and adding four new species (monkey, dog, budding yeast and fission yeast), along with a new human microarray platform. A reliability scoring function was also implemented, based on coexpression conservation to filter out coexpression with low reliability. Second, the network drawing function was updated, to implement automatic cluster analyses with enrichment analyses in Gene Ontology and in cis elements, along with interactive network analyses with Cytoscape Web. With these updates, COXPRESdb will become a more powerful tool for analyses of functional and regulatory networks of genes in a variety of animal species.