Transcriptome-wide analysis of a baculovirus using nanopore sequencing

DEMINERS enables clinical metagenomics and comparative transcriptomic analysis by increasing throughput and accuracy of nanopore direct RNA sequencing

Nanopore direct RNA sequencing (DRS) is a powerful tool for RNA biology but suffers from low basecalling accuracy, low throughput, and high input requirements. We present DEMINERS, a novel DRS toolkit combining an RNA multiplexing workflow, a Random Forest-based barcode classifier, and an optimized convolutional neural network basecaller with species-specific training. DEMINERS enables accurate demultiplexing of up to 24 samples, reducing RNA input and runtime. Applications include clinical metagenomics, cancer transcriptomics, and parallel transcriptomic comparisons, uncovering microbial diversity in COVID-19 and m6A's role in malaria and glioma. DEMINERS offers a robust, high-throughput solution for precise transcript and RNA modification analysis.

KSHV 3.0: a state-of-the-art annotation of the Kaposi's sarcoma-associated herpesvirus transcriptome using cross-platform sequencing

Kaposi's sarcoma-associated herpesvirus (KSHV) is a large, oncogenic DNA virus belonging to the gammaherpesvirus subfamily. KSHV has been extensively studied with various high-throughput RNA-sequencing approaches to map the transcription start and end sites, the splice junctions, and the translation initiation sites. Despite these efforts, the comprehensive annotation of the viral transcriptome remains incomplete. In the present study, we generated a long-read sequencing data set of the lytic and latent KSHV transcriptome using native RNA and direct cDNA-sequencing methods. This was supplemented with Cap Analysis of Gene Expression sequencing based on a short-read platform. We also utilized data sets from previous publications for our analysis. As a result of this combined approach, we have identified a number of novel viral transcripts and RNA isoforms and have either corroborated or improved the annotation of previously identified viral RNA molecules, thereby notably enhancing our comprehension of the transcriptomic architecture of the KSHV genome. We also evaluated the coding capability of transcripts previously thought to be non-coding by integrating our data on the viral transcripts with translatomic information from other publications.IMPORTANCEDeciphering the viral transcriptome of Kaposi's sarcoma-associated herpesvirus is of great importance because we can gain insight into the molecular mechanism of viral replication and pathogenesis, which can help develop potential targets for antiviral interventions. Specifically, the identification of substantial transcriptional overlaps by this work suggests the existence of a genome-wide interference between transcriptional machineries. This finding indicates the presence of a novel regulatory layer, potentially controlling the expression of viral genes.

KSHV 3.0: A State-of-the-Art Annotation of the Kaposi's Sarcoma-Associated Herpesvirus Transcriptome Using Cross-Platform Sequencing

Kaposi's sarcoma-associated herpesvirus (KSHV) is a large, oncogenic DNA virus belonging to the gammaherpesvirus subfamily. KSHV has been extensively studied with various high-throughput RNA-sequencing approaches to map the transcription start and end sites, the splice junctions, and the translation initiation sites. Despite these efforts, the comprehensive annotation of the viral transcriptome remains incomplete. In the present study, we generated a long-read sequencing dataset of the lytic and latent KSHV transcriptome using native RNA and direct cDNA sequencing methods. This was supplemented with CAGE sequencing based on a short-read platform. We also utilized datasets from previous publications for our analysis. As a result of this combined approach, we have identified a number of novel viral transcripts and RNA isoforms and have either corroborated or improved the annotation of previously identified viral RNA molecules, thereby notably enhancing our comprehension of the transcriptomic architecture of the KSHV genome. We also evaluated the coding capability of transcripts previously thought to be non-coding, by integrating our data on the viral transcripts with translatomic information from other publications.

Discovering novel reproductive genes in a non-model fly using de novo GridION transcriptomics

Gene discovery has important implications for investigating phenotypic trait evolution, adaptation, and speciation. Male reproductive tissues, such as accessory glands (AGs), are hotspots for recruitment of novel genes that diverge rapidly even among closely related species/populations. These genes synthesize seminal fluid proteins that often affect post-copulatory sexual selection-they can mediate male-male sperm competition, ejaculate-female interactions that modify female remating and even influence reproductive incompatibilities among diverging species/populations. Although de novo transcriptomics has facilitated gene discovery in non-model organisms, reproductive gene discovery is still challenging without a reference database as they are often novel and bear no homology to known proteins. Here, we use reference-free GridION long-read transcriptomics, from Oxford Nanopore Technologies (ONT), to discover novel AG genes and characterize their expression in the widespread dung fly, Sepsis punctum. Despite stark population differences in male reproductive traits (e.g.: Body size, testes size, and sperm length) as well as female re-mating, the male AG genes and their secretions of S. punctum are still unknown. We implement a de novo ONT transcriptome pipeline incorporating quality-filtering and rigorous error-correction procedures, and we evaluate gene sequence and gene expression results against high-quality Illumina short-read data. We discover highly-expressed reproductive genes in AG transcriptomes of S. punctum consisting of 40 high-quality and high-confidence ONT genes that cross-verify against Illumina genes, among which 26 are novel and specific to S. punctum. Novel genes account for an average of 81% of total gene expression and may be functionally relevant in seminal fluid protein production. For instance, 80% of genes encoding secretory proteins account for 74% total gene expression. In addition, median sequence similarities of ONT nucleotide and protein sequences match within-Illumina sequence similarities. Read-count based expression quantification in ONT is congruent with Illumina's Transcript per Million (TPM), both in overall pattern and within functional categories. Rapid genomic innovation followed by recruitment of de novo genes for high expression in S. punctum AG tissue, a pattern observed in other insects, could be a likely mechanism of evolution of these genes. The study also demonstrates the feasibility of adapting ONT transcriptomics for gene discovery in non-model systems.

Baculovirus Display of Varicella-Zoster Virus Glycoprotein E Induces Robust Humoral and Cellular Immune Responses in Mice

Varicella-zoster virus (VZV) is the causative agent of varicella and herpes zoster (HZ) and can pose a significant challenge to human health globally. The initial VZV infection-more common in children-causes a self-limiting chicken pox. However, in later life, the latent VZV can become reactivated in these patients, causing HZ and postherpetic neuralgia (PHN), a serious and painful complication. VZV glycoprotein E (gE) has been developed into a licensed subunit vaccine against HZ (Shingrix). However, its efficacy relies on the concomitant delivery of a robust adjuvant (AS01B). Here, we sought to create a new immunogen for vaccine design by displaying the VZV-gE on the baculovirus surface (Bac-gE). Correct localization and display of gE on the engineered baculovirus was verified by flow cytometry and immune electron microscopy. We show that Bac-gE provides excellent antigenicity against VZV and induces not only stronger gE-specific CD4+ and CD8+ T cell responses but also higher levels of VZV-specific neutralizing antibodies as compared with other vaccine strategies in mice. Collectively, we show that the baculovirus display of VZV-gE confers ideal humoral and cellular immune responses required for HZ vaccine development, paving the way for a baculovirus-based vaccine design.

Dual isoform sequencing reveals complex transcriptomic and epitranscriptomic landscapes of a prototype baculovirus

In this study, two long-read sequencing (LRS) techniques, MinION from Oxford Nanopore Technologies and Sequel from the Pacific Biosciences, were used for the transcriptional characterization of a prototype baculovirus, Autographa californica multiple nucleopolyhedrovirus. LRS is able to read full-length RNA molecules, and thereby distinguish between transcript isoforms, mono- and polycistronic RNAs, and overlapping transcripts. Altogether, we detected 875 transcript species, of which 759 were novel and 116 were annotated previously. These RNA molecules include 41 novel putative protein coding transcripts [each containing 5'-truncated in-frame open reading frames (ORFs), 14 monocistronic transcripts, 99 polygenic RNAs, 101 non-coding RNAs, and 504 untranslated region isoforms. This work also identified novel replication origin-associated transcripts, upstream ORFs, cis-regulatory sequences and poly(A) sites. We also detected RNA methylation in 99 viral genes and RNA hyper-editing in the longer 5'-UTR transcript isoform of the canonical ORF 19 transcript.

MinION nanopore sequencing and assembly of a complete human papillomavirus genome

BACKGROUND

The MinION sequencer belongs to the third generation of sequencing technology that allows for the generation of ultra-long reads, representing a potentially more effective approach to characterize entire viral genome sequences than other time-consuming and low-throughput methodologies.

METHODS

We report the use of the MinION nanopore sequencer to sequence the full-length genome of human papillomavirus (HPV)-ICB2 (7441 bp), which was previously characterized in our laboratory. Three independent MinION libraries were prepared and sequenced using either three consecutive 12 -h runs (Protocol A) or a single run of 48 h starting from a pool of three barcoded DNA libraries (Protocol B). A fully automated bioinformatics pipeline was developed for the reconstruction of the viral genome.

RESULTS

Protocols A and B generated 9,354,933 and 3,255,879 reads, respectively. Read length N50 values ranged between 6976 and 7360 nucleotides over the four sequencing runs. Bioinformatics analysis showed that both protocols allowed for the reconstruction of the whole viral genome, with pairwise percentages of identity to HPV-ICB2 of 100 % for protocol A and 99.98 % for protocol B.

CONCLUSION

Our results show that the use of the MinION nanopore sequencer represents an effective strategy for whole-genome sequencing of HPVs with a minimal error rate.

TALC: Transcript-level Aware Long-read Correction

MOTIVATION

Long-read sequencing technologies are invaluable for determining complex RNA transcript architectures but are error-prone. Numerous 'hybrid correction' algorithms have been developed for genomic data that correct long reads by exploiting the accuracy and depth of short reads sequenced from the same sample. These algorithms are not suited for correcting more complex transcriptome sequencing data.

RESULTS

We have created a novel reference-free algorithm called Transcript-level Aware Long-Read Correction (TALC) which models changes in RNA expression and isoform representation in a weighted De Bruijn graph to correct long reads from transcriptome studies. We show that transcript-level aware correction by TALC improves the accuracy of the whole spectrum of downstream RNA-seq applications and is thus necessary for transcriptome analyses that use long read technology.

AVAILABILITY AND IMPLEMENTATION

TALC is implemented in C++ and available at https://github.com/lbroseus/TALC.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Entomopathogenic Viruses in the Neotropics: Current Status and Recently Discovered Species

The market for biological control of insect pests in the world and in Brazil has grown in recent years due to the unwanted ecological and human health impacts of chemical insecticides. Therefore, research on biological control agents for pest management has also increased. For instance, insect viruses have been used to protect crops and forests around the world for decades. Among insect viruses, the baculoviruses are the most studied and used viral biocontrol agent. More than 700 species of insects have been found to be naturally infected by baculoviruses, with 90% isolated from lepidopteran insects. In this review, some basic aspects of baculovirus infection in vivo and in vitro infection, gene content, viral replication will be discussed. Furthermore, we provide examples of the use of insect viruses for biological pest control and recently characterized baculoviruses in Brazil.

The how and why of lncRNA function: An innate immune perspective

Next-generation sequencing has provided a more complete picture of the composition of the human transcriptome indicating that much of the "blueprint" is a vastness of poorly understood non-protein-coding transcripts. This includes a newly identified class of genes called long noncoding RNAs (lncRNAs). The lack of sequence conservation for lncRNAs across species meant that their biological importance was initially met with some skepticism. LncRNAs mediate their functions through interactions with proteins, RNA, DNA, or a combination of these. Their functions can often be dictated by their localization, sequence, and/or secondary structure. Here we provide a review of the approaches typically adopted to study the complexity of these genes with an emphasis on recent discoveries within the innate immune field. Finally, we discuss the challenges, as well as the emergence of new technologies that will continue to move this field forward and provide greater insight into the biological importance of this class of genes. This article is part of a Special Issue entitled: ncRNA in control of gene expression edited by Kotb Abdelmohsen.