Published papers

Comparative whole-genome sequencing enables the identification of specific mutations during adaptation of bacteria to new environments and allelic replacement can establish their causality. However, the mechanisms of action are hard to decipher and little has been achieved for epistatic mutations, especially at the metabolic level. Here we show that a strain of Escherichia coli carrying mutations in the rpoC and glpK genes, derived from adaptation in glycerol, uses two distinct metabolic strategies to gain growth advantage. A 27-bp deletion in the rpoC gene first increases metabolic efficiency. Then, a point mutation in the glpK gene promotes growth by improving glycerol utilization but results in increased carbon wasting as overflow metabolism. In a strain carrying both mutations, these contrasting carbon/energy saving and wasting mechanisms work together to give an 89% increase in growth rate. This study provides insight into metabolic reprogramming during adaptive laboratory evolution for fast cellular growth.

Background: Due to their bi-directional replication machinery starting from a single finite origin, bacterial genomes show characteristic nucleotide compositional bias between the two replichores, which can be visualised through GC skew or (C-G)/(C+G). Although this polarisation is used for computational prediction of replication origins in many bacterial genomes, the degree of GC skew visibility varies widely among different species, necessitating a quantitative measurement of GC skew strength in order to provide confidence measures for GC skew-based predictions of replication origins.Results: Here we discuss a quantitative index for the measurement of GC skew strength, named the generalised GC skew index (gGCSI), which is applicable to genomes of any length, including bacterial chromosomes and plasmids. We demonstrate that gGCSI is independent of the window size and can thus be used to compare genomes with different sizes, such as bacterial chromosomes and plasmids. It can suggest the existence of different replication mechanisms in archaea and of rolling-circle replication in plasmids. Correlation of gGCSI values between plasmids and their corresponding host chromosomes suggests that within the same strain, these replicons have reproduced using the same replication machinery and thus exhibit similar strengths of replication strand skew.Conclusions: gGCSI can be applied to genomes of any length and thus allows comparative study of replication-related mutation and selection pressures in genomes of different lengths such as bacterial chromosomes and plasmids. Using gGCSI, we showed that replication-related mutation or selection pressure is similar for replicons with similar machinery. © 2009 Arakawa et al; licensee BioMed Central Ltd.

Background: Biochemical pathways provide an essential context for understanding comprehensive experimental data and the systematic workings of a cell. Therefore, the availability of online pathway browsers will facilitate post-genomic research, just as genome browsers have contributed to genomics. Many pathway maps have been provided online as part of public pathway databases. Most of these maps, however, function as the gateway interface to a specific database, and the comprehensiveness of their represented entities, data mapping capabilities, and user interfaces are not always sufficient for generic usage. Methodology/Principal Findings: We have identified five central requirements for a pathway browser: (1) availability of large integrated maps showing genes, enzymes, and metabolites; (2) comprehensive search features and data access; (3) data mapping for transcriptomic, proteomic, and metabolomic experiments, as well as the ability to edit and annotate pathway maps; (4) easy exchange of pathway data; and (5) intuitive user experience without the requirement for installation and regular maintenance. According to these requirements, we have evaluated existing pathway databases and tools and implemented a web-based pathway browser named Pathway Projector as a solution. Conclusions/Significance: Pathway Projector provides integrated pathway maps that are based upon the KEGG Atlas, with the addition of nodes for genes and enzymes, and is implemented as a scalable, zoomable map utilizing the Google Maps API. Users can search pathway-related data using keywords, molecular weights, nucleotide sequences, and amino acid sequences, or as possible routes between compounds. In addition, experimental data from transcriptomic, proteomic, and metabolomic analyses can be readily mapped. Pathway Projector is freely available for academic users at http://www.g-language.org/PathwayProjector/. © 2009 Kono et al.

We developed a sample preparation protocol for rapid and unbiased analysis of the membrane proteome using an alimentary canal-mimicking system in which proteases are activated in the presence of bile salts. In this rapid and unbiased protocol, immobilized trypsin is used in the presence of deoxycholate and lauroylsarcosine to increase digestion efficiency as well as to increase the solubility of the membrane proteins. Using 22.5 μg of Escherichia coli whole cell lysate, we quantitatively demonstrated that membrane proteins were extracted and digested at the same level as soluble proteins without any solubility-related bias. The recovery of membrane proteins was independent of the number of transmembrane domains per protein. In the analysis of the membrane-enriched fraction from 22.5 μg of E. coli cell lysate, the abundance distribution of the membrane proteins was in agreement with that of the membrane protein-coding genes when this protocol, coupled with strong cation exchange prefractionation prior to nano-LC-MS/MS analysis, was used. Because this protocol allows unbiased sample preparation, protein abundance estimation based on the number of observed peptides per protein was applied to both soluble and membrane proteins simultaneously, and the copy numbers per cell for 1,453 E. coli proteins, including 545 membrane proteins, were successfully obtained. Finally, this protocol was applied to quantitative analysis of guanosine tetra- and pentaphosphate-dependent signaling in E. coli wild-type and relA knock-out strains. © 2009 by The American Society for Biochemistry and Molecular Biology, Inc.

Most bacterial genomes have biased nucleotide composition, and the asymmetry is considered to be caused by a single-stranded DNA (ssDNA) deamination arising from the bacterial replication machinery. In order to evaluate the relationship experimentally, the position and frequency of ssDNA formed during replication must be verified clearly. Although many ssDNA detection technologies exist, almost all methods have been developed for eukaryotic genomes. To apply these to bacterial genomes, which harbor a smaller amount of DNA than those of eukaryotes, more efficient, new methods are required. Therefore, we developed a novel strand-specific ssDNA sequencing method, called 4S-seq, for the bacterial genome. The 4S-seq method enriches ssDNA in the extracted genomic DNA by a dsDNA-specific nuclease and implements a strand-specific library using a biotin label with a customized tag. As a result, the 4S-seq is able to calculate the ssDNA content in each strand (Watson/Crick) at each position of the genome efficiently.

As the worldwide prevalence of colorectal cancer (CRC) increases, it is vital to reduce its morbidity and mortality through early detection. Saliva-based tests are an ideal noninvasive tool for CRC detection. Here, we explored and validated salivary biomarkers to distinguish patients with CRC from those with adenoma (AD) and healthy controls (HC). Saliva samples were collected from patients with CRC, AD, and HC. Untargeted salivary hydrophilic metabolite profiling was conducted using capillary electrophoresis–mass spectrometry and liquid chromatography–mass spectrometry. An alternative decision tree (ADTree)-based machine learning (ML) method was used to assess the discrimination abilities of the quantified metabolites. A total of 2602 unstimulated saliva samples were collected from subjects with CRC (n = 235), AD (n = 50), and HC (n = 2317). Data were randomly divided into training (n = 1301) and validation datasets (n = 1301). The clustering analysis showed a clear consistency of aberrant metabolites between the two groups. The ADTree model was optimized through cross-validation (CV) using the training dataset, and the developed model was validated using the validation dataset. The model discriminating CRC + AD from HC showed area under the receiver-operating characteristic curves (AUC) of 0.860 (95% confidence interval [CI]: 0.828-0.891) for CV and 0.870 (95% CI: 0.837-0.903) for the validation dataset. The other model discriminating CRC from AD + HC showed an AUC of 0.879 (95% CI: 0.851-0.907) and 0.870 (95% CI: 0.838-0.902), respectively. Salivary metabolomics combined with ML demonstrated high accuracy and versatility in detecting CRC.

Cancer-derived small extracellular vesicles (sEVs) are multifunctional particles with a lipid bilayer structure that are involved in cancer progression, such as malignant proliferation, distant metastasis, and cancer immunity evasion. The separation protocol used to isolate sEVs is an important process and thus, several have been developed, including ultracentrifugation (UC), size exclusion chromatography (SEC), and affinity purification using antibodies against sEV surface antigens. However, the effects of different separation methods on sEV components have not been adequately examined. Here, we developed a semi-automated system for collecting sEVs by combining SEC and preparative high-performance liquid chromatography and applied it to metabolome analysis. The developed SEC system could recover sEVs more efficiently and non-destructively than UC, suggesting that it is an appropriate recovery method for metabolic analysis and reflects biological conditions. Furthermore, using the developed SEC system, we performed metabolome analysis of sEVs from isocitrate dehydrogenase 1 (IDH)-mutated human colon HCT116 cells, which produce the oncogenic metabolite, 2-hydroxyglutaric acid (2-HG). IDH1-mutated HCT116 cells released significantly more sEVs than wild-type (WT) cells. The metabolomic profiles of IDH1 mutant and WT cells showed distinct differences between the cells and their sEVs. Notably, in IDH mutant cells, large amounts of 2-HG were detected not only in cells, but also in sEVs. These results indicate that the SEC system we developed has wide potential applications in sEVs research.

High-throughput metabolomics has enabled the development of large-scale cohort studies. Long-term studies require multiple batch-based measurements, which require sophisticated quality control (QC) to eliminate unexpected bias to obtain biologically meaningful quantified metabolomic profiles. Liquid chromatography–mass spectrometry was used to analyze 10,833 samples in 279 batch measurements. The quantified profile included 147 lipids including acylcarnitine, fatty acids, glucosylceramide, lactosylceramide, lysophosphatidic acid, and progesterone. Each batch included 40 samples, and 5 QC samples were measured for 10 samples of each. The quantified data from the QC samples were used to normalize the quantified profiles of the sample data. The intra- and inter-batch median coefficients of variation (CV) among the 147 lipids were 44.3% and 20.8%, respectively. After normalization, the CV values decreased by 42.0% and 14.7%, respectively. The effect of this normalization on the subsequent analyses was also evaluated. The demonstrated analyses will contribute to obtaining unbiased, quantified data for large-scale metabolomics.

The Clp1 family proteins, consisting of the Clp1 and Nol9/Grc3 groups, have polynucleotide kinase (PNK) activity at the 5′ end of RNA strands and are important enzymes in the processing of some precursor RNAs. However, it remains unclear how this enzyme family diversified in the eukaryotes. We performed a large-scale molecular evolutionary analysis of the full-length genomes of 358 eukaryotic species to classify the diverse Clp1 family proteins. The average number of Clp1 family proteins in eukaryotes was 2.3 ± 1.0, and most representative species had both Clp1 and Nol9/Grc3 proteins, suggesting that the Clp1 and Nol9/Grc3 groups were already formed in the eukaryotic ancestor by gene duplication. We also detected an average of 4.1 ± 0.4 Clp1 family proteins in members of the protist phylum Euglenozoa. For example, in Trypanosoma brucei, there are three genes of the Clp1 group and one gene of the Nol9/Grc3 group. In the Clp1 group proteins encoded by these three genes, the C-terminal domains have been replaced by unique characteristics domains, so we designated these proteins Tb-Clp1-t1, Tb-Clp1-t2, and Tb-Clp1-t3. Experimental validation showed that only Tb-Clp1-t2 has PNK activity against RNA strands. As in this example, N-terminal and C-terminal domain replacement also contributed to the diversification of the Clp1 family proteins in other eukaryotic species. Our analysis also revealed that the Clp1 family proteins in humans and plants diversified through isoforms created by alternative splicing.

Vascular endothelial growth factor-A (VEGF-A) is one of the primary factors promoting angiogenesis in endothelial cells. Although defects in VEGF-A signaling are linked to diverse pathophysiological conditions, the early phosphorylation-dependent signaling events pertinent to VEGF-A signaling remain poorly defined. Hence, a temporal quantitative phosphoproteomic analysis was performed in human umbilical vein endothelial cells (HUVECs) treated with VEGF-A-165 for 1, 5 and 10 min. This led to the identification and quantification of 1971 unique phosphopeptides corresponding to 961 phosphoproteins and 2771 phosphorylation sites in total. Specifically, 69, 153, and 133 phosphopeptides corresponding to 62, 125, and 110 phosphoproteins respectively, were temporally phosphorylated at 1, 5, and 10 min upon addition of VEGF-A. These phosphopeptides included 14 kinases, among others. This study also captured the phosphosignaling events directed through RAC, FAK, PI3K-AKT-MTOR, ERK, and P38 MAPK modules with reference to our previously assembled VEGF-A/VEGFR2 signaling pathway map in HUVECs. Apart from a significant enrichment of biological processes such as cytoskeleton organization and actin filament binding, our results also suggest a role of AAK1-AP2M1 in the regulation of VEGFR endocytosis. Taken together, the temporal quantitative phosphoproteomics analysis of VEGF signaling in HUVECs revealed early signaling events and we believe that this analysis will serve as a starting point for the analysis of differential signaling across VEGF members toward the full elucidation of their role in the angiogenesis processes. Graphical abstract: [Figure not available: see fulltext.]

Despite increasing global prevalence, the precise pathogenesis and terms for objective diagnosis of neurodegenerative dementias remain controversial, and comprehensive understanding of the disease remains lacking. Here, we conducted metabolomic analysis of serum and saliva obtained from patients with neurodegenerative dementias (n = 10), including Alzheimer's disease, frontotemporal lobe dementia, and Lewy body disease, as well as from age-matched healthy controls (n = 9). Using CE-TOF-MS, six metabolites in serum (β-alanine, creatinine, hydroxyproline, glutamine, iso-citrate, and cytidine) and two in saliva (arginine and tyrosine) were significantly different between dementias and controls. Using multivariate analysis, serum was confirmed as a more efficient biological fluid for diagnosis compared to saliva; additionally, 45 metabolites in total were identified as candidate markers that could discriminate at least one pair of diagnostic groups from the healthy control group. These metabolites possibly provide an objective method for diagnosing dementia-type by multiphase screening. Moreover, diagnostic-type-dependent differences were observed in several tricarboxylic acid cycle compounds detected in serum, indicating that some pathways in glucose metabolism may be altered in dementia patients. This pilot study revealed novel alterations in metabolomic profiles between various neurodegenerative dementias, which would contribute to etiological investigations. © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Background: Histone modifications play an important role in gene regulation. Acetylation of histone 3 lysine 9 (H3K9ac) is generally associated with transcription initiation and unfolded chromatin, thereby positively influencing gene expression. Deep sequencing of the 5' ends of gene transcripts using DeepCAGE delivers detailed information about the architecture and expression level of gene promoters. The combination of H3K9ac ChIP-chip and DeepCAGE in a myeloid leukemia cell line (THP-1) allowed us to study the spatial distribution of H3K9ac around promoters using a novel clustering approach. The promoter classes were analyzed for association with relevant genomic sequence features.Results: We performed a clustering of 4,481 promoters according to their surrounding H3K9ac signal and analyzed the clustered promoters for association with different sequence features. The clustering revealed three groups with major H3K9ac signal upstream, centered and downstream of the promoter. Narrow single peak promoters tend to have a concentrated activity of H3K9ac in the upstream region, while broad promoters tend to have a concentrated activity of H3K9ac and RNA polymerase II binding in the centered and downstream regions. A subset of promoters with high gene expression level, compared to subsets with low and medium gene expression, shows dramatic increase in H3K9ac activity in the upstream cluster only; this may indicate that promoters in the centered and downstream clusters are predominantly regulated at post-initiation steps. Furthermore, the upstream cluster is depleted in CpG islands and more likely to regulate un-annotated genes.Conclusions: Clustering core promoters according to their surrounding acetylation signal is a promising approach for the study of histone modifications. When examining promoters clustered into groups according to their surrounding H3K9 acetylation signal, we find that the relative localization and intensity of H3K9ac is very specific depending on characteristic sequence features of the promoter. Experimental data from DeepCAGE and ChIP-chip experiments using undifferentiated (monocyte) and differentiated (macrophage) THP-1 cells leads us to the same conclusions. © 2010 Kratz et al; licensee BioMed Central Ltd.

The structural dynamics of chromatin have been implicated in the regulation of fundamental eukaryotic processes, such as DNA transcription, replication and repair. Although previous studies have revealed that the chromatin landscape, nucleosome remodeling and histone modification events are intimately tied into cellular energetics and redox state, few studies undertake defined time-resolved measurements of these state variables. Here, we use metabolically synchronous, continuously-grown yeast cultures to measure DNA occupancy and track global patterns with respect to the metabolic state of the culture. Combined with transcriptome analyses and ChIP-qPCR experiments, these paint an intriguing picture where genome-wide nucleosome focusing occurs during the recovery of energy charge, followed by clearance of the promoter regions and global transcriptional slow-down, thus indicating a nucleosome-mediated “reset point” for the cycle. The reset begins at the end of the catabolic and stress-response transcriptional programs and ends prior to the start of the anabolic and cell-growth transcriptional program, and the histones on genes from both the catabolic and anabolic superclusters are deacetylated.

Acute kidney injury (AKI) is defined as a rapid decline in kidney function. The associated syndromes may lead to increased morbidity and mortality, but its early detection remains difficult. Using capillary electrophoresis time-of-flight mass spectrometry (CE-TOFMS), we analyzed the urinary metabolomic profile of patients admitted to the intensive care unit (ICU) after invasive surgery. Urine samples were collected at six time points: before surgery, at ICU admission and 6, 12, 24 and 48 h after. First, urine samples from 61 initial patients (non-AKI: 23, mild AKI: 24, severe AKI: 14) were measured, followed by the measurement of urine samples from 60 additional patients (non-AKI: 40, mild AKI: 20). Glycine and ethanolamine were decreased in patients with AKI compared with non-AKI patients at 6–24 h in the two groups. The linear statistical model constructed at each time point by machine learning achieved the best performance at 24 h (median AUC: 89%, cross-validated) for the 1st group. When cross-validated between the two groups, the AUC showed the best value of 70% at 12 h. These results identified metabolites and time points that show patterns specific to subjects who develop AKI, paving the way for the development of better biomarkers.

Malignant pleural mesothelioma (MPM) is an invasive malignancy that develops in the pleural cavity, and antifolates are used as chemotherapeutics for treating. The majority of antifolates, including pemetrexed (PMX), inhibit enzymes involved in purine and pyrimidine synthesis. MPM patients frequently develop drug resistance in clinical practice, however the associated drug-resistance mechanism is not well understood. This study was aimed to elucidate the mechanism underlying resistance to PMX in MPM cell lines. We found that among the differentially expressed genes associated with drug resistance (determined by RNA sequencing), TYMS expression was higher in the established resistant cell lines than in the parental cell lines. Knocking down TYMS expression significantly reduced drug resistance in the resistant cell lines. Conversely, TYMS overexpression significantly increased drug resistance in the parental cells. Metabolomics analysis revealed that the levels of dTMP were higher in the resistant cell lines than in the parental cell lines; however, resistant cells showed no changes in dTTP levels after PMX treatment. We found that the nucleic acid-biosynthetic pathway is important for predicting the efficacy of PMX in MPM cells. The results of chromatin immunoprecipitation-quantitative polymerase chain reaction (ChIP-qPCR) assays suggested that H3K27 acetylation in the 5′-UTR of TYMS may promote its expression in drug-resistant cells. Our findings indicate that the intracellular levels of dTMP are potential biomarkers for the effective treatment of patients with MPM and suggest the importance of regulatory mechanisms of TYMS expression in the disease.

Oscillations play a significant role in biological systems, with many examples in the fast, ultradian, circadian, circalunar, and yearly time domains. However, determining periodicity in such data can be problematic. There are a number of computational methods to identify the periodic components in large datasets, such as signal-to-noise based Fourier decomposition, Fisher's g-test and autocorrelation. However, the available methods assume a sinusoidal model and do not attempt to quantify the waveform shape and the presence of multiple periodicities, which provide vital clues in determining the underlying dynamics. Here, we developed a Fourier based measure that generates a de-noised waveform from multiple significant frequencies. This waveform is then correlated with the raw data from the respiratory oscillation found in yeast, to provide oscillation statistics including waveform metrics and multi-periods. The method is compared and contrasted to commonly used statistics. Moreover, we show the utility of the program in the analysis of noisy datasets and other high-throughput analyses, such as metabolomics and flow cytometry, respectively.

Artificially designed nucleases such as zinc-finger nucleases (ZFNs) and transcription activator-like effector nucleases (TALENs) can induce a targeted DNA double-strand break at the specific target genomic locus, leading to the frameshiftmediated gene disruption. However, the assays for their activity on the endogenous genomic loci remain limited. Herein, we describe a versatile modified lacZ assay to detect frameshifts in the nuclease target site. Short fragments of the genome DNA at the target or putative off-target loci were amplified from the genomic DNA of TALEN-treated or control embryos, and were inserted into the lacZa sequence for the conventional blue-white selection. The frequency of the frameshifts in the fragment can be estimated from the numbers of blue and white colonies. Insertions and/or deletions were easily determined by sequencing the plasmid DNAs recovered from the positive colonies. Our technique should offer broad application to the artificial nucleases for genome editing in various types of model organisms.

The cell cycle-dependent nucleocytoplasmic transport of proteins is predominantly regulated by CDK kinase activities; however, it is currently difficult to predict the proteins thus regulated, largely because of the low prediction efficiency of the motifs involved. Here, we report the successful prediction of CDK1-regulated nucleocytoplasmic shuttling proteins using a prediction system for nuclear localization signals (NLSs). By systematic amino acid replacement analyses in budding yeast, we created activity-based profiles for different classes of importin-α-dependent NLSs that represent the functional contributions of different amino acids at each position within an NLS class. We then developed a computer program for prediction of the classical importin-α/β pathway-specific NLSs (cNLS Mapper, available at http//nls-mapper.iab.keio.ac.jp/) that calculates NLS activities by using these profiles and an additivity-based motif scoring algorithm. This calculation method achieved significantly higher prediction accuracy in terms of both sensitivity and specificity than did current methods. The search for NLSs that overlap the consensus CDK1 phosphorylation site by using cNLS Mapper identified all previously reported and 5 previously uncharacterized yeast proteins (Yen1, Psy4, Pds1, Msa1, and Dna2) displaying CDK1- and cell cycle-regulated nuclear transport. CDK1 activated or repressed their nuclear import activity, depending on the position of CDK1-phosphorylation sites within NLSs. The application of this strategy to other functional linear motifs should be useful in systematic studies of protein-protein networks.

Background The hoverfly Microdon (Chymophila) katsurai Maruyama et Hironaga 2004 was speculated to be a myrmecophilous species associated with the ant Polyrhachis lamellidens based on observations of adults near the ant nest. However, there have been no reports regarding the observation of immature stages of this species in association with P. lamellidens. New information For the first time, we found three M. katsurai larvae inside a P. lamellidens nest and conducted rearing experiments on the larval M. katsurai. P. lamellidens workers did not show any inspection or attack behavior against the M. katsurai larvae under rearing conditions, suggesting that M. katsurai larvae can survive inside a P. lamellidens nest. Although no predatory behavior by the M. katsurai larvae was observed, all the M. katsurai larvae pupated and emerged in a rearing environment. The dorsal surface of the larval M. katsurai has a distinct pale green color with a uniform reticular structure. The puparium of M. katsurai shows several morphological features that are characteristic of the subgenus Chymophila. We conclude that M. katsurai is likely a myrmecophilous species that utilizes P. lamellidens as a specific host and that classification of M. katsurai based on puparium morphology is concordant with that based on adult morphology.