Transmission of SARS-CoV-2 on mink farms between humans and mink and back to humans

Two-way transmission on mink farms Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a zoonotic virus—one that spilled over from another species to infect and transmit among humans. We know that humans can infect other animals with SARS-CoV-2, such as domestic cats and even tigers in zoos. Oude Munnink et al. used whole-genome sequencing to show that SARS-CoV-2 infections were rife among mink farms in the southeastern Netherlands, all of which are destined to be closed by March 2021 (see the Perspective by Zhou and Shi). Toward the end of June 2020, 68% of mink farm workers tested positive for the virus or had antibodies to SARS-CoV-2. These large clusters of infection were initiated by human COVID-19 cases with viruses that bear the D614G mutation. Sequencing has subsequently shown that mink-to-human transmission also occurred. More work must be done to understand whether there is a risk that mustelids may become a reservoir for SARS-CoV-2. Science, this issue p. 172; see also p. 120

I n late December 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was identified as the cause of a viral pneumonia outbreak, possibly related to a seafood and live animal market in Wuhan, China (1). Since then, SARS-CoV-2 has spread across the world. By 8 October 2020, more than 36.1 million people had been infected with SARS-CoV-2, resulting in more than 1 million deaths (2). In the Netherlands, more than 155,000 infections have been confirmed, more than 6500 SARS-CoV-2-related deaths have been reported, and nonpharmaceutical interventions have been put into place to prevent further spread of SARS-CoV-2 (3).
In view of the similarities of the new virus and SARS-CoV-1, a zoonotic origin of the outbreak was suspected and linked to a Wuhan fresh market where various animals-including fish, shellfish, poultry, wild birds, and exotic animals-were sold. However, other cases with onset well before the period correlated with the Wuhan market-associated cluster were observed, which suggests the possibility of other sources (4). Although closely related coronaviruses found in bats (5,6) and pangolins (7,8) have the greatest sequence identity to SARS-CoV-2, the most likely divergence of SARS-CoV-2 from the most closely related bat sequence is estimated to have occurred between 1948 and 1982 (9). Therefore, the animal reservoir(s) of SARS-CoV-2 is(are) yet to be identified.
In response to the outbreaks in mink farms, the Dutch national response system for zoonotic diseases was activated. Although the public health risk of exposure to animals with SARS-CoV-2 was determined to be low, increased awareness of animals' possible involvement in the COVID-19 epidemic was needed. Therefore, from 20 May 2020 onward, mink farmers, veterinarians, and laboratories were obliged to report symptoms in mink (family Mustelidae) to the Netherlands Food and Consumer Product Safety Authority (NFCPSA), and an extensive surveillance system was established (35).
Whole-genome sequencing (WGS) can be used to monitor the emergence and spread of pathogens (36)(37)(38)(39). As part of the surveillance effort in the Netherlands, more than 1750 SARS-CoV-2 viruses have been sequenced to date from patients from different parts of the country (40). Here, we describe an indepth investigation into the SARS-CoV-2 outbreak in mink farms and mink farm employees in the Netherlands, combining epidemiological information, surveillance data, and WGS on the human-animal interface.
SARS-CoV-2 was first diagnosed on two mink farms (designated NB1 and NB2) in the Netherlands on 23 and 25 April 2020, respectively. After the initial detection of SARS-CoV-2 on these farms, a thorough investigation was initiated to identify potential transmission routes and to perform an environmental and occupational risk assessment. Here, we describe the results of the investigation of the first 16 SARS-CoV-2-affected mink farms by combining SARS-CoV-2 diagnostics, WGS, and in-depth interviews.
Owners and employees of the 16 mink farms with SARS-CoV-2-positive animals were included in the contact tracing investigation by the Netherlands Municipal Health Services and were tested according to national protocol. Ninety-seven individuals were tested by either serological assays and/or RT-PCR. Fortythree of 88 (49%) upper respiratory tract samples tested positive by RT-PCR, whereas 38 of 75 (51%) serum samples tested positive for SARS-CoV-2-specific antibodies. In total, 66 of 97 (68%) tested individuals had evidence for SARS-CoV-2 infection (Table 1).
During the interview on 28 April, four of five employees from NB1 reported that they had experienced respiratory symptoms before the outbreak was detected in minks, but none of them had been tested for SARS-CoV-2. Their symptom-onset dates ranged from 1 April to 9 May. For 16 of the mink (sampled on 28 April) and one farm employee (sampled on 4 May), a whole-genome sequence was obtained (hCov-19/Netherlands/NoordBrabant_177/2020). The human sequence clusters within the mink sequences, although it differs from the closest mink sequence by 7 nucleotides (nts) ( Fig. 1 and cluster A in Figs. 2 and 3). On farm NB2, SARS-CoV-2 was diagnosed on 25 April. Retrospective analysis showed that one employee from NB2 had been hospitalized with SARS-CoV-2 on 31 March. All samples from the eight employees taken on 30 April were negative by RT-PCR but tested positive for SARS-CoV-2 antibodies. The virus sequence obtained from NB2 animals was distinct from that of NB1 animals, indicating a separate introduction (Figs. 2 and 3, cluster B).
On mink farm NB3, SARS-CoV-2 infection was diagnosed on 7 May. Initially, all seven employees tested negative for SARS-CoV-2, but when retested between 19 and 26 May after developing COVID-19-related symptoms, five of seven individuals working or living on the farm tested positive for SARS-CoV-2 RNA. Whole-genome sequences were obtained from these five individuals. The clustering of these sequences with the sequences derived from NB3 mink, together with initial negative test results and subsequent symptom onset, indicate that the employees were infected with SARS-CoV-2 after mink on the farm became infected. An additional infection was identified from contact tracing: An individual who did not visit the farm but had close contact with one of the employees became infected with the SARS-CoV-2 strain found on farm NB3. Animal and human sequences from farm NB3 were close to those from farm NB1, and both were categorized in cluster A.
Similarly, on mink farm NB7, zoonotic transmission from mink to humans likely occurred. On this farm, SARS-CoV-2 infection in mink was diagnosed on 31 May. NB7 employees initially tested negative for SARS-CoV-2, but several of them subsequently began to show symptoms. Between 10 June and 1 July, samples were taken from 10 employees, 8 of which tested positive for SARS-CoV-2 RNA. WGS of two NB7 employee samples showed that their virus sequences clustered with the sequences from the mink at this farm.
The sequences generated from animals and employees of mink farms were compared with~1775 whole-genome sequences in the national database. To discriminate between community-acquired and mink farm-           infections and to determine the potential risk for people who live near mink farms, WGS was also performed on 34 SARS-CoV-2-positive samples (collected from 4 March until 29 April 2020) from individuals who live in the same geographic area (as determined by four-digit postal code) in which farms NB1 to NB4 are located. These local sequences, sampled in a proxy of 19 km 2 , reflected the general diversity of SARS-CoV-2 seen in the Netherlands and were not related to the clusters of mink sequences found on the mink farms, thereby indicating that no spillover to people living in close proximity to mink farms had occurred and that the sequences from SARS-CoV-2-infected animals and farm workers clustered by farm ( Fig. 2; sequences from community shown in magenta). The sequences from the mink farm investi-gation were also compared with sequences from Poland (n = 65) because many of the mink farm workers were seasonal migrants from Poland, but the Polish sequences were more divergent.
Phylogenetic analysis of the mink SARS-CoV-2 genomes showed that mink sequences of 16 farms were grouped into five different clusters (Figs. 2 and 3). Viruses from farms NB1, NB3, NB4, NB8, NB12, NB13, and NB16 belonged to cluster A; sequences from NB2 formed a distinct cluster (cluster B); those from farms NB6, NB7, NB9, and NB14 formed cluster C; those from NB5, NB8, NB10, and NB15 formed cluster D; and those from NB11 were designated as cluster E. On farm NB8, SARS-CoV-2 viruses from clusters A and D were found. A detailed inventory of possible common characteristics-including farm owner, shared personnel, feed supplier, and veterinary service provider-was made. Multiple farms within a cluster shared the same owner; however, in most cases no common factor could be identified for most farms, and clustering could not be explained by geographic distance (Table 2 and Fig. 4).
In total, 18 sequences from mink farm employees or their close contacts were generated from seven different farms. In most cases, these human sequences were nearly identical to the mink sequences from the same farm. For NB1, the situation was different: The human sequence clusters deeply within the sequences derived from mink (Fig. 1), differing from the closest related mink sequence by 7 nts. This was also the case on farm NB14, with a 4-nt difference from the closest related mink sequence. Sequences of employees at mink farm NB8 clustered with those of animals from NB12, likely because personnel were exchanged between these two farms. SARS-CoV-2 was detected on mink farms NB1 and NB4 after reports of respiratory symptoms and increased mortality in mink. The sequences from farm NB1 showed a difference of 0 to 9 single-nucleotide polymorphisms (SNPs) (average of 3.9 nts) and sequences from NB4 had a difference of 0 to 8 SNPs (average of 3.6 nts), which is more than is generally observed in outbreaks in human settings. In addition, two deletions (one of 12 nts and one of 134 nts) were observed in a single sequence from NB1. The sequences of mink at NB6 had differences of 0 to 12 SNPs, and a deletion of 9 nts was observed in one sequence, whereas diversity was lower for the subsequent farm sequences ( Table 2). After the initial detection of SARS-CoV-2, farms were screened weekly. The first, second, fifth, and sixth weekly screening yielded new positive cases.
Several nonsynonymous mutations were identified among the mink sequences when compared with the Wuhan reference sequence NC_045512.2. However, no particular amino acid substitutions were found in all-mink samples ( fig. S1). Of note, three clusters had the position 614G variant (clusters A, C, and E), and two had the original variant. On the basis of the data available at this stage, there were no obvious disease presentation differences in animals or humans between clusters, but data collection and analysis are ongoing for cases after NB16. The mutations we observed can also be found in the general human population, and the same mutations were found in human cases related to the mink farms.
Our work shows evidence of ongoing SARS-CoV-2 transmission in mink farms and spillover events to humans. More research in minks and other mustelid species will be needed to determine whether these species are at risk of becoming reservoirs of SARS-CoV-2. After the detection of SARS-CoV-2 on mink farms, 68% of the tested farm workers and/or relatives or contacts would later become or had been infected with SARS-CoV-2, indicating that contact with SARS-CoV-2-infected mink is a risk factor for contracting COVID-19. Recently, an eightfold increase in cytidine-to-uridine (C→U) substitutions compared with uridineto-cytidine (U→C) substitutions was described, suggestive of host adaptation (41). In the mink sequences, we observed a 3.5-fold increase in C→U compared with U→C substitutions, but the number of substitutions (185) was limited.
A high diversity was observed in the sequences from some mink farms, which is likely explained by multiple generations of viral infections in animals before the increase in mortality was detected. Current estimates         indicate that the substitution rate of SARS-CoV-2 is around 1.16 × 10 −3 substitutions per site per year in the human population (42), which corresponds to around one mutation per 2 weeks. This could mean that the virus was already circulating in mink farms for some time before it was identified. However, a relatively high sequence diversity was also observed in farmed mink that tested negative week before a positive test, hinting toward a faster evolutionary rate of the virus in the mink population. Mink farms have large populations of animals, living at high density, which may promote virus transmission. However, the moment of viral introduction was not known, which makes it difficult to draw definite conclusions about the substitution rate in mink farms. Our sequencing did not reveal any systematic mutations that would need to be assessed for potential phenotypic effects. The generation interval for SARS-CoV-2 in humans has been estimated to be around 4 to 5 days (43), but with high-dose exposure in a farm with a high number and density of animals, this interval could potentially be shorter. Further evidence that animals were the most likely source of human infection was provided by the clear phylogenetic separation between mink farm-related human and animal sequences and sequences from human cases within the same geographic area (as determined by four-digit postal code). However, some of the farm-related humans may have been infected within their household, not directly from mink. Spillback into the local community was not observed in our sequence data.
Thus far, the investigation has failed to identify common factors that might explain farm-to-farm spread, possibly via temporary workers who were not included in testing. Since our observations were made, SARS-CoV-2 infections have also been described in mink farms elsewhere (44)(45)(46). Additional research efforts are needed, as it is imperative that the fur production and trading sector not become a reservoir for future spillover of SARS-CoV-2 to humans. code, and materials used in this work are publicly available. All sequences are publicly available in the GISAID database under accession numbers EPI_ISL_461190 to EPI_ISL_461192, EPI_ISL_461200, EPI_ISL_461202, EPI_ISL_461203, EPI_ISL_522987 to EPI_ISL_523034, EPI_ISL_523040, EPI_ISL_523046, EPI_ISL_523068, EPI_ISL_523070, EPI_ISL_523072, EPI_ISL_523073, EPI_ISL_523075, EPI_ISL_523085, EPI_ISL_523089 to EPI_ISL_523120, EPI_ISL_523282 to EPI_ISL_523286, EPI_ISL_523301, EPI_ISL_523310 to EPI_ISL_523312, EPI_ISL_523333, EPI_ISL_523493, and EPI_ISL_523494. Ethical approval was not required for this study, as anonymous aggregated data were used, and no medical interventions were made on animal or human subjects. This work is licensed under a Creative Commons Attribution 4.0 International (CC BY 4.0) license, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. To view a copy of this license, visit https://creativecommons.org/licenses/by/4.0/. This license does not apply to figures/photos/artwork or other content included in the article that is credited to a third party; obtain authorization from the rights holder before using such material.