Comparative analysis of European bat lyssavirus 1 pathogenicity in the mouse model

Transcriptomic responses of bat cells to European bat lyssavirus 1 infection under conditions simulating euthermia and hibernation

BACKGROUND

Coevolution between pathogens and their hosts decreases host morbidity and mortality. Bats host and can tolerate viruses which can be lethal to other vertebrate orders, including humans. Bat adaptations to infection include localized immune response, early pathogen sensing, high interferon expression without pathogen stimulation, and regulated inflammatory response. The immune reaction is costly, and bats suppress high-cost metabolism during torpor. In the temperate zone, bats hibernate in winter, utilizing a specific behavioural adaptation to survive detrimental environmental conditions and lack of energy resources. Hibernation torpor involves major physiological changes that pose an additional challenge to bat-pathogen coexistence. Here, we compared bat cellular reaction to viral challenge under conditions simulating hibernation, evaluating the changes between torpor and euthermia.

RESULTS

We infected the olfactory nerve-derived cell culture of Myotis myotis with an endemic bat pathogen, European bat lyssavirus 1 (EBLV-1). After infection, the bat cells were cultivated at two different temperatures, 37 °C and 5 °C, to examine the cell response during conditions simulating euthermia and torpor, respectively. The mRNA isolated from the cells was sequenced and analysed for differential gene expression attributable to the temperature and/or infection treatment. In conditions simulating euthermia, infected bat cells produce an excess signalling by multitude of pathways involved in apoptosis and immune regulation influencing proliferation of regulatory cell types which can, in synergy with other produced cytokines, contribute to viral tolerance. We found no up- or down-regulated genes expressed in infected cells cultivated at conditions simulating torpor compared to non-infected cells cultivated under the same conditions. When studying the reaction of uninfected cells to the temperature treatment, bat cells show an increased production of heat shock proteins (HSPs) with chaperone activity, improving the bat's ability to repair molecular structures damaged due to the stress related to the temperature change.

CONCLUSIONS

The lack of bat cell reaction to infection in conditions simulating hibernation may contribute to the virus tolerance or persistence in bats. Together with the cell damage repair mechanisms induced in response to hibernation, the immune regulation may promote bats' ability to act as reservoirs of zoonotic viruses such as lyssaviruses.

Oral Rabies Vaccine Strain SPBN GASGAS: Genetic Stability after Serial In Vitro and In Vivo Passaging

Oral vaccination of wildlife has shown to be a very effective management tool in rabies control. Evaluation of the genetic stability of vaccine viruses before distributing vaccine baits in the environment is essential because all available oral rabies vaccines, including the genetically engineered rabies virus vaccine strain SPBN GASGAS (Rabitec), are based on replication-competent viruses. To evaluate the genetic stability of this vaccine strain, five serial passages of the Master Seed Virus (MSV) in the production cell line BHK21 Cl13 were performed. Furthermore, to test possible reversion to virulence, a back-passage study in suckling mouse brain (SMB) was performed. Subsequently, the pooled 5th SMB passage was inoculated intracerebrally (i.c.) in adult and suckling mice. The full genome sequences of the isolated 5th passage, in vivo and in vitro, were compared at both the consensus and the quasispecies level with the MSV. Additionally, the full genome sequence of the 6th SMB passage from the individual animals was determined and compared. Full-length integration of the double glycoprotein and modified base substitutions at amino acid position 194 and 333 of the glycoprotein could be verified in all 5th and 6th passage samples. Overall, 11 single nucleotide polymorphisms (SNPs) were detected in the 5th pooled SMB passage, 4 with frequency between 10 and 20%, and 7 with between 2.5 and 10%. SNPs that resulted in amino acid exchange were found in genes: N (one SNP), G (four SNPs), and L (three SNPs). However, none of these SNPs were associated with reversion to virulence since all adult mice inoculated i.c. with this material survived. In the individual samples of the 6th SMB passage 24 additional SNPs (>2.5%) were found, of which only 1 SNP (L-gene, position 6969) had a prevalence of >50% in 3 of 17 samples. The obtained results confirmed the stable expression of genetic modifications and the genetic stability of the consensus strain after serial in vivo and in vitro passaging.

Comparative pathogenesis of different phylogroup I bat lyssaviruses in a standardized mouse model

A plethora of bat-associated lyssaviruses potentially capable of causing the fatal disease rabies are known today. Transmitted via infectious saliva, occasionally-reported spillover infections from bats to other mammals demonstrate the permeability of the species-barrier and highlight the zoonotic potential of bat-related lyssaviruses. However, it is still unknown whether and, if so, to what extent, viruses from different lyssavirus species vary in their pathogenic potential. In order to characterize and systematically compare a broader group of lyssavirus isolates for their viral replication kinetics, pathogenicity, and virus release through saliva-associated virus shedding, we used a mouse infection model comprising a low (10² TCID₅₀) and a high (10⁵ TCID₅₀) inoculation dose as well as three different inoculation routes (intramuscular, intranasal, intracranial). Clinical signs, incubation periods, and survival were investigated. Based on the latter two parameters, a novel pathogenicity matrix was introduced to classify lyssavirus isolates. Using a total of 13 isolates from ten different virus species, this pathogenicity index varied within and between virus species. Interestingly, Irkut virus (IRKV) and Bokeloh bat lyssavirus (BBLV) obtained higher pathogenicity scores (1.14 for IRKV and 1.06 for BBLV) compared to rabies virus (RABV) isolates ranging between 0.19 and 0.85. Also, clinical signs differed significantly between RABV and other bat lyssaviruses. Altogether, our findings suggest a high diversity among lyssavirus isolates concerning survival, incubation period, and clinical signs. Virus shedding significantly differed between RABVs and other lyssaviruses. Our results demonstrated that active shedding of infectious virus was exclusively associated with two RABV isolates (92% for RABV-DogA and 67% for RABV-Insectbat), thus providing a potential explanation as to why sustained spillovers are solely attributed to RABVs. Interestingly, 3D imaging of a selected panel of brain samples from bat-associated lyssaviruses demonstrated a significantly increased percentage of infected astrocytes in mice inoculated with IRKV (10.03%; SD±7.39) compared to RABV-Vampbat (2.23%; SD±2.4), and BBLV (0.78%; SD±1.51), while only individual infected cells were identified in mice infected with Duvenhage virus (DUVV). These results corroborate previous studies on RABV that suggest a role of astrocyte infection in the pathogenicity of lyssaviruses.

Globally, there are at present 17 different officially recognized lyssavirus species posing a potential threat for human and animal health. Bats have been identified as carriers for the vast majority of those zoonotic viruses, which cause the fatal disease rabies and are transmitted through infectious saliva. The occurrence of sporadic spillover events where lyssaviruses are spread from bats to other mammalian species highlights the importance of studying pathogenicity and virus shedding in regard to a potentially sustained onward cross-species transmission. Therefore, as part of this study, we compared 13 different isolates from ten lyssavirus species in a standardized mouse infection model, focusing on clinical signs, incubation periods, and survival. Based on the latter two, a novel pathogenicity index to classify different lyssavirus species was established. This pathogenicity index varied within and between different lyssavirus species and revealed a higher ranking of other bat-related lyssaviruses in comparison to the tested Rabies virus (RABV) isolates. Altogether, our results demonstrate a high diversity among the investigated isolates concerning pathogenicity and clinical picture. Furthermore, we comparatively analyzed virus shedding via saliva and while there was no indication towards a reduced pathogenicity of bat-associated lyssaviruses as opposed to RABV, shedding was increased in RABV isolates. Additionally, we investigated neuronal cell tropism and revealed that bat lyssaviruses are not only capable of infecting neurons but also astrocytes.

First case of lethal encephalitis in Western Europe due to European bat lyssavirus type 1

BACKGROUND

Inaccurate diagnosis of encephalitis is a major issue as immunosuppressive treatments can be deleterious in case of viral infection. The European bat lyssavirus type 1, a virus related to rabies virus, is endemic in European bats. No human case has yet been reported in Western Europe. A 59 year-old patient without specific past medical history died from encephalitis. A colony of bats lived in an outbuilding of his house. No diagnosis was made using standard procedures.

METHODS

We used a Next Generation Sequencing (NGS) based transcriptomic protocol to search for pathogens in autopsy samples (meninges and brain frontal lobe). Results were confirmed by PCR and by antibody testing in serum. Immunochemistry was used to characterize inflammatory cells and viral antigens in brain lesions. Cells and mice were inoculated with brain extracts for virus isolation.

RESULTS

The patient's brain lesions were severe and diffuse in white and gray matter. Perivascular inflammatory infiltrates were abundant and rich in plasma cells. NGS identified European bat lyssavirus type 1a in brain, which was confirmed by PCR. A high titer of neutralizing antibodies was found in serum. No viral antigen was detected and the virus could not be isolated by cell culture or by mouse inoculation.

CONCLUSIONS

The patient died from European bat lyssavirus type 1a infection. NGS was key to identifying this unexpected viral etiology in an epidemiological context that did not suggest rabies. People exposed to bats should be strongly advised to be vaccinated with rabies vaccines, which are effective against EBLV-1.

Full-Genome Sequences and Phylogenetic Analysis of Archived Danish European Bat Lyssavirus 1 (EBLV-1) Emphasize a Higher Genetic Resolution and Spatial Segregation for Sublineage 1a

European bat lyssavirus type 1 (EBLV-1) is the causative agent for almost all reported rabies cases found in European bats. In recent years, increasing numbers of available EBLV-1 full genomes and their phylogenetic analyses helped to further elucidate the distribution and genetic characteristics of EBLV-1 and its two subtypes, namely EBLV-1a and EBLV-1b. Nonetheless, the absence of full-genome sequences from regions with known detections of EBLV-1 still limit the understanding of the phylogeographic relations between viruses from different European regions. In this study, a set of 21 archived Danish EBLV-1 samples from the years 1985 to 2009 was processed for the acquisition of full-genome sequences using a high-throughput sequencing approach. Subsequent phylogenetic analysis encompassing all available EBLV-1 full genomes from databases revealed the Danish sequences belong to the EBLV-1a subtype and further highlighted the distinct, close phylogenetic relationship of Danish, Dutch and German isolates in this region. In addition, the formation of five putative groups nearly exclusively formed by Danish isolates and the overall increased resolution of the EBLV-1a branch indicate a higher genetic diversity and spatial segregation for this sublineage than was previously known. These results emphasize the importance of phylogenetic analyses of full-genome sequences of lyssaviruses for genetic geography.

Bat-Borne Coronaviruses in Jordan and Saudi Arabia: A Threat to Public Health?

Emerging infectious diseases are of great concern to public health, as highlighted by the ongoing coronavirus disease 2019 (COVID-19) pandemic. Such diseases are of particular danger during mass gathering and mass influx events, as large crowds of people in close proximity to each other creates optimal opportunities for disease transmission. The Hashemite Kingdom of Jordan and the Kingdom of Saudi Arabia are two countries that have witnessed mass gatherings due to the arrival of Syrian refugees and the annual Hajj season. The mass migration of people not only brings exotic diseases to these regions but also brings new diseases back to their own countries, e.g., the outbreak of MERS in South Korea. Many emerging pathogens originate in bats, and more than 30 bat species have been identified in these two countries. Some of those bat species are known to carry viruses that cause deadly diseases in other parts of the world, such as the rabies virus and coronaviruses. However, little is known about bats and the pathogens they carry in Jordan and Saudi Arabia. Here, the importance of enhanced surveillance of bat-borne infections in Jordan and Saudi Arabia is emphasized, promoting the awareness of bat-borne diseases among the general public and building up infrastructure and capability to fill the gaps in public health preparedness to prevent future pandemics.

Astrocyte Infection during Rabies Encephalitis Depends on the Virus Strain and Infection Route as Demonstrated by Novel Quantitative 3D Analysis of Cell Tropism

Although conventional immunohistochemistry for neurotropic rabies virus (RABV) usually shows high preference for neurons, non-neuronal cells are also potential targets, and abortive astrocyte infection is considered a main trigger of innate immunity in the CNS. While in vitro studies indicated differences between field and less virulent lab-adapted RABVs, a systematic, quantitative comparison of astrocyte tropism in vivo is lacking. Here, solvent-based tissue clearing was used to measure RABV cell tropism in infected brains. Immunofluorescence analysis of 1 mm-thick tissue slices enabled 3D-segmentation and quantification of astrocyte and neuron infection frequencies. Comparison of three highly virulent field virus clones from fox, dog, and raccoon with three lab-adapted strains revealed remarkable differences in the ability to infect astrocytes in vivo. While all viruses and infection routes led to neuron infection frequencies between 7–19%, striking differences appeared for astrocytes. Whereas astrocyte infection by field viruses was detected independent of the inoculation route (8–27%), only one lab-adapted strain infected astrocytes route-dependently [0% after intramuscular (i.m.) and 13% after intracerebral (i.c.) inoculation]. Two lab-adapted vaccine viruses lacked astrocyte infection altogether (0%, i.c. and i.m.). This suggests a model in which the ability to establish productive astrocyte infection in vivo functionally distinguishes field and attenuated lab RABV strains.

In-depth genome analyses of viruses from vaccine-derived rabies cases and corresponding live-attenuated oral rabies vaccines

Live-attenuated rabies virus strains such as those derived from the field isolate Street Alabama Dufferin (SAD) have been used extensively and very effectively as oral rabies vaccines for the control of fox rabies in both Europe and Canada. Although these vaccines are safe, some cases of vaccine-derived rabies have been detected during rabies surveillance accompanying these campaigns. In recent analysis it was shown that some commercial SAD vaccines consist of diverse viral populations, rather than clonal genotypes. For cases of vaccine-derived rabies, only consensus sequence data have been available to date and information concerning their population diversity was thus lacking. In our study, we used high-throughput sequencing to analyze 11 cases of vaccine-derived rabies, and compared their viral population diversity to the related oral rabies vaccines using pairwise Manhattan distances. This extensive deep sequencing analysis of vaccine-derived rabies cases observed during oral vaccination programs provided deeper insights into the effect of accidental in vivo replication of genetically diverse vaccine strains in the central nervous system of target and non-target species under field conditions. The viral population in vaccine-derived cases appeared to be clonal in contrast to their parental vaccines. The change from a state of high population diversity present in the vaccine batches to a clonal genotype in the affected animal may indicate the presence of a strong bottleneck during infection. In conclusion, it is very likely that these few cases are the consequence of host factors and not the result of the selection of a more virulent genotype. Furthermore, this type of vaccine-derived rabies leads to the selection of clonal genotypes and the selected variants were genetically very similar to potent SAD vaccines that have undergone a history of in vitro selection.

Defining objective clusters for rabies virus sequences using affinity propagation clustering

Rabies is caused by lyssaviruses, and is one of the oldest known zoonoses. In recent years, more than 21,000 nucleotide sequences of rabies viruses (RABV), from the prototype species rabies lyssavirus, have been deposited in public databases. Subsequent phylogenetic analyses in combination with metadata suggest geographic distributions of RABV. However, these analyses somewhat experience technical difficulties in defining verifiable criteria for cluster allocations in phylogenetic trees inviting for a more rational approach. Therefore, we applied a relatively new mathematical clustering algorythm named ‘affinity propagation clustering’ (AP) to propose a standardized sub-species classification utilizing full-genome RABV sequences. Because AP has the advantage that it is computationally fast and works for any meaningful measure of similarity between data samples, it has previously been applied successfully in bioinformatics, for analysis of microarray and gene expression data, however, cluster analysis of sequences is still in its infancy. Existing (516) and original (46) full genome RABV sequences were used to demonstrate the application of AP for RABV clustering. On a global scale, AP proposed four clusters, i.e. New World cluster, Arctic/Arctic-like, Cosmopolitan, and Asian as previously assigned by phylogenetic studies. By combining AP with established phylogenetic analyses, it is possible to resolve phylogenetic relationships between verifiably determined clusters and sequences. This workflow will be useful in confirming cluster distributions in a uniform transparent manner, not only for RABV, but also for other comparative sequence analyses.

Rabies is one of the oldest known zoonoses, caused by lyssaviruses. In recent years, more than 21,000 nucleotide sequences for rabies viruses (RABV) have been deposited in public databases. In this study, a novel mathematical approach called affinity propagation (AP) clustering, a highly powerful tool, to verifiably divide full genome RABV sequences into genetic clusters, was used. A panel of existing and novel RABV full genome sequences was used to demonstrate the application of AP for RABV clustering. Using a combination of AP with established phylogenetic analyses is useful in resolving phylogenetic relationships between more objectively determined clusters and sequences. This workflow will help to substantiate a transparent cluster distribution, not only for RABV, but also for other comparative sequence analyses.

Bat Biology, Genomes, and the Bat1K Project: To Generate Chromosome-Level Genomes for All Living Bat Species

Bats are unique among mammals, possessing some of the rarest mammalian adaptations, including true self-powered flight, laryngeal echolocation, exceptional longevity, unique immunity, contracted genomes, and vocal learning. They provide key ecosystem services, pollinating tropical plants, dispersing seeds, and controlling insect pest populations, thus driving healthy ecosystems. They account for more than 20% of all living mammalian diversity, and their crown-group evolutionary history dates back to the Eocene. Despite their great numbers and diversity, many species are threatened and endangered. Here we announce Bat1K, an initiative to sequence the genomes of all living bat species (n∼1,300) to chromosome-level assembly. The Bat1K genome consortium unites bat biologists (>148 members as of writing), computational scientists, conservation organizations, genome technologists, and any interested individuals committed to a better understanding of the genetic and evolutionary mechanisms that underlie the unique adaptations of bats. Our aim is to catalog the unique genetic diversity present in all living bats to better understand the molecular basis of their unique adaptations; uncover their evolutionary history; link genotype with phenotype; and ultimately better understand, promote, and conserve bats. Here we review the unique adaptations of bats and highlight how chromosome-level genome assemblies can uncover the molecular basis of these traits. We present a novel sequencing and assembly strategy and review the striking societal and scientific benefits that will result from the Bat1K initiative.

Host Genetic Variation Does Not Determine Spatio-Temporal Patterns of European Bat 1 Lyssavirus

The majority of bat rabies cases in Europe are attributed to European bat 1 lyssavirus (EBLV-1), circulating mainly in serotine bats (Eptesicus serotinus). Two subtypes have been defined (EBLV-1a and EBLV-1b), each associated with a different geographical distribution. In this study, we undertake a comprehensive sequence analysis based on 80 newly obtained EBLV-1 nearly complete genome sequences from nine European countries over a 45-year period to infer selection pressures, rates of nucleotide substitution, and evolutionary time scale of these two subtypes in Europe. Our results suggest that the current lineage of EBLV-1 arose in Europe ∼600 years ago and the virus has evolved at an estimated average substitution rate of ∼4.19×10-5 subs/site/year, which is among the lowest recorded for RNA viruses. In parallel, we investigate the genetic structure of French serotine bats at both the nuclear and mitochondrial level and find that they constitute a single genetic cluster. Furthermore, Mantel tests based on interindividual distances reveal the absence of correlation between genetic distances estimated between viruses and between host individuals. Taken together, this indicates that the genetic diversity observed in our E. serotinus samples does not account for EBLV-1a and -1b segregation and dispersal in Europe.