An easy-to-use Decoy Database Builder software tool, implementing different decoy strategies for false discovery rate calculation in automated MS/MS protein identifications

Observations from the Proteomics Bench

Many challenges in proteomics result from the high-throughput nature of the experiments. This paper first presents pre-analytical problems, which still occur, although the call for standardization in omics has been ongoing for many years. This article also discusses aspects that affect bioinformatic analysis based on three sets of reference data measured with different orbitrap instruments. Despite continuous advances in mass spectrometer technology as well as analysis software, data-set-wise quality control is still necessary, and decoy-based estimation, although challenged by modern instruments, should be utilized. We draw attention to the fact that numerous young researchers perceive proteomics as a mature, readily applicable technology. However, it is important to emphasize that the maximum potential of the technology can only be realized by an educated handling of its limitations.

Target formation in muscle fibres indicates reinnervation - A proteomic study in muscle samples from peripheral neuropathies

AIMS

Target skeletal muscle fibres - defined by different concentric areas in oxidative enzyme staining - can occur in patients with neurogenic muscular atrophy. Here, we used our established hypothesis-free proteomic approach with the aim of deciphering the protein composition of targets. We also searched for potential novel interactions between target proteins.

METHODS

Targets and control areas were laser microdissected from skeletal muscle sections of 20 patients with neurogenic muscular atrophy. Samples were analysed by a highly sensitive mass spectrometry approach, enabling relative protein quantification. The results were validated by immunofluorescence studies. Protein interactions were investigated by yeast two-hybrid assays, coimmunoprecipitation experiments and bimolecular fluorescence complementation.

RESULTS

More than 1000 proteins were identified. Among these, 55 proteins were significantly over-represented and 40 proteins were significantly under-represented in targets compared to intraindividual control samples. The majority of over-represented proteins were associated with the myofibrillar Z-disc and actin dynamics, followed by myosin and myosin-associated proteins, proteins involved in protein biosynthesis and chaperones. Under-represented proteins were mainly mitochondrial proteins. Functional studies revealed that the LIM domain of the over-represented protein LIMCH1 interacts with isoform A of Xin actin-binding repeat-containing protein 1 (XinA).

CONCLUSIONS

In particular, proteins involved in myofibrillogenesis are over-represented in target structures, which indicate an ongoing process of sarcomere assembly and/or remodelling within this specific area of the muscle fibres. We speculate that target structures are the result of reinnervation processes in which filamin C-associated myofibrillogenesis is tightly regulated by the BAG3-associated protein quality system.

The first versatile human iPSC-based model of ectopic virus induction allows new insights in RNA-virus disease

A detailed description of pathophysiological effects that viruses exert on their host is still challenging. For the first time, we report a highly controllable viral expression model based on an iPS-cell line from a healthy human donor. The established viral model system enables a dose-dependent and highly localized RNA-virus expression in a fully controllable environment, giving rise for new applications for the scientific community.

Homozygous expression of the myofibrillar myopathy-associated p.W2710X filamin C variant reveals major pathomechanisms of sarcomeric lesion formation

Filamin C (FLNc) is mainly expressed in striated muscle cells where it localizes to Z-discs, myotendinous junctions and intercalated discs. Recent studies have revealed numerous mutations in the FLNC gene causing familial and sporadic myopathies and cardiomyopathies with marked clinical variability. The most frequent myopathic mutation, p.W2710X, which is associated with myofibrillar myopathy, deletes the carboxy-terminal 16 amino acids from FLNc and abolishes the dimerization property of Ig-like domain 24. We previously characterized "knock-in" mice heterozygous for this mutation (p.W2711X), and have now investigated homozygous mice using protein and mRNA expression analyses, mass spectrometry, and extensive immunolocalization and ultrastructural studies. Although the latter mice display a relatively mild myopathy under normal conditions, our analyses identified major mechanisms causing the pathophysiology of this disease: in comparison to wildtype animals (i) the expression level of FLNc protein is drastically reduced; (ii) mutant FLNc is relocalized from Z-discs to particularly mechanically strained parts of muscle cells, i.e. myotendinous junctions and myofibrillar lesions; (iii) the number of lesions is greatly increased and these lesions lack Bcl2-associated athanogene 3 (BAG3) protein; (iv) the expression of heat shock protein beta-7 (HSPB7) is almost completely abolished. These findings indicate grave disturbances of BAG3-dependent and -independent autophagy pathways that are required for efficient lesion repair. In addition, our studies reveal general mechanisms of lesion formation and demonstrate that defective FLNc dimerization via its carboxy-terminal domain does not disturb assembly and basic function of myofibrils. An alternative, more amino-terminally located dimerization site might compensate for that loss. Since filamins function as stress sensors, our data further substantiate that FLNc is important for mechanosensing in the context of Z-disc stabilization and maintenance.

A metastable subproteome underlies inclusion formation in muscle proteinopathies

Protein aggregation is a pathological feature of neurodegenerative disorders. We previously demonstrated that protein inclusions in the brain are composed of supersaturated proteins, which are abundant and aggregation-prone, and form a metastable subproteome. It is not yet clear, however, whether this phenomenon is also associated with non-neuronal protein conformational disorders. To respond to this question, we analyzed proteomic datasets from biopsies of patients with genetic and acquired protein aggregate myopathy (PAM) by quantifying the changes in composition, concentration and aggregation propensity of proteins in the fibers containing inclusions and those surrounding them. We found that a metastable subproteome is present in skeletal muscle from healthy patients. The expression of this subproteome escalate as proteomic samples are taken more proximal to the pathologic inclusion, eventually exceeding its solubility limits and aggregating. While most supersaturated proteins decrease or maintain steady abundance across healthy fibers and inclusion-containing fibers, proteins within the metastable subproteome rise in abundance, suggesting that they escape regulation. Taken together, our results show in the context of a human conformational disorder that the supersaturation of a metastable subproteome underlies widespread aggregation and correlates with the histopathological state of the tissue.

pValid: Validation Beyond the Target-Decoy Approach for Peptide Identification in Shotgun Proteomics

As the de facto validation method in mass spectrometry-based proteomics, the target-decoy approach determines a threshold to estimate the false discovery rate and then filters those identifications beyond the threshold. However, the incorrect identifications within the threshold are still unknown and further validation methods are needed. In this study, we characterized a framework of validation and investigated a number of common and novel validation methods. We first defined the accuracy of a validation method by its false-positive rate (FPR) and false-negative rate (FNR) and, further, proved that a validation method with lower FPR and FNR led to identifications with higher sensitivity and precision. Then we proposed a validation method named pValid that incorporated an open database search and a theoretical spectrum prediction strategy via a machine-learning technology. pValid was compared with four common validation methods as well as a synthetic peptide validation method. Tests on three benchmark data sets indicated that pValid had an FPR of 0.03% and an FNR of 1.79% on average, both superior to the other four common validation methods. Tests on a synthetic peptide data set also indicated that the FPR and FNR of pValid were better than those of the synthetic peptide validation method. Tests on a large-scale human proteome data set indicated that pValid successfully flagged the highest number of incorrect identifications among all five methods. Further considering its cost-effectiveness, pValid has the potential to be a feasible validation tool for peptide identification.

Effects of 12 Weeks of Hypertrophy Resistance Exercise Training Combined with Collagen Peptide Supplementation on the Skeletal Muscle Proteome in Recreationally Active Men

Evidence has shown that protein supplementation following resistance exercise training (RET) helps to further enhance muscle mass and strength. Studies have demonstrated that collagen peptides containing mostly non-essential amino acids increase fat-free mass (FFM) and strength in sarcopenic men. The aim of this study was to investigate whether collagen peptide supplementation in combination with RET influences the protein composition of skeletal muscle. Twenty-five young men (age: 24.2 ± 2.6 years, body mass (BM): 79.6 ± 5.6 kg, height: 185.0 ± 5.0 cm, fat mass (FM): 11.5% ± 3.4%) completed body composition and strength measurements and vastus lateralis biopsies were taken before and after a 12-week training intervention. In a double-blind, randomized design, subjects consumed either 15 g of specific collagen peptides (COL) or a non-caloric placebo (PLA) every day within 60 min after their training session. A full-body hypertrophy workout was completed three times per week and included four exercises using barbells. Muscle proteome analysis was performed by liquid chromatography tandem mass spectrometry (LC-MS/MS). BM and FFM increased significantly in COL compared with PLA, whereas no differences in FM were detected between the two groups. Both groups improved in strength levels, with a slightly higher increase in COL compared with PLA. In COL, 221 higher abundant proteins were identified. In contrast, only 44 proteins were of higher abundance in PLA. In contrast to PLA, the upregulated proteins in COL were mostly associated with the protein metabolism of the contractile fibers. In conclusion, the use of RET in combination with collagen peptide supplementation results in a more pronounced increase in BM, FFM, and muscle strength than RET alone. More proteins were upregulated in the COL intervention most of which were associated with contractile fibers.

Defining the RNA interactome by total RNA-associated protein purification

The RNA binding proteome (RBPome) was previously investigated using UV crosslinking and purification of poly(A)-associated proteins. However, most cellular transcripts are not polyadenylated. We therefore developed total RNA-associated protein purification (TRAPP) based on 254 nm UV crosslinking and purification of all RNA-protein complexes using silica beads. In a variant approach (PAR-TRAPP), RNAs were labelled with 4-thiouracil prior to 350 nm crosslinking. PAR-TRAPP in yeast identified hundreds of RNA binding proteins, strongly enriched for canonical RBPs. In comparison, TRAPP identified many more proteins not expected to bind RNA, and this correlated strongly with protein abundance. Comparing TRAPP in yeast and E. coli showed apparent conservation of RNA binding by metabolic enzymes. Illustrating the value of total RBP purification, we discovered that the glycolytic enzyme enolase interacts with tRNAs. Exploiting PAR-TRAPP to determine the effects of brief exposure to weak acid stress revealed specific changes in late 60S ribosome biogenesis. Furthermore, we identified the precise sites of crosslinking for hundreds of RNA-peptide conjugates, using iTRAPP, providing insights into potential regulation. We conclude that TRAPP is a widely applicable tool for RBPome characterization.

A new estimation of protein-level false discovery rate

Background

In mass spectrometry-based proteomics, protein identification is an essential task. Evaluating the statistical significance of the protein identification result is critical to the success of proteomics studies. Controlling the false discovery rate (FDR) is the most common method for assuring the overall quality of the set of identifications. Existing FDR estimation methods either rely on specific assumptions or rely on the two-stage calculation process of first estimating the error rates at the peptide-level, and then combining them somehow at the protein-level. We propose to estimate the FDR in a non-parametric way with less assumptions and to avoid the two-stage calculation process.

Results

We propose a new protein-level FDR estimation framework. The framework contains two major components: the Permutation+BH (Benjamini–Hochberg) FDR estimation method and the logistic regression-based null inference method. In Permutation+BH, the null distribution of a sample is generated by searching data against a large number of permuted random protein database and therefore does not rely on specific assumptions. Then, p-values of proteins are calculated from the null distribution and the BH procedure is applied to the p-values to achieve the relationship of the FDR and the number of protein identifications. The Permutation+BH method generates the null distribution by the permutation method, which is inefficient for online identification. The logistic regression model is proposed to infer the null distribution of a new sample based on existing null distributions obtained from the Permutation+BH method.

Conclusions

In our experiment based on three public available datasets, our Permutation+BH method achieves consistently better performance than MAYU, which is chosen as the benchmark FDR calculation method for this study. The null distribution inference result shows that the logistic regression model achieves a reasonable result both in the shape of the null distribution and the corresponding FDR estimation result.

An sRNA and Cold Shock Protein Homolog-Based Feedforward Loop Post-transcriptionally Controls Cell Cycle Master Regulator CtrA

Adjustment of cell cycle progression is crucial for bacterial survival and adaptation under adverse conditions. However, the understanding of modulation of cell cycle control in response to environmental changes is rather incomplete. In α-proteobacteria, the broadly conserved cell cycle master regulator CtrA underlies multiple levels of control, including coupling of cell cycle and cell differentiation. CtrA levels are known to be tightly controlled through diverse transcriptional and post-translational mechanisms. Here, small RNA (sRNA)-mediated post-transcriptional regulation is uncovered as an additional level of CtrA fine-tuning. Computational predictions as well as transcriptome and proteome studies consistently suggested targeting of ctrA and the putative cold shock chaperone cspA5 mRNAs by the trans-encoded sRNA (trans-sRNA) GspR (formerly SmelC775) in several Sinorhizobium species. GspR strongly accumulated in the stationary growth phase, especially in minimal medium (MM) cultures. Lack of the gspR locus confers a fitness disadvantage in competition with the wild type, while its overproduction hampers cell growth, suggesting that this riboregulator interferes with cell cycle progression. An eGFP-based reporter in vivo assay, involving wild-type and mutant sRNA and mRNA pairs, experimentally confirmed GspR-dependent post-transcriptional down-regulation of ctrA and cspA5 expression, which most likely occurs through base-pairing to the respective mRNA. The energetically favored secondary structure of GspR is predicted to comprise three stem-loop domains, with stem-loop 1 and stem-loop 3 targeting ctrA and cspA5 mRNA, respectively. Moreover, this work reports evidence for post-transcriptional control of ctrA by CspA5. Thus, this regulation and GspR-mediated post-transcriptional repression of ctrA and cspA5 expression constitute a coherent feed-forward loop, which may enhance the negative effect of GspR on CtrA levels. This novel regulatory circuit involving the riboregulator GspR, CtrA, and a cold shock chaperone may contribute to fine-tuning of ctrA expression.

Loquacious-PD facilitates Drosophila Dicer-2 cleavage through interactions with the helicase domain and dsRNA

Loquacious-PD (Loqs-PD) is required for biogenesis of many endogenous siRNAs in Drosophila In vitro, Loqs-PD enhances the rate of dsRNA cleavage by Dicer-2 and also enables processing of substrates normally refractory to cleavage. Using purified components, and Loqs-PD truncations, we provide a mechanistic basis for Loqs-PD functions. Our studies indicate that the 22 amino acids at the C terminus of Loqs-PD, including an FDF-like motif, directly interact with the Hel2 subdomain of Dicer-2's helicase domain. This interaction is RNA-independent, but we find that modulation of Dicer-2 cleavage also requires dsRNA binding by Loqs-PD. Furthermore, while the first dsRNA-binding motif of Loqs-PD is dispensable for enhancing cleavage of optimal substrates, it is essential for enhancing cleavage of suboptimal substrates. Finally, our studies define a previously unrecognized Dicer interaction interface and suggest that Loqs-PD is well positioned to recruit substrates into the helicase domain of Dicer-2.

Improving a Synechocystis-based photoautotrophic chassis through systematic genome mapping and validation of neutral sites

The use of microorganisms as cell factories frequently requires extensive molecular manipulation. Therefore, the identification of genomic neutral sites for the stable integration of ectopic DNA is required to ensure a successful outcome. Here we describe the genome mapping and validation of five neutral sites in the chromosome of Synechocystis sp. PCC 6803, foreseeing the use of this cyanobacterium as a photoautotrophic chassis. To evaluate the neutrality of these loci, insertion/deletion mutants were produced, and to assess their functionality, a synthetic green fluorescent reporter module was introduced. The constructed integrative vectors include a BioBrick-compatible multiple cloning site insulated by transcription terminators, constituting robust cloning interfaces for synthetic biology approaches. Moreover, Synechocystis mutants (chassis) ready to receive purpose-built synthetic modules/circuits are also available. This work presents a systematic approach to map and validate chromosomal neutral sites in cyanobacteria, and that can be extended to other organisms.

New glycoproteomics software, GlycoPep Evaluator, generates decoy glycopeptides de novo and enables accurate false discovery rate analysis for small data sets

Glycoproteins are biologically significant large molecules that participate in numerous cellular activities. In order to obtain site-specific protein glycosylation information, intact glycopeptides, with the glycan attached to the peptide sequence, are characterized by tandem mass spectrometry (MS/MS) methods such as collision-induced dissociation (CID) and electron transfer dissociation (ETD). While several emerging automated tools are developed, no consensus is present in the field about the best way to determine the reliability of the tools and/or provide the false discovery rate (FDR). A common approach to calculate FDRs for glycopeptide analysis, adopted from the target-decoy strategy in proteomics, employs a decoy database that is created based on the target protein sequence database. Nonetheless, this approach is not optimal in measuring the confidence of N-linked glycopeptide matches, because the glycopeptide data set is considerably smaller compared to that of peptides, and the requirement of a consensus sequence for N-glycosylation further limits the number of possible decoy glycopeptides tested in a database search. To address the need to accurately determine FDRs for automated glycopeptide assignments, we developed GlycoPep Evaluator (GPE), a tool that helps to measure FDRs in identifying glycopeptides without using a decoy database. GPE generates decoy glycopeptides de novo for every target glycopeptide, in a 1:20 target-to-decoy ratio. The decoys, along with target glycopeptides, are scored against the ETD data, from which FDRs can be calculated accurately based on the number of decoy matches and the ratio of the number of targets to decoys, for small data sets. GPE is freely accessible for download and can work with any search engine that interprets ETD data of N-linked glycopeptides. The software is provided at https://desairegroup.ku.edu/research.

Effects of zinc on particulate methane monooxygenase activity and structure

Particulate methane monooxygenase (pMMO) is a membrane-bound metalloenzyme that oxidizes methane to methanol in methanotrophic bacteria. Zinc is a known inhibitor of pMMO, but the details of zinc binding and the mechanism of inhibition are not understood. Metal binding and activity assays on membrane-bound pMMO from Methylococcus capsulatus (Bath) reveal that zinc inhibits pMMO at two sites that are distinct from the copper active site. The 2.6 Å resolution crystal structure of Methylocystis species strain Rockwell pMMO reveals two previously undetected bound lipids, and metal soaking experiments identify likely locations for the two zinc inhibition sites. The first is the crystallographic zinc site in the pmoC subunit, and zinc binding here leads to the ordering of 10 previously unobserved residues. A second zinc site is present on the cytoplasmic side of the pmoC subunit. Parallels between these results and zinc inhibition studies of several respiratory complexes suggest that zinc might inhibit proton transfer in pMMO.

Comparison of alternative MS/MS and bioinformatics approaches for confident phosphorylation site localization

Over the past years, phosphoproteomics has advanced to a prime tool in signaling research. Since then, an enormous amount of information about in vivo protein phosphorylation events has been collected providing a treasure trove for gaining a better understanding of the molecular processes involved in cell signaling. Yet, we still face the problem of how to achieve correct modification site localization. Here we use alternative fragmentation and different bioinformatics approaches for the identification and confident localization of phosphorylation sites. Phosphopeptide-enriched fractions were analyzed by multistage activation, collision-induced dissociation and electron transfer dissociation (ETD), yielding complementary phosphopeptide identifications. We further found that MASCOT, OMSSA and Andromeda each identified a distinct set of phosphopeptides allowing the number of site assignments to be increased. The postsearch engine SLoMo provided confident phosphorylation site localization, whereas different versions of PTM-Score integrated in MaxQuant differed in performance. Based on high-resolution ETD and higher collisional dissociation (HCD) data sets from a large synthetic peptide and phosphopeptide reference library reported by Marx et al. [Nat. Biotechnol. 2013, 31 (6), 557-564], we show that an Andromeda/PTM-Score probability of 1 is required to provide an false localization rate (FLR) of 1% for HCD data, while 0.55 is sufficient for high-resolution ETD spectra. Additional analyses of HCD data demonstrated that for phosphotyrosine peptides and phosphopeptides containing two potential phosphorylation sites, PTM-Score probability cutoff values of <1 can be applied to ensure an FLR of 1%. Proper adjustment of localization probability cutoffs allowed us to significantly increase the number of confident sites with an FLR of <1%.Our findings underscore the need for the systematic assessment of FLRs for different score values to report confident modification site localization.

High-resolution mass spectrometry driven discovery of peptidic danger signals in insect immunity

The 'danger model' is an alternative concept for immune response postulating that the immune system reacts to entities that do damage (danger associated molecular patterns, DAMP) and not only to entities that are foreign (pathogen-associated molecular patterns, PAMP) as proposed by classical immunology concepts. In this study we used Galleria mellonella to validate the danger model in insects. Hemolymph of G. mellonella was digested with thermolysin (as a representative for virulence-associated metalloproteinases produced by humanpathogens) followed by chromatographic fractionation. Immune-stimulatory activity was tested by measuring lysozyme activity with the lytic zone assays against Micrococcus luteus cell wall components. Peptides were analyzed by nano-scale liquid chromatography coupled to high-resolution Fourier transform mass spectrometers. Addressing the lack of a genome sequence we complemented the rudimentary NCBI protein database with a recently established transcriptome and de novo sequencing methods for peptide identification. This approach led to identification of 127 peptides, 9 of which were identified in bioactive fractions. Detailed MS/MS experiments in comparison with synthetic analogues confirmed the amino acid sequence of all 9 peptides. To test the potential of these putative danger signals to induce immune responses we injected the synthetic analogues into G. mellonella and monitored the anti-bacterial activity against living Micrococcus luteus. Six out of 9 peptides identified in the bioactive fractions exhibited immune-stimulatory activity when injected. Hence, we provide evidence that small peptides resulting from thermolysin-mediated digestion of hemolymph proteins function as endogenous danger signals which can set the immune system into alarm. Consequently, our study indicates that the danger model also plays a role in insect immunity.

Proteomic-based insight into Malpighian tubules of silkworm Bombyx mori

Malpighian tubules (MTs) are highly specific organs of arthropods (Insecta, Myriapoda and Arachnida) for excretion and osmoregulation. In order to highlight the important genes and pathways involved in multi-functions of MTs, we performed a systematic proteomic analysis of silkworm MTs in the present work. Totally, 1,367 proteins were identified by one-dimensional gel electrophoresis coupled with liquid chromatography-tandem mass spectrometry, and as well as by Trans Proteomic Pipeline (TPP) and Absolute protein expression (APEX) analyses. Forty-one proteins were further identified by two-dimensional gel electrophoresis. Some proteins were revealed to be significantly associated with various metabolic processes, organic solute transport, detoxification and innate immunity. Our results might lay a good foundation for future functional studies of MTs in silkworm and other lepidoptera.

Restoration of cytoskeleton homeostasis after gigaxonin gene transfer for giant axonal neuropathy

Giant axonal neuropathy (GAN) is caused by loss of function of the gigaxonin protein. On a cellular level GAN is characterized by intermediate filament (IF) aggregation, leading to a progressive and fatal peripheral neuropathy in humans. This study sought to determine if re-introduction of the GAN gene into GAN-deficient cells and mice would restore proper cytoskeleton IF homeostasis. Treatment of primary skin fibroblast cultures from three different GAN patients with an adeno-associated virus type 2 (AAV2) vector containing a normal human GAN transgene significantly reduced the number of cells displaying vimentin IF aggregates. A proteomic analysis of these treated cells was also performed, wherein the abundance of 32 of 780 identified proteins significantly changed in response to gigaxonin gene transfer. While 29 of these responding proteins have not been directly described in association with gigaxonin, three were previously identified as being disregulated in GAN and were now shifted toward normal levels. To assess the potential application of this approach in vivo and eventually in humans, GAN mice received an intracisternal injection of an AAV9/GAN vector to globally deliver the GAN gene to the brainstem and spinal cord. The treated mice showed a nearly complete clearance of peripherin IF accumulations at 3 weeks post-injection. These studies demonstrate that gigaxonin gene transfer can reverse the cellular IF aggregate pathology associated with GAN.

Shotgun analysis on the peritrophic membrane of the silkworm Bombyx mori

The insect midgut epithelium is generally lined with a unique chitin and protein structure, the peritrophic membrane (PM), which facilitates food digestion and protects the gut epithelium. We used gel electrophoresis and mass spectrometry to identify the extracted proteins from the silkworm PM to obtain an in-depth understanding of the biological function of the silkworm PM components. A total of 305 proteins, with molecular weights ranging from 8.02 kDa to 788.52 kDa and the isoelectric points ranging from 3.39 to 12.91, were successfully identified. We also found several major classes of PM proteins, i.e. PM chitin-binding protein, invertebrate intestinal mucin, and chitin deacetylase. The protein profile provides a basis for further study of the physiological events in the PM of Bombyx mori. [BMB Reports 2012; 45(11): 665-670]

Differential proteomic analysis of abnormal intramyoplasmic aggregates in desminopathy

Desminopathy is a subtype of myofibrillar myopathy caused by desmin mutations and characterized by protein aggregates accumulating in muscle fibers. The aim of this study was to assess the protein composition of these aggregates. Aggregates and intact myofiber sections were obtained from skeletal muscle biopsies of five desminopathy patients by laser microdissection and analyzed by a label-free spectral count-based proteomic approach. We identified 397 proteins with 22 showing significantly higher spectral indices in aggregates (ratio >1.8, p<0.05). Fifteen of these proteins not previously reported as specific aggregate components provide new insights regarding pathomechanisms of desminopathy. Results of proteomic analysis were supported by immunolocalization studies and parallel reaction monitoring. Three mutant desmin variants were detected directly on the protein level as components of the aggregates, suggesting their direct involvement in aggregate-formation and demonstrating for the first time that proteomic analysis can be used for direct identification of a disease-causing mutation in myofibrillar myopathy. Comparison of the proteomic results in desminopathy with our previous analysis of aggregate composition in filaminopathy, another myofibrillar myopathy subtype, allows to determine subtype-specific proteomic profile that facilitates identification of the specific disorder.

BIOLOGICAL SIGNIFICANCE

Our proteomic analysis provides essential new insights in the composition of pathological protein aggregates in skeletal muscle fibers of desminopathy patients. The results contribute to a better understanding of pathomechanisms in myofibrillar myopathies and provide the basis for hypothesis-driven studies. The detection of specific proteomic profiles in different myofibrillar myopathy subtypes indicates that proteomic analysis may become a useful tool in differential diagnosis of protein aggregate myopathies.

Analysis of the membrane proteome of ciprofloxacin-resistant macrophages by stable isotope labeling with amino acids in cell culture (SILAC)

Overexpression of multidrug transporters is a well-established mechanism of resistance to chemotherapy, but other changes may be co-selected upon exposure to drugs that contribute to resistance. Using a model of J774 macrophages made resistant to the fluoroquinolone antibiotic ciprofloxacin and comparing it with the wild-type parent cell line, we performed a quantitative proteomic analysis using the stable isotope labeling with amino acids in cell culture technology coupled with liquid chromatography electrospray ionization Fourier transform tandem mass spectrometry (LC-ESI-FT-MS/MS) on 2 samples enriched in membrane proteins (fractions F1 and F2 collected from discontinuous sucrose gradient). Nine hundred proteins were identified with at least 3 unique peptides in these 2 pooled fractions among which 61 (F1) and 69 (F2) showed a significantly modified abundance among the 2 cell lines. The multidrug resistance associated protein Abcc4, known as the ciprofloxacin efflux transporter in these cells, was the most upregulated, together with Dnajc3, a protein encoded by a gene located downstream of Abcc4. The other modulated proteins are involved in transport functions, cell adhesion and cytoskeleton organization, immune response, signal transduction, and metabolism. This indicates that the antibiotic ciprofloxacin is able to trigger a pleiotropic adaptative response in macrophages that includes the overexpression of its efflux transporter.

The proteome of human liver peroxisomes: identification of five new peroxisomal constituents by a label-free quantitative proteomics survey

The peroxisome is a key organelle of low abundance that fulfils various functions essential for human cell metabolism. Severe genetic diseases in humans are caused by defects in peroxisome biogenesis or deficiencies in the function of single peroxisomal proteins. To improve our knowledge of this important cellular structure, we studied for the first time human liver peroxisomes by quantitative proteomics. Peroxisomes were isolated by differential and Nycodenz density gradient centrifugation. A label-free quantitative study of 314 proteins across the density gradient was accomplished using high resolution mass spectrometry. By pairing statistical data evaluation, cDNA cloning and in vivo colocalization studies, we report the association of five new proteins with human liver peroxisomes. Among these, isochorismatase domain containing 1 protein points to the existence of a new metabolic pathway and hydroxysteroid dehydrogenase like 2 protein is likely involved in the transport or β-oxidation of fatty acids in human peroxisomes. The detection of alcohol dehydrogenase 1A suggests the presence of an alternative alcohol-oxidizing system in hepatic peroxisomes. In addition, lactate dehydrogenase A and malate dehydrogenase 1 partially associate with human liver peroxisomes and enzyme activity profiles support the idea that NAD⁺ becomes regenerated during fatty acid β-oxidation by alternative shuttling processes in human peroxisomes involving lactate dehydrogenase and/or malate dehydrogenase. Taken together, our data represent a valuable resource for future studies of peroxisome biochemistry that will advance research of human peroxisomes in health and disease.

Crowdsourcing in proteomics: public resources lead to better experiments

With the growing interest in the field of proteomics, the amount of publicly available proteome resources has also increased dramatically. This means that there are many useful resources available for almost all aspects of a proteomics experiment. However, it remains vital to use the right resource, for the right purpose, at the right time. This review is therefore meant to aid the reader in obtaining an overview of the available resources and their application, thus providing the necessary background to choose the appropriate resources for the experiment at hand. Many of the resources are also taking advantage of so-called crowdsourcing to maximize the potential of the resource. What this means and how this can improve future experiments will also be discussed. The text roughly follows the steps involved in a proteomics experiment, starting with the planning of the experiment, via the processing of the data and the analysis of the results, to the community-wide sharing of the produced data.

FE65 regulates and interacts with the Bloom syndrome protein in dynamic nuclear spheres – potential relevance to Alzheimer's disease

The intracellular domain of the amyloid precursor protein (AICD) is generated following cleavage of the precursor by the γ-secretase complex and is involved in membrane to nucleus signaling, for which the binding of AICD to the adapter protein FE65 is essential. Here we show that FE65 knockdown causes a down regulation of the protein BLM and the MCM protein family and that elevated nuclear levels of FE65 result in stabilization of the BLM protein in nuclear mobile spheres. These spheres are able to grow and fuse, and potentially correspond to the nuclear domain 10. BLM plays a role in DNA replication and repair mechanisms and FE65 was also shown to play a role in the cell's response to DNA damage. A set of proliferation assays in our work revealed that FE65 knockdown cells exhibit reduced cell replication in HEK293T cells. On the basis of these results, we hypothesize that nuclear FE65 levels (nuclear FE65/BLM containing spheres) may regulate cell cycle re-entry in neurons due to increased interaction of FE65 with BLM and/or an increase in MCM protein levels. Thus, FE65 interactions with BLM and MCM proteins may contribute to the neuronal cell cycle re-entry observed in Alzheimer disease brains.

Proteomic analysis of alterations induced by perinatal hypoxic-ischemic brain injury

Perinatal hypoxic-ischemic brain injury is an important cause of neurological deficits still causing mortality and morbidity in the early period of life. As efficient clinical or pharmaceutical strategies to prevent or reduce the outcome of perinatal hypoxic-ischemic brain damage are limited, the development of new therapies is of utmost importance. To evolve innovative therapeutic concepts, elucidation of the mechanisms contributing to the neurological impairments upon hypoxic-ischemic brain injury is necessary. Therefore, we aimed for the identification of proteins that are affected by hypoxic-ischemic brain injury in neonatal rats. To assess changes in protein expression two days after induction of brain damage, a 2D-DIGE based proteome analysis was performed. Among the proteins altered after hypoxic-ischemic brain injury, Calcineurin A, Coronin-1A, as well as GFAP were identified, showing higher expression in lesioned hemispheres. Validation of the changes in Calcineurin A expression by Western Blot analysis demonstrated several truncated forms of this protein generated by limited proteolysis after hypoxia-ischemia. Further analysis revealed activation of calpain, which is involved in the limited proteolysis of Calcineurin. Active forms of Calcineurin are associated with the dephosphorylation of Darpp-32, an effect that was also demonstrated in lesioned hemispheres after perinatal brain injury.

A combined laser microdissection and mass spectrometry approach reveals new disease relevant proteins accumulating in aggregates of filaminopathy patients

Filaminopathy is a subtype of myofibrillar myopathy caused by mutations in FLNC, the gene encoding filamin C, and histologically characterized by pathologic accumulation of several proteins within skeletal muscle fibers. With the aim to get new insights in aggregate composition, we collected aggregates and control tissue from skeletal muscle biopsies of six myofibrillar myopathy patients harboring three different FLNC mutations by laser microdissection and analyzed the samples by a label-free mass spectrometry approach. A total of 390 proteins were identified, and 31 of those showed significantly higher spectral indices in aggregates compared with patient controls with a ratio >1.8. These proteins included filamin C, other known myofibrillar myopathy associated proteins, and a striking number of filamin C binding partners. Across the patients the patterns were extremely homogeneous. Xin actin-binding repeat containing protein 2, heat shock protein 27, nebulin-related-anchoring protein, and Rab35 could be verified as new filaminopathy biomarker candidates. In addition, further experiments identified heat shock protein 27 and Xin actin-binding repeat containing protein 2 as novel filamin C interaction partners and we could show that Xin actin-binding repeat containing protein 2 and the known interaction partner Xin actin-binding repeat containing protein 1 simultaneously associate with filamin C. Ten proteins showed significant lower spectral indices in aggregate samples compared with patient controls (ratio <0.56) including M-band proteins myomesin-1 and myomesin-2. Proteomic findings were consistent with previous and novel immunolocalization data. Our findings suggest that aggregates in filaminopathy have a largely organized structure of proteins also interacting under physiological conditions. Different filamin C mutations seem to lead to almost identical aggregate compositions. The finding that filamin C was detected as highly abundant protein in aggregates in filaminopathy indicates that our proteomic approach may be suitable to identify new candidate genes among the many MFM patients with so far unknown mutation.

Proteomic analysis in giant axonal neuropathy: new insights into disease mechanisms

INTRODUCTION

Giant axonal neuropathy (GAN) is a progressive hereditary disease that affects the peripheral and central nervous systems. It is characterized morphologically by aggregates of intermediate filaments in different tissues. Mutations have been reported in the gene that codes for gigaxonin. Nevertheless, the underlying molecular mechanism remains obscure.

METHODS

Cell lines from 4 GAN patients and 4 controls were analyzed by iTRAQ.

RESULTS

Among the dysregulated proteins were ribosomal protein L29, ribosomal protein L37, galectin-1, glia-derived nexin, and aminopeptidase N. Also, nuclear proteins linked to formin-binding proteins were found to be dysregulated. Although the major role of gigaxonin is reported to be degradation of cytoskeleton-associated proteins, the amount of 76 structural cytoskeletal proteins was unaltered.

CONCLUSIONS

Several of the dysregulated proteins play a role in cytoskeletal reorganization. Based on these findings, we speculate that disturbed cytoskeletal regulation is responsible for formation of aggregates of intermediate filaments.

The amyloid precursor protein (APP) family members are key players in S-adenosylmethionine formation by MAT2A and modify BACE1 and PSEN1 gene expression-relevance for Alzheimer's disease

Central hallmark of Alzheimer's disease are senile plaques mainly composed of β-amyloid, which is a cleavage product of the amyloid precursor protein (APP). The physiological function of APP and its family members APLP1 and APLP2 is poorly understood. In order to fill this gap, we established a cell-culture based model with simultaneous knockdown of all members of the family. A comprehensive proteome study of the APP/APLP1/APLP2 knockdown cell lysates versus controls revealed significant protein abundance changes of more than 30 proteins. Targeted validation of selected candidates by immunoblotting supported the significant down-regulation of the methionine adenosyltransferase II, alpha (MAT2A) as well as of peroxiredoxin 4 in the knockdown cells. Moreover, MAT2A was significantly down-regulated at the mRNA level as well. MAT2A catalyzes the production of S-adenosylmethionine from methionine and ATP, which plays a pivotal role in the methylation of neurotransmitters, DNA, proteins, and lipids. MAT2A-dependent significant up-regulation of S-adenosylmethionine was also detectable in the knockdown cells compared with controls. Our results point to a role of the APP family proteins in cellular methylation mechanisms and fit to findings of disturbed S-adenosylmethionine levels in tissue and CSF of Alzheimer disease patients versus controls. Importantly, methylation plays a central role for neurotransmitter generation like acetylcholine pointing to a crucial relevance of our findings for Alzheimer's disease. In addition, we identified differential gene expression of BACE1 and PSEN1 in the knockdown cells, which is possibly a consequence of MAT2A deregulation and may indicate a self regulatory mechanism.

Important issues in planning a proteomics experiment: statistical considerations of quantitative proteomic data

Mass spectrometry is frequently used in quantitative proteomics to detect differentially regulated proteins. A very important but unfortunately oftentimes neglected part in detecting differential proteins is the statistical analysis. Data from proteomics experiments are usually high-dimensional and hence require profound statistical methods. It is especially important to already correctly design a proteomic experiment before it is conducted in the laboratory. Only this can ensure that the statistical analysis is capable of detecting truly differential proteins afterwards. This chapter thus covers aspects of both statistical planning and the actual analysis of quantitative proteomic experiments.

Search and decoy: the automatic identification of mass spectra

In recent years, the generation and interpretation of MS/MS spectra for the identification of peptides and proteins has matured to a frequently used automatic workflow in Proteomics. Several software solutions for the automated analysis of MS/MS spectra allow for high-throughput/high-performance analyses of complex samples. Related to MS/MS searches, target-decoy approaches have gained more and more popularity: in a "decoy" part of the search database nonexistent sequences mimic real sequences (the "target" sequences). With their help, the number of falsely identified peptides/proteins can be estimated after a search and the resulting protein list can be cut at a specified false discovery rate (FDR). This is an essential prerequisite for all quantitative approaches, as they rely on correct identifications. Especially the label-free approach "spectral counting"-gaining more and more popularity due to low costs and simplicity-depends directly on the correctness of peptide-spectrum matches (PSMs). This work's aim is to describe five popular search engines-especially their general properties regarding protein identification, but also their quantification abilities, if those go beyond spectral counting. By doing so, Proteomics researchers are enabled to compare their features and to choose an appropriate solution for their specific question. Furthermore, the search engines are applied to a spectrum data set generated from a complex sample with a Thermo LTQ Velos OrbiTrap (Thermo Fisher Scientific, Waltham, MA, USA). The results of the search engines are compared, e.g., regarding time requirements, peptides and proteins found, and the search engines' behavior using the decoy approach.

Construction of a chassis for hydrogen production: physiological and molecular characterization of a Synechocystis sp. PCC 6803 mutant lacking a functional bidirectional hydrogenase

Cyanobacteria are photosynthetic prokaryotes that are promising 'low-cost' microbial cell factories due to their simple nutritional requirements and metabolic plasticity, and the availability of tools for their genetic manipulation. The unicellular non-nitrogen-fixing Synechocystis sp. PCC 6803 is the best studied cyanobacterial strain and its genome was the first to be sequenced. The vast amount of physiological and molecular data available, together with a relatively small genome, makes Synechocystis suitable for computational metabolic modelling and to be used as a photoautotrophic chassis in synthetic biology applications. To prepare it for the introduction of a synthetic hydrogen producing device, a Synechocystis sp. PCC 6803 deletion mutant lacking an active bidirectional hydrogenase (ΔhoxYH) was produced and characterized at different levels: physiological, proteomic and transcriptional. The results showed that, under conditions favouring hydrogenase activity, 17 of the 210 identified proteins had significant differential fold changes in comparisons of the mutant with the wild-type. Most of these proteins are related to the redox and energy state of the cell. Transcriptional studies revealed that only six genes encoding those proteins exhibited significant differences in transcript levels. Moreover, the mutant exhibits similar growth behaviour compared with the wild-type, reflecting Synechocystis plasticity and metabolic adaptability. Overall, this study reveals that the Synechocystis ΔhoxYH mutant is robust and can be used as a photoautotrophic chassis for the integration of synthetic constructs, i.e. molecular constructs assembled from well characterized biological and/or synthetic parts (e.g. promoters, regulators, coding regions, terminators) designed for a specific purpose.

Expression of PTRF in PC-3 Cells modulates cholesterol dynamics and the actin cytoskeleton impacting secretion pathways

Expression of caveolin-1 is up-regulated in prostate cancer metastasis and is associated with aggressive recurrence of the disease. Intriguingly, caveolin-1 is also secreted from prostate cancer cell lines and has been identified in secreted prostasomes. Caveolin-1 is the major structural component of the plasma membrane invaginations called caveolae. Co-expression of the coat protein Polymerase I and transcript release factor (PTRF) is required for caveolae formation. We recently found that expression of caveolin-1 in the aggressive prostate cancer cell line PC-3 is not accompanied by PTRF, leading to noncaveolar caveolin-1 lipid rafts. Moreover, ectopic expression of PTRF in PC-3 cells sequesters caveolin-1 into caveolae. Here we quantitatively analyzed the effect of PTRF expression on the PC-3 proteome using stable isotope labeling by amino acids in culture and subcellular proteomics. We show that PTRF reduced the secretion of a subset of proteins including secreted proteases, cytokines, and growth regulatory proteins, partly via a reduction in prostasome secretion. To determine the cellular mechanism accounting for the observed reduction in secreted proteins we analyzed total membrane and the detergent-resistant membrane fractions. Our data show that PTRF expression selectively impaired the recruitment of actin cytoskeletal proteins to the detergent-resistant membrane, which correlated with altered cholesterol distribution in PC-3 cells expressing PTRF. Consistent with this, modulating cellular cholesterol altered the actin cytoskeleton and protein secretion in PC-3 cells. Intriguingly, several proteins that function in ER to Golgi trafficking were reduced by PTRF expression. Taken together, these results suggest that the noncaveolar caveolin-1 found in prostate cancer cells generates a lipid raft microenvironment that accentuates secretion pathways, possibly at the step of ER sorting/exit. Importantly, these effects could be modulated by PTRF expression.

The cDNA sequence of three hemocyanin subunits from the garden snail Helix lucorum

Hemocyanins are blue copper containing respiratory proteins residing in the hemolymph of many molluscs and arthropods. They can have different molecular masses and quaternary structures. Moreover, several molluscan hemocyanins are isolated with one, two or three isoforms occurring as decameric, didecameric, multidecameric or tubule aggregates. We could recently isolate three different hemocyanin isopolypeptides from the hemolymph of the garden snail Helix lucorum (HlH). These three structural subunits were named α(D)-HlH, α(N)-HlH and β-HlH. We have cloned and sequenced their cDNA which is the first result ever reported for three isoforms of a molluscan hemocyanin. Whereas the complete gene sequence of α(D)-HlH and β-HlH was obtained, including the 5' and 3' UTR, 180bp of the 5' end and around 900bp at the 3' end are missing for the third subunit. The subunits α(D)-HlH and β-HlH comprise a signal sequence of 19 amino acids plus a polypeptide of 3409 and 3414 amino acids, respectively. We could determine 3031 residues of the α(N)-HLH subunit. Sequence comparison with other molluscan hemocyanins shows that α(D)-HlH is more related to Aplysia californicum hemocyanin than to each of its own isopolypeptides. The structural subunits comprise 8 different functional units (FUs: a, b, c, d, e, f, g, h) and each functional unit possesses a highly conserved copper-A and copper-B site for reversible oxygen binding. Potential N-glycosylation sites are present in all three structural subunits. We confirmed that all three different isoforms are effectively produced and secreted in the hemolymph of H. lucorum by analyzing a tryptic digest of the purified native hemocyanin by MALDI-TOF and LC-FTICR mass spectrometry.

A guide through the computational analysis of isotope-labeled mass spectrometry-based quantitative proteomics data: an application study

BACKGROUND

Mass spectrometry-based proteomics has reached a stage where it is possible to comprehensively analyze the whole proteome of a cell in one experiment. Here, the employment of stable isotopes has become a standard technique to yield relative abundance values of proteins. In recent times, more and more experiments are conducted that depict not only a static image of the up- or down-regulated proteins at a distinct time point but instead compare developmental stages of an organism or varying experimental conditions.

RESULTS

Although the scientific questions behind these experiments are of course manifold, there are, nevertheless, two questions that commonly arise: 1) which proteins are differentially regulated regarding the selected experimental conditions, and 2) are there groups of proteins that show similar abundance ratios, indicating that they have a similar turnover? We give advice on how these two questions can be answered and comprehensively compare a variety of commonly applied computational methods and their outcomes.

CONCLUSIONS

This work provides guidance through the jungle of computational methods to analyze mass spectrometry-based isotope-labeled datasets and recommends an effective and easy-to-use evaluation strategy. We demonstrate our approach with three recently published datasets on Bacillus subtilis 12 and Corynebacterium glutamicum 3. Special focus is placed on the application and validation of cluster analysis methods. All applied methods were implemented within the rich internet application QuPE 4. Results can be found at http://qupe.cebitec.uni-bielefeld.de.

Proteomic and immunoproteomic characterization of a DIVA subunit vaccine against Actinobacillus pleuropneumoniae

Background

Protection of pigs by vaccination against Actinobacillus pleuropneumoniae, the causative agent of porcine pleuropneumonia, is hampered by the presence of 15 different serotypes. A DIVA subunit vaccine comprised of detergent-released proteins from A. pleuropneumoniae serotypes 1, 2 and 5 has been developed and shown to protect pigs from clinical symptoms upon homologous and heterologous challenge. This vaccine has not been characterized in-depth so far. Thus we performed i) mass spectrometry in order to identify the exact protein content of the vaccine and ii) cross-serotype 2-D immunoblotting in order to discover cross-reactive antigens. By these approaches we expected to gain results enabling us to argue about the reasons for the efficacy of the analyzed vaccine.

Results

We identified 75 different proteins in the vaccine. Using the PSORTb algorithm these proteins were classified according to their cellular localization. Highly enriched proteins are outer membrane-associated lipoproteins like OmlA and TbpB, integral outer membrane proteins like FrpB, TbpA, OmpA1, OmpA2, HgbA and OmpP2, and secreted Apx toxins. The subunit vaccine also contained large amounts of the ApxIVA toxin so far thought to be expressed only during infection. Applying two-dimensional difference gel electrophoresis (2-D DIGE) we showed different isoforms and variations in expression levels of several proteins among the strains used for vaccine production. For detection of cross-reactive antigens we used detergent released proteins of serotype 7. Sera of pigs vaccinated with the detergent-released proteins of serotypes 1, 2, and 5 detected seven different proteins of serotype 7, and convalescent sera of pigs surviving experimental infection with serotype 7 reacted with 13 different proteins of the detergent-released proteins of A. pleuropneumoniae serotypes 1, 2, and 5.

Conclusions

A detergent extraction-based subunit vaccine of A. pleuropneumoniae was characterized by mass spectrometry. It contained a large variety of immunogenic and virulence associated proteins, among them the ApxIVA toxin. The identification of differences in expression as well as isoform variation between the serotypes implied the importance of combining proteins of different serotypes for vaccine generation. This finding was supported by immunoblotting showing the induction of cross-reactive antibodies against several surface associated proteins in immunized animals.

Laboratory data and sample management for proteomics

Proteomic experiments can be difficult to handle because of the large amount of data in different formats that is generated. Samples need to be managed and generated, data needs to be integrated with samples and annotation information. A laboratory information management system (LIMS) can be used to overcome some of the data handling problems. In this chapter, we discuss the role of a LIMS in the proteomics laboratory, and show two step-by-step examples of usage of the Proteios Software Environment (ProSE) to handle two different proteomics workflows.

mzIdentML: an open community-built standard format for the results of proteomics spectrum identification algorithms

To deal with the data flood of current mass spectrometry methods, standard data formats are needed. The Proteomics Standards Initiative (PSI) of the Human Proteome Organisation (HUPO) develops open storage and transfer standards for and with the community. The Proteomics Informatics work group of the PSI has recently released an XML-based format to store the parameters and results of spectrum identification algorithms (the so-called search engines), which identify peptides and/or proteins from mass spectra. Here, this format called "mzIdentML" is described by giving principle design concepts and presenting examples of important use cases.

Bioinformatics challenges in mass spectrometry-driven proteomics

Mass spectrometry-based proteomics has become an essential part of the analytical toolbox of the life sciences. With the ability to identify and quantify hundreds to thousands of proteins in high throughput, the field has contributed its fair share to the data avalanche coming from the so-called omics fields. As a result, the challenges involved in processing and managing this flood of data have grown as well. This chapter will point out and discuss these challenges, starting from the processing of raw mass spectrometry data into peaks, over the identification of peptides and proteins, to the quantification of the identified molecules. Finally, the informatics aspects of the nascent field of targeted proteomics are outlined as well.

Using Laboratory Information Management Systems as central part of a proteomics data workflow

The organization and storage of proteomics data are challenging issues today and even more for the rising amount of information in the future. This review article describes the advantages of using Laboratory Information Management Systems (LIMS) in proteomics laboratories. Seven typical LIMS are explored in detail to describe their role in an even bigger interrelation. They are a central part of the proteomics data workflow, starting with data generation and ending with the publication in journals and repositories. Therefore, they enable community-wide data utilization and further Systems Biology discoveries.

Quantitative proteomics analysis of chondrogenic differentiation of C3H10T1/2 mesenchymal stem cells by iTRAQ labeling coupled with on-line two-dimensional LC/MS/MS

The chondrogenic potential of multipotent mesenchymal stem cells (MSCs) makes them a promising source for cell-based therapy of cartilage defects; however, the exact intracellular molecular mechanisms of chondrogenesis as well as self-renewal of MSCs remain largely unknown. To gain more insight into the underlying molecular mechanisms, we applied isobaric tag for relative and absolute quantitation (iTRAQ) labeling coupled with on-line two-dimensional LC/MS/MS technology to identify proteins differentially expressed in an in vitro model for chondrogenesis: chondrogenic differentiation of C3H10T1/2 cells, a murine embryonic mesenchymal cell line, was induced by micromass culture and 100 ng/ml bone morphogenetic protein 2 treatment for 6 days. A total of 1756 proteins were identified with an average false discovery rate <0.21%. Linear regression analysis of the quantitative data gave strong correlation coefficients: 0.948 and 0.923 for two replicate two-dimensional LC/MS/MS analyses and 0.881, 0.869, and 0.927 for three independent iTRAQ experiments, respectively (p < 0.0001). Among 1753 quantified proteins, 100 were significantly altered (95% confidence interval), and six of them were further validated by Western blotting. Functional categorization revealed that the 17 up-regulated proteins mainly comprised hallmarks of mature chondrocytes and enzymes participating in cartilage extracellular matrix synthesis, whereas the 83 down-regulated were predominantly involved in energy metabolism, chromatin organization, transcription, mRNA processing, signaling transduction, and cytoskeleton; except for a number of well documented proteins, the majority of these altered proteins were novel for chondrogenesis. Finally, the biological roles of BTF3l4 and fibulin-5, two novel chondrogenesis-related proteins identified in the present study, were verified in the context of chondrogenic differentiation. These data will provide valuable clues for our better understanding of the underlying mechanisms that modulate these complex biological processes and assist in the application of MSCs in cell-based therapy for cartilage regeneration.