Proteogenomic mapping of Mycoplasma hyopneumoniae virulent strain 232

Background

Mycoplasma hyopneumoniae causes respiratory disease in swine and contributes to the porcine respiratory disease complex, a major disease problem in the swine industry. The M. hyopneumoniae strain 232 genome is one of the smallest and best annotated microbial genomes, containing only 728 annotated genes and 691 known proteins. Standard protein databases for mass spectrometry only allow for the identification of known and predicted proteins, which if incorrect can limit our understanding of the biological processes at work. Proteogenomic mapping is a methodology which allows the entire 6-frame genome translation of an organism to be used as a mass spectrometry database to help identify unknown proteins as well as correct and confirm existing annotations. This methodology will be employed to perform an in-depth analysis of the M. hyopneumoniae proteome.

Results

Proteomic analysis indicates 483 of 691 (70%) known M. hyopneumoniae strain 232 proteins are expressed under the culture conditions given in this study. Furthermore, 171 of 328 (52%) hypothetical proteins have been confirmed. Proteogenomic mapping resulted in the identification of previously unannotated genes gatC and rpmF and 5-prime extensions to genes mhp063, mhp073, and mhp451, all conserved and annotated in other M. hyopneumoniae strains and Mycoplasma species. Gene prediction with Prodigal, a prokaryotic gene predicting program, completely supports the new genomic coordinates calculated using proteogenomic mapping.

Conclusions

Proteogenomic mapping showed that the protein coding genes of the M. hyopneumoniae strain 232 identified in this study are well annotated. Only 1.8% of mapped peptides did not correspond to genes defined by the current genome annotation. This study also illustrates how proteogenomic mapping can be an important tool to help confirm, correct and append known gene models when using a genome sequence as search space for peptide mass spectra. Using a gene prediction program which scans for a wide variety of promoters can help ensure genes are accurately predicted or not missed completely. Furthermore, protein extraction using differential detergent fractionation effectively increases the number of membrane and cytoplasmic proteins identifiable my mass spectrometry.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-576) contains supplementary material, which is available to authorized users.     

Assignment of Putative Functions to Membrane “Hypothetical Proteins” from the <Emphasis Type="Italic">Trypanosoma cruzi</Emphasis> Genome

Protozoan parasites cause thousands of deaths each year in developing countries. The genome projects of these parasites opened a new era in the identification of therapeutic targets. However, the putative function could be predicted for fewer than half of the protein-coding genes. In this work, all Trypanosoma cruzi proteins containing predicted transmembrane spans were processed through an automated computational routine and further analyzed in order to assign the most probable function. The analysis consisted of dissecting the whole predicted protein in different regions. More than 5,000 sequences were processed, and the predicted biological functions were grouped into 19 categories according to the hits obtained after analysis. One focus of interest, due to the scarce information available on trypanosomatids, is the proteins involved in signal-transduction processes. In the present work, we identified 54 proteins belonging to this group, which were individually analyzed. The results show that by means of a simple pipeline it was possible to attribute probable functions to sequences annotated as coding for “hypothetical proteins.” Also, we successfully identified the majority of candidates participating in the signal-transduction pathways in T. cruzi.     

Experimental validation of in silico model‐predicted isocitrate dehydrogenase and phosphomannose isomerase from Dehalococcoides mccartyi

Gene sequences annotated as proteins of unknown or non‐specific function and hypothetical proteins account for a large fraction of most genomes. In the strictly anaerobic and organohalide respiring Dehalococcoides mccartyi, this lack of annotation plagues almost half the genome. Using a combination of bioinformatics analyses and genome‐wide metabolic modelling, new or more specific annotations were proposed for about 80 of these poorly annotated genes in previous investigations of D. mccartyi metabolism. Herein, we report the experimental validation of the proposed reannotations for two such genes (KB1_0495 and KB1_0553) from D. mccartyi strains in the KB‐1 community. KB1_0495 or DmIDH was originally annotated as an NAD⁺‐dependent isocitrate dehydrogenase, but biochemical assays revealed its activity primarily with NADP⁺ as a cofactor. KB1_0553, also denoted as DmPMI, was originally annotated as a hypothetical protein/sugar isomerase domain protein. We previously proposed that it was a bifunctional phosphoglucose isomerase/phosphomannose isomerase, but only phosphomannose isomerase activity was identified and confirmed experimentally. Further bioinformatics analyses of these two protein sequences suggest their affiliation to potentially novel enzyme families within their respective larger enzyme super families.     

Survey of Surface Proteins from the Pathogenic Mycoplasma hyopneumoniae Strain 7448 Using a Biotin Cell Surface Labeling Approach

The characterization of the repertoire of proteins exposed on the cell surface by Mycoplasma hyopneumoniae (M. hyopneumoniae), the etiological agent of enzootic pneumonia in pigs, is critical to understand physiological processes associated with bacterial infection capacity, survival and pathogenesis. Previous in silico studies predicted that about a third of the genes in the M. hyopneumoniae genome code for surface proteins, but so far, just a few of them have experimental confirmation of their expression and surface localization. In this work, M. hyopneumoniae surface proteins were labeled in intact cells with biotin, and affinity-captured biotin-labeled proteins were identified by a gel-based liquid chromatography-tandem mass spectrometry approach. A total of 20 gel slices were separately analyzed by mass spectrometry, resulting in 165 protein identifications corresponding to 59 different protein species. The identified surface exposed proteins better defined the set of M. hyopneumoniae proteins exposed to the host and added confidence to in silico predictions. Several proteins potentially related to pathogenesis, were identified, including known adhesins and also hypothetical proteins with adhesin-like topologies, consisting of a transmembrane helix and a large tail exposed at the cell surface. The results provided a better picture of the M. hyopneumoniae cell surface that will help in the understanding of processes important for bacterial pathogenesis. Considering the experimental demonstration of surface exposure, adhesion-like topology predictions and absence of orthologs in the closely related, non-pathogenic species Mycoplasma flocculare, several proteins could be proposed as potential targets for the development of drugs, vaccines and/or immunodiagnostic tests for enzootic pneumonia.     

Structural Genomics of Minimal Organisms and Protein Fold Space

The initial aim of the Berkeley Structural Genomics Center is to obtain a near-complete structural complement of two minimal organisms, closely related pathogens Mycoplasma genitalium and M. pneumoniae. The former has fewer than 500 genes and the latter fewer than 700 genes. To achieve this goal, the current protein targets have been selected starting with those predicted to be most tractable and likely to yield new structural and functional information. During the past 3 years, the semi-automated structural genomics pipeline has been set up from cloning, expression, purification, and ultimately to structural determination. The results from the pipeline substantially increased the coverage of the protein fold space of M. pneumoniae and M. genitalium. Furthermore, about 1/2 of the structures of ‘unique’ protein sequences revealed new and novel folds, and over 2/3 of the structures of previously annotated ‘hypothetical proteins’ inferred their molecular functions.     

Hydrogen peroxide production and myo‐inositol metabolism as important traits for virulence of Mycoplasma hyopneumoniae

Mycoplasma hyopneumoniae is the causative agent of enzootic pneumonia. In our previous work, we reconstructed the metabolic models of this species along with two other mycoplasmas from the respiratory tract of swine: Mycoplasma hyorhinis, considered less pathogenic but which nonetheless causes disease and Mycoplasma flocculare, a commensal bacterium. We identified metabolic differences that partially explained their different levels of pathogenicity. One important trait was the production of hydrogen peroxide from the glycerol metabolism only in the pathogenic species. Another important feature was a pathway for the metabolism of myo‐inositol in M. hyopneumoniae. Here, we tested these traits to understand their relation to the different levels of pathogenicity, comparing not only the species but also pathogenic and attenuated strains of M. hyopneumoniae. Regarding the myo‐inositol metabolism, we show that only M. hyopneumoniae assimilated this carbohydrate and remained viable when myo‐inositol was the primary energy source. Strikingly, only the two pathogenic strains of M. hyopneumoniae produced hydrogen peroxide in complex medium. We also show that this production was dependent on the presence of glycerol. Although further functional tests are needed, we present in this work two interesting metabolic traits of M. hyopneumoniae that might be directly related to its enhanced virulence.     

A comprehensive analysis of Trypanosoma brucei mitochondrial proteome

The composition of the large, single, mitochondrion (mt) of Trypanosoma brucei was characterized by MS (2‐D LC‐MS/MS and gel‐LC‐MS/MS) analyses. A total of 2897 proteins representing a substantial proportion of procyclic form cellular proteome were identified, which confirmed the validity of the vast majority of gene predictions. The data also showed that the genes annotated as hypothetical (species specific) were overpredicted and that virtually all genes annotated as hypothetical, unlikely are not expressed. By comparing the MS data with genome sequence, 40 genes were identified that were not previously predicted. The data are placed in a publicly available web‐based database (www.TrypsProteome.org). The total mitochondrial proteome is estimated at 1008 proteins, with 401, 196, and 283 assigned to the mt with high, moderate, and lower confidence, respectively. The remaining mitochondrial proteins were estimated by statistical methods although individual assignments could not be made. The identified proteins have predicted roles in macromolecular, metabolic, energy generating, and transport processes providing a comprehensive profile of the protein content and function of the T. brucei mt.     

Analysing the outer membrane subproteome of Methylococcus capsulatus (Bath) using proteomics and novel biocomputing tools

High-resolution two-dimensional gel electrophoresis and mass spectrometry has been used to identify the outer membrane (OM) subproteome of the Gram-negative bacterium Methylococcus capsulatus (Bath). Twenty-eight unique polypeptide sequences were identified from protein samples enriched in OMs. Only six of these polypeptides had previously been identified. The predictions from novel bioinformatic methods predicting β-barrel outer membrane proteins (OMPs) and OM lipoproteins were compared to proteins identified experimentally. BOMP () predicted 43 β-barrel OMPs (1.45%) from the 2,959 annotated open reading frames. This was a lower percentage than predicted from other Gram-negative proteomes (1.8–3%). More than half of the predicted BOMPs in M. capsulatus were annotated as (conserved) hypothetical proteins with significant similarity to very few sequences in Swiss-Prot or TrEMBL. The experimental data and the computer predictions indicated that the protein composition of the M. capsulatus OM subproteome was different from that of other Gram-negative bacteria studied in a similar manner. A new program, Lipo, was developed that can analyse entire predicted proteomes and give a list of recognised lipoproteins categorised according to their lipo-box similarity to known Gram-negative lipoproteins (). This report is the first using a proteomics and bioinformatics approach to identify the OM subproteome of an obligate methanotroph.     

Metabolic modeling of energy balances in Mycoplasma hyopneumoniae shows that pyruvate addition increases growth rate

Mycoplasma hyopneumoniae is cultured on large‐scale to produce antigen for inactivated whole‐cell vaccines against respiratory disease in pigs. However, the fastidious nutrient requirements of this minimal bacterium and the low growth rate make it challenging to reach sufficient biomass yield for antigen production. In this study, we sequenced the genome of M. hyopneumoniae strain 11 and constructed a high quality constraint‐based genome‐scale metabolic model of 284 chemical reactions and 298 metabolites. We validated the model with time‐series data of duplicate fermentation cultures to aim for an integrated model describing the dynamic profiles measured in fermentations. The model predicted that 84% of cellular energy in a standard M. hyopneumoniae cultivation was used for non‐growth associated maintenance and only 16% of cellular energy was used for growth and growth associated maintenance. Following a cycle of model‐driven experimentation in dedicated fermentation experiments, we were able to increase the fraction of cellular energy used for growth through pyruvate addition to the medium. This increase in turn led to an increase in growth rate and a 2.3 times increase in the total biomass concentration reached after 3–4 days of fermentation, enhancing the productivity of the overall process. The model presented provides a solid basis to understand and further improve M. hyopneumoniae fermentation processes. Biotechnol. Bioeng. 2017;114: 2339–2347. © 2017 The Authors. Biotechnology and Bioengineering published by Wiley Periodicals, Inc.     

The genome sequence of Mycoplasma hyopneumoniae strain 232, the agent of swine mycoplasmosis

We present the complete genome sequence of Mycoplasma hyopneumoniae, an important member of the porcine respiratory disease complex. The genome is composed of 892,758 bp and has an average G+C content of 28.6 mol%. There are 692 predicted protein coding sequences, the average protein size is 388 amino acids, and the mean coding density is 91%. Functions have been assigned to 304 (44%) of the predicted protein coding sequences, while 261 (38%) of the proteins are conserved hypothetical proteins and 127 (18%) are unique hypothetical proteins. There is a single 16S-23S rRNA operon, and there are 30 tRNA coding sequences. The cilium adhesin gene has six paralogs in the genome, only one of which contains the cilium binding site. The companion gene, P102, also has six paralogs. Gene families constitute 26.3% of the total coding sequences, and the largest family is the 34-member ABC transporter family. Protein secretion occurs through a truncated pathway consisting of SecA, SecY, SecD, PrsA, DnaK, Tig, and LepA. Some highly conserved eubacterial proteins, such as GroEL and GroES, are notably absent. The DnaK-DnaJ-GrpR complex is intact, providing the only control over protein folding. There are several proteases that might serve as virulence factors, and there are 53 coding sequences with prokaryotic lipoprotein lipid attachment sites. Unlike other mycoplasmas, M. hyopneumoniae contains few genes with tandem repeat sequences that could be involved in phase switching or antigenic variation. Thus, it is not clear how M. hyopneumoniae evades the immune response and establishes a chronic infection.     

Post-translational processing targets functionally diverse proteins in Mycoplasma hyopneumoniae

Mycoplasma hyopneumoniae is a genome-reduced, cell wall-less, bacterial pathogen with a predicted coding capacity of less than 700 proteins and is one of the smallest self-replicating pathogens. The cell surface of M. hyopneumoniae is extensively modified by processing events that target the P97 and P102 adhesin families. Here, we present analyses of the proteome of M. hyopneumoniae-type strain J using protein-centric approaches (one- and two-dimensional GeLC–MS/MS) that enabled us to focus on global processing events in this species. While these approaches only identified 52% of the predicted proteome (347 proteins), our analyses identified 35 surface-associated proteins with widely divergent functions that were targets of unusual endoproteolytic processing events, including cell adhesins, lipoproteins and proteins with canonical functions in the cytosol that moonlight on the cell surface. Affinity chromatography assays that separately used heparin, fibronectin, actin and host epithelial cell surface proteins as bait recovered cleavage products derived from these processed proteins, suggesting these fragments interact directly with the bait proteins and display previously unrecognized adhesive functions. We hypothesize that protein processing is underestimated as a post-translational modification in genome-reduced bacteria and prokaryotes more broadly, and represents an important mechanism for creating cell surface protein diversity.     

Identification of coding sequences from a freshly prepared Trypanosoma brucei brucei expression library by polymerase chain reaction

Animal African trypanosomiasis (AAT) also known as Nagana is a devastating disease among domestic animals in large parts of Sub-Saharan Africa causing loses in milk and meat production as well as traction power. However, there is currently no commercial vaccine against AAT. The parasites have also developed resistance to some of the drugs in use. Moreover, the use of affordable computer-aided wet bench methods in the search for vaccine and/or new drug targets against this disease have not yet been fully explored in developing countries. This study, therefore, explored the use of PCR to screen a freshly prepared bloodstream form Trypanosoma brucei brucei (T. b. brucei) expression library for coding sequences followed by bioinformatics analyses specifying the functions and importance of these proteins to parasite survival. Eleven protein coding sequences were identified from twenty nine purified clones. The putative retro transposon hot spot protein 4 (RHSP 4) was the only protein with a fully annotated DNA sequence. All the others were hypothetical or had partial or unqualified annotations. RHSP 4 and pyruvate dehydrogenase E1 component, alpha sub-unit (PDE1α) are involved in aerobic respiration whereas succinyl-Co A-3-ketoacid-coenzyme A transferase mitochondrial precursor (SKTMP) is predicted to be involved in ketone body catabolism. Cystathionine beta-synthase (CBS) and alpha-1,3-mannosyltransferase (αMT) have been predicted in cysteine biosynthesis and vesicular transport respectively. The functions of the hypothetical proteins encountered have neither been experimentally determined nor predicted. We hypothesize that both CBS and PDE1α are good drug targets. Overall, about 300 plates are required to PCR screen the entire Trypanosoma brucei genome in approximately eight months. This method is therefore, applicable and affordable in the search for new drug targets under conditions of limited resources among developing countries.     

Updating the proteome of the uncultivable hemotrophic Mycoplasma suis in experimentally infected pigs

Mycoplasma suis belongs to the hemotrophic mycoplasmas that are associated with acute and chronic anemia in a wide range of livestock and wild animals. The inability to culture M. suis in vitro has hindered its characterization at the molecular level. Since the publication of M. suis genome sequences in 2011 only one proteome study has been published. Aim of the presented study was to significantly extend the proteome coverage of M. suis strain KI_3806 during acute infection by applying three different protein extraction methods followed by 1D SDS‐PAGE and LC‐MS/MS. A total of 404 of 795 M. suis KI_3806 proteins (50.8%) were identified. Data analysis revealed the expression of 83.7% of the predicted ORFs with assigned functions but also highlights the expression of 179 of 523 (34.2%) hypothetical proteins with unknown functions. Computational analyses identified expressed membrane‐associated hypothetical proteins that might be involved in adhesion or host–pathogen interaction. Furthermore, analyses of the expressed proteins indicated the existence of a hexose‐6‐phosphate‐transporter and an ECF transporter. In conclusion, our proteome study provides a further step toward the elucidation of the unique life cycle of M. suis and the establishment of an in vitro culture. All MS data have been deposited in the ProteomeXchange with identifier PXD002294 ( http://proteomecentral.proteomexchange.org/dataset/PXD002294 ).     

Analysis of two large functionally uncharacterized regions in the <Emphasis Type="Italic">Methanopyrus kandleri</Emphasis> AV19 genome

Background

For most sequenced prokaryotic genomes, about a third of the protein coding genes annotated are "orphan proteins", that is, they lack homology to known proteins. These hypothetical genes are typically short and randomly scattered throughout the genome. This trend is seen for most of the bacterial and archaeal genomes published to date.

Results

In contrast we have found that a large fraction of the genes coding for such orphan proteins in the Methanopyrus kandleri AV19 genome occur within two large regions. These genes have no known homologs except from other M. kandleri genes. However, analysis of their lengths, codon usage, and Ribosomal Binding Site (RBS) sequences shows that they are most likely true protein coding genes and not random open reading frames.

Conclusions

Although these regions can be considered as candidates for massive lateral gene transfer, our bioinformatics analysis suggests that this is not the case. We predict many of the organism specific proteins to be transmembrane and belong to protein families that are non-randomly distributed between the regions. Consistent with this, we suggest that the two regions are most likely unrelated, and that they may be integrated plasmids.

     

Swine and poultry pathogens: the complete genome sequences of two strains of Mycoplasma hyopneumoniae and a strain of Mycoplasma synoviae

This work reports the results of analyses of three complete mycoplasma genomes, a pathogenic (7448) and a nonpathogenic (J) strain of the swine pathogen Mycoplasma hyopneumoniae and a strain of the avian pathogen Mycoplasma synoviae; the genome sizes of the three strains were 920,079 bp, 897,405 bp, and 799,476 bp, respectively. These genomes were compared with other sequenced mycoplasma genomes reported in the literature to examine several aspects of mycoplasma evolution. Strain-specific regions, including integrative and conjugal elements, and genome rearrangements and alterations in adhesin sequences were observed in the M. hyopneumoniae strains, and all of these were potentially related to pathogenicity. Genomic comparisons revealed that reduction in genome size implied loss of redundant metabolic pathways, with maintenance of alternative routes in different species. Horizontal gene transfer was consistently observed between M. synoviae and Mycoplasma gallisepticum. Our analyses indicated a likely transfer event of hemagglutinin-coding DNA sequences from M. gallisepticum to M. synoviae.     

Gene content and density in banana (<Emphasis Type="Italic">Musa acuminata</Emphasis>) as revealed by genomic sequencing of BAC clones

The complete sequence of Musa acuminata bacterial artificial chromosome (BAC) clones is presented and, consequently, the first analysis of the banana genome organization. One clone (MuH9) is 82,723 bp long with an overall G+C content of 38.2%. Twelve putative protein-coding sequences were identified, representing a gene density of one per 6.9 kb, which is slightly less than that previously reported for Arabidopsis but similar to rice. One coding sequence was identified as a partial M. acuminata malate synthase, while the remaining sequences showed a similarity to predicted or hypothetical proteins identified in genome sequence data. A second BAC clone (MuG9) is 73,268 bp long with an overall G+C content of 38.5%. Only seven putative coding regions were discovered, representing a gene density of only one gene per 10.5 kb, which is strikingly lower than that of the first BAC. One coding sequence showed significant homology to the soybean ribonucleotide reductase (large subunit). A transition point between coding regions and repeated sequences was found at approximately 45 kb, separating the coding upstream BAC end from its downstream end that mainly contained transposon-like sequences and regions similar to known repetitive sequences of M. acuminata. This gene organization resembles Gramineae genome sequences, where genes are clustered in gene-rich regions separated by gene-poor DNA containing abundant transposons.Communicated by J.S. Heslop-Harrison     

Structural Bioinformatics of <Emphasis Type="Italic">Neisseria meningitidis</Emphasis> LD-Carboxypeptidase: Implications for Substrate Binding and Specificity

Neisseria meningitidis, a gram negative bacterium, is the leading cause of bacterial meningitis and severe sepsis. Neisseria meningitidis genome contains 2,160 predicted coding regions including 1,000 hypothetical genes. Re-annotation of N. meningitidis hypothetical proteins identified nine putative peptidases. Among them, the NMB1620 protein was annotated as LD-carboxypeptidase involved in peptidoglycan recycling. Structural bioinformatics studies of NMB1620 protein using homology modeling and ligand docking were carried out. Structural comparison of substrate binding site of LD-carboxypeptidase was performed based on binding of tetrapeptide substrate ‘l-alanyl-d-glutamyl-meso-diaminopimelyl-d-alanine’. Inspection of different subsite-forming residues showed changeability in the S1 subsite across different bacterial species. This variability was predicted to provide a structural basis to S1-subsite for accommodating different amino acid residues at P1 position of the tetrapeptide substrate ‘l-alanyl-d-glutamyl-meso-diaminopimelyl-d-alanine’.     

Sequence of Leptospira santarosai serovar Shermani genome and prediction of virulence-associated genes

Leptospirosis, a widespread zoonosis, is a re-emerging infectious disease caused by pathogenic Leptospira species. In Taiwan, Leptospira santarosai serovar Shermani is the most frequently isolated serovar, causing both renal and systemic infections. This study aimed to generate a L. santarosai serovar Shermani genome sequence and categorize its hypothetical genes, particularly those associated with virulence. The genome sequence consists of 3,936,333 nucleotides and 4033 predicted genes. Additionally, 2244 coding sequences could be placed into clusters of orthologous groups and the number of genes involving cell wall/membrane/envelope biogenesis and defense mechanisms was higher than that of other Leptospira spp. Comparative genetic analysis based on BLASTX data revealed that about 73% and 68.8% of all coding sequences have matches to pathogenic L. interrogans and L. borgpetersenii, respectively, and about 57.6% to saprophyte L. biflexa. Among the hypothetical proteins, 421 have a transmembrane region, 172 have a signal peptide and 17 possess a lipoprotein signature. According to PFAM prediction, 32 hypothetical proteins have properties of toxins and surface proteins mediated bacterial attachment, suggesting they may have roles associated with virulence. The availability of the genome sequence of L. santarosai serovar Shermani and the bioinformatics re-annotation of leptospiral hypothetical proteins will facilitate further functional genomic studies to elucidate the pathogenesis of leptospirosis and develop leptospiral vaccines.