Impact of RNA degradation on gene expression profiling

Background

Cancer is a heterogeneous disease caused by genomic aberrations and characterized by significant variability in clinical outcomes and response to therapies. Several subtypes of common cancers have been identified based on alterations of individual cancer genes, such as HER2, EGFR, and others. However, cancer is a complex disease driven by the interaction of multiple genes, so the copy number status of individual genes is not sufficient to define cancer subtypes and predict responses to treatments. A classification based on genome-wide copy number patterns would be better suited for this purpose.

Method

To develop a more comprehensive cancer taxonomy based on genome-wide patterns of copy number abnormalities, we designed an unsupervised classification algorithm that identifies genomic subgroups of tumors. This algorithm is based on a modified genomic Non-negative Matrix Factorization (gNMF) algorithm and includes several additional components, namely a pilot hierarchical clustering procedure to determine the number of clusters, a multiple random initiation scheme, a new stop criterion for the core gNMF, as well as a 10-fold cross-validation stability test for quality assessment.

Result

We applied our algorithm to identify genomic subgroups of three major cancer types: non-small cell lung carcinoma (NSCLC), colorectal cancer (CRC), and malignant melanoma. High-density SNP array datasets for patient tumors and established cell lines were used to define genomic subclasses of the diseases and identify cell lines representative of each genomic subtype. The algorithm was compared with several traditional clustering methods and showed improved performance. To validate our genomic taxonomy of NSCLC, we correlated the genomic classification with disease outcomes. Overall survival time and time to recurrence were shown to differ significantly between the genomic subtypes.

Conclusions

We developed an algorithm for cancer classification based on genome-wide patterns of copy number aberrations and demonstrated its superiority to existing clustering methods. The algorithm was applied to define genomic subgroups of three cancer types and identify cell lines representative of these subgroups. Our data enabled the assembly of representative cell line panels for testing drug candidates.     

Functional expression and characterization of a purine nucleobase transporter gene from Leishmania major

Leishmania major, like all the other kinetoplastid protozoa, are unable to synthesize purines and rely on purine nucleobase and nucleoside acquisition across the parasite plasma membrane by specific permeases. Although, several genes have been cloned that encode nucleoside transporters in Leishmania and Trypanosoma brucei, much less progress has been made on nucleobase transporters, especially at the molecular level. The studies reported here have cloned and expressed the first gene for a L. major nucleobase transporter, designated LmaNT3. The LmaNT3 permease shows 33% identity to L. donovani nucleoside transporter 1.1 (LdNT1.1) and is, thus, a member of the equilibrative nucleoside transporter (ENT) family. ENT family members identified to date are nucleoside transporters, some of which also transport one or several nucleobases. Functional expression studies in Xenopus laevis oocytes revealed that LmaNT3 mediates high levels of uptake of hypoxanthine, xanthine, adenine and guanine. Moreover, LmaNT3 is an high affinity transporter with K(m) values for hypoxanthine, xanthine, adenine and guanine of 16.5 +/- 1.5, 8.5 +/- 0.6, 8.5 +/- 1.1, and 8.8 +/- 4.0 microM, respectively. LmaNT3 is, thus, the first member of the ENT family identified in any organism that functions as a nucleobase rather than nucleoside or nucleoside/nucleobase transporter.     

Genetic variation in RPOIILS gene encoding RNA polymerase II largest subunit from Leishmania major

Leishmaniasis is a geographically widespread severe disease which includes visceral leishmaniasis, cutaneous leishmaniasis (CL). There are 350 million people at risk in over 80 countries. In the Old World, CL is usually caused by Leishmania major, Leishmania tropica, and Leishmania aethiopica complex which 90 % of cases occurring in Afghanistan, Algeria, Iran, Iraq, Saudi Arabia, Syria, Brazil, and Peru. Recently, some reports showed that some strains of L. major have internal transcribed space (ITS-1) with differential size exhibiting homology with the related gene in a divergent genus of kinetoplastida, the Crithidia. This prompted us to analyze the mentioned gene in 100 isolates obtained from patients with suspected CL. After obtaining samples from 100 patients, DNA extraction was performed and ITS-1 was analyzed using PCR–RFLP. These samples were sequenced for verifying their homology. Then, RPOIILS gene was analyzed in the samples that their ITS-1 gene exhibiting homology with the related gene in Crithidia. Results showed that 10 % of the isolates have ITS-1 exhibiting different size with the routine ones. Sequencing of them showed their similarity to the one from Crithidia fasciculata. RPOIILS gene encoding RNA polymerase II largest subunit analysis showed genetic diversity. This study might also help in solving the problems concerning Leishmaniasis outbreak currently facing in Iran and some other endemic regions of the world.     

General suppression of macrophage gene expression during Leishmania donovani infection

Within the mammalian host, Leishmania donovani is an obligatory intracellular protozoan that resides and multiplies exclusively in the phagolysosomes of macrophages. The outcome of this infection is governed by the interaction between Leishmania and macrophage molecules that ultimately effect the expression of genes within both cells. To explore the effect of this intracellular infection on macrophage gene expression, a cDNA expression array analysis was performed to compare gene expression profiles in noninfected and L. donovani-infected macrophages. In this manner, it was possible to examine the effect of infection on the expression of several hundred well-characterized host cell genes in an unbiased manner. Interestingly, approximately 40% of the genes whose expression was detected in macrophages were down-regulated during infection with L. donovani. However, several genes were also induced during the infection process, some of which could play a role in recruitment of additional macrophages to the site of infection. Taken together, the general suppression of gene expression in addition to the selective induction of key genes is likely to play an important role in allowing the parasite to survive and proliferate within its host macrophage cell.     

Fibronectin gene expression during limb cartilage differentiation

A critical event in limb cartilage differentiation is a transient cellular condensation process in which prechondrogenic mesenchymal cells become closely juxtaposed and interact with one another prior to initiating cartilage matrix deposition. Fibronectin (FN) has been suggested to be involved in regulating the onset of condensation and chondrogenesis by actively promoting prechondrogenic aggregate formation during the process. We have performed a systematic quantitative study of the expression of the FN gene during the progression of chondrogenesis in vitro and in vivo. In high-density micromass cultures of limb mesenchymal cells, FN mRNA levels increase about 5-fold coincident with the crucial condensation process, and remain relatively high during the initial deposition of cartilage matrix by the cells. Thereafter, FN mRNA levels progressively decline to relatively low levels as the cultures form a virtually uniform mass of cartilage. The changes in FN mRNA levels in vitro are paralleled closely by changes in the relative rate of FN synthesis as determined by pulse-labeling and immunoprecipitation analysis. The relative rate of FN synthesis increases 4- to 5-fold at condensation and the onset of chondrogenesis, after which it progressively declines to low levels as cartilage matrix accumulates. High levels of FN gene expression also occur at the onset of chondrogenesis in vivo. In the proximal central core regions of the limb bud in which condensation and cartilage matrix deposition are being initiated, FN mRNA levels and the relative rates of FN synthesis become progressively about 4-fold higher than in the distal subridge region, which consists of undifferentiated mesenchymal cells that have not yet initiated condensation.(ABSTRACT TRUNCATED AT 250 WORDS)     

Integrated genomic and gene expression profiling identifies two major genomic circuits in urothelial carcinoma

Similar to other malignancies, urothelial carcinoma (UC) is characterized by specific recurrent chromosomal aberrations and gene mutations. However, the interconnection between specific genomic alterations, and how patterns of chromosomal alterations adhere to different molecular subgroups of UC, is less clear. We applied tiling resolution array CGH to 146 cases of UC and identified a number of regions harboring recurrent focal genomic amplifications and deletions. Several potential oncogenes were included in the amplified regions, including known oncogenes like E2F3, CCND1, and CCNE1, as well as new candidate genes, such as SETDB1 (1q21), and BCL2L1 (20q11). We next combined genome profiling with global gene expression, gene mutation, and protein expression data and identified two major genomic circuits operating in urothelial carcinoma. The first circuit was characterized by FGFR3 alterations, overexpression of CCND1, and 9q and CDKN2A deletions. The second circuit was defined by E3F3 amplifications and RB1 deletions, as well as gains of 5p, deletions at PTEN and 2q36, 16q, 20q, and elevated CDKN2A levels. TP53/MDM2 alterations were common for advanced tumors within the two circuits. Our data also suggest a possible RAS/RAF circuit. The tumors with worst prognosis showed a gene expression profile that indicated a keratinized phenotype. Taken together, our integrative approach revealed at least two separate networks of genomic alterations linked to the molecular diversity seen in UC, and that these circuits may reflect distinct pathways of tumor development.