Long-Read Sequencing Revealed an Extensive Transcript Complexity in Herpesviruses

EST-SSR Marker Development and Full-Length Transcriptome Sequence Analysis of Tiger Lily (Lilium lancifolium Thunb)

The fast advancement and deployment of sequencing technologies after the Human Genome Project have greatly increased our knowledge of the eukaryotic genome sequences. However, due to technological concerns, high-quality genomic data has been confined to a few key organisms. Moreover, our understanding of which portions of genomes make up genes and which transcript isoforms synthesize these genes is scarce. Therefore, the current study has been designed to explore the reliability of the tiger lily (Lilium lancifolium Thunb) transcriptome. The PacBio-SMRT was used for attaining the complete transcriptomic profile. We obtained a total of 815,624 CCS (Circular Consensus Sequence) reads with an average length of 1295 bp. The tiger lily transcriptome has been sequenced for the first time using third-generation long-read technology. Furthermore, unigenes (38,707), lncRNAs (6852), and TF members (768) were determined based on the transcriptome data, followed by evaluating SSRs (3319). It has also been revealed that 105 out of 128 primer pairs effectively amplified PCR products. Around 15,608 transcripts were allocated to 25 distinct KOG Clusters, and 10,706 unigenes were grouped into 52 functional categories in the annotated transcripts. Until now, no tiger lily lncRNAs have been discovered. Results of this study may serve as an extensive set of reference transcripts and help us learn more about the transcriptomes of tiger lilies and pave the path for further research.

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.

Combined Short and Long-Read Sequencing Reveals a Complex Transcriptomic Architecture of African Swine Fever Virus

African swine fever virus (ASFV) is a large DNA virus belonging to the Asfarviridae family. Despite its agricultural importance, little is known about the fundamental molecular mechanisms of this pathogen. Short-read sequencing (SRS) can produce a huge amount of high-precision sequencing reads for transcriptomic profiling, but it is inefficient for comprehensively annotating transcriptomes. Long-read sequencing (LRS) can overcome some of SRS's limitations, but it also has drawbacks, such as low-coverage and high error rate. The limitations of the two approaches can be surmounted by the combined use of these techniques. In this study, we used Illumina SRS and Oxford Nanopore Technologies LRS platforms with multiple library preparation methods (amplified and direct cDNA sequencings and native RNA sequencing) for constructing the ASFV transcriptomic atlas. This work identified many novel transcripts and transcript isoforms and annotated the precise termini of previously described RNAs. This study identified a novel species of ASFV transcripts, the replication origin-associated RNAs. Additionally, we discovered several nested genes embedded into larger canonical genes. In contrast to the current view that the ASFV transcripts are monocistronic, we detected a significant extent of polycistronism. A multifaceted meshwork of transcriptional overlaps was also discovered.

New approach of genetic characterization of group A rotaviruses by the nanopore sequencing method

Nanopore sequencing of virus genomes represented by segmented RNA (e.g. rotaviruses) requires the development of specific approaches. Due to the massive use of rotavirus vaccines, the relevance of monitoring the genetic diversity of circulating strains of group A rotaviruses (RVA) increased. The WHO recommended method of multiplex type-specific PCR does not allow genotyping of all clinically significant strains of RVA and identifying inter-strain differences within the genotype. We have described a new principle of amplification of RVA gene segments using six primers for reverse transcription and one universal primer for PCR for nanopore sequencing. The amplification of RVA genome was tested on clinical samples and three phylogenetically distant laboratory RVA strains, Wa (G1P[8]), DS-1 (G2P[4]) and 568 (G3P[3]). The developed protocol of sample preparation and nanopore sequencing allowed obtaining full-length sequences for gene segments of RVA, including the diagnostically significant segments 9 (VP7), 4 (VP4) and 6 (VP6) with high accuracy and coverage. The accuracy of sequencing of the rotavirus genome exceeded 99.5 %, and the genome coverage varied for different strains from 59.0 to 99.6 % (on average 86 %). The developed approach of nanopore sequencing of RVA genome could be a prospective tool for epidemiological studies and surveillance of rotavirus infection.

Long-read assays shed new light on the transcriptome complexity of a viral pathogen

Characterization of global transcriptomes using conventional short-read sequencing is challenging due to the insensitivity of these platforms to transcripts isoforms, multigenic RNA molecules, and transcriptional overlaps. Long-read sequencing (LRS) can overcome these limitations by reading full-length transcripts. Employment of these technologies has led to the redefinition of transcriptional complexities in reported organisms. In this study, we applied LRS platforms from Pacific Biosciences and Oxford Nanopore Technologies to profile the vaccinia virus (VACV) transcriptome. We performed cDNA and direct RNA sequencing analyses and revealed an extremely complex transcriptional landscape of this virus. In particular, VACV genes produce large numbers of transcript isoforms that vary in their start and termination sites. A significant fraction of VACV transcripts start or end within coding regions of neighbouring genes. This study provides new insights into the transcriptomic profile of this viral pathogen.

Meta-analytic approach for transcriptome profiling of herpes simplex virus type 1

In this meta-analysis, we re-analysed and compared herpes simplex virus type 1 transcriptomic data generated by eight studies using various short- and long-read sequencing techniques and different library preparation methods. We identified a large number of novel mRNAs, non-coding RNAs and transcript isoforms, and validated many previously published transcripts. Here, we present the most complete HSV-1 transcriptome to date. Furthermore, we also demonstrate that various sequencing techniques, including both cDNA and direct RNA sequencing approaches, are error-prone, which can be circumvented by using integrated approaches. This work draws attention to the need for using multiple sequencing approaches and meta-analyses in transcriptome profiling studies to obtain reliable results.

The Progress of Multi-Omics Technologies: Determining Function in Lactic Acid Bacteria Using a Systems Level Approach

Lactic Acid Bacteria (LAB) have long been recognized as having a significant impact ranging from commercial to health domains. A vast amount of research has been carried out on these microbes, deciphering many of the pathways and components responsible for these desirable effects. However, a large proportion of this functional information has been derived from a reductionist approach working with pure culture strains. This provides limited insight into understanding the impact of LAB within intricate systems such as the gut microbiome or multi strain starter cultures. Whole genome sequencing of strains and shotgun metagenomics of entire systems are powerful techniques that are currently widely used to decipher function in microbes, but they also have their limitations. An available genome or metagenome can provide an image of what a strain or microbiome, respectively, is potentially capable of and the functions that they may carry out. A top-down, multi-omics approach has the power to resolve the functional potential of an ecosystem into an image of what is being expressed, translated and produced. With this image, it is possible to see the real functions that members of a system are performing and allow more accurate and impactful predictions of the effects of these microorganisms. This review will discuss how technological advances have the potential to increase the yield of information from genomics, transcriptomics, proteomics and metabolomics. The potential for integrated omics to resolve the role of LAB in complex systems will also be assessed. Finally, the current software approaches for managing these omics data sets will be discussed.

Physiological roles of antisense RNAs in prokaryotes

Prokaryotes encounter constant and often brutal modifications to their environment. In order to survive, they need to maintain fitness, which includes adapting their protein expression patterns. Many factors control gene expression but this review focuses on just one, namely antisense RNAs (asRNAs), a class of non-coding RNAs (ncRNAs) characterized by their location in cis and their perfect complementarity with their targets. asRNAs were considered for a long time to be trivial and only to be found on mobile genetic elements. However, recent advances in methodology have revealed that their abundance and potential activities have been underestimated. This review aims to illustrate the role of asRNA in various physiologically crucial functions in both archaea and bacteria, which can be regrouped in three categories: cell maintenance, horizontal gene transfer and virulence. A literature survey of asRNAs demonstrates the difficulties to characterize and assign a role to asRNAs. With the aim of facilitating this task, we describe recent technological advances that could be of interest to identify new asRNAs and to discover their function.

MinION sequencing to genotype US strains of infectious laryngotracheitis virus

Over the last decade the US broiler industry has fought long-lasting outbreaks of infectious laryngotracheitis (ILTV). Previously, nine genotypes (I-IX) of ILTVs have been recognized using the polymerase chain reaction-restriction fragment length polymorphisms (PCR-RFLP) method with three viral alleles (gB, gM and UL47/gG). In this study, the genotyping system was simplified to six genotypes by amplicon sequencing and examining discriminating single nucleotide polymorphisms (SNPs) within these open reading frames. Using phylogenomic analysis of 27 full genomes of ILTV, a single allele (ORF A/ORF B) was identified containing SNPs that could differentiate ILTVs into genotypes congruent with the phylogenetic partitioning. The allelic variations allowed for the cataloging of the 27 strains into 5 genotypes: vaccinal TCO, vaccinal CEO, virulent CEO-like, virulent US and virulent US backyard flocks from 1980 to 1990, correlating with the PCR-RFLP genotypes I/ II/ III (TCO), IV (CEO), V (virulent CEO-like), VI (virulent US) and VII/VIII/IX (virulent US backyard flock isolates). With the unique capabilities of third generation sequencing, we investigated the application of Oxford Nanopore MinION technology for rapid sequencing of the amplicons generated in the single-allele assay. This technology was an improvement over Sanger-based sequencing of the single allele amplicons due to a booster amplification step in the MinION sequencing protocol. Overall, there was a 90% correlation between the genotyping results of the single-allele assay and the multi-allele assay. Surveillance of emerging ILTV strains could greatly benefit from real-time amplicon sequencing using the single-allele assay and MinION sequencing. RESEARCH HIGHLIGHTS A multi-allelic assay identified nine ILTV genotypes circulating in the US Single-allele genotyping is congruent with whole genome phylogenetic partitioning US ILTV strains can be grouped into five genotypes using the single-allele assay The single-allele assay can be done using MinION sequencing of barcoded amplicons.

RNA splicing analysis in genomic medicine

High-throughput next-generation sequencing technologies have led to a rapid increase in the number of sequence variants identified in clinical practice via diagnostic genetic tests. Current bioinformatic analysis pipelines fail to take adequate account of the possible splicing effects of such variants, particularly where variants fall outwith canonical splice site sequences, and consequently the pathogenicity of such variants may often be missed. The regulation of splicing is highly complex and as a result, in silico prediction tools lack sufficient sensitivity and specificity for reliable use. Variants of all kinds can be linked to aberrant splicing in disease and the need for correct identification and diagnosis grows ever more crucial as novel splice-switching antisense oligonucleotide therapies start to enter clinical usage. RT-PCR provides a useful targeted assay of the splicing effects of identified variants, while minigene assays, massive parallel reporter assays and animal models can also be used for more detailed study of a particular splicing system, given enough time and resources. However, RNA-sequencing (RNA-seq) has the potential to be used as a rapid diagnostic tool in genomic medicine. By utilising data science approaches and machine learning, it may prove possible to finally understand and interpret the 'splicing code' and apply this knowledge in human disease diagnostics.

Proteogenomic Identification of a Novel Protein-Encoding Gene in Bovine Herpesvirus 1 That Is Expressed during Productive Infection

Bovine herpesvirus 1 (BoHV-1) is one of several microbes that contributes to the development of the bovine respiratory disease (BRD) and can also induce abortions in cattle. As other alpha-herpesvirinae subfamily members, BoHV-1 efficiently replicates in many cell types and subsequently establishes a life-long latent infection in sensory neurons. BoHV-1 encodes more than 70 proteins that are expressed in a well-defined manner during productive infection. However, in silico open reading frame (ORF) prediction of the BoHV-1 genome suggests that the virus may encode more than one hundred proteins. In this study we used mass spectrometry followed by proteogenomic mapping to reveal the existence of 92 peptides that map to previously un-annotated regions of the viral genome. Twenty-one of the newly termed “intergenic peptides” were predicted to have a viable ORF around them. Twelve of these produced an mRNA transcript as demonstrated by strand-specific RT-PCR. We further characterized the 5′ and 3′ termini of one mRNA transcript, ORF-A, and detected a 55 kDa protein produced during active infection using a custom-synthesized antibody. We conclude that the coding potential of BoHV-1 is underestimated.