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.

Variants of Flavin-Containing Monooxygenase 3 Found in Subjects in an Updated Database of Genome Resources.

Single-nucleotide substitutions of human flavin-containing monooxygenase 3 (FMO3) identified in the whole-genome sequences of the updated Japanese population reference panel (now containing 38,000 subjects) were investigated. In this study, two stop codon mutations, two frameshifts, and 43 amino-acid-substituted FMO3 variants were identified. Among these 47 variants, one stop codon mutation, one frameshift, and 24 substituted variants were already recorded in the National Center for Biotechnology Information database. Functionally impaired FMO3 variants are known to be associated with the metabolic disorder trimethylaminuria; consequently, the enzymatic activities of the 43 substituted FMO3 variants were investigated. Twenty-seven recombinant FMO3 variants expressed in bacterial membranes had similar activities toward trimethylamine N-oxygenation (∼75%-125%) to that of wild-type FMO3 (98 minutes-1). However, six recombinant FMO3 variants (Arg51Gly, Val283Ala, Asp286His, Val382Ala, Arg387His, and Phe451Leu) had moderately decreased (∼50%) activities toward trimethylamine N-oxygenation, and 10 recombinant FMO3 variants (Gly11Asp, Gly39Val, Met66Lys, Asn80Lys, Val151Glu, Gly193Arg, Arg387Cys, Thr453Pro, Leu457Trp, and Met497Arg) showed severely decreased FMO3 catalytic activity (<10%). Because of the known deleterious effects of FMO3 C-terminal stop codons, the four truncated FMO3 variants (Val187SerfsTer25, Arg238Ter, Lys416SerfsTer72, and Gln427Ter) were suspected to be inactive with respect to trimethylamine N-oxygenation. The FMO3 p.Gly11Asp and p.Gly193Arg variants were located within the conserved sequences of flavin adenine dinucleotide (positions 9-14) and NADPH (positions 191-196) binding sites, which are important for FMO3 catalytic function. Whole-genome sequence data and kinetic analyses revealed that 20 of the 47 nonsense or missense FMO3 variants had moderately or severely impaired activity toward N-oxygenation of trimethylaminuria. SIGNIFICANCE STATEMENT: The number of single-nucleotide substitutions in human flavin-containing monooxygenase 3 (FMO3) recorded in the expanded Japanese population reference panel database was updated. One stop mutation, FMO3 p.Gln427Ter; one frameshift (p.Lys416SerfsTer72); and 19 novel amino-acid-substituted FMO3 variants were identified, along with p.Arg238Ter, p.Val187SerfsTer25, and 24 amino-acid-substituted variants already recorded with reference SNP (rs) numbers. Recombinant FMO3 Gly11Asp, Gly39Val, Met66Lys, Asn80Lys, Val151Glu, Gly193Arg, Arg387Cys, Thr453Pro, Leu457Trp, and Met497Arg variants showed severely decreased FMO3 catalytic activity, possibly associated with the trimethylaminuria.

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, Km and Vmax , and intrinsic clearance (CLint = Vmax /Km ) of the wild-type DHPase and eight variants were measured. Five of these variants (R118Q, H295R, T418I, Y448H, and T513A) showed significantly reduced CLint 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.

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.

Risk analysis of COVID-19 hospitalization and critical care by race and region in the United States: a cohort study.

BACKGROUND: This study aimed to identify the current risk factors for coronavirus disease 2019 severity and examine its association with medication use. METHODS: We used data from a large United States electronic health record database to conduct an anonymized cohort study of 171,491 patients with coronavirus disease 2019. The study was conducted from January 1, 2020, to August 27, 2021. Data on age, race, sex, history of diseases, and history of medication prescriptions were analyzed using the Cox proportional hazards model analysis to calculate hazard ratios for hospitalization and severe risk. RESULTS: Factors that increased the risk of hospitalization and critical care were age ≥ 65 years, male sex, type 2 diabetes, hypertension, interstitial pneumonia, and cardiovascular disease. In particular, age ≥ 65 years significantly increased the risk of hospitalization (hazard ratio, 2.81 [95% confidence interval, 2.58-3.07]; P < 0.001) and critical care (hazard ratio, 3.45 [2.88-4.14]; P < 0.001). In contrast, patients with hyperlipidemia had a reduced risk. However, patients with hyperlipidemia who were not taking statins had a significantly increased risk of hospitalization (hazard ratio, 1.24 [1.16-1.34]; P < 0.001). Sodium-glucose cotransporter-2 inhibitors, angiotensin-converting enzyme inhibitors, angiotensin II receptor blockers, glucocorticoids, and statins significantly reduced the risk of hospitalization and critical care. The risk of hospitalization and critical care increased in patients of all ethnicities with type 2 diabetes. The factors that significantly increased the risk of hospitalization in all regions were older age, hypertension, chronic obstructive pulmonary disease, and cardiovascular disease. CONCLUSION: This study identified factors that increase or reduce the risk of severe coronavirus disease. The provision of appropriate drug treatment and modification of lifestyle-related risk factors may reduce coronavirus disease severity.

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.

Impacts of NRF2 activation in non-small-cell lung cancer cell lines on extracellular metabolites.

Aberrant activation of NRF2 is as a critical prognostic factor that drives the malignant progression of various cancers. Cancer cells with persistent NRF2 activation heavily rely on NRF2 activity for therapeutic resistance and aggressive tumorigenic capacity. To clarify the metabolic features of NRF2-activated lung cancers, we conducted targeted metabolomic (T-Met) and global metabolomic (G-Met) analyses of non-small-cell lung cancer (NSCLC) cell lines in combination with exome and transcriptome analyses. Exome analysis of 88 cell lines (49 adenocarcinoma, 14 large cell carcinoma, 15 squamous cell carcinoma and 10 others) identified non-synonymous mutations in the KEAP1, NRF2 and CUL3 genes. Judging from the elevated expression of NRF2 target genes, these mutations are expected to result in the constitutive stabilization of NRF2. Out of the 88 cell lines, 52 NSCLC cell lines (29 adenocarcinoma, 10 large cell carcinoma, 9 squamous cell carcinoma and 4 others) were subjected to T-Met analysis. Classification of the 52 cell lines into three groups according to the NRF2 target gene expression enabled us to draw typical metabolomic signatures induced by NRF2 activation. From the 52 cell lines, 18 NSCLC cell lines (14 adenocarcinoma, 2 large cell carcinoma, 1 squamous cell carcinoma and 1 others) were further chosen for G-Met and detailed transcriptome analyses. G-Met analysis of their culture supernatants revealed novel metabolites associated with NRF2 activity, which may be potential diagnostic biomarkers of NRF2 activation. This study also provides useful information for the exploration of new metabolic nodes for selective toxicity towards NRF2-activated NSCLC.

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.

Identification of somatic genetic alterations in ovarian clear cell carcinoma with next generation sequencing.

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

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.

Iron attenuates erythropoietin production by decreasing hypoxia-inducible transcription factor 2α concentrations in renal interstitial fibroblasts.

Iron is an essential mineral for oxygen delivery and for a variety of enzymatic activities, but excessive iron results in oxidative cytotoxicity. Because iron is primarily used in red blood cells, defective erythropoiesis caused by loss of the erythroid growth factor erythropoietin (Epo) elevates iron storage levels in serum and tissues. Here, we investigated the effects of iron in a mouse model of Epo-deficiency anemia, in which serum iron concentration was significantly elevated. We found that intraperitoneal injection of iron-dextran caused severe iron deposition in renal interstitial fibroblasts, the site of Epo production. Iron overload induced by either intraperitoneal injection or feeding decreased activity of endogenous Epo gene expression by reducing levels of hypoxia-inducible transcription factor 2α (HIF2α), the major transcriptional activator of the Epo gene. Administration of an iron-deficient diet to the anemic mice reduced serum iron to normal concentration and enhanced the ability of renal Epo production. These results demonstrate that iron overload due to Epo deficiency attenuates endogenous Epo gene expression in the kidneys. Thus, iron suppresses Epo production by reducing HIF2α concentration in renal interstitial fibroblasts.

Evaluation of reported pathogenic variants and their frequencies in a Japanese population based on a whole-genome reference panel of 2049 individuals.

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

Functional Characterization of 21 Allelic Variants of Dihydropyrimidine Dehydrogenase Identified in 1070 Japanese Individuals.

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

HLA-VBSeq: accurate HLA typing at full resolution from whole-genome sequencing data.

BACKGROUND: Human leucocyte antigen (HLA) genes play an important role in determining the outcome of organ transplantation and are linked to many human diseases. Because of the diversity and polymorphisms of HLA loci, HLA typing at high resolution is challenging even with whole-genome sequencing data. RESULTS: We have developed a computational tool, HLA-VBSeq, to estimate the most probable HLA alleles at full (8-digit) resolution from whole-genome sequence data. HLA-VBSeq simultaneously optimizes read alignments to HLA allele sequences and abundance of reads on HLA alleles by variational Bayesian inference. We show the effectiveness of the proposed method over other methods through the analysis of predicting HLA types for HLA class I (HLA-A, -B and -C) and class II (HLA-DQA1,-DQB1 and -DRB1) loci from the simulation data of various depth of coverage, and real sequencing data of human trio samples. CONCLUSIONS: HLA-VBSeq is an efficient and accurate HLA typing method using high-throughput sequencing data without the need of primer design for HLA loci. Moreover, it does not assume any prior knowledge about HLA allele frequencies, and hence HLA-VBSeq is broadly applicable to human samples obtained from a genetically diverse population.

Comprehensive transcriptomic analysis of molecularly targeted drugs in cancer for target pathway evaluation.

Targeted therapy is a rational and promising strategy for the treatment of advanced cancer. For the development of clinical agents targeting oncogenic signaling pathways, it is important to define the specificity of compounds to the target molecular pathway. Genome-wide transcriptomic analysis is an unbiased approach to evaluate the compound mode of action, but it is still unknown whether the analysis could be widely applicable to classify molecularly targeted anticancer agents. We comprehensively obtained and analyzed 129 transcriptomic datasets of cancer cells treated with 83 anticancer drugs or related agents, covering most clinically used, molecularly targeted drugs alongside promising inhibitors of molecular cancer targets. Hierarchical clustering and principal component analysis revealed that compounds targeting similar target molecules or pathways were clustered together. These results confirmed that the gene signatures of these drugs reflected their modes of action. Of note, inhibitors of oncogenic kinase pathways formed a large unique cluster, showing that these agents affect a shared molecular pathway distinct from classical antitumor agents and other classes of agents. The gene signature analysis further classified kinome-targeting agents depending on their target signaling pathways, and we identified target pathway-selective signature gene sets. The gene expression analysis was also valuable in uncovering unexpected target pathways of some anticancer agents. These results indicate that comprehensive transcriptomic analysis with our database (http://scads.jfcr.or.jp/db/cs/) is a powerful strategy to validate and re-evaluate the target pathways of anticancer compounds.

Validation of multiple single nucleotide variation calls by additional exome analysis with a semiconductor sequencer to supplement data of whole-genome sequencing of a human population.

BACKGROUND: Validation of single nucleotide variations in whole-genome sequencing is critical for studying disease-related variations in large populations. A combination of different types of next-generation sequencers for analyzing individual genomes may be an efficient means of validating multiple single nucleotide variations calls simultaneously. RESULTS: Here, we analyzed 12 independent Japanese genomes using two next-generation sequencing platforms: the Illumina HiSeq 2500 platform for whole-genome sequencing (average depth 32.4×), and the Ion Proton semiconductor sequencer for whole exome sequencing (average depth 109×). Single nucleotide polymorphism (SNP) calls based on the Illumina Human Omni 2.5-8 SNP chip data were used as the reference. We compared the variant calls for the 12 samples, and found that the concordance between the two next-generation sequencing platforms varied between 83% and 97%. CONCLUSIONS: Our results show the versatility and usefulness of the combination of exome sequencing with whole-genome sequencing in studies of human population genetics and demonstrate that combining data from multiple sequencing platforms is an efficient approach to validate and supplement SNP calls.

Rare but impaired flavin-containing monooxygenase 3 (FMO3) variants reported in a recently updated Japanese mega-databank of genome resources.

Genetic variants of human flavin-containing monooxygenase 3 (FMO3) were investigated using an updated Japanese population panel containing 54,000 subjects (the previous panel contained 38,000 subjects). One stop codon mutation and six amino acid-substituted FMO3 variants were newly identified in the updated databank. Of these, two substituted variants (p.Thr329Ala and p.Arg492Trp) were previously identified in compound haplotypes with p.[(Glu158Lys; Glu308Gly)] and were associated with the metabolic disorder trimethylaminuria. Three recombinant FMO3 protein variants (p.Ser137Leu, p.Ala334Val, and p.Ile426Val) expressed in bacterial membranes had similar activities toward trimethylamine N-oxygenation (∼75-125 %) as wild-type FMO3 (117 min-1); however, the recombinant novel FMO3 variant Phe313Ile showed moderately decreased FMO3 catalytic activity (∼20 % of wild-type). Because of the known deleterious effects of FMO3 C-terminal stop codons, the novel truncated FMO3 Gly184Ter variant was suspected to be inactive. To easily identify the four impaired FMO3 variants (one stop codon mutation and three amino-acid substitutions) in the clinical setting, simple confirmation methods for these FMO3 variants are proposed using polymerase chain reaction/restriction fragment length polymorphism or allele-specific PCR methods. The updated whole-genome sequence data and kinetic analyses revealed that four of the seven single-nucleotide nonsense or missense FMO3 variants had moderately or severely impaired activity toward trimethylamine N-oxygenation.

Development of a gene expression database and related analysis programs for evaluation of anticancer compounds.

Genome-wide transcriptional expression analysis is a powerful strategy for characterizing the biological activity of anticancer compounds. It is often instructive to identify gene sets involved in the activity of a given drug compound for comparison with different compounds. Currently, however, there is no comprehensive gene expression database and related application system that is; (i) specialized in anticancer agents; (ii) easy to use; and (iii) open to the public. To develop a public gene expression database of antitumor agents, we first examined gene expression profiles in human cancer cells after exposure to 35 compounds including 25 clinically used anticancer agents. Gene signatures were extracted that were classified as upregulated or downregulated after exposure to the drug. Hierarchical clustering showed that drugs with similar mechanisms of action, such as genotoxic drugs, were clustered. Connectivity map analysis further revealed that our gene signature data reflected modes of action of the respective agents. Together with the database, we developed analysis programs that calculate scores for ranking changes in gene expression and for searching statistically significant pathways from the Kyoto Encyclopedia of Genes and Genomes database in order to analyze the datasets more easily. Our database and the analysis programs are available online at our website (http://scads.jfcr.or.jp/db/cs/). Using these systems, we successfully showed that proteasome inhibitors are selectively classified as endoplasmic reticulum stress inducers and induce atypical endoplasmic reticulum stress. Thus, our public access database and related analysis programs constitute a set of efficient tools to evaluate the mode of action of novel compounds and identify promising anticancer lead compounds.

Detection and Correction of Sample Misidentifications in a Biobank Using the MassARRAY System and Genomic Information.

With the number of samples increasing in many biobanks, one of the most pressing tasks is recording the correct relationships between information and the specimens. Genomic information is useful in determining the identity of these specimens. The Tohoku Medical Megabank Organization is running one of the largest biobanks in Japan. Here, we introduce a management system, which includes the development of a new probe set for the MassARRAY system for use during the production of proliferating T cells (T cells) and lymphoblastoid cell lines (LCLs). We selected single nucleotide variants that could be detected by next-generation sequencing and showed high resolution with ∼0.5 minor allele frequencies. After checking the set of probes against 96 samples from 48 people, we obtained no contradictory results in comparison with our genome sequence information. When we applied the set to our 3035 LCLs and 2256 T cells, the result showed 98.93% consistency with the corresponding genomic information. We surveyed the handling records of the 1.07% of samples that showed inconsistencies, and found that most had resulted from human errors (ID swapping between samples) during manual operations. After improving a few error-prone protocols, the error rate dropped to 0.47% for LCLs and 0% for T cells. Overall, the system that we developed shows high accuracy with easy and fast operability, and provides a good opportunity to improve the validation procedure to facilitate high-quality banking, especially in cases involving genomic information.