Deep RNA sequencing reveals complex transcriptional landscape of a bat adenovirus

IRF7 from the black flying fox induces a STAT1-independent ISG response in unstimulated cell lines that protects against diverse RNA viruses

Bats are considered unique in their ability to harbor large numbers of viruses and serve as reservoirs for zoonotic viruses that have the potential to spill over into humans. However, these animals appear relatively resistant to the pathogenic effects of many viruses. Mounting evidence suggests that bats may tolerate viral infections due to unique immune features. These include evolutionary innovations in inflammatory pathways and in the molecules involved in viral sensing, interferon induction, and downstream interferon-induced antiviral effectors. We sought to determine whether interferon-stimulated genes (ISGs) from the black flying fox ( Pteropus alecto ) encoded proteins with unique antiviral activity relative to their human orthologs. Accordingly, we compared the antiviral activity of over 50 ISG human-bat ortholog pairs to identify differences in individual effector functions. We identified IRF7 from Pteropus alecto (Pa.IRF7) as a potent and broad-acting antiviral molecule that provides robust antiviral protection without prior activation. We show that Pa.IRF7 uniquely induces a subset of protective ISGs independent of canonical IFN signaling, which leads to protection from alphaviruses, a flavivirus, a rhabdovirus, and a paramyxovirus. In uninfected cells, Pa.IRF7 partially localizes to the nucleus and can directly bind interferon-sensitive regulatory elements (ISREs). Compared to human IRF7, Pa.IRF7 also has additional serines in its C terminal domain that contribute to antiviral activity and may serve as unique phosphorylation hubs for activation. These properties constitute major differences between bat and human IRF7 that offer additional insight into the potential uniqueness of the black flying fox immune system.

Testing bats in rehabilitation for SARS-CoV-2 before release into the wild

Several studies have suggested SARS-CoV-2 originated from a viral ancestor in bats, but whether transmission occurred directly or via an intermediary host to humans remains unknown. Concerns of spillover of SARS-CoV-2 into wild bat populations are hindering bat rehabilitation and conservation efforts in the United Kingdom and elsewhere. Current protocols state that animals cared for by individuals who have tested positive for SARS-CoV-2 cannot be released into the wild and must be isolated to reduce the risk of transmission to wild populations. Here, we propose a reverse transcription-quantitative polymerase chain reaction (RT-qPCR)-based protocol for detection of SARS-CoV-2 in bats, using fecal sampling. Bats from the United Kingdom were tested following suspected exposure to SARS-CoV-2 and tested negative for the virus. With current UK and international legislation, the identification of SARS-CoV-2 infection in wild animals is becoming increasingly important, and protocols such as the one developed here will help improve understanding and mitigation of SARS-CoV-2 in the future.

Early in vivo transcriptome of Trichoplusia ni ascovirus core genes

Ascoviruses are large double-stranded DNA insect viruses that destroy the nucleus and transform each cell into 20 or more viral vesicles for replication. In the present study we used RNA-sequencing to compare the expression of Trichoplusia ni ascovirus 6a1 (TnAV-6a1) core genes during the first week of infection, with emphasis on the first 48 h, comparing transcript levels in major somatic tissues (epidermis, tracheal matrix and fat body), the sites infected initially, with those of the haemolymph, where viral vesicles circulate and most replication occurs. By 48 h post-infection (p.i.), only 26 genes were expressed in somatic tissues at ≥5 log2 reads per kilobase per million, whereas in the haemolymph 48 genes were expressed at a similar level by the same time. Early and high expression of TnAV caspase-2-like gene occurred in all tissues, implying it is required for replication, but that it is probably not associated with apoptosis induction, which occurs in infections of Spodoptera frugiperda ascovirus 1 a (SfAV-1a), the ascovirus type species. Other highly expressed viral genes at 48 h p.i. in viral vesicles included a dynein-like beta chain and lipid-modifying enzymes, suggesting their importance to vesicle formation and growth as well as virion synthesis. Finally, as occurs in SfAV expression, we found bicistronic and tricistronic mRNA messages produced by TnAV.

The Experience of Greece as a Model to Contain COVID-19 Infection Spread

The severe acute respiratory syndrome corona virus 2 (SARS-CoV-2) emerged in late 2019 and has caused a pandemic known as corona virus disease 2019 (COVID-19), responsible for the death of more than 2 million people worldwide. The outbreak of COVID-19 has posed an unprecedented threat on human lives and public safety. The aim of this review is to describe key aspects of the bio-pathology of the novel disease, and discuss aspects of its spread, as well as targeted protective strategies that can help shape the outcome of the present and future health crises. Greece is used as a model to inhibit SARS-COV-2 spread, since it is one of the countries with the lowest fatality rates among nations of the European Union (E.U.), following two consecutive waves of COVID-19 pandemic. Furthermore, niche research technological approaches and scientific recommendations that emerged during the COVID-19 era are discussed.

Design, validation and evaluation of a SYBR green-based quantitative PCR array for comprehensive analysis of adenovirus type 5 transcriptional patterns

The adenoviral genome encodes coordinately expressed early and late gene transcriptional units that specify a complex collection of extensively spliced overlapping mRNAs. These complexities confound the generation of compatible, validated and optimized qPCR assays that permit comprehensive evaluation of adenoviral transcription. We have developed and evaluated a compilation of qPCR assays that represent the majority of the human adenovirus 5 (hAdV5) genome and allow for absolute and relative quantification of transcriptional activity. A panel of specific adenovirus gene primer pairs was designed through computational modeling to be compatible under a single reaction condition, precisely amplify spliced transcript products within each gene class, and not result in cellular or viral RNA/DNA background amplification. Primer pairs and reaction conditions were optimized to generate a single amplification product that was specific for its target amplicon with minimal intra-assay variability. The specificity of target amplicons was confirmed by dissociation curve analysis, gel electrophoresis and sequencing. In all, thirty-two primer sets representing specific gene products, as well as, pan early and late gene regions were validated under identical amplification conditions, thereby enabling a comprehensive assessment of adenoviral transcription within a single plate array. In order to generate positive control templates and to facilitate absolute quantification of gene expression, all target amplicons were cloned to create gene target-specific standards. These plasmid amplicon controls demonstrated that the SYBR qPCR assays exhibited optimal amplification efficiencies with a high sensitivity of detection to less than 10 copies and a linear amplification across at least eight orders of magnitude. The effectiveness and utility of the comprehensive adenoviral transcriptional array was assessed by investigating the changes in Ad5Wt gene expression at 72 versus 24 h post infection. Predictably, overall gene expression was globally increased at 72 h post infection; however, levels of E2 and Late transcripts exhibited the greatest increased expression, reflecting their necessity at this time point for genomic replication and virion assembly. Taken together, these data demonstrate that the adenoviral qPCR transcriptional array is a modular, scalable, and cost-effective method to comprehensively and accurately assess hAdV5 gene transcription. This array is broadly applicable to facilitate: adenoviral vector development; assessment of cell complementation of knockout viruses; antiviral mechanism of action evaluation; next-generation sequencing data validation.

Virus- and Interferon Alpha-Induced Transcriptomes of Cells from the Microbat Myotis daubentonii

Antiviral interferons (IFN-alpha/beta) are possibly responsible for the high tolerance of bats to zoonotic viruses. Previous studies focused on the IFN system of megabats (suborder Yinpterochiroptera). We present statistically robust RNA sequencing (RNA-seq) data on transcriptomes of cells from the “microbat” Myotis daubentonii (suborder Yangochiroptera) responding at 6 and 24 h to either an IFN-inducing virus or treatment with IFN. Our data reveal genes triggered only by virus, either in both humans and Myotis (CCL4, IFNL3, CH25H), or exclusively in Myotis (STEAP4). Myotis cells also express a series of conserved IFN-stimulated genes (ISGs) and an unusually high paralog number of the antiviral ISG BST2 (tetherin) but lack several ISGs that were described for megabats (EMC2, FILIP1, IL17RC, OTOGL, SLC24A1). Also, in contrast to megabats, we detected neither different IFN-alpha subtypes nor an unusually high baseline expression of IFNs. Thus, Yangochiroptera microbats, represented by Myotis, may possess an IFN system with distinctive features.

•

Virus- and IFN-responsive transcriptomes of the microbat Myotis daubentonii

•

CCL4, IFNL3, CH25H, STEAP4 are IFNB-like genes triggered by virus only

•

Microbats encode more paralogs of BST2 (tetherin) than any other mammal

•

Clear differences between the IFN systems of microbats and megabats

A metaanalysis of bat phylogenetics and positive selection based on genomes and transcriptomes from 18 species

This work represents a large, order-wide evolutionary analysis of the order Chiroptera (bats). Our pipeline for assembling sequence data and curating orthologous multiple sequence alignments includes methods for improving results when combining genomic and transcriptomic data sources. The resulting phylogenetic tree divides the order Chiroptera into Yinpterochiroptera and Yangochiroptera, in disagreement with the previous division into Megachiroptera and Microchiroptera and in agreement with some other recent molecular studies, and also provides evidence for other contested branch placements. We also performed a genome-wide analysis of positive selection and found 181 genes with signatures of positive selection. Enrichment analysis shows these positively selected genes to be primarily related to immune responses but also, surprisingly, collagen formation.

Historically, the evolution of bats has been analyzed using a small number of genetic loci for many species or many genetic loci for a few species. Here we present a phylogeny of 18 bat species, each of which is represented in 1,107 orthologous gene alignments used to build the tree. We generated a transcriptome sequence of Hypsignathus monstrosus, the African hammer-headed bat, and additional transcriptome sequence for Rousettus aegyptiacus, the Egyptian fruit bat. We then combined these data with existing genomic and transcriptomic data from 16 other bat species. In the analysis of such datasets, there is no clear consensus on the most reliable computational methods for the curation of quality multiple sequence alignments since these public datasets represent multiple investigators and methods, including different source materials (chromosomal DNA or expressed RNA). Here we lay out a systematic analysis of parameters and produce an advanced pipeline for curating orthologous gene alignments from combined transcriptomic and genomic data, including a software package: the Mismatching Isoform eXon Remover (MIXR). Using this method, we created alignments of 11,677 bat genes, 1,107 of which contain orthologs from all 18 species. Using the orthologous gene alignments created, we assessed bat phylogeny and also performed a holistic analysis of positive selection acting in bat genomes. We found that 181 genes have been subject to positive natural selection. This list is dominated by genes involved in immune responses and genes involved in the production of collagens.

Adenoviromics: Mining the Human Adenovirus Species D Genome

Human adenovirus (HAdV) infections cause disease world-wide. Whole genome sequencing has now distinguished 90 distinct genotypes in 7 species (A-G). Over half of these 90 HAdVs fall within species D, with essentially all of the HAdV-D whole genome sequences generated in the last decade. Herein, we describe recent new findings made possible by mining of this expanded genome database, and propose future directions to elucidate new functional elements and new functions for previously known viral components.

The IFN Response in Bats Displays Distinctive IFN-Stimulated Gene Expression Kinetics with Atypical RNASEL Induction

Bats host a large number of zoonotic viruses, including several viruses that are highly pathogenic to other mammals. The mechanisms underlying this rich viral diversity are unknown, but they may be linked to unique immunological features that allow bats to act as asymptomatic viral reservoirs. Vertebrates respond to viral infection by inducing IFNs, which trigger antiviral defenses through IFN-stimulated gene (ISG) expression. Although the IFN system of several bats is characterized at the genomic level, less is known about bat IFN-mediated transcriptional responses. In this article, we show that IFN signaling in bat cells from the black flying fox (Pteropus alecto) consists of conserved and unique ISG expression profiles. In IFN-stimulated cells, bat ISGs comprise two unique temporal subclusters with similar early induction kinetics but distinct late-phase declines. In contrast, human ISGs lack this decline phase and remained elevated for longer periods. Notably, in unstimulated cells, bat ISGs were expressed more highly than their human counterparts. We also found that the antiviral effector 2-5A-dependent endoribonuclease, which is not an ISG in humans, is highly IFN inducible in black flying fox cells and contributes to cell-intrinsic control of viral infection. These studies reveal distinctive innate immune features that may underlie a unique virus-host relationship in bats.

Evolution and Cryo-electron Microscopy Capsid Structure of a North American Bat Adenovirus and Its Relationship to Other Mastadenoviruses

Since the first description of adenoviruses in bats in 2006, a number of micro- and megabat species in Europe, Africa, and Asia have been shown to carry a wide diversity of adenoviruses. Here, we report on the evolutionary, biological, and structural characterization of a novel bat adenovirus (BtAdV) recovered from a Rafinesque's big-eared bat (Corynorhinus rafinesquii) in Kentucky, USA, which is the first adenovirus isolated from North American bats. This virus (BtAdV 250-A) exhibits a close phylogenetic relationship with Canine mastadenovirus A (CAdV A), as previously observed with other BtAdVs. To further investigate the relationships between BtAdVs and CAdVs, we conducted mass spectrometric analysis and single-particle cryo-electron microscopy reconstructions of the BtAdV 250-A capsid and also analyzed the in vitro host ranges of both viruses. Our results demonstrate that BtAdV 250-A represents a new mastadenovirus species that, in contrast to CAdV, has a unique capsid morphology that contains more prominent extensions of protein IX and can replicate efficiently in a phylogenetically diverse range of species. These findings, in addition to the recognition that both the genetic diversity of BtAdVs and the number of different bat species from disparate geographic regions infected with BtAdVs appears to be extensive, tentatively suggest that bats may have served as a potential reservoir for the cross-species transfer of adenoviruses to other hosts, as theorized for CAdV.

IMPORTANCE

Although many adenoviruses are host specific and likely codiverged with their hosts over millions of years, other adenoviruses appear to have emerged through successful cross-species transmission events on more recent time scales. The wide geographic distribution and genetic diversity of adenoviruses in bats and their close phylogenetic relationship to Canine mastadenovirus A (CAdV A) has raised important questions about how CAdV A, and possibly other mammalian adenoviruses, may have emerged. Although most adenoviruses tend to cause limited disease in their natural hosts, CAdV A is unusual in that it may cause high morbidity and sometimes fatal infections in immunocompetent hosts and is thus an important pathogen of carnivores. Here, we performed a comparative evolutionary and structural study of representative bat and canine adenoviruses to better understand the relationship between these two viral groups.

First discovery of novel 3-hydroxy-quinazoline-2,4(1H,3H)-diones as specific anti-vaccinia and adenovirus agents via 'privileged scaffold' refining approach

A series of 1,2,3-triazolyl 3-hydroxy-quinazoline-2,4(1H,3H)-diones was constructed utilizing Cu(I)-catalyzed azide-alkyne 1,3-dipolar cycloaddition (CuAAC) method. The biological significance of the novel synthesized quinazolines was highlighted by evaluating them in vitro for antiviral activity, wherein several compounds exhibited excellent activity specifically against vaccinia and adenovirus. Especially, 24b11 displayed the most potent inhibitory activity against vaccinia with an EC50 value of 1.7μM, which was 15 fold than that of the reference drug Cidofovir (EC50=25μM). 24b13 was the most potent compound against adenovirus-2 with an EC50 value of 6.2μM, which proved lower than all the reference drugs. Preliminary structure-activity relationships were also discussed. To the best of our knowledge, no data are present in the literature on antiviral activity of 3-hydroxy-quinazoline-2,4(1H,3H)-diones against DNA-viruses. Thus, these findings warrant further investigations (library expansion and compound refinement) on this novel class of antiviral agents.

Novel bat adenoviruses with an extremely large E3 gene

Bats carry diverse RNA viruses, some of which are responsible for human diseases. Compared to bat-borne RNA viruses, relatively little information is known regarding bat-borne DNA viruses. In this study, we isolated and characterized three novel bat adenoviruses (BtAdV WIV9-11) from Rhinolophus sinicus. Their genomes, which are highly similar to each other but distinct from those of previously sequenced adenoviruses (AdVs), are 37 545, 37 566 and 38 073 bp in size, respectively. An unusually large E3 gene was identified in their genomes. Phylogenetic and taxonomic analyses suggested that these isolates represent a distinct species of the genus Mastadenovirus. Cell susceptibility assays revealed a broad cell tropism for these isolates, indicating that they have a potentially wide host range. Our results expand the understanding of genetic diversity of bat AdVs.

Comparison of the Life Cycles of Genetically Distant Species C and Species D Human Adenoviruses Ad6 and Ad26 in Human Cells

Our understanding of adenovirus (Ad) biology is largely extrapolated from human species C Ad5. Most humans are immune to Ad5, so lower-seroprevalence viruses like human Ad6 and Ad26 are being tested as therapeutic vectors. Ad6 and Ad26 differ at the DNA level by 34%. To better understand how this might impact their biology, we examined the life cycle of the two viruses in human lung cells in vitro. Both viruses infected A549 cells with similar efficiencies, executed DNA replication with identical kinetics within 12 h, and began killing cells within 72 h. While Ad6-infected cells remained adherent until death, Ad26-infected cells detached within 12 h of infection but remained viable. Next-generation sequencing (NGS) of mRNA from infected cells demonstrated that viral transcripts constituted 1% of cellular mRNAs within 6 h and 8 to 16% within 12 h. Quantitative PCR and NGS revealed the activation of key early genes at 6 h and transition to late gene activation by 12 h by both viruses. There were marked differences in the balance of E1A and E1B activation by the two viruses and in the expression of E3 immune evasion mRNAs. Ad6 was markedly more effective at suppressing major histocompatibility complex class I (MHC I) display on the cell surface and in evading TRAIL-mediated apoptosis than was Ad26. These data demonstrate shared as well as divergent life cycles in these genetically distant human adenoviruses. An understanding of these differences expands the knowledge of alternative Ad species and may inform the selection of related Ads for therapeutic development.

IMPORTANCE

A burgeoning number of adenoviruses (Ads) are being harnessed as therapeutics, yet the biology of these viruses is generally extrapolated from Ad2 and Ad5. Here, we are the first to compare the transcriptional programs of two genetically distant Ads by mRNA next-generation sequencing (NGS). Species C Ad6 and Ad26 are being pursued as lower-seroprevalence Ad vectors but differ at the DNA level by 34%. Head-to-head comparison in human lung cells by NGS revealed that the two viruses generally conform to our general understanding of the Ad transcriptional program. However, fine mapping revealed subtle and strong differences in how these two viruses execute these programs, including differences in the balance of E1A and E1B mRNAs and in E3 immune evasion genes. This suggests that not all adenoviruses behave like Ad2 and Ad5 and that they may have unique strategies to infect cells and evade the immune system.