Divergent clonal selection dominates medulloblastoma at recurrence.

The development of targeted anti-cancer therapies through the study of cancer genomes is intended to increase survival rates and decrease treatment-related toxicity. We treated a transposon-driven, functional genomic mouse model of medulloblastoma with 'humanized' in vivo therapy (microneurosurgical tumour resection followed by multi-fractionated, image-guided radiotherapy). Genetic events in recurrent murine medulloblastoma exhibit a very poor overlap with those in matched murine diagnostic samples (<5%). Whole-genome sequencing of 33 pairs of human diagnostic and post-therapy medulloblastomas demonstrated substantial genetic divergence of the dominant clone after therapy (<12% diagnostic events were retained at recurrence). In both mice and humans, the dominant clone at recurrence arose through clonal selection of a pre-existing minor clone present at diagnosis. Targeted therapy is unlikely to be effective in the absence of the target, therefore our results offer a simple, proximal, and remediable explanation for the failure of prior clinical trials of targeted therapy.

Enhancer hijacking activates GFI1 family oncogenes in medulloblastoma.

Medulloblastoma is a highly malignant paediatric brain tumour currently treated with a combination of surgery, radiation and chemotherapy, posing a considerable burden of toxicity to the developing child. Genomics has illuminated the extensive intertumoral heterogeneity of medulloblastoma, identifying four distinct molecular subgroups. Group 3 and group 4 subgroup medulloblastomas account for most paediatric cases; yet, oncogenic drivers for these subtypes remain largely unidentified. Here we describe a series of prevalent, highly disparate genomic structural variants, restricted to groups 3 and 4, resulting in specific and mutually exclusive activation of the growth factor independent 1 family proto-oncogenes, GFI1 and GFI1B. Somatic structural variants juxtapose GFI1 or GFI1B coding sequences proximal to active enhancer elements, including super-enhancers, instigating oncogenic activity. Our results, supported by evidence from mouse models, identify GFI1 and GFI1B as prominent medulloblastoma oncogenes and implicate 'enhancer hijacking' as an efficient mechanism driving oncogene activation in a childhood cancer.

Correspondence regarding "Effect of active smoking on the human bronchial epithelium transcriptome".

BACKGROUND: In the work of Chari et al. entitled "Effect of active smoking on the human bronchial epithelium transcriptome" the authors use SAGE to identify candidate gene expression changes in bronchial brushings from never, former, and current smokers. These gene expression changes are categorized into those that are reversible or irreversible upon smoking cessation. A subset of these identified genes is validated on an independent cohort using RT-PCR. The authors conclude that their results support the notion of gene expression changes in the lungs of smokers which persist even after an individual has quit. RESULTS: This correspondence raises questions about the validity of the approach used by the authors to analyze their data. The majority of the reported results suffer deficiencies due to the methods used. The most fundamental of these are explained in detail: biases introduced during data processing, lack of correction for multiple testing, and an incorrect use of clustering for gene discovery. A randomly generated "null" dataset is used to show the consequences of these shortcomings. CONCLUSION: Most of Chari et al.'s findings are consistent with what would be expected by chance alone. Although there is clear evidence of reversible changes in gene expression, the majority of those identified appear to be false positives. However, contrary to the authors' claims, no irreversible changes were identified. There is a broad consensus that genetic change due to smoking persists once an individual has quit smoking; unfortunately, this study lacks sufficient scientific rigour to support or refute this hypothesis or identify any specific candidate genes. The pitfalls of large-scale analysis, as exemplified here, may not be unique to Chari et al.

Whole transcriptome analysis reveals differential gene expression profile reflecting macrophage polarization in response to influenza A H5N1 virus infection.

BACKGROUND: Avian influenza A H5N1 virus can cause lethal disease in humans. The virus can trigger severe pneumonia and lead to acute respiratory distress syndrome. Data from clinical, in vitro and in vivo suggest that virus-induced cytokine dysregulation could be a contributory factor to the pathogenesis of human H5N1 disease. However, the precise mechanism of H5N1 infection eliciting the unique host response are still not well understood. METHODS: To obtain a better understanding of the molecular events at the earliest time points, we used RNA-Seq to quantify and compare the host mRNA and miRNA transcriptomes induced by the highly pathogenic influenza A H5N1 (A/Vietnam/3212/04) or low virulent H1N1 (A/Hong Kong/54/98) viruses in human monocyte-derived macrophages at 1-, 3-, and 6-h post infection. RESULTS: Our data reveals that two macrophage populations corresponding to M1 (classically activated) and M2 (alternatively activated) macrophage subtypes respond distinctly to H5N1 virus infection when compared to H1N1 virus or mock infection, a distinction that could not be made from previous microarray studies. When this confounding variable is considered in our statistical model, a clear set of dysregulated genes and pathways emerges specifically in H5N1 virus-infected macrophages at 6-h post infection, whilst was not found with H1N1 virus infection. Furthermore, altered expression of genes in these pathways, which have been previously implicated in viral host response, occurs specifically in the M1 subtype. We observe a significant up-regulation of genes in the RIG-I-like receptor signaling pathway. In particular, interferons, and interferon-stimulated genes are broadly affected. The negative regulators of interferon signaling, the suppressors of cytokine signaling, SOCS-1 and SOCS-3, were found to be markedly up-regulated in the initial round of H5N1 virus replication. Elevated levels of these suppressors could lead to the eventual suppression of cellular antiviral genes, contributing to pathophysiology of H5N1 virus infection. CONCLUSIONS: Our study provides important mechanistic insights into the understanding of H5N1 viral pathogenesis and the multi-faceted host immune responses. The dysregulated genes could be potential candidates as therapeutic targets for treating H5N1 disease.

Statistical analysis and significance testing of serial analysis of gene expression data using a Poisson mixture model.

BACKGROUND: Serial analysis of gene expression (SAGE) is used to obtain quantitative snapshots of the transcriptome. These profiles are count-based and are assumed to follow a Binomial or Poisson distribution. However, tag counts observed across multiple libraries (for example, one or more groups of biological replicates) have additional variance that cannot be accommodated by this assumption alone. Several models have been proposed to account for this effect, all of which utilize a continuous prior distribution to explain the excess variance. Here, a Poisson mixture model, which assumes excess variability arises from sampling a mixture of distinct components, is proposed and the merits of this model are discussed and evaluated. RESULTS: The goodness of fit of the Poisson mixture model on 15 sets of biological SAGE replicates is compared to the previously proposed hierarchical gamma-Poisson (negative binomial) model, and a substantial improvement is seen. In further support of the mixture model, there is observed: 1) an increase in the number of mixture components needed to fit the expression of tags representing more than one transcript; and 2) a tendency for components to cluster libraries into the same groups. A confidence score is presented that can identify tags that are differentially expressed between groups of SAGE libraries. Several examples where this test outperforms those previously proposed are highlighted. CONCLUSION: The Poisson mixture model performs well as a) a method to represent SAGE data from biological replicates, and b) a basis to assign significance when testing for differential expression between multiple groups of replicates. Code for the R statistical software package is included to assist investigators in applying this model to their own data.

A SAGE approach to discovery of genes involved in autophagic cell death.

Programmed cell death (PCD), important in normal animal physiology and disease, can be divided into at least two morphological subtypes, including type I, or apoptosis, and type II, or autophagic cell death. While many molecules involved in apoptosis have been discovered and studied intensively during the past decade, autophagic cell death is not well characterized molecularly. Here we report the first comprehensive identification of molecules associated with autophagic cell death during normal metazoan development in vivo. During Drosophila metamorphosis, the larval salivary glands undergo autophagic cell death regulated by a hormonally induced transcriptional cascade. To identify and analyze the genes expressed, we examined wild-type patterns of gene expression in three predeath stages of Drosophila salivary glands using serial analysis of gene expression (SAGE) [7]. 1244 transcripts, including genes involved in autophagy, defense response, cytoskeleton remodeling, noncaspase proteolysis, and apoptosis, were expressed differentially prior to salivary gland death. Mutant expression analysis indicated that several of these genes were regulated by E93, a gene required for salivary gland cell death. Our analyses strongly support both the emerging notion that there is overlap with respect to the molecules involved in autophagic cell death and apoptosis, and that there are important differences.

An efficient strategy for large-scale high-throughput transposon-mediated sequencing of cDNA clones.

We describe an efficient high-throughput method for accurate DNA sequencing of entire cDNA clones. Developed as part of our involvement in the Mammalian Gene Collection full-length cDNA sequencing initiative, the method has been used and refined in our laboratory since September 2000. Amenable to large scale projects, we have used the method to generate >7 Mb of accurate sequence from 3695 candidate full-length cDNAs. Sequencing is accomplished through the insertion of Mu transposon into cDNAs, followed by sequencing reactions primed with Mu-specific sequencing primers. Transposon insertion reactions are not performed with individual cDNAs but rather on pools of up to 96 clones. This pooling strategy reduces the number of transposon insertion sequencing libraries that would otherwise be required, reducing the costs and enhancing the efficiency of the transposon library construction procedure. Sequences generated using transposon-specific sequencing primers are assembled to yield the full-length cDNA sequence, with sequence editing and other sequence finishing activities performed as required to resolve sequence ambiguities. Although analysis of the many thousands (22 785) of sequenced Mu transposon insertion events revealed a weak sequence preference for Mu insertion, we observed insertion of the Mu transposon into 1015 of the possible 1024 5mer candidate insertion sites.

A knowledge discovery object model API for Java.

BACKGROUND: Biological data resources have become heterogeneous and derive from multiple sources. This introduces challenges in the management and utilization of this data in software development. Although efforts are underway to create a standard format for the transmission and storage of biological data, this objective has yet to be fully realized. RESULTS: This work describes an application programming interface (API) that provides a framework for developing an effective biological knowledge ontology for Java-based software projects. The API provides a robust framework for the data acquisition and management needs of an ontology implementation. In addition, the API contains classes to assist in creating GUIs to represent this data visually. CONCLUSIONS: The Knowledge Discovery Object Model (KDOM) API is particularly useful for medium to large applications, or for a number of smaller software projects with common characteristics or objectives. KDOM can be coupled effectively with other biologically relevant APIs and classes. Source code, libraries, documentation and examples are available at http://www.bcgsc.ca/bioinfo/software.

DiscoverySpace: an interactive data analysis application.

DiscoverySpace is a graphical application for bioinformatics data analysis. Users can seamlessly traverse references between biological databases and draw together annotations in an intuitive tabular interface. Datasets can be compared using a suite of novel tools to aid in the identification of significant patterns. DiscoverySpace is of broad utility and its particular strength is in the analysis of serial analysis of gene expression (SAGE) data. The application is freely available online.

Large-scale production of SAGE libraries from microdissected tissues, flow-sorted cells, and cell lines.

We describe the details of a serial analysis of gene expression (SAGE) library construction and analysis platform that has enabled the generation of >298 high-quality SAGE libraries and >30 million SAGE tags primarily from sub-microgram amounts of total RNA purified from samples acquired by microdissection. Several RNA isolation methods were used to handle the diversity of samples processed, and various measures were applied to minimize ditag PCR carryover contamination. Modifications in the SAGE protocol resulted in improved cloning and DNA sequencing efficiencies. Bioinformatic measures to automatically assess DNA sequencing results were implemented to analyze the integrity of ditag structure, linker or cross-species ditag contamination, and yield of high-quality tags per sequence read. Our analysis of singleton tag errors resulted in a method for correcting such errors to statistically determine tag accuracy. From the libraries generated, we produced an essentially complete mapping of reliable 21-base-pair tags to the mouse reference genome sequence for a meta-library of approximately 5 million tags. Our analyses led us to reject the commonly held notion that duplicate ditags are artifacts. Rather than the usual practice of discarding such tags, we conclude that they should be retained to avoid introducing bias into the results and thereby maintain the quantitative nature of the data, which is a major theoretical advantage of SAGE as a tool for global transcriptional profiling.

Iron-regulated transcription and capsule formation in the fungal pathogen Cryptococcus neoformans.

Cryptococcus neoformans is the leading cause of fungal meningitis in humans. Production of a polysaccharide capsule is a key virulence property for the fungus and capsule synthesis is regulated by iron levels. Given that iron acquisition is an important aspect of virulence for many pathogens, we employed serial analysis of gene expression (SAGE) to examine the transcriptome under iron-limiting and iron-replete conditions. Initially, we demonstrated by SAGE and Northern analysis that iron limitation results in an elevated transcript level for the CAP60 gene that is required for capsule production. We also identified genes encoding putative components for iron transport and homeostasis, including the FTR1 (iron permease) gene, with higher transcript levels in the low-iron condition. An FTR1 disruption mutant grows more slowly than wild-type cells in low-iron medium, and shows delayed growth and altered capsule regulation in iron-replete medium. Iron deprivation also resulted in elevated SAGE tags for putative extracellular mannoproteins and the GPI8 gene encoding a glycosylphosphatidylinositol (GPI) transamidase. The GPI8 gene appears to be essential while disruption of the CIG1 gene encoding a mannoprotein resulted in impaired growth in low-iron medium and altered capsule response to the iron-replete condition. Additionally, we found that iron-replete conditions led to elevated transcripts for genes for iron storage, nitrogen metabolism, glycolysis, mitochondrial function, lipid metabolism and calmodulin-calcineurin signalling. Overall, these studies provide the first view of the C. neoformans transcriptional response to different iron levels.

Effect of TERT and ATM on gene expression profiles in human fibroblasts.

Telomeres protect chromosomes from degradation, end-to-end fusion, and illegitimate recombination. Loss of telomeres may lead to cell death or senescence or may cause genomic instability, leading to tumor formation. Expression of human telomerase reverse transcriptase (TERT) in human fibroblast cells elongates their telomeres and extends their lifespan. Ataxia telangiectasia mutated (ATM) deficiency in A-T human fibroblasts results in accelerated telomere shortening, abnormal cell-cycle response to DNA damage, and early senescence. Gene expression profiling was performed by serial analysis of gene expression (SAGE) on BJ normal human skin fibroblasts, A-T cells, and BJ and A-T cells transduced with TERT cDNA and expressing telomerase activity. In the four SAGE libraries, 36,921 unique SAGE tags were detected. Pairwise comparisons between the libraries showed differential expression levels of 1%-8% of the tags. Transcripts affected by both TERT and ATM were identified according to expression patterns, making them good candidates for further studies of pathways affected by both TERT and ATM. These include MT2A, P4HB, LGALS1, CFL1, LDHA, S100A10, EIF3S8, RANBP9, and SEC63. These genes are involved in apoptosis or processes related to cell growth, and most have been found to be deregulated in cancer. Our results have provided further insight into the roles of TERT and ATM by identifying genes likely to be involved in their function. Supplementary material for this article can be found on the Genes, Chromosomes and Cancer website at http://www.interscience.wiley.com/jpages/1045-2257/suppmat/index.html.

Temperature-regulated transcription in the pathogenic fungus Cryptococcus neoformans.

The basidiomycete fungus Cryptococcus neoformans is an opportunistic pathogen of worldwide importance that causes meningitis, leading to death in immunocompromised individuals. Unlike many basidiomycete fungi, C. neoformans is thermotolerant, and its ability to grow at 37 degrees C is considered to be a virulence factor. We used serial analysis of gene expression (SAGE) to characterize the transcriptomes of C. neoformans strains that represent two varieties with different polysaccharide capsule serotypes. These include a serotype D strain of the C. neoformans variety neoformans and a serotype A strain of variety grubii. In this report, we describe the construction and characterization of SAGE libraries from each strain grown at 25 degrees C and 37 degrees C. The SAGE data reveal transcriptome differences between the two strains, even at this early stage of analysis, and identify sets of genes with higher transcript levels at 25 degrees C or 37 degrees C. Notably, growth at the lower temperature increased transcript levels for histone genes, indicating a general influence of temperature on chromatin structure. At 37 degrees C, we noted elevated transcript levels for several genes encoding heat shock proteins and translation machinery. Some of these genes may play a role in temperature-regulated phenotypes in C. neoformans, such as the adaptation of the fungus to growth in the host and the dimorphic transition between budding and filamentous growth. Overall, this work provides the most comprehensive gene expression data available for C. neoformans; this information will be a critical resource both for gene discovery and genome annotation in this pathogen.