Susceptibility of Domestic Goat (Capra aegagrus hircus) to Experimental Infection with Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) B.1.351/Beta Variant

A wide range of animal species are susceptible to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Natural and/or experimental infections have been reported in pet, zoo, farmed and wild animals. Interestingly, some SARS-CoV-2 variants, such as B.1.1.7/Alpha, B.1.351/Beta, and B.1.1.529/Omicron, were demonstrated to infect some animal species not susceptible to classical viral variants. The present study aimed to elucidate if goats (Capra aegagrus hircus) are susceptible to the B.1.351/Beta variant. First, an in silico approach was used to predict the affinity between the receptor-binding domain of the spike protein of SARS-CoV-2 B.1.351/Beta variant and angiotensin-converting enzyme 2 from goats. Moreover, we performed an experimental inoculation with this variant in domestic goat and showed evidence of infection. SARS-CoV-2 was detected in nasal swabs and tissues by RT-qPCR and/or immunohistochemistry, and seroneutralisation was confirmed via ELISA and live virus neutralisation assays. However, the viral amount and tissue distribution suggest a low susceptibility of goats to the B.1.351/Beta variant. Therefore, although monitoring livestock is advisable, it is unlikely that goats play a role as SARS-CoV-2 reservoir species, and they are not useful surrogates to study SARS-CoV-2 infection in farmed animals.


Introduction
The causative agent of the coronavirus infectious disease 2019 (COVID-19) pandemic, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is mainly transmitted via respiratory droplets or aerosols from human to human [1]. Additionally, SARS-CoV-2 can infect several animal species that may act as intermediate hosts and play a role in the transmission and maintenance of the virus. Considering that the origin of this pandemic is

In Silico Study
We ran MODELLER v10.1 [35] to generate ten different homology models of goat ACE2 in complex with the spike RBD of the SARS-CoV-2 ancestral variant. The crystallographic structure of human ACE2 in complex with the spike RBD of the ancestral variant (PDB ID 6M0J) was used as a template. Then, using the MODELLER models as input, we used FoldX v5 [36] to predict the changes in binding affinity induced by the mutations of the variants. We considered the following mutations located at the RBD: K417N, E484K, N501Y for B.1.351/Beta, L452R, T478K for B.1.617.2/Delta, K417T, E484K, N501Y for P.1/Gamma, and G339D, S371L, S373P, S375F, K417N, N440K, G446S, S477N, T478K, E484A, Q493R, G496S, Q498R, N501Y, Y505H for B.1.1.529/Omicron (BA.1) variants. We ran FoldX modules repairPDB followed by BuildModel. We fixed all FoldX parameters to their default values except parameter vdwDesign, which was set to zero. As a positive control, we performed the same steps using the human ACE2 sequence to predict the changes in binding affinity caused by the mutations of the SARS-CoV-2 VOC in humans.
The infectious titre of the SARS-CoV-2 stock was calculated by determining the dilution that caused cytopathic effect (CPE) in 50% of the inoculated Vero E6 cells (50% tissue culture infectious dose endpoint, TCID 50 ).

Experimental Study Design
This experimental study was approved by the Institutional Animal Welfare Committee of the Institute of Agrifood Research and Technology (CEEA-IRTA) and by the Ethical Commission of Animal Experimentation of the Autonomous Government of Catalonia (CEA-OH/11586/3) and conducted by certified staff. Experiments involving SARS-CoV-2 were performed at the Biosafety Level-3 (BSL-3) facilities of the Biocontainment Unit of IRTA-CReSA (Barcelona, Spain).
A total of 18 goats of 2-3 months of age were acquired from a Spanish commercial farm (La botiga d'Ullastrell S.L, Barcelona, Spain) and included in the study. Three animals were used as non-inoculated controls and were necropsied on the day of animal arrival. The remaining 15 animals were allocated inside an experimental box of the animal BSL-3 facilities, and after 7 days of acclimation, they were intranasally (IN) inoculated with 2 mL of 10 6 TCID 50 /mL of SARS-CoV-2 B.1.351/Beta variant (1 mL in each nostril) using a diffusor device (LMA™ MAD ® Nasal, Teleflex LLC). At 2, 4, 7, 10 and 18 days post-inoculation (dpi), three goats/day were euthanised and complete necropsies were performed (Figure 1). At necropsy, nasal and rectal swabs (NS/RS), blood, and the following tissues were collected: nasal turbinate (NT; caudal and cranial portions), olfactory bulb, parotid gland, trachea, pharynx, tonsil, lung (apical, medial, and caudal portions), lymph nodes (LN; cervical, mesenteric, and mediastinal ones), kidney, spleen, jejunum, colon, and heart. Clinical signs and rectal temperatures were recorded daily during the whole study. Figure 1. Schematic representation of the study design. Three out of eighteen goats were necropsied on the day of arrival and served as non-inoculated controls. Nasal and rectal swabs, blood, and tissue samples were collected. After acclimatisation, the remaining fifteen goats were intranasally inoculated with the SARS-CoV-2 B.1.351/Beta variant (2 mL of 1 × 10 6 TCID50/animal). At 2, 4, 7,10, and 18 days post-infection, three goats were necropsied each day and nasal and rectal swabs, blood, and tissues were collected. Clinical signs and rectal temperature were recorded daily.

Virus Detection and Isolation
Tissues, NS, and RS were collected during necropsy and were transferred into cryotubes containing 500 μL DMEM supplemented with 100 U/mL penicillin, 100 μg/mL streptomycin, and 2 mM glutamine, vortexed and stored at −75 °C until use. Blood was centrifuged (1800× g, 10 min at 4 °C) and sera were collected and kept at −20 °C until further use. Detection of SARS-CoV-2 RNA was performed via genomic and subgenomic reverse transcription-quantitative PCR (RT-qPCR) in all collected samples, as previously described [38,39]. Briefly, viral RNA was extracted using the Indimag Pathogen Kit (Indical Biosciences, Leipzig, Germany) on a Biosprint 96 workstation (Qiagen, Hilden, Germany) in accordance with the manufacturers' instructions. Then, RT-qPCR was carried out using AgPath-IDTM One-Step RT-PCR Reagents (Applied Biosystems, Life Technologies), and amplification was performed using a 7500 Fast Real-Time PCR System (Applied Biosystems, Life Technologies, Waltham, MA, USA). Samples with a Ct value < 40 were considered positive for SARS-CoV-2 RNA.
RT-qPCR positive samples with Ct < 30 were evaluated for the presence of infectious virus by titration in Vero E6 cells, as previously reported [37]. Briefly, ten-fold dilutions from each sample were transferred to Vero E6 monolayers and incubated at 37 °C and 5% CO2. Plates were monitored under a light microscope and wells were assessed for the presence of CPE for 6 days. The amount of infectious virus was calculated by determining the TCID50 using the Reed and Muench method [40].

Neutralising Antibodies
Blood samples of the 18 experimental goats were centrifuged at 1800 × g for 10 min at 4 °C, and then the obtained sera were inactivated at 56 °C for 1 h. Neutralising antibodies targeting SARS-CoV-2 RBD from all samples were measured with the GenScript cPass TM SARS-CoV-2 Neutralization Antibody Detection Kit (Genscript, the Netherlands), following the manufacturer's protocol. The percentage of inhibition of each sample was calculated by the following formula: % Inhibition = (1 − (OD450 sample/OD450 negative control)) × 100. Samples and controls were included in duplicate (SD ≤ 10%). Inhibition of ≥ 30% was considered a positive neutralisation.
In addition, live virus neutralisation assay was performed as previously described [37]. Pre-inactivated serum samples were initially diluted 1:20 and then 2-fold serially diluted in supplemented DMEM, mixed 1:1 with the SARS-CoV-2 isolate Figure 1. Schematic representation of the study design. Three out of eighteen goats were necropsied on the day of arrival and served as non-inoculated controls. Nasal and rectal swabs, blood, and tissue samples were collected. After acclimatisation, the remaining fifteen goats were intranasally inoculated with the SARS-CoV-2 B.1.351/Beta variant (2 mL of 1 × 10 6 TCID 50 /animal). At 2, 4, 7,10, and 18 days post-infection, three goats were necropsied each day and nasal and rectal swabs, blood, and tissues were collected. Clinical signs and rectal temperature were recorded daily.

Virus Detection and Isolation
Tissues, NS, and RS were collected during necropsy and were transferred into cryotubes containing 500 µL DMEM supplemented with 100 U/mL penicillin, 100 µg/mL streptomycin, and 2 mM glutamine, vortexed and stored at −75 • C until use. Blood was centrifuged (1800× g, 10 min at 4 • C) and sera were collected and kept at −20 • C until further use. Detection of SARS-CoV-2 RNA was performed via genomic and subgenomic reverse transcription-quantitative PCR (RT-qPCR) in all collected samples, as previously described [38,39]. Briefly, viral RNA was extracted using the Indimag Pathogen Kit (Indical Biosciences, Leipzig, Germany) on a Biosprint 96 workstation (Qiagen, Hilden, Germany) in accordance with the manufacturers' instructions. Then, RT-qPCR was carried out using AgPath-IDTM One-Step RT-PCR Reagents (Applied Biosystems, Life Technologies), and amplification was performed using a 7500 Fast Real-Time PCR System (Applied Biosystems, Life Technologies, Waltham, MA, USA). Samples with a Ct value < 40 were considered positive for SARS-CoV-2 RNA.
RT-qPCR positive samples with Ct < 30 were evaluated for the presence of infectious virus by titration in Vero E6 cells, as previously reported [37]. Briefly, ten-fold dilutions from each sample were transferred to Vero E6 monolayers and incubated at 37 • C and 5% CO 2 . Plates were monitored under a light microscope and wells were assessed for the presence of CPE for 6 days. The amount of infectious virus was calculated by determining the TCID 50 using the Reed and Muench method [40].

Neutralising Antibodies
Blood samples of the 18 experimental goats were centrifuged at 1800× g for 10 min at 4 • C, and then the obtained sera were inactivated at 56 • C for 1 h. Neutralising antibodies targeting SARS-CoV-2 RBD from all samples were measured with the GenScript cPass TM SARS-CoV-2 Neutralization Antibody Detection Kit (Genscript, the Netherlands), following the manufacturer's protocol. The percentage of inhibition of each sample was calculated by the following formula: % Inhibition = (1 − (OD 450 sample/OD 450 negative control)) × 100. Samples and controls were included in duplicate (SD ≤ 10%). Inhibition of ≥ 30% was considered a positive neutralisation.
In addition, live virus neutralisation assay was performed as previously described [37]. Pre-inactivated serum samples were initially diluted 1:20 and then 2-fold serially diluted in supplemented DMEM, mixed 1:1 with the SARS-CoV-2 isolate (EPI_ISL_3164134), and incubated for 1 h at 37 • C. Then, each dilution mixture (in duplicates) was transferred to Vero E6 monolayers containing 100 TCID 50 of SARS-CoV-2 per well and cultured for 3 days at 37 • C and 5% CO 2 . Then, CPE was measured using the CellTiter-Glo luminescent cell viability assay (Promega, Madison, WI, USA) in accordance with the manufacturer's protocol. Luminescence was measured as luminescence units in a Fluoroskan Ascent FL luminometer (ThermoFisher Scientific, Waltham, MA, USA). Serum virus neutralisation titre 50 (SNT 50 ) was defined as the reciprocal value of the sample dilution that showed 50% of the SARS-CoV-2-induced CPE in Vero E6 cells.

Histopathology
Upper (NT) respiratory tract and tonsil were harvested and fixed by immersion in 10% buffered formalin solution. Fixed tissues were embedded in paraffin, and cut sections (3 µm thick) were stained with haematoxylin and eosin (H/E) for histopathology studies using an optical microscope.
A previously described immunohistochemistry (IHC) technique to detect SARS-CoV-2 NP antigen [41] using a Rabbit monoclonal antibody (40143-R019, SinoBiological, Beijing, China) at dilution 1:10,000 was applied on NT and tonsil. The amount of viral antigen in tissue samples was semi-quantitatively scored: lack of antigen detection, low, moderate, and high amount following a previously published scoring system [37].

Statistical Analyses
Dose-response curves of neutralisation assay in serum samples were adjusted to a non-linear fit regression model calculated with a normalised logistic curve with variable slope. For data normalisation, uninfected cells and untreated virus-infected cells were used as negative and positive control of infection (% Neutralisation = (RLUmax − RLUexperimental)/(RLUmax − RLUmin) × 100), respectively. All SNT 50 values were determined as the concentration of sera blocking 50% of the SARS-CoV-2-induced CPE in Vero E6 cells and were expressed as reciprocal dilution. All statistical analyses were performed with GraphPad Prism 8.4.3 (GraphPad Software, Inc., San Diego, CA, USA).

529/Omicron Variants
We ran FoldX to predict the changes in binding affinity (∆∆G = ∆G variant − ∆G ancestral ) for the spike RBD in complex with goat ACE2, induced by the mutations of the variants we studied (see Figure 2). As a control, we also computed the binding affinity changes for the spike RBD in complex with human ACE2. Results for goat ACE2 indicated that Beta ∆∆G was significantly lower (higher binding affinity changes) than Delta, Gamma, and Omicron ∆∆Gs (Mann-Whitney-Wilcoxon test p-values equal to 1.06 × 10 −2 , 3.64 × 10 −3 , and 3.20 × 10 −2 , respectively). We found a similar trend for human ACE2 (Mann-Whitney-Wilcoxon test p-values equal to 6.58 × 10 −4 , 2.91 × 10 −4 , and 1.06 × 10 −2 for Delta, Gamma, and Omicron, respectively). Positive/negative ∆∆Gs values indicated if the binding affinity was lower/higher for the evaluated variant in comparison with the ancestral spike RBD in complex with goat ACE2. The average ∆∆Gs for the Beta variant was −2.18, with a 95% confidence interval (CI) of (−3.67, −0.69), for the Delta variant −0.06, with a 95% CI of (−0.58, 0.46), for the Gamma variant 0.25 with a 95% CI of (0.97, 1.48), and for the Omicron variant 0.20, with a 95% CI of (−1.86, 2.26). We did not find significant differences between the predicted ∆∆Gs for goat and human for any of the variants (Mann-Whitney-Wilcoxon test p-values of 4.85 × 10 −1 , 5.15 × 10 −1 , 7.64 × 10 −1 , and 6.61 × 10 −1 for Delta, Beta, Gamma, and Omicron, respectively).

Clinical Signs and Pathological Findings
No clinical signs or increase in rectal temperature from baseline were observed after challenge in any of the goats along the study. Moreover, no gross or microscopic lesions attributable to SARS-CoV-2 infection were recorded in any of the studied animals.

Virus Detection and Replicative Viral Isolation in Goat Samples
Fourteen out of the fifteen goats tested positive for SARS-CoV-2 by genomic RT-qPCR in NS at 2 dpi and ten out of twelve at 4 dpi with low viral loads (Ct ≥ 27). Only two out of the nine remaining goats were positive with low viral loads (Ct > 30) at 7 dpi ( Figure  3a). Moreover, low levels of viral RNA (Ct > 23) were observed in caudal and cranial NT (Figure 3b), lymphoid tissues (tonsil and cervical and/or mediastinal LN) (Figure 3c,d), and respiratory tract (trachea and lung) (Figure 3e,f). SARS-CoV-2 RNA from tonsil and LN samples were detected until 18 dpi in goats (Figure 3c,d). Regarding subgenomic RNA, tonsil samples were positive at 2, 4, and 7 dpi (Ct > 23) as well as cranial and caudal NT at 2 dpi (Ct > 25). Blood and all other samples (RS, olfactory bulb, parotid gland, trachea, spleen, kidney, jejunum, and heart) tested negative (Ct ≥ 40) via genomic and subgenomic RT-qPCR. FoldX predicted ∆∆G for B.1.351/Beta, B.1.617.2/Delta, P.1/Gamma, and B.1.1.529/Omicron variants for human and goat ACE2 (blue and orange boxes, respectively). For goat, the computed ∆∆G was significantly lower for the Beta variant than for the Delta, Gamma, and Omicron variants, with Mann-Whitney-Wilcoxon test p-values of 1.06 × 10 −2 , 3.64 × 10 −3 , and 3.20 × 10 −2 , respectively. The average ∆∆G predicted from the 10 Modeller models is plotted as a white dot for each variant. We found no significant differences between the predicted ∆∆G values for human and goat of the different variants (Mann-Whitney-Wilcoxon test p-values of 4.85 × 10 −1 , 5.15 × 10 −1 , 7.64 × 10 −1 , 6.61 × 10 −1 for Delta, Beta, Gamma, and Omicron, respectively). In this boxplot representation, the lower and upper ends of each box represent the first (Q1) and third (Q3) quartiles of the ∆∆G predicted values, respectively. The horizontal line, inside each box, represents the median, or second quartile (Q2), and the mean is plotted as a white dot for each variant. The box "whiskers" extend to values that are 1.5 times the size of the interquartile range (IQR = Q3 − Q1). Values that fall outside this range are displayed independently as black diamonds. Mann-Whitney-Wilcoxon test p-values are annotated according to the following criteria: ns (0.05 < p-value ≤ 1), * (0.02 < p-value ≤ 0.05), ** (0.001 < p-value ≤ 0.02).

Clinical Signs and Pathological Findings
No clinical signs or increase in rectal temperature from baseline were observed after challenge in any of the goats along the study. Moreover, no gross or microscopic lesions attributable to SARS-CoV-2 infection were recorded in any of the studied animals.

Virus Detection and Replicative Viral Isolation in Goat Samples
Fourteen out of the fifteen goats tested positive for SARS-CoV-2 by genomic RT-qPCR in NS at 2 dpi and ten out of twelve at 4 dpi with low viral loads (Ct ≥ 27). Only two out of the nine remaining goats were positive with low viral loads (Ct > 30) at 7 dpi (Figure 3a). Moreover, low levels of viral RNA (Ct > 23) were observed in caudal and cranial NT (Figure 3b), lymphoid tissues (tonsil and cervical and/or mediastinal LN) (Figure 3c,d), and respiratory tract (trachea and lung) (Figure 3e,f). SARS-CoV-2 RNA from tonsil and LN samples were detected until 18 dpi in goats (Figure 3c,d). Regarding subgenomic RNA, tonsil samples were positive at 2, 4, and 7 dpi (Ct > 23) as well as cranial and caudal NT at 2 dpi (Ct > 25). Blood and all other samples (RS, olfactory bulb, parotid gland, trachea, spleen, kidney, jejunum, and heart) tested negative (Ct ≥ 40) via genomic and subgenomic RT-qPCR. Only samples with a Ct < 30 by RT-qPCR were tested for virus titration on cell culture (Table S1). We isolated infectious virus on Vero E6 cells at 2 dpi from one tonsil (3.1 TCID50/mL) and one cranial NT (1.9 TCID50/mL) from two different goats. SARS-CoV-2 could not be isolated from all other tested samples.

Detection of SARS-CoV-2 Nucleoprotein in Tissues by Immunohistochemistry
Tonsils from inoculated goats collected at 2, 4, and 7 dpi (one goat per day) were positive for IHC staining ( Figure 4); labelling was mainly found in dendritic-like cells mostly located around tonsillar crypts. The amount of viral antigen was low in all animals. All other samples were negative. Only samples with a Ct < 30 by RT-qPCR were tested for virus titration on cell culture (Table S1). We isolated infectious virus on Vero E6 cells at 2 dpi from one tonsil (3.1 TCID 50 /mL) and one cranial NT (1.9 TCID 50 /mL) from two different goats. SARS-CoV-2 could not be isolated from all other tested samples.

Detection of SARS-CoV-2 Nucleoprotein in Tissues by Immunohistochemistry
Tonsils from inoculated goats collected at 2, 4, and 7 dpi (one goat per day) were positive for IHC staining ( Figure 4); labelling was mainly found in dendritic-like cells mostly located around tonsillar crypts. The amount of viral antigen was low in all animals. All other samples were negative.

Discussion
Goat is an important livestock animal that has already been suggested to be susceptible to SARS-CoV-2 by previous studies, both in vitro and in vivo [31,42]. Similarity between goat and human ACE2 receptors has been demonstrated via comparative genomic analysis as well as binding affinity between goat ACE2 and SARS-CoV-2 S protein [42]. Other previous in vitro studies also reported that goat ACE2 supports cell entry of SARS-CoV-2 [43,44], and its replication has been demonstrated in DF1 cells expressing both goat ACE2 and transmembrane serine protease 2 (TMPRSS2) [43]. Since the susceptibility of animal species could change after the appearance of SARS-CoV-2 variants and our in silico study predicted potential affinity between the B.1.351/Beta variant and ACE2 goat receptor, the aim of the present study was to investigate the susceptibility of goat to this variant of SARS-CoV-2 by means of an in vivo experimental inoculation.
Previously, Bosco-Lauth et al. [31] described SARS-CoV-2 RT-qPCR positive results in two out of three inoculated goats at 3 days after experimental challenge with the ancestral variant. However, none of the animals shed detectable infectious virus, nor did other ruminants such as sheep and cattle. The aforementioned study together with other previous studies reported that goat, sheep, or cattle developed limited subclinical infection after SARS-CoV-2 inoculation [20,32]. In sheep, Gaudreault et al. [32] isolated SARS-CoV-2 from trachea samples at 4 dpi, while RNA was detected in respiratory and lymphoid tissues at 4, 8 and 21 dpi. Meanwhile, Ulrich et al. [45] did not perform virus isolation from RT-qPCR positive samples obtained from SARS-CoV-2 inoculated cattle. Although our study performed in goats was carried out with a higher sample size than the mentioned experimental infections in livestock ruminants, comparable results regarding overall susceptibility in goats were observed. Low amounts of viral genome and antigen in tissues were found in the inoculated goats. Although the detected nucleic acid may be considered remnants of the inoculum, the detection of subgenomic RNA in NT and tonsil samples at 2 dpi suggests recent virus replication, at least in those tissues. In addition, the virus isolation from these tissues supports the presence of infectious virus in this species after challenge with the B.1.351/Beta variant. On the other hand, results from IHC exhibited presence of nucleocapsid protein of SARS-CoV-2 in the tonsil of a single goat at each of the time points (2, 4, and 7 dpi) in cells with morphology compatible with dendritic cells

Discussion
Goat is an important livestock animal that has already been suggested to be susceptible to SARS-CoV-2 by previous studies, both in vitro and in vivo [31,42]. Similarity between goat and human ACE2 receptors has been demonstrated via comparative genomic analysis as well as binding affinity between goat ACE2 and SARS-CoV-2 S protein [42]. Other previous in vitro studies also reported that goat ACE2 supports cell entry of SARS-CoV-2 [43,44], and its replication has been demonstrated in DF1 cells expressing both goat ACE2 and transmembrane serine protease 2 (TMPRSS2) [43]. Since the susceptibility of animal species could change after the appearance of SARS-CoV-2 variants and our in silico study predicted potential affinity between the B.1.351/Beta variant and ACE2 goat receptor, the aim of the present study was to investigate the susceptibility of goat to this variant of SARS-CoV-2 by means of an in vivo experimental inoculation.
Previously, Bosco-Lauth et al. [31] described SARS-CoV-2 RT-qPCR positive results in two out of three inoculated goats at 3 days after experimental challenge with the ancestral variant. However, none of the animals shed detectable infectious virus, nor did other ruminants such as sheep and cattle. The aforementioned study together with other previous studies reported that goat, sheep, or cattle developed limited subclinical infection after SARS-CoV-2 inoculation [20,32]. In sheep, Gaudreault et al. [32] isolated SARS-CoV-2 from trachea samples at 4 dpi, while RNA was detected in respiratory and lymphoid tissues at 4, 8 and 21 dpi. Meanwhile, Ulrich et al. [45] did not perform virus isolation from RT-qPCR positive samples obtained from SARS-CoV-2 inoculated cattle. Although our study performed in goats was carried out with a higher sample size than the mentioned experimental infections in livestock ruminants, comparable results regarding overall susceptibility in goats were observed. Low amounts of viral genome and antigen in tissues were found in the inoculated goats. Although the detected nucleic acid may be considered remnants of the inoculum, the detection of subgenomic RNA in NT and tonsil samples at 2 dpi suggests recent virus replication, at least in those tissues. In addition, the virus isolation from these tissues supports the presence of infectious virus in this species after challenge with the B.1.351/Beta variant. On the other hand, results from IHC exhibited presence of nucleocapsid protein of SARS-CoV-2 in the tonsil of a single goat at each of the time points (2, 4, and 7 dpi) in cells with morphology compatible with dendritic cells (DCs). It is likely that this immunolabelling represents viral internalisation and accumulation in these cells, since at least human DCs have been shown to trap SARS-CoV-2 even though not showing a productive infection [46,47]. Further studies would be needed to understand whether viral replication or DC internalisation occur in tonsils in this or other potentially susceptible species.
Differences in viral loads observed in our study compared to those of previous experimental infections in goats, as well as in the other ruminant species, may be due to the higher affinity predicted for B.1.351/Beta variant in goats than for ancestral variants [31,32]. This may confirm that the emergence of new variants raises concerns regarding the risk of new mutations in SARS-CoV-2 that could increase the range of susceptibility of zoonotic reservoirs.
In our study, the dose of challenge may be higher than the one that goats could be exposed to under natural conditions. Thus, it is unknown whether inoculation with a lower dose of virus would have shown evidence of infection or not. In a parallel study, we tested 208 goats that had been in contact with COVID-19 positive farmers, yet none of the goats developed humoral immune responses against SARS-CoV-2 (unpublished work) supporting that they were not exposed to the virus albeit the circumstances.
Overall, our results suggest that domestic goat has low susceptibility to SARS-CoV-2 B.1.351/Beta variant infection with low amounts of viral genome and antigen in tissues and evidence of seroconversion from 7 dpi onwards. Moreover, challenged goats did not show clinical signs or gross and microscopic lesions consistent with SARS-CoV-2 infection. Thus, the domestic goat seems to be a poorly competent host for SARS-CoV-2 (at least for the Beta variant) and probably has a negligible role in virus transmission, in agreement with studies carried out in other ruminant livestock species [31,32,46]. However, it is important to continue monitoring potential susceptible wild and domestic species to minimize the risk of transmission at the human-animal interface. Moreover, the circulating variants such as Omicron and the emergence of new ones raise concerns regarding the risk of new mutations in SARS-CoV-2 that could expand host range susceptibility and generate new zoonotic reservoirs. In fact, a recent study with 27 ACE2 orthologues including goat suggests a broader species receptor binding of Omicron compared to other variants such as Delta [48]. Thus, the infectivity of new viral variants in different animal species should be continuously monitored as improved viral replication cannot be ruled out in species currently defined as poorly susceptible.