Protein abundance ratios for global studies of prokaryotes

The use of multidimensional capillary HPLC combined with MS/MS has allowed high qualitative and quantitative proteome coverage of prokaryotic organisms. The determination of protein abundance change between two or more conditions has matured to the point that false discovery rates can be very low and for smaller proteomes coverage is sufficiently high to explicitly consider false negative error. Selected aspects of using these methods for global protein abundance assessments are reviewed. These include instrumental issues that influence the reliability of abundance ratios; a comparison of sources of nonlinearity, errors, and data compression in proteomics and spotted cDNA arrays; strengths and weaknesses of spectral counting versus stable isotope metabolic labeling; and a survey of microbiological applications of global abundance analysis at the protein level. Proteomic results for two organisms that have been studied extensively using these methods are reviewed in greater detail. Spectral counting and metabolic labeling data are compared and the utility of proteomics for global gene regulation studies are discussed for the methanogenic Archaeon Methanococcus maripaludis. The oral pathogen Porphyromonas gingivalis is discussed as an example of an organism where a large percentage of the proteome differs in relative abundance between the intracellular and extracellular phenotype.     

Analysis of gene expression profiles associated with cisplatin resistance in human ovarian cancer cell lines and tissues using cDNA microarray

Gene expression profiles were analyzed by using cDNA microarray for a cisplatin-sensitive cell line (KF), and three- and thirty-fold cisplatin-resistant ovarian cancer cell lines (KFr and KFrP200) both showing no p53 mutation within exon 5, 6, 7, 8 and no pglycoprotein overexpression. Expression of GST-pi mRNA increased as the level of resistance to cisplatin became high. Microarray analysis revealed that DNA repair associated genes, i.e., XRCC5, XRCC6, ERCC5, hMLH1 were over-expressed in three-fold cisplatin-resistant cell line, KFr as compared to cisplatin-sensitive parental cell line, KF. Apoptosis inhibitors, i.e., IGFR type I and II were over-expressed, and apoptosis inducer, i.e., caspase 3 and BAK were underexpressed in highly cisplatin-resistant cell line, KFrP200 as compared to KFr. As for clinical cases, cDNA microarray was used to compare gene expression profiles directly between two groups, i.e., the chemotherapy (CAP) sensitive group (n = 2) and the resistant group (n = 2). Six genes such as beta tubulin, high-mobility group (nonhistone chromosomal) protein 1, connective tissue growth factor, insulin-like growth factor binding protein 2, alpha tubulin, and RAS-related gene were overexpressed in CAP therapy resistance group, whereas seven genes such as CD9 antigen, alpha-2-macroglobulin, caveolin 2, interleukin 1 receptor antagonist, Rho GTPase activating protein 1, reticulon 3, cyclin-dependent kinase 10, keratin 7 were underexpressed in CAP therapy resistance group. By increasing clinical case number and gene number of microarray to be used in the analysis of expression profile of gene cluster affecting anticancer drug resistance and sensitivity of the ovarian cancer, it would be possible to apply microarray analysis to personalization of chemotherapy such as selection of effective chemotherapy protocol and prediction of therapeutic effect in the near future.     

Differential responses of Coffea arabica L. leaves and roots to chemically induced systemic acquired resistance.

Coffea arabica is susceptible to several pests and diseases, some of which affect the leaves and roots. Systemic acquired resistance (SAR) is the main defence mechanism activated in plants in response to pathogen attack. Here, we report the effects of benzo(1,2,3)thiadiazole-7-carbothioic acid-S-methyl ester (BTH), a SAR chemical inducer, on the expression profile of C. arabica. Two cDNA libraries were constructed from the mRNA isolated from leaves and embryonic roots to create 1587 nonredundant expressed sequence tags (ESTs). We developed a cDNA microarray containing 1506 ESTs from the leaves and embryonic roots, and 48 NBS-LRR (nucleotide-binding site leucine-rich repeat) gene fragments derived from 2 specific genomic libraries. Competitive hybridization between untreated and BTH-treated leaves resulted in 55 genes that were significantly overexpressed and 16 genes that were significantly underexpressed. In the roots, 37 and 42 genes were over and underexpressed, respectively. A general shift in metabolism from housekeeping to defence occurred in the leaves and roots after BTH treatment. We observed a systemic increase in pathogenesis-related protein synthesis, in the oxidative burst, and in the cell wall strengthening processes. Moreover, responses in the roots and leaves varied significantly.     

Expression profiling by oligonucleotide microarrays spotted on coated polymer slides

We have developed a ready-to-spot polymer microarray slide, which is coated with a uniform layer of reactive electrophilic groups using anthraquinone-mediated photo-coupling chemistry. The slide coating reduces the hydrophobicity of the native polymer significantly, thereby enabling robust and efficient one-step coupling of spotted 5' amino-linked oligonucleotides onto the polymer slide. The utility of the coated polymer slide in gene expression profiling was assessed by fabrication of spotted oligonucleotide microarrays using a collection of 5' amino-linked 70-mer oligonucleotide probes representing 96 yeast genes from Operon. Two-colour hybridizations with labelled cDNA target pools derived from standard grown and heat-shocked wild type yeast cells could reproducibly measure heat shock induced expression of seven different heat shock protein (HSP) genes. Moreover, the observed fold changes were comparable to those reported previously using spotted cDNA arrays and high-density 25-mer oligonucleotide arrays from Affymetrix. The low hybridization signals obtained from the DeltaSSA4 mutant cDNA target, together with the high signal detected in two-colour hybridizations with heat-shocked wild type yeast relative to the DeltaSSA4 mutant strain implies that unspecific binding of cDNA target to the SSA4-specific 70-mer oligonucleotide probes is negligible. Combined, our results indicate that the coated polymer microarray slide represents a robust and cost-effective array platform for pre-spotted oligonucleotide arrays.     

Current issues for DNA microarrays: platform comparison, double linear amplification, and universal RNA reference

DNA microarray technology has been widely used to simultaneously determine the expression levels of thousands of genes. A variety of approaches have been used, both in the implementation of this technology and in the analysis of the large amount of expression data. However, several practical issues still have not been resolved in a satisfactory manner, and among the most critical is the lack of agreement in the results obtained in different array platforms. In this study, we present a comparison of several microarray platforms [Affymetrix oligonucleotide arrays, custom complementary DNA (cDNA) arrays, and custom oligo arrays printed with oligonucleotides from three different sources] as well as analysis of various methods used for microarray target preparation and the reference design. The results indicate that the pairwise correlations of expression levels between platforms are relative low overall but that the log ratios of the highly expressed genes are strongly correlated, especially between Affymetrix and cDNA arrays. The microarray measurements were compared with quantitative real-time-polymerase chain reaction (QRT-PCR) results for 23 genes, and the varying degrees of agreement for each platform were characterized. We have also developed and tested a double amplification method which allows the use of smaller amounts of starting material. The added round of amplification produced reproducible results as compared to the arrays hybridized with single round amplified targets. Finally, the reliability of using a universal RNA reference for two-channel microarrays was tested and the results suggest that comparisons of multiple experimental conditions using the same control can be accurate.     

tsORFdb: Theoretical Small Open Reading Frames (ORFs) database and massProphet: Peptide Mass Fingerprinting (PMF) tool for unknown small functional ORFs

Peptide mass fingerprinting (PMF) has become one of the most widely used methods for rapid identification of proteins in proteomics research. Many peaks, however, remain unassigned after PMF analysis, partly because of post-translational modification and the limited scope of protein sequences. Almost all PMF tools employ only known or predicted protein sequences and do not include open reading frames (ORFs) in the genome, which eliminates the chance of finding novel functional peptides. Unlike most tools that search protein sequences from known coding sequences, the tool we developed uses a database for theoretical small ORFs (tsORFs) and a PMF application using a tsORFs database (tsORFdb). The tsORFdb is a database for ORFeome that encompasses all potential tsORFs derived from whole genome sequences as well as the predicted ones. The massProphet system tries to extend the search scope to include the ORFeome using the tsORFdb. The tsORFdb and massProphet should be useful for proteomics research to give information about unknown small ORFs as well as predicted and registered proteins.     

Quantitative proteomics of intracellular Porphyromonas gingivalis

Whole-cell quantitative proteomic analyses were conducted to investigate the change from an extracellular to intracellular lifestyle for Porphyromonas gingivalis, a Gram-negative intracellular pathogen associated with periodontal disease. Global protein abundance data for P. gingivalis strain ATCC 33277 internalized for 18 h within human gingival epithelial cells and controls exposed to gingival cell culture medium were obtained at sufficient coverage to provide strong evidence that these changes are profound. A total of 385 proteins were overexpressed in internalized P. gingivalis relative to controls; 240 proteins were shown to be underexpressed. This represented in total about 28% of the protein encoding ORFs annotated for this organism, and slightly less than half of the proteins that were observed experimentally. Production of several proteases, including the classical virulence factors RgpA, RgpB, and Kgp, was decreased. A separate validation study was carried out in which a 16-fold dilution of the P. gingivalis proteome was compared to the undiluted sample in order to assess the quantitative false negative rate (all ratios truly alternative). Truly null (no change) abundance ratios from technical replicates were used to assess the rate of quantitative false positives over the entire proteome. A global comparison between the direction of abundance change observed and previously published bioinformatic gene pair predictions for P. gingivalis will assist with future studies of P. gingivalis gene regulation and operon prediction.     

Proteogenomic mapping as a complementary method to perform genome annotation

The accelerated rate of genomic sequencing has led to an abundance of completely sequenced genomes. Annotation of the open reading frames (ORFs) (i.e., gene prediction) in these genomes is an important task and is most often performed computationally based on features in the nucleic acid sequence. Using recent advances in proteomics, we set out to predict the set of ORFs for an organism based principally on expressed protein-based evidence. Using a novel search strategy, we mapped peptides detected in a whole-cell lysate of Mycoplasma pneumoniae onto a genomic scaffold and extended these "hits" into ORFs bound by traditional genetic signals to generate a "proteogenomic map". We were able to generate an ORF model for M. pneumoniae strain FH using proteomic data with a high correlation to models based on sequence features. Ultimately, we detected over 81% of the genomically predicted ORFs in M. pneumoniae strain M129 (the originally sequenced strain). We were also able to detect several new ORFs not originally predicted by genomic methods, various N-terminal extensions, and some evidence that would suggest that certain predicted ORFs are bogus. Some of these differences may be a result of the strain analyzed but demonstrate the robustness of protein analysis across closely related genomes. This technique is a cost-effective means to add value to genome annotation, and a prerequisite for proteome quantitation and in vivo interaction measures.     

Complete genome sequence of the genetically tractable hydrogenotrophic methanogen Methanococcus maripaludis

The genome sequence of the genetically tractable, mesophilic, hydrogenotrophic methanogen Methanococcus maripaludis contains 1,722 protein-coding genes in a single circular chromosome of 1,661,137 bp. Of the protein-coding genes (open reading frames [ORFs]), 44% were assigned a function, 48% were conserved but had unknown or uncertain functions, and 7.5% (129 ORFs) were unique to M. maripaludis. Of the unique ORFs, 27 were confirmed to encode proteins by the mass spectrometric identification of unique peptides. Genes for most known functions and pathways were identified. For example, a full complement of hydrogenases and methanogenesis enzymes was identified, including eight selenocysteine-containing proteins, with each being paralogous to a cysteine-containing counterpart. At least 59 proteins were predicted to contain iron-sulfur centers, including ferredoxins, polyferredoxins, and subunits of enzymes with various redox functions. Unusual features included the absence of a Cdc6 homolog, implying a variation in replication initiation, and the presence of a bacterial-like RNase HI as well as an RNase HII typical of the Archaea. The presence of alanine dehydrogenase and alanine racemase, which are uniquely present among the Archaea, explained the ability of the organism to use L- and D-alanine as nitrogen sources. Features that contrasted with the related organism Methanocaldococcus jannaschii included the absence of inteins, even though close homologs of most intein-containing proteins were encoded. Although two-thirds of the ORFs had their highest Blastp hits in Methanocaldococcus jannaschii, lateral gene transfer or gene loss has apparently resulted in genes, which are often clustered, with top Blastp hits in more distantly related groups.     

Affinity proteomics for systematic protein profiling of chromosome 21 gene products in human tissues

Here we show that an affinity proteomics strategy using affinity-purified antibodies raised against recombinant human protein fragments can be used for chromosome-wide protein profiling. The approach is based on affinity reagents raised toward bioinformatics-designed protein epitope signature tags corresponding to unique regions of individual gene loci. The genes of human chromosome 21 identified by the genome efforts were investigated, and the success rates for de novo cloning, protein production, and antibody generation were 85, 76, and 56%, respectively. Using human tissue arrays, a systematic profiling of protein expression and subcellular localization was undertaken for the putative gene products. The results suggest that this affinity proteomics strategy can be used to produce a proteome atlas, describing distribution and expression of proteins in normal tissues as well as in common cancers and other forms of diseased tissues.     

Correlation of protein expression, Gleason score and DNA ploidy in prostate cancer

The prognosis of prostate cancer correlates with tumor differentiation. Gleason score and DNA ploidy are two prognostic factors that correlate with prognosis. We analyzed differences in protein expression in prostate cancer of high and low aggressiveness according to these measures. From 35 prostatectomy specimens, 29 cancer samples and 10 benign samples were harvested by scraping cells from cut surfaces. DNA ploidy was assessed by image cytometry. Protein preparations from cell suspensions were examined by 2-DE. Protein spots that differed quantitatively between sample groups were identified by MS fingerprinting of tryptic fragments and MS/MS sequence analysis. We found 39 protein spots with expression levels that were raised or lowered in correlation with Gleason score and/or DNA ploidy pattern (31 overexpressed in high-malignant cancer, 8 underexpressed). Of these, 30 were identified by MS. Among overexpressed proteins were heat-shock, structural and membrane proteins and enzymes involved in gene silencing, protein synthesis/degradation, mitochondrial protein import (metaxin 2), detoxification (GST-pi) and energy metabolism. Stroma-associated proteins were generally underexpressed. The protein expression of prostate cancer correlates with tumor differentiation. Potential prognostic markers may be found among proteins that are differentially expressed and the clinical value of these should be validated.     

Cross-species microarray hybridization to identify developmentally regulated genes in the filamentous fungus <Emphasis Type="Italic">Sordaria macrospora</Emphasis>

The filamentous fungus Sordaria macrospora forms complex three-dimensional fruiting bodies that protect the developing ascospores and ensure their proper discharge. Several regulatory genes essential for fruiting body development were previously isolated by complementation of the sterile mutants pro1, pro11 and pro22. To establish the genetic relationships between these genes and to identify downstream targets, we have conducted cross-species microarray hybridizations using cDNA arrays derived from the closely related fungus Neurospora crassa and RNA probes prepared from wild-type S. macrospora and the three developmental mutants. Of the 1,420 genes which gave a signal with the probes from all the strains used, 172 (12%) were regulated differently in at least one of the three mutants compared to the wild type, and 17 (1.2%) were regulated differently in all three mutant strains. Microarray data were verified by Northern analysis or quantitative real time PCR. Among the genes that are up- or down-regulated in the mutant strains are genes encoding the pheromone precursors, enzymes involved in melanin biosynthesis and a lectin-like protein. Analysis of gene expression in double mutants revealed a complex network of interaction between the pro gene products.     

Genomic and proteomic expression patterns in HPV-16 E6 gene transfected stable human carcinoma cell lines

cDNA microarray and proteomics studies were performed to analyze the genomic and proteomic expression patterns in HPV-16 E6 gene transfected stable human carcinoma cell lines. Among 1024 known genes and ESTs tested by cDNA microarray, we found 50 upregulated and 35 downregulated genes in RC10.1 HPV-16 E6 transfected human colon adenocarcinoma cells compared to RKO cells, and 27 upregulated and 43 downregulated genes in A549E6 HPV-16 E6 transfected human lung adenocarcinoma cells compared to A549 cells. Employing two dimensional gel electrophoresis and MALDI-TOF-MS, the global pattern of protein expressions in RC10.1 human colon adenocarcinoma and A549E6 human lung adenocarcinoma cell lines stably expressing the HPV 16-E6 gene were compared with those of RKO and A549 cell lines to generate a differential protein expression catalog. We found 13 upregulated and 13 downregulated proteins in RC10.1 (E6-expressing RKO) cells compared to RKO cells and 12 upregulated and 14 downregulated proteins in A549E6 (E6-expressing A549) cells compared to A549 cells. The identified genes and proteins were classified into several groups according to the subcellular function. Expressing pattern of three genes and proteins (CDK5, Bak, and I-TRAF) were matched in both analyses of cDNA microarray and proteomics. These powerful approaches using cDNA microarray and proteomics could provide in-depth information on the impact of HPV-16 E6-related genes and proteins. Differential gene and protein expression patterns by transfection of HPV-16 E6 will provide the nucleus of valuable resource for investigation of the biochemical basis of cervical carcinogenesis. Further understanding of this data base may provide valuable resources for developing novel diagnostic markers and therapeutic targets of cervical cancer.