Volume 27, Number 9—September 2021
Severe Acute Respiratory Syndrome Coronavirus 2 in Farmed Mink (Neovison vison), Poland
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the etiologic agent of coronavirus disease and has been spreading worldwide since December 2019. The virus can infect different animal species under experimental conditions, and mink on fur farms in Europe and other areas are susceptible to SARS-CoV-2 infection. We investigated SARS-CoV-2 infection in 91 mink from a farm in northern Poland. Using reverse transcription PCR, antigen detection, and next-generation sequencing, we confirmed that 15 animals were positive for SARS-CoV-2. We verified this finding by sequencing full viral genomes and confirmed a virus variant that has sporadic mutations through the full genome sequence in the spike protein (G75V and C1247F). We were unable to find other SARS-CoV-2 sequences simultaneously containing these 2 mutations. Country-scale monitoring by veterinary inspection should be implemented to detect SARS-CoV-2 in other mink farms.
Identifying possible pathogen hosts and studying transmission dynamics of hosts in their populations are crucial steps in controlling zoonotic diseases. The origin of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is probably bats (1), but the potential intermediate host has not yet been confirmed. SARS-CoV-2 seems to readily jump from humans to other animal species, particularly carnivores (i.e., dogs, cats, ferrets, lions, pumas) (2,3), raising concerns about new animal sources of coronavirus disease (COVID-19) (4,5).
SARS-CoV-2 infections in mink were reported from farms in Denmark and the Netherlands and later in other regions (6–9) (Figure 1). Because of SARS-CoV-2 outbreaks in mink farms and their appearance in the surrounding communities, the European Centre for Disease Prevention and Control and the World Health Organization have emphasized the need for surveying the host‒animal interface by collaboration among virologists and epidemiologists to track and characterize viral mutations (11). After SARS-CoV-2 infections in mink in the Netherlands, the Dutch Ministry of Agriculture decided to cull the mink from all farms. In Denmark, the Danish National Institute of Public Health announced the culling of all 17 million mink in the country after the virus had passed back from the mink farms into the human community. The data available for Denmark on these mink-associated SARS-CoV-2 variants suggest that these variants can spread rapidly on mink farms and in nearby human communities (12). However, humans infected with the mink-related variants do not appear to have more severe clinical symptoms than those infected with non–mink-related variants (9).
Poland is the second-largest producer of mink pelts in Europe. Poland has 810 fur animal farms, including those for foxes, mink, raccoons, dogs, and chinchillas. The 354 active mink farms in Poland contain ≈6.3 million mink. During 2019, mink farmers in Poland sold 8.5 million mink pelts (13,14).
As of May 5, 2021, Poland had recorded 2,838,180 COVID-19 cases and 70,336 total related deaths (15). Considering the recent reports of SARS-CoV-2 infections in mink in other countries in Europe and the high incidence of human SARS-CoV-2 infections in Poland, we monitored SARS-CoV-2 in mink on 1 farm in Pomorskie Voivodeship in northern Poland.
We collected throat swab (BIOCOMA, http://www.biocomma.com) specimens from 91 mink culled for pelting at a mink farm in Pomorskie Voivodeship in northern Poland on November 17, 2020. The farm owner reported no respiratory symptoms in the animals. We collected blood samples directly from the heart by using cardiac puncture and a sterile 5-mL syringe immediately after death of the mink. After separation of the blood clot, we centrifuged samples at 5,000 rpm for 10 min by using an MPW High-Speed Brushless Centrifuge (MPW Med Instruments, https://mpw.pl). We collected serum and stored it at −80°C until the samples could be analyzed.
We added 150 μL of each swab specimen sample in inactivation buffer to 300 μL of RLT lysis buffer (RNeasy Mini Kit; QIAGEN, https://www.qiagen.com). We then mixed samples vortexing and incubated them for 10 min at room temperature. After incubation, we added 400 μL of 70% ethanol to each sample and mixed them by pipetting. We transferred lysates to RNeasy Mini spin columns with collection tubes (QIAGEN) and centrifuged them at 13,000 rpm for 1 min. We then washed the columns once with 700 μL of RW1 buffer and twice with 500 μL of RPE buffer. Between every wash, we centrifuged the columns and discarded the flow-through. We performed elution by adding 50 μL of PCR-grade water to the columns and incubating them for 2 min. We placed the columns into new tubes and centrifuged them at 13,000 rpm for 1 min. After isolation, we stored the samples for <2 hours at 4°C. No human-origin samples were processed at the same time.
SARS-CoV-2 Case Definition
We defined SARS-CoV-2‒positive animals as suggested by the World Organisation for Animal Health (10). We considered mink to be SARS-CoV-2 positive if SARS-CoV-2 was isolated from a sample taken directly from an animal (nasal or oropharyngeal swab sample) or if viral nucleic acid was identified in a sample taken directly from an animal, giving cause for suspicion of previous association or contact with SARS-CoV-2, by targeting >2 specific genomic regions at a level indicating presence of infectious virus or by targeting 1 genomic region, followed by sequencing of a secondary target.
Real-Time Reverse Transcription PCR
For each sample, we prepared the reaction mixture by using a TaqPath 1-Step quantitative RT-PCR (reverse transcription PCR) Master Mix (ThermoFisher Scientific, https://www.thermofisher.com), polymerase, diethyl pyrocarbonate‒treated water (EURx, https://eurx.com.pl) and primers and probes for the RNA-dependent RNA polymerase (RdRp) and envelope (E) genes (16) in white, 8-well, quantitative PCR strips with optical clear caps (Applied Biosystems, https://www.thermofisher.com). We also prepared positive control plasmids made in-house containing the RdRp and E genes and a no-template control containing diethyl pyrocarbonate‒treated water instead of template reactions. We mixed reactions and loaded them into a Light Cycler 480 (Applied Biosystems, https://www.thermofisher.com). Cycling conditions were incubation with uracil N-glycosylase for 2 min at 25°C, reverse transcription for 15 min at 50°C, and enzyme activation for 2 min at 95°C, followed by 40 amplification cycles consisting of 3 s at 95°C and 30 s at 60°C. After each amplification cycle, we measured the signal from each sample in both the FAM (RdRp gene) and HEX (E gene) channels. Samples with a crossing point (Cp) <35 for any gene were considered positive for SARS-CoV-2.
SARS-CoV-2 Antigen Detection in Mink
We performed 2 different antigen tests to confirm the presence of viral antigen in either the swab or serum samples. Antigen tests were conducted on samples positive by quantitative PCR. We used 3 negative samples from the same batch as the positive samples as controls. A total of 150 μL of transport medium from each swab specimen was transferred to tubes from a COVID-19 Antigen Detection Kit (Zhuhai Lituo Biotechnology Co., Ltd., https://www.lituo.com.cn) containing extraction buffer. This antigen test was in the form of a cassette with a lateral flow assay that targets the nucleocapsid protein of SARS-CoV-2. Samples were mixed and incubated for 1 min at room temperature. We added 2 drops of each sample to the sample window on the test cassettes. Results consisted of 2 bands: 1 for the control and 1 for the target. If both bands showed a burgundy line, the test result was considered positive. We read results after 12 min.
We tested 91 mink serum samples by using a SARS-CoV-2 antigen ELISA (COV-04-S; Salofa Oy, https://www.salofa.com) according to the kit instructions. This test is a double-antibody sandwich ELISA. The results were obtained according to the formula based on the concentration standards provided in the kit. The cutoff value for this test was 2.97 pg/mL. The tests were repeated twice, and additional dilutions were performed to determine the final concentration as suggested in the kit instructions.
Full SARS-CoV-2 Genome Sequencing and Classification
We performed SARS-CoV-2 genome sequencing at University of Gdansk in Poland and the University of Helsinki in Finland by using only samples containing RNA isolated from virus-positive swab specimens (amplification of target gene in an RT-PCR; this gene has a Cp <35) or that were inconclusive (only 1 target gene amplification with Cp <35). Samples with positive results in the SARS-CoV-2 antigen detection assays were also sequenced.
At Gdansk, 2 independent protocols were used for SARS-CoV-2 genome sequencing. The first protocol was an Illumina (https://www.illumina.com) RNA preparation with enrichment for respiratory virus oligos panel V2, followed by an Illumina MiniSeq medium output run that produced 150-nt paired-end reads. The second protocol was an ARTIC version 3 amplicon generation (https://bmcgenomics.biomedcentral.com), followed by an Oxford Nanopore Technology MinIONB run (https://nanoporetech.com). No human origin samples were processed at the same time. No DNA/rRNA depletion methods were used. Reads were base called, de-barcoded, and trimmed to remove adaptor, barcode, and PCR primer sequences. Oxford Nanopore Technology reads were used for SARS-CoV-2 genome assembly in ARTIC-nCoV-bioinformaticsSOP-v1.1.0 (https://artic.network/ncov-2019).
Illumina paired reads were used to prepare contigs by de novo assembly by using the Geneious Prime 2020.2 (https://www.geneious.com) software suite with integrated tools for variant calling (minimum coverage = 100, minimum variant frequency = 0.25) and consensus sequence generation. The fasta files generated by the Illumina procedure were further analyzed in Kraken2 2.1.1 software (https://www.ccb.jhu.edu/software/kraken2) to classify every read to the reference database containing viral and American mink genomes (17).
In Helsinki, the sequencing libraries were prepared by using the Illumina DNA Prep Kit (New England BioLabs, https://www.neb.com). We measured library fragment sizes by using agarose gel electrophoresis and concentrations by using the Qubit dsDNA HS Assay Kit (Life Technologies, https://www.thermofisher.com) and the NEBNext Library Quant Kit for Illumina (New England BioLabs). Sequencing was conducted by using the MiSeq V3 Reagent Kit (Illumina) with 250-bp reads. We trimmed raw sequence reads and removed low-quality (quality score <30) and short (<50 nt) sequences by using Trimmomatic (18). Trimmed sequence reads were mapped against the SARS-CoV-2 reference sequence (GenBank accession no. NC_045512.2) by using BWA-MEM version 0.7.17 (19), followed by sorting and removal of duplicate reads by using SAMTools version 1.10 (20).
Phylogenetic Analysis of SARS-CoV-2 Isolates
The dataset consisted of all genetic sequences of SARS-CoV-2 from this study (Poland, Germany, Lithuania, Latvia, Estonia, Russia, and Ukraine), which was completed as the representative pool Europe by Nexstrain (https://nextstrain.org/ncov/europe) and resulted in 5,778 entries. We performed phylogenetic analysis by using the procedure recommended by Nextstrain with modifications in the subsampling region filtering procedure, in which the number of sequences per country was 40 (21). We used Augur toolkit version 10.1.1 (Nextstrain) for phylogenetic analysis and Auspice version 2.10.1 (Nextstrain) for visualization. Possible time of divergence for samples was inferred by using the TreeTime pipeline (https://www.treetime.com) implemented in the Nextstrain analysis and presented in the phylogenetic tree (22).
Statistical Analysis and Ethics
We calculated 95% CIs by using published procedures (23,24). This study was conducted with due regard for European Union Principles and the Polish Law on Animal Protection. No permit from the Local Bioethical Committee for Animal Experimentation was obtained because animals were culled by the owner for production of pelts. Samples were collected postmortem.
Prevalence of SARS-CoV-2
We confirmed that 15 mink (16.5%, 95% CI 8.4%–28.6%) were positive for SARS-CoV-2. We summarized and provide the diagnostic results (Table 1).
SARS-CoV-2 Antigen Detection in Mink
Samples mink_4, mink_5, mink_48, mink_50, mink_77, and mink_83 had highly visible signals in both the control and test lines. Samples mink_20, mink_36, mink_42, mink_46, mink_49, mink_67, mink_76, and mink_88 had a highly visible control line and a much less pronounced test line. In all other samples, only the test line was visible. All 8 real-time RT-PCR‒positive samples were also positive by the antigen test. In addition, 5 E gene‒positive samples were also positive in the antigen test, but 4 were negative. The sample from mink_20 was positive in the antigen test, but SARS-CoV-2 RNA was not detected by RT-PCR.
Read Classification and SARS-CoV-2 Genome Sequences
The final validation of detection of SARS-CoV-2 in the mink was classification of the next-generation sequencing reads on the basis of the database containing reference viral, human, and American mink genomes. We used 3 independent approaches to obtain full viral genomes (Table 2). Only samples that had a complete SARS-CoV-2 genome sequence are shown. The number of Illumina reads generated for samples mink_4, mink_42, mink_49, mink_76, and mink_88 was not enough to produce full SARS-CoV-2 genomes. For these samples, the genomes were obtained by using the ARTIC procedure.
Phylogenetic Analysis of SARS-CoV-2 in Farmed Mink
We checked the 12 mink-originated SARS-CoV-2 sequences for mink-specific mutations detected earlier in mink from the Netherlands and Denmark, but found none. This finding suggested recent and separate introduction of SARS-CoV2 into mink from Poland (Figure 2). Alignment of full-genome sequences from 12 samples showed multiple polymorphisms at different nucleotide sites. Many of these polymorphisms gave rise to changes in amino acids when compared with the reference Wuhan-Hu-1 sequence (GenBank accession no. MN908947). Two specific mutations present in all samples were found in the spike protein: G75V and C1247F. The G75V mutation is present in 199 isolates published in GISAID (https://www.gisaid.org), and the C1257F mutation in 83 isolates. Other rare amino acid variants present in every SARS-CoV-2 isolate from mink in Poland were found in 5 additional proteins: nonstructural protein (nsp) 2, nsp3, nsp14, nsp15, and nucleocapsid protein.
On the basis of the dataset, we inferred a phylogenetic relationship by estimating the divergence times between each isolate (Figure 2). The analysis estimated that the most recent common ancestor for SARS-CoV-2 from mink in Poland and the 2 most similar sequences (German/NW-HHU-340/2020 and Norway/4235/2020) diverged on approximately September 31, 2020. We recognized mutations in amino acid sequences. If the molecular evolution started after virus introduction into the farm, this incident is estimated to have occurred on approximately October 4, 2020. Complete genome sequences of SARS-CoV-2 isolated from farmed mink in Poland have been deposited in GISAID (accession nos. EPI_ISL_732948‒59).
Identifying new species that can serve as animal sources of SARS-CoV-2 and predicting where novel outbreaks are most likely to occur are crucial steps for preventing and minimizing the extent of SARS-CoV-2 infections among humans (25). Recent reports confirmed the presence of SARS-CoV-2 in different animal species, including fur animals (i.e., mink and racoon dogs) (26,27). We report a 16.5% prevalence of SARS-CoV-2 in mink tested at a fur farm in northern Poland, confirming the presence of SARS-CoV-2 in farmed mink in Poland.
During our study, we used 2 different sequencing technologies to sequence the SARS-CoV-2 genomes. We found that amplicon-based nanopore sequencing gave better results than the bed-based enrichment Illumina approach. Conversely, Illumina reads showed a broader context because we were able to classify background reads that do not map to the SARS-CoV-2 genome as being of host origin. Therefore, these reads can be used as proof of sample origin. We also showed that the RdRp target for the quantitative PCR is less effective than the E gene in our experiment settings. The full genome of SARS-CoV-2 was assembled when both target genes were detected or the E gene was detected by a pair of positive signals in the antigenic assay.
Poland is one of the largest fur producers in Europe. Considering the number of farmed mink in the country and the large number of persons employed in this sector, we seek to increase awareness in the scientific community and mink industry that mink are susceptible to SARS-CoV-2 infection. Previous studies reported viral RNA detection in airborne inhalable dust in mink farms (8). Moreover, close contact of farmworkers with animals during feeding, culling, and dehiding increases their risk for exposure. We believe that a country-scale biomonitoring program should be activated as soon as possible to prevent the fur production sector from being a reservoir for future spillover of SARS-CoV-2 to humans. Samples for molecular diagnostics should be obtained from all farms in Poland following the highest standards for material collection, sample handling, and molecular detection of SARS-CoV-2.
We report a possible new genotype of SARS-CoV-2 that has sporadic mutations throughout the full genome sequence. Two mutations located in the spike protein (G75V and C1247F) were present in all isolates reported in this study. The G75V mutation is localized in the N terminal domain and could be responsible for interactions with host receptors or stabilizing the spike protein in a constrained prefusion state (28). To date, no other SARS-CoV-2 sequences deposited in GISAID have these 2 mutations simultaneously (29). We have recently detected possible zoonotic spillover of SARS-CoV-2 in worker employed at the farm described in this study (L. Rabalski et al., unpub. data). Preliminary genome analysis showed that the newly described isolates carry the combination of mutations typical of viruses isolated in November 2020, but additional new changes have accumulated since that time. We believe that wide monitoring of humans living near the mink farm should be performed to search for possible spillover and presence of new virus variants. Constant epizootiology monitoring is a crucial step in preventing new outbreaks of zoonotic diseases.
Dr. Rabalski is a virologist and assistant professor in the Laboratory of Recombinant Vaccines, Intercollegiate Faculty of Biotechnology of University of Gdansk and Medical University of Gdansk, Gdansk, Poland. His research interests include virology, epidemiology of viral pathogens, virus genomics, novel tools in molecular virology, and next-generation sequencing.
We thank the originating laboratories for obtaining the specimens; the submitting laboratories for generating genetic sequence data and sharing it through the GISAID Initiative; Alicja Rost, Ewa Zieliniewicz, and Karolina Baranowicz for laboratory assistance; Bartosz Wasąg for assistance with sequencing; the veterinary surgeons for assistance with sample collection; and the University of Gdańsk, the Medical University of Gdańsk, and the University of Helsinki for their support.
This research was funded through the 2018–2019 BiodivERsA joint call for research proposals, under the BiodivERsA3 ERA-Net COFUND program; the funding organizations ANR (France), DFG (Germany), EPA (Ireland), FWO (Belgium), and NCN (Poland); the Versatile Emerging Infections European Union Horizon 2020 (grant no. 874735); and the Jane and Aatos Erkko Foundation. M.G. was supported by the National Science Centre, Poland, under the BiodivERsA3 program (2019/31/Z/NZ8/04028), and the Young Scientists Program at the Medical University of Gdansk (664/772/61/71-1413). L.R. was supported by the Ministry of Science and Higher Education (decision no. 54/WFS 2020: Co-infections with SARS-CoV-2, database of COVID-19 accompanying infections).
- Zhou P, Yang X-L, Wang X-G, Hu B, Zhang L, Zhang W, et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature. 2020;579:270–3.
- Jo WK, de Oliveira-Filho EF, Rasche A, Greenwood AD, Osterrieder K, Drexler JF. Potential zoonotic sources of SARS-CoV-2 infections. Transbound Emerg Dis. 2020;10:1111.
- Mallapaty S. The search for animals harbouring coronavirus - and why it matters. Nature. 2021;591:26–8.
- Alexander MR, Schoeder CT, Brown JA, Smart CD, Moth C, Wikswo JP, et al. Predicting susceptibility to SARS-CoV-2 infection based on structural differences in ACE2 across species. FASEB J. 2020;34:15946–60.
- Damas J, Hughes GM, Keough KC, Painter CA, Persky NS, Corbo M, et al. Broad host range of SARS-CoV-2 predicted by comparative and structural analysis of ACE2 in vertebrates. Proc Natl Acad Sci U S A. 2020;117:22311–22.
- Hammer AS, Quaade ML, Rasmussen TB, Fonager J, Rasmussen M, Mundbjerg K, et al. SARS-CoV-2 transmission between mink (Neovison vison) and humans, Denmark. Emerg Infect Dis. 2021;27:547–51.
- Koopmans M. SARS-CoV-2 and the human-animal interface: outbreaks on mink farms. Lancet Infect Dis. 2021;21:18–9.
- Oreshkova N, Molenaar RJ, Vreman S, Harders F, Oude Munnink BB, Hakze-van der Honing RW, et al. SARS-CoV-2 infection in farmed minks, the Netherlands, April and May 2020. Euro Surveill. 2020;25:
- Oude Munnink BB, Sikkema RS, Nieuwenhuijse DF, Molenaar RJ, Munger E, Molenkamp R, et al. Transmission of SARS-CoV-2 on mink farms between humans and mink and back to humans. Science. 2021;371:172–7.
- World Organisation for Animal Health. COVID-19, 2021 [cited 2021 Jan 12]. https://www.oie.int/en/scientific-expertise/specific-information-and-recommendations/questions-and-answers-on-2019novel-coronavirus/events-in-animals
- World Health Organization. WHO-convened global study of the origins of SARS-CoV-2, 2020 [cited 2020 Dec 13]. https://www.who.int/publications/m/item/who-convened-global-study-of-the-origins-of-sars-cov-2
- Boklund A, Hammer AS, Quaade ML, Rasmussen TB, Lohse L, Strandbygaard B, et al. SARS-CoV-2 in Danish mink farms: course of the epidemic and a descriptive analysis of the outbreaks in 2020. Animals (Basel). 2021;11:164.
- European Union Fur Association. Sustainable fur, 2020 [cited 2020 Dec 13]. https://www.sustainablefur.com
- Związku Przedsiębiorców i Pracodawców. Facts about breeding fur animals in Poland [in Polish]. Warsaw; 2020 [cited 2021 Jun 17]. https://zpp.net.pl/wp-content/uploads/2020/09/Raport-ZPP_Fakty-o-hodowli-zwierząt-futerkowych-w-Polsce.pdf
- European Centre for Disease Prevention and Control. COVID-19 situation update for the EU/EEA and the UK, 2021 [cited 2021 May 12]. https://www.ecdc.europa.eu/en/cases-2019-ncov-eueea
- Corman VM, Landt O, Kaiser M, Molenkamp R, Meijer A, Chu DK, et al. Detection of 2019 novel coronavirus (2019-nCoV) by real-time RT-PCR. Euro Surveill. 2020;25:
- Wood DE, Lu J, Langmead B. Improved metagenomic analysis with Kraken 2. Genome Biol. 2019;20:257.
- Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.
- Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM, 2013 [cited 2021 May 13]. https://github.com/lh3/bwa
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al.; 1000 Genome Project Data Processing Subgroup. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078–9.
- Hadfield J, Megill C, Bell SM, Huddleston J, Potter B, Callender C, et al. Nextstrain: real-time tracking of pathogen evolution. Bioinformatics. 2018;34:4121–3.
- Sagulenko P, Puller V, Neher RA. TreeTime: Maximum-likelihood phylodynamic analysis. Virus Evol. 2018;4:
- Sokal RR, Rohlf FJ. Biometry: the principles and practices of statistics in biological research. New York: W. H. Freeman; 1995.
- Grzybek M, Bajer A, Bednarska M, Al-Sarraf M, Behnke-Borowczyk J, Harris PD, et al. Long-term spatiotemporal stability and dynamic changes in helminth infracommunities of bank voles (Myodes glareolus) in NE Poland. Parasitology. 2015;142:1722–43.
- Hemida MG, Ba Abduallah MM. The SARS-CoV-2 outbreak from a one health perspective. One Health. 2020;10:
- Abdel-Moneim AS, Abdelwhab EM. Evidence for SARS-CoV-2 infection of animal hosts. Pathogens. 2020;9:529.
- Freuling CM, Breithaupt A, Müller T, Sehl J, Balkema-Buschmann A, Rissmann M, et al. Susceptibility of raccoon dogs for experimental SARS-CoV-2 infection. Emerg Infect Dis. 2020;26:2982–5.
- Arya R, Kumari S, Pandey B, Mistry H, Bihani SC, Das A, et al. Structural insights into SARS-CoV-2 proteins. J Mol Biol. 2021;433:
- Shu Y, McCauley J. GISAID: Global initiative on sharing all influenza data - from vision to reality. Euro Surveill. 2017;22:30494.
TablesCite This Article
Original Publication Date: August 09, 2021