Hepatobiliary organoids derived from leporids support the replication of hepatotropic lagoviruses.

The genus Lagovirus of the family Caliciviridae contains some of the most virulent vertebrate viruses known. Lagoviruses infect leporids, such as rabbits, hares and cottontails. Highly pathogenic viruses such as Rabbit haemorrhagic disease virus 1 (RHDV1) cause a fulminant hepatitis that typically leads to disseminated intravascular coagulation within 24-72 h of infection, killing over 95 % of susceptible animals. Research into the pathophysiological mechanisms that are responsible for this extreme phenotype has been hampered by the lack of a reliable culture system. Here, we report on a new ex vivo model for the cultivation of lagoviruses in cells derived from the European rabbit (Oryctolagus cuniculus) and European brown hare (Lepus europaeus). We show that three different lagoviruses, RHDV1, RHDV2 and RHDVa-K5, replicate in monolayer cultures derived from rabbit hepatobiliary organoids, but not in monolayer cultures derived from cat (Felis catus) or mouse (Mus musculus) organoids. Virus multiplication was demonstrated by (i) an increase in viral RNA levels, (ii) the accumulation of dsRNA viral replication intermediates and (iii) the expression of viral structural and non-structural proteins. The establishment of an organoid culture system for lagoviruses will facilitate studies with considerable implications for the conservation of endangered leporid species in Europe and North America, and the biocontrol of overabundant rabbit populations in Australia and New Zealand.

What a Difference a Gene Makes: Identification of Virulence Factors of Cowpox Virus.

Cowpox virus (CPXV) is a zoonotic orthopoxvirus (OPV) that causes spillover infections from its animal hosts to humans. In 2009, several human CPXV cases occurred through transmission from pet rats. An isolate from a diseased rat, RatPox09, exhibited significantly increased virulence in Wistar rats and caused high mortality compared to that caused by the mildly virulent laboratory strain Brighton Red (BR). The RatPox09 genome encodes four genes which are absent in the BR genome. We hypothesized that their gene products could be major factors influencing the high virulence of RatPox09. To address this hypothesis, we employed several BR-RatPox09 chimeric viruses. Using Red-mediated mutagenesis, we generated BR-based knock-in mutants with single or multiple insertions of the respective RatPox09 genes. High-throughput sequencing was used to verify the genomic integrity of all recombinant viruses, and transcriptomic analyses confirmed that the expression profiles of the genes that were adjacent to the modified ones were unaltered. While the in vitro growth kinetics were comparable to those of BR and RatPox09, we discovered that a knock-in BR mutant containing the four RatPox09-specific genes was as virulent as the RatPox09 isolate, causing death in over 75% of infected Wistar rats. Unexpectedly, the insertion of gCPXV0030 (g7tGP) alone into the BR genome resulted in significantly higher clinical scores and lower survival rates matching the rate for rats infected with RatPox09. The insertion of gCPXV0284, encoding the BTB (broad-complex, tramtrack, and bric-à-brac) domain protein D7L, also increased the virulence of BR, while the other two open reading frames failed to rescue virulence independently. In summary, our results confirmed our hypothesis that a relatively small set of four genes can contribute significantly to CPXV virulence in the natural rat animal model.IMPORTANCE With the cessation of vaccination against smallpox and its assumed cross-protectivity against other OPV infections, waning immunity could open up new niches for related poxviruses. Therefore, the identification of virulence mechanisms in CPXV is of general interest. Here, we aimed to identify virulence markers in an experimental rodent CPXV infection model using bacterial artificial chromosome (BAC)-based virus recombineering. We focused our work on the recent zoonotic CPXV isolate RatPox09, which is highly pathogenic in Wistar rats, unlike the avirulent BR reference strain. In several animal studies, we were able to identify a novel set of CPXV virulence genes. Two of the identified virulence genes, encoding a putative BTB/POZ protein (CPXVD7L) and a B22R-family protein (CPXV7tGP), respectively, have not yet been described to be involved in CPXV virulence. Our results also show that single genes can significantly affect virulence, thus facilitating adaptation to other hosts.

Investigations into the presence of nidoviruses in pythons.

BACKGROUND: Pneumonia and stomatitis represent severe and often fatal diseases in different captive snakes. Apart from bacterial infections, paramyxo-, adeno-, reo- and arenaviruses cause these diseases. In 2014, new viruses emerged as the cause of pneumonia in pythons. In a few publications, nidoviruses have been reported in association with pneumonia in ball pythons and a tiger python. The viruses were found using new sequencing methods from the organ tissue of dead animals. METHODS: Severe pneumonia and stomatitis resulted in a high mortality rate in a captive breeding collection of green tree pythons. Unbiased deep sequencing lead to the detection of nidoviral sequences. A developed RT-qPCR was used to confirm the metagenome results and to determine the importance of this virus. A total of 1554 different boid snakes, including animals suffering from respiratory diseases as well as healthy controls, were screened for nidoviruses. Furthermore, in addition to two full-length sequences, partial sequences were generated from different snake species. RESULTS: The assembled full-length snake nidovirus genomes share only an overall genome sequence identity of less than 66.9% to other published snake nidoviruses and new partial sequences vary between 99.89 and 79.4%. Highest viral loads were detected in lung samples. The snake nidovirus was not only present in diseased animals, but also in snakes showing no typical clinical signs. CONCLUSIONS: Our findings further highlight the possible importance of snake nidoviruses in respiratory diseases and proof multiple circulating strains with varying disease potential. Nidovirus detection in clinical healthy individuals might represent testing during the incubation period or reconvalescence. Our investigations show new aspects of nidovirus infections in pythons. Nidoviruses should be included in routine diagnostic workup of diseased reptiles.

A Variegated Squirrel Bornavirus Associated with Fatal Human Encephalitis.

Between 2011 and 2013, three breeders of variegated squirrels (Sciurus variegatoides) had encephalitis with similar clinical signs and died 2 to 4 months after onset of the clinical symptoms. With the use of a metagenomic approach that incorporated next-generation sequencing and real-time reverse-transcriptase quantitative polymerase chain reaction (RT-qPCR), the presence of a previously unknown bornavirus was detected in a contact squirrel and in brain samples from the three patients. Phylogenetic analyses showed that this virus, tentatively named variegated squirrel 1 bornavirus (VSBV-1), forms a lineage separate from that of the known bornavirus species. (Funded by the Federal Ministry of Food and Agriculture [Germany] and others.).

Variegated Squirrel Bornavirus 1 in Squirrels, Germany and the Netherlands.

We screened squirrels in Germany and the Netherlands for the novel zoonotic variegated squirrel bornavirus 1 (VSBV-1). The detection of VSBV-1 in 11 squirrels indicates a considerable risk for transmission to humans handling those animals. Therefore, squirrels in contact with humans should routinely be tested for VSBV-1.

Quasispecies composition and diversity do not reveal any predictors for chronic classical swine fever virus infection.

Classical swine fever (CSF) can run acute, chronic, and prenatal courses in both domestic pigs and wild boar. Although chronic infections are rare events, their epidemiological impact is very high due to the long-term shedding of virus. So far, little is known about the factors that influence disease course and outcome from either the host or virus's perspective. To elucidate the viral determinants, we analyzed the role of the viral populations for the development of chronic CSF virus (CSFV) infections. Three different animal trials that had led to both chronic and acute infections were chosen for a detailed analysis by deep sequencing. The three inocula represented sub-genogroups 2.1 and 2.3, and two viruses were wild-type CSFV, one derived from an infectious cDNA clone. These viruses and samples derived from acutely and chronically infected animals were subjected to next-generation sequencing. Subsequently, the derived full-length genomes were compared at both the consensus and the quasispecies level. At consensus level, no differences were observed between the parental viruses and the viruses obtained from chronically infected animals. Despite a considerable level of variability at the quasispecies level, no indications were found for any predictive pattern with regard to the chronicity of the CSFV infections. While there might be no direct marker for chronicity, moderate virulence of some CSFV strains in itself seems to be a crucial prerequisite for the establishment of long-term infections which does not need further genetic adaption. Thus, general host and virus factors need further investigation.

Bluetongue virus serotype 27: detection and characterization of two novel variants in Corsica, France.

During the compulsory vaccination programme against bluetongue virus serotype 1 (BTV-1) in Corsica (France) in 2014, a BTV strain belonging to a previously uncharacterized serotype (BTV-27) was isolated from asymptomatic goats. The present study describes the detection and molecular characterization of two additional distinct BTV-27 variants found in goats in Corsica in 2014 and 2015. The full coding genome of these two novel BTV-27 variants show high homology (90-93 % nucleotide/93-95 % amino acid) with the originally described BTV-27 isolate from Corsican goats in 2014. These three variants constitute the novel serotype BTV-27 ('BTV-27/FRA2014/v01 to v03'). Phylogenetic analyses with the 26 other established BTV serotypes revealed the closest relationship to BTV-25 (SWI2008/01) (80 % nucleotide/86 % amino acid) and to BTV-26 (KUW2010/02) (73-74 % nucleotide/80-81 % amino acid). However, highest sequence homologies between individual segments of BTV-27/FRA2014/v01-v03 with BTV-25 and BTV-26 vary. All three variants share the same segment 2 nucleotype with BTV-25. Neutralization assays of anti-BTV27/FRA2014/v01-v03 sera with a reassortant virus containing the outer capsid proteins of BTV-25 (BTV1VP2/VP5 BTV25) further confirmed that BTV-27 represents a distinct BTV serotype. Relationships between the variants and with BTV-25 and BTV-26, hypotheses about their origin, reassortment events and evolution are discussed.

Out of the Reservoir: Phenotypic and Genotypic Characterization of a Novel Cowpox Virus Isolated from a Common Vole.

UNLABELLED: The incidence of human cowpox virus (CPXV) infections has increased significantly in recent years. Serological surveys have suggested wild rodents as the main CPXV reservoir. We characterized a CPXV isolated during a large-scale screening from a feral common vole. A comparison of the full-length DNA sequence of this CPXV strain with a highly virulent pet rat CPXV isolate showed a sequence identity of 96%, including a large additional open reading frame (ORF) of about 6,000 nucleotides which is absent in the reference CPXV strain Brighton Red. Electron microscopy analysis demonstrated that the vole isolate, in contrast to the rat strain, forms A-type inclusion (ATI) bodies with incorporated virions, consistent with the presence of complete ati and p4c genes. Experimental infections showed that the vole CPXV strain caused only mild clinical symptoms in its natural host, while all rats developed severe respiratory symptoms followed by a systemic rash. In contrast, common voles infected with a high dose of the rat CPXV showed severe signs of respiratory disease but no skin lesions, whereas infection with a low dose led to virus excretion with only mild clinical signs. We concluded that the common vole is susceptible to infection with different CPXV strains. The spectrum ranges from well-adapted viruses causing limited clinical symptoms to highly virulent strains causing severe respiratory symptoms. In addition, the low pathogenicity of the vole isolate in its eponymous host suggests a role of common voles as a major CPXV reservoir, and future research will focus on the correlation between viral genotype and phenotype/pathotype in accidental and reservoir species. IMPORTANCE: We report on the first detection and isolation of CPXV from a putative reservoir host, which enables comparative analyses to understand the infection cycle of these zoonotic orthopox viruses and the relevant genes involved. In vitro studies, including whole-genome sequencing as well as in vivo experiments using the Wistar rat model and the vole reservoir host allowed us to establish links between genomic sequences and the in vivo properties (virulence) of the novel vole isolate in comparison to those of a recent zoonotic CPXV isolated from pet rats in 2009. Furthermore, the role of genes present only in a reservoir isolate can now be further analyzed. These studies therefore allow unique insights and conclusions about the role of the rodent reservoir in CPXV epidemiology and transmission and about the zoonotic threat that these viruses represent.

Characterization of Shuni viruses detected in Israel.

Shuni virus (SHUV) was recently identified in Israel in several brains of ovine, bovine, and goat fetuses and newborn animals with congenital arthrogryposis-hydranencephaly syndrome. In the present study, the sequences of several Israeli SHUV strains were analyzed in detail; based on the small genome segment which encodes the nucleocapsid protein and the small nonstructural protein (NSs), a very high similarity of 99-100 % among each other was found. In contrast to the highly conserved N protein, several mutations were found within the NSs-coding sequence of SHUVs present in brain samples of malformed fetuses, resulting in a considerably frequent appearance of stop codons. Interferon alpha/beta production was demonstrated in an in-vitro interferon bioassay; hence, the virus isolated from the brain of a malformed sheep fetus acquired mutations, resulting in the loss of its NSs protein function.

Comparison of porcine epidemic diarrhea viruses from Germany and the United States, 2014.

Since 2013, highly virulent porcine epidemic diarrhea virus has caused considerable economic losses in the United States. To determine the relation of US strains to those recently causing disease in Germany, we compared genomes and found that the strain from Germany is closely related to variants in the United States.

Mixed triple: allied viruses in unique recent isolates of highly virulent type 2 bovine viral diarrhea virus detected by deep sequencing.

UNLABELLED: In February 2013, very severe acute clinical symptoms were observed in calves, heifers, and dairy cattle in several farms in North Rhine Westphalia and Lower Saxony, Germany. Deep sequencing revealed the coexistence of three distinct genome variants within recent highly virulent bovine viral diarrhea virus type 2 (BVDV-2) isolates. While the major portion (ca. 95%) of the population harbored a duplication of a 222-nucleotide (nt) segment within the p7-NS2-encoding region, the minority reflected the standard structure of a BVDV-2 genome. Additionally, unusual mutations were found in both variants, within the highly conserved p7 protein and close to the p7-NS2 cleavage site. Using a reverse genetic system with a BVDV-2a strain harboring a similar duplication, it could be demonstrated that during replication, genomes without duplication are generated de novo from genomes with duplication. The major variant with duplication is compulsorily escorted by the minor variant without duplication. RNA secondary structure prediction allowed the analysis of the unique but stable mixture of three BVDV variants and also provided the explanation for their generation. Finally, our results suggest that the variant with duplication plays the major role in the highly virulent phenotype. IMPORTANCE: This study emphasizes the importance of full-genome deep sequencing in combination with manual in-depth data analysis for the investigation of viruses in basic research and diagnostics. Here we investigated recent highly virulent bovine viral diarrhea virus isolates from a 2013 series of outbreaks. We discovered a unique special feature of the viral genome, an unstable duplication of 222 nucleotides which is eventually deleted by viral polymerase activity, leading to an unexpectedly mixed population of viral genomes for all investigated isolates. Our study is of high importance to the field because we demonstrate that these insertion/deletion events allow another level of genome plasticity of plus-strand RNA viruses, in addition to the well-known polymerase-induced single nucleotide variations which are generally considered the main basis for viral adaptation and evolution.

Novel bluetongue virus in goats, Corsica, France, 2014.

During 2000-2013, 4 genotypes of bluetongue virus (BTV) were detected in Corsica, France. At the end of 2013, a compulsory BTV-1 vaccination campaign was initiated among domestic ruminants; biological samples from goats were tested as part of a corresponding monitoring program. A BTV strain with nucleotide sequences suggestive of a novel serotype was detected.

Preparation of modified long-mer RNAs and analysis of FMN binding to the ypaA aptamer from B. subtilis.

In recent years, RNA has been shown to fulfil a number of cellular functions. This has led to much interest in elucidation of the structure of functional RNA molecules, and thus, in the preparation of suitably functionalized RNAs. The chemical synthesis of RNAs allows for the site-specific modification; however, is limited to sequences of about 60-70 nucleotides in length. At the example of the flavine mononucleotide (FMN) responsive aptamer of the ypaA riboswitch from B. subtilis, we demonstrate the highly efficient preparation of site-specifically modified long-mer RNAs. Our strategy consists of the chemical synthesis of fragments followed by enzymatic or chemical ligation. Splint ligation with T4 RNA ligase turned out to be most successful among the enyzymatic protocols. Highly efficient chemical ligation was performed by azide-alkyne cycloaddition of suitably modified RNA fragments. Wild-type and 2-aminopurine (2-AP)-modified variants of the ypaA aptamer were prepared. FMN binding to all synthesized ypaA aptamer variants is demonstrated. However, dissociation of FMN from its binding site by reduction of the isoalloxazin unit as demonstrated before for a small-hairpin-derived aptazyme could not be shown. This implies that either FMN is less accessible to reduction when it is bound to its natural aptamer; that reduced FMN remains bound to the aptamer; or that FMN upon reduction indeed is released from its binding site, without the aptamer folding back in the natural ligand-free state. The results of this study are of general interest to the preparation of site-specifically modified RNAs for investigation into structure and function.

First Detection and Circulation of RHDV2 in New Zealand.

Rabbit haemorrhage disease virus 2 (RHDV2) is a highly pathogenic lagovirus that causes lethal disease in rabbits and hares (lagomorphs). Since its first detection in Europe in 2010, RHDV2 has spread worldwide and has been detected in over 35 countries so far. Here, we provide the first detailed report of the detection and subsequent circulation of RHDV2 in New Zealand. RHDV2 was first detected in New Zealand in 2018, with positive samples retrospectively identified in December 2017. Subsequent time-resolved phylogenetic analysis suggested a single introduction into the North Island between March and November 2016. Genetic analysis identified a GI.3P-GI.2 variant supporting a non-Australian origin for the incursion; however, more accurate identification of the source of the incursion remains challenging due to the wide global distribution of the GI.3P-GI.2 variant. Furthermore, our analysis suggests the spread of the virus between the North and South Islands of New Zealand at least twice, dated to mid-2017 and around 2018. Further phylogenetic analysis also revealed a strong phylogeographic pattern. So far, no recombination events with endemic benign New Zealand rabbit caliciviruses have been identified. This study highlights the need for further research and surveillance to monitor the distribution and diversity of lagoviruses in New Zealand and to detect incursions of novel variants.

First Detection of Benign Rabbit Caliciviruses in Chile.

Pathogenic lagoviruses (Rabbit hemorrhagic disease virus, RHDV) are widely spread across the world and are used in Australia and New Zealand to control populations of feral European rabbits. The spread of the non-pathogenic lagoviruses, e.g., rabbit calicivirus (RCV), is less well studied as the infection results in no clinical signs. Nonetheless, RCV has important implications for the spread of RHDV and rabbit biocontrol as it can provide varying levels of cross-protection against fatal infection with pathogenic lagoviruses. In Chile, where European rabbits are also an introduced species, myxoma virus was used for localised biocontrol of rabbits in the 1950s. To date, there have been no studies investigating the presence of lagoviruses in the Chilean feral rabbit population. In this study, liver and duodenum rabbit samples from central Chile were tested for the presence of lagoviruses and positive samples were subject to whole RNA sequencing and subsequent data analysis. Phylogenetic analysis revealed a novel RCV variant in duodenal samples that likely originated from European RCVs. Sequencing analysis also detected the presence of a rabbit astrovirus in one of the lagovirus-positive samples.

Utilizing Molecular Epidemiology and Citizen Science for the Surveillance of Lagoviruses in Australia.

Australia has multiple lagoviruses with differing pathogenicity. The circulation of these viruses was traditionally determined through opportunistic sampling events. In the lead up to the nationwide release of RHDVa-K5 (GI.1aP-GI.1a) in 2017, an existing citizen science program, RabbitScan, was augmented to allow members of the public to submit samples collected from dead leporids for lagovirus testing. This study describes the information obtained from the increased number of leporid samples received between 2015 and 2022 and focuses on the recent epidemiological interactions and evolutionary trajectory of circulating lagoviruses in Australia between October 2020 and December 2022. A total of 2771 samples were tested from January 2015 to December 2022, of which 1643 were lagovirus-positive. Notable changes in the distribution of lagovirus variants were observed, predominantly in Western Australia, where RHDV2-4c (GI.4cP-GI.2) was detected again in 2021 after initially being reported to be present in 2018. Interestingly, we found evidence that the deliberately released RHDVa-K5 was able to establish and circulate in wild rabbit populations in WA. Overall, the incorporation of citizen science approaches proved to be a cost-efficient method to increase the sampling area and enable an in-depth analysis of lagovirus distribution, genetic diversity, and interactions. The maintenance of such programs is essential to enable continued investigations of the critical parameters affecting the biocontrol of feral rabbit populations in Australia, as well as to enable the detection of any potential future incursions.

Pathogen profiling of Australian rabbits by metatranscriptomic sequencing.

Australia is known for its long history of using biocontrol agents, such as myxoma virus (MYXV) and rabbit haemorrhagic disease virus (RHDV), to manage wild European rabbit populations. Interestingly, while undertaking RHDV surveillance of rabbits that were found dead, we observed that approximately 40% of samples were negative for RHDV. To investigate whether other infectious agents are responsible for killing rabbits in Australia, we subjected a subset of these RHDV-negative liver samples to metatranscriptomic sequencing. In addition, we investigated whether the host transcriptome data could provide additional differentiation between likely infectious versus non-infectious causes of death. We identified transcripts from several Clostridia species, Pasteurella multocida, Pseudomonas spp., and Eimeria stiedae, in liver samples of several rabbits that had died suddenly, all of which are known to infect rabbits and are capable of causing disease and mortality. In addition, we identified Hepatitis E virus and Cyniclomyces yeast in some samples, both of which are not usually associated with severe disease. In one-third of the sequenced total liver RNAs, no infectious agent could be identified. While metatranscriptomic sequencing cannot provide definitive evidence of causation, additional host transcriptome analysis provided further insights to distinguish between pathogenic microbes and commensals or environmental contaminants. Interestingly, three samples where no pathogen could be identified showed evidence of up-regulated host immune responses, while immune response pathways were not up-regulated when E. stiedae, Pseudomonas, or yeast were detected. In summary, although no new putative rabbit pathogens were identified, this study provides a robust workflow for future investigations into rabbit mortality events.

Lagovirus Non-structural Protein p23: A Putative Viroporin That Interacts With Heat Shock Proteins and Uses a Disulfide Bond for Dimerization.

The exact function(s) of the lagovirus non-structural protein p23 is unknown as robust cell culture systems for the Rabbit haemorrhagic disease virus (RHDV) and other lagoviruses have not been established. Instead, a range of in vitro and in silico models have been used to study p23, revealing that p23 oligomerizes, accumulates in the cytoplasm, and possesses a conserved C-terminal region with two amphipathic helices. Furthermore, the positional homologs of p23 in other caliciviruses have been shown to possess viroporin activity. Here, we report on the mechanistic details of p23 oligomerization. Site-directed mutagenesis revealed the importance of an N-terminal cysteine for dimerization. Furthermore, we identified cellular interactors of p23 using stable isotope labeling with amino acids in cell culture (SILAC)-based proteomics; heat shock proteins Hsp70 and 110 interact with p23 in transfected cells, suggesting that they 'chaperone' p23 proteins before their integration into cellular membranes. We investigated changes to the global transcriptome and proteome that occurred in infected rabbit liver tissue and observed changes to the misfolded protein response, calcium signaling, and the regulation of the endoplasmic reticulum (ER) network. Finally, flow cytometry studies indicate slightly elevated calcium concentrations in the cytoplasm of p23-transfected cells. Taken together, accumulating evidence suggests that p23 is a viroporin that might form calcium-conducting channels in the ER membranes.

Distribution and Genetic Diversity of Hepatitis E Virus in Wild and Domestic Rabbits in Australia.

In 2020, Hepatitis E virus (HEV) was detected for the first time in Australian rabbits. To improve our understanding of the genetic diversity and distribution of the virus, 1635 rabbit liver samples from locations across Australia were screened via RT-qPCR for HEV. HEV genomes were amplified and sequenced from 48 positive samples. Furthermore, we tested 380 serum samples from 11 locations across Australia for antibodies against HEV. HEV was detected in rabbits from all states and territories, except the Northern Territory. Seroprevalence varied between locations (from 0% to 22%), demonstrating that HEV is widely distributed in rabbit populations across Australia. Phylogenetic analyses showed that Australian HEV sequences are genetically diverse and that HEV was likely introduced into Australia independently on several occasions. In summary, this study broadens our understanding of the genetic diversity of rabbit HEV globally and shows that the virus is endemic in both domestic and wild rabbit populations in Australia.