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

Exploring the transcriptomic profile of human monkeypox virus via CAGE and native RNA sequencing approaches

In this study, we employed short- and long-read sequencing technologies to delineate the transcriptional architecture of the human monkeypox virus and to identify key regulatory elements that govern its gene expression. Specifically, we conducted a transcriptomic analysis to annotate the transcription start sites (TSSs) and transcription end sites (TESs) of the virus by utilizing Cap Analysis of gene expression sequencing on the Illumina platform and direct RNA sequencing on the Oxford Nanopore technology device. Our investigations uncovered significant complexity in the use of alternative TSSs and TESs in viral genes. In this research, we also detected the promoter elements and poly(A) signals associated with the viral genes. Additionally, we identified novel genes in both the left and right variable regions of the viral genome.IMPORTANCEGenerally, gaining insight into how the transcription of a virus is regulated offers insights into the key mechanisms that control its life cycle. The recent outbreak of the human monkeypox virus has underscored the necessity of understanding the basic biology of its causative agent. Our results are pivotal for constructing a comprehensive transcriptomic atlas of the human monkeypox virus, providing valuable resources for future studies.

African Swine Fever Virus Protein-Protein Interaction Prediction

The African swine fever virus (ASFV) is an often deadly disease in swine and poses a threat to swine livestock and swine producers. With its complex genome containing more than 150 coding regions, developing effective vaccines for this virus remains a challenge due to a lack of basic knowledge about viral protein function and protein-protein interactions between viral proteins and between viral and host proteins. In this work, we identified ASFV-ASFV protein-protein interactions (PPIs) using artificial intelligence-powered protein structure prediction tools. We benchmarked our PPI identification workflow on the Vaccinia virus, a widely studied nucleocytoplasmic large DNA virus, and found that it could identify gold-standard PPIs that have been validated in vitro in a genome-wide computational screening. We applied this workflow to more than 18,000 pairwise combinations of ASFV proteins and were able to identify seventeen novel PPIs, many of which have corroborating experimental or bioinformatic evidence for their protein-protein interactions, further validating their relevance. Two protein-protein interactions, I267L and I8L, I267L__I8L, and B175L and DP79L, B175L__DP79L, are novel PPIs involving viral proteins known to modulate host immune response.

Bacteriophages avoid autoimmunity from cognate immune systems as an intrinsic part of their life cycles

Dormant prophages protect lysogenic cells by expressing diverse immune systems, which must avoid targeting their cognate prophages upon activation. Here we report that multiple Staphylococcus aureus prophages encode Tha (tail-activated, HEPN (higher eukaryotes and prokaryotes nucleotide-binding) domain-containing anti-phage system), a defence system activated by structural tail proteins of incoming phages. We demonstrate the function of two Tha systems, Tha-1 and Tha-2, activated by distinct tail proteins. Interestingly, Tha systems can also block reproduction of the induced tha-positive prophages. To prevent autoimmunity after prophage induction, these systems are inhibited by the product of a small overlapping antisense gene previously believed to encode an excisionase. This genetic organization, conserved in S. aureus prophages, allows Tha systems to protect prophages and their bacterial hosts against phage predation and to be turned off during prophage induction, balancing immunity and autoimmunity. Our results show that the fine regulation of these processes is essential for the correct development of prophages' life cycle.

Transcription termination and readthrough in African swine fever virus

Introduction

African swine fever virus (ASFV) is a nucleocytoplasmic large DNA virus (NCLDV) that encodes its own host-like RNA polymerase (RNAP) and factors required to produce mature mRNA. The formation of accurate mRNA 3' ends by ASFV RNAP depends on transcription termination, likely enabled by a combination of sequence motifs and transcription factors, although these are poorly understood. The termination of any RNAP is rarely 100% efficient, and the transcriptional "readthrough" at terminators can generate long mRNAs which may interfere with the expression of downstream genes. ASFV transcriptome analyses reveal a landscape of heterogeneous mRNA 3' termini, likely a combination of bona fide termination sites and the result of mRNA degradation and processing. While short-read sequencing (SRS) like 3' RNA-seq indicates an accumulation of mRNA 3' ends at specific sites, it cannot inform about which promoters and transcription start sites (TSSs) directed their synthesis, i.e., information about the complete and unprocessed mRNAs at nucleotide resolution.

Methods

Here, we report a rigorous analysis of full-length ASFV transcripts using long-read sequencing (LRS). We systematically compared transcription termination sites predicted from SRS 3' RNA-seq with 3' ends mapped by LRS during early and late infection.

Results

Using in-vitro transcription assays, we show that recombinant ASFV RNAP terminates transcription at polyT stretches in the non-template strand, similar to the archaeal RNAP or eukaryotic RNAPIII, unaided by secondary RNA structures or predicted viral termination factors. Our results cement this T-rich motif (U-rich in the RNA) as a universal transcription termination signal in ASFV. Many genes share the usage of the same terminators, while genes can also use a range of terminators to generate transcript isoforms varying enormously in length. A key factor in the latter phenomenon is the highly abundant terminator readthrough we observed, which is more prevalent during late compared with early infection.

Discussion

This indicates that ASFV mRNAs under the control of late gene promoters utilize different termination mechanisms and factors to early promoters and/or that cellular factors influence the viral transcriptome landscape differently during the late stages of infection.

The transcriptomic insight into the differential susceptibility of African Swine Fever in inbred pigs

African swine fever (ASF) is a global threat to animal health and food security. ASF is typically controlled by strict biosecurity, rapid diagnosis, and culling of affected herds. Much progress has been made in developing modified live virus vaccines against ASF. There is host variation in response to ASF infection in the field and under controlled conditions. To better understand the dynamics underlying this host differential morbidity, whole transcriptome profiling was carried out in twelve immunized and five sham immunized pigs. Seventeen MHC homozygous inbred Large white Babraham pigs were sampled at three time points before and after the challenge. The changes in the transcriptome profiles of infected animals were surveyed over time. In addition, the immunization effect on the host response was studied as well among the contrasts of all protection subgroups. The results showed two promising candidate genes to distinguish between recovered and non-recovered pigs after infection with a virulent African swine fever virus (ASFV) pre-infection: HTRA3 and GFPT2 (padj < 0.05). Variant calling on the transcriptome assemblies showed a two-base pair insertion into the ACOX3 gene closely located to HTRA3 that may regulate its expression as a putative genomic variant for ASF. Several significant DGEs, enriched gene ontology (GO) terms, and KEGG pathways at 1 day and 7 days post-infection, compared to the pre-infection, indicate a significant inflammation response immediately after ASF infection. The presence of the virus was confirmed by the mapping of RNA-Seq reads on two whole viral genome sequences. This was concordant with a higher virus load in the non-recovered animals 7 days post-infection. There was no transcriptome signature on the immunization at pre-infection and 1 day post-infection. More samples and data from additional clinical trials may support these findings.

In-depth Temporal Transcriptome Profiling of Monkeypox and Host Cells using Nanopore Sequencing

The recent human Monkeypox outbreak underlined the importance of studying basic biology of orthopoxviruses. However, the transcriptome of its causative agent has not been investigated before neither with short-, nor with long-read sequencing approaches. This Oxford Nanopore long-read RNA-Sequencing dataset fills this gap. It will enable the in-depth characterization of the transcriptomic architecture of the monkeypox virus, and may even make possible to annotate novel host transcripts. Moreover, our direct cDNA and native RNA sequencing reads will allow the estimation of gene expression changes of both the virus and the host cells during the infection. Overall, our study will lead to a deeper understanding of the alterations caused by the viral infection on a transcriptome level.

Metagenomics and metabarcoding experimental choices and their impact on microbial community characterization in freshwater recirculating aquaculture systems

BACKGROUND

Microbial communities in recirculating aquaculture systems (RAS) play a role in system success, nutrient cycling, and water quality. Considering the increasing socio-economic role of fish farming, e.g., regarding food security, an in-depth understanding of aquaculture microbial communities is also relevant from a management perspective, especially regarding the growth, development, and welfare of the farmed animal. However, the current data on the composition of microbial communities within RAS is patchy, which is partly attributable to diverging method choices that render comparative analyses challenging. Therefore, there is a need for accurate, standardized, and user-friendly methods to study microbial communities in aquaculture systems.

RESULTS

We compared sequencing approach performances (3 types of 16S short amplicon sequencing, PacBio long-read amplicon sequencing, and amplification-free shotgun metagenomics) in the characterization of microbial communities in two commercial RAS fish farms. Results showed that 16S primer choice and amplicon length affect some values (e.g., diversity measures, number of assigned taxa or distinguishing ASVs) but have no impact on spatio-temporal patterns between sample types, farms and time points. This implies that 16S rRNA approaches are adequate for community studies. The long-read amplicons underperformed regarding the quantitative resolution of spatio-temporal patterns but were suited to identify functional services, e.g., nitrification cycling and the detection of pathogens. Finally, shotgun metagenomics extended the picture to fungi, viruses, and bacteriophages, opening avenues for exploring inter-domain interactions. All sequencing datasets agreed on major prokaryotic players, such as Actinobacteriota, Bacteroidota, Nitrospirota, and Proteobacteria.

CONCLUSION

The different sequencing approaches yielded overlapping and highly complementary results, with each contributing unique data not obtainable with the other approaches. We conclude that a tiered approach constitutes a strategy for obtaining the maximum amount of information on aquaculture microbial communities and can inform basic research on community evolution dynamics. For specific and/or applied questions, single-method approaches are more practical and cost-effective and could lead to better farm management practices.

Analysis of gene expression in monocytes of immunized pigs after infection with homologous or heterologous African swine fever virus

African swine fever is a deadly disease of pigs caused by the large DNA virus (ASFV). Despite intensive research, little is known about the molecular mechanisms of ASFV pathogenesis. Transcriptome analysis of host and viral genes in infected macrophages revealed changes in expression of genes involved in various biological processes, including immune response, inflammatory response and apoptosis. To understand the mechanisms of virus pathogenesis, we used transcriptome analysis to identify the differences in gene expression between peripheral blood monocytes (PBMCs) isolated from pigs immunized with attenuated Congo ASFV strain (KK262), and then infected in vitro with virulent homologous Congo strain (K49) or heterologous Mozambique strain (M78). We found that overexpression of IFN-γ was detected only in cells infected with M78, although the expression of interferon-stimulated genes was increased in both types of cells. In addition, up-regulation of pro-inflammatory cytokines and chemokines was found in PBMCs infected with the heterologous strain M78, in contrast to the cells infected with K49. These data may indicate the beginning of an early immune response in cells infected with a heterologous, but not homologous strain. Transcriptome analysis revealed down-regulation of genes involved in endocytosis and phagocytosis in cells infected with the K49 strain, but not in PBMCs infected with M78. On the contrary, we detected activation of endoplasmic reticulum stress response genes in cells infected with a homologous strain, but not in cells infected with a heterologous strain. This study is the first attempt to determine the differences in the response to ASF infection between homologous and heterologous strains at the cellular level. Our results showed that not only genes of the immune response, but also genes involved in endocytosis and cellular stress response may be important for the formation of cross-protective immunity. This data may be useful for vaccine development or testing of candidate vaccines.

In-Depth Temporal Transcriptome Profiling of an Alphaherpesvirus Using Nanopore Sequencing

In this work, a long-read sequencing (LRS) technique based on the Oxford Nanopore Technology MinION platform was used for quantifying and kinetic characterization of the poly(A) fraction of bovine alphaherpesvirus type 1 (BoHV-1) lytic transcriptome across a 12-h infection period. Amplification-based LRS techniques frequently generate artefactual transcription reads and are biased towards the production of shorter amplicons. To avoid these undesired effects, we applied direct cDNA sequencing, an amplification-free technique. Here, we show that a single promoter can produce multiple transcription start sites whose distribution patterns differ among the viral genes but are similar in the same gene at different timepoints. Our investigations revealed that the circ gene is expressed with immediate–early (IE) kinetics by utilizing a special mechanism based on the use of the promoter of another IE gene (bicp4) for the transcriptional control. Furthermore, we detected an overlap between the initiation of DNA replication and the transcription from the bicp22 gene, which suggests an interaction between the two molecular machineries. This study developed a generally applicable LRS-based method for the time-course characterization of transcriptomes of any organism.

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.

African Swine Fever Virus and host response - transcriptome profiling of the Georgia 2007/1 strain and porcine macrophages

African swine fever virus (ASFV) has a major global economic impact. With a case fatality in domestic pigs approaching 100%, it currently presents the largest threat to animal farming. Although genomic differences between attenuated and highly virulent ASFV strains have been identified, the molecular determinants for virulence at the level of gene expression have remained opaque. Here, we characterize the transcriptome of ASFV genotype II Georgia 2007/1 (GRG) during infection of the physiologically relevant host cells, porcine macrophages. In this study, we applied cap analysis gene expression sequencing (CAGE-seq) to map th0e 5′ ends of viral mRNAs at 5 and 16 h postinfection. A bioinformatics analysis of the sequence context surrounding the transcription start sites (TSSs) enabled us to characterize the global early and late promoter landscape of GRG. We compared transcriptome maps of the GRG isolate and the lab-attenuated BA71V strain that highlighted GRG virulence-specific transcripts belonging to multigene families, including two predicted MGF 100 genes, I7L and I8L. In parallel, we monitored transcriptome changes in the infected host macrophage cells. Of the 9,384 macrophage genes studied, transcripts for 652 host genes were differentially regulated between 5 and 16 h postinfection compared with only 25 between uninfected cells and 5 h postinfection. NF-κB activated genes and lysosome components such as S100 were upregulated, and chemokines such as CCL24, CXCL2, CXCL5, and CXCL8 were downregulated.

IMPORTANCE African swine fever virus (ASFV) causes hemorrhagic fever in domestic pigs, with case fatality rates approaching 100% and no approved vaccines or antivirals. The highly virulent ASFV Georgia 2007/1 strain (GRG) was the first isolated when ASFV spread from Africa to the Caucasus region in 2007, then spreading through Eastern Europe and, more recently, across Asia. We used an RNA-based next-generation sequencing technique called CAGE-seq to map the starts of viral genes across the GRG DNA genome. This has allowed us to investigate which viral genes are expressed during early or late stages of infection and how this is controlled, comparing their expression to the nonvirulent ASFV-BA71V strain to identify key genes that play a role in virulence. In parallel, we investigated how host cells respond to infection, which revealed how the ASFV suppresses components of the host immune response to ultimately win the arms race against its porcine host.

Integrative profiling of Epstein-Barr virus transcriptome using a multiplatform approach

BACKGROUND

Epstein-Barr virus (EBV) is an important human pathogenic gammaherpesvirus with carcinogenic potential. The EBV transcriptome has previously been analyzed using both Illumina-based short read-sequencing and Pacific Biosciences RS II-based long-read sequencing technologies. Since the various sequencing methods have distinct strengths and limitations, the use of multiplatform approaches have proven to be valuable. The aim of this study is to provide a more complete picture on the transcriptomic architecture of EBV.

METHODS

In this work, we apply the Oxford Nanopore Technologies MinION (long-read sequencing) platform for the generation of novel transcriptomic data, and integrate these with other's data generated by another LRS approach, Pacific BioSciences RSII sequencing and Illumina CAGE-Seq and Poly(A)-Seq approaches. Both amplified and non-amplified cDNA sequencings were applied for the generation of sequencing reads, including both oligo-d(T) and random oligonucleotide-primed reverse transcription. EBV transcripts are identified and annotated using the LoRTIA software suite developed in our laboratory.

RESULTS

This study detected novel genes embedded into longer host genes containing 5'-truncated in-frame open reading frames, which potentially encode N-terminally truncated proteins. We also detected a number of novel non-coding RNAs and transcript length isoforms encoded by the same genes but differing in their start and/or end sites. This study also reports the discovery of novel splice isoforms, many of which may represent altered coding potential, and of novel replication-origin-associated transcripts. Additionally, novel mono- and multigenic transcripts were identified. An intricate meshwork of transcriptional overlaps was revealed.

CONCLUSIONS

An integrative approach applying multi-technique sequencing technologies is suitable for reliable identification of complex transcriptomes because each techniques has different advantages and limitations, and the they can be used for the validation of the results obtained by a particular approach.