A Ten-microRNA Expression Signature Predicts Survival in Glioblastoma

Glioblastoma (GBM) is the most common and aggressive primary brain tumor with very poor patient median survival. To identify a microRNA (miRNA) expression signature that can predict GBM patient survival, we analyzed the miRNA expression data of GBM patients (n = 222) derived from The Cancer Genome Atlas (TCGA) dataset. We divided the patients randomly into training and testing sets with equal number in each group. We identified 10 significant miRNAs using Cox regression analysis on the training set and formulated a risk score based on the expression signature of these miRNAs that segregated the patients into high and low risk groups with significantly different survival times (hazard ratio HR] = 2.4; 95% CI = 1.4-3.8; p < 0.0001). Of these 10 miRNAs, 7 were found to be risky miRNAs and 3 were found to be protective. This signature was independently validated in the testing set (HR = 1.7; 95% CI = 1.1-2.8; p = 0.002). GBM patients with high risk scores had overall poor survival compared to the patients with low risk scores. Overall survival among the entire patient set was 35.0% at 2 years, 21.5% at 3 years, 18.5% at 4 years and 11.8% at 5 years in the low risk group, versus 11.0%, 5.5%, 0.0 and 0.0% respectively in the high risk group (HR = 2.0; 95% CI = 1.4-2.8; p < 0.0001). Cox multivariate analysis with patient age as a covariate on the entire patient set identified risk score based on the 10 miRNA expression signature to be an independent predictor of patient survival (HR = 1.120; 95% CI = 1.04-1.20; p = 0.003). Thus we have identified a miRNA expression signature that can predict GBM patient survival. These findings may have implications in the understanding of gliomagenesis, development of targeted therapy and selection of high risk cancer patients for adjuvant therapy.

A DNA methylation prognostic signature of glioblastoma: identification of NPTX2-PTEN-NF-kappa B nexus

Glioblastoma (GBM) is the most common, malignant adult primary tumor with dismal patient survival, yet the molecular determinants of patient survival are poorly characterized. Global methylation profile of GBM samples (our cohort; n = 44) using high-resolution methylation microarrays was carried out. Cox regression analysis identified a 9-gene methylation signature that predicted survival in GBM patients. A risk-score derived from methylation signature predicted survival in univariate analysis in our and The Cancer Genome Atlas (TCGA) cohort. Multivariate analysis identified methylation risk score as an independent survival predictor in TCGA cohort. Methylation risk score stratified the patients into low-risk and high-risk groups with significant survival difference. Network analysis revealed an activated NF-kappa B pathway association with high-risk group. NF-kappa B inhibition reversed glioma chemoresistance, and RNA interference studies identified interleukin-6 and intercellular adhesion molecule-1 as key NF-kappa B targets in imparting chemoresistance. Promoter hypermethylation of neuronal pentraxin II (NPTX2), a risky methylated gene, was confirmed by bisulfite sequencing in GBMs. GBMs and glioma cell lines had low levels of NPTX2 transcripts, which could be reversed upon methylation inhibitor treatment. NPTX2 overexpression induced apoptosis, inhibited proliferation and anchorage-independent growth, and rendered glioma cells chemosensitive. Furthermore, NPTX2 repressed NF-kappa B activity by inhibiting AKT through a p53-PTEN-dependent pathway, thus explaining the hypermethylation and downregulation of NPTX2 in NF-kappa B-activated high-risk GBMs. Taken together, a 9-gene methylation signature was identified as an independent GBM prognosticator and could be used for GBM risk stratification. Prosurvival NF-kappa B pathway activation characterized high-risk patients with poor prognosis, indicating it to be a therapeutic target. (C) 2013 AACR.

The repertoire of protein kinases encoded in the draft version of the human genome: atypical variations and uncommon domain combinations

Background: Phosphorylation by protein kinases is central to cellular signal transduction. Abnormal functioning of kinases has been implicated in developmental disorders and malignancies. Their activity is regulated by second messengers and by the binding of associated domains, which are also influential in translocating the catalytic component to their substrate sites, in mediating interaction with other proteins and carrying out their biological roles. Results: Using sensitive profile-search methods and manual analysis, the human genome has been surveyed for protein kinases. A set of 448 sequences, which show significant similarity to protein kinases and contain the critical residues essential for kinase function, have been selected for an analysis of domain combinations after classifying the kinase domains into subfamilies. The unusual domain combinations in particular kinases suggest their involvement in ubiquitination pathways and alternative modes of regulation for mitogen-activated protein kinase kinases (MAPKKs) and cyclin-dependent kinase (CDK)-like kinases. Previously unexplored kinases have been implicated in osteoblast differentiation and embryonic development on the basis of homology with kinases of known functions from other organisms. Kinases potentially unique to vertebrates are involved in highly evolved processes such as apoptosis, protein translation and tyrosine kinase signaling. In addition to coevolution with the kinase domain, duplication and recruitment of non-catalytic domains is apparent in signaling domains such as the PH, DAG-PE, SH2 and SH3 domains. Conclusions: Expansion of the functional repertoire and possible existence of alternative modes of regulation of certain kinases is suggested by their uncommon domain combinations. Experimental verification of the predicted implications of these kinases could enhance our understanding of their biological roles.

Elucidating the cancer-specific genetic alteration spectrum of glioblastoma derived cell lines from whole exome and RNA sequencing

Cell lines derived from tumor tissues have been used as a valuable system to study gene regulation and cancer development. Comprehensive characterization of the genetic background of cell lines could provide clues on novel genes responsible for carcinogenesis and help in choosing cell lines for particular studies. Here, we have carried out whole exome and RNA sequencing of commonly used glioblastoma (GBM) cell lines (U87, T98G, LN229, U343, U373 and LN18) to unearth single nucleotide variations (SNVs), indels, differential gene expression, gene fusions and RNA editing events. We obtained an average of 41,071 SNVs out of which 1,594 (3.88%) were potentially cancer-specific. The cell lines showed frequent SNVs and indels in some of the genes that are known to be altered in GBM-EGFR, TP53, PTEN, SPTA1 and NF1. Chromatin modifying genes-ATRX, MLL3, MLL4, SETD2 and SRCAP also showed alterations. While no cell line carried IDH1 mutations, five cell lines showed hTERT promoter activating mutations with a concomitant increase in hTERT transcript levels. Five significant gene fusions were found of which NUP93-CYB5B was validated. An average of 18,949 RNA editing events was also obtained. Thus we have generated a comprehensive catalogue of genetic alterations for six GBM cell lines.

Elucidating the cancer-specific genetic alteration spectrum of glioblastoma derived cell lines from whole exome and RNA sequencing

Cell lines derived from tumor tissues have been used as a valuable system to study gene regulation and cancer development. Comprehensive characterization of the genetic background of cell lines could provide clues on novel genes responsible for carcinogenesis and help in choosing cell lines for particular studies. Here, we have carried out whole exome and RNA sequencing of commonly used glioblastoma (GBM) cell lines (U87, T98G, LN229, U343, U373 and LN18) to unearth single nucleotide variations (SNVs), indels, differential gene expression, gene fusions and RNA editing events. We obtained an average of 41,071 SNVs out of which 1,594 (3.88%) were potentially cancer-specific. The cell lines showed frequent SNVs and indels in some of the genes that are known to be altered in GBM-EGFR, TP53, PTEN, SPTA1 and NF1. Chromatin modifying genes-ATRX, MLL3, MLL4, SETD2 and SRCAP also showed alterations. While no cell line carried IDH1 mutations, five cell lines showed hTERT promoter activating mutations with a concomitant increase in hTERT transcript levels. Five significant gene fusions were found of which NUP93-CYB5B was validated. An average of 18,949 RNA editing events was also obtained. Thus we have generated a comprehensive catalogue of genetic alterations for six GBM cell lines.

Repertoire of Protein Kinases Encoded in the Genome of Takifugu rubripes

Takifugu rubripes is teleost fish widely used in comparative genomics to understand the human system better due to its similarities both in number of genes and structure of genes. In this work we survey the fugu genome, and, using sensitive computational approaches, we identify the repertoire of putative protein kinases and classify them into groups and subfamilies. The fugu genome encodes 519 protein kinase-like sequences and this number of putative protein kinases is comparable closely to that of human. However, in spite of its similarities to human kinases at the group level, there are differences at the subfamily level as noted in the case of KIS and DYRK subfamilies which contribute to differences which are specific to the adaptation of the organism. Also, certain unique domain combination of galectin domain and YkA domain suggests alternate mechanisms for immune response and binding to lipoproteins. Lastly, an overall similarity with the MAPK pathway of humans suggests its importance to understand signaling mechanisms in humans. Overall the fugu serves as a good model organism to understand roles of human kinases as far as kinases such as LRRK and IRAK and their associated pathways are concerned.

Sequence and structural basis for chromosomal fragility during translocations in cancer

Chromosomal aberration is considered to be one of the major characteristic features in many cancers. Chromosomal translocation, one type of genomic abnormality, can lead to deregulation of critical genes involved in regulating important physiological functions such as cell proliferation and DNA repair. Although chromosomal translocations were thought to be random events, recent findings suggest that certain regions in the human genome are more susceptible to breakage than others. The possibility of deviation from the usual B-DNA conformation in such fragile regions has been an active area of investigation. This review summarizes the factors that contribute towards the fragility of these regions in the chromosomes, such as DNA sequences and the role of different forms of DNA structures. Proteins responsible for chromosomal fragility, and their mechanism of action are also discussed. The effect of positioning of chromosomes within the nucleus favoring chromosomal translocations and the role of repair mechanisms are also addressed.

Effects of Withania somnifera and Tinospora cordifolia Extracts on the Side Population Phenotype of Human Epithelial Cancer Cells: Toward Targeting Multidrug Resistance in Cancer

Recent reports suggest the existence of a subpopulation of stem-like cancer cells, termed as cancer stem cells (CSCs), which bear functional and phenotypic resemblance with the adult, tissue-resident stem cells. Side population (SP) assay based on differential efflux of Hoechst 33342 has been effectively used for the isolation of CSCs. The drug resistance properties of SP cells are typically due to the increased expression of ABC transporters leading to drug efflux. Conventionally used chemotherapeutic drugs may often leads to an enrichment of SP, revealing their inability to target the drug-resistant SP and CSCs. Thus, identification of agents that can reduce the SP phenotype is currently in vogue in cancer therapeutics. Withania somnifera (WS) and Tinospora cordifolia (TC) have been used in Ayurveda for treating various diseases, including cancer. In the current study, we have investigated the effects of ethanolic (ET) extracts of WS and TC on the cancer SP phenotype. Interestingly, we found significant decrease in SP on treatment with TC-ET, but not with WS-ET. The SP-inhibitory TC-ET was further fractionated into petroleum ether (TC-PET), dichloromethane (TC-DCM), and n-butyl alcohol (TC-nBT) fractions using bioactivity-guided fractionation. Our data revealed that TC-PET and TC-DCM, but not TC-nBT, significantly inhibited SP in a dose-dependent manner. Furthermore, flow cytometry-based functional assays revealed that TC-PET and TC-DCM significantly inhibited ABC-B1 and ABC-G2 transporters and sensitized cancer cells toward chemotherapeutic drug-mediated cytotoxicity. Thus, the TC-PET and TC-DCM may harbor phytochemicals with the potential to reverse the drug-resistant phenotype, thus improving the efficacy of cancer chemotherapy.

Timing the Emergence of Resistance to Anti-HIV Drugs with Large Genetic Barriers

New antiretroviral drugs that offer large genetic barriers to resistance, such as the recently approved inhibitors of HIV-1 protease, tipranavir and darunavir, present promising weapons to avert the failure of current therapies for HIV infection. Optimal treatment strategies with the new drugs, however, are yet to be established. A key limitation is the poor understanding of the process by which HIV surmounts large genetic barriers to resistance. Extant models of HIV dynamics are predicated on the predominance of deterministic forces underlying the emergence of resistant genomes. In contrast, stochastic forces may dominate, especially when the genetic barrier is large, and delay the emergence of resistant genomes. We develop a mathematical model of HIV dynamics under the influence of an antiretroviral drug to predict the waiting time for the emergence of genomes that carry the requisite mutations to overcome the genetic barrier of the drug. We apply our model to describe the development of resistance to tipranavir in in vitro serial passage experiments. Model predictions of the times of emergence of different mutant genomes with increasing resistance to tipranavir are in quantitative agreement with experiments, indicating that our model captures the dynamics of the development of resistance to antiretroviral drugs accurately. Further, model predictions provide insights into the influence of underlying evolutionary processes such as recombination on the development of resistance, and suggest guidelines for drug design: drugs that offer large genetic barriers to resistance with resistance sites tightly localized on the viral genome and exhibiting positive epistatic interactions maximally inhibit the emergence of resistant genomes.

Positional Preferences of Polypurine/Polypyrimidine Tracts in Saccharomyces cerevisiae Genome: Implications for cis Regulation of Gene Expression

The complete genome of the baker's yeast S. cerevisiae was analyzed for the presence of polypurine/polypyrimidine (poly[pu/py]) repeats and their occurrences were classified on the basis of their location within and outside open reading frames (ORFs). The analysis reveals that such sequence motifs are present abundantly both in coding as well as noncoding regions. Clear positional preferences are seen when these tracts occur in noncoding regions. These motifs appear to occur predominantly at a unit nucleosomal length both upstream and downstream of ORFs. Moreover, there is a biased distribution of polypurines in the coding strands when these motifs occur within open reading frames. The significance of the biased distribution is discussed with reference to the occurrence of these motifs in other known mRNA sequences and expressed sequence tags. A model for cis regulation of gene expression is proposed based on the ability of these motifs to form an intermolecular triple helix structure when present within the coding region and/or to modulate nucleosome positioning via enhanced histone affinity when present outside coding regions.

A male-specific nuclease-resistant chromatin fraction in the mealybug Planococcus lilacinus

In mealybugs, chromatin condensation is related to both genomic imprinting and sex determination. The paternal chromosomal complement is condensed and genetically inactive in sons but not in daughters. During a study of chromatin organization in Planococcus lilacinus, digestion with micrococcal nuclease showed that 3% to 5% of the male genome is resistant to the enzyme. This Nuclease Resistant Chromatin (NRC) apparently has a nucleosomal organization. Southern hybridization of genomic DNA suggests that NRC sequences are present in both sexes and occur throughout the genome. Cloned NRC DNA is A+T-rich with stretches of adenines similar to those present in mouse alpha-satellite sequences. NRC DNA also contains sequence motifs that are typically associated with the nuclear matrix. Salt-fractionation experiments showed that NRC sequences are matrix associated. These observations are discussed in relation to the unusual cytological features of mealybug chromosomes, including the possible existence of multiple centres of inactivation.

Structural Studies on the Second Mycobacterium smegmatis Dps: Invariant and Variable Features of Structure, Assembly and Function

A second DNA binding protein from stationary-phase cells of Mycobacterium smegmatis (MsDps2) has been identified from the bacterial genome. It was cloned, expressed and characterised and its crystal structure was determined. The core dodecameric structure of MsDps2 is the same as that of the Dps from the organism described earlier (MsDps1). However, MsDps2 possesses a long N-terminal tail instead of the C-terminal tail in MsDps1. This tail appears to be involved in DNA binding. It is also intimately involved in stabilizing the dodecamer. Partly on account of this factor, MsDps2 assembles straightway into the dodecamer, while MsDps1 does so on incubation after going through an intermediate trimeric stage. The ferroxidation centre is similar in the two proteins, while the pores leading to it exhibit some difference. The mode of sequestration of DNA in the crystalline array of molecules, as evidenced by the crystal structures, appears to be different in MsDps1 and MsDps2, highlighting the variability in the mode of Dps–DNA complexation. A sequence search led to the identification of 300 Dps molecules in bacteria with known genome sequences. Fifty bacteria contain two or more types of Dps molecules each, while 195 contain only one type. Some bacteria, notably some pathogenic ones, do not contain Dps. A sequence signature for Dps could also be derived from the analysis.

Genome-wide analysis and experimentation of plant serine/threonine/tyrosine-specific protein kinases

Protein tyrosine phosphorylation plays an important role in cell growth, development and oncogenesis. No classical protein tyrosine kinase has hitherto been cloned from plants. Does protein tyrosine kinase exist in plants? To address this, we have performed a genomic survey of protein tyrosine kinase motifs in plants using the delineated tyrosine phosphorylation motifs from the animal system. The Arabidopsis thaliana genome encodes 57 different protein kinases that have tyrosine kinase motifs. Animal non-receptor tyrosine kinases, SRC, ABL, LYN, FES, SEK, KIN and RAS have structural relationship with putative plant tyrosine kinases. In an extended analysis, animal receptor and non-receptor kinases, Raf and Ras kinases, mixed lineage kinases and plant serine/threonine/tyrosine (STY) protein kinases, form a well-supported group sharing a common origin within the superfamily of STY kinases. We report that plants lack bona fide tyrosine kinases, which raise an intriguing possibility that tyrosine phosphorylation is carried out by dual-specificity STY protein kinases in plants. The distribution pattern of STY protein kinase families on Arabidopsis chromosomes indicates that this gene family is partly a consequence of duplication and reshuffling of the Arabidopsis genome and of the generation of tandem repeats. Genome-wide analysis is supported by the functional expression and characterization of At2g24360 and phosphoproteomics of Arabidopsis. Evidence for tyrosine phosphorylated proteins is provided by alkaline hydrolysis, anti-phosphotyrosine immunoblotting, phosphoamino acid analysis and peptide mass fingerprinting. These results report the first comprehensive survey of genome-wide and tyrosine phosphoproteome analysis of plant STY protein kinases.

Effect of isonicotinic acid hydrazide-copper complex on Rous sarcoma virus and its genome RNA

The copper complex of the antituberculous drug, isonicotinic acid hydrazide (INH), inhibits the RNA-dependent DNA polymerase of Rous sarcoma virus and inactivates its ability to malignantly transform chick embryo cells. The INH-copper complex binds to the 70S genome RNA of Rous sarcoma virus (RSV), which may account for its ability to inhibit the RNA-dependent DNA polymerase. The complex binds RNA more effectively than DNA in contrast to M-IBT-copper complexes, which bind both types of nucleic acids equally. The homopolymers, poly rA and poly rU, are bound by the INH-copper complex to a greater extent than poly rC. Isonicotinic acid hydrazide alone and CuSO4 alone bind neither DNA, RNA, poly (rA), poly (rU), nor poly (rC). However, CuSO4 alone binds poly (rI); INH alone does not. In addition to viral DNA synthesis, chick-embryo cell DNA synthesis is inhibited by the INH-copper complex. The extent of inhibition of cellular DNA synthesis is greater than that of cellular RNA and protein synthesis. No selective inhibition of transformation in cells previously infected with Rous sarcoma virus is observed.