LFQRatio: A Normalization Method to Decipher Quantitative Proteome Changes in Microbial Coculture Systems

The value of synthetic microbial communities in biotechnology is gaining traction due to their ability to undertake more complex metabolic tasks than monocultures. However, a thorough understanding of strain interactions, productivity, and stability is often required to optimize growth and scale up cultivation. Quantitative proteomics can provide valuable insights into how microbial strains adapt to changing conditions in biomanufacturing. However, current workflows and methodologies are not suitable for simple artificial coculture systems where strain ratios are dynamic. Here, we established a workflow for coculture proteomics using an exemplar system containing two members, Azotobacter vinelandii and Synechococcus elongatus. Factors affecting the quantitative accuracy of coculture proteomics were investigated, including peptide physicochemical characteristics such as molecular weight, isoelectric point, hydrophobicity, and dynamic range as well as factors relating to protein identification such as varying proteome size and shared peptides between species. Different quantification methods based on spectral counts and intensity were evaluated at the protein and cell level. We propose a new normalization method, named "LFQRatio", to reflect the relative contributions of two distinct cell types emerging from cell ratio changes during cocultivation. LFQRatio can be applied to real coculture proteomics experiments, providing accurate insights into quantitative proteome changes in each strain.

The HLA-B*57:01 allele corresponds to a very large MHC haploblock likely explaining its massive effect for HIV-1 elite control

Introduction

We have reanalyzed the genomic data of the International Collaboration for the Genomics of HIV (ICGH), centering on HIV-1 Elite Controllers.

Methods

We performed a genome-wide Association Study comparing 543 HIV Elite Controllers with 3,272 uninfected controls of European descent. Using the latest database for imputation, we analyzed 35,552 Single Nucleotide Polymorphisms (SNPs) within the Major Histocompatibility Complex (MHC) region.

Results

Our analysis identified 2,626 SNPs significantly associated (p<5. 10-8) with elite control of HIV-1 infection, including well-established MHC signals such as the rs2395029-G allele which tags HLA-B*57:01. A thorough investigation of SNPs in linkage disequilibrium with rs2395029 revealed an extensive haploblock spanning 1.9 megabases in the MHC region tagging HLA-B*57:01, comprising 379 SNP alleles impacting 72 genes. This haploblock contains damaging variations in proteins like NOTCH4 and DXO and is also associated with a strong differential pattern of expression of multiple MHC genes such as HLA-B, MICB, and ZBTB12. The study was expanded to include two cohorts of seropositive African-American individuals, where a haploblock tagging the HLA-B*57:03 allele was similarly associated with control of viral load. The mRNA expression profile of this haploblock in African Americans closely mirrored that in the European cohort.

Discussion

These findings suggest that additional molecular mechanisms beyond the conventional antigen-presenting role of class I HLA molecules may contribute to the observed influence of HLA-B*57:01/B*57:03 alleles on HIV-1 elite control. Overall, this study has uncovered a large haploblock associated with HLA-B*57 alleles, providing novel insights into their massive effect on HIV-1 elite control.

Gene expression profiling of peripheral blood mononuclear cells from women with cervical lesions reveals new markers of cancer

Cervical cancer (CC) is a multifactorial disease of which human papillomavirus (HPV) is the main etiological agent. Despite cervical Pap smear screening and anti‑HPV vaccination, CC remains a major public health issue. Identification of specific gene expression signatures in the blood could allow better insight into the immune response of CC and could provide valuable information for the development of novel biomarkers. The present study performed a transcriptomic analysis of peripheral blood mononuclear cells (PBMCs) from Senegalese patients with CC (n=31), low‑grade cervical intraepithelial neoplasia (CIN1; n=27) and from healthy control (CTR) subjects (n=29). Individuals in the CIN1 and CTR groups exhibited similar patterns in gene expression. A total of 182 genes were revealed to be differentially expressed in patients with CC compared with individuals in the CIN1 and CTR groups. The IL1R2, IL18R1, MMP9 and FKBP5 genes were the most upregulated, whereas the T‑cell receptor α gene TRA was the most downregulated in the CC group compared with in the CIN1 and CTR groups. The pathway enrichment analysis of the differentially expressed genes revealed pathways directly and indirectly linked to inflammation. To the best of our knowledge, the present study is the first large transcriptomic study on CC performed using PBMCs from African women; the results revealed the involvement of genes and pathways related to inflammation, most notably the IL‑1 pathway, and the involvement of downregulation of the T‑cell receptor α, a key component of the immune response. Several of the stated genes have already been reported in other cancer studies as putative blood biomarkers, thus reinforcing the requirement for deeper investigation. These findings may aid in the development of innovative clinical biomarkers for CC prevention and should be further replicated in other populations.

Linear epitope mapping of the humoral response against SARS-CoV-2 in two independent African cohorts

Profiling of the antibody responses to severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) proteins in African populations is scarce. Here, we performed a detailed IgM and IgG epitope mapping study against 487 peptides covering SARS-CoV-2 wild-type structural proteins. A panel of 41 pre-pandemic and 82 COVID-19 RT-PCR confirmed sera from Madagascar and Senegal were used. We found that the main 36 immunodominant linear epitopes identified were (i) similar in both countries, (ii) distributed mainly in the Spike and the Nucleocapsid proteins, (iii) located outside the RBD and NTD regions where most of the reported SARS-CoV-2 variant mutations occur, and (iv) identical to those reported in European, North American, and Asian studies. Within the severe group, antibody levels were inversely correlated with the viral load. This first antibody epitope mapping study performed in patients from two African countries may be helpful to guide rational peptide-based diagnostic assays or vaccine development.

Different Pigmentation Risk Loci for High-Risk Monosomy 3 and Low-Risk Disomy 3 Uveal Melanomas

BACKGROUND

Uveal melanoma (UM), a rare malignant tumor of the eye, is predominantly observed in populations of European ancestry. UMs carrying a monosomy 3 (M3) frequently relapse mainly in the liver, whereas UMs with disomy 3 (D3) are associated with more favorable outcome. Here, we explored the UM genetic predisposition factors in a large genome-wide association study (GWAS) of 1142 European UM patients and 882 healthy controls .

METHODS

We combined 2 independent datasets (Global Screening Array) with the dataset described in a previously published GWAS in UM (Omni5 array), which were imputed separately and subsequently merged. Patients were stratified according to their chromosome 3 status, and identified UM risk loci were tested for differential association with M3 or D3 subgroups. All statistical tests were 2-sided.

RESULTS

We recapitulated the previously identified risk locus on chromosome 5 on CLPTM1L (rs421284: odds ratio [OR] =1.58, 95% confidence interval [CI] = 1.35 to 1.86; P = 1.98 × 10-8) and identified 2 additional risk loci involved in eye pigmentation: IRF4 locus on chromosome 6 (rs12203592: OR = 1.76, 95% CI = 1.44 to 2.16; P = 3.55 × 10-8) and HERC2 locus on chromosome 15 (rs12913832: OR= 0.57, 95% CI = 0.48 to 0.67; P = 1.88 × 10-11). The IRF4 rs12203592 single-nucleotide polymorphism was found to be exclusively associated with risk for the D3 UM subtype (ORD3 = 2.73, 95% CI = 1.87 to 3.97; P = 1.78 × 10-7), and the HERC2 rs12913832 single-nucleotide polymorphism was exclusively associated with risk for the M3 UM subtype (ORM3 = 2.43, 95% CI = 1.79 to 3.29; P = 1.13 × 10-8). However, the CLPTM1L risk locus was equally statistically significant in both subgroups.

CONCLUSIONS

This work identified 2 additional UM risk loci known for their role in pigmentation. Importantly, we demonstrate that UM tumor biology and metastatic potential are influenced by patients' genetic backgrounds.

HLA-check: evaluating HLA data from SNP information

Background

The major histocompatibility complex (MHC) region of the human genome, and specifically the human leukocyte antigen (HLA) genes, play a major role in numerous human diseases. With the recent progress of sequencing methods (eg, Next-Generation Sequencing, NGS), the accurate genotyping of this region has become possible but remains relatively costly. In order to obtain the HLA information for the millions of samples already genotyped by chips in the past ten years, efficient bioinformatics tools, such as SNP2HLA or HIBAG, have been developed that infer HLA information from the linkage disequilibrium existing between HLA alleles and SNP markers in the MHC region.

Results

In this study, we first used ShapeIT and Impute2 to implement an imputation method akin to SNP2HLA and found a comparable quality of imputation on a European dataset. More importantly, we developed a new tool, HLA-check, that allows for the detection of aberrant HLA allele calling with regard to the SNP genotypes in the region. Adding this tool to the HLA imputation software increases dramatically their accuracy, especially for HLA class I genes.

Conclusion

Overall, HLA-check was able to identify a limited number of implausible HLA typings (less than 10%) in a population, and these samples can then either be removed or be retyped by NGS for HLA association analysis.

A new 3p25 locus is associated with liver fibrosis progression in human immunodeficiency virus/hepatitis C virus-coinfected patients

There is growing evidence that human genetic variants contribute to liver fibrosis in subjects with hepatitis C virus (HCV) monoinfection, but this aspect has been little investigated in patients coinfected with HCV and human immunodeficiency virus (HIV). We performed the first genome-wide association study of liver fibrosis progression in patients coinfected with HCV and HIV, using the well-characterized French National Agency for Research on AIDS and Viral Hepatitis CO13 HEPAVIH cohort. Liver fibrosis was assessed by elastography (FibroScan), providing a quantitative fibrosis score. After quality control, a genome-wide association study was conducted on 289 Caucasian patients, for a total of 8,426,597 genotyped (Illumina Omni2.5 BeadChip) or reliably imputed single-nucleotide polymorphisms. Single-nucleotide polymorphisms with P values <10^-6 were investigated in two independent replication cohorts of European patients infected with HCV alone. Two signals of genome-wide significance (P < 5 × 10^-8 ) were obtained. The first, on chromosome 3p25 and corresponding to rs61183828 (P = 3.8 × 10^-9 ), was replicated in the two independent cohorts of patients with HCV monoinfection. The cluster of single-nucleotide polymorphisms in linkage disequilibrium with rs61183828 was located close to two genes involved in mechanisms affecting both cell signaling and cell structure (CAV3) or HCV replication (RAD18). The second signal, obtained with rs11790131 (P = 9.3 × 10^-9 ) on chromosome region 9p22, was not replicated.

CONCLUSION

This genome-wide association study identified a new locus associated with liver fibrosis severity in patients with HIV/HCV coinfection, on chromosome 3p25, a finding that was replicated in patients with HCV monoinfection; these results provide new relevant hypotheses for the pathogenesis of liver fibrosis in patients with HIV/HCV coinfection that may help define new targets for drug development or new prognostic tests, to improve patient care. (Hepatology 2016;64:1462-1472).

Quantitative proteomic analysis of the influence of lignin on biofuel production by Clostridium acetobutylicum ATCC 824

BACKGROUND

Clostridium acetobutylicum has been a focus of research because of its ability to produce high-value compounds that can be used as biofuels. Lignocellulose is a promising feedstock, but the lignin-cellulose-hemicellulose biomass complex requires chemical pre-treatment to yield fermentable saccharides, including cellulose-derived cellobiose, prior to bioproduction of acetone-butanol-ethanol (ABE) and hydrogen. Fermentation capability is limited by lignin and thus process optimization requires knowledge of lignin inhibition. The effects of lignin on cellular metabolism were evaluated for C. acetobutylicum grown on medium containing either cellobiose only or cellobiose plus lignin. Microscopy, gas chromatography and 8-plex iTRAQ-based quantitative proteomic technologies were applied to interrogate the effect of lignin on cellular morphology, fermentation and the proteome.

RESULTS

Our results demonstrate that C. acetobutylicum has reduced performance for solvent production when lignin is present in the medium. Medium supplemented with 1 g L(-1) of lignin led to delay and decreased solvents production (ethanol; 0.47 g L(-1) for cellobiose and 0.27 g L(-1) for cellobiose plus lignin and butanol; 0.13 g L(-1) for cellobiose and 0.04 g L(-1) for cellobiose plus lignin) at 20 and 48 h, respectively, resulting in the accumulation of acetic acid and butyric acid. Of 583 identified proteins (FDR < 1 %), 328 proteins were quantified with at least two unique peptides. Up- or down-regulation of protein expression was determined by comparison of exponential and stationary phases of cellobiose in the presence and absence of lignin. Of relevance, glycolysis and fermentative pathways were mostly down-regulated, during exponential and stationary growth phases in presence of lignin. Moreover, proteins involved in DNA repair, transcription/translation and GTP/ATP-dependent activities were also significantly affected and these changes were associated with altered cell morphology.

CONCLUSIONS

This is the first comprehensive analysis of the cellular responses of C. acetobutylicum to lignin at metabolic and physiological levels. These data will enable targeted metabolic engineering strategies to optimize biofuel production from biomass by overcoming limitations imposed by the presence of lignin.

Identification of Genes Whose Expression Profile Is Associated with Non-Progression towards AIDS Using eQTLs

BACKGROUND

Many genome-wide association studies have been performed on progression towards the acquired immune deficiency syndrome (AIDS) and they mainly identified associations within the HLA loci. In this study, we demonstrate that the integration of biological information, namely gene expression data, can enhance the sensitivity of genetic studies to unravel new genetic associations relevant to AIDS.

METHODS

We collated the biological information compiled from three databases of expression quantitative trait loci (eQTLs) involved in cells of the immune system. We derived a list of single nucleotide polymorphisms (SNPs) that are functional in that they correlate with differential expression of genes in at least two of the databases. We tested the association of those SNPs with AIDS progression in two cohorts, GRIV and ACS. Tests on permuted phenotypes of the GRIV and ACS cohorts or on randomised sets of equivalent SNPs allowed us to assess the statistical robustness of this method and to estimate the true positive rate.

RESULTS

Eight genes were identified with high confidence (p = 0.001, rate of true positives 75%). Some of those genes had previously been linked with HIV infection. Notably, ENTPD4 belongs to the same family as CD39, whose expression has already been associated with AIDS progression; while DNAJB12 is part of the HSP90 pathway, which is involved in the control of HIV latency. Our study also drew our attention to lesser-known functions such as mitochondrial ribosomal proteins and a zinc finger protein, ZFP57, which could be central to the effectiveness of HIV infection. Interestingly, for six out of those eight genes, down-regulation is associated with non-progression, which makes them appealing targets to develop drugs against HIV.

Reduced keratin expression in colorectal neoplasia and associated fields is reversible by diet and resection

BACKGROUND

Patients with adenomatous colonic polyps are at increased risk of developing further polyps suggesting field-wide alterations in cancer predisposition. The current study aimed to identify molecular alterations in the normal mucosa in the proximity of adenomatous polyps and to assess the modulating effect of butyrate, a chemopreventive compound produced by fermentation of dietary residues.

METHODS

A cross-sectional study was undertaken in patients with adenomatous polyps: biopsy samples were taken from the adenoma, and from macroscopically normal mucosa on the contralateral wall to the adenoma and from the mid-sigmoid colon. In normal subjects biopsies were taken from the mid-sigmoid colon. Biopsies were frozen for proteomic analysis or formalin-fixed for immunohistochemistry. Proteomic analysis was undertaken using iTRAQ workflows followed by bioinformatics analyses. A second dietary fibre intervention study arm used the same endpoints and sampling strategy at the beginning and end of a high-fibre intervention.

RESULTS

Key findings were that keratins 8, 18 and 19 were reduced in expression level with progressive proximity to the lesion. Lesional tissue exhibited multiple K8 immunoreactive bands and overall reduced levels of keratin. Biopsies from normal subjects with low faecal butyrate also showed depressed keratin expression. Resection of the lesion and elevation of dietary fibre intake both appeared to restore keratin expression level.

CONCLUSION

Changes in keratin expression associate with progression towards neoplasia, but remain modifiable risk factors. Dietary strategies may improve secondary chemoprevention.

TRIAL REGISTRATION NUMBER

ISRCTN90852168.

Evidence after imputation for a role of MICA variants in nonprogression and elite control of HIV type 1 infection

Past genome-wide association studies (GWAS) involving individuals with AIDS have mainly identified associations in the HLA region. Using the latest software, we imputed 7 million single-nucleotide polymorphisms (SNPs)/indels of the 1000 Genomes Project from the GWAS-determined genotypes of individuals in the Genomics of Resistance to Immunodeficiency Virus AIDS nonprogression cohort and compared them with those of control cohorts. The strongest signals were in MICA, the gene encoding major histocompatibility class I polypeptide-related sequence A (P = 3.31 × 10(-12)), with a particular exonic deletion (P = 1.59 × 10(-8)) in full linkage disequilibrium with the reference HCP5 rs2395029 SNP. Haplotype analysis also revealed an additive effect between HLA-C, HLA-B, and MICA variants. These data suggest a role for MICA in progression and elite control of human immunodeficiency virus type 1 infection.

Archaeal signal transduction: impact of protein phosphatase deletions on cell size, motility, and energy metabolism in Sulfolobus acidocaldarius

In this study, the in vitro and in vivo functions of the only two identified protein phosphatases, Saci-PTP and Saci-PP2A, in the crenarchaeal model organism Sulfolobus acidocaldarius were investigated. Biochemical characterization revealed that Saci-PTP is a dual-specific phosphatase (against pSer/pThr and pTyr), whereas Saci-PP2A exhibited specific pSer/pThr activity and inhibition by okadaic acid. Deletion of saci_pp2a resulted in pronounced alterations in growth, cell shape and cell size, which could be partially complemented. Transcriptome analysis of the three strains (Δsaci_ptp, Δsaci_pp2a and the MW001 parental strain) revealed 155 genes that were differentially expressed in the deletion mutants, and showed significant changes in expression of genes encoding the archaella (archaeal motility structure), components of the respiratory chain and transcriptional regulators. Phosphoproteome studies revealed 801 unique phosphoproteins in total, with an increase in identified phosphopeptides in the deletion mutants. Proteins from most functional categories were affected by phosphorylation, including components of the motility system, the respiratory chain, and regulatory proteins. In the saci_pp2a deletion mutant the up-regulation at the transcript level, as well as the observed phosphorylation pattern, resembled starvation stress responses. Hypermotility was also observed in the saci_pp2a deletion mutant. The results highlight the importance of protein phosphorylation in regulating essential cellular processes in the crenarchaeon S. acidocaldarius.

An insight into iTRAQ: where do we stand now?

The iTRAQ (isobaric tags for relative and absolute quantification) technique is widely employed in proteomic workflows requiring relative quantification. Here, we review the iTRAQ literature; in particular, we focus on iTRAQ usage in relation to other commonly used quantitative techniques e.g. stable isotope labelling in culture (SILAC), label-free methods and selected reaction monitoring (SRM). As a result, we identify several issues arising with respect to iTRAQ. Perhaps frustratingly, iTRAQ's attractiveness has been undermined by a number of technical and analytical limitations: it may not be truly quantitative, as the changes in abundance reported will generally be underestimated. We discuss weaknesses and strengths of iTRAQ as a methodology for relative quantification in the light of this and other technical issues. We focus on technical developments targeted at iTRAQ accuracy and precision, use of 4-plex over 8-plex reagents and application of iTRAQ to post-translational modification (PTM) workflows. We also discuss iTRAQ in relation to label-free approaches, to which iTRAQ is losing ground.

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.

A proteomic analysis of differential cellular responses to the short-chain fatty acids butyrate, valerate and propionate in colon epithelial cancer cells

The short chain fatty acids (SCFAs) are inhibitors of histone deacetylases (HDACi); they are produced naturally in the colon by fermentation. They affect cellular processes at a molecular and transcriptional level, the mechanisms of which may involve large numbers of proteins and integrated pathways. Butyrate is the most biologically potent of the SCFAs in colon epithelial cells, inhibiting human colon carcinoma cell proliferation and inducing apoptosis in vitro. In order to investigate the hypothesis that propionate and valerate possess unique and independent actions from butyrate, we combined proteomic and cellomic approaches for large-scale comparative analysis. Proteomic evaluation was undertaken using an iTRAQ tandem mass-spectrometry workflow and high-throughput High-content Analysis microscopy (HCA) was applied to generate cellomic information on the cell cycle and the cytoskeletal structure. Our results show that these SCFAs possess specific effects. Butyrate was shown to have more pronounced effects on the keratins and intermediate filaments (IFs); while valerate altered the β-tubulin isotypes' expression and the microtubules (MTs); propionate was involved in both mechanisms, displaying intermediate effects. These data suggest distinct physiological roles for SCFAs in colon epithelial function, offering new possibilities for cancer therapeutics.

Proteomic evaluation and validation of cathepsin D regulated proteins in macrophages exposed to Streptococcus pneumoniae

Macrophages are central effectors of innate immune responses to bacteria. We have investigated how activation of the abundant macrophage lysosomal protease, cathepsin D, regulates the macrophage proteome during killing of Streptococcus pneumoniae. Using the cathepsin D inhibitor pepstatin A, we demonstrate that cathepsin D differentially regulates multiple targets out of 679 proteins identified and quantified by eight-plex isobaric tag for relative and absolute quantitation. Our statistical analysis identified 18 differentially expressed proteins that passed all paired t-tests (α = 0.05). This dataset was enriched for proteins regulating the mitochondrial pathway of apoptosis or inhibiting competing death programs. Five proteins were selected for further analysis. Western blotting, followed by pharmacological inhibition or genetic manipulation of cathepsin D, verified cathepsin D-dependent regulation of these proteins, after exposure to S. pneumoniae. Superoxide dismutase-2 up-regulation was temporally related to increased reactive oxygen species generation. Gelsolin, a known regulator of mitochondrial outer membrane permeabilization, was down-regulated in association with cytochrome c release from mitochondria. Eukaryotic elongation factor (eEF2), a regulator of protein translation, was also down-regulated by cathepsin D. Using absence of the negative regulator of eEF2, eEF2 kinase, we confirm that eEF2 function is required to maintain expression of the anti-apoptotic protein Mcl-1, delaying macrophage apoptosis and confirm using a murine model that maintaining eEF2 function is associated with impaired macrophage apoptosis-associated killing of Streptococcus pneumoniae. These findings demonstrate that cathepsin D regulates multiple proteins controlling the mitochondrial pathway of macrophage apoptosis or competing death processes, facilitating intracellular bacterial killing.

Minimising iTRAQ ratio compression through understanding LC-MS elution dependence and high-resolution HILIC fractionation

Application of iTRAQ-based workflows for protein profiling has become widespread. Concomitantly, the idiosyncratic limitations of iTRAQ, such as its tendency to underestimate quantifications, have been studied and recognised. This report shows that the influence of ratio compression and limiting transmission in iTRAQ MS/MS in high-complexity mixtures (iTRAQ-labelled lysates) can be partly alleviated using high-resolution sample fractionation. Here, we also investigate in greater detail the dependency of iTRAQ quantification on the dynamics of online chromatography in low-complexity mixtures (iTRAQ-labelled standards). These findings will allow more efficient strategies to be designed for iTRAQ proteomics, alleviating iTRAQ underestimation and thus facilitating the detection of subtle abundance changes.

Quantitative protein expression and cell surface characteristics of Escherichia coli MG1655 biofilms

Cell surface physicochemical characterization techniques were combined with quantitative changes in protein expression, to investigate the biological and biophysical changes of Escherichia coli MG1655 cells when grown as a biofilm (BIO). The overall surface charge of BIO cells was found to be less negative, highlighting the need for a lower electrophoretic mobility for attachment to occur. Comparison of the chemical functional groups on the cell surface showed similar profiles, with the absorbance intensity higher for proteins and carbohydrates in the BIO cells. Quantitative proteomic analysis demonstrated that 3 proteins were significantly increased, and 9 proteins significantly decreased in abundance, in cells grown as a BIO compared to their planktonic counterparts, with 7 of these total 12 proteins unique to this study. Proteins showing significant increased or decreased abundance include proteins involved in acid resistance, DNA protection and binding and ABC transporters. Further predictive analysis of the metabolic pathways showed an increased abundance of the amino acid metabolism and tricarboxylic acid (TCA) cycle, with a decrease in expression within the pentose phosphate and glycolysis pathways. It is therefore hypothesized that cells grown as a BIO are still energetically viable potentially using amino acids as an indirect carbon backbone source into the TCA cycle.

Improving N-glycosylation efficiency in Escherichia coli using shotgun proteomics, metabolic network analysis, and selective reaction monitoring

Recently, the prospect of using Escherichia coli as a host for human glycoprotein production has increased due to detailed characterization of the prokaryotic N-glycosylation process and the ability to transfer the system into this bacterium. Although functionality of the native Campylobacter jejuni N-glycosylation system in E. coli has been demonstrated, the efficiency of the process using the well-characterized C. jejuni glycoprotein AcrA, was found to be low at 13.4±0.9% of total extracted protein. A combined approach using isobaric labeling of peptides and probability-based network analysis of metabolic changes was applied to forward engineer E. coli to improve glycosylation efficiency of AcrA. Enhancing flux through the glyoxylate cycle was identified as a potential metabolic manipulation to improve modification efficiency and was achieved by increasing the expression of isocitrate lyase. While the overall recombinant protein titre did not change significantly, the amount of glycosylated protein increased by approximately 300%.

Eight-plex iTRAQ analysis of variant metastatic human prostate cancer cells identifies candidate biomarkers of progression: An exploratory study

BACKGROUND

Due to the heterogeneity in the biological behavior of prostate cancer, biomarkers that can reliably distinguish indolent from aggressive disease are urgently needed to inform treatment choices.

METHODS

We employed 8-plex isobaric Tags for Relative and Absolute Quantitation (iTRAQ), to profile the proteomes of two distinct panels of isogenic prostate cancer cells with varying growth and metastatic potentials, in order to identify novel biomarkers associated with progression. The LNCaP, LNCaP-Pro5, and LNCaP-LN3 panel of cells represent a model of androgen-responsive prostate cancer, while the PC-3, PC-3M, and PC-3M-LN4 panel represent a model of androgen-insensitive disease.

RESULTS

Of the 245 unique proteins identified and quantified (>or=95% confidence; >or=2 peptides/protein), 17 showed significant differential expression (>or=+/-1.5), in at least one of the variant LNCaP cells relative to parental cells. Similarly, comparisons within the PC-3 panel identified 45 proteins to show significant differential expression in at least one of the variant PC-3 cells compared with parental cells. Differential expression of selected candidates was verified by Western blotting or immunocytochemistry, and corresponding mRNA expression was determined by quantitative real-time PCR (qRT-PCR). Immunostaining of prostate tissue microarrays for ERp5, one of the candidates identified, showed a significant higher immunoexpression in pre-malignant lesions compared with non-malignant epithelium (P < 0.0001, Mann-Whitney U-test), and in high Gleason grade (4-5) versus low grade (2-3) cancers (P < 0.05).

CONCLUSIONS

Our study provides proof of principle for the application of an 8-plex iTRAQ approach to uncover clinically relevant candidate biomarkers for prostate cancer progression.

A systems biology approach to investigate the response of Synechocystis sp. PCC6803 to a high salt environment

Background

Salt overloading during agricultural processes is causing a decrease in crop productivity due to saline sensitivity. Salt tolerant cyanobacteria share many cellular characteristics with higher plants and therefore make ideal model systems for studying salinity stress. Here, the response of fully adapted Synechocystis sp. PCC6803 cells to the addition of 6% w/v NaCl was investigated using proteomics combined with targeted analysis of transcripts.

Results

Isobaric mass tagging of peptides led to accurate relative quantitation and identification of 378 proteins, and approximately 40% of these were differentially expressed after incubation in BG-11 media supplemented with 6% salt for 9 days. Protein abundance changes were related to essential cellular functional alterations. Differentially expressed proteins involved in metabolic responses were also analysed using the probabilitistic tool Mixed Model on Graphs (MMG), where the role of energy conversion through glycolysis and reducing power through pentose phosphate pathway were highlighted. Temporal RT-qPCR experiments were also run to investigate protein expression changes at the transcript level, for 14 non-metabolic proteins. In 9 out of 14 cases the mRNA changes were in accordance with the proteins.

Conclusion

Synechocystis sp. PCC6803 has the ability to regulate essential metabolic processes to enable survival in high salt environments. This adaptation strategy is assisted by further regulation of proteins involved in non-metabolic cellular processes, supported by transcriptional and post-transcriptional control. This study demonstrates the effectiveness of using a systems biology approach in answering environmental, and in particular, salt adaptation questions in Synechocystis sp. PCC6803

Quantitative overview of N2 fixation in Nostoc punctiforme ATCC 29133 through cellular enrichments and iTRAQ shotgun proteomics

Nostoc punctiforme ATCC 29133 is a photoautotrophic cyanobacterium with the capacity to fix atmospheric N 2. Its ability to mediate this process is similar to that described for Nostoc sp. PCC 7120, where vegetative cells differentiate into heterocysts. Quantitative proteomic investigations at both the filament level and the heterocyst level are presented using isobaric tagging technology (iTRAQ), with 721 proteins at the 95% confidence interval quantified across both studies. Observations from both experiments yielded findings confirmatory of both transcriptional studies, and published Nostoc sp. PCC 7120 iTRAQ data. N. punctiforme exhibits similar metabolic trends, though changes in a number of metabolic pathways are less pronounced than in Nostoc sp. PCC 7120. Results also suggest a number of proteins that may benefit from future investigations. These include ATP dependent Zn-proteases, N-reserve degraders and also redox balance proteins. Complementary proteomic data sets from both organisms present key precursor knowledge that is important for future cyanobacterial biohydrogen research.

Neutral evolution of proteins: The superfunnel in sequence space and its relation to mutational robustness

Following Kimura's neutral theory of molecular evolution [M. Kimura, The Neutral Theory of Molecular Evolution (Cambridge University Press, Cambridge, 1983) (reprinted in 1986)], it has become common to assume that the vast majority of viable mutations of a gene confer little or no functional advantage. Yet, in silico models of protein evolution have shown that mutational robustness of sequences could be selected for, even in the context of neutral evolution. The evolution of a biological population can be seen as a diffusion on the network of viable sequences. This network is called a "neutral network." Depending on the mutation rate mu and the population size N, the biological population can evolve purely randomly (muN<<1) or it can evolve in such a way as to select for sequences of higher mutational robustness (muN>>1). The stringency of the selection depends not only on the product muN but also on the exact topology of the neutral network, the special arrangement of which was named "superfunnel." Even though the relation between mutation rate, population size, and selection was thoroughly investigated, a study of the salient topological features of the superfunnel that could affect the strength of the selection was wanting. This question is addressed in this study. We use two different models of proteins: on lattice and off lattice. We compare neutral networks computed using these models to random networks. From this, we identify two important factors of the topology that determine the stringency of the selection for mutationally robust sequences. First, the presence of highly connected nodes ("hubs") in the network increases the selection for mutationally robust sequences. Second, the stringency of the selection increases when the correlation between a sequence's mutational robustness and its neighbors' increases. The latter finding relates a global characteristic of the neutral network to a local one, which is attainable through experiments or molecular modeling.

Identifying differentially expressed subnetworks with MMG

BACKGROUND

Mixture model on graphs (MMG) is a probabilistic model that integrates network topology with (gene, protein) expression data to predict the regulation state of genes and proteins. It is remarkably robust to missing data, a feature particularly important for its use in quantitative proteomics. A new implementation in C and interfaced with R makes MMG extremely fast and easy to use and to extend.

AVAILABILITY

The original implementation (Matlab) is still available from http://www.dcs.shef.ac.uk/~guido/; the new implementation is available from http://wrightlab.group.shef.ac.uk/people_noirel.htm, from CRAN, and has been submitted to BioConductor, http://www.bioconductor.org/.

Automated extraction of meaningful pathways from quantitative proteomics data

Technological developments in the life sciences have resulted in an ever-accelerating pace of data production. Systems Biology tries to shed light upon these data by building complex models describing the interactions between biological components. However, extracting information from this morass of data requires the use of sophisticated computational techniques. Here, we propose a method suitable to integrate data drawn from quantitative proteomics into a metabolic scaffold and identify the metabolic pathways which are collectively up-regulated or down-regulated. The availability of such a tool is highly desirable as the extracted information could then be taken as a starting point for in-depth analyses, in particular in fields like Synthetic Biology, where datasets need be characterized routinely.

Neutral evolution of protein-protein interactions: a computational study using simple models

Background

Protein-protein interactions are central to cellular organization, and must have appeared at an early stage of evolution. To understand better their role, we consider a simple model of protein evolution and determine the effect of an explicit selection for Protein-protein interactions.

Results

In the model, viable sequences all have the same fitness, following the neutral evolution theory. A very simple, two-dimensional lattice representation of the protein structures is used, and the model only considers two kinds of amino acids: hydrophobic and polar. With these approximations, exact calculations are performed. The results do not depend too strongly on these assumptions, since a model using a 3D, off-lattice representation of the proteins gives results in qualitative agreement with the 2D one. With both models, the evolutionary dynamics lead to a steady state population that is enriched in sequences that dimerize with a high affinity, well beyond the minimal level needed to survive. Correspondingly, sequences close to the viability threshold are less abundant in the steady state, being subject to a larger proportion of lethal mutations. The set of viable sequences has a "funnel" shape, consistent with earlier studies: sequences that are highly populated in the steady state are "close" to each other (with proximity being measured by the number of amino acids that differ).

Conclusion

This bias in the the steady state sequences should lead to an increased resistance of the population to environmental change and an increased ability to evolve.