Single-Molecule RNA Sequencing Reveals IFNγ-Induced Differential Expression of Immune Escape Genes in Merkel Cell Polyomavirus–Positive MCC Cell Lines

Merkel cell carcinoma (MCC) is a rare and highly aggressive cancer, which is mainly caused by genomic integration of the Merkel cell polyomavirus and subsequent expression of a truncated form of its large T antigen. The resulting primary tumor is known to be immunogenic and under constant pressure to escape immune surveillance. Because interferon gamma (IFNγ), a key player of immune response, is secreted by many immune effector cells and has been shown to exert both anti-tumoral and pro-tumoral effects, we studied the transcriptomic response of MCC cells to IFNγ. In particular, immune modulatory effects that may help the tumor evade immune surveillance were of high interest to our investigation. The effect of IFNγ treatment on the transcriptomic program of three MCC cell lines (WaGa, MKL-1, and MKL-2) was analyzed using single-molecule sequencing via the Oxford Nanopore platform. A significant differential expression of several genes was detected across all three cell lines. Subsequent pathway analysis and manual annotation showed a clear upregulation of genes involved in the immune escape of tumor due to IFNγ treatment. The analysis of selected genes on protein level underlined our sequencing results. These findings contribute to a better understanding of immune escape of MCC and may help in clinical treatment of MCC patients. Furthermore, we demonstrate that single-molecule sequencing can be used to assess characteristics of large eukaryotic transcriptomes and thus contribute to a broader access to sequencing data in the community due to its low cost of entry.

A whole genome atlas of 81 Psilocybe genomes as a resource for psilocybin production.

The Psilocybe genus is well known for the synthesis of valuable psychoactive compounds such as Psilocybin, Psilocin, Baeocystin and Aeruginascin. The ubiquity of Psilocybin synthesis in Psilocybe has been attributed to a horizontal gene transfer mechanism of a ~20Kb gene cluster. A recently published highly contiguous reference genome derived from long read single molecule sequencing has underscored interesting variation in this Psilocybin synthesis gene cluster. This reference genome has also enabled the shotgun sequencing of spores from many Psilocybe strains to better catalog the genomic diversity in the Psilocybin synthesis pathway. Here we present the de novo assembly of 81 Psilocybe genomes compared to the P.envy reference genome. Surprisingly, the genomes of Psilocybe galindoi, Psilocybe tampanensis and Psilocybe azurescens lack sequence coverage over the previously described Psilocybin synthesis pathway but do demonstrate amino acid sequence homology to a less contiguous gene cluster and may illuminate the previously proposed evolution of psilocybin synthesis.

Contingency and selection in mitochondrial genome dynamics

Eukaryotic cells contain numerous copies of mitochondrial DNA (mtDNA), allowing for the coexistence of mutant and wild-type mtDNA in individual cells. The fate of mutant mtDNA depends on their relative replicative fitness within cells and the resulting cellular fitness within populations of cells. Yet the dynamics of the generation of mutant mtDNA and features that inform their fitness remain unaddressed. Here we utilize long read single-molecule sequencing to track mtDNA mutational trajectories in Saccharomyces cerevisiae. We show a previously unseen pattern that constrains subsequent excision events in mtDNA fragmentation. We also provide evidence for the generation of rare and contentious non-periodic mtDNA structures that lead to persistent diversity within individual cells. Finally, we show that measurements of relative fitness of mtDNA fit a phenomenological model that highlights important biophysical parameters governing mtDNA fitness. Altogether, our study provides techniques and insights into the dynamics of large structural changes in genomes that may be applicable in more complex organisms.

Genome resource for Elsinoë batatas, the causal agent of stem and foliage scab disease of sweet potato

Elsinoë batatas is a phytopathogenic fungus causing stem and foliage scab disease of sweet potato. At present, there is no reference genome available for E. batatas, limiting basic research for the pathogen. The present study applied the nanopore single molecule sequencing technology to sequence the E. batatas genome. This study thus reports the first high-quality genome sequence of E. batatas, with a total contig size of 26.49 Mb, 50.8% GC content and an N50 of 2,546,814bp. The sequences obtained serve as a reference for analysis of E. batatas isolates and provide a resource to better understand the biology of stem and foliage scab disease of sweet potato.

A whole genome atlas of 81 Psilocybe genomes as a resource for psilocybin production.

The Psilocybe genus is well known for the synthesis of valuable psychoactive compounds such as Psilocybin, Psilocin, Baeocystin and Aeruginascin. The ubiquity of Psilocybin synthesis in Psilocybe has been attributed to a horizontal gene transfer mechanism of a ~20Kb gene cassette. A recently published highly contiguous reference genome derived from long read single molecule sequencing has underscored interesting variation in this Psilocybin synthesis gene cassette. This reference genome has also enabled the shotgun sequencing of spores from many Psilocybe strains to better catalog the genomic diversity in the Psilocybin synthesis pathway. Here we present the de novo assembly of genomes of 81 Psilocybe genomes compared to the P.envy reference genome. Surprisingly, the genomes of Psilocybe galindoi, Psilocybe tampanensis and Psilocybe azurescens lack sequence coverage over the previously described Psilocybin synthesis pathway but do demonstrate amino acid sequence homology to an alternative pathway and may illuminate previously proposed convergent evolution of Psilocybin synthesis.

A high resolution single molecule sequencing-based Arabidopsis transcriptome using novel methods of Iso-seq analysis

Background Accurate and comprehensive annotation of transcript sequences is essential for transcript quantification and differential gene and transcript expression analysis. Single molecule long read sequencing technologies provide improved integrity of transcript structures including alternative splicing, and transcription start and polyadenylation sites. However, accuracy is significantly affected by sequencing errors, mRNA degradation or incomplete cDNA synthesis. Results We present a new and comprehensive Arabidopsis thaliana Reference Transcript Dataset 3 (AtRTD3). AtRTD3 contains over 160k transcripts - twice that of the best current Arabidopsis transcriptome and including over 1,500 novel genes. 79% of transcripts are from Iso-seq with accurately defined splice junctions and transcription start and end sites. We developed novel methods to determine splice junctions and transcription start and end sites accurately. Mis-match profiles around splice junctions provided a powerful feature to distinguish correct splice junctions and remove false splice junctions. Stratified approaches identified high confidence transcription start/end sites and removed fragmentary transcripts due to degradation. AtRTD3 is a major improvement over existing transcriptomes as demonstrated by analysis of an Arabidopsis cold response RNA-seq time-series. AtRTD3 provided higher resolution of transcript expression profiling and identified cold- and light-induced differential transcription start and polyadenylation site usage. Conclusions AtRTD3 is the most comprehensive Arabidopsis transcriptome currently available. It improves the precision of differential gene and transcript expression, differential alternative splicing, and transcription start/end site usage from RNA-seq data. The novel methods for identifying accurate splice junctions and transcription start/end sites are widely applicable and will improve single molecule sequencing analysis from any species.

lra: A long read aligner for sequences and contigs

It is computationally challenging to detect variation by aligning single-molecule sequencing (SMS) reads, or contigs from SMS assemblies. One approach to efficiently align SMS reads is sparse dynamic programming (SDP), where optimal chains of exact matches are found between the sequence and the genome. While straightforward implementations of SDP penalize gaps with a cost that is a linear function of gap length, biological variation is more accurately represented when gap cost is a concave function of gap length. We have developed a method, lra, that uses SDP with a concave-cost gap penalty, and used lra to align long-read sequences from PacBio and Oxford Nanopore (ONT) instruments as well as de novo assembly contigs. This alignment approach increases sensitivity and specificity for SV discovery, particularly for variants above 1kb and when discovering variation from ONT reads, while having runtime that are comparable (1.05-3.76×) to current methods. When applied to calling variation from de novo assembly contigs, there is a 3.2% increase in Truvari F1 score compared to minimap2+htsbox. lra is available in bioconda (https://anaconda.org/bioconda/lra) and github (https://github.com/ChaissonLab/LRA).

A draft reference assembly of the Psilocybe cubensis genome

We describe the use of high-fidelity single molecule sequencing to assemble the genome of the psychoactive Psilocybe cubensis mushroom. The genome is 46.6Mb, 46% GC, and in 32 contigs with an N50 of 3.3Mb. The BUSCO completeness scores are 97.6% with 1.2% duplicates. The Psilocybin synthesis cluster exists in a single 3.2Mb contig. The dataset is available from NCBI BioProject with accessions PRJNA687911 and PRJNA700437.