A. Sorana Morrissy

Next-generation tag sequencing for cancer gene expression profiling

We describe a new method, Tag-seq, which employs ultra high-throughput sequencing of 21 base pair cDNA tags for sensitive and cost-effective gene expression profiling. We compared Tag-seq data to LongSAGE data and observed improved representation of several classes of rare transcripts, including transcription factors, antisense transcripts, and intronic sequences, the latter possibly representing novel exons or genes. We observed increases in the diversity, abundance, and dynamic range of such rare transcripts and took advantage of the greater dynamic range of expression to identify, in cancers and normal libraries, altered expression ratios of alternative transcript isoforms. The strand-specific information of Tag-seq reads further allowed us to detect altered expression ratios of sense and antisense (S-AS) transcripts between cancer and normal libraries. S-AS transcripts were enriched in known cancer genes, while transcript isoforms were enriched in miRNA targeting sites. We found that transcript abundance had a stronger GC-bias in LongSAGE than Tag-seq, such that AT-rich tags were less abundant than GC-rich tags in LongSAGE. Tag-seq also performed better in gene discovery, identifying >98% of genes detected by LongSAGE and profiling a distinct subset of the transcriptome characterized by AT-rich genes, which was expressed at levels below those detectable by LongSAGE. Overall, Tag-seq is sensitive to rare transcripts, has less sequence composition bias relative to LongSAGE, and allows differential expression analysis for a greater range of transcripts, including transcripts encoding important regulatory molecules.

Alternative expression analysis by RNA sequencing

In alternative expression analysis by sequencing (ALEXA-seq), we developed a method to analyze massively parallel RNA sequence data to catalog transcripts and assess differential and alternative expression of known and predicted mRNA isoforms in cells and tissues. As proof of principle, we used the approach to compare fluorouracil-resistant and -nonresistant human colorectal cancer cell lines. We assessed the sensitivity and specificity of the approach by comparison to exon tiling and splicing microarrays and validated the results with reverse transcription–PCR, quantitative PCR and Sanger sequencing. We observed global disruption of splicing in fluorouracil-resistant cells characterized by expression of new mRNA isoforms resulting from exon skipping, alternative splice site usage and intron retention. Alternative expression annotation databases, source code, a data viewer and other resources to facilitate analysis are available at http://www.alexaplatform.org/alexa_seq/.

Cytogenetic Prognostication Within Medulloblastoma Subgroups

PURPOSE Medulloblastoma comprises four distinct molecular subgroups: WNT, SHH, Group 3, and Group 4. Current medulloblastoma protocols stratify patients based on clinical features: patient age, metastatic stage, extent of resection, and histologic variant. Stark prognostic and genetic differences among the four subgroups suggest that subgroup-specific molecular biomarkers could improve patient prognostication. PATIENTS AND METHODS Molecular biomarkers were identified from a discovery set of 673 medulloblastomas from 43 cities around the world. Combined risk stratification models were designed based on clinical and cytogenetic biomarkers identified by multivariable Cox proportional hazards analyses. Identified biomarkers were tested using fluorescent in situ hybridization (FISH) on a nonoverlapping medulloblastoma tissue microarray (n = 453), with subsequent validation of the risk stratification models. RESULTS Subgroup information improves the predictive accuracy of a multivariable survival model compared with clinical biomarkers alone. Most previously published cytogenetic biomarkers are only prognostic within a single medulloblastoma subgroup. Profiling six FISH biomarkers (GLI2, MYC, chromosome 11 [chr11], chr14, 17p, and 17q) on formalin-fixed paraffin-embedded tissues, we can reliably and reproducibly identify very low-risk and very high-risk patients within SHH, Group 3, and Group 4 medulloblastomas. CONCLUSION Combining subgroup and cytogenetic biomarkers with established clinical biomarkers substantially improves patient prognostication, even in the context of heterogeneous clinical therapies. The prognostic significance of most molecular biomarkers is restricted to a specific subgroup. We have identified a small panel of cytogenetic biomarkers that reliably identifies very high-risk and very low-risk groups of patients, making it an excellent tool for selecting patients for therapy intensification and therapy de-escalation in future clinical trials.

Intertumoral Heterogeneity within Medulloblastoma Subgroups

While molecular subgrouping has revolutionized medulloblastoma classification, the extent of heterogeneity within subgroups is unknown. Similarity network fusion (SNF) applied to genome-wide DNA methylation and gene expression data across 763 primary samples identifies very homogeneous clusters of patients, supporting the presence of medulloblastoma subtypes. After integration of somatic copy-number alterations, and clinical features specific to each cluster, we identify 12 different subtypes of medulloblastoma. Integrative analysis using SNF further delineates group 3 from group 4 medulloblastoma, which is not as readily apparent through analyses of individual data types. Two clear subtypes of infants with Sonic Hedgehog medulloblastoma with disparate outcomes and biology are identified. Medulloblastoma subtypes identified through integrative clustering have important implications for stratification of future clinical trials.

The whole-genome landscape of medulloblastoma subtypes

Current therapies for medulloblastoma, a highly malignant childhood brain tumour, impose debilitating effects on the developing child, and highlight the need for molecularly targeted treatments with reduced toxicity. Previous studies have been unable to identify the full spectrum of driver genes and molecular processes that operate in medulloblastoma subgroups. Here we analyse the somatic landscape across 491 sequenced medulloblastoma samples and the molecular heterogeneity among 1,256 epigenetically analysed cases, and identify subgroup-specific driver alterations that include previously undiscovered actionable targets. Driver mutations were confidently assigned to most patients belonging to Group 3 and Group 4 medulloblastoma subgroups, greatly enhancing previous knowledge. New molecular subtypes were differentially enriched for specific driver events, including hotspot in-frame insertions that target KBTBD4 and ‘enhancer hijacking’ events that activate PRDM6. Thus, the application of integrative genomics to an extensive cohort of clinical samples derived from a single childhood cancer entity revealed a series of cancer genes and biologically relevant subtype diversity that represent attractive therapeutic targets for the treatment of patients with medulloblastoma.

Extensive relationship between antisense transcription and alternative splicing in the human genome

To analyze the relationship between antisense transcription and alternative splicing, we developed a computational approach for the detection of antisense-correlated exon splicing events using Affymetrix exon array data. Our analysis of expression data from 176 lymphoblastoid cell lines revealed that the majority of expressed sense-antisense genes exhibited alternative splicing events that were correlated to the expression of the antisense gene. Most of these events occurred in areas of sense-antisense (SAS) gene overlap, which were significantly enriched in both exons and nucleosome occupancy levels relative to nonoverlapping regions of the same genes. Nucleosome occupancy was highly correlated with Pol II abundance across overlapping regions and with concomitant increases in local alternative exon usage. These results are consistent with an antisense transcription-mediated mechanism of splicing regulation in normal human cells. A comparison of the prevalence of antisense-correlated splicing events between individuals of Mormon versus African descent revealed population-specific events that may indicate the continued evolution of new SAS loci. Furthermore, the presence of antisense transcription was correlated to alternative splicing across multiple metazoan species, suggesting that it may be a conserved mechanism contributing to splicing regulation.

Hippo Signaling Influences HNF4A and FOXA2 Enhancer Switching during Hepatocyte Differentiation

Summary Cell fate acquisition is heavily influenced by direct interactions between master regulators and tissue-specific enhancers. However, it remains unclear how lineage-specifying transcription factors, which are often expressed in both progenitor and mature cell populations, influence cell differentiation. Using in vivo mouse liver development as a model, we identified thousands of enhancers that are bound by the master regulators HNF4A and FOXA2 in a differentiation-dependent manner, subject to chromatin remodeling, and associated with differentially expressed target genes. Enhancers exclusively occupied in the embryo were found to be responsive to developmentally regulated TEAD2 and coactivator YAP1. Our data suggest that Hippo signaling may affect hepatocyte differentiation by influencing HNF4A and FOXA2 interactions with temporal enhancers. In summary, transcription factor-enhancer interactions are not only tissue specific but also differentiation dependent, which is an important consideration for researchers studying cancer biology or mammalian development and/or using transformed cell lines.

Spatial heterogeneity in medulloblastoma

Spatial heterogeneity of transcriptional and genetic markers between physically isolated biopsies of a single tumor poses major barriers to the identification of biomarkers and the development of targeted therapies that will be effective against the entire tumor. We analyzed the spatial heterogeneity of multiregional biopsies from 35 patients, using a combination of transcriptomic and genomic profiles. Medulloblastomas (MBs), but not high-grade gliomas (HGGs), demonstrated spatially homogeneous transcriptomes, which allowed for accurate subgrouping of tumors from a single biopsy. Conversely, somatic mutations that affect genes suitable for targeted therapeutics demonstrated high levels of spatial heterogeneity in MB, malignant glioma, and renal cell carcinoma (RCC). Actionable targets found in a single MB biopsy were seldom clonal across the entire tumor, which brings the efficacy of monotherapies against a single target into question. Clinical trials of targeted therapies for MB should first ensure the spatially ubiquitous nature of the target mutation.

Tissue-specific alternative polyadenylation at the imprinted gene Mest regulates allelic usage at Copg2

The gene Mest (also known as Peg1) is regulated by genomic imprinting in the mouse and only the paternal allele is active for transcription. MEST is similarly imprinted in humans, where it is a candidate for the growth retardation Silver-Russell syndrome. The MEST protein belongs to an ancient family of hydrolases but its function is still unknown. It is highly conserved in vertebrates although imprinted expression is only observed in marsupials and eutherians, thus a recent evolutionary event. Here we describe the identification of new imprinted RNA products at the Mest locus, longer variants of the RNA, called MestXL, transcribed >10 kb into the downstream antisense gene Copg2. During development MestXL is produced exclusively in the developing central nervous system (CNS) by alternative polyadenylation. Copg2 is biallelically expressed in the embryo except in MestXL-expressing tissues, where we observed preferential expression from the maternal allele. To analyze the function of the MestXL transcripts in Copg2 regulation, we studied the effects of a targeted allele at Mest introducing a truncation in the mRNA. We show that both the formation of the MestXL isoforms and the allelic bias at Copg2 are lost in the CNS of mutants embryos. Our results propose a new mechanism to regulate allelic usage in the mammalian genome, via tissue-specific alternative polyadenylation and transcriptional interference in sense–antisense pairs at imprinted loci.

Sleeping Beauty Mouse Models Identify Candidate Genes Involved in Gliomagenesis

Genomic studies of human high-grade gliomas have discovered known and candidate tumor drivers. Studies in both cell culture and mouse models have complemented these approaches and have identified additional genes and processes important for gliomagenesis. Previously, we found that mobilization of Sleeping Beauty transposons in mice ubiquitously throughout the body from the Rosa26 locus led to gliomagenesis with low penetrance. Here we report the characterization of mice in which transposons are mobilized in the Glial Fibrillary Acidic Protein (GFAP) compartment. Glioma formation in these mice did not occur on an otherwise wild-type genetic background, but rare gliomas were observed when mobilization occurred in a p19Arf heterozygous background. Through cloning insertions from additional gliomas generated by transposon mobilization in the Rosa26 compartment, several candidate glioma genes were identified. Comparisons to genetic, epigenetic and mRNA expression data from human gliomas implicates several of these genes as tumor suppressor genes and oncogenes in human glioblastoma.