From proteomic multimarker profiling to interesting proteins: thymosin-β(4) and kininogen-1 as new potential biomarkers for inflammatory hepatic lesions

Despite tremendous efforts in disclosing the pathophysiological and epidemiological factors associated with liver fibrogenesis, non-invasive diagnostic measures to estimate the clinical outcome and progression of liver fibrogenesis are presently limited. Therefore, there is a mandatory need for methodologies allowing the reasonable and reliable assessment of the severity and/or progression of hepatic fibrogenesis. We here performed proteomic serum profiling by matrix-assisted laser desorption ionization time-of-flight mass spectrometry in 179 samples of patients chronically infected with hepatitis C virus and 195 control sera. Multidimensional analysis of spectra allowed the definition of algorithms capable to distinguish class-specific protein expression profiles in serum samples. Overall about 100 peaks could be detected per single spectrum. Different algorithms including protein peaks in the range of 2000 and 10,000 Da were generated after pre-fractionation on a weak cation exchange surface. A specificity of 93% with a sensitivity of 86% as mean of the test set results was found, respectively. The nature of three of these protein peaks that belonged to kininogen-1 and thymosin-β(4) was further analysed by tandem mass spectrometry (MS)/MS. We further found that kininogen-1 mRNA was significantly down-regulated in cirrhotic livers. We have identified kininogen-1 and thymosin-β(4) as potential new biomarkers for human chronic hepatitis C and conclude that serum profiling is a reliable technique to identify hepatitis-associated expression patterns. Based on the high throughput capability, the identified differential protein panel may serve as a diagnostic marker and warrants further validation in larger cohorts.     

Enrichment of low molecular weight fraction of serum for MS analysis of peptides associated with hepatocellular carcinoma

A challenging aspect of biomarker discovery in serum is the interference of abundant proteins with identification of disease-related proteins and peptides. This study describes enrichment of serum by denaturing ultrafiltration, which enables an efficient profiling and identification of peptides up to 5 kDa. We consistently detect several hundred peptide-peaks in MALDI-TOF and SELDI-TOF spectra of enriched serum. The sample preparation is fast and reproducible with an average CV for all 276 peaks in the MALDI-TOF spectrum of 11%. Compared to unenriched serum, the number of peaks in enriched spectra is 4 times higher at an S/N ratio of 5 and 20 times higher at an S/N ratio of 10. To demonstrate utility of the methods, we compared 20 enriched sera of patients with hepatocellular carcinoma (HCC) and 20 age-matched controls using MALDI-TOF. The comparison of 332 peaks at p < 0.001 identified 45 differentially abundant peaks that classified HCC with 90% accuracy in this small pilot study. Direct TOF/TOF sequencing of the most abundant peptide matches with high probability des-Ala-fibrinopeptide A. This study shows that enrichment of the low molecular weight fraction of serum facilitates an efficient discovery of peptides that could serve as biomarkers for detection of HCC as well as other diseases.     

Identification of potential serum biomarkers for Wilms tumor after excluding confounding effects of common systemic inflammatory factors

Wilms tumor is the most common pediatric tumor of the kidney. Previous studies have identified several serum biomarkers for Wilms tumor; however, they lack sufficient specificity and may not adequately distinguish Wilms tumor from confounding conditions. To date, no specific protein biomarker has been confirmed for this pediatric tumor. To identify novel serum biomarkers for Wilms tumor, we used proteomic technologies to perform protein profiling of serum samples from pre-surgery and post-surgery patients with Wilms tumor and healthy controls. Some common systemic inflammatory factors were included to control for systemic inflammation. By comparing protein peaks among the three groups of sera, we identified two peaks (11,526 and 4,756 Da) showing significant differential expression not only between pre-surgery and control sera but also between pre-surgery and post-surgery sera. These two peaks were identified as serum amyloid A1 (SAA1) and apolipoprotein C-III (APO C-III). Western blot analysis confirmed that both proteins were expressed at higher levels in pre-surgery sera than in post-surgery and control sera. Using the method of leave-1-out for cross detection, we demonstrate that detection of these two candidate biomarkers had high sensitivity and specificity in discriminating pre-surgery sera from post-surgery and normal control sera. Taken together, these findings suggest that SAA1 and APO C-III are two potential biomarkers for Wilms tumor.     

Characterizing the reproducibility of a protein profiling method for the analysis of mouse bronchoalveolar lavage fluid

The detection of biomarkers in biological fluids has been advanced by the introduction of mass spectrometry screening methods such as matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOFMS), which enables the detection of the presence and the molecular mass of proteins in unfractionated mixtures. The generation of reproducible mass spectra over the course of an experiment is vital in obtaining data in which differences in protein profiles between diseased and healthy states can be assessed correctly. We have developed a protocol to automate the collection of protein profiling data from a large number of samples using MALDI-TOFMS, and we used these samples to characterize the technical reproducibility of the method. This protocol has been used for the analysis of proteins found in bronchoalveolar lavage fluid samples from mice with the ultimate goal of enabling the discovery of differential expression patterns predictive of the development of chronic obstructive pulmonary disease. Samples were purified using magnetic bead-based technology and analyzed on an AnchorChip target plate. Our results demonstrate that the number of peaks detected reproducibly decreases significantly as sample size increases, which motivates the need for technical replicates to be explicitly included in the analysis of MALDI-TOF-based protein profiling studies.     

Novel statistical approaches to identify risk factors for soil-transmitted helminth infection in Timor-Leste

Soil-transmitted helminths (STHs) are parasitic intestinal worms that infect almost a fifth of the global population. Sustainable control of STHs requires understanding the complex interaction of factors contributing to transmission. Identifying risk factors has mainly relied on logistic regression models where the underlying assumption of independence between variables is not always satisfied. Previously demonstrated risk factors including water, sanitation and hygiene (WASH) access and behaviours, and socioeconomic status are intrinsically linked. Similarly, environmental factors including climate, soil and land attributes are often strongly correlated. Alternative methods such as recursive partitioning and Bayesian networks can handle correlated variables, but there are no published studies comparing these methods with logistic regression in the context of STH risk factor analysis. Baseline cross-sectional data from school-aged children in the (S)WASH-D for Worms study were used to compare risk factors identified from modelling the same data using three different statistical techniques. Outcomes of interest were infection with Ascaris spp. and any hookworm species (Necator americanus, Ancylostoma duodenale, and Ancylostoma ceylanicum). Mixed-effects logistic regression identified the fewest risk factors. Recursive partitioning identified the most WASH and demographic risk factors, while Bayesian networks identified the most environmental risk factors. Recursive partitioning produced classification trees that visualised potentially at-risk population sub-groups. Bayesian networks helped visualise relationships between variables and enabled interactive modelling of outcomes based on different scenarios for the predictor variables of interest. Model performance was similar across all techniques. Risk factors identified across all techniques were vegetation for Ascaris spp., and cleaning oneself with water after defecating for hookworm. This study adds to the limited body of evidence exploring alternative data modelling approaches in identifying risk factors for STH infections. Our findings suggest these approaches can provide novel insights for more robust interpretation.     

Analysis of Human Proteome Organization Plasma Proteome Project (HUPO PPP) reference specimens using surface enhanced laser desorption/ionization-time of flight (SELDI-TOF) mass spectrometry: multi-institution correlation of spectra and identification of biomarkers

We report on a multicenter analysis of HUPO reference specimens using SELDI-TOF MS. Eight sites submitted data obtained from serum and plasma reference specimen analysis. Spectra from five sites passed preliminary quality assurance tests and were subjected to further analysis. Intralaboratory CVs varied from 15 to 43%. A correlation coefficient matrix generated using data from these five sites demonstrated high level of correlation, with values >0.7 on 37 of 42 spectra. More than 50 peaks were differentially present among the various sample types, as observed on three chip surfaces. Additionally, peaks at approximately 9200 and approximately 15,950 m/z were present only in select reference specimens. Chromatographic fractionation using anion-exchange, membrane cutoff, and reverse phase chromatography, was employed for protein purification of the approximately 9200 m/z peak. It was identified as the haptoglobin alpha subunit after peptide mass fingerprinting and high-resolution MS/MS analysis. The differential expression of this protein was confirmed by Western blot analysis. These pilot studies demonstrate the potential of the SELDI platform for reproducible and consistent analysis of serum/plasma across multiple sites and also for targeted biomarker discovery and protein identification. This approach could be exploited for population-based studies in all phases of the HUPO PPP.     

Characterization of data analysis methods for information recovery from metabolic 1H NMR spectra using artificial complex mixtures

The assessment of data analysis methods in ¹H NMR based metabolic profiling is hampered owing to a lack of knowledge of the exact sample composition. In this study, an artificial complex mixture design comprising two artificially defined groups designated normal and disease, each containing 30 samples, was implemented using 21 metabolites at concentrations typically found in human urine and having a realistic distribution of inter-metabolite correlations. These artificial mixtures were profiled by ¹H NMR spectroscopy and used to assess data analytical methods in the task of differentiating the two conditions. When metabolites were individually quantified, volcano plots provided an excellent method to track the effect size and significance of the change between conditions. Interestingly, the Welch t test detected a similar set of metabolites changing between classes in both quantified and spectral data, suggesting that differential analysis of ¹H NMR spectra using a false discovery rate correction, taking into account fold changes, is a reliable approach to detect differential metabolites in complex mixture studies. Various multivariate regression methods based on partial least squares (PLS) were applied in discriminant analysis mode. The most reliable methods in quantified and spectral ¹H NMR data were PLS and RPLS linear and logistic regression respectively. A jackknife based strategy for variable selection was assessed on both quantified and spectral data and results indicate that it may be possible to improve on the conventional Orthogonal-PLS methodology in terms of accuracy and sensitivity. A key improvement of our approach consists of objective criteria to select significant signals associated with a condition that provides a confidence level on the discoveries made, which can be implemented in metabolic profiling studies.     

Evaluation of matrix-assisted laser desorption/ionization-time of flight mass spectrometry proteomic profiling: identification of alpha 2-HS glycoprotein B-chain as a biomarker of diet

Biomarkers have the potential to impact a wide range of public health concerns, including early detection of diseases, drug discovery, and improved accuracy of monitoring effects of interventions. Given new technological developments, broad-based screening approaches will likely advance biomarker discovery at an accelerated pace. Matrix-assisted laser desorption/ionization-time of flight mass spectrometry (MALDI-TOF MS) allows for the elucidation of individual protein masses from a complex mixture with high throughput. We have developed a method for identifying serum biomarkers using MALDI-TOF and statistical analysis. However, before applying this approach to screening of complex diseases, we evaluated the approach in a controlled dietary intervention study. In this study, MALDI-TOF spectra were generated using samples from a randomized controlled trial. During separate feeding periods, 38 participants ate a basal diet devoid of fruits and vegetables and a basal diet supplemented with cruciferous (broccoli) family vegetables. Serum samples were obtained at the end of each 7-day feeding period and treated to remove large, abundant proteins. MALDI-TOF spectra were analyzed using peak picking algorithms and logistic regression models. Our bioinformatics methods identified two significant peaks at m/z values of 2740 and 1847 that could classify participants based on diet (basal vs. cruciferous) with 76% accuracy. The 2740 m/z peak was identified as the B-chain of alpha 2-HS glycoprotein, a serum protein previously found to vary with diet and be involved in insulin resistance and immune function.     

Combined experimental and statistical strategy for mass spectrometry based serum protein profiling for diagnosis of breast cancer: a case-control study

Serum protein profiling by mass spectrometry is a promising method for early detection of cancer. We have implemented a combined strategy based on matrix-assisted laser desorption ionization mass spectrometry (MALDI MS) and statistical data analysis for serum protein profiling and applied it in a well-described breast cancer case-control study. A rigorous sample collection protocol ensured high quality specimen and reduced bias from preanalytical factors. Preoperative serum samples obtained from 48 breast cancer patients and 28 controls were used to generate MALDI MS protein profiles. A total of nine mass spectrometric protein profiles were obtained for each serum sample. A total of 533 common peaks were defined and represented a 'reference protein profile'. Among these 533 common peaks, we identified 72 peaks exhibiting statistically significant intensity differences ( p < 0.01) between cases and controls. A diagnostic rule based on these 72 mass values was constructed and exhibited a cross-validated sensitivity and specificity of approximately 85% for the detection of breast cancer. With this method, it was possible to distinguish early stage cancers from controls without major loss of sensitivity and specificity. We conclude that optimized serum sample handling and mass spectrometry data acquisition strategies in combination with statistical analysis provide a viable platform for serum protein profiling in cancer diagnosis.     

Serum Biomarkers of Keshan Disease Assessed Using a Protein Profiling Approach Based on ClinProt Technique

The etiology of Keshan disease (KD), an endemic myocardiopathy in regions of China, is largely unknown. To show the protein changes in serum from KD patients versus controls and idiopathic dilated cardiomyopathy (IDCM) and to search specific biological markers for differential diagnosis for KD. Serum of 65 patients with KD was compared with 29 patients with IDCM, 62 controls from KD areas and 28 controls from non-KD areas by ClinProt/MALDI-ToF technique. The genetic algorithm, quick classifier algorithm and supervised neural network algorithm methods were used to screen marker proteins and establish diagnostic model. Thirty-four differential peaks were identified in KD patients compared with the healthy controls from non-KD areas. Thirty-eight differentially peaks were identified in KD patients and controls from KD areas; and sixty-seven differentially peaks were identified in patients with KD and patients with IDCM. We believe that marker protein peaks screened in KD patients, healthy controls and IDCM patients may provide clues for the differential diagnosis and treatment of KD.     

Circulating miR-19a and miR-205 in Serum May Predict the Sensitivity of Luminal A Subtype of Breast Cancer Patients to Neoadjuvant Chemotherapy with Epirubicin Plus Paclitaxel

Background

The luminal A subtype of breast cancer has a good prognosis and is sensitive to endocrine therapy but is less sensitive to chemotherapy. It is necessary to identify biomarkers to predict chemosensitivity and avoid over-treatment. We hypothesized that miRNAs in the serum might be associated with chemosensitivity.

Methods

Sixty-eight breast cancer patients received neoadjuvant chemotherapy with epirubicin plus paclitaxel. The serum of the patients was collected before chemotherapy and stored at −80°C. The samples were classified into two groups in term of the chemosensitivity. We identified the differential expression patterns of miRNAs between the chemotherapy sensitive and resistant groups using microRNA profiling. Four miRNAs that were differentially expressed between the two groups were further validated in another 56 samples. We created a model fitting formula and a receiver operating characteristics (ROC) curve using logistic regression analysis to evaluate the prediction potency.

Results

We identified 8 miRNAs differentially expressed between the two groups: 6 miRNAs were up-regulated, and 2 miRNAs were down-regulated in the resistant group compared with the sensitive group. The expression of miR-19a and miR-205 were determined to have significant differences between the two groups (P<0.05). A predictive model of these two miRNAs was created by the logistic regression analysis. The probability of this model was 89.71%. Based on the ROC curve, the specificity was 75.00%, and the sensitivity was 81.25%.

Conclusions

The combination of miR-19a and miR-205 in the serum may predict the chemosensitivity of luminal A subtype of breast cancer to epirubicin plus paclitaxel neoadjuvant chemotherapy.     

Analysis of MALDI-TOF mass spectrometry data for discovery of peptide and glycan biomarkers of hepatocellular carcinoma

This paper presents computational methods to analyze MALDI-TOF mass spectrometry data for quantitative comparison of peptides and glycans in serum. The methods are applied to identify candidate biomarkers in serum samples of 203 participants from Egypt; 73 hepatocellular carcinoma (HCC) cases, 52 patients with chronic liver disease (CLD) consisting of cirrhosis and fibrosis cases, and 78 population controls. Two complementary sample preparation methods were applied prior to generating mass spectra: (1) low molecular weight (LMW) enrichment of each serum sample was carried out for MALDI-TOF quantification of peptides, and (2) glycans were enzymatically released from proteins in each serum sample and permethylated for MALDI-TOF quantification of glycans. A peak selection algorithm was applied to identify the most useful peptide and glycan peaks for accurate detection of HCC cases from high-risk population of patients with CLD. In addition to global peaks selected by the whole population based approach, where identically labeled patients are treated as a single group, subgroup-specific peaks were identified by searching for peaks that are differentially abundant in a subgroup of patients only. The peak selection process was preceded by peak screening, where we eliminated peaks that have significant association with covariates such as age, gender, and viral infection based on the peptide and glycan spectra from population controls. The performance of the selected peptide and glycan peaks was evaluated in terms of their ability in detecting HCC cases from patients with CLD in a blinded validation set and through the cross-validation method. Finally, we investigated the possibility of using both peptides and glycans in a panel to enhance the diagnostic capability of these candidate markers. Further evaluation is needed to examine the potential clinical utility of the candidate peptide and glycan markers identified in this study.     

Factors influencing soay sheep survival: a Bayesian analysis

This article presents a Bayesian analysis of mark-recapture-recovery data on Soay sheep. A reversible jump Markov chain Monte Carlo technique is used to determine age classes of common survival, and to model the survival probabilities in those classes using logistic regression. This involves environmental and individual covariates, as well as random effects. Auxiliary variables are used to impute missing covariates measured on individual sheep. The Bayesian approach suggests different models from those previously obtained using classical statistical methods. Following model averaging, features that were not previously detected, and which are of ecological importance, are identified.     

Multimarker Proteomic Profiling for the Prediction of Cardiovascular Mortality in Patients with Chronic Heart Failure

Risk stratification of patients with systolic chronic heart failure (HF) is critical to better identify those who may benefit from invasive therapeutic strategies such as cardiac transplantation. Proteomics has been used to provide prognostic information in various diseases. Our aim was to investigate the potential value of plasma proteomic profiling for risk stratification in HF. A proteomic profiling using surface enhanced laser desorption ionization - time of flight - mass spectrometry was performed in a case/control discovery population of 198 patients with systolic HF (left ventricular ejection fraction <45%): 99 patients who died from cardiovascular cause within 3 years and 99 patients alive at 3 years. Proteomic scores predicting cardiovascular death were developed using 3 regression methods: support vector machine, sparse partial least square discriminant analysis, and lasso logistic regression. Forty two ion m/z peaks were differentially intense between cases and controls in the discovery population and were used to develop proteomic scores. In the validation population, score levels were higher in patients who subsequently died within 3 years. Similar areas under the curves (0.66 – 0.68) were observed for the 3 methods. After adjustment on confounders, proteomic scores remained significantly associated with cardiovascular mortality. Use of the proteomic scores allowed a significant improvement in discrimination of HF patients as determined by integrated discrimination improvement and net reclassification improvement indexes. In conclusion, proteomic analysis of plasma proteins may help to improve risk prediction in HF patients.     

Proteomic serum biomarkers for neuromuscular diseases

Introduction: The clinical evaluation of neuromuscular symptoms often includes the assessment of altered blood proteins or changed enzyme activities. However, the blood concentration of many muscle-derived serum markers is not specific for different neuromuscular disorders and also shows alterations in the course of these diseases. Thus, the establishment of more reliable biomarker signatures for improved muscle diagnostics is required.

Areas covered: To address the lack of muscle disease-specific marker molecules, mass spectrometry-based proteomics was applied to the systematic identification and biochemical characterization of new serum biomarker candidates. This article outlines serum proteomics in relation to neuromuscular disorders and reviews the bioanalytical results from recent proteomic profiling studies of representative neuromuscular disorders, including motor neuron disease, muscular dystrophies and sarcopenia of old age. Pathophysiological changes in the skeletal muscle proteome are reflected by serum alterations in a variety of sarcomeric proteins, metabolic enzymes and signaling proteins.

Expert commentary: Based on the proteomic identification of actively secreted or passively released skeletal muscle proteins following pathophysiological insults, new biomarker candidates can now be used to develop liquid biopsy procedures for superior diagnostic approaches, design novel prognostic tools and establish more reliable methods for the systematic evaluation of experimental therapies to treat neuromuscular disease.     

Ryanodine receptors, a family of intracellular calcium ion channels, are expressed throughout early vertebrate development

Background

We developed a new version of the open source software package Peptrix that can yet compare large numbers of Orbitrap? LC-MS data. The peptide profiling results for Peptrix on MS1 spectra were compared with those obtained from a small selection of open source and commercial software packages: msInspect, Sieve? and Progenesis?. The properties compared in these packages were speed, total number of detected masses, redundancy of masses, reproducibility in numbers and CV of intensity, overlap of masses, and differences in peptide peak intensities. Reproducibility measurements were taken for the different MS1 software applications by measuring in triplicate a complex peptide mixture of immunoglobulin on the Orbitrap? mass spectrometer. Values of peptide masses detected from the high intensity peaks of the MS1 spectra by peptide profiling were verified with values of the MS2 fragmented and sequenced masses that resulted in protein identifications with a significant score.

Findings

Peptrix finds about the same number of peptide features as the other packages, but peptide masses are in some cases approximately 5 to 10 times less redundant present in the peptide profile matrix. The Peptrix profile matrix displays the largest overlap when comparing the number of masses in a pair between two software applications. The overlap of peptide masses between software packages of low intensity peaks in the spectra is remarkably low with about 50% of the detected masses in the individual packages. Peptrix does not differ from the other packages in detecting 96% of the masses that relate to highly abundant sequenced proteins. MS1 peak intensities vary between the applications in a non linear way as they are not processed using the same method.

Conclusions

Peptrix is capable of peptide profiling using Orbitrap? files and finding differential expressed peptides in body fluid and tissue samples. The number of peptide masses detected in Orbitrap? files can be increased by using more MS1 peptide profiling applications, including Peptrix, since it appears from the comparison of Peptrix with the other applications that all software packages have likely a high false negative rate of low intensity peptide peaks (missing peptides).     

A well‐characterised peak identification list of MALDI MS profile peaks for human blood serum

MALDI MS profiling, using easily available body fluids such as blood serum, has attracted considerable interest for its potential in clinical applications. Despite the numerous reports on MALDI MS profiling of human serum, there is only scarce information on the identity of the species making up these profiles, particularly in the mass range of larger peptides. Here, we provide a list of more than 90 entries of MALDI MS profile peak identities up to 10 kDa obtained from human blood serum. Various modifications such as phosphorylation were detected among the peptide identifications. The overlap with the few other MALDI MS peak lists published so far was found to be limited and hence our list significantly extends the number of identified peaks commonly found in MALDI MS profiling of human blood serum.     

Toward the proteomic identification of biomarkers for the prediction of HBV related hepatocellular carcinoma

Early detection is a key step for effective intervention of hepatocellular carcinoma (HCC), the lack of sensitive and specific biomarkers is a major reason for the high rate of HCC-related mortality. This report described an integrated strategy by combining SELDI-ProteinChip, sophisticated algorithm analysis, acetonitrile (ACN) pre-treatment and two-dimensional electrophoresis (2DE)-peptide mass fingerprinting (PMF) techniques to identify serological markers for the prediction of HBV-related HCC. Proteomic profiling of three groups of serum specimens from HBV-related HCC (50 cases), HBV infection (45 cases), and normal subjects (30 cases) was conducted by using SELDI-ProteinChip system and the resulting different protein peaks were subjected to stepwise statistical analyses. Three most discriminatory peaks at 5890, 11615, and 11724 Da, respectively, were screened out from the statistical algorithm and a predictive model based on the three peaks was constructed and tested using the newly enrolled serum samples. 2DE was applied to separate and compare the serum samples that were pre-treated by ACN precipitation. The protein spots obviously intensified in HCC sera in the 2DE region of 12 kDa were identified by PMF to be serum SAA, which was validated by SELDI-TOF spectra of HCC sera after immunoprecipitation using anti-SAA antibody and by Western blot experiments. Given the fact that SAA is not a specific biomarker, further attempt is being made to identify the other two most discriminatory peaks to realize the possibility of using the predictive model for HCC surveillance and prediction.     

A bayesian mixture model for comparative spectral count data in shotgun proteomics

Recent developments in mass-spectrometry-based shotgun proteomics, especially methods using spectral counting, have enabled large-scale identification and differential profiling of complex proteomes. Most such proteomic studies are interested in identifying proteins, the abundance of which is different under various conditions. Several quantitative methods have recently been proposed and implemented for this purpose. Building on some techniques that are now widely accepted in the microarray literature, we developed and implemented a new method using a Bayesian model to calculate posterior probabilities of differential abundance for thousands of proteins in a given experiment simultaneously. Our Bayesian model is shown to deliver uniformly superior performance when compared with several existing methods.     

parallelMCMCcombine: An R Package for Bayesian Methods for Big Data and Analytics

Recent advances in big data and analytics research have provided a wealth of large data sets that are too big to be analyzed in their entirety, due to restrictions on computer memory or storage size. New Bayesian methods have been developed for data sets that are large only due to large sample sizes. These methods partition big data sets into subsets and perform independent Bayesian Markov chain Monte Carlo analyses on the subsets. The methods then combine the independent subset posterior samples to estimate a posterior density given the full data set. These approaches were shown to be effective for Bayesian models including logistic regression models, Gaussian mixture models and hierarchical models. Here, we introduce the R package parallelMCMCcombine which carries out four of these techniques for combining independent subset posterior samples. We illustrate each of the methods using a Bayesian logistic regression model for simulation data and a Bayesian Gamma model for real data; we also demonstrate features and capabilities of the R package. The package assumes the user has carried out the Bayesian analysis and has produced the independent subposterior samples outside of the package. The methods are primarily suited to models with unknown parameters of fixed dimension that exist in continuous parameter spaces. We envision this tool will allow researchers to explore the various methods for their specific applications and will assist future progress in this rapidly developing field.