Emergence and spread of SARS-CoV-2 variants from farmed mink to humans and back during the epidemic in Denmark, June-November 2020.

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) not only caused the COVID-19 pandemic but also had a major impact on farmed mink production in several European countries. In Denmark, the entire population of farmed mink (over 15 million animals) was culled in late 2020. During the period of June to November 2020, mink on 290 farms (out of about 1100 in the country) were shown to be infected with SARS-CoV-2. Genome sequencing identified changes in the virus within the mink and it is estimated that about 4000 people in Denmark became infected with these mink virus variants. However, the routes of transmission of the virus to, and from, the mink have been unclear. Phylogenetic analysis revealed the generation of multiple clusters of the virus within the mink. Detailed analysis of changes in the virus during replication in mink and, in parallel, in the human population in Denmark, during the same time period, has been performed here. The majority of cases in mink involved variants with the Y453F substitution and the H69/V70 deletion within the Spike (S) protein; these changes emerged early in the outbreak. However, further introductions of the virus, by variants lacking these changes, from the human population into mink also occurred. Based on phylogenetic analysis of viral genome data, we estimate, using a conservative approach, that about 17 separate examples of mink to human transmission occurred in Denmark but up to 59 such events (90% credible interval: (39-77)) were identified using parsimony to count cross-species jumps on transmission trees inferred using Bayesian methods. Using the latter approach, 136 jumps (90% credible interval: (117-164)) from humans to mink were found, which may underlie the farm-to-farm spread. Thus, transmission of SARS-CoV-2 from humans to mink, mink to mink, from mink to humans and between humans were all observed.

Infection, recovery and re-infection of farmed mink with SARS-CoV-2.

Mink, on a farm with about 15,000 animals, became infected with SARS-CoV-2. Over 75% of tested animals were positive for SARS-CoV-2 RNA in throat swabs and 100% of tested animals were seropositive. The virus responsible had a deletion of nucleotides encoding residues H69 and V70 within the spike protein gene as well as the A22920T mutation, resulting in the Y453F substitution within this protein, seen previously in mink. The infected mink recovered and after free-testing of 300 mink (a level giving 93% confidence of detecting a 1% prevalence), the animals remained seropositive. During further follow-up studies, after a period of more than 2 months without any virus detection, over 75% of tested animals again scored positive for SARS-CoV-2 RNA. Whole genome sequencing showed that the viruses circulating during this re-infection were most closely related to those identified in the first outbreak on this farm but additional sequence changes had occurred. Animals had much higher levels of anti-SARS-CoV-2 antibodies in serum samples after the second round of infection than at free-testing or during recovery from initial infection, consistent with a boosted immune response. Thus, it was concluded that following recovery from an initial infection, seropositive mink were readily re-infected by SARS-CoV-2.

SARS-CoV-2 Transmission between Mink (Neovison vison) and Humans, Denmark.

Severe acute respiratory syndrome coronavirus 2 has caused a pandemic in humans. Farmed mink (Neovison vison) are also susceptible. In Denmark, this virus has spread rapidly among farmed mink, resulting in some respiratory disease. Full-length virus genome sequencing revealed novel virus variants in mink. These variants subsequently appeared within the local human community.

Analysis of Virus Population Profiles within Pigs Infected with Virulent Classical Swine Fever Viruses: Evidence for Bottlenecks in Transmission but Absence of Tissue-Specific Virus Variants.

Classical swine fever virus (CSFV) contains a specific motif within the E2 glycoprotein that differs between strains of different virulence. In the highly virulent CSFV strain Koslov, this motif comprises residues S763/L764 in the polyprotein. However, L763/P764 represent the predominant alleles in published CSFV genomes. In this study, changes were introduced into the CSFV strain Koslov (here called vKos_SL) to generate modified CSFVs with substitutions at residues 763 and/or 764 (vKos_LL, vKos_SP, and vKos_LP). The properties of these mutant viruses, in comparison to those of vKos_SL, were determined in pigs. Each of the viruses was virulent and induced typical clinical signs of CSF, but the vKos_LP strain produced them significantly earlier. Full-length CSFV cDNA amplicons (12.3 kb) derived from sera of infected pigs were deep sequenced and cloned to reveal the individual haplotypes that contributed to the single-nucleotide polymorphism (SNP) profiles observed in the virus population. The SNP profiles for vKos_SL and vKos_LL displayed low-level heterogeneity across the entire genome, whereas vKos_SP and vKos_LP displayed limited diversity with a few high-frequency SNPs. This indicated that vKos_SL and vKos_LL exhibited a higher level of fitness in the host and more stability at the consensus level, whereas several consensus changes were observed in the vKos_SP and vKos_LP sequences, pointing to adaptation. For each virus, only a subset of the variants present within the virus inoculums were maintained in the infected pigs. No clear tissue-dependent quasispecies differentiation occurred within inoculated pigs; however, clear evidence for transmission bottlenecks to contact animals was observed, with subsequent loss of sequence diversity.IMPORTANCE The surface-exposed E2 protein of classical swine fever virus is required for its interaction with host cells. A short motif within this protein varies between strains of different virulence. The importance of two particular amino acid residues in determining the properties of a highly virulent strain of the virus has been analyzed. Each of the different viruses tested proved highly virulent, but one of them produced earlier, but not more severe, disease. By analyzing the virus genomes present within infected pigs, it was found that the viruses which replicated within inoculated animals were only a subset of those within the virus inoculum. Furthermore, following contact transmission, it was shown that a very restricted set of viruses had transferred between animals. There were no significant differences in the virus populations present in various tissues of the infected animals. These results indicate mechanisms of virus population change during transmission between animals.

Modifications of the nucleoprotein of viral haemorrhagic septicaemia virus showed gain of virulence in intraperitoneally infected rainbow trout.

Viral haemorrhagic septicaemia virus (VHSV) is the cause of an important listed disease in European rainbow trout (Oncorhynchus mykiss) aquaculture and can be present in a wide range of fish species, including marine fish, which can act as viral reservoir. Recent studies revealed putative genetic virulence markers of VHSV to rainbow trout highlighting the roles of the nucleoprotein, phosphoprotein and non-virion protein. Using reverse genetics, we produced recombinant viruses by introducing parts of or the entire nucleoprotein from a high-virulent isolate VHSV into a low-virulent backbone. Furthermore, we also made recombinant viruses by introducing residue modifications in the nucleoprotein that seem to play a role in virulence. Rainbow trout challenged with these recombinant viruses (rVHSVs) by intraperitoneal injection (IP) developed clinical signs and showed lower survival when compared to the parental rVHSV whereas fish challenged by immersion did not show clinical signs except for the high-virulent control. The mutations did not influence the viral growth in cell culture. The recombinant viruses and parental recombinant were unable to replicate and show cytopathic effect in EPC cells whereas the high-virulent control was well adapted in all the fish cell lines tested. We showed evidence that corroborates with the hypothesis that the nucleoprotein has virulence motifs associated with VHSV virulence in rainbow trout.

Preliminary report of an outbreak of SARS-CoV-2 in mink and mink farmers associated with community spread, Denmark, June to November 2020.

In June-November 2020, SARS-CoV-2-infected mink were detected in 290 of 1,147 Danish mink farms. In North Denmark Region, 30% (324/1,092) of people found connected to mink farms tested SARS-CoV-2-PCR-positive and approximately 27% (95% confidence interval (CI): 25-30) of SARS-CoV-2-strains from humans in the community were mink-associated. Measures proved insufficient to mitigate spread. On 4 November, the government ordered culling of all Danish mink. Farmed mink constitute a potential virus reservoir challenging pandemic control.

Identification of specific amino acid residues in the border disease virus glycoprotein E2 that modify virus growth in pig cells but not in sheep cells.

Border disease virus (BDV) envelope glycoprotein E2 is required for entry into cells and is a determinant of host tropism for sheep and pig cells. Here, we describe adaptive changes in the BDV E2 protein that modify virus replication in pig cells. To achieve this, two BDV isolates, initially collected from a pig and a sheep on the same farm, were passaged in primary sheep and pig cells in parallel with a rescued variant of the pig virus derived from a cloned full-length BDV cDNA. The pig isolate and the rescued virus shared the same amino acid sequence, but the sheep isolate differed at ten residues, including two substitutions in E2 (K771E and Y925H). During serial passage in cells, the viruses displayed clear selectivity for growth in sheep cells; only the cDNA-derived virus adapted to grow in pig cells. Sequencing revealed an amino acid substitution (Q739R) in the E2 domain DA of this rescued virus. Adaptation at the same residue (Q739K/Q739R) was also observed after passaging of the pig isolate in sheep cells. Use of reverse genetics confirmed that changing residue Q739 to R or K (each positively charged) was sufficient to achieve adaptation to pig cells. Furthermore, this change in host tropism was suppressed if Q739R was combined with K771E. Another substitution (Q728R), conferring an additional positive charge, acquired during passaging, restored the growth of the Q739R/K771E variant. Overall, this study provided evidence that specific, positively charged, residues in the E2 domain DA are crucial for pig-cell tropism of BDV.

Distinct roles for the IIId2 sub-domain in pestivirus and picornavirus internal ribosome entry sites.

Viral internal ribosomes entry site (IRES) elements coordinate the recruitment of the host translation machinery to direct the initiation of viral protein synthesis. Within hepatitis C virus (HCV)-like IRES elements, the sub-domain IIId(1) is crucial for recruiting the 40S ribosomal subunit. However, some HCV-like IRES elements possess an additional sub-domain, termed IIId2, whose function remains unclear. Herein, we show that IIId2 sub-domains from divergent viruses have different functions. The IIId2 sub-domain present in Seneca valley virus (SVV), a picornavirus, is dispensable for IRES activity, while the IIId2 sub-domains of two pestiviruses, classical swine fever virus (CSFV) and border disease virus (BDV), are required for 80S ribosomes assembly and IRES activity. Unlike in SVV, the deletion of IIId2 from the CSFV and BDV IRES elements impairs initiation of translation by inhibiting the assembly of 80S ribosomes. Consequently, this negatively affects the replication of CSFV and BDV. Finally, we show that the SVV IIId2 sub-domain is required for efficient viral RNA synthesis and growth of SVV, but not for IRES function. This study sheds light on the molecular evolution of viruses by clearly demonstrating that conserved RNA structures, within distantly related RNA viruses, have acquired different roles in the virus life cycles.

Strategy for efficient generation of numerous full-length cDNA clones of classical swine fever virus for haplotyping.

BACKGROUND: Direct molecular cloning of full-length cDNAs derived from viral RNA is an approach to identify the individual viral genomes within a virus population. This enables characterization of distinct viral haplotypes present during infection. RESULTS: In this study, we recover individual genomes of classical swine fever virus (CSFV), present in a pig infected with vKos that was rescued from a cDNA clone corresponding to the highly virulent CSFV Koslov strain. Full-length cDNA amplicons (ca. 12.3 kb) were made by long RT-PCR, using RNA extracted from serum, and inserted directly into a cloning vector prior to detailed characterization of the individual viral genome sequences. The amplicons used for cloning were deep sequenced, which revealed low level sequence variation (< 5%) scattered across the genome consistent with the clone-derived origin of vKos. Numerous full-length cDNA clones were generated using these amplicons and full-genome sequencing of individual cDNA clones revealed insights into the virus diversity and the haplotypes present during infection. Most cDNA clones were unique, containing several single-nucleotide polymorphisms, and phylogenetic reconstruction revealed a low degree of order. CONCLUSIONS: This optimized methodology enables highly efficient construction of full-length cDNA clones corresponding to individual viral genomes present within RNA virus populations.

Determinants of the VP1/2A junction cleavage by the 3C protease in foot-and-mouth disease virus-infected cells.

The foot-and-mouth disease virus (FMDV) capsid precursor, P1-2A, is cleaved by FMDV 3C protease to yield VP0, VP3, VP1 and 2A. Cleavage of the VP1/2A junction is the slowest. Serotype O FMDVs with uncleaved VP1-2A (having a K210E substitution in VP1; at position P2 in cleavage site) have been described previously and acquired a second site substitution (VP1 E83K) during virus rescue. Furthermore, introduction of the VP1 E83K substitution alone generated a second site change at the VP1/2A junction (2A L2P, position P2' in cleavage site). These virus adaptations have now been analysed using next-generation sequencing to determine sub-consensus level changes in the virus; this revealed other variants within the E83K mutant virus population that changed residue VP1 K210. The construction of serotype A viruses with a blocked VP1/2A cleavage site (containing K210E) has now been achieved. A collection of alternative amino acid substitutions was made at this site, and the properties of the mutant viruses were determined. Only the presence of a positively charged residue at position P2 in the cleavage site permitted efficient cleavage of the VP1/2A junction, consistent with analyses of diverse FMDV genome sequences. Interestingly, in contrast to the serotype O virus results, no second site mutations occurred within the VP1 coding region of serotype A viruses with the blocked VP1/2A cleavage site. However, some of these viruses acquired changes in the 2C protein that is involved in enterovirus morphogenesis. These results have implications for the testing of potential antiviral agents targeting the FMDV 3C protease.

Identification of rhabdoviral sequences in oropharyngeal swabs from German and Danish bats.

BACKGROUND: In the frame of active lyssavirus surveillance in bats, oropharyngeal swabs from German (N = 2297) and Danish (N = 134) insectivorous bats were investigated using a newly developed generic pan-lyssavirus real-time reverse transcriptase PCR (RT-qPCR). FINDINGS: In total, 15 RT-qPCR positive swabs were detected. Remarkably, sequencing of positive samples did not confirm the presence of bat associated lyssaviruses but revealed nine distinct novel rhabdovirus-related sequences. CONCLUSIONS: Several novel rhabdovirus-related sequences were detected both in German and Danish insectivorous bats. The results also prove that the novel generic pan-lyssavirus RT-qPCR offers a very broad detection range that allows the collection of further valuable data concerning the broad and complex diversity within the family Rhabdoviridae.

Studies on genetic diversity of bovine viral diarrhea viruses in Danish cattle herds.

Scandinavian countries have successfully pursued bovine viral diarrhea virus (BVDV) eradication without the use of vaccines. In Denmark, control and eradication of BVDV were achieved during the last two decades, but occasionally new BVDV infections are detected in some Danish cattle herds. The aim of this study was to determine recent BVDV subtypes isolated from 4 Danish herds (A, B, C, and D) isolated in 2009-2012 and to analyze the genetic variation of these isolates within the same herd and its relation with those of other herds. The results showed that three herds (B, C, D) were BVDV 1-b and only one herd (herd A) was BVDV 1-d, no other subtypes were detected. The deduced E2 amino acids result showed a high identity percent (99-100 %) between isolates originating from the same herd, but with higher variation compared to isolates of the other herds. Some of these new Danish strains have closer relationship to BVDVs from outside Denmark than to older Danish strains indicating that these are new introductions to Denmark. In conclusion, BVDV-1 subtypes recently detected in Denmark were only subtypes 1b and 1d, and BVDV infections established in a herd is genetically stable over a long time period.

Porcine respiratory coronavirus genome sequences; comparisons and relationships to transmissible gastroenteritis viruses.

Porcine respiratory coronavirus (PRCV) was initially detected in Europe, and later in the United States of America (US), in the 1980s. In this study we obtained and compared PRCV sequences from Europe and the US, and investigated how these are related to transmissible gastroenteritis virus (TGEV) sequences. The whole genome sequences of Danish (1/90-DK), Italian (PRCV15087/12 III NPTV Parma), and Belgian PRCV (91V44) strains are presented. These sequences were aligned with nine other PRCV sequences from Europe and the US, and 43 TGEV sequences. Following alignment of the PRCV sequences, it was apparent that multiple amino acid variations in the structural proteins were distinct between the European and US strains. The alignments were used to build phylogenetic trees to infer the evolutionary relationships between the strains. In these trees, the European PRCV strains clustered as a separate group, whereas the US strains of PRCV all clustered with TGEVs.

A Deep Sequencing Strategy for Investigation of Virus Variants within African Swine Fever Virus-Infected Pigs.

African swine fever virus (ASFV) is the causative agent of African swine fever, an economically important disease of pigs, often with a high case fatality rate. ASFV has demonstrated low genetic diversity among isolates collected within Eurasia. To explore the influence of viral variants on clinical outcomes and infection dynamics in pigs experimentally infected with ASFV, we have designed a deep sequencing strategy. The variant analysis revealed unique SNPs at <10% frequency in several infected pigs as well as some SNPs that were found in more than one pig. In addition, a deletion of 10,487 bp (resulting in the complete loss of 21 genes) was present at a nearly 100% frequency in the ASFV DNA from one pig at position 6362-16849. This deletion was also found to be present at low levels in the virus inoculum and in two other infected pigs. The current methodology can be used for the currently circulating Eurasian ASFVs and also adapted to other ASFV strains and genotypes. Comprehensive deep sequencing is critical for following ASFV molecular evolution, especially for the identification of modifications that affect virus virulence.

Microbiological and decomposition analysis of mass mink burial sites during the COVID-19 pandemic.

In 2020, Denmark buried approximately four million culled, farmed mink in mass graves treated with slaked lime due to widespread SARS-CoV-2 infections. After six months, environmental concerns prompted the exhumation of these cadavers. Our analysis encompassed visual inspections, soil pH measurements, and gas emission assessments of the grave environment. Additionally, we evaluated carcasses for decay status, cadaverine content, and the presence of various pathogens, including SARS-CoV-2 and mink coronavirus. Our findings revealed minimal microbial activity and limited carcass decomposition. Although viral RNA from SARS-CoV-2 and mink coronavirus, along with DNA from Aleutian mink disease virus, were detected, the absence of infectious SARS-CoV-2 in cell culture assays suggests slow natural degradation processes. This study provides critical insights for future considerations in managing mass burial scenarios during outbreaks of livestock-associated zoonotic pathogens.

Efficient generation of recombinant RNA viruses using targeted recombination-mediated mutagenesis of bacterial artificial chromosomes containing full-length cDNA.

BACKGROUND: Infectious cDNA clones are a prerequisite for directed genetic manipulation of RNA viruses. Here, a strategy to facilitate manipulation and rescue of classical swine fever viruses (CSFVs) from full-length cDNAs present within bacterial artificial chromosomes (BACs) is described. This strategy allows manipulation of viral cDNA by targeted recombination-mediated mutagenesis within bacteria. RESULTS: A new CSFV-BAC (pBeloR26) derived from the Riems vaccine strain has been constructed and subsequently modified in the E2 coding sequence, using the targeted recombination strategy to enable rescue of chimeric pestiviruses (vR26_E2gif and vR26_TAV) with potential as new marker vaccine candidates. Sequencing of the BACs revealed a high genetic stability during passages within bacteria. The complete genome sequences of rescued viruses, after extensive passages in mammalian cells showed that modifications in the E2 protein coding sequence were stably maintained. A single amino acid substitution (D3431G) in the RNA dependent RNA polymerase was observed in the rescued viruses vR26_E2gif and vR26, which was reversion to the parental Riems sequence. CONCLUSIONS: These results show that targeted recombination-mediated mutagenesis provides a powerful tool for expediting the construction of novel RNA genomes and should be applicable to the manipulation of other RNA viruses.

Analysis of classical swine fever virus RNA replication determinants using replicons.

Self-replicating RNAs (replicons), with or without reporter gene sequences, derived from the genome of the Paderborn strain of classical swine fever virus (CSFV) have been produced. The full-length viral cDNA, propagated within a bacterial artificial chromosome, was modified by targeted recombination within Escherichia coli. RNA transcripts were produced in vitro and introduced into cells by electroporation. The translation and replication of the replicon RNAs could be followed by the accumulation of luciferase (from Renilla reniformis or Gaussia princeps) protein expression (where appropriate), as well as by detection of CSFV NS3 protein production within the cells. Inclusion of the viral E2 coding region within the replicon was advantageous for replication efficiency. Production of chimeric RNAs, substituting the NS2 and NS3 coding regions (as a unit) from the Paderborn strain with the equivalent sequences from the highly virulent Koslov strain or the vaccine strain Riems, blocked replication. However, replacing the Paderborn NS5B coding sequence with the RNA polymerase coding sequence from the Koslov strain greatly enhanced expression of the reporter protein from the replicon. In contrast, replacement with the Riems NS5B sequence significantly impaired replication efficiency. Thus, these replicons provide a system for determining specific regions of the CSFV genome required for genome replication without the constraints of maintaining infectivity.

Modulation of translation initiation efficiency in classical swine fever virus.

Modulation of translation initiation efficiency on classical swine fever virus (CSFV) RNA can be achieved by targeted mutations within the internal ribosome entry site (IRES). In this study, cDNAs corresponding to the wild-type (wt) or mutant forms of the IRES of CSFV strain Paderborn were amplified and inserted into dicistronic reporter plasmids encoding Fluc and Rluc under the control of a T7 promoter. The mutations were within domains II, IIId(1), and IIIf of the IRES. The plasmids were transfected into baby hamster kidney (BHK) cells infected with recombinant vaccinia virus vTF7-3, which expresses the T7 RNA polymerase. IRES mutants with different levels of IRES activity were identified and then introduced by homologous recombination into bacterial artificial chromosomes (BACs) containing CSFV Paderborn cDNA downstream of a T7 promoter. From the wt and mutant BACs, full-length CSFV RNA transcripts were produced in vitro and electroporated into porcine PK15 cells. Rescued mutant viruses were obtained from RNAs that contained mutations within domain IIIf which retained more than 75% of the wt translation efficiency. Sequencing of cDNA generated from these rescued viruses verified the maintenance of the introduced changes within the IRES. The growth characteristics of each rescued mutant virus were compared to those of the wt virus. It was shown that viable mutant viruses with reduced translation initiation efficiency can be designed and generated and that viruses containing mutations within domain IIIf of the IRES have reduced growth in cell culture compared to the wt virus.

Transplacental transmission of field and rescued strains of BTV-2 and BTV-8 in experimentally infected sheep.

Transplacental transmission of bluetongue virus has been shown previously for the North European strain of serotype 8 (BTV-8) and for tissue culture or chicken egg-adapted vaccine strains but not for field strains of other serotypes. In this study, pregnant ewes (6 per group) were inoculated with either field or rescued strains of BTV-2 and BTV-8 in order to determine the ability of these viruses to cross the placental barrier. The field BTV-2 and BTV-8 strains was passaged once in Culicoides KC cells and once in mammalian cells. All virus inoculated sheep became infected and seroconverted against the different BTV strains used in this study. BTV RNA was detectable in the blood of all but two ewes for over 28 days but infectious virus could only be detected in the blood for a much shorter period. Interestingly, transplacental transmission of BTV-2 (both field and rescued strains) was demonstrated at high efficiency (6 out of 13 lambs born to BTV-2 infected ewes) while only 1 lamb of 12 born to BTV-8 infected ewes showed evidence of in utero infection. In addition, evidence for horizontal transmission of BTV-2 between ewes was observed. As expected, the parental BTV-2 and BTV-8 viruses and the viruses rescued by reverse genetics showed very similar properties to each other. This study showed, for the first time, that transplacental transmission of BTV-2, which had been minimally passaged in cell culture, can occur; hence such transmission might be more frequent than previously thought.

Increased Presence of Circulating Cell-Free, Fragmented, Host DNA in Pigs Infected with Virulent African Swine Fever Virus.

African swine fever virus (ASFV) causes severe hemorrhagic disease in domestic pigs and wild boar, often with high case fatality rates. The virus replicates in the circulating cells of the monocyte-macrophage lineage and within lymphoid tissues. The infection leads to high fever and a variety of clinical signs. In this study, it was observed that ASFV infection in pigs resulted in a >1000-fold increase in the level of circulating cell-free DNA (cfDNA), derived from the nuclei of host cells in the serum. This change occurred in parallel with the increase in circulating ASFV DNA. In addition, elevated levels (about 30-fold higher) of host mitochondrial DNA (mtDNA) were detected in the serum from ASFV-infected pigs. For comparison, the release of the cellular enzyme, lactate dehydrogenase (LDH), a commonly used marker of cellular damage, was also found to be elevated during ASFV infection, but later and less consistently. The sera from pigs infected with classical swine fever virus (CSFV), which causes a clinically similar disease to ASFV, were also tested but, surprisingly, this infection did not result in the release of cfDNA, mtDNA, or LDH. It was concluded that the level of cfDNA in the serum is a sensitive host marker of virulent ASFV infection.