BackgroundMass spectrometry (MS) coupled with online separation methods is commonly applied for differential and quantitative profiling of biological samples in metabolomic as well as proteomic research. Such approaches are used for systems biology, functional genomics, and biomarker discovery, among others. An ongoing challenge of these molecular profiling approaches, however, is the development of better data processing methods. Here we introduce a new generation of a popular open-source data processing toolbox, MZmine 2.ResultsA key concept of the MZmine 2 software design is the strict separation of core functionality and data processing modules, with emphasis on easy usability and support for high-resolution spectra processing. Data processing modules take advantage of embedded visualization tools, allowing for immediate previews of parameter settings. Newly introduced functionality includes the identification of peaks using online databases, MSn data support, improved isotope pattern support, scatter plot visualization, and a new method for peak list alignment based on the random sample consensus (RANSAC) algorithm. The performance of the RANSAC alignment was evaluated using synthetic datasets as well as actual experimental data, and the results were compared to those obtained using other alignment algorithms.ConclusionsMZmine 2 is freely available under a GNU GPL license and can be obtained from the project website at: http://mzmine.sourceforge.net/. The current version of MZmine 2 is suitable for processing large batches of data and has been applied to both targeted and non-targeted metabolomic analyses.

BackgroundAdvances in DNA sequencing and proteomics have facilitated quantitative comparisons of snake venom composition. Most studies have employed one approach or the other. Here, both Illumina cDNA sequencing and LC/MS were used to compare the transcriptomes and proteomes of two pit vipers, Protobothrops flavoviridis and Ovophis okinavensis, which differ greatly in their biology.ResultsSequencing of venom gland cDNA produced 104,830 transcripts. The Protobothrops transcriptome contained transcripts for 103 venom-related proteins, while the Ovophis transcriptome contained 95. In both, transcript abundances spanned six orders of magnitude. Mass spectrometry identified peptides from 100% of transcripts that occurred at higher than contaminant (e.g. human keratin) levels, including a number of proteins never before sequenced from snakes. These transcriptomes reveal fundamentally different envenomation strategies. Adult Protobothrops venom promotes hemorrhage, hypotension, incoagulable blood, and prey digestion, consistent with mammalian predation. Ovophis venom composition is less readily interpreted, owing to insufficient pharmacological data for venom serine and metalloproteases, which comprise more than 97.3% of Ovophis transcripts, but only 38.0% of Protobothrops transcripts. Ovophis venom apparently represents a hybrid strategy optimized for frogs and small mammals.ConclusionsThis study illustrates the power of cDNA sequencing combined with MS profiling. The former quantifies transcript composition, allowing detection of novel proteins, but cannot indicate which proteins are actually secreted, as does MS. We show, for the first time, that transcript and peptide abundances are correlated. This means that MS can be used for quantitative, non-invasive venom profiling, which will be beneficial for studies of endangered species.

Metabolomics is a rapidly growing branch of post-genomic chemical biology. The fission yeast Schizosaccharomyces pombe is an excellent eukaryotic model organism. Although the entire S. pombe genome has been sequenced and detailed transcriptomic analyses were performed, little metabolic profiling has been done. Here we report the first global semi-quantitative analysis of the S. pombe metabolome using liquid chromatography high-resolution mass spectrometry. Procedures to obtain metabolic compounds from S. pombe extracts were established. One hundred and twenty-three distinct metabolites were identified while approximately 1900 peaks from the approximately 6000 observed were assigned. A software system (MZviewer) was developed to visualize semi-quantitative metabolome data using a dynamically generated scatter plot. We examined the metabolome of S. pombe cells exponentially grown in synthetic culture medium (EMM2) at two different temperatures, 26 degrees C and 36 degrees C. The profiles were similar except for varying amounts of certain amino acids and a significant increase in several compounds at 36 degrees C, such as trehalose (200-fold), glycerophosphoethanolamine (50-fold), arabitol (16-fold), ribulose (8-fold), and ophthalmic acid (5-fold). Reproducibility was demonstrated using a deletion mutant sib1Delta that lacked ferrichrome synthetase and showed no significant metabolic effects except the disappearance of the hexapeptide ferrichrome and the appearance of a putative dipeptide precursor. Taking advantage of the metabolic profile similarity at 26 degrees C and 36 degrees C, we analyzed the metabolome of a temperature-sensitive hcs1-143 mutant defective in the HMG-CoA synthase. As expected, HMG-CoA was decreased. In addition, extensive secondary metabolic effects, including a decrease in urea cycle intermediates and an increase in acetylated compounds, were observed. These findings confirm that S. pombe can be applied as an appropriate model to monitor metabolic responses to environmental conditions as well as genetic perturbations.

Molluscan shells, mainly composed of calcium carbonate, also contain organic components such as proteins and polysaccharides. Shell organic matrices construct frameworks of shell structures and regulate crystallization processes during shell formation. To date, a number of shell matrix proteins (SMPs) have been identified, and their functions in shell formation have been studied. However, previous studies focused only on SMPs extracted from adult shells, secreted after metamorphosis. Using proteomic analyses combined with genomic and transcriptomic analyses, we have identified 31 SMPs from larval shells of the pearl oyster, Pinctada fucata, and 111 from the Pacific oyster, Crassostrea gigas. Larval SMPs are almost entirely different from those of adults in both species. RNA-seq data also confirm that gene expression profiles for larval and adult shell formation are nearly completely different. Therefore, bivalves have two repertoires of SMP genes to construct larval and adult shells. Despite considerable differences in larval and adult SMPs, some functional domains are shared by both SMP repertoires. Conserved domains include von Willebrand factor type A (VWA), chitin-binding (CB), carbonic anhydrase (CA), and acidic domains. These conserved domains are thought to play crucial roles in shell formation. Furthermore, a comprehensive survey of animal genomes revealed that the CA and VWA–CB domain-containing protein families expanded in molluscs after their separation from other Lophotrochozoan linages such as the Brachiopoda. After gene expansion, some family members were co-opted for molluscan SMPs that may have triggered to develop mineralized shells from ancestral, nonmineralized chitinous exoskeletons.

Current proteomic studies clarified canonical synaptic proteins that are common to many types of synapses. However, proteins of diversified functions in a subset of synapses are largely hidden because of their low abundance or structural similarities to abundant proteins. To overcome this limitation, we have developed an “ultra-definition” (UD) subcellular proteomic workflow. Using purified synaptic vesicle (SV) fraction from rat brain, we identified 1,466 proteins, three times more than reported previously. This refined proteome includes all canonical SV proteins, as well as numerous proteins of low abundance, many of which were hitherto undetected. Comparison of UD quantifications between SV and synaptosomal fractions has enabled us to distinguish SV-resident proteins from potential SV-visitor proteins. We found 134 SV residents, of which 86 are present in an average copy number per SV of less than one, including vesicular transporters of nonubiquitous neurotransmitters in the brain. We provide a fully annotated resource of all categorized SV-resident and potential SV-visitor proteins, which can be utilized to drive novel functional studies, as we characterized here Aak1 as a regulator of synaptic transmission. Moreover, proteins in the SV fraction are associated with more than 200 distinct brain diseases. Remarkably, a majority of these proteins was found in the low-abundance proteome range, highlighting its pathological significance. Our deep SV proteome will provide a fundamental resource for a variety of future investigations on the function of synapses in health and disease.

Regulations of proliferation and quiescence in response to nutritional cues are important for medicine and basic biology. The fission yeast Schizosaccharomyces pombe serves as a model, owing to the shift of proliferating cells to the metabolically active quiescence (designate G0 phase hereafter) by responding to low nitrogen source. S. pombe G0 phase cells keep alive for months without growth and division. Nitrogen replenishment reinstates vegetative proliferation phase (designate VEG). Some 40 genes required for G0 maintenance were identified, but many more remain to be identified. We here show, using mutants, that the proteasome is required for maintaining G0 quiescence. Functional outcomes of proteasome in G0 and VEG phases appear to be distinct. Upon proteasome dysfunction, a number of antioxidant proteins and compounds responsive to ROS (reactive oxygen species) are produced. In addition, autophagy-mediated destruction of mitochondria occurs, which suppresses the loss of viability by eliminating ROS-generating mitochondria. These defensive responses are found in G0 but not in VEG, suggesting that the main function of proteasome in G0 phase homeostasis is to minimize ROS. Proteasome and autophagy are thus collaborative to support the lifespan of S. pombe G0 phase.cell cycle | cellular quiescence | proteolysis | mitochondria

BackgroundWhile many studies have shown that extracellular proteins evolve rapidly, how selection acts on them remains poorly understood. We used snake venoms to understand the interaction between ecology, expression level, and evolutionary rate in secreted protein systems. Venomous snakes employ well-integrated systems of proteins and organic constituents to immobilize prey. Venoms are generally optimized to subdue preferred prey more effectively than non-prey, and many venom protein families manifest positive selection and rapid gene family diversification. Although previous studies have illuminated how individual venom protein families evolve, how selection acts on venoms as integrated systems, is unknown.ResultsUsing next-generation transcriptome sequencing and mass spectrometry, we examined microevolution in two pitvipers, allopatrically separated for at least 1.6 million years, and their hybrids. Transcriptomes of parental species had generally similar compositions in regard to protein families, but for a given protein family, the homologs present and concentrations thereof sometimes differed dramatically. For instance, a phospholipase A2 transcript comprising 73.4 % of the Protobothrops elegans transcriptome, was barely present in the P. flavoviridis transcriptome (<0.05 %). Hybrids produced most proteins found in both parental venoms. Protein evolutionary rates were positively correlated with transcriptomic and proteomic abundances, and the most abundant proteins showed positive selection. This pattern holds with the addition of four other published crotaline transcriptomes, from two more genera, and also for the recently published king cobra genome, suggesting that rapid evolution of abundant proteins may be generally true for snake venoms. Looking more broadly at Protobothrops, we show that rapid evolution of the most abundant components is due to positive selection, suggesting an interplay between abundance and adaptation.ConclusionsGiven log-scale differences in toxin abundance, which are likely correlated with biosynthetic costs, we hypothesize that as a result of natural selection, snakes optimize return on energetic investment by producing more of venom proteins that increase their fitness. Natural selection then acts on the additive genetic variance of these components, in proportion to their contributions to overall fitness. Adaptive evolution of venoms may occur most rapidly through changes in expression levels that alter fitness contributions, and thus the strength of selection acting on specific secretome components.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-1832-6) contains supplementary material, which is available to authorized users.

Venom gland transcriptomes and proteomes of six Micrurus taxa (M. corallinus, M. lemniscatus carvalhoi, M. lemniscatus lemniscatus, M. paraensis, M. spixii spixii, and M. surinamensis) were investigated, providing the most comprehensive, quantitative data on Micrurus venom composition to date, and more than tripling the number of Micrurus venom protein sequences previously available. The six venomes differ dramatically. All are dominated by 2–6 toxin classes that account for 91–99% of the toxin transcripts. The M. s. spixii venome is compositionally the simplest. In it, three-finger toxins (3FTxs) and phospholipases A2 (PLA2s) comprise >99% of the toxin transcripts, which include only four additional toxin families at levels ≥0.1%. Micrurus l. lemniscatus venom is the most complex, with at least 17 toxin families. However, in each venome, multiple structural subclasses of 3FTXs and PLA2s are present. These almost certainly differ in pharmacology as well. All venoms also contain phospholipase B and vascular endothelial growth factors. Minor components (0.1–2.0%) are found in all venoms except that of M. s. spixii. Other toxin families are present in all six venoms at trace levels (<0.005%). Minor and trace venom components differ in each venom. Numerous novel toxin chemistries include 3FTxs with previously unknown 8- and 10-cysteine arrangements, resulting in new 3D structures and target specificities. 9-cysteine toxins raise the possibility of covalent, homodimeric 3FTxs or heterodimeric toxins with unknown pharmacologies. Probable muscarinic sequences may be reptile-specific homologs that promote hypotension via vascular mAChRs. The first complete sequences are presented for 3FTxs putatively responsible for liberating glutamate from rat brain synaptosomes. Micrurus C-type lectin-like proteins may have 6–9 cysteine residues and may be monomers, or homo- or heterodimers of unknown pharmacology. Novel KSPIs, 3× longer than any seen previously, appear to have arisen in three species by gene duplication and fusion. Four species have transcripts homologous to the nociceptive toxin, (MitTx) α-subunit, but all six species had homologs to the β-subunit. The first non-neurotoxic, non-catalytic elapid phospholipase A2s are reported. All are probably myonecrotic. Phylogenetic analysis indicates that the six taxa diverged 15–35 million years ago and that they split from their last common ancestor with Old World elapines nearly 55 million years ago. Given their early diversification, many cryptic micrurine taxa are anticipated.