Mitogen-induced defective mitosis transforms neural progenitor cells.

BACKGROUND: Chromosome instability (CIN) with recurrent copy number alterations is a feature of many solid tumors, including glioblastoma (GBM), yet the genes that regulate cell division are rarely mutated in cancers. Here, we show that the brain-abundant mitogen, platelet-derived growth factor-A (PDGFA) fails to induce the expression of kinetochore and spindle assembly checkpoint genes leading to defective mitosis in neural progenitor cells (NPCs). METHODS: Using a recently reported in vitro model of the initiation of high-grade gliomas from murine NPCs, we investigated the immediate effects of PDGFA exposure on the nuclear and mitotic phenotypes and patterns of gene and protein expression in NPCs, a putative GBM cell of origin. RESULTS: NPCs divided abnormally in defined media containing PDGFA with P53-dependent effects. In wild-type cells, defective mitosis was associated with P53 activation and cell death, but in some null cells, defective mitosis was tolerated. Surviving cells had unstable genomes and proliferated in the presence of PDGFA accumulating random and clonal chromosomal rearrangements. The outcome of this process was a population of tumorigenic NPCs with recurrent gains and losses of chromosomal regions that were syntenic to those recurrently gained and lost in human GBM. By stimulating proliferation without setting the stage for successful mitosis, PDGFA-transformed NPCs lacking P53 function. CONCLUSIONS: Our work describes a mechanism of transformation of NPCs by a brain-associated mitogen, raising the possibility that the unique genomic architecture of GBM is an adaptation to defective mitosis that ensures the survival of affected cells.

The evolution of pre-mRNA splicing and its machinery revealed by reduced extremophilic red algae.

The Cyanidiales are a group of mostly thermophilic and acidophilic red algae that thrive near volcanic vents. Despite their phylogenetic relationship, the reduced genomes of Cyanidioschyzon merolae and Galdieria sulphuraria are strikingly different with respect to pre-mRNA splicing, a ubiquitous eukaryotic feature. Introns are rare and spliceosomal machinery is extremely reduced in C. merolae, in contrast to G. sulphuraria. Previous studies also revealed divergent spliceosomes in the mesophilic red alga Porphyridium purpureum and the red algal derived plastid of Guillardia theta (Cryptophyta), along with unusually high levels of unspliced transcripts. To further examine the evolution of splicing in red algae, we compared C. merolae and G. sulphuraria, investigating splicing levels, intron position, intron sequence features, and the composition of the spliceosome. In addition to identifying 11 additional introns in C. merolae, our transcriptomic analysis also revealed typical eukaryotic splicing in G. sulphuraria, whereas most transcripts in C. merolae remain unspliced. The distribution of intron positions within their host genes was examined to provide insight into patterns of intron loss in red algae. We observed increasing variability of 5' splice sites and branch donor regions with increasing intron richness. We also found these relationships to be connected to reductions in and losses of corresponding parts of the spliceosome. Our findings highlight patterns of intron and spliceosome evolution in related red algae under the pressures of genome reduction.

CIViCdb 2022: evolution of an open-access cancer variant interpretation knowledgebase.

CIViC (Clinical Interpretation of Variants in Cancer; civicdb.org) is a crowd-sourced, public domain knowledgebase composed of literature-derived evidence characterizing the clinical utility of cancer variants. As clinical sequencing becomes more prevalent in cancer management, the need for cancer variant interpretation has grown beyond the capability of any single institution. CIViC contains peer-reviewed, published literature curated and expertly-moderated into structured data units (Evidence Items) that can be accessed globally and in real time, reducing barriers to clinical variant knowledge sharing. We have extended CIViC's functionality to support emergent variant interpretation guidelines, increase interoperability with other variant resources, and promote widespread dissemination of structured curated data. To support the full breadth of variant interpretation from basic to translational, including integration of somatic and germline variant knowledge and inference of drug response, we have enabled curation of three new Evidence Types (Predisposing, Oncogenic and Functional). The growing CIViC knowledgebase has over 300 contributors and distributes clinically-relevant cancer variant data currently representing >3200 variants in >470 genes from >3100 publications.

A platform for oncogenomic reporting and interpretation.

Manual interpretation of variants remains rate limiting in precision oncology. The increasing scale and complexity of molecular data generated from comprehensive sequencing of cancer samples requires advanced interpretative platforms as precision oncology expands beyond individual patients to entire populations. To address this unmet need, we introduce a Platform for Oncogenomic Reporting and Interpretation (PORI), comprising an analytic framework that facilitates the interpretation and reporting of somatic variants in cancer. PORI integrates reporting and graph knowledge base tools combined with support for manual curation at the reporting stage. PORI represents an open-source platform alternative to commercial reporting solutions suitable for comprehensive genomic data sets in precision oncology. We demonstrate the utility of PORI by matching 9,961 pan-cancer genome atlas tumours to the graph knowledge base, calculating therapeutically informative alterations, and making available reports describing select individual samples.

Rearrangement-mediated cis-regulatory alterations in advanced patient tumors reveal interactions with therapy.

The global impact of somatic structural variants (SVs) on gene regulation in advanced tumors with complex treatment histories has been mostly uncharacterized. Here, using whole-genome and RNA sequencing from 570 recurrent or metastatic tumors, we report the altered expression of hundreds of genes in association with nearby SV breakpoints, including oncogenes and G-protein-coupled receptor-related genes such as PLEKHG2. A significant fraction of genes with SV-expression associations correlate with worse patient survival in primary and advanced cancers, including SRD5A1. In many instances, SV-expression associations involve retrotransposons being translocated near genes. High overall SV burden is associated with treatment with DNA alkylating agents or taxanes and altered expression of metabolism-associated genes. SV-expression associations within tumors from topoisomerase I inhibitor-treated patients include chromatin-related genes. Within anthracycline-treated tumors, SV breakpoints near chromosome 1p genes include PDE4B. Patient treatment and history can help understand the widespread SV-mediated cis-regulatory alterations found in cancer.

Uncovering Clinically Relevant Gene Fusions with Integrated Genomic and Transcriptomic Profiling of Metastatic Cancers.

PURPOSE: Gene fusions are important oncogenic drivers and many are actionable. Whole-genome and transcriptome (WGS and RNA-seq, respectively) sequencing can discover novel clinically relevant fusions. EXPERIMENTAL DESIGN: Using WGS and RNA-seq, we reviewed the prevalence of fusions in a cohort of 570 patients with cancer, and compared prevalence to that predicted with commercially available panels. Fusions were annotated using a consensus variant calling pipeline (MAVIS) and required that a contig of the breakpoint could be constructed and supported from ≥2 structural variant detection approaches. RESULTS: In 570 patients with advanced cancer, MAVIS identified 81 recurrent fusions by WGS and 111 by RNA-seq, of which 18 fusions by WGS and 19 by RNA-seq were noted in at least 3 separate patients. The most common fusions were EML4-ALK in thoracic malignancies (9/69, 13%), and CMTM8-CMTM7 in colorectal cancer (4/73, 5.5%). Combined genomic and transcriptomic analysis identified novel fusion partners for clinically relevant genes, such as NTRK2 (novel partners: SHC3, DAPK1), and NTRK3 (novel partners: POLG, PIBF1). CONCLUSIONS: Utilizing WGS/RNA-seq facilitates identification of novel fusions in clinically relevant genes, and detected a greater proportion than commercially available panels are expected to find. A significant benefit of WGS and RNA-seq is the innate ability to retrospectively identify variants that becomes clinically relevant over time, without the need for additional testing, which is not possible with panel-based approaches.

Genome and Transcriptome Biomarkers of Response to Immune Checkpoint Inhibitors in Advanced Solid Tumors.

PURPOSE: Immune checkpoint inhibitors (ICI) have revolutionized the treatment of solid tumors with dramatic and durable responses seen across multiple tumor types. However, identifying patients who will respond to these drugs remains challenging, particularly in the context of advanced and previously treated cancers. EXPERIMENTAL DESIGN: We characterized fresh tumor biopsies from a heterogeneous pan-cancer cohort of 98 patients with metastatic predominantly pretreated disease through the Personalized OncoGenomics program at BC Cancer (Vancouver, Canada) using whole genome and transcriptome analysis (WGTA). Baseline characteristics and follow-up data were collected retrospectively. RESULTS: We found that tumor mutation burden, independent of mismatch repair status, was the most predictive marker of time to progression (P = 0.007), but immune-related CD8+ T-cell and M1-M2 macrophage ratio scores were more predictive for overall survival (OS; P = 0.0014 and 0.0012, respectively). While CD274 [programmed death-ligand 1 (PD-L1)] gene expression is comparable with protein levels detected by IHC, we did not observe a clinical benefit for patients with this marker. We demonstrate that a combination of markers based on WGTA provides the best stratification of patients (P = 0.00071, OS), and also present a case study of possible acquired resistance to pembrolizumab in a patient with non-small cell lung cancer. CONCLUSIONS: Interpreting the tumor-immune interface to predict ICI efficacy remains challenging. WGTA allows for identification of multiple biomarkers simultaneously that in combination may help to identify responders, particularly in the context of a heterogeneous population of advanced and previously treated cancers, thus precluding tumor type-specific testing.

Glioma-derived IL-33 orchestrates an inflammatory brain tumor microenvironment that accelerates glioma progression.

Despite a deeper molecular understanding, human glioblastoma remains one of the most treatment refractory and fatal cancers. It is known that the presence of macrophages and microglia impact glioblastoma tumorigenesis and prevent durable response. Herein we identify the dual function cytokine IL-33 as an orchestrator of the glioblastoma microenvironment that contributes to tumorigenesis. We find that IL-33 expression in a large subset of human glioma specimens and murine models correlates with increased tumor-associated macrophages/monocytes/microglia. In addition, nuclear and secreted functions of IL-33 regulate chemokines that collectively recruit and activate circulating and resident innate immune cells creating a pro-tumorigenic environment. Conversely, loss of nuclear IL-33 cripples recruitment, dramatically suppresses glioma growth, and increases survival. Our data supports the paradigm that recruitment and activation of immune cells, when instructed appropriately, offer a therapeutic strategy that switches the focus from the cancer cell alone to one that includes the normal host environment.

Pan-cancer analysis of advanced patient tumors reveals interactions between therapy and genomic landscapes.

Advanced and metastatic tumors with complex treatment histories drive cancer mortality. Here we describe the POG570 cohort, a comprehensive whole-genome, transcriptome and clinical dataset, amenable for exploration of the impacts of therapies on genomic landscapes. Previous exposure to DNA-damaging chemotherapies and mutations affecting DNA repair genes, including POLQ and genes encoding Polζ, were associated with genome-wide, therapy-induced mutagenesis. Exposure to platinum therapies coincided with signatures SBS31 and DSB5 and, when combined with DNA synthesis inhibitors, signature SBS17b. Alterations in ESR1, EGFR, CTNNB1, FGFR1, VEGFA and DPYD were consistent with drug resistance and sensitivity. Recurrent noncoding events were found in regulatory region hotspots of genes including TERT, PLEKHS1, AP2A1 and ADGRG6. Mutation burden and immune signatures corresponded with overall survival and response to immunotherapy. Our data offer a rich resource for investigation of advanced cancers and interpretation of whole-genome and transcriptome sequencing in the context of a cancer clinic.

Standard operating procedure for curation and clinical interpretation of variants in cancer.

Manually curated variant knowledgebases and their associated knowledge models are serving an increasingly important role in distributing and interpreting variants in cancer. These knowledgebases vary in their level of public accessibility, and the complexity of the models used to capture clinical knowledge. CIViC (Clinical Interpretation of Variants in Cancer - www.civicdb.org) is a fully open, free-to-use cancer variant interpretation knowledgebase that incorporates highly detailed curation of evidence obtained from peer-reviewed publications and meeting abstracts, and currently holds over 6300 Evidence Items for over 2300 variants derived from over 400 genes. CIViC has seen increased adoption by, and also undertaken collaboration with, a wide range of users and organizations involved in research. To enhance CIViC's clinical value, regular submission to the ClinVar database and pursuit of other regulatory approvals is necessary. For this reason, a formal peer reviewed curation guideline and discussion of the underlying principles of curation is needed. We present here the CIViC knowledge model, standard operating procedures (SOP) for variant curation, and detailed examples to support community-driven curation of cancer variants.

Comprehensive genomic profiling of glioblastoma tumors, BTICs, and xenografts reveals stability and adaptation to growth environments.

Glioblastoma multiforme (GBM) is the most deadly brain tumor, and currently lacks effective treatment options. Brain tumor-initiating cells (BTICs) and orthotopic xenografts are widely used in investigating GBM biology and new therapies for this aggressive disease. However, the genomic characteristics and molecular resemblance of these models to GBM tumors remain undetermined. We used massively parallel sequencing technology to decode the genomes and transcriptomes of BTICs and xenografts and their matched tumors in order to delineate the potential impacts of the distinct growth environments. Using data generated from whole-genome sequencing of 201 samples and RNA sequencing of 118 samples, we show that BTICs and xenografts resemble their parental tumor at the genomic level but differ at the mRNA expression and epigenomic levels, likely due to the different growth environment for each sample type. These findings suggest that a comprehensive genomic understanding of in vitro and in vivo GBM model systems is crucial for interpreting data from drug screens, and can help control for biases introduced by cell-culture conditions and the microenvironment in mouse models. We also found that lack of MGMT expression in pretreated GBM is linked to hypermutation, which in turn contributes to increased genomic heterogeneity and requires new strategies for GBM treatment.

Relative Mutation Rates in Nucleomorph-Bearing Algae.

Chlorarachniophyte and cryptophyte algae are unique among plastid-containing species in that they have a nucleomorph genome: a compact, highly reduced nuclear genome from a photosynthetic eukaryotic endosymbiont. Despite their independent origins, the nucleomorph genomes of these two lineages have similar genomic architectures, but little is known about the evolutionary pressures impacting nucleomorph DNA, particularly how their rates of evolution compare to those of the neighboring genetic compartments (the mitochondrion, plastid, and nucleus). Here, we use synonymous substitution rates to estimate relative mutation rates in the four genomes of nucleomorph-bearing algae. We show that the relative mutation rates of the host versus endosymbiont nuclear genomes are similar in both chlorarachniophytes and cryptophytes, despite the fact that nucleomorph gene sequences are notoriously highly divergent. There is some evidence, however, for slightly elevated mutation rates in the nucleomorph DNA of chlorarachniophytes-a feature not observed in that of cryptophytes. For both lineages, relative mutation rates in the plastid appear to be lower than those in the nucleus and nucleomorph (and, in one case, the mitochondrion), which is consistent with studies of other plastid-bearing protists. Given the divergent nature of nucleomorph genes, our finding of relatively low evolutionary rates in these genomes suggests that for both lineages a burst of evolutionary change and/or decreased selection pressures likely occurred early in the integration of the secondary endosymbiont.

Evolution and Diversity of Pre-mRNA Splicing in Highly Reduced Nucleomorph Genomes.

Eukaryotic genes are interrupted by introns that are removed in a conserved process known as pre-mRNA splicing. Though well-studied in select model organisms, we are only beginning to understand the variation and diversity of this process across the tree of eukaryotes. We explored pre-mRNA splicing and other features of transcription in nucleomorphs, the highly reduced remnant nuclei of secondary endosymbionts. Strand-specific transcriptomes were sequenced from the cryptophyte Guillardia theta and the chlorarachniophyte Bigelowiella natans, whose plastids are derived from red and green algae, respectively. Both organisms exhibited elevated nucleomorph antisense transcription and gene expression relative to their respective nuclei, suggesting unique properties of gene regulation and transcriptional control in nucleomorphs. Marked differences in splicing were observed between the two nucleomorphs: the few introns of the G. theta nucleomorph were largely retained in mature transcripts, whereas the many short introns of the B. natans nucleomorph are spliced at typical eukaryotic levels (>90%). These differences in splicing levels could be reflecting the ancestries of the respective plastids, the different intron densities due to independent genome reduction events, or a combination of both. In addition to extending our understanding of the diversity of pre-mRNA splicing across eukaryotes, our study also indicates potential links between splicing, antisense transcription, and gene regulation in reduced genomes.

The New Red Algal Subphylum Proteorhodophytina Comprises the Largest and Most Divergent Plastid Genomes Known.

Red algal plastid genomes are often considered ancestral and evolutionarily stable, and thus more closely resembling the last common ancestral plastid genome of all photosynthetic eukaryotes [1, 2]. However, sampling of red algal diversity is still quite limited (e.g., [2-5]). We aimed to remedy this problem. To this end, we sequenced six new plastid genomes from four undersampled and phylogenetically disparate red algal classes (Porphyridiophyceae, Stylonematophyceae, Compsopogonophyceae, and Rhodellophyceae) and discovered an unprecedented degree of genomic diversity among them. These genomes are rich in introns, enlarged intergenic regions, and transposable elements (in the rhodellophycean Bulboplastis apyrenoidosa), and include the largest and most intron-rich plastid genomes ever sequenced (that of the rhodellophycean Corynoplastis japonica; 1.13 Mbp). Sophisticated phylogenetic analyses accounting for compositional heterogeneity show that these four "basal" red algal classes form a larger monophyletic group, Proteorhodophytina subphylum nov., and confidently resolve the large-scale relationships in the Rhodophyta. Our analyses also suggest that secondary red plastids originated before the diversification of all mesophilic red algae. Our genomic survey has challenged the current paradigmatic view of red algal plastid genomes as "living fossils" [1, 2, 6] by revealing an astonishing degree of divergence in size, organization, and non-coding DNA content. A closer look at red algae shows that they comprise the most ancestral (e.g., [2, 7, 8]) as well as some of the most divergent plastid genomes known.

Dramatically reduced spliceosome in Cyanidioschyzon merolae.

The human spliceosome is a large ribonucleoprotein complex that catalyzes pre-mRNA splicing. It consists of five snRNAs and more than 200 proteins. Because of this complexity, much work has focused on the Saccharomyces cerevisiae spliceosome, viewed as a highly simplified system with fewer than half as many splicing factors as humans. Nevertheless, it has been difficult to ascribe a mechanistic function to individual splicing factors or even to discern which are critical for catalyzing the splicing reaction. We have identified and characterized the splicing machinery from the red alga Cyanidioschyzon merolae, which has been reported to harbor only 26 intron-containing genes. The U2, U4, U5, and U6 snRNAs contain expected conserved sequences and have the ability to adopt secondary structures and form intermolecular base-pairing interactions, as in other organisms. C. merolae has a highly reduced set of 43 identifiable core splicing proteins, compared with ∼90 in budding yeast and ∼140 in humans. Strikingly, we have been unable to find a U1 snRNA candidate or any predicted U1-associated proteins, suggesting that splicing in C. merolae may occur without the U1 small nuclear ribonucleoprotein particle. In addition, based on mapping the identified proteins onto the known splicing cycle, we propose that there is far less compositional variability during splicing in C. merolae than in other organisms. The observed reduction in splicing factors is consistent with the elimination of spliceosomal components that play a peripheral or modulatory role in splicing, presumably retaining those with a more central role in organization and catalysis.

Transcriptome analysis of the parasite Encephalitozoon cuniculi: an in-depth examination of pre-mRNA splicing in a reduced eukaryote.

BACKGROUND: The microsporidian Encephalitozoon cuniculi possesses one of the most reduced and compacted eukaryotic genomes. Reduction in this intracellular parasite has affected major cellular machinery, including the loss of over fifty core spliceosomal components compared to S. cerevisiae. To identify expression changes throughout the parasite's life cycle and also to assess splicing in the context of this reduced system, we examined the transcriptome of E. cuniculi using Illumina RNA-seq. RESULTS: We observed that nearly all genes are expressed at three post-infection time-points examined. A large fraction of genes are differentially expressed between the first and second (37.7%) and first and third (43.8%) time-points, while only four genes are differentially expressed between the latter two. Levels of intron splicing are very low, with 81% of junctions spliced at levels below 50%. This is dramatically lower than splicing levels found in two other fungal species examined. We also describe the first case of alternative splicing in a microsporidian, an unexpected complexity given the reduction in spliceosomal components. CONCLUSIONS: Low levels of splicing observed are likely the result of an inefficient spliceosome; however, at least in one case, splicing appears to be playing a functional role. Although several RNA decay genes are encoded in E. cuniculi, the lack of a few key players could be reducing decay levels and therefore increasing the proportion of unspliced transcripts. Significant proportions of genes are differentially expressed in the first forty-eight hours but not after, indicative of genetic changes that precede the intracellular to infective stage transition.

Algal genomes reveal evolutionary mosaicism and the fate of nucleomorphs.

Cryptophyte and chlorarachniophyte algae are transitional forms in the widespread secondary endosymbiotic acquisition of photosynthesis by engulfment of eukaryotic algae. Unlike most secondary plastid-bearing algae, miniaturized versions of the endosymbiont nuclei (nucleomorphs) persist in cryptophytes and chlorarachniophytes. To determine why, and to address other fundamental questions about eukaryote-eukaryote endosymbiosis, we sequenced the nuclear genomes of the cryptophyte Guillardia theta and the chlorarachniophyte Bigelowiella natans. Both genomes have >21,000 protein genes and are intron rich, and B. natans exhibits unprecedented alternative splicing for a single-celled organism. Phylogenomic analyses and subcellular targeting predictions reveal extensive genetic and biochemical mosaicism, with both host- and endosymbiont-derived genes servicing the mitochondrion, the host cell cytosol, the plastid and the remnant endosymbiont cytosol of both algae. Mitochondrion-to-nucleus gene transfer still occurs in both organisms but plastid-to-nucleus and nucleomorph-to-nucleus transfers do not, which explains why a small residue of essential genes remains locked in each nucleomorph.

Patterns of 5' untranslated region length distribution in Encephalitozoon cuniculi: implications for gene regulation and potential links between transcription and splicing.

Encephalitozoon cuniculi, a eukaryotic intracellular parasite belonging to the group Microsporidia, has a highly reduced and compacted genome. Its mRNA transcripts have been found to differ between the two life stages, the spore and meront, of the parasite. Spore transcripts generally have more transcription start sites, longer 5' untranslated regions (UTRs), and overlap more frequently with upstream genes than those of meronts. A previous analysis of 31 meront gene transcripts showed that most have short 5'UTRs, and intron-containing genes, mostly ribosomal protein genes, exclusively have very short 5'UTRs. Here we analyzed a larger set of transcripts from meronts, and we find a pattern of 5'UTR length distribution similar to other reduced genomes. There is an abundance of very short 5'UTRs that are <20 bp in length, and very few 5'UTRs that are much longer. We also find a relationship between gene categories and 5'UTR length: intron-containing genes and ribosomal protein genes have exclusively short 5'UTRs. We suggest that the abundance of short 5'UTRs may be related to a class of highly expressed genes that benefit the parasite's growth cycle. Also, the longer 5'UTRs may be playing a role in down-regulating expression of genes that require temporal or environment-induced expression.

Splicing and transcription differ between spore and intracellular life stages in the parasitic microsporidia.

Microsporidia are a diverse group of highly derived fungal relatives that are intracellular parasites of many animals. Both transcription and introns have been shown to be unusual in microsporidia: The complete genome of the human parasite Encephalitozoon cuniculi has only a few very short introns, and two distantly related microsporidian spores have been shown to harbor transcripts encoding several genes that overlap on different strands. However, microsporidia alternate between two life stages: the intracellular proliferative stage and the extracellular and largely metabolically dormant infectious spore. To date, most studies have focused on the spore. Here, we have compared transcription profiles for a number of genes from both life stages of microsporidia and found major differences in both the prevalence of overlapping transcription and splicing. Specifically, spore transcripts in E. cuniculi have longer 5' untranslated regions, overlap more frequently with upstream genes, and have a significantly higher number of transcription initiation sites compared with intracellular transcripts from the same species. In addition, we demonstrate that splicing occurs exclusively in the intracellular stage and not in spore messenger RNAs (mRNAs) in both E. cuniculi and the distantly related Antonospora locustae. These differences between the microsporidian life stages raise questions about the functional importance of transcripts in the spore. We hypothesize that at least some transcripts in spores are a product of the cell's transition into a dormant state and that these unusual mRNAs could play a structural role rather than an informational one.