Gene Rearrangement and Expression of PRKACA and PRKACB Govern Morphobiology of Pancreatobiliary Oncocytic Neoplasms.

Intraductal oncocytic papillary neoplasms (IOPNs) are distinct from intraductal papillary mucinous neoplasms based on characteristic morphologic and genetic features represented by fusion genes involving PRKACA or PRKACB (PRKACA/B). However, pancreatic and biliary tumors with partial oncocytic features are often encountered clinically, and their molecular features are yet to be clarified. This study included 80 intraductal papillary neoplasms: 32 tumors with mature IOPN morphology (typical), 28 with partial or subclonal oncocytic features (atypical), and 20 without oncocytic features (control). We analyzed PRKACA/B fusion genes, including ATP1B1::PRKACA, DNAJB1::PRKACA, and ATP1B1::PRKACB, by reverse-transcription PCR; mRNA expression of fusion genes and nonrearranged PRKACA/B genes by quantitative reverse-transcription PCR; mutations in KRAS, BRAF, and GNAS by targeted sequencing or droplet digital PCR; and the expression of cyclic adenosine monophosphate (cAMP)-dependent protein kinase catalytic subunits α (PRKACA) and ß (PRKACB), phosphorylated cAMP response element-binding protein, and aberrations of p16, p53, SMAD4, STK11, and ß-catenin by immunohistochemistry. PRKACA/B fusion genes were detected in 100% (32/32) of typical, 46% (13/28) of atypical, and 0% (0/20) of control (P < .05). Expression of PRKACA, PRKACB, and phosphorylated cAMP response element-binding protein was upregulated in neoplasms with PRKACA/B fusion genes (P < .05). mRNA expression of the PRKACA/B fusion genes and protein expression of PRKACA or PRKACB tended to be higher in typical than in atypical cases (mRNA, P = .002; protein expression, P = .054). In some atypical neoplasms with mixed subtypes, PRKACA/B fusion genes were superimposed exclusively on oncocytic components. Typical IOPNs harbored fewer KRAS and GNAS mutations than control samples and fewer alterations in p53 and STK11 than atypical samples (P < .05). In conclusion, PRKACA/B fusion genes not only are the characteristic drivers of IOPNs but also play a crucial role in the development of subclonal oncocytic neoplasms. Moreover, oncocytic morphology is strongly associated with upregulation of PRKACA/B, which may provide clues for potential therapeutic options.

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.

Morphometric analysis of nuclear shape irregularity as a novel predictor of programmed death-ligand 1 expression in lung squamous cell carcinoma.

Immune checkpoint inhibitor (ICI) therapy has been established as one of the key treatment strategies for lung squamous cell carcinoma (LUSQ). The status of programmed death-ligand 1 (PD-L1) in tumor cells and/or immune cells using immunohistochemistry has been primarily used as a surrogate marker for determining ICI treatment; however, when the tissues to be examined are small, false-negative results could be unavoidable due to the heterogeneity of PD-L1 immunoreactivity. To overcome this practical limitation, we attempted to explore the status of nuclear atypia evaluated using morphometry as a potential predictor of PD-L1 status in LUSQ. We correlated the parameters related to nuclear atypia with PD-L1 status using two different cohorts of LUSQ patients (95 cases from The Cancer Genome Atlas database and 30 cases from the Miyagi Cancer Center). Furthermore, we studied the gene mutation status to elucidate the genetic profile of PD-L1 predictable cases. The results revealed that nuclear atypia, especially morphometric parameters related to nuclear shape irregularity, including aspect ratio, circularity, roundness, and solidity, were all significantly associated with PD-L1 status. Additionally, LUSQ cases with high PD-L1 expression and pronounced nuclear atypia were significantly associated with C10orf71 and COL14A1 mutations compared with those with low PD-L1 expression and mild nuclear atypia. We demonstrated for the first time that nuclear shape irregularity could represent a novel predictor of PD-L1 expression in LUSQ. Including the morphometric parameters related to nuclear atypia in conjunction with PD-L1 status could help determine an effective ICI therapeutic strategy; however, further investigation is required.

Salivary gland cancer organoids are valid for preclinical genotype-oriented medical precision trials.

Salivary gland cancers (SGCs) are heterogeneous tumors, and precision oncology represents a promising therapeutic approach; however, its impact on SGCs remains obscure. This study aimed to establish a translational model for testing molecular-targeted therapies by combining patient-derived organoids and genomic analyses of SGCs. We enrolled 29 patients, including 24 with SGCs and 5 with benign tumors. Resected tumors were subjected to organoid and monolayer cultures, as well as whole-exome sequencing. Organoid and monolayer cultures of SGCs were successfully established in 70.8% and 62.5% of cases, respectively. Organoids retained most histopathological and genetic profiles of their original tumors. In contrast, 40% of the monolayer-cultured cells did not harbor somatic mutations of their original tumors. The efficacy of molecular-targeted drugs tested on organoids depended on their oncogenic features. Organoids recapitulated the primary tumors and were useful for testing genotype-oriented molecular targeted therapy, which is valuable for precision medicine in patients with SGCs.

Plasma metabolic disturbances during pregnancy and postpartum in women with depression.

Examining plasma metabolic profiling during pregnancy and postpartum could help clinicians understand the risk factors for postpartum depression (PPD) development. This analysis targeted paired plasma metabolites in mid-late gestational and 1 month postpartum periods in women with (n = 209) or without (n = 222) PPD. Gas chromatogram-mass spectrometry was used to analyze plasma metabolites at these two time points. Among the 170 objected plasma metabolites, principal component analysis distinguished pregnancy and postpartum metabolites but failed to discriminate women with and without PPD. Compared to women without PPD, those with PPD exhibited 37 metabolites with disparate changes during pregnancy and the 1-month postpartum period and an enriched citrate cycle. Machine learning and multivariate statistical analysis identified two or three compounds that could be potential biomarkers for PPD prediction during pregnancy. Our findings suggest metabolic disturbances in women with depression and may help to elucidate metabolic processes associated with PPD development.

X chromosome aneuploidies and schizophrenia: association analysis and phenotypic characterization.

AIM: The aims of the present study were: (i) to examine the association between schizophrenia (SCZ) and 47, XXY or 47, XXX in a large case-control sample; and (ii) to characterize the clinical features of patients with SCZ with these X chromosome aneuploidies. METHODS: To identify 47, XXY and 47, XXX, array comparative genomic hybridization (aCGH) was performed in 3188 patients with SCZ and 3586 controls. We examined the association between 47, XXY and 47, XXX and SCZ in males and females separately using exact conditional tests to control for platform effects. Clinical data were retrospectively examined for patients with SCZ with X chromosome aneuploidies. RESULTS: Of the analyzed samples, 3117 patients (97.8%) and 3519 controls (98.1%) passed our quality control. X chromosome aneuploidies were exclusively identified in patients: 47, XXY in seven patients (0.56%), 47, XXX in six patients (0.42%). Statistical analysis revealed a significant association between SCZ and 47, XXY (P = 0.028) and 47, XXX (P = 0.011). Phenotypic data were available from 12 patients. Treatment-resistance to antipsychotics and manic symptoms were observed in six patients each (four with 47, XXY and two with 47, XXX for both), respectively. Statistical analysis revealed that treatment-resistance to antipsychotics, mood stabilizer use, and manic symptoms were significantly more common in patients with 47, XXY than in male patients without pathogenic copy number variations. CONCLUSION: These findings indicate that both 47, XXY and 47, XXX are significantly associated with risk for SCZ. Patients with SCZ with 47, XXY may be characterized by treatment-resistance and manic symptoms.

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.

The association between ERK inhibitor sensitivity and molecular characteristics in colorectal cancer.

The mitogen-activated protein kinase (MAPK) pathway plays an important role in the colorectal cancer (CRC) progression, being supposed to be activated by the gene mutations, such as BRAF or KRAS. Although the inhibitors of extracellular signal-regulated kinase (ERK) have demonstrated efficacy in the cells with the BRAF or KRAS mutations, a clinical response is not always associated with the molecular signature. The patient-derived organoids (PDO) have emerged as a powerful in vitro model system to study cancer, and it has been widely applied for the drug screening. The present study aims to analyze the association between the molecular characteristics which analyzed by next-generation sequencing (NGS) and sensitivity to the ERK inhibitor (i.e., SCH772984) in PDO derived from CRC specimens. A drug sensitivity test for the SCH772984 was conducted using 14 CRC cell lines, and the results demonstrated that the sensitivity was in agreement with the BRAF mutation, but was not completely consistent with the KRAS status. In the drug sensitivity test for PDO, 6 out of 7 cases with either BRAF or KRAS mutations showed sensitivity to the SCH772984, while 5 out of 6 cases of both BRAF and KRAS wild-types were resistant. The results of this study suggested that the molecular status of the clinical specimens are likely to represent the sensitivity in the PDOs but is not necessarily absolutely overlapping. PDO might be able to complement the limitations of the gene panel and have the potential to provide a novel precision medicine.

Development of a system combining comprehensive genotyping and organoid cultures for identifying and testing genotype-oriented personalised medicine for pancreatobiliary cancers.

BACKGROUND: Pancreatobiliary cancer is a highly aggressive tumour with a dismal prognosis. Personalised medicine represents a promising and effective therapeutic approach for this intractable disease. In this study, we aimed to establish a system for identifying and testing genotype-oriented targeted drugs for pancreatobiliary cancers by combining exome sequencing and organoid culture of primary tumours. METHODS: Tumour cells isolated from resected tumours were subjected to organoid cultures based on published protocols with modifications. Exome sequencing was performed on the primary tumours. Histopathological and molecular features of the primary tumours were validated in the corresponding organoids. Genotype-oriented candidate targeted drugs were identified from exome sequencing, and their efficacies were tested in the organoids. RESULTS: Organoid cultures succeeded in 30 of 54 (55.6%) cases. Six primary cancers of the biliary tract and gall bladder were subjected to exome sequencing, which revealed a variety of somatic mutations of genes involved in signalling pathways, epigenetic modifiers, genome maintenance and metabolic enzymes. Most of the organoids of these 6 cases showed identical histopathological features and genomic aberrations as those of the primary tumours. Some of the aberrations were candidates for targeted therapies. Integrin-linked kinase (ILK) was one such candidate target, and an ILK inhibitor was confirmed to suppress proliferation of patient-derived organoids. CONCLUSIONS: By combining exome sequencing and organoid culture, our model enabled to identify genotype-oriented targets for personalised medicine and to test efficacies of candidate targeted drugs in the organoids. The current proof-of-concept approach could increase therapeutic opportunities for patients with pancreatobiliary cancers.

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.

Whole exome sequencing and establishment of an organoid culture of the carcinoma showing thymus-like differentiation (CASTLE) of the parotid gland.

Carcinoma showing thymus-like differentiation (CASTLE) is a rare tumor, especially in the parotid gland. We encountered a CASTLE of the parotid gland and analyzed its clinicopathological features, as well as the genotype using whole exome sequencing (WES). Moreover, we successfully established an organoid culture cell line from the primary tumor tissue. The patient was a 23-year-old woman who underwent superficial parotidectomy with peripheral neck dissection, followed by radiotherapy. Pathologically, the resected specimen showed atypical epithelioid nests and trabeculae with squamous differentiation, separated by thick fibrous septa, accompanied by dense lymphocytes and plasma cell infiltration. Immunohistochemistry revealed that the tumor cells were positive for AE1/AE3, p40, p63, p16, CK5/6, and CD5, and the background lymphocytes were positive for CD5 and CD99. Based on these findings, the tumor was diagnosed as CASTLE. WES uncovered five nonsynonymous and splicing somatic mutations, namely, FREM2 p.Val861Phe, CLK3 p.Phe376Leu, DLGAP1 p.Lys294Asn, NOX1 p.Val165Met, and PSG9 c.430 + 4A > T. Organoid culture cells preserved the histopathological characteristics of the epithelioid component of CASTLE and harbored all five somatic mutations detected in the primary tumor. In conclusion, for the first time to the best of our knowledge, we successfully analyzed a comprehensive genotype and established an organoid culture cell line of a parotid gland CASTLE, which should serve for analyzing the nature of this rare tumor.

Identification of Dominant Transcripts in Oxidative Stress Response by a Full-Length Transcriptome Analysis.

Our body responds to environmental stress by changing the expression levels of a series of cytoprotective enzymes/proteins through multilayered regulatory mechanisms, including the KEAP1-NRF2 system. While NRF2 upregulates the expression of many cytoprotective genes, there are fundamental limitations in short-read RNA sequencing (RNA-Seq), resulting in confusion regarding interpreting the effectiveness of cytoprotective gene induction at the transcript level. To precisely delineate isoform usage in the stress response, we conducted independent full-length transcriptome profiling (isoform sequencing; Iso-Seq) analyses of lymphoblastoid cells from three volunteers under normal and electrophilic stress-induced conditions. We first determined the first exon usage in KEAP1 and NFE2L2 (encoding NRF2) and found the presence of transcript diversity. We then examined changes in isoform usage of NRF2 target genes under stress conditions and identified a few isoforms dominantly expressed in the majority of NRF2 target genes. The expression levels of isoforms determined by Iso-Seq analyses showed striking differences from those determined by short-read RNA-Seq; the latter could be misleading concerning the abundance of transcripts. These results support that transcript usage is tightly regulated to produce functional proteins under electrophilic stress. Our present study strongly argues that there are important benefits that can be achieved by long-read transcriptome sequencing.

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.

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.

Matataki: an ultrafast mRNA quantification method for large-scale reanalysis of RNA-Seq data.

BACKGROUND: Data generated by RNA sequencing (RNA-Seq) is now accumulating in vast amounts in public repositories, especially for human and mouse genomes. Reanalyzing these data has emerged as a promising approach to identify gene modules or pathways. Although meta-analyses of gene expression data are frequently performed using microarray data, meta-analyses using RNA-Seq data are still rare. This lag is partly due to the limitations in reanalyzing RNA-Seq data, which requires extensive computational resources. Moreover, it is nearly impossible to calculate the gene expression levels of all samples in a public repository using currently available methods. Here, we propose a novel method, Matataki, for rapidly estimating gene expression levels from RNA-Seq data. RESULTS: The proposed method uses k-mers that are unique to each gene for the mapping of fragments to genes. Since aligning fragments to reference sequences requires high computational costs, our method could reduce the calculation cost by focusing on k-mers that are unique to each gene and by skipping uninformative regions. Indeed, Matataki outperformed conventional methods with regards to speed while demonstrating sufficient accuracy. CONCLUSIONS: The development of Matataki can overcome current limitations in reanalyzing RNA-Seq data toward improving the potential for discovering genes and pathways associated with disease at reduced computational cost. Thus, the main bottleneck of RNA-Seq analyses has shifted to achieving the decompression of sequenced data. The implementation of Matataki is available at https://github.com/informationsea/Matataki .

Genomic testing for pancreatic cancer in clinical practice as real-world evidence.

BACKGROUND: Precision medicine guided by comprehensive genome sequencing represents a potential treatment strategy for pancreatic cancer. However, clinical sequencing for pancreatic cancer entails several practical difficulties. We have launched an in-house clinical sequencing system and started genomic testing for patients with cancer in clinical practice. We have analyzed the clinical utility of this system in pancreatic cancer. METHODS: We retrospectively reviewed 20 patients with pancreatic cancer who visited our division. Genomic DNA was extracted from both tumor tissue and peripheral blood mononuclear cells obtained from the patients. We performed a comprehensive genomic testing using targeted amplicon sequencing for 160 cancer-related genes. The primary endpoints were the detection rates of potential actionable and druggable gene alterations. The secondary endpoints were the detection rate of secondary germline findings, the rate of re-biopsy required for genome sequencing, survival time after the initial visit (post-sequencing survival time), and turnaround time. RESULTS: Although re-biopsy was required for 25% (5/20) of all patients, genomic testing was performed in all patients. Actionable and druggable gene alterations were detected in 100% (20/20) and 35% (7/20) of patients, respectively, whereas secondary germline findings were detected in 5% (1/20) of patients. The median turnaround times for physicians and patients were 20 and 26 days, respectively. The median post-sequencing survival time was 10.3 months. Only 10% (2/20) of all patients were treated with therapeutic agents based on the outcomes of genomic testing. CONCLUSIONS: The clinical application of comprehensive genomic testing for pancreatic cancer was feasible and promising in clinical practice.

Zinc finger-IRF composite elements bound by Ikaros/IRF4 complexes function as gene repression in plasma cell.

The transcription factor (TF) interferon regulatory factor-4 (IRF4) promotes both germinal center (GC) reactions and plasma cell (PC) differentiation by binding to alternative DNA motifs including AP-1-IRF composite elements, Ets-IRF composite elements (EICEs), and interferon sequence response elements (ISREs). Although all of these motifs mediate transcriptional activation by IRF4, it is still unknown how some of the IRF4 target genes are downregulated upon PC differentiation. Here, we revealed a molecular mechanism of IRF4-mediated gene downregulation during PC differentiation. By combining IRF4 chromatin immunoprecipitation sequence and gene expression analysis, we identified zinc finger-IRF composite elements (ZICEs) in IRF4 binding regions aligned with genes whose expression was downregulated in PCs. The zinc finger TFs Ikaros and Aiolos were identified as IRF4 binding partners in PCs, and Ikaros but not Aiolos was essential for IRF4 binding to the ZICE sequence and for PC differentiation. The Ebf1 gene, which positively controls B-cell activation and GC reactions, was identified as one of the Ikaros/IRF4 target genes. Importantly, while the ZICE embeds the ISRE motif, IRF4 bound the ZICE motif as heterodimers with Ikaros for repression of target genes, which include Ebf1 In contrast, if the zinc finger motif is juxtaposed to the EICE motif, the Ikaros/PU.1/IRF4 complex functioned to activate target gene expression. Our findings revealed a novel mode of IRF4 activity upon PC differentiation where upon forming an Ikaros/IRF4 DNA-bound complex, a subset of genes is repressed.

Inflammatory responses induce an identity crisis of alveolar macrophages, leading to pulmonary alveolar proteinosis.

Pulmonary alveolar proteinosis (PAP) is a severe respiratory disease characterized by dyspnea caused by accumulation of surfactant protein. Dysfunction of alveolar macrophages (AMs), which regulate the homeostasis of surfactant protein, leads to the development of PAP; for example, in mice lacking BTB and CNC homology 2 (Bach2). However, how Bach2 helps prevent PAP is unknown, and the cell-specific effects of Bach2 are undefined. Using mice lacking Bach2 in specific cell types, we found that the PAP phenotype of Bach2-deficient mice is due to Bach2 deficiency in more than two types of immune cells. Depletion of hyperactivated T cells in Bach2-deficient mice restored normal function of AMs and ameliorated PAP. We also found that, in Bach2-deficient mice, hyperactivated T cells induced gene expression patterns that are specific to other tissue-resident macrophages and dendritic cells. Moreover, Bach2-deficient AMs exhibited a reduction in cell cycle progression. IFN-γ released from T cells induced Bach2 expression in AMs, in which Bach2 then bound to regulatory regions of inflammation-associated genes in myeloid cells. Of note, in AMs, Bach2 restricted aberrant responses to excessive T cell-induced inflammation, whereas, in T cells, Bach2 puts a brake on T cell activation. Moreover, Bach2 stimulated the expression of multiple histone genes in AMs, suggesting a role of Bach2 in proper histone expression. We conclude that Bach2 is critical for the maintenance of AM identity and self-renewal in inflammatory environments. Treatments targeting T cells may offer new therapeutic strategies for managing secondary PAP.

Establishment of Protocols for Global Metabolomics by LC-MS for Biomarker Discovery.

Metabolomics is a promising avenue for biomarker discovery. Although the quality of metabolomic analyses, especially global metabolomics (G-Met) using mass spectrometry (MS), largely depends on the instrumentation, potential bottlenecks still exist at several basic levels in the metabolomics workflow. Therefore, we established a precise protocol initially for the G-Met analyses of human blood plasma to overcome some these difficulties. In our protocol, samples are deproteinized in a 96-well plate using an automated liquid-handling system, and conducted either using a UHPLC-QTOF/MS system equipped with a reverse phase column or a LC-FTMS system equipped with a normal phase column. A normalization protocol of G-Met data was also developed to compensate for intra- and inter-batch differences, and the variations were significantly reduced along with our normalization, especially for the UHPLC-QTOF/MS data with a C18 reverse-phase column for positive ions. Secondly, we examined the changes in metabolomic profiles caused by the storage of EDTA-blood specimens to identify quality markers for the evaluation of the specimens' pre-analytical conditions. Forty quality markers, including lysophospholipids, dipeptides, fatty acids, succinic acid, amino acids, glucose, and uric acid were identified by G-Met for the evaluation of plasma sample quality and established the equation of calculating the quality score. We applied our quality markers to a small-scale study to evaluate the quality of clinical samples. The G-Met protocols and quality markers established here should prove useful for the discovery and development of biomarkers for a wider range of diseases.