Tymon Rubel

Modeling oncogenic signaling in colon tumors by multidirectional analyses of microarray data directed for maximization of analytical reliability.

Background Clinical progression of colorectal cancers (CRC) may occur in parallel with distinctive signaling alterations. We designed multidirectional analyses integrating microarray-based data with biostatistics and bioinformatics to elucidate the signaling and metabolic alterations underlying CRC development in the adenoma-carcinoma sequence. Methodology/Principal Findings Studies were performed on normal mucosa, adenoma, and carcinoma samples obtained during surgery or colonoscopy. Collections of cryostat sections prepared from the tissue samples were evaluated by a pathologist to control the relative cell type content. The measurements were done using Affymetrix GeneChip HG-U133plus2, and probe set data was generated using two normalization algorithms: MAS5.0 and GCRMA with least-variant set (LVS). The data was evaluated using pair-wise comparisons and data decomposition into singular value decomposition (SVD) modes. The method selected for the functional analysis used the Kolmogorov-Smirnov test. Expressional profiles obtained in 105 samples of whole tissue sections were used to establish oncogenic signaling alterations in progression of CRC, while those representing 40 microdissected specimens were used to select differences in KEGG pathways between epithelium and mucosa. Based on a consensus of the results obtained by two normalization algorithms, and two probe set sorting criteria, we identified 14 and 17 KEGG signaling and metabolic pathways that are significantly altered between normal and tumor samples and between benign and malignant tumors, respectively. Several of them were also selected from the raw microarray data of 2 recently published studies (GSE4183 and GSE8671). Conclusion/Significance Although the proposed strategy is computationally complex and labor–intensive, it may reduce the number of false results.

Integrating proteomic and transcriptomic high-throughput surveys for search of new biomarkers of colon tumors

To the search of new colon tumor biomarkers in the transition from normal colon (NC) mucosa to adenoma (AD) and adenocarcinoma (AC), we integrated microarray data with the results of a high-throughput proteomic workflow. In proteomic study, we used a modified isoelectric focusing protocol on strips with an immobilized pH gradient to separate peptides labeled with iTRAQ (isobaric tags for relative and absolute quantitation) tags followed by liquid chromatography–tandem mass spectrometry analysis. Gene expression measurements were done using Affymetrix GeneChip HG-U133plus2 microarrays and quantitative reverse transcriptase PCR (q-RT-PCR). We identified 3,886 proteins with at least two peptides. Of them, 1,061 proteins were differentially expressed [FC ≥ 1.5; FDR ≤ 0.01] in two pair-wise comparisons: AD vs. NC and AC vs. AD while 15 and 23 proteins were progressively up-regulated and down-regulated in the NC/AD/AC sequence, respectively. The quantitative proteomic information was subsequently correlated with microarray data. For a collection of genes with the same direction of changes of both mRNA and protein levels, we obtained 785/853/795 genes in AD vs. NC/AC vs. NC/AC vs. AD comparison, respectively. Further evaluation of sequentially altered gene expression by q-RT-PCR on individual samples of 24 NCs, 42 ADs, and 26 ACs confirmed progressive expression of six genes: biglycan, calumenin, collagen type XII, alpha 1 (COL12A1), monoamine oxidase A (MAOA), ectonucleoside triphosphate diphosphohydrolase 5 (ENTPD5), and MOCO sulphurase C-terminal domain-containing 2 (MOSC2). Among them, three continuously down-regulated (MAOA, ENTPD5, and MOSC2) and one continuously overexpressed (COL12A1) are reported, to our best knowledge, for the first time in a connection to colon cancer onset.

Molecular defense mechanisms of Barrett’s metaplasia estimated by an integrative genomics

Barrett’s esophagus is characterized by the replacement of squamous epithelium with specialized intestinal metaplastic mucosa. The exact mechanisms of initiation and development of Barrett’s metaplasia remain unknown, but a hypothesis of “successful adaptation” against noxious reflux components has been proposed. To search for the repertoire of adaptation mechanisms of Barrett’s metaplasia, we employed high-throughput functional genomic and proteomic methods that defined the molecular background of metaplastic mucosa resistance to reflux. Transcriptional profiling was established for 23 pairs of esophageal squamous epithelium and Barrett’s metaplasia tissue samples using Affymetrix U133A 2.0 GeneChips and validated by quantitative real-time polymerase chain reaction. Differences in protein composition were assessed by electrophoretic and mass-spectrometry-based methods. Among 2,822 genes differentially expressed between Barrett’s metaplasia and squamous epithelium, we observed significantly overexpressed metaplastic mucosa genes that encode cytokines and growth factors, constituents of extracellular matrix, basement membrane and tight junctions, and proteins involved in prostaglandin and phosphoinositol metabolism, nitric oxide production, and bioenergetics. Their expression likely reflects defense and repair responses of metaplastic mucosa, whereas overexpression of genes encoding heat shock proteins and several protein kinases in squamous epithelium may reflect lower resistance of normal esophageal epithelium than Barrett’s metaplasia to reflux components. Despite the methodological and interpretative difficulties in data analyses discussed in this paper, our studies confirm that Barrett’s metaplasia may be regarded as a specific microevolution allowing for accumulation of mucosal morphological and physiological changes that better protect against reflux injury.

Comprehensive Analysis of the Palindromic Motif TCTCGCGAGA: A Regulatory Element of the HNRNPK Promoter

Definitive identification of promoters, their cis-regulatory motifs, and their trans-acting proteins requires experimental analysis. To define the HNRNPK promoter and its cognate DNA–protein interactions, we performed a comprehensive study combining experimental approaches, including luciferase reporter gene assays, chromatin immunoprecipitations (ChIP), electrophoretic mobility shift assays (EMSA), and mass spectrometry (MS). We discovered that out of the four potential HNRNPK promoters tested, the one containing the palindromic motif TCTCGCGAGA exhibited the highest activity in a reporter system assay. Although further EMSA and MS analyses, performed to uncover the identity of the palindrome-binding transcription factor, did identify a complex of DNA-binding proteins, neither method unambiguously identified the pertinent direct trans-acting protein(s). ChIP revealed similar chromatin states at the promoters with the palindromic motif and at housekeeping gene promoters. A ChIP survey showed significantly higher recruitment of PARP1, a protein identified by MS as ubiquitously attached to DNA probes, within heterochromatin sites. Computational analyses indicated that this palindrome displays features that mark nucleosome boundaries, causing the surrounding DNA landscape to be constitutively open. Our strategy of diverse approaches facilitated the direct characterization of various molecular properties of HNRNPK promoter bearing the palindromic motif TCTCGCGAGA, despite the obstacles that accompany in vitro methods.

Histone H3 lysine 27 acetylation is altered in colon cancer

BackgroundHistone post-translational modifications (PTMs) play an important role in the regulation of the expression of genes, including those involved in cancer development and progression. However, our knowledge of PTM patterns in human tumours is limited.MethodsMS-based analyses were used to quantify global alterations of histone PTMs in colorectal cancer (CRC) samples. Histones isolated from 12 CRCs and their corresponding normal mucosa by acidic extraction were separated by SDS-PAGE and analysed by liquid chromatography-mass spectrometry.ResultsAmong 96 modified peptides, 41 distinct PTM sites were identified, of which 7, 13, 11, and 10 were located within the H2A, H2B, H3, and H4 sequences, respectively, and distributed among the amino-terminal tails and the globular domain of the four histones. Modification intensities were quantified for 33 sites, of which 4 showed significant (p-value ≤ 0.05) differences between CRC tissues and healthy mucosa samples. We identified histone H3 lysine 27 acetylation (H3K27Ac) as a modification upregulated in CRC, which had not been shown previously.ConclusionsThe present results indicate the usefulness of a bottom-up proteomic approach for the detection of histone modifications at a global scale. The differential abundance of H3K27Ac mark in CRC, a PTM associated with active enhancers, suggests its role in regulating genes whose expression changes in CRC.

An integrated LC-ESI-MS platform for quantitation of serum peptide ladders. Application for colon carcinoma study.

Mounting evidence indicates that MS analysis of the human blood peptidome allows to distinguish between cancer and non‐cancer samples, giving promise for a new MS‐based diagnostic tool. However, several aspects of already published work have been criticized and demand for more methodical approach has been formulated. Motivated by this we undertook a systematic study of the plasma and serum peptidome using an integrated ESI‐LC‐MS‐based platform, equipped with new data analysis tools for relative and absolute peptide quantitation. We used a high resolution LC‐ESI‐MS to analyze well‐separated MS signals corresponding to peptides, and measured the variability of >1000 peptide signal amplitudes across a set of plasma and serum samples from healthy individuals. By spiking serum samples with known amounts of isotopically labeled versions of a selected set of peptides we measured the variability of their absolute concentration in this sample set and demonstrated a strong influence of clotting time on the concentration of these peptides in serum. Finally, we used this new LC‐ESI‐MS analytical platform for the differential analysis of healthy versus colon cancer serum samples and found that it was possible to distinguish the two groups with 89.8% sensitivity and 94.6% specificity.

Diffprot - software for non-parametric statistical analysis of differential proteomics data.

Mass spectrometry-based global proteomics experiments generate large sets of data that can be converted into useful information only with an appropriate statistical approach. We present Diffprot - a software tool for statistical analysis of MS-derived quantitative data. With implemented resampling-based statistical test and local variance estimate, Diffprot allows to draw significant results from small scale experiments and effectively eliminates false positive results. To demonstrate the advantages of this software, we performed two spike-in tests with complex biological matrices, one label-free and one based on iTRAQ quantification; in addition, we performed an iTRAQ experiment on bacterial samples. In the spike-in tests, protein ratios were estimated and were in good agreement with theoretical values; statistical significance was assigned to spiked proteins and single or no false positive results were obtained with Diffprot. We compared the performance of Diffprot with other statistical tests - widely used t-test and non-parametric Wilcoxon test. In contrast to Diffprot, both generated many false positive hits in the spike-in experiment. This proved the superiority of the resampling-based method in terms of specificity, making Diffprot a rational choice for small scale high-throughput experiments, when the need to control the false positive rate is particularly pressing.

Comparative kinome analysis to identify putative colon tumor biomarkers

Kinase domains are the type of protein domain most commonly found in genes associated with tumorigenesis. Because of this, the human kinome (the protein kinase component of the genome) represents a promising source of cancer biomarkers and potential targets for novel anti-cancer therapies. Alterations in the human colon kinome during the progression from normal colon (NC) through adenoma (AD) to adenocarcinoma (AC) were investigated using integrated transcriptomic and proteomic datasets. Two hundred thirty kinase genes and 42 kinase proteins showed differential expression patterns (fold change ≥ 1.5) in at least one tissue pair-wise comparison (AD vs. NC, AC vs. NC, and/or AC vs. AD). Kinases that exhibited similar trends in expression at both the mRNA and protein levels were further analyzed in individual samples of NC (n = 20), AD (n = 39), and AC (n = 24) by quantitative reverse transcriptase PCR. Individual samples of NC and tumor tissue were distinguishable based on the mRNA levels of a set of 20 kinases. Altered expression of several of these kinases, including chaperone activity of bc1 complex-like (CABC1) kinase, bromodomain adjacent to zinc finger domain protein 1B (BAZ1B) kinase, calcium/calmodulin-dependent protein kinase type II subunit delta (CAMK2D), serine/threonine-protein kinase 24 (STK24), vaccinia-related kinase 3 (VRK3), and TAO kinase 3 (TAOK3), has not been previously reported in tumor tissue. These findings may have diagnostic potential and may lead to the development of novel targeted therapeutic interventions for colorectal cancer.

Limited prolonged effects of rifaximin treatment on irritable bowel syndrome-related differences in the fecal microbiome and metabolome

ABSTRACT Irritable bowel syndrome (IBS) is a chronic functional disorder and its development may be linked, directly and indirectly, to intestinal dysbiosis. Here we investigated the interactions between IBS symptoms and the gut microbiome, including the relation to rifaximin (1200 mg daily; 11.2 g per a treatment). We recruited 72 patients, including 31 with IBS-D (diarrhea), 11 with IBS-C (constipation), and 30 with IBS-M (mixed constipation and diarrhea) and 30 healthy controls (HCs). Of them, 68%, 64%, and 53% patients with IBS-D, IBS-C, and IBS-M, respectively, achieved 10–12 week-term improvement after the rifaximin treatment. Stool samples were collected before and after the treatment, and fecal microbiotic profiles were analyzed by deep sequencing of 16S rRNA, while stool metabolic profiles were studied by hydrogen 1-nuclear magnetic resonance (1H-NMR) and gas chromatography–mass spectrometry (GC-MS). Of 26 identified phyla, only Bacteroidetes, Firmicutes, Proteobacteria, and Actinobacteria were consistently found in all samples. Bacteroidetes was predominant in fecal samples from HCs and IBS-D and IBS-M subjects, whereas Firmicutes was predominant in samples from IBS-C subjects. Species richness, but not community diversity, differentiated all IBS patients from HCs. Metabolic fingerprinting, using NMR spectra, distinguished HCs from all IBS patients. Thirteen metabolites identified by GC-MS differed HCs and IBS patients. However, neither metagenomics nor metabolomics analyses identified significant differences between patients with and without improvement after treatment.

Functional features of gene expression profiles differentiating gastrointestinal stromal tumours according to KIT mutations and expression

BackgroundGastrointestinal stromal tumours (GISTs) represent a heterogeneous group of tumours of mesenchymal origin characterized by gain-of-function mutations in KIT or PDGFRA of the type III receptor tyrosine kinase family. Although mutations in either receptor are thought to drive an early oncogenic event through similar pathways, two previous studies reported the mutation-specific gene expression profiles. However, their further conclusions were rather discordant. To clarify the molecular characteristics of differentially expressed genes according to GIST receptor mutations, we combined microarray-based analysis with detailed functional annotations.MethodsTotal RNA was isolated from 29 frozen gastric GISTs and processed for hybridization on GENECHIP® HG-U133 Plus 2.0 microarrays (Affymetrix). KIT and PDGFRA were analyzed by sequencing, while related mRNA levels were analyzed by quantitative RT-PCR.ResultsFifteen and eleven tumours possessed mutations in KIT and PDGFRA, respectively; no mutation was found in three tumours. Gene expression analysis identified no discriminative profiles associated with clinical or pathological parameters, even though expression of hundreds of genes differentiated tumour receptor mutation and expression status. Functional features of genes differentially expressed between the two groups of GISTs suggested alterations in angiogenesis and G-protein-related and calcium signalling.ConclusionOur study has identified novel molecular elements likely to be involved in receptor-dependent GIST development and allowed confirmation of previously published results. These elements may be potential therapeutic targets and novel markers of KIT mutation status.