TRIB2 in human AML: a biological and clinical investigation

Acute myeloid leukemia (AML) involves the proliferation, abnormal survival and arrest of cells at a very early stage of myeloid cell differentiation. The biological and clinical heterogeneity of this disease complicates treatment and highlights the significance of understanding the underlying causes of AML, which may constitute potential therapeutic targets, as well as offer prognostic information. Tribbles homolog 2 (Trib2) is a potent murine oncogene capable of inducing transplantable AML with complete penetrance. The pathogenicity of Trib2 is attributed to its ability to induce proteasomal degradation of the full length isoform of the transcription factor CCAAT/enhancer-binding protein alpha (C/EBPα p42). The role of TRIB2 in human AML cells, however, has not been systematically investigated or targeted. Across human cancers, TRIB2 oncogenic activity was found to be associated with its elevated expression. In the context of AML, TRIB2 overexpression was suggested to be associated with the large and heterogeneous subset of cytogenetically normal AML patients. Based upon the observation that overexpression of TRIB2 has a role in cellular transformation, the effect of modulating its expression in human AML was examined in a human AML cell line that expresses high levels of TRIB2, U937 cells. Specific suppression of TRIB2 led to impaired cell growth, as a consequence of both an increase in apoptosis and a decrease in cell proliferation. Consistent with these in vitro results, TRIB2 silencing strongly reduced progression of the U937 in vivo xenografts, accompanied by detection of a lower spleen weight when compared with mice transplanted with TRIB2- expressing control cells. Gene expression analysis suggested that TRIB2 modulates apoptosis and cell-cycle sensitivity by influencing the expression of a subset of genes known to have implications on these phenotypes. Furthermore, TRIB2 was found to be expressed in a significant subset of AML patient samples analysed. To investigate whether increased expression of this gene could be afforded prognostic significance, primary AML cells with dichotomized levels of TRIB2 transcripts were evaluated in terms of their xenoengraftment potential, an assay reported to correlate with disease aggressiveness observed in humans. A small cohort of analysed samples with higher TRIB2 expression did not associate with preferential leukaemic cell engraftment in highly immune-deficient mice, hence, not predicting for an adverse prognosis. However, further experiments including a larger cohort of well characterized AML patients would be needed to clarify TRIB2 significance in the diagnostic setting. Collectively, these data support a functional role for TRIB2 in the maintenance of the oncogenic properties of human AML cells and suggest TRIB2 can be considered a rational therapeutic target. Proteasome inhibition has emerged as an attractive target for the development of novel anti-cancer therapies and results from translational research and clinical trials support the idea that proteasome inhibitors should be considered in the treatment of AML. The present study argued that proteasome inhibition would effectively inhibit the function of TRIB2 by abrogating C/EBPα p42 protein degradation and that it would be an effective pharmacological targeting strategy in TRIB2-positive AMLs. Here, a number of cell models expressing high levels of TRIB2 were successfully targeted by treatment with proteasome inhibitors, as demonstrated by multiple measurements that included increased cytotoxicity, inhibition of clonogenic growth and anti-AML activity in vivo. Mechanistically, it was shown that block of the TRIB2 degradative function led to an increase of C/EBPα p42 and that response was specific to the TRIB2-C/EBPα axis. Specificity was addressed by a panel of experiments showing that U937 cells (express detectable levels of endogenous TRIB2 and C/EBPα) treated with the proteasome inhibitor bortezomib (Brtz) displayed a higher cytotoxic response upon TRIB2 overexpression and that ectopic expression of C/EBPα rescued cell death. Additionally, in C/EBPα-negative leukaemia cells, K562 and Kasumi 1, Brtz-induced toxicity was not increased following TRIB2 overexpression supporting the specificity of the compound on the TRIB2-C/EBPα axis. Together these findings provide pre-clinical evidence that TRIB2- expressing AML cells can be pharmacologically targeted with proteasome inhibition due, in part, to blockage of the TRIB2 proteolytic function on C/EBPα p42. A large body of evidence indicates that AML arises through the stepwise acquisition of genetic and epigenetic changes. Mass spectrometry data has identified an interaction between TRIB2 and the epigenetic regulator Protein Arginine Methyltransferase 5 (PRMT5). Following assessment of TRIB2‟s role in AML cell survival and effective targeting of the TRIB2-C/EBPα degradation pathway, a putative TRIB2/PRMT5 cooperation was investigated in order to gain a deeper understanding of the molecular network in which TRIB2 acts as a potent myeloid oncogene. First, a microarray data set was interrogated for PRMT5 expression levels and the primary enzyme responsible for symmetric dimethylation was found to be transcribed at significantly higher levels in AML patients when compared to healthy controls. Next, depletion of PRMT5 in the U937 cell line was shown to reduce the transformative phenotype in the high expressing TRIB2 AML cells, which suggests that PRMT5 and TRIB2 may cooperate to maintain the leukaemogenic potential. Importantly, PRMT5 was identified as a TRIB2-interacting protein by means of a protein tagging approach to purify TRIB2 complexes from 293T cells. These findings trigger further research aimed at understanding the underlying mechanism and the functional significance of this interplay. In summary, the present study provides experimental evidence that TRIB2 has an important oncogenic role in human AML maintenance and, importantly in such a molecularly heterogeneous disease, provides the rational basis to consider proteasome inhibition as an effective targeting strategy for AML patients with high TRIB2 expression. Finally, the identification of PRMT5 as a TRIB2-interacting protein opens a new level of regulation to consider in AML. This work may contribute to our further understanding and therapeutic strategies in acute leukaemias.

Transcriptional profiling in Saccharomyces cerevisiae relevant for predicting alachlor mechanisms of toxicity

Alachlor has been a commonly applied herbicide and is a substance of ecotoxicological concern. The present study aims to identify molecular biomarkers in the eukaryotic model Saccharomyces cerevisiae that can be used to predict potential cytotoxic effects of alachlor, while providing new mechanistic clues with possible relevance for experimentally less accessible eukaryotes. It focuses on genome-wide expression profiling in a yeast population in response to two exposure scenarios exerting effects from slight to moderate magnitude at phenotypic level. In particular, 100 and 264 genes, respectively, were found as differentially expressed on a 2-h exposure of yeast cells to the lowest observed effect concentration (110 mg/L) and the 20% inhibitory concentration (200 mg/L) of alachlor, in comparison with cells not exposed to the herbicide. The datasets of alachlor-responsive genes showed functional enrichment in diverse metabolic, transmembrane transport, cell defense, and detoxification categories. In general, the modifications in transcript levels of selected candidate biomarkers, assessed by quantitative reverse transcriptase polymerase chain reaction, confirmed the microarray data and varied consistently with the growth inhibitory effects of alachlor. Approximately 16% of the proteins encoded by alachlor-differentially expressed genes were found to share significant homology with proteins from ecologically relevant eukaryotic species. The biological relevance of these results is discussed in relation to new insights into the potential adverse effects of alachlor in health of organisms from ecosystems, particularly in worst-case situations such as accidental spills or careless storage, usage, and disposal.

Unraveling the function of bacterial-type phosphoenolpyruvate carboxylase in vascular plants

Thesis (Master, Biology) -- Queen's University, 2016-09-29 20:09:46.997

Autoinducer-2 plays a crucial role in gut colonization and probiotic functionality of Bifidobacterium breve UCC2003

In the present study we show that luxS of Bifidobacterium breve UCC2003 is involved in the production of the interspecies signaling molecule autoinducer-2 (AI-2), and that this gene is essential for gastrointestinal colonization of a murine host, while it is also involved in providing protection against Salmonella infection in Caenorhabditis elegans. We demonstrate that a B. breve luxS-insertion mutant is significantly more susceptible to iron chelators than the WT strain and that this sensitivity can be partially reverted in the presence of the AI-2 precursor DPD. Furthermore, we show that several genes of an iron starvation-induced gene cluster, which are downregulated in the luxS-insertion mutant and which encodes a presumed iron-uptake system, are transcriptionally upregulated under in vivo conditions. Mutation of two genes of this cluster in B. breve UCC2003 renders the derived mutant strains sensitive to iron chelators while deficient in their ability to confer gut pathogen protection to Salmonella-infected nematodes. Since a functional luxS gene is present in all tested members of the genus Bifidobacterium, we conclude that bifidobacteria operate a LuxS-mediated system for gut colonization and pathogen protection that is correlated with iron acquisition.

Bifidobacterium breve UCC2003 metabolises the human milk oligosaccharides lacto-N-tetraose and lacto-N-neo-tetraose through overlapping, yet distinct pathways

In this study, we demonstrate that the prototype B. breve strain UCC2003 possesses specific metabolic pathways for the utilisation of lacto-N-tetraose (LNT) and lacto-N-neotetraose (LNnT), which represent the central moieties of Type I and Type II human milk oligosaccharides (HMOs), respectively. Using a combination of experimental approaches, the enzymatic machinery involved in the metabolism of LNT and LNnT was identified and characterised. Homologs of the key genetic loci involved in the utilisation of these HMO substrates were identified in B. breve, B. bifidum, B. longum subsp. infantis and B. longum subsp. longum using bioinformatic analyses, and were shown to be variably present among other members of the Bifidobacterium genus, with a distinct pattern of conservation among human-associated bifidobacterial species.

Selection bias in gene extraction on the basis of microarray gene-expression data

In the context of cancer diagnosis and treatment, we consider the problem of constructing an accurate prediction rule on the basis of a relatively small number of tumor tissue samples of known type containing the expression data on very many (possibly thousands) genes. Recently, results have been presented in the literature suggesting that it is possible to construct a prediction rule from only a few genes such that it has a negligible prediction error rate. However, in these results the test error or the leave-one-out cross-validated error is calculated without allowance for the selection bias. There is no allowance because the rule is either tested on tissue samples that were used in the first instance to select the genes being used in the rule or because the cross-validation of the rule is not external to the selection process; that is, gene selection is not performed in training the rule at each stage of the cross-validation process. We describe how in practice the selection bias can be assessed and corrected for by either performing a cross-validation or applying the bootstrap external to the selection process. We recommend using 10-fold rather than leave-one-out cross-validation, and concerning the bootstrap, we suggest using the so-called. 632+ bootstrap error estimate designed to handle overfitted prediction rules. Using two published data sets, we demonstrate that when correction is made for the selection bias, the cross-validated error is no longer zero for a subset of only a few genes.

A mixture model-based approach to the clustering of microarray expression data

Motivation: This paper introduces the software EMMIX-GENE that has been developed for the specific purpose of a model-based approach to the clustering of microarray expression data, in particular, of tissue samples on a very large number of genes. The latter is a nonstandard problem in parametric cluster analysis because the dimension of the feature space (the number of genes) is typically much greater than the number of tissues. A feasible approach is provided by first selecting a subset of the genes relevant for the clustering of the tissue samples by fitting mixtures of t distributions to rank the genes in order of increasing size of the likelihood ratio statistic for the test of one versus two components in the mixture model. The imposition of a threshold on the likelihood ratio statistic used in conjunction with a threshold on the size of a cluster allows the selection of a relevant set of genes. However, even this reduced set of genes will usually be too large for a normal mixture model to be fitted directly to the tissues, and so the use of mixtures of factor analyzers is exploited to reduce effectively the dimension of the feature space of genes. Results: The usefulness of the EMMIX-GENE approach for the clustering of tissue samples is demonstrated on two well-known data sets on colon and leukaemia tissues. For both data sets, relevant subsets of the genes are able to be selected that reveal interesting clusterings of the tissues that are either consistent with the external classification of the tissues or with background and biological knowledge of these sets.

Integrating data from heterogeneous DNA microarray platforms

DNA microarrays are one of the most used technologies for gene expression measurement. However, there are several distinct microarray platforms, from different manufacturers, each with its own measurement protocol, resulting in data that can hardly be compared or directly integrated. Data integration from multiple sources aims to improve the assertiveness of statistical tests, reducing the data dimensionality problem. The integration of heterogeneous DNA microarray platforms comprehends a set of tasks that range from the re-annotation of the features used on gene expression, to data normalization and batch effect elimination. In this work, a complete methodology for gene expression data integration and application is proposed, which comprehends a transcript-based re-annotation process and several methods for batch effect attenuation. The integrated data will be used to select the best feature set and learning algorithm for a brain tumor classification case study. The integration will consider data from heterogeneous Agilent and Affymetrix platforms, collected from public gene expression databases, such as The Cancer Genome Atlas and Gene Expression Omnibus.

Test of four colon cancer risk-scores in formalin fixed paraffin embedded microarray gene expression data.

BACKGROUND: Prognosis prediction for resected primary colon cancer is based on the T-stage Node Metastasis (TNM) staging system. We investigated if four well-documented gene expression risk scores can improve patient stratification. METHODS: Microarray-based versions of risk-scores were applied to a large independent cohort of 688 stage II/III tumors from the PETACC-3 trial. Prognostic value for relapse-free survival (RFS), survival after relapse (SAR), and overall survival (OS) was assessed by regression analysis. To assess improvement over a reference, prognostic model was assessed with the area under curve (AUC) of receiver operating characteristic (ROC) curves. All statistical tests were two-sided, except the AUC increase. RESULTS: All four risk scores (RSs) showed a statistically significant association (single-test, P < .0167) with OS or RFS in univariate models, but with HRs below 1.38 per interquartile range. Three scores were predictors of shorter RFS, one of shorter SAR. Each RS could only marginally improve an RFS or OS model with the known factors T-stage, N-stage, and microsatellite instability (MSI) status (AUC gains < 0.025 units). The pairwise interscore discordance was never high (maximal Spearman correlation = 0.563) A combined score showed a trend to higher prognostic value and higher AUC increase for OS (HR = 1.74, 95% confidence interval [CI] = 1.44 to 2.10, P < .001, AUC from 0.6918 to 0.7321) and RFS (HR = 1.56, 95% CI = 1.33 to 1.84, P < .001, AUC from 0.6723 to 0.6945) than any single score. CONCLUSIONS: The four tested gene expression-based risk scores provide prognostic information but contribute only marginally to improving models based on established risk factors. A combination of the risk scores might provide more robust information. Predictors of RFS and SAR might need to be different.

GSVA: gene set variation analysis for microarray and RNA-seq data

Gene set enrichment (GSE) analysis is a popular framework for condensing information from gene expression profiles into a pathway or signature summary. The strengths of this approach over single gene analysis include noise and dimension reduction, as well as greater biological interpretability. As molecular profiling experiments move beyond simple case-control studies, robust and flexible GSE methodologies are needed that can model pathway activity within highly heterogeneous data sets. To address this challenge, we introduce Gene Set Variation Analysis (GSVA), a GSE method that estimates variation of pathway activity over a sample population in an unsupervised manner. We demonstrate the robustness of GSVA in a comparison with current state of the art sample-wise enrichment methods. Further, we provide examples of its utility in differential pathway activity and survival analysis. Lastly, we show how GSVA works analogously with data from both microarray and RNA-seq experiments. GSVA provides increased power to detect subtle pathway activity changes over a sample population in comparison to corresponding methods. While GSE methods are generally regarded as end points of a bioinformatic analysis, GSVA constitutes a starting point to build pathway-centric models of biology. Moreover, GSVA contributes to the current need of GSE methods for RNA-seq data. GSVA is an open source software package for R which forms part of the Bioconductor project and can be downloaded at http://www.bioconductor.org.

Data analysis tools and methods for DNA microarray and high-throughput sequencing data

The recent rapid development of biotechnological approaches has enabled the production of large whole genome level biological data sets. In order to handle thesedata sets, reliable and efficient automated tools and methods for data processingand result interpretation are required. Bioinformatics, as the field of studying andprocessing biological data, tries to answer this need by combining methods and approaches across computer science, statistics, mathematics and engineering to studyand process biological data. The need is also increasing for tools that can be used by the biological researchers themselves who may not have a strong statistical or computational background, which requires creating tools and pipelines with intuitive user interfaces, robust analysis workflows and strong emphasis on result reportingand visualization. Within this thesis, several data analysis tools and methods have been developed for analyzing high-throughput biological data sets. These approaches, coveringseveral aspects of high-throughput data analysis, are specifically aimed for gene expression and genotyping data although in principle they are suitable for analyzing other data types as well. Coherent handling of the data across the various data analysis steps is highly important in order to ensure robust and reliable results. Thus,robust data analysis workflows are also described, putting the developed tools andmethods into a wider context. The choice of the correct analysis method may also depend on the properties of the specific data setandthereforeguidelinesforchoosing an optimal method are given. The data analysis tools, methods and workflows developed within this thesis have been applied to several research studies, of which two representative examplesare included in the thesis. The first study focuses on spermatogenesis in murinetestis and the second one examines cell lineage specification in mouse embryonicstem cells.

Knowledge-based analysis of microarray gene expression data by using support vector machines

We introduce a method of functionally classifying genes by using gene expression data from DNA microarray hybridization experiments. The method is based on the theory of support vector machines (SVMs). SVMs are considered a supervised computer learning method because they exploit prior knowledge of gene function to identify unknown genes of similar function from expression data. SVMs avoid several problems associated with unsupervised clustering methods, such as hierarchical clustering and self-organizing maps. SVMs have many mathematical features that make them attractive for gene expression analysis, including their flexibility in choosing a similarity function, sparseness of solution when dealing with large data sets, the ability to handle large feature spaces, and the ability to identify outliers. We test several SVMs that use different similarity metrics, as well as some other supervised learning methods, and find that the SVMs best identify sets of genes with a common function using expression data. Finally, we use SVMs to predict functional roles for uncharacterized yeast ORFs based on their expression data.

A Comparative Analysis of Predictive Learning Algorithms on High-Dimensional Microarray Cancer Data

This research evaluates pattern recognition techniques on a subclass of big data where the dimensionality of the input space (p) is much larger than the number of observations (n). Specifically, we evaluate massive gene expression microarray cancer data where the ratio κ is less than one. We explore the statistical and computational challenges inherent in these high dimensional low sample size (HDLSS) problems and present statistical machine learning methods used to tackle and circumvent these difficulties. Regularization and kernel algorithms were explored in this research using seven datasets where κ < 1. These techniques require special attention to tuning necessitating several extensions of cross-validation to be investigated to support better predictive performance. While no single algorithm was universally the best predictor, the regularization technique produced lower test errors in five of the seven datasets studied.

Production of DNA microarray and expression analysis of genes from Xylella fastidiosa in different culture media

DNA Microarray was developed to monitor the expression of many genes from Xylella fastidiosa, allowing the side by-side comparison of two situations in a single experiment. The experiments were performed using X. fastidiosa cells grown in two culture media: BCYE and XDM2. The primers were synthesized, spotted onto glass slides and the array was hybridized against fluorescently labeled cDNAs. The emitted signals were quantified, normalized and the data were statistically analyzed to verify the differentially expressed genes. According to the data, 104 genes were differentially expressed in XDM2 and 30 genes in BCYE media. The present study showed that DNA microarray technique efficiently differentiate the expressed genes under different conditions.