A functionally conserved STORR gene fusion in Papaver species that diverged 16.8 million years ago

The STORR gene fusion event is considered a key step in the evolution of benzylisoquinoline alkaloid (BIA) metabolism in opium poppy as the resulting bi-modular protein performs the isomerization of (S)- to (R)- reticuline which is required for morphinan biosynthesis. Our previous analysis of the opium poppy genome suggested the STORR gene fusion event occurred before a whole genome duplication event 7.2 million years ago. Here we use a combination of phylogenetic, transcriptomic, metabolomic, biochemical and genomic analysis to investigate the origin of the STORR gene fusion across the Papaveraceae family. The pro-morphinan/morphinan subclass of BIAs was present in a subset of 10 Papaver species including P. somniferum (opium poppy) and this correlated with the presence of the STORR gene fusion with one important exception. P. californicum does not produce morphinans but it does contain a STORR gene fusion that epimerizes (S)- to (R)- reticuline when heterologously expressed in yeast. The high similarity of the amino acid sequence linking the two modules of STORR along with phylogenetic gene tree analysis strongly suggests the gene fusion occurred only once and between 17-25 million years ago before the separation of P. californicum from the other Papaver species. We discovered that the most abundant BIA in P. californicum is (R)- glaucine, a member of the aporphine subclass of BIAs. Only the (S) isomer of this compound has previously been reported from nature. These results lead us to conclude that the function of the STORR gene fusion is not exclusive to morphinan production in the Papaveraceae.

Pan-genomic characterization of high-risk pediatric papillary thyroid carcinoma

Pediatric papillary thyroid carcinomas (pPTCs) are often indolent tumors with excellent long-term outcome, although subsets of cases are clinically troublesome and recur. Although it is generally thought to exhibit similar molecular aberrancies as their counterpart tumors in adults, the pan-genomic landscape of clinically aggressive pPTCs has not been previously described. In this study, five pairs of primary and synchronously metastatic pPTC from patients with high-risk phenotypes were characterized using parallel whole-genome and -transcriptome sequencing. Primary tumors and their metastatic components displayed an exceedingly low number of coding somatic mutations and gross chromosomal alterations overall, with surprisingly few shared mutational events. Two cases exhibited one established gene fusion event each (SQSTM1-NTRK3 and NCOA4-RET) in both primary and metastatic tissues, and one case each was positive for a BRAF V600E mutation and a germline truncating CHEK2 mutation, respectively. One single case was without apparent driver events and was considered as a genetic orphan. Non-coding mutations in cancer-associated regions were generally not present. By expressional analyses, fusion-driven primary and metastatic pPTC clustered separately from the mutation-driven cases and the sole genetic orphan. We conclude that pPTCs are genetically indolent tumors with exceedingly stable genomes. Several mutations found exclusively in the metastatic samples which may represent novel genetic events that drive the metastatic behavior, and the differences in mutational compositions suggest early clonal divergence between primary tumors and metastases. Moreover, an overrepresentation of mutational and expressional dysregulation of immune regulatory pathways was noted among fusion-positive pPTC metastases, suggesting that these tumors might facilitate spread through immune evasive mechanisms.

Detection of a Novel, and Likely Ancestral, Tn916-Like Element from a Human Saliva Metagenomic Library

Tn916 is a conjugative transposon (CTn) and the first reported and most well characterised of the Tn916/Tn1545 family of CTns. Tn916-like elements have a characteristic modular structure and different members of this family have been identified based on similarities and variations in these modules. In addition to carrying genes encoding proteins required for their conjugation, Tn916-like elements also carry accessory, antimicrobial resistance genes; most commonly the tetracycline resistance gene, tet(M). Our study aimed to identify and characterise tetracycline resistance genes from the human saliva metagenome using a functional metagenomic approach. We identified a tetracycline-resistant clone, TT31, the sequencing of which revealed it to encode both tet(M) and tet(L). Comparison of the TT31 sequence with the accessory, regulation, and recombination modules of other Tn916-like elements indicated that a partial Tn916-like element encoding a truncated orf9 was cloned in TT31. Analysis indicated that a previous insertion within the truncated orf9 created the full length orf9 found in most Tn916-like transposons; demonstrating that orf9 is, in fact, the result of a gene fusion event. Thus, we hypothesise that the Tn916-like element cloned in TT31 likely represents an ancestral Tn916.

154 Genomic Landscape of Radiation-Induced Meningiomas

INTRODUCTION Majority of pediatric cancers require the irradiation of the central nervous system (CNS), and as more patients survive into adulthood from improved oncological therapy the sequelae of brain radiation are increasing in prevalence. Radiation-induced meningiomas (RIMs), one such secondary effect, demonstrate a clinically more aggressive behaviour than sporadic meningiomas (SMs). We aimed to describe the genomic mutational landscape of RIMs METHODS We analyzed a principal cohort of 18 RIMs, with 31 RIMs overall, from patients who received childhood radiation therapy and 30 SMs, as a comparator population. We performed a multiplatform integrative genomic analysis; including methylation, whole exome and RNA sequencing. RESULTS >RIMs exhibited a five-fold increase in copy number alterations, commonly the loss of chromosome 1p (17/18 RIMs) and 22q (17/18 RIMs), which was significantly more than observed in sporadic meningiomas. Furthermore, RNA sequencing data revealed an NF2 gene fusion event in 35.3% of RIMs In all 6 cases, there was a complete NF2 exon spliced into a complete exon of a reciprocal gene, suggesting that the breakpoints of genomic rearrangement are intronic. All tumours with the NF2 fusion also possessed monosomy of chromosome 22q, rendering the cells with homozygous disruption of NF2. Clinically, RIMs with the NF2 fusion exhibited ill-defined borders and a tendency to develop in anatomic frontal location. The NF2 fusion RIMs, also, had a significantly faster growth rate compared to non-fusion RIMS (P < 0.05). Also, targeted sequencing panel confirmed that RIMs had fewer nonsynonymous NF2 mutations (6.5% vs. 30% in SM) and absence of mutations in TRAF7, SMO, KLF4, PIK3CA and AKT1, genes traditionally involved in SMs. CONCLUSION Our study demonstrates that RIMs have distinct genomic drivers of oncogenesis as compared to SMs, specifically NF2 inactivation through fusion event. Radiation therapy possibly triggers genomic structural rearrangements through error-prone repair of double-stranded DNA breaks.

Acquisition of genome-wide copy number alterations in monozygotic twins with acute lymphoblastic leukemia

Chimeric fusion genes are highly prevalent in childhood acute lymphoblastic leukemia (ALL) and are mostly prenatal, early genetic events in the evolutionary trajectory of this cancer. ETV6-RUNX1–positive ALL also has multiple (∼ 6 per case) copy number alterations (CNAs) as revealed by genome-wide single-nucleotide polymorphism arrays. Recurrent CNAs are probably “driver” events contributing critically to clonal diversification and selection, but at diagnosis, their developmental timing is “buried” in the leukemia's covert natural history. This conundrum can be resolved with twin pairs. We identified and compared CNAs in 5 pairs of monozygotic twins with concordant ETV6-RUNX1–positive ALL and 1 pair discordant for ETV6-RUNX1 positive ALL. We compared, within each pair, CNAs classified as potential “driver” or “passenger” mutations based upon recurrency and, where known, gene function. An average of 5.1 (range 3-11) CNAs (excluding immunoglobulin/T-cell receptor alterations) were identified per case. All “driver” CNAs (total of 32) were distinct within each of the 5 twin pairs with concordant ALL. “Driver” CNAs in another twin with ALL were all absent in the shared ETV6-RUNX1–positive preleukemic clone of her healthy co-twin. These data place all “driver” CNAs secondary to the prenatal gene fusion event and most probably postnatal in the sequential, molecular pathogenesis of ALL.

Evolutionary Origins of the Eukaryotic Shikimate Pathway: Gene Fusions, Horizontal Gene Transfer, and Endosymbiotic Replacements

ABSTRACT Currently the shikimate pathway is reported as a metabolic feature of prokaryotes, ascomycete fungi, apicomplexans, and plants. The plant shikimate pathway enzymes have similarities to prokaryote homologues and are largely active in chloroplasts, suggesting ancestry from the plastid progenitor genome. Toxoplasma gondii, which also possesses an alga-derived plastid organelle, encodes a shikimate pathway with similarities to ascomycete genes, including a five-enzyme pentafunctional arom. These data suggests that the shikimate pathway and the pentafunctional arom either had an ancient origin in the eukaryotes or was conveyed by eukaryote-to-eukaryote horizontal gene transfer (HGT). We expand sampling and analyses of the shikimate pathway genes to include the oomycetes, ciliates, diatoms, basidiomycetes, zygomycetes, and the green and red algae. Sequencing of cDNA from Tetrahymena thermophila confirmed the presence of a pentafused arom, as in fungi and T. gondii. Phylogenies and taxon distribution suggest that the arom gene fusion event may be an ancient eukaryotic innovation. Conversely, the Plantae lineage (represented here by both Viridaeplantae and the red algae) acquired different prokaryotic genes for all seven steps of the shikimate pathway. Two of the phylogenies suggest a derivation of the Plantae genes from the cyanobacterial plastid progenitor genome, but if the full Plantae pathway was originally of cyanobacterial origin, then the five other shikimate pathway genes were obtained from a minimum of two other eubacterial genomes. Thus, the phylogenies demonstrate both separate HGTs and shared derived HGTs within the Plantae clade either by primary HGT transfer or secondarily via the plastid progenitor genome. The shared derived characters support the holophyly of the Plantae lineage and a single ancestral primary plastid endosymbiosis. Our analyses also pinpoints a minimum of 50 gene/domain loss events, demonstrating that loss and replacement events have been an important process in eukaryote genome evolution.