A community challenge to evaluate RNA-seq, fusion detection, and isoform quantification methods for cancer discovery.

The accurate identification and quantitation of RNA isoforms present in the cancer transcriptome is key for analyses ranging from the inference of the impacts of somatic variants to pathway analysis to biomarker development and subtype discovery. The ICGC-TCGA DREAM Somatic Mutation Calling in RNA (SMC-RNA) challenge was a crowd-sourced effort to benchmark methods for RNA isoform quantification and fusion detection from bulk cancer RNA sequencing (RNA-seq) data. It concluded in 2018 with a comparison of 77 fusion detection entries and 65 isoform quantification entries on 51 synthetic tumors and 32 cell lines with spiked-in fusion constructs. We report the entries used to build this benchmark, the leaderboard results, and the experimental features associated with the accurate prediction of RNA species. This challenge required submissions to be in the form of containerized workflows, meaning each of the entries described is easily reusable through CWL and Docker containers at https://github.com/SMC-RNA-challenge. A record of this paper's transparent peer review process is included in the supplemental information.

Full-length transcript characterization of SF3B1 mutation in chronic lymphocytic leukemia reveals downregulation of retained introns.

While splicing changes caused by somatic mutations in SF3B1 are known, identifying full-length isoform changes may better elucidate the functional consequences of these mutations. We report nanopore sequencing of full-length cDNA from CLL samples with and without SF3B1 mutation, as well as normal B cell samples, giving a total of 149 million pass reads. We present FLAIR (Full-Length Alternative Isoform analysis of RNA), a computational workflow to identify high-confidence transcripts, perform differential splicing event analysis, and differential isoform analysis. Using nanopore reads, we demonstrate differential 3' splice site changes associated with SF3B1 mutation, agreeing with previous studies. We also observe a strong downregulation of intron retention events associated with SF3B1 mutation. Full-length transcript analysis links multiple alternative splicing events together and allows for better estimates of the abundance of productive versus unproductive isoforms. Our work demonstrates the potential utility of nanopore sequencing for cancer and splicing research.

Responses of the picoprasinophyte Micromonas commoda to light and ultraviolet stress.

Micromonas is a unicellular marine green alga that thrives from tropical to polar ecosystems. We investigated the growth and cellular characteristics of acclimated mid-exponential phase Micromonas commoda RCC299 over multiple light levels and over the diel cycle (14:10 hour light:dark). We also exposed the light:dark acclimated M. commoda to experimental shifts from moderate to high light (HL), and to HL plus ultraviolet radiation (HL+UV), 4.5 hours into the light period. Cellular responses of this prasinophyte were quantified by flow cytometry and changes in gene expression by qPCR and RNA-seq. While proxies for chlorophyll a content and cell size exhibited similar diel variations in HL and controls, with progressive increases during day and decreases at night, both parameters sharply decreased after the HL+UV shift. Two distinct transcriptional responses were observed among chloroplast genes in the light shift experiments: i) expression of transcription and translation-related genes decreased over the time course, and this transition occurred earlier in treatments than controls; ii) expression of several photosystem I and II genes increased in HL relative to controls, as did the growth rate within the same diel period. However, expression of these genes decreased in HL+UV, likely as a photoprotective mechanism. RNA-seq also revealed two genes in the chloroplast genome, ycf2-like and ycf1-like, that had not previously been reported. The latter encodes the second largest chloroplast protein in Micromonas and has weak homology to plant Ycf1, an essential component of the plant protein translocon. Analysis of several nuclear genes showed that the expression of LHCSR2, which is involved in non-photochemical quenching, and five light-harvesting-like genes, increased 30 to >50-fold in HL+UV, but was largely unchanged in HL and controls. Under HL alone, a gene encoding a novel nitrite reductase fusion protein (NIRFU) increased, possibly reflecting enhanced N-assimilation under the 625 µmol photons m-2 s-1 supplied in the HL treatment. NIRFU's domain structure suggests it may have more efficient electron transfer than plant NIR proteins. Our analyses indicate that Micromonas can readily respond to abrupt environmental changes, such that strong photoinhibition was provoked by combined exposure to HL and UV, but a ca. 6-fold increase in light was stimulatory.

Identifying Aspects of the Post-Transcriptional Program Governing the Proteome of the Green Alga Micromonas pusilla.

Micromonas is a unicellular motile alga within the Prasinophyceae, a green algal group that is related to land plants. This picoeukaryote (<2 µm diameter) is widespread in the marine environment but is not well understood at the cellular level. Here, we examine shifts in mRNA and protein expression over the course of the day-night cycle using triplicated mid-exponential, nutrient replete cultures of Micromonas pusilla CCMP1545. Samples were collected at key transition points during the diel cycle for evaluation using high-throughput LC-MS proteomics. In conjunction, matched mRNA samples from the same time points were sequenced using pair-ended directional Illumina RNA-Seq to investigate the dynamics and relationship between the mRNA and protein expression programs of M. pusilla. Similar to a prior study of the marine cyanobacterium Prochlorococcus, we found significant divergence in the mRNA and proteomics expression dynamics in response to the light:dark cycle. Additionally, expressional responses of genes and the proteins they encoded could also be variable within the same metabolic pathway, such as we observed in the oxygenic photosynthesis pathway. A regression framework was used to predict protein levels from both mRNA expression and gene-specific sequence-based features. Several features in the genome sequence were found to influence protein abundance including codon usage as well as 3' UTR length and structure. Collectively, our studies provide insights into the regulation of the proteome over a diel cycle as well as the relationships between transcriptional and translational programs in the widespread marine green alga Micromonas.

Evidence-based green algal genomics reveals marine diversity and ancestral characteristics of land plants.

BACKGROUND: Prasinophytes are widespread marine green algae that are related to plants. Cellular abundance of the prasinophyte Micromonas has reportedly increased in the Arctic due to climate-induced changes. Thus, studies of these unicellular eukaryotes are important for marine ecology and for understanding Viridiplantae evolution and diversification. RESULTS: We generated evidence-based Micromonas gene models using proteomics and RNA-Seq to improve prasinophyte genomic resources. First, sequences of four chromosomes in the 22 Mb Micromonas pusilla (CCMP1545) genome were finished. Comparison with the finished 21 Mb genome of Micromonas commoda (RCC299; named herein) shows they share ≤8,141 of ~10,000 protein-encoding genes, depending on the analysis method. Unlike RCC299 and other sequenced eukaryotes, CCMP1545 has two abundant repetitive intron types and a high percent (26 %) GC splice donors. Micromonas has more genus-specific protein families (19 %) than other genome sequenced prasinophytes (11 %). Comparative analyses using predicted proteomes from other prasinophytes reveal proteins likely related to scale formation and ancestral photosynthesis. Our studies also indicate that peptidoglycan (PG) biosynthesis enzymes have been lost in multiple independent events in select prasinophytes and plants. However, CCMP1545, polar Micromonas CCMP2099 and prasinophytes from other classes retain the entire PG pathway, like moss and glaucophyte algae. Surprisingly, multiple vascular plants also have the PG pathway, except the Penicillin-Binding Protein, and share a unique bi-domain protein potentially associated with the pathway. Alongside Micromonas experiments using antibiotics that halt bacterial PG biosynthesis, the findings highlight unrecognized phylogenetic complexity in PG-pathway retention and implicate a role in chloroplast structure or division in several extant Viridiplantae lineages. CONCLUSIONS: Extensive differences in gene loss and architecture between related prasinophytes underscore their divergence. PG biosynthesis genes from the cyanobacterial endosymbiont that became the plastid, have been selectively retained in multiple plants and algae, implying a biological function. Our studies provide robust genomic resources for emerging model algae, advancing knowledge of marine phytoplankton and plant evolution.

Intron Invasions Trace Algal Speciation and Reveal Nearly Identical Arctic and Antarctic Micromonas Populations.

Spliceosomal introns are a hallmark of eukaryotic genes that are hypothesized to play important roles in genome evolution but have poorly understood origins. Although most introns lack sequence homology to each other, new families of spliceosomal introns that are repeated hundreds of times in individual genomes have recently been discovered in a few organisms. The prevalence and conservation of these introner elements (IEs) or introner-like elements in other taxa, as well as their evolutionary relationships to regular spliceosomal introns, are still unknown. Here, we systematically investigate introns in the widespread marine green alga Micromonas and report new families of IEs, numerous intron presence-absence polymorphisms, and potential intron insertion hot-spots. The new families enabled identification of conserved IE secondary structure features and establishment of a novel general model for repetitive intron proliferation across genomes. Despite shared secondary structure, the IE families from each Micromonas lineage bear no obvious sequence similarity to those in the other lineages, suggesting that their appearance is intimately linked with the process of speciation. Two of the new IE families come from an Arctic culture (Micromonas Clade E2) isolated from a polar region where abundance of this alga is increasing due to climate induced changes. The same two families were detected in metagenomic data from Antarctica--a system where Micromonas has never before been reported. Strikingly high identity between the Arctic isolate and Antarctic coding sequences that flank the IEs suggests connectivity between populations in the two polar systems that we postulate occurs through deep-sea currents. Recovery of Clade E2 sequences in North Atlantic Deep Waters beneath the Gulf Stream supports this hypothesis. Our research illuminates the dynamic relationships between an unusual class of repetitive introns, genome evolution, speciation, and global distribution of this sentinel marine alga.

Marine algae and land plants share conserved phytochrome signaling systems.

Phytochrome photosensors control a vast gene network in streptophyte plants, acting as master regulators of diverse growth and developmental processes throughout the life cycle. In contrast with their absence in known chlorophyte algal genomes and most sequenced prasinophyte algal genomes, a phytochrome is found in Micromonas pusilla, a widely distributed marine picoprasinophyte (<2 µm cell diameter). Together with phytochromes identified from other prasinophyte lineages, we establish that prasinophyte and streptophyte phytochromes share core light-input and signaling-output domain architectures except for the loss of C-terminal response regulator receiver domains in the streptophyte phytochrome lineage. Phylogenetic reconstructions robustly support the presence of phytochrome in the common progenitor of green algae and land plants. These analyses reveal a monophyletic clade containing streptophyte, prasinophyte, cryptophyte, and glaucophyte phytochromes implying an origin in the eukaryotic ancestor of the Archaeplastida. Transcriptomic measurements reveal diurnal regulation of phytochrome and bilin chromophore biosynthetic genes in Micromonas. Expression of these genes precedes both light-mediated phytochrome redistribution from the cytoplasm to the nucleus and increased expression of photosynthesis-associated genes. Prasinophyte phytochromes perceive wavelengths of light transmitted farther through seawater than the red/far-red light sensed by land plant phytochromes. Prasinophyte phytochromes also retain light-regulated histidine kinase activity lost in the streptophyte phytochrome lineage. Our studies demonstrate that light-mediated nuclear translocation of phytochrome predates the emergence of land plants and likely represents a widespread signaling mechanism in unicellular algae.

Estimating the nucleotide diversity in Ceratodon purpureus (Ditrichaceae) from 218 conserved exon-primed, intron-spanning nuclear loci.

PREMISE OF THE STUDY: We developed and tested primers for 218 nuclear loci for studying population genetics, phylogeography, and genome evolution in bryophytes. • METHODS AND RESULTS: We aligned expressed sequence tags (ESTs) from Ceratodon purpureus to the Physcomitrella patens genome sequence, and designed primers that are homologous to conserved exons but span introns in the P. patens genome. We tested these primers on four isolates from New York, USA; Otavalo, Ecuador; and two laboratory isolates from Austria (WT4 and GG1). The median genome-wide nucleotide diversity was 0.008 substitutions/site, but the range was large (0-0.14), illustrating the among-locus heterogeneity in the species. • CONCLUSIONS: These loci provide a valuable resource for finely resolved, genome-wide population genetic and species-level phylogenetic analyses of C. purpureus and its relatives.

The transcriptional diversity of 25 Drosophila cell lines.

Drosophila melanogaster cell lines are important resources for cell biologists. Here, we catalog the expression of exons, genes, and unannotated transcriptional signals for 25 lines. Unannotated transcription is substantial (typically 19% of euchromatic signal). Conservatively, we identify 1405 novel transcribed regions; 684 of these appear to be new exons of neighboring, often distant, genes. Sixty-four percent of genes are expressed detectably in at least one line, but only 21% are detected in all lines. Each cell line expresses, on average, 5885 genes, including a common set of 3109. Expression levels vary over several orders of magnitude. Major signaling pathways are well represented: most differentiation pathways are "off" and survival/growth pathways "on." Roughly 50% of the genes expressed by each line are not part of the common set, and these show considerable individuality. Thirty-one percent are expressed at a higher level in at least one cell line than in any single developmental stage, suggesting that each line is enriched for genes characteristic of small sets of cells. Most remarkable is that imaginal disc-derived lines can generally be assigned, on the basis of expression, to small territories within developing discs. These mappings reveal unexpected stability of even fine-grained spatial determination. No two cell lines show identical transcription factor expression. We conclude that each line has retained features of an individual founder cell superimposed on a common "cell line" gene expression pattern.

The developmental transcriptome of Drosophila melanogaster.

Drosophila melanogaster is one of the most well studied genetic model organisms; nonetheless, its genome still contains unannotated coding and non-coding genes, transcripts, exons and RNA editing sites. Full discovery and annotation are pre-requisites for understanding how the regulation of transcription, splicing and RNA editing directs the development of this complex organism. Here we used RNA-Seq, tiling microarrays and cDNA sequencing to explore the transcriptome in 30 distinct developmental stages. We identified 111,195 new elements, including thousands of genes, coding and non-coding transcripts, exons, splicing and editing events, and inferred protein isoforms that previously eluded discovery using established experimental, prediction and conservation-based approaches. These data substantially expand the number of known transcribed elements in the Drosophila genome and provide a high-resolution view of transcriptome dynamics throughout development.

Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation.

High-throughput mRNA sequencing (RNA-Seq) promises simultaneous transcript discovery and abundance estimation. However, this would require algorithms that are not restricted by prior gene annotations and that account for alternative transcription and splicing. Here we introduce such algorithms in an open-source software program called Cufflinks. To test Cufflinks, we sequenced and analyzed >430 million paired 75-bp RNA-Seq reads from a mouse myoblast cell line over a differentiation time series. We detected 13,692 known transcripts and 3,724 previously unannotated ones, 62% of which are supported by independent expression data or by homologous genes in other species. Over the time series, 330 genes showed complete switches in the dominant transcription start site (TSS) or splice isoform, and we observed more subtle shifts in 1,304 other genes. These results suggest that Cufflinks can illuminate the substantial regulatory flexibility and complexity in even this well-studied model of muscle development and that it can improve transcriptome-based genome annotation.

The completion of the Mammalian Gene Collection (MGC).

Since its start, the Mammalian Gene Collection (MGC) has sought to provide at least one full-protein-coding sequence cDNA clone for every human and mouse gene with a RefSeq transcript, and at least 6200 rat genes. The MGC cloning effort initially relied on random expressed sequence tag screening of cDNA libraries. Here, we summarize our recent progress using directed RT-PCR cloning and DNA synthesis. The MGC now contains clones with the entire protein-coding sequence for 92% of human and 89% of mouse genes with curated RefSeq (NM-accession) transcripts, and for 97% of human and 96% of mouse genes with curated RefSeq transcripts that have one or more PubMed publications, in addition to clones for more than 6300 rat genes. These high-quality MGC clones and their sequences are accessible without restriction to researchers worldwide.

Using N-SCAN or TWINSCAN to predict gene structures in genomic DNA sequences.

N-SCAN is a gene-prediction system that combines the methods of ab initio predictors like GENSCAN with information derived from genome comparison. It is the latest in the TWINSCAN series of programs. This unit describes the use of N-SCAN to identify gene structures in eukaryotic genomic sequences. Protocols for using N-SCAN through its Web interface and from the command line in a Linux environment are provided. Detailed discussion about the appropriate parameter settings, input-sequence processing, and choice of genome for comparison are included.

Iterative gene prediction and pseudogene removal improves genome annotation.

Correct gene prediction is impaired by the presence of processed pseudogenes: nonfunctional, intronless copies of real genes found elsewhere in the genome. Gene prediction programs frequently mistake processed pseudogenes for real genes or exons, leading to biologically irrelevant gene predictions. While methods exist to identify processed pseudogenes in genomes, no attempt has been made to integrate pseudogene removal with gene prediction, or even to provide a freestanding tool that identifies such erroneous gene predictions. We have created PPFINDER (for Processed Pseudogene finder), a program that integrates several methods of processed pseudogene finding in mammalian gene annotations. We used PPFINDER to remove pseudogenes from N-SCAN gene predictions, and show that gene prediction improves substantially when gene prediction and pseudogene masking are interleaved. In addition, we used PPFINDER with gene predictions as a parent database, eliminating the need for libraries of known genes. This allows us to run the gene prediction/PPFINDER procedure on newly sequenced genomes for which few genes are known.

The PCR suite.

The web application PCR Suite is an extension of the primer design program Primer3. It allows the design of primer sets encompassing single nucleotide polymorphisms, all exons of a single gene, all open reading frames in a list of cDNAs or the creation of overlapping PCR products.

Mutations in the DJ-1 gene associated with autosomal recessive early-onset parkinsonism.

The DJ-1 gene encodes a ubiquitous, highly conserved protein. Here, we show that DJ-1 mutations are associated with PARK7, a monogenic form of human parkinsonism. The function of the DJ-1 protein remains unknown, but evidence suggests its involvement in the oxidative stress response. Our findings indicate that loss of DJ-1 function leads to neurodegeneration. Elucidating the physiological role of DJ-1 protein may promote understanding of the mechanisms of brain neuronal maintenance and pathogenesis of Parkinson's disease.

Disruption of a long-range cis-acting regulator for Shh causes preaxial polydactyly.

Preaxial polydactyly (PPD) is a common limb malformation in human. A number of polydactylous mouse mutants indicate that misexpression of Shh is a common requirement for generating extra digits. Here we identify a translocation breakpoint in a PPD patient and a transgenic insertion site in the polydactylous mouse mutant sasquatch (Ssq). The genetic lesions in both lie within the same respective intron of the LMBR1/Lmbr1 gene, which resides approximately 1 Mb away from Shh. Genetic analysis of Ssq reveals that the Lmbr1 gene is incidental to the phenotype and that the mutation directly interrupts a cis-acting regulator of Shh. This regulator is most likely the target for generating PPD mutations in human.

A double RING-H2 domain in RNF32, a gene expressed during sperm formation.

The RING domain is a cysteine-rich zinc-binding motif, which is found in a wide variety of proteins, among which are several proto-oncogenes and the gene implicated in autosomal recessive juvenile parkinsonism, Parkin. The domain mediates binding to other proteins, either via their RING domains or other motifs. In several proteins, RING domains are found in combination with other cysteine-rich binding motifs and some proteins contain two RING domains. Recent evidence suggests that RING finger proteins function in the ubiquitin pathway as E3 ligases. A variant of the RING domain is the RING-H2 domain, in which one of the cysteines is replaced by a histidine. We have cloned and characterized a novel gene, RNF32, located on chromosome 7q36. RNF32 is contained in 37 kb of genomic DNA and consists of 9 constitutive and 8 alternatively spliced exons, most of which are alternative first exons. A long and a short transcript of the gene are expressed; the short transcript containing exons 1-4 only. This gene encodes two RING-H2 domains separated by an IQ domain of unknown function. This is the first reported gene with a double RING-H2 domain. In humans, RNF32 overlaps with a processed retroposon located on the opposite strand, C7orf13. RNF32 is specifically expressed in testis and ovary, whereas C7orf13 is testis-specific, suggesting that its expression may be regulated by elements in the RNF32 promoter region. RNF32 is expressed during spermatogenesis, most likely in spermatocytes and/or in spermatids, suggesting a possible role in sperm formation.