Inorganic polyphosphate as an energy source in tumorigenesis

Cancer cells have high demands for energy to maintain their exceedingly proliferative growth. However, the mechanism of energy expenditure in cancer is not well understood. We hypothesize that cancer cells might utilize energy-rich inorganic polyphosphate (polyP), as energetic reserve. PolyP is comprised of orthophosphates linked by phosphoanhydride bonds, as in ATP. Here, we show that polyP is highly abundant in several types of cancer cells, including brain tumor-initiating cells (BTICs), i.e., stem-like cells derived from a mouse brain tumor model that we have previously described. The polymer is avidly consumed during starvation of the BTICs. Depletion of ATP by inhibiting glycolysis and mitochondrial ATP-synthase (OXPHOS) further decreases the levels of polyP and alters morphology of the cells. Moreover, enzymatic hydrolysis of the polymer impairs the viability of cancer cells and significantly deprives ATP stores. These results suggest that polyP might be utilized as a source of phosphate energy in cancer. While the role of polyP as an energy source is established for bacteria, this finding is the first demonstration that polyP may play a similar role in the metabolism of cancer cells.

Global Reduction of H3K4me3 Improves Chemotherapeutic Efficacy for Pediatric Ependymomas

BACKGROUND

Ependymomas (EPNs) are the third most common brain tumor in children. These tumors are resistant to available chemotherapeutic treatments, therefore new effective targeted therapeutics must be identified. Increasing evidence shows epigenetic alterations including histone posttranslational modifications (PTMs), are associated with malignancy, chemotherapeutic resistance and prognosis for pediatric EPNs. In this study we examined histone PTMs in EPNs and identified potential targets to improve chemotherapeutic efficacy.

METHODS

Global histone H3 lysine 4 trimethylation (H3K4me3) levels were detected in pediatric EPN tumor samples with immunohistochemistry and immunoblots. Candidate genes conferring therapeutic resistance were profiled in pediatric EPN tumor samples with micro-array. Promoter H3K4me3 was examined for two candidate genes, CCND1 and ERBB2, with chromatin-immunoprecipitation coupled with real-time PCR (ChIP-PCR). These methods and MTS assay were used to verify a relationship between H3K4me3 levels and CCND1 and ERBB2, and to investigate cell viability in response to chemotherapeutic drugs in primary cultured pediatric EPN cells.

RESULTS

H3K4me3 levels positively correlate with WHO grade malignancy in pediatric EPNs and are associated with progression free survival in patients with posterior fossa group A EPNs (PF-EPN-A). Reduction of H3K4me3 by silencing its methyltransferase SETD1A, in primary cultured EPN cells increased cell response to chemotherapy.

CONCLUSIONS

Our results support the development of a novel treatment that targets H3K4me3 to increase chemotherapeutic efficacy in pediatric PF-EPN-A tumors.

Stabilization of HIF-1α and HIF-2α, up-regulation of MYCC and accumulation of stabilized p53 constitute hallmarks of CNS-PNET animal model

Recently, we described a new animal model of CNS primitive neuroectodermal tumors (CNS-PNET), which was generated by orthotopic transplantation of human Radial Glial (RG) cells into NOD-SCID mice's brain sub-ventricular zone. In the current study we conducted comprehensive RNA-Seq analyses to gain insights on the mechanisms underlying tumorigenesis in this mouse model of CNS-PNET. Here we show that the RNA-Seq profiles derived from these tumors cluster with those reported for patients' PNETs. Moreover, we found that (i) stabilization of HIF-1α and HIF-2α, which are involved in mediation of the hypoxic responses in the majority of cell types, (ii) up-regulation of MYCC, a key onco-protein whose dysregulation occurs in ~70% of human tumors, and (iii) accumulation of stabilized p53, which is commonly altered in human cancers, constitute hallmarks of our tumor model, and might represent the basis for CNS-PNET tumorigenesis in this model. We discuss the possibility that these three events might be interconnected. These results indicate that our model may prove invaluable to uncover the molecular events leading to MYCC and TP53 alterations, which would be of broader interest considering their relevance to many human malignancies. Lastly, this mouse model might prove useful for drug screening targeting MYCC and related members of its protein interaction network.

Expression of miR-18a and miR-210 in Normal Breast Tissue as Candidate Biomarkers of Breast Cancer Risk

miRNAs are noncoding RNAs with abnormal expression in breast cancer; their expression in high-risk benign breast tissue may relate to breast cancer risk. We examined miRNA profiles in contralateral unaffected breasts (CUB) of patients with breast cancer and validated resulting candidates in two additional sample sets. Expression profiles of 754 mature miRNAs were examined using TaqMan Low Density Arrays in 30 breast cancer samples [15 estrogen receptor (ER)-positive and 15 ER-negative] and paired CUBs and 15 reduction mammoplasty controls. Pairwise comparisons identified miRNAs with significantly differential expression. Seven candidate miRNAs were examined using qRT-PCR in a second CUB sample set (40 cases, 20 ER+, 20 ER-) and 20 reduction mammoplasty controls. Further validation was performed in 80 benign breast biopsy (BBB) samples; 40 from cases who subsequently developed breast cancer and 40 from controls who did not. Logistic regression, using tertiles of miRNA expression, was used to discriminate cases from controls. Seven miRNAs were differentially expressed in tumors and CUBs versus reduction mammoplasty samples. Among them, miR-18a and miR-210 were validated in the second CUB set, showing significantly higher expression in tumor and CUBs than in reduction mammoplasty controls. The expression of miR-18a and miR-210 was also significantly higher in BBB cases than in BBB controls. When both miR-18a and miR-210 were expressed in the upper tertiles in BBB, OR for subsequent cancer was 3.20, P = 0.023. miR-18a and miR-210 are expressed at higher levels in CUBs of patients with breast cancer, and in BBB prior to cancer development, and are therefore candidate breast cancer risk biomarkers. Cancer Prev Res; 10(1); 89-97. ©2016 AACR.

The P4 Health Spectrum - A Predictive, Preventive, Personalized and Participatory Continuum for Promoting Healthspan

Chronic diseases (i.e., noncommunicable diseases), mainly cardiovascular disease, cancer, respiratory diseases and type-2-diabetes, are now the leading cause of death, disability and diminished quality of life on the planet. Moreover, these diseases are also a major financial burden worldwide, significantly impacting the economy of many countries. Healthcare systems and medicine have progressively improved upon the ability to address infectious diseases and react to adverse health events through both surgical interventions and pharmacology; we have become efficient in delivering reactive care (i.e., initiating interventions once an individual is on the verge of or has actually suffered a negative health event). However, with slowly progressing and often 'silent' chronic diseases now being the main cause of illness, healthcare and medicine must evolve into a proactive system, moving away from a merely reactive approach to care. Minimal interactions among the specialists and limited information to the general practitioner and to the individual receiving care lead to a fragmented health approach, non-concerted prescriptions, a scattered follow-up and a suboptimal cost-effectiveness ratio. A new approach in medicine that is predictive, preventive, personalized and participatory, which we label here as "P4" holds great promise to reduce the burden of chronic diseases by harnessing technology and an increasingly better understanding of environment-biology interactions, evidence-based interventions and the underlying mechanisms of chronic diseases. In this concept paper, we propose a 'P4 Health Continuum' model as a framework to promote and facilitate multi-stakeholder collaboration with an orchestrated common language and an integrated care model to increase the healthspan.

MicroRNA-10a Regulation of Proinflammatory Mediators: An Important Component of Untreated Juvenile Dermatomyositis

Objective.

To identify differentially expressed microRNA (miRNA) in muscle biopsies (MBx) from 15 untreated children with juvenile dermatomyositis (JDM) compared with 5 controls.

Methods.

Following MBx miRNA profiling, differentially expressed miRNA and their protein targets were validated by quantitative real-time PCR (qRT-PCR) and immunological assay. The association of miRNA-10a and miRNA-10b with clinical data was evaluated, including Disease Activity Score (DAS), von Willebrand factor antigen (vWF:Ag), nailfold capillary end row loops, duration of untreated disease, and tumor necrosis factor (TNF)-α-308A allele.

Results.

In JDM, 16/362 miRNA were significantly differentially expressed [false discovery rate (FDR) < 0.05]. Among these, miRNA-10a was the most downregulated miRNA in both FDR and ranking of fold change: miRNA-10a = −2.27-fold, miRNA-10b = −1.80-fold. Decreased miRNA-10a and miRNA-10b expressions were confirmed using qRT-PCR: −4.16 and −2.59 fold, respectively. The qRT-PCR documented that decreased miRNA-10a expression was related to increased vascular cell adhesion molecule 1 in 13 of these JDM cases (correlation −0.67, p = 0.012), unlike miRNA-10b data (not significant). Concurrent JDM plasma contained increased levels of interleukin (IL) 6 (p = 0.0363), IL-8 (p = 0.0005), TNF-α (p = 0.0011), and monocyte chemoattractant proteins 1 (p = 0.0139). Decreased miRNA-10a, but not miRNA-10b, was associated with the TNF-α-308A allele (p = 0.015). In the 15 JDM, a trend of association of miRNA-10a (but not miRNA-10b) with vWF:Ag and DAS was observed.

Conclusion.

MiRNA-10a downregulation is an important element in untreated JDM muscle pathophysiology. We speculate that muscle miRNA expression in adult dermatomyositis differs from muscle miRNA expression in untreated childhood JDM.

CD84 is markedly up-regulated in Kawasaki disease arteriopathy

The major goals of Kawasaki disease (KD) therapy are to reduce inflammation and prevent thrombosis in the coronary arteries (CA), but some children do not respond to currently available non-specific therapies. New treatments have been difficult to develop because the molecular pathogenesis is unknown. In order to identify dysregulated gene expression in KD CA, we performed high-throughput RNA sequencing on KD and control CA, validated potentially dysregulated genes by real-time reverse transcription-polymerase chain reaction (RT-PCR) and localized protein expression by immunohistochemistry. Signalling lymphocyte activation molecule CD84 was up-regulated 16-fold (P < 0·01) in acute KD CA (within 2 months of onset) and 32-fold (P < 0·01) in chronic CA (5 months to years after onset). CD84 was localized to inflammatory cells in KD tissues. Genes associated with cellular proliferation, motility and survival were also up-regulated in KD CA, and immune activation molecules MX2 and SP140 were up-regulated in chronic KD. CD84, which facilitates immune responses and stabilizes platelet aggregates, is markedly up-regulated in KD CA in patients with acute and chronic arterial disease. We provide the first molecular evidence of dysregulated inflammatory responses persisting for months to years in CA significantly damaged by KD.

Human rhinovirus infection causes different DNA methylation changes in nasal epithelial cells from healthy and asthmatic subjects

Background

Mechanisms underlying the development of virus-induced asthma exacerbations remain unclear. To investigate if epigenetic mechanisms could be involved in virus-induced asthma exacerbations, we undertook DNA methylation profiling in asthmatic and healthy nasal epithelial cells (NECs) during Human Rhinovirus (HRV) infection in vitro.

Methods

Global and loci-specific methylation profiles were determined via Alu element and Infinium Human Methylation 450 K microarray, respectively. Principal components analysis identified the genomic loci influenced the most by disease-status and infection. Real-time PCR and pyrosequencing were used to confirm gene expression and DNA methylation, respectively.

Results

HRV infection significantly increased global DNA methylation in cells from asthmatic subjects only (43.6% to 44.1%, p = 0.04). Microarray analysis revealed 389 differentially methylated loci either based on disease status, or caused by virus infection. There were disease-associated DNA methylation patterns that were not affected by HRV infection as well as HRV-induced DNA methylation changes that were unique to each group. A common methylation locus stood out in response to HRV infection in both groups, where the small nucleolar RNA, H/ACA box 12 (SNORA12) is located. Further analysis indicated that a relationship existed between SNORA12 DNA methylation and gene expression in response to HRV infection.

Conclusions

We describe for the first time that Human rhinovirus infection causes DNA methylation changes in airway epithelial cells that differ between asthmatic and healthy subjects. These epigenetic differences may possibly explain the mechanism by which respiratory viruses cause asthma exacerbations.

Histone deacetylase 1 deficiency impairs differentiation and electrophysiological properties of cardiomyocytes derived from induced pluripotent cells

Epigenetic and chromatin modifications play particularly important roles in embryonic and induced pluripotent stem cells (ESCs and iPSCs) allowing for the cells to both differentiate and dedifferentiate back to a pluripotent state. We analyzed how the loss of a key chromatin-modifying enzyme, histone deacetylase 1 (HDAC1), affects early and cardiovascular differentiation of both ESCs and iPSCs. We also investigated potential differences between these two cell types when differentiation is induced. Our data indicate an essential role for HDAC1 in deacetylating regulatory regions of key pluripotency-associated genes during early differentiation. Although HDAC1 functions primarily as a HDAC, its loss also affects DNA methylation in ESCs and iPSCs both during pluripotency and differentiation. We show that HDAC1 plays a crucial, nonredundant role in cardiomyocyte differentiation and maturation. Our data also elucidate important differences between ESCs and iPSCs, when levels of this enzyme are reduced, that affect their ability to differentiate into functional cardiomyocytes. As varying levels of chromatin-modifying enzymes are likely to exist in patient-derived iPSCs, understanding the molecular circuitry of these enzymes in ESCs and iPSCs is critical for their potential use in cardiovascular therapeutic applications

Differentially expressed miRNAs in retinoblastoma

MicroRNAs (miRNAs) are short non-coding RNA transcripts that have the ability to regulate the expression of target genes, and have been shown to influence the development of various tumors. The purpose of our study is to identify aberrantly expressed miRNAs in retinoblastoma for the discovery of potential therapeutic targets for this disease, and to gain a greater understanding of the mechanisms driving retinoblastoma progression. We report 41 differentially expressed miRNAs (p<0.05) in 12 retinoblastomas as compared to three normal human retinae. Of these miRNAs, many are newly identified as being differentially expressed in retinoblastoma. Further, we report the validations of five of the most downregulated miRNAs in primary human retinoblastomas (p<0.05), human retinoblastoma cell lines, and mouse retinoblastoma cell lines. This serves as the largest and most comprehensive retinoblastoma miRNA analysis to date with corresponding clinical and pathological characteristics. This is an essential step in the discovery of miRNAs associated with retinoblastoma progression, and in the identification of potential therapeutic targets for this disease.

ANALYSIS OF ALEXANDRIUM TAMARENSE (DINOPHYCEAE) GENES REVEALS THE COMPLEX EVOLUTIONARY HISTORY OF A MICROBIAL EUKARYOTE()

Microbial eukaryotes may extinguish much of their nuclear phylogenetic history due to endosymbiotic/horizontal gene transfer (E/HGT). We studied E/HGT in 32,110 contigs of expressed sequence tags (ESTs) from the dinoflagellate Alexandrium tamarense (Dinophyceae) using a conservative phylogenomic approach. The vast majority of predicted proteins (86.4%) in this alga are novel or dinoflagellate-specific. We searched for putative homologs of these predicted proteins against a taxonomically broadly sampled protein database that includes all currently available data from algae and protists and reconstructed a phylogeny from each of the putative homologous protein sets. Of the 2,523 resulting phylogenies, 14-17% are potentially impacted by E/HGT involving both prokaryote and eukaryote lineages, with 2-4% showing clear evidence of reticulate evolution. The complex evolutionary histories of the remaining proteins, many of which may also have been affected by E/HGT, cannot be interpreted using our approach with currently available gene data. We present empirical evidence of reticulate genome evolution that combined with inadequate or highly complex phylogenetic signal in many proteins may impede genome-wide approaches to infer the tree of microbial eukaryotes.

TGF-β mediated DNA methylation in prostate cancer

Almost all tumors harbor a defective negative feedback loop of signaling by transforming growth factor-β (TGF-β). Epigenetic mechanisms of gene regulation, including DNA methylation, are fundamental to normal cellular function and also play a major role in carcinogenesis. Recent evidence demonstrated that TGF-β signaling mediates cancer development and progression. Many key events in TGF-β signaling in cancer included auto-induction of TGF-β1 and increased expression of DNA methyltransferases (DNMTs), suggesting that DNA methylation plays a significant role in cancer development and progression. In this review, we performed an extensive survey of the literature linking TGF-β signaling to DNA methylation in prostate cancer. It appeared that almost all DNA methylated genes detected in prostate cancer are directly or indirectly related to TGF-β signaling. This knowledge has provided a basis for our future directions of prostate cancer research and strategies for prevention and therapy for prostate cancer.

Epigenetically reprogramming metastatic tumor cells with an embryonic microenvironment

We have previously shown that the microenvironment of human embryonic stem cells (hESCs) is able to change and reprogram aggressive cancer cells to a less aggressive state. Some mechanisms implicated in the phenotypic changes observed after this exposure are mainly associated with the Nodal signaling pathway, which plays a key role in tumor cell plasticity. However, several other molecular mechanisms might be related directly and/or indirectly to these changes, including microRNA (miRNA) regulation and DNA methylation.

AIM

To further explore the epigenetic mechanisms potentially underlying the phenotypic changes that occur after exposing metastatic melanoma cells to a hESC microenvironment.

MATERIALS & METHODS

A total of 365 miRNAs were screened using the TaqMan® Low Density Arrays. We also evaluated whether DNA methylation could be one of the factors regulating the expression of the inhibitor of Nodal, Lefty, in hESCs (where it is highly expressed) vs melanoma cells (where it is not expressed).

RESULTS

Using these experimental approaches, we identified miRNAs that are up- and down-regulated in melanoma cells exposed to a hESC microenvironment, such as miR-302a and miR-27b, respectively. We also demonstrate that Notch4 is one of the targets of miR-302a, which is upstream of Nodal. Additionally, one of the mechanisms that might explain the absence of the inhibitor of Nodal, Lefty, in cancer cells is silencing by DNA methylation, which provides new insights into the unregulated expression of Nodal in melanoma.

CONCLUSION

These findings suggest that epigenetic changes such as DNA methylation and regulation by microRNAs might play a significant role in tumor cell plasticity and the metastatic phenotype.

Swarm rat chondrosarcoma cells as an in vivo model: lung colonization and effects of tissue environment on tumor growth

Swarm rat chondrosarcoma cells have been used extensively for biochemical studies of extra-cellular matrix metabolism in cartilage. However, these cells also possess tumor-like behavior in vivo and are useful in investigation of chondrosarcoma biology. the current study was designed to develop a metastatic model using swarm rat chondrosarcoma cells, and to assess the effect of tissue-environment on tumor behavior in vivo. Tumors were implanted subcutaneously or into bone, and animals were assessed radiographically and microscopically for tumor growth and metastasis. The subcutaneous tumor grew to an average mass of 35 g, while tumor implanted into bone grew 75 mg. Transplantation of the cells into the bone led to extensive bone remodeling with invasion of the medullary cavity and destruction of the bone cortex. Light microscopy demonstrated no significant differences in the number of mitoses, cellular atypia or extracellular matrix staining between the two sites of tumor implantation. Interestingly, lung colonization was observed in none of the animals in the subcutaneous tumor injection group, while tumors colonized the lungs in 95% of the rats with tumor injected into bone. Analysis of cDNA libraries from subcutaneous and bone-transplanted tumors demonstrated a complex and diverse array of expressed transcripts, and there were significant differences in gene expression between tumors at different sites. The results of this study suggest swarm rat chondrosarcoma is a model that resembles human chondrosarcoma mimicking its ability to infiltrate and remodel local bone and to colonize the lungs. Furthermore, the interaction between host-tissue and tumor cells plays a major role in the tumor behavior in this model. Identifying these interactions will lead to further understanding of chondrosarcoma and contribute to therapeutic targets in the future.

Genetic and epigenetic variations contributed by Alu retrotransposition

Background

De novo retrotransposition of Alu elements has been recognized as a major driver for insertion polymorphisms in human populations. In this study, we exploited Alu-anchored bisulfite PCR libraries to identify evolutionarily recent Alu element insertions, and to investigate their genetic and epigenetic variation.

Results

A total of 327 putatively recent Alu insertions were identified, altogether represented by 1,762 sequence reads. Nearly all such de novo retrotransposition events (316/327) were novel. Forty-seven out of forty-nine randomly selected events, corresponding to nineteen genomic loci, were sequence-verified. Alu element insertions remained hemizygous in one or more individuals in sixteen of the nineteen genomic loci. The Alu elements were found to be enriched for young Alu families with characteristic sequence features, such as the presence of a longer poly(A) tail. In addition, we documented the occurrence of a duplication of the AT-rich target site in their immediate flanking sequences, a hallmark of retrotransposition. Furthermore, we found the sequence motif (TT/AAAA) that is recognized by the ORF2P protein encoded by LINE-1 in their 5'-flanking regions, consistent with the fact that Alu retrotransposition is facilitated by LINE-1 elements. While most of these Alu elements were heavily methylated, we identified an Alu localized 1.5 kb downstream of TOMM5 that exhibited a completely unmethylated left arm. Interestingly, we observed differential methylation of its immediate 5' and 3' flanking CpG dinucleotides, in concordance with the unmethylated and methylated statuses of its internal 5' and 3' sequences, respectively. Importantly, TOMM5's CpG island and the 3 Alu repeats and 1 MIR element localized upstream of this newly inserted Alu were also found to be unmethylated. Methylation analyses of two additional genomic loci revealed no methylation differences in CpG dinucleotides flanking the Alu insertion sites in the two homologous chromosomes, irrespective of the presence or absence of the insertion.

Conclusions

We anticipate that the combination of methodologies utilized in this study, which included repeat-anchored bisulfite PCR sequencing and the computational analysis pipeline herein reported, will prove invaluable for the generation of genetic and epigenetic variation maps.

Identification of microRNAs as potential prognostic markers in ependymoma

Introduction

We have examined expression of microRNAs (miRNAs) in ependymomas to identify molecular markers of value for clinical management. miRNAs are non-coding RNAs that can block mRNA translation and affect mRNA stability. Changes in the expression of miRNAs have been correlated with many human cancers.

Materials and Methods

We have utilized TaqMan Low Density Arrays to evaluate the expression of 365 miRNAs in ependymomas and normal brain tissue. We first demonstrated the similarity of expression profiles of paired frozen tissue (FT) and paraffin-embedded specimens (FFPE). We compared the miRNA expression profiles of 34 FFPE ependymoma samples with 8 microdissected normal brain tissue specimens enriched for ependymal cells. miRNA expression profiles were then correlated with tumor location, histology and other clinicopathological features.

Results

We have identified miRNAs that are over-expressed in ependymomas, such as miR-135a and miR-17-5p, and down-regulated, such as miR-383 and miR-485-5p. We have also uncovered associations between expression of specific miRNAs which portend a worse prognosis. For example, we have identified a cluster of miRNAs on human chromosome 14q32 that is associated with time to relapse. We also found that miR-203 is an independent marker for relapse compared to the parameters that are currently used. Additionally, we have identified three miRNAs (let-7d, miR-596 and miR-367) that strongly correlate to overall survival.

Conclusion

We have identified miRNAs that are differentially expressed in ependymomas compared with normal ependymal tissue. We have also uncovered significant associations of miRNAs with clinical behavior. This is the first report of clinically relevant miRNAs in ependymomas.

Efficacy of interstitial continuous vincristine infusion in a bioluminescent rodent intracranial tumor model

Interstitial chemotherapeutic drug infusion can bypass the blood-brain barrier, and provide high regional drug concentrations without systemic exposure. However, toxicity and efficacy for drugs administered via interstitial continuous (i.c.) infusion have not been characterized. In the current study, vincristine (VIN) was infused into the right frontal lobes of healthy Fisher 344 rats at 30, 45, 60, and 120 μg/ml over a period of 7 days at 1 μl/h, using an Alzet osmotic pump to evaluate toxicity. C6 rat glioblastoma cells transduced with a luciferase gene were inoculated into the right frontal lobe of a second group of rats. VIN was administered to tumor bearing rats via i.c. infusion 7 days later and tumor growth was monitored by bioluminescence intensity (BLI) to assess VIN efficacy, intravenous (i.v.) drug administration was used as a comparison drug delivery method. The results suggested that VIN toxicity is dose-dependent. Efficacy studies showed increased BLI, which correlates with histopathological tumor size, in saline-infused and i.v.-treated tumor-bearing rats. These rats survived an average of 28 ± 0.85 days and 33 ± 1.38 days, respectively. Both groups had large tumors at the time of death. Animals treated with VIN via i.c. infusion survived until day 90, the observation endpoint for this study. This was significantly longer than average survival times in the previous two groups. These results demonstrate that VIN via i.c. infusion is effective in reducing C6 glioblastoma tumors and prolonging rodent survival time compared to i.v. injection and suggest that chemotherapeutic drug administration via i.c. infusion may be a promising strategy for treating malignant brain tumors.

Abstract 3009: Genetic and epigenetic variations resulting from Alu retrotransposition

Alu repeats are the most prevalent retrotransposons, comprising approximately 10% of the human genome. Alu retrotransposons have amplified during primate evolution by an RNA-mediated copy and paste mechanism. Among the Alu subfamilies, AluY elements and its variants exhibit the highest rates of retrotransposition. De novo retrotransposition of active Alu elements have given rise to insertion polymorphisms in human populations. While generating methylation maps of Alu elements in human normal cerebellum and in ependymoma tumors, Xie and colleagues (2009, 2010) also identified putative new Alu element insertions. In this study, we explored these Alu libraries to identify recent Alu element insertions and to understand the (epi)genetic variations that result from such events. This study may also enable identification of DNA structural variation that may occur at the onset of tumorigenesis and/or during tumor progression.

A total of 327 putative recent Alu element insertions represented by 1,762 sequence reads were identified, including 316 events that have not yet been documented in either dbSNP or dbRIP. Forty-seven out of a total of forty-nine randomly selected events representing nineteen genomic loci were sequence-verified. Alu element insertions remained heterozygous for sixteen out of nineteen genomic loci in one or more individuals. The polymorphic Alu elements are enriched for young Alu families and they present some characteristic sequence features, such as longer poly(A) tails. In addition, they also present the TT/AAAA consensus sequence at their 3’-flanking regions, and AT-rich target site duplications. Despite being homozygous or heterozygous insertions, these Alu elements are heavily methylated. Further methylation analysis on two genomic loci revealed no methylation difference in two heterozygous alleles and also no methylation difference at CpG dinucleotides flanking the Alu insertion sites. These results indicate that the recently integrated Alu elements have not influenced the methylation status of their neighboring CpG dinucleotides.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 3009. doi:10.1158/1538-7445.AM2011-3009

Abstract 168: Differential microRNA expression profiles of breast cancer and matched benign breast samples stratified by ER status

Background: The identification of women at high risk for breast cancer remains a challenge. Robust risk markers that better stratify women for their probability of developing breast cancer would represent a significant advance for breast cancer prevention research. microRNAs (miRNAs) are a large group of non-coding RNAs, generally 18-25 nt long, which block messenger RNA (mRNA) translation and affect mRNA stability. Since the molecular signature of ER-positive breast cancer differs from that of ER-negative breast cancer, we hypothesized that miRNA profile differences may be reflected in the contralateral normal breast tissue of women with breast cancer.Methods:we used a well-annotated set of samples comprised of matched contra-lateral benign breast tissue (MCBB) from 30 women (15 ER positive and 15 ER negative tumor with their corresponding MCBBs matched for age and race) to generate expression profiles of 754 mature miRNAs by TaqMan Low Density Arrays (TLDAs). Following background subtraction, miRNA levels were normalized against controls and compared across groups using tests of contrast on estimated parameters from a full ANOVA model with Empirical Bayes variance correction and False Discovery Rate (FDR) control at p≪0.05 using the Benjamini-Hochberg method. Following this strategy, we have performed different pair-wise comparisons as follows: 1) ER+ to ER- tumors, 2) ER+ tumor to MCBB tissue, 3) ER- tumor to MCBB, and 4) MCBB tissue from women with ER+ versus ER- tumors.Results: We observe clear patterns of differential expression between ER+ and ER- tumors. Four miRNAs were significantly down-regulated in ER+ tumors compared to ER-negative, whereas miR190b was significantly up-regulated in ER- lesions. More striking were the patterns of differential expression of 23 miRNAs in ER- breast cancer tissue compared to MCCB. A similar analysis in ER+ breast cancer and MCBB yeilded at least 29 miRNAs with significantly differential expression. Among these differentially expressed miRNAs, nine miRNAs were common in the ER- and ER+ cancers, and the remainder were distinct. Notably, we found that miR-18a and miR-18b expression were more than 50-fold higher in ER- cancer samples compared to MCBB, but when we compared ER+ tumor to MCBB, we found no significant difference in miR-18b expression. On the other hand, miR-96 and miR-183 are over-expressed in both ER+ and ER- tumors when compared to the MCBB set. Comparing MCBB from women with ER-and ER+ breast cancer yielded no significant differences. Conclusion: we see striking differences in the patterns of miRNA expression differences between breast cancer samples and matched contralateral benign tissue when stratified by ER status.The majority of the miRNAs identified in our analyses are novel. These findings have significant potential for the identification of new ER- status specific and biologically significant molecular biomarkers.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 168. doi:10.1158/1538-7445.AM2011-168

Epigenetic reprogramming as a key contributor to melanocyte malignant transformation

Melanoma progression requires deregulation of gene expression by currently uncharacterized epigenetic mechanisms. A mouse model based on changes in cell microenvironment was developed by our group to study melanocyte malignant transformation. Melanoma cell lines (4C11- and 4C11+) were obtained as result of 5 sequential anchorage blockades of non-tumorigenic melan-a melanocytes. Melan-a cells submitted to 4 de-adhesion cycles were also established (4C), are non-tumorigenic and represent an intermediary phase of tumor progression. The aim of this work was to identify factors contributing to epigenetic modifications in early and later phases of malignant transformation induced by anchorage impediment. Epigenetic alterations occur early in tumorigenesis; 4C cell line shows changes in global and gene-specific DNA methylation and histone marks. Many histone modifications differ between melan-a, 4C, 4C11- (non-metastatic melanoma cell line) and 4C11+ (metastatic melanoma cell line) which could be associated with changes in gene and microRNA expression. These epigenetic alterations seem to play a key role in malignant transformation since melanocytes treated with 5-Aza-2'-deoxycytidine before each anchorage blockade do not transform. Some epigenetic changes seem to be also responsible for the maintenance of malignant phenotype, since melanoma cell lines (4C11- and 4C11+) treated in vitro with 5-Aza-2'-deoxycytidine or Trichostatin A showed reduction of tumor growth in vivo. Changes in gene expression reflecting cell adaptation to new environment were also observed. We propose a model in which sustained microenvironmental stress in melanocytes results in epigenetic reprogramming. Thus, after adaptation, cells may acquire epigenetic marks that could contribute to the establishment of a malignant phenotype.

Genome-wide quantitative assessment of variation in DNA methylation patterns

Genomic DNA methylation contributes substantively to transcriptional regulations that underlie mammalian development and cellular differentiation. Much effort has been made to decipher the molecular mechanisms governing the establishment and maintenance of DNA methylation patterns. However, little is known about genome-wide variation of DNA methylation patterns. In this study, we introduced the concept of methylation entropy, a measure of the randomness of DNA methylation patterns in a cell population, and exploited it to assess the variability in DNA methylation patterns of Alu repeats and promoters. A few interesting observations were made: (i) within a cell population, methylation entropy varies among genomic loci; (ii) among cell populations, the methylation entropies of most genomic loci remain constant; (iii) compared to normal tissue controls, some tumors exhibit greater methylation entropies; (iv) Alu elements with high methylation entropy are associated with high GC content but depletion of CpG dinucleotides and (v) Alu elements in the intronic regions or far from CpG islands are associated with low methylation entropy. We further identified 12 putative allelic-specific methylated genomic loci, including four Alu elements and eight promoters. Lastly, using subcloned normal fibroblast cells, we demonstrated the highly variable methylation patterns are resulted from low fidelity of DNA methylation inheritance.

Cancer hallmarks in induced pluripotent cells: new insights

Studies are beginning to emerge that demonstrate intriguing differences between human-induced pluripotent stem cells (hiPSCs) and human embryonic stem cells (hESCs). Here, we investigated the expression of key members of the Nodal embryonic signaling pathway, critical to the maintenance of pluripotency in hESCs. Western blot and real-time RT-PCR analyses reveal slightly lower levels of Nodal (a TGF-beta family member) and Cripto-1 (Nodal's co-receptor) and a dramatic decrease in Lefty (Nodal's inhibitor and TGF-beta family member) in hiPSCs compared with hESCs. The noteworthy drop in hiPSC's Lefty expression correlated with an increase in the methylation of Lefty B CpG island. Based on these findings, we addressed a more fundamental question related to the consequences of epigenetically reprogramming hiPSCs, especially with respect to maintaining a stable ESC phenotype. A global comparative analysis of 365 microRNAs (miRs) in two hiPSC versus four hESC lines ultimately identified 10 highly expressed miRs in hiPCSs with >10-fold difference, which have been shown to be cancer related. These data demonstrate cancer hallmarks expressed by hiPSCs, which will require further assessment for their impact on future therapies..

Targeting an embryonic signaling pathway to suppress the metastatic phenotype

The diagnosis of melanoma is becoming more frequent. Although surgical excision of early lesions is associated with relatively significant high cure rates, treatment modalities are largely unsuccessful for advanced disease. Characteristics such as cellular heterogeneity, phenotypic plasticity and the expression of embryonic related molecules and signaling pathways that are important for self renewal and pluripotency, present a challenge in targeting the most aggressive tumor cells. However, the absence of major regulatory checkpoints in these tumor cells allows aberrant activation of embryonic signaling pathways. In particular, the embryonic morphogen Nodal, belonging to the TGF-β superfamily, is an important regulator of embryonic stem cell fate. We have recently demonstrated that Nodal is also expressed in aggressive melanoma, although unlike embryonic stem cells, melanoma does not express Lefty, a natural inhibitor of Nodal, thus allowing Nodal to be expressed in an unregulated manner. Further analysis revealed that Lefty is silenced in these tumor cells by DNA methylation. Down-regulation of Nodal expression in melanoma cells results in a significant reduction in clonogenicity and suppression of tumorigenesis. Additional translational studies indicate that antibodies against Nodal are capable of targeting and inducing apoptosis in human melanoma tumors lodged in the lungs of nude mice. Furthermore, our data show another stem cell signaling pathway, Notch, is expressed in these aggressive melanoma cells, and there is direct evidence demonstrating the molecular cross-talk between Nodal and Notch. Thus, these two prominent stem cell pathways contributing to melanoma cell plasticity may serve as promising dual therapeutic targets for clinical intervention.

Folic acid remodels chromatin on Hes1 and Neurog2 promoters during caudal neural tube development

The mechanism(s) behind folate rescue of neural tube closure are not well understood. In this study we show that maternal intake of folate prior to conception reverses the proliferation potential of neural crest stem cells in homozygous Splotch embryos (Sp(-/-)) via epigenetic mechanisms. It is also shown that the pattern of differentiation seen in these cells is similar to wild-type (WT). Cells from open caudal neural tubes of Sp(-/-) embryos exhibit increased H3K27 methylation and decreased expression of KDM6B possibly due to up-regulation of KDM6B targeting micro-RNAs such as miR-138, miR-148a, miR-185, and miR-339-5p. In our model, folate reversed these epigenetic marks in folate-rescued Sp(-/-) embryos. Using tissue from caudal neural tubes of murine embryos we also examined H3K27me2 and KDM6B association with Hes1 and Neurog2 promoters at embryonic day E10.5, the proliferative stage, and E12.5, when neural differentiation begins. In Sp(-/-) embryos compared with WT, levels of H3K27me2 associated with the Hes1 promoter were increased at E10.5, and levels associated with the Neurog2 promoter were increased at E12.5. KDM6B association with Hes1 and Neurog2 promoters was inversely related to H3K27me2 levels. These epigenetic changes were reversed in folate-rescued Sp(-/-) embryos. Thus, one of the mechanisms by which folate may rescue the Sp(-/-) phenotype is by increasing the expression of KDM6B, which in turn decreases H3K27 methylation marks on Hes1 and Neurog2 promoters thereby affecting gene transcription.

Microenvironment alters epigenetic and gene expression profiles in Swarm rat chondrosarcoma tumors

BACKGROUND

Chondrosarcomas are malignant cartilage tumors that do not respond to traditional chemotherapy or radiation. The 5-year survival rate of histologic grade III chondrosarcoma is less than 30%. An animal model of chondrosarcoma has been established--namely, the Swarm Rat Chondrosarcoma (SRC)--and shown to resemble the human disease. Previous studies with this model revealed that tumor microenvironment could significantly influence chondrosarcoma malignancy.

METHODS

To examine the effect of the microenvironment, SRC tumors were initiated at different transplantation sites. Pyrosequencing assays were utilized to assess the DNA methylation of the tumors, and SAGE libraries were constructed and sequenced to determine the gene expression profiles of the tumors. Based on the gene expression analysis, subsequent functional assays were designed to determine the relevancy of the specific genes in the development and progression of the SRC.

RESULTS

The site of transplantation had a significant impact on the epigenetic and gene expression profiles of SRC tumors. Our analyses revealed that SRC tumors were hypomethylated compared to control tissue, and that tumors at each transplantation site had a unique expression profile. Subsequent functional analysis of differentially expressed genes, albeit preliminary, provided some insight into the role that thymosin-β4, c-fos, and CTGF may play in chondrosarcoma development and progression.

CONCLUSION

This report describes the first global molecular characterization of the SRC model, and it demonstrates that the tumor microenvironment can induce epigenetic alterations and changes in gene expression in the SRC tumors. We documented changes in gene expression that accompany changes in tumor phenotype, and these gene expression changes provide insight into the pathways that may play a role in the development and progression of chondrosarcoma. Furthermore, specific functional analysis indicates that thymosin-β4 may have a role in chondrosarcoma metastasis.

Transcriptional profiling of polycythemia vera identifies gene expression patterns both dependent and independent from the action of JAK2V617F

PURPOSE

To understand the changes in gene expression in polycythemia vera (PV) progenitor cells and their relationship to JAK2V617F.

EXPERIMENTAL DESIGN

Messenger RNA isolated from CD34(+) cells from nine PV patients and normal controls was profiled using Affymetrix arrays. Gene expression change mediated by JAK2V617F was determined by profiling CD34(+) cells transduced with the kinase and by analysis of leukemia cell lines harboring JAK2V617F, treated with an inhibitor.

RESULTS

A PV expression signature was enriched for genes involved in hematopoietic development, inflammatory responses, and cell proliferation. By quantitative reverse transcription-PCR, 23 genes were consistently deregulated in all patient samples. Several of these genes such as WT1 and KLF4 were regulated by JAK2, whereas others such as NFIB and EVI1 seemed to be deregulated in PV by a JAK2-independent mechanism. Using cell line models and comparing gene expression profiles of cell lines and PV CD34(+) PV specimens, we have identified panels of 14 JAK2-dependent genes and 12 JAK2-independent genes. These two 14- and 12-gene sets could separate not only PV from normal CD34(+) specimens, but also other MPN such as essential thrombocytosis and primary myelofibrosis from their normal counterparts.

CONCLUSIONS

A subset of the aberrant gene expression in PV progenitor cells can be attributed to the action of the mutant kinase, but there remain a significant number of genes characteristic of the disease but deregulated by as yet unknown mechanisms. Genes deregulated in PV as a result of the action of JAK2V617F or independent of the kinase may represent other targets for therapy.

Construction of a medicinal leech transcriptome database and its application to the identification of leech homologs of neural and innate immune genes

Background

The medicinal leech, Hirudo medicinalis, is an important model system for the study of nervous system structure, function, development, regeneration and repair. It is also a unique species in being presently approved for use in medical procedures, such as clearing of pooled blood following certain surgical procedures. It is a current, and potentially also future, source of medically useful molecular factors, such as anticoagulants and antibacterial peptides, which may have evolved as a result of its parasitizing large mammals, including humans. Despite the broad focus of research on this system, little has been done at the genomic or transcriptomic levels and there is a paucity of openly available sequence data. To begin to address this problem, we constructed whole embryo and adult central nervous system (CNS) EST libraries and created a clustered sequence database of the Hirudo transcriptome that is available to the scientific community.

Results

A total of ~133,000 EST clones from two directionally-cloned cDNA libraries, one constructed from mRNA derived from whole embryos at several developmental stages and the other from adult CNS cords, were sequenced in one or both directions by three different groups: Genoscope (French National Sequencing Center), the University of Iowa Sequencing Facility and the DOE Joint Genome Institute. These were assembled using the phrap software package into 31,232 unique contigs and singletons, with an average length of 827 nt. The assembled transcripts were then translated in all six frames and compared to proteins in NCBI's non-redundant (NR) and to the Gene Ontology (GO) protein sequence databases, resulting in 15,565 matches to 11,236 proteins in NR and 13,935 matches to 8,073 proteins in GO. Searching the database for transcripts of genes homologous to those thought to be involved in the innate immune responses of vertebrates and other invertebrates yielded a set of nearly one hundred evolutionarily conserved sequences, representing all known pathways involved in these important functions.

Conclusions

The sequences obtained for Hirudo transcripts represent the first major database of genes expressed in this important model system. Comparison of translated open reading frames (ORFs) with the other openly available leech datasets, the genome and transcriptome of Helobdella robusta, shows an average identity at the amino acid level of 58% in matched sequences. Interestingly, comparison with other available Lophotrochozoans shows similar high levels of amino acid identity, where sequences match, for example, 64% with Capitella capitata (a polychaete) and 56% with Aplysia californica (a mollusk), as well as 58% with Schistosoma mansoni (a platyhelminth). Phylogenetic comparisons of putative Hirudo innate immune response genes present within the Hirudo transcriptome database herein described show a strong resemblance to the corresponding mammalian genes, indicating that this important physiological response may have older origins than what has been previously proposed.

An insight into the sialome of Glossina morsitans morsitans

BACKGROUND

Blood feeding evolved independently in worms, arthropods and mammals. Among the adaptations to this peculiar diet, these animals developed an armament of salivary molecules that disarm their host's anti-bleeding defenses (hemostasis), inflammatory and immune reactions. Recent sialotranscriptome analyses (from the Greek sialo = saliva) of blood feeding insects and ticks have revealed that the saliva contains hundreds of polypeptides, many unique to their genus or family. Adult tsetse flies feed exclusively on vertebrate blood and are important vectors of human and animal diseases. Thus far, only limited information exists regarding the Glossina sialome, or any other fly belonging to the Hippoboscidae.

RESULTS

As part of the effort to sequence the genome of Glossina morsitans morsitans, several organ specific, high quality normalized cDNA libraries have been constructed, from which over 20,000 ESTs from an adult salivary gland library were sequenced. These ESTs have been assembled using previously described ESTs from the fat body and midgut libraries of the same fly, thus totaling 62,251 ESTs, which have been assembled into 16,743 clusters (8,506 of which had one or more EST from the salivary gland library). Coding sequences were obtained for 2,509 novel proteins, 1,792 of which had at least one EST expressed in the salivary glands. Despite library normalization, 59 transcripts were overrepresented in the salivary library indicating high levels of expression. This work presents a detailed analysis of the salivary protein families identified. Protein expression was confirmed by 2D gel electrophoresis, enzymatic digestion and mass spectrometry. Concurrently, an initial attempt to determine the immunogenic properties of selected salivary proteins was undertaken.

CONCLUSIONS

The sialome of G. m. morsitans contains over 250 proteins that are possibly associated with blood feeding. This set includes alleles of previously described gene products, reveals new evidence that several salivary proteins are multigenic and identifies at least seven new polypeptide families unique to Glossina. Most of these proteins have no known function and thus, provide a discovery platform for the identification of novel pharmacologically active compounds, innovative vector-based vaccine targets, and immunological markers of vector exposure.

Global demethylation of rat chondrosarcoma cells after treatment with 5-aza-2'-deoxycytidine results in increased tumorigenicity

Abnormal patterns of DNA methylation are observed in several types of human cancer. While localized DNA methylation of CpG islands has been associated with gene silencing, the effect that genome-wide loss of methylation has on tumorigenesis is not completely known. To examine its effect on tumorigenesis, we induced DNA demethylation in a rat model of human chondrosarcoma using 5-aza-2-deoxycytidine. Rat specific pyrosequencing assays were utilized to assess the methylation levels in both LINEs and satellite DNA sequences following 5-aza-2-deoxycytidine treatment. Loss of DNA methylation was accompanied by an increase in invasiveness of the rat chondrosarcoma cells, in vitro, as well as by an increase in tumor growth in vivo. Subsequent microarray analysis provided insight into the gene expression changes that result from 5-aza-2-deoxycytidine induced DNA demethylation. In particular, two genes that may function in tumorigenesis, sox-2 and midkine, were expressed at low levels in control cells but upon 5-aza-2-deoxycytidine treatment these genes became overexpressed. Promoter region DNA analysis revealed that these genes were methylated in control cells but became demethylated following 5-aza-2-deoxycytidine treatment. Following withdrawal of 5-aza-2-deoxycytidine, the rat chondrosarcoma cells reestablished global DNA methylation levels that were comparable to that of control cells. Concurrently, invasiveness of the rat chondrosarcoma cells, in vitro, decreased to a level indistinguishable to that of control cells. Taken together these experiments demonstrate that global DNA hypomethylation induced by 5-aza-2-deoxycytidine may promote specific aspects of tumorigenesis in rat chondrosarcoma cells.

Longitudinal assessment of regional directed delivery in a rodent malignant glioma model

OBJECT

Direct delivery of chemotherapeutic agents for the treatment of brain tumors is an area of focus in the development of therapeutic paradigms because this method of delivery circumvents the blood-brain barrier without causing adverse systemic side effects. Few studies have investigated longitudinal tumor response to this type of therapy. In this study, the authors examined the time course of tumor response to direct delivery of a chemotherapeutic agent in a rodent malignant glioma model.

METHODS

To visualize tumor response to chemotherapy, the authors used bioluminescence imaging in a rodent model. Rat 9L gliosarcoma cells expressing a luciferase gene were inoculated into adult male rat striata. Ten days following surgery the animals were randomly divided into 4 groups. Groups 1 and 2 received 20 and 40 microl carboplatin (1 mg/ml), respectively, via convection-enhanced delivery (CED); Group 3 received 60 mg/kg carboplatin intraperitoneally; and Group 4 received no treatment. Tumor growth was correlated with luminescence levels twice weekly.

RESULTS

Differential growth curves were observed for the 4 groups. Systemically treated rats showed decreasing photon flux emission at 15.0 + or - 4.7 days; rats treated with 20- or 40-microl CED showed decreased emissions at 4.0 + or - 2.0 and 3.2 + or - 1.3 days after treatment, respectively. Histopathologically, 6 of 12 CED-treated animals exhibited no residual tumor at the end point of the study.

CONCLUSIONS

Direct and systemic delivery of carboplatin was examined to determine how the method of drug delivery affects tumor growth. The present report is one of the first in vivo studies to examine the time course of tumor response to direct drug delivery. The results indicate that direct drug delivery may be a promising option for treating gliomas.

Multiple mechanisms influence regulation of the cystic fibrosis transmembrane conductance regulator gene promoter

The cystic fibrosis transmembrane conductance regulator (CFTR) gene is driven by a promoter that cannot alone account for the temporal and tissue-specific regulation of the gene. This has led to the search for additional regulatory elements that cooperate with the basal promoter to achieve coordinated expression. We previously identified two alternative upstream exons of the gene that were mutually exclusive of the first exon, and one of which showed temporal regulation in the human and sheep lung. We now demonstrate that this alternative splice product generates a stable protein, which initiates translation at an ATG in exon 4, and thus lacks the N terminus of CFTR. The other splice variant inhibits translation of the protein. In a search for the promoter used by the upstream exons, we identified a novel element that contributes to the activity of the basal CFTR promoter in airway epithelial cells, but does not function independently. Finally, we demonstrate that, in primary airway cells, skin fibroblasts, and both airway and intestinal cell lines, the CFTR promoter is unmethylated, irrespective of CFTR expression status. Thus, methylation is not the main cause of inactivation of CFTR transcription.

Monosomy of chromosome 10 associated with dysregulation of epidermal growth factor signaling in glioblastomas

CONTEXT

Glioblastomas--uniformly fatal brain tumors--often have both monosomy of chromosome 10 and gains of the epidermal growth factor receptor (EGFR) gene locus on chromosome 7, an association for which the mechanism is poorly understood.

OBJECTIVES

To assess whether coselection of EGFR gains on 7p12 and monosomy 10 in glioblastomas promotes tumorigenic epidermal growth factor (EGF) signaling through loss of the annexin A7 (ANXA7) gene on 10q21.1-q21.2 and whether ANXA7 acts as a tumor suppressor gene by regulating EGFR in glioblastomas.

DESIGN, SETTING, AND PATIENTS

Multidimensional analysis of gene, coding sequence, promoter methylation, messenger RNA (mRNA) transcript, protein data for ANXA7 (and EGFR), and clinical patient data profiles of 543 high-grade gliomas from US medical centers and The Cancer Genome Atlas pilot project (made public 2006-2008; and unpublished, tumors collected 2001-2008). Functional analyses using LN229 and U87 glioblastoma cells.

MAIN OUTCOME MEASURES

Associations among ANXA7 gene dosage, coding sequence, promoter methylation, mRNA transcript, and protein expression. Effect of ANXA7 haploinsufficiency on EGFR signaling and patient survival. Joint effects of loss of ANXA7 and gain of EGFR expression on tumorigenesis.

RESULTS

Heterozygous ANXA7 gene deletion is associated with significant loss of ANXA7 mRNA transcript expression (P = 1 x 10(-15); linear regression) and a reduction (mean [SEM]) of 91.5% (2.3%) of ANXA7 protein expression compared with ANXA7 wild-type glioblastomas (P = .004; unpaired t test). ANXA7 loss of function stabilizes the EGFR protein (72%-744% increase in EGFR protein abundance) and augments EGFR transforming signaling in glioblastoma cells. ANXA7 haploinsufficiency doubles tumorigenic potential of glioblastoma cells, and combined ANXA7 knockdown and EGFR overexpression promotes tumorigenicity synergistically. The heterozygous loss of ANXA7 in approximately 75% of glioblastomas in the The Cancer Genome Atlas plus infrequency of ANXA7 mutation (approximately 6% of tumors) indicates its role as a haploinsufficiency gene. ANXA7 mRNA transcript expression, dichotomized at the median, associates with patient survival in 191 glioblastomas (log-rank P = .008; hazard ratio [HR], 0.667; 95% confidence interval [CI], 0.493-0.902; 46.9 vs 74.8 deaths/100 person-years for high vs low ANXA7 mRNA expression) and with a separate group of 180 high-grade gliomas (log-rank P = .00003; HR, 0.476; 95% CI, 0.333-0.680; 21.8 vs 50.0 deaths/100 person-years for high vs low ANXA7 mRNA expression). Deletion of the ANXA7 gene associates with poor patient survival in 189 glioblastomas (log-rank P = .042; HR, 0.686; 95% CI, 0.476-0.989; 54.0 vs 80.1 deaths/100 person-years for wild-type ANXA7 vs ANXA7 deletion).

CONCLUSION

Haploinsufficiency of the tumor suppressor ANXA7 due to monosomy of chromosome 10 provides a clinically relevant mechanism to augment EGFR signaling in glioblastomas beyond that resulting from amplification of the EGFR gene.

Analysis of expressed sequence tags from the four main developmental stages of Trypanosoma congolense

Trypanosoma congolense is one of the most economically important pathogens of livestock in Africa. Culture-derived parasites of each of the three main insect stages of the T. congolense life cycle, i.e., the procyclic, epimastigote and metacyclic stages, and bloodstream stage parasites isolated from infected mice, were used to construct stage-specific cDNA libraries and expressed sequence tags (ESTs or cDNA clones) in each library were sequenced. Thirteen EST clusters encoding different variant surface glycoproteins (VSGs) were detected in the metacyclic library and 26 VSG EST clusters were found in the bloodstream library, 6 of which are shared by the metacyclic library. Rare VSG ESTs are present in the epimastigote library, and none were detected in the procyclic library. ESTs encoding enzymes that catalyze oxidative phosphorylation and amino acid metabolism are about twice as abundant in the procyclic and epimastigote stages as in the metacyclic and bloodstream stages. In contrast, ESTs encoding enzymes involved in glycolysis, the citric acid cycle and nucleotide metabolism are about the same in all four developmental stages. Cysteine proteases, kinases and phosphatases are the most abundant enzyme groups represented by the ESTs. All four libraries contain T. congolense-specific expressed sequences not present in the Trypanosoma brucei and Trypanosoma cruzi genomes. Normalized cDNA libraries were constructed from the metacyclic and bloodstream stages, and found to be further enriched for T. congolense-specific ESTs. Given that cultured T. congolense offers an experimental advantage over other African trypanosome species, these ESTs provide a basis for further investigation of the molecular properties of these four developmental stages, especially the epimastigote and metacyclic stages for which it is difficult to obtain large quantities of organisms. The T. congolense EST databases are available at: http://www.sanger.ac.uk/Projects/T_congolense/EST_index.shtml. The sequence data have been submitted to EMBL under the following accession numbers: FN263376-FN292969.

Interstitial continuous infusion therapy in a malignant glioma model in rats

PURPOSE

Local direct delivery of chemotherapeutic agents for the treatment of brain tumors is an area of focus in the development of new therapeutic paradigms. These techniques need improvement, especially in terms of drug retention in brain tissue.

MATERIALS AND METHODS

In this study, we used a rat glioma model to examine carboplatin distribution, as measured by platinum penetration, after delivery via interstitial continuous (i.c.) infusion. We also examined rat survival times in response to carboplatin and oxaliplatin. I.C. infusion, a modified version of convection-enhanced delivery (CED) for local drug delivery, uses low volume (1 microl per hour) continuous infusion directly into the tumor.

RESULTS

I.C. infusion produced a nearly 360-fold higher concentration of platinum in tumor tissue and significantly prolonged rodent survival time compared to intraperitoneal (i.p.) infusion.

CONCLUSIONS

We showed i.c. infusion allows for circumvention of the blood-brain barrier, focused drug distribution, and sustained drug delivery. This method could be a promising strategy for treating brain tumors.

High-throughput sequence-based epigenomic analysis of Alu repeats in human cerebellum

DNA methylation, the only known covalent modification of mammalian DNA, occurs primarily in CpG dinucleotides. 51% of CpGs in the human genome reside within repeats, and 25% within Alu elements. Despite that, no method has been reported for large-scale ascertainment of CpG methylation in repeats. Here we describe a sequencing-based strategy for parallel determination of the CpG-methylation status of thousands of Alu repeats, and a computation algorithm to design primers that enable their specific amplification from bisulfite converted genomic DNA. Using a single primer pair, we generated amplicons of high sequence complexity, and derived CpG-methylation data from 31 178 Alu elements and their 5′ flanking sequences, altogether representing over 4 Mb of a human cerebellum epigenome. The analysis of the Alu methylome revealed that the methylation level of Alu elements is high in the intronic and intergenic regions, but low in the regions close to transcription start sites. Several hypomethylated Alu elements were identified and their hypomethylated status verified by pyrosequencing. Interestingly, some Alu elements exhibited a strikingly tissue-specific pattern of methylation. We anticipate the amplicons herein described to prove invaluable as epigenome representations, to monitor epigenomic alterations during normal development, in aging and in diseases such as cancer.

An experimental brainstem tumor model using in vivo bioluminescence imaging in rat

PURPOSE

Currently, there is no conclusive treatment for brainstem tumor. To facilitate the development of new treatments, it is essential to establish predictive preclinical in vivo models in which therapeutic modalities can be evaluated. Although a few rodent models have been reported, there is no novel approach that can monitor tumor response qualitatively and quantitatively.

MATERIALS AND METHODS

Bioluminescence imaging was used to characterize a rat brainstem tumor model. In this model, 9L gliosarcoma cells, transduced with an onco-retroviral vector containing the luciferase coding sequence, were inoculated into Fisher 344 rats.

RESULT

Histopathological assessment showed successful cell implantation into the brainstem. There was a strong correlation between pathological tumor volume and luminescence strength. Longitudinal quantitative responses of the tumor after application of a therapeutic agent were also demonstrated.

CONCLUSION

This study demonstrates a robust rodent model with the ability to monitor brainstem tumor growth and response to chemotherapeutic agents.

Gene discovery in an invasive tephritid model pest species, the Mediterranean fruit fly, Ceratitis capitata

Background

The medfly, Ceratitis capitata, is a highly invasive agricultural pest that has become a model insect for the development of biological control programs. Despite research into the behavior and classical and population genetics of this organism, the quantity of sequence data available is limited. We have utilized an expressed sequence tag (EST) approach to obtain detailed information on transcriptome signatures that relate to a variety of physiological systems in the medfly; this information emphasizes on reproduction, sex determination, and chemosensory perception, since the study was based on normalized cDNA libraries from embryos and adult heads.

Results

A total of 21,253 high-quality ESTs were obtained from the embryo and head libraries. Clustering analyses performed separately for each library resulted in 5201 embryo and 6684 head transcripts. Considering an estimated 19% overlap in the transcriptomes of the two libraries, they represent about 9614 unique transcripts involved in a wide range of biological processes and molecular functions. Of particular interest are the sequences that share homology with Drosophila genes involved in sex determination, olfaction, and reproductive behavior. The medfly transformer2 (tra2) homolog was identified among the embryonic sequences, and its genomic organization and expression were characterized.

Conclusion

The sequences obtained in this study represent the first major dataset of expressed genes in a tephritid species of agricultural importance. This resource provides essential information to support the investigation of numerous questions regarding the biology of the medfly and other related species and also constitutes an invaluable tool for the annotation of complete genome sequences. Our study has revealed intriguing findings regarding the transcript regulation of tra2 and other sex determination genes, as well as insights into the comparative genomics of genes implicated in chemosensory reception and reproduction.

A pig multi-tissue normalised cDNA library: large-scale sequencing, cluster analysis and 9K micro-array resource generation

Background

Domestic animal breeding and product quality improvement require the control of reproduction, nutrition, health and welfare in these animals. It is thus necessary to improve our knowledge of the major physiological functions and their interactions. This would be greatly enhanced by the availability of expressed gene sequences in the databases and by cDNA arrays allowing the transcriptome analysis of any function.

The objective within the AGENAE French program was to initiate a high-throughput cDNA sequencing program of a 38-tissue normalised library and generate a diverse microarray for transcriptome analysis in pig species.

Results

We constructed a multi-tissue cDNA library, which was normalised and subtracted to reduce the redundancy of the clones. Expressed Sequence Tags were produced and 24449 high-quality sequences were released in EMBL database. The assembly of all the public ESTs (available through SIGENAE website) resulted in 40786 contigs and 54653 singletons. At least one Agenae sequence is present in 11969 contigs (12.5%) and in 9291 of the deeper-than-one-contigs (22.8%). Sequence analysis showed that both normalisation and subtraction processes were successful and that the initial tissue complexity was maintained in the final libraries. A 9K nylon cDNA microarray was produced and is available through CRB-GADIE. It will allow high sensitivity transcriptome analyses in pigs.

Conclusion

In the present work, a pig multi-tissue cDNA library was constructed and a 9K cDNA microarray designed. It contributes to the Expressed Sequence Tags pig data, and offers a valuable tool for transcriptome analysis.

Systematic approaches for incorporating control spots and data quality information to improve normalization of cDNA microarray data

BACKGROUND

Normalization and data quality control are two important aspects in microarray data analysis. Proper normalization and data quality control ensure that intensity ratios provide meaningful and accurate measurement of relative gene expression values. Control spots such as spikes and housekeeping genes with known concentrations in two channels are often used for calibrating experimental parameters. They provide valuable information about experimental variation which can be utilized for better normalization. They are also needed for proper normalization in cases that the most of the spots tend to change in one direction. In addition, it is desirable to include information on spot quality. Such information is available in a typical microarray data set, but is not fully utilized by existing normalization methods.

RESULTS

We propose two extensions of the two-way semi-linear model (TW-SLM) for appropriately combining control genes and spot quality information in normalization. The first extension (TW-SLMC) is designed to systematically incorporate control spots in a semi-parametric model to calibrate estimated normalization curves so that the relative fold changes of gene expressions are accurately estimated. Extrapolation is not required in this approach. The second extension (TW-SLMQ) is proposed to incorporate spot quality measure into normalization. This approach down-weights spots with lower quality scores in normalization. These two extensions can be used simultaneously for normalizing a data set. Two microarray data sets are used to demonstrate the proposed methods.

AVAILABILITY

An R based computing package is developed for the proposed methods and available from the corresponding authors.

Genome sequence of Aedes aegypti, a major arbovirus vector

We present a draft sequence of the genome of Aedes aegypti, the primary vector for yellow fever and dengue fever, which at approximately 1376 million base pairs is about 5 times the size of the genome of the malaria vector Anopheles gambiae. Nearly 50% of the Ae. aegypti genome consists of transposable elements. These contribute to a factor of approximately 4 to 6 increase in average gene length and in sizes of intergenic regions relative to An. gambiae and Drosophila melanogaster. Nonetheless, chromosomal synteny is generally maintained among all three insects, although conservation of orthologous gene order is higher (by a factor of approximately 2) between the mosquito species than between either of them and the fruit fly. An increase in genes encoding odorant binding, cytochrome P450, and cuticle domains relative to An. gambiae suggests that members of these protein families underpin some of the biological differences between the two mosquito species.

Life cycle transcriptome of the malaria mosquito Anopheles gambiae and comparison with the fruitfly Drosophila melanogaster

The African mosquito Anopheles gambiae is the major vector of human malaria. We report a genome-wide survey of mosquito gene expression profiles clustered temporally into developmental programs and spatially into adult tissue-specific patterns. Global expression analysis shows that genes that belong to related functional categories or that encode the same or functionally linked protein domains are associated with characteristic developmental programs or tissue patterns. Comparative analysis of our data together with data published from Drosophila melanogaster reveal an overall strong and positive correlation of developmental expression between orthologous genes. The degree of correlation varies, depending on association of orthologs with certain developmental programs or functional groups. Interestingly, the similarity of gene expression is not correlated with the coding sequence similarity of orthologs, indicating that expression profiles and coding sequences evolve independently. In addition to providing a comprehensive view of temporal and spatial gene expression during the A. gambiae life cycle, this large-scale comparative transcriptomic analysis has detected important evolutionary features of insect transcriptomes.

Cyanobacterial contribution to algal nuclear genomes is primarily limited to plastid functions

A single cyanobacterial primary endosymbiosis that occurred approximately 1.5 billion years ago is believed to have given rise to the plastid in the common ancestor of the Plantae or Archaeplastida--the eukaryotic supergroup comprising red, green (including land plants), and glaucophyte algae. Critical to plastid establishment was the transfer of endosymbiont genes to the host nucleus (i.e., endosymbiotic gene transfer [EGT]). It has been postulated that plastid-derived EGT played a significant role in plant nuclear-genome evolution, with 18% (or 4,500) of all nuclear genes in Arabidopsis thaliana having a cyanobacterial origin with about one-half of these recruited for nonplastid functions. Here, we determine whether the level of cyanobacterial gene recruitment proposed for Arabidopsis is of the same magnitude in the algal sisters of plants by analyzing expressed-sequence tag (EST) data from the glaucophyte alga Cyanophora paradoxa. Bioinformatic analysis of 3,576 Cyanophora nuclear genes shows that 10.8% of these with significant database hits are of cyanobacterial origin and one-ninth of these have nonplastid functions. Our data indicate that unlike plants, early-diverging algal groups appear to retain a smaller number of endosymbiont genes in their nucleus, with only a minor proportion of these recruited for nonplastid functions.

Analysis of fat body transcriptome from the adult tsetse fly, Glossina morsitans morsitans

Tsetse flies (Diptera: Glossinidia) are vectors of pathogenic African trypanosomes. To develop a foundation for tsetse physiology, a normalized expressed sequence tag (EST) library was constructed from fat body tissue of immune-stimulated Glossina morsitans morsitans. Analysis of 20,257 high-quality ESTs yielded 6372 unique genes comprised of 3059 tentative consensus (TC) sequences and 3313 singletons (available at http://aksoylab.yale.edu). We analysed the putative fat body transcriptome based on homology to other gene products with known functions available in the public domain. In particular, we describe the immune-related products, reproductive function related yolk proteins and milk-gland protein, iron metabolism regulating ferritins and transferrin, and tsetse's major energy source proline biosynthesis. Expression analysis of the three yolk proteins indicates that all are detected in females, while only the yolk protein with similarity to lipases, is expressed in males. Milk gland protein, apparently important for larval nutrition, however, is primarily synthesized by accessory milk gland tissue.

Discovery of immune-related genes expressed in hemocytes of the tarantula spider Acanthoscurria gomesiana

The present study reports the identification of immune related transcripts from hemocytes of the spider Acanthoscurria gomesiana by high throughput sequencing of expressed sequence tags (ESTs). To generate ESTs from hemocytes, two cDNA libraries were prepared: one by directional cloning (primary) and the other by the normalization of the first (normalized). A total of 7584 clones were sequenced and the identical ESTs were clustered, resulting in 3723 assembled sequences (AS). At least 20% of these sequences are putative novel genes. The automatic functional annotation of AS based on Gene Ontology revealed several abundant transcripts related to the following functional classes: hemocyanin, lectin, and structural constituents of ribosome and cytoskeleton. From this annotation, 73 transcripts possibly involved in immune response were also identified, suggesting the existence of several molecular processes not previously described for spiders, such as: pathogen recognition, coagulation, complement activation, cell adhesion and intracellular signaling pathway for the activation of cellular defenses.

Insights into a dinoflagellate genome through expressed sequence tag analysis

Background

Dinoflagellates are important marine primary producers and grazers and cause toxic "red tides". These taxa are characterized by many unique features such as immense genomes, the absence of nucleosomes, and photosynthetic organelles (plastids) that have been gained and lost multiple times. We generated EST sequences from non-normalized and normalized cDNA libraries from a culture of the toxic species Alexandrium tamarense to elucidate dinoflagellate evolution. Previous analyses of these data have clarified plastid origin and here we study the gene content, annotate the ESTs, and analyze the genes that are putatively involved in DNA packaging.

Results

Approximately 20% of the 6,723 unique (11,171 total 3'-reads) ESTs data could be annotated using Blast searches against GenBank. Several putative dinoflagellate-specific mRNAs were identified, including one novel plastid protein. Dinoflagellate genes, similar to other eukaryotes, have a high GC-content that is reflected in the amino acid codon usage. Highly represented transcripts include histone-like (HLP) and luciferin binding proteins and several genes occur in families that encode nearly identical proteins. We also identified rare transcripts encoding a predicted protein highly similar to histone H2A.X. We speculate this histone may be retained for its role in DNA double-strand break repair.

Conclusion

This is the most extensive collection to date of ESTs from a toxic dinoflagellate. These data will be instrumental to future research to understand the unique and complex cell biology of these organisms and for potentially identifying the genes involved in toxin production.

1274 full-open reading frames of transcripts expressed in the developing mouse nervous system

As part of the trans-National Institutes of Health (NIH) Mouse Brain Molecular Anatomy Project (BMAP), and in close coordination with the NIH Mammalian Gene Collection Program (MGC), we initiated a large-scale project to clone, identify, and sequence the complete open reading frame (ORF) of transcripts expressed in the developing mouse nervous system. Here we report the analysis of the ORF sequence of 1274 cDNAs, obtained from 47 full-length-enriched cDNA libraries, constructed by using a novel approach, herein described. cDNA libraries were derived from size-fractionated cytoplasmic mRNA isolated from brain and eye tissues obtained at several embryonic stages and postnatal days. Altogether, including the full-ORF MGC sequences derived from these libraries by the MGC sequencing team, NIH_BMAP full-ORF sequences correspond to approximately 20% of all transcripts currently represented in mouse MGC. We show that NIH_BMAP clones comprise 68% of mouse MGC cDNAs > or =5 kb, and 54% of those > or =4 kb, as of March 15, 2004. Importantly, we identified transcripts, among the 1274 full-ORF sequences, that are exclusively or predominantly expressed in brain and eye tissues, many of which encode yet uncharacterized proteins.

Aedes aegypti genomics

The mosquito, Aedes aegypti, is the primary, worldwide arthropod vector for the yellow fever and dengue viruses. As it is also one of the most tractable mosquito species for laboratory studies, it has been and remains one of the most intensively studied arthropod species. This has resulted in the development of detailed genetic and physical maps for Ae. aegypti and considerable insight into its genome organization. The research community is well-advanced in developing important molecular tools that will facilitate a whole genome sequencing effort. This includes generation of BAC clone end sequences, physical mapping of selected BAC clones and generation of EST sequences. Whole genome sequence information for Ae. aegypti will provide important insight into mosquito chromosome evolution and allow for the identification of genes and gene function. These functions may be common to all mosquitoes or perhaps unique to individual species, possibly specific to host-seeking and blood-feeding behaviors, as well as the innate immune response to pathogens encountered during blood-feeding. This information will be invaluable to the global effort to develop novel strategies for preventing arthropod-borne disease transmission.