# Patterns of genital tract mustelid gammaherpesvirus 1 (Musghv-1) reactivation are linked to stressors in european badgers (Meles Meles)

Ming-Shan Tsai, Sarah François, Chris Newman, David W. Macdonald, Christina D. Buesching

► To cite this version:

Ming-Shan Tsai, Sarah François, Chris Newman, David W. Macdonald, Christina D. Buesching. Patterns of genital tract mustelid gammaherpesvirus 1 (Musghv-1) reactivation are linked to stressors in european badgers (Meles Meles). Biomolecules, 2021, 11 (5), pp.716. 10.3390/biom11050716. hal-04338374

HAL Id: hal-04338374

https://hal.inrae.fr/hal-04338374

Submitted on 12 Dec 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Patterns of Genital Tract *Mustelid Gammaherpesvirus 1* (Musghv-1) Reactivation Are Linked to Stressors in European Badgers (*Meles Meles*)

Ming-shan Tsai <fnref n="1" />,*, Sarah François <fnref n="2" />, Chris Newman <fnref n="1" />,<fnref n="3" />, David W. Macdonald <fnref n="1" /> and Christina D. Buesching <fnref n="3" />,<fnref n="4" />

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Tsai, M.-s.; François, S.; Newman, C.; Macdonald, D.W.; Buesching, C.D. Patterns of Genital Tract *Mustelid Gammaherpesvirus 1* (Musghv-1) Reactivation Are Linked to Stressors in European Badgers (*Meles Meles*. *Biomolecules* **2021**, *11*, 716. https://doi.org/10.3390/biom11050716

Academic Editor: Marshall Williams

Received: 6 April 2021

Accepted: 6 May 2021

Published: 11 May 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Gammaherpesvirus reactivation can promote diseases or impair reproduction. Understanding reactivation patterns and associated risks of different stressors is therefore important. Nevertheless, outside the laboratory or captive environment, studies on the effects of stress on gammaherpesvirus reactivation in wild mammals are lacking. Here we used *Mustelid gammaherpesvirus 1* (MusGHV-1) infection in European badgers (*Meles meles*) as a host-pathogen wildlife model to study the effects of a variety of demographic, physiological and environmental stressors on virus shedding in the genital tract. We collected 251 genital swabs from 150 free-ranging individuals across three seasons and screened them for the presence of MusGHV-1 DNA using PCR targeting the DNA polymerase gene. We explored possible links between MusGHV-1 DNA presence and seven variables reflecting stressors, using logistic regression analysis. The results reveal different sets of risk factors between juveniles and adults, likely reflecting primary infection and reactivation. In adults, virus shedding was more likely in badgers in poorer body condition and younger than 5 years or older than 7; while in juveniles, virus shedding is more likely in females and individuals in better body condition. However, living in social groups with more cubs was a risk factor for all badgers. We discuss possible explanations for these risk factors and their links to stress in badgers.

**Keywords:** sexually transmittable infections; wildlife disease; epidemiology; seasonal effect; weather conditions

## 1. Introduction

Alongside the *Alphaherpesvirinae*, *Betaherpesvirinae* and *Deltaherpesvirinae*, the Gammaherpesvirinae is a subfamily of double-stranded DNA viruses belonging to the family Herpesviridae [1,2] containing four genera, the *Macavirus*, *Rhadinovirus*, *Lymphocryptovirus* and *Percavirus*. After primary infection, most known gammaherpesviruses establish a latent stage in the B lymphocytes, and can then be reactivated repeatedly throughout life. Numerous factors are known to trigger reactivation, including stress [3,4], trauma (e.g., surgery: [5]) or primary co-infection with other pathogens [6]. Reactivation is a process of viral lytic infection, which involves virus replication within the host cell, eventually causing the cells to burst and release infectious virions. Reactivation of gammaherpesviruses occurs in plasma cells and epithelial cells of mucosa that function as portals for external contact (e.g., mouth, nose, eyes and genital tract), thus facilitating transmission [7]. Reactivation is typically asymptomatic or induces only mild disease, but, depending on

<footnote n="1">¹ Recanati-Kaplan Centre, Wildlife Conservation Research Unit, Department of Zoology, University of Oxford, Abingdon Road, Tubney House, Tubney, Oxfordshire OX13 5QL, UK; chris.newman@lmh.ox.ac.uk (C.N.); david.macdonald@zoo.ox.ac.uk (D.W.M.)</footnote>

<footnote n="2">² Evolve.Zoo, Peter Medawar Building for Pathogen Research, Department of Zoology, University of Oxford, South Park Road, Oxford OX1 3SY, UK; sarah.francois@zoo.ox.ac.uk</footnote>

<footnote n="3">³ Cook's Lake Farming Forestry and Wildlife Inc. (Ecological Consultancy), Queens County, NS B0J 2H0, Canada; christina.buesching@ubc.ca</footnote>

<footnote n="4">⁴ Department of Biology, Irving K. Barber Faculty of Science, University of British Columbia, Kelowna, BC V1V 1V7, Canada</footnote>

<footnote>* Correspondence: mingshan.tsai@zoo.ox.ac.uk</footnote>
strain pathogenicity [8], can also contribute to the development of other severe diseases such as cancer [9], and is associated with a higher risk of contracting co-infections with additional highly virulent pathogens (e.g., herpesvirus DNA detection is a risk factor for the presence of *Chlamydia pecorum* DNA in koalas [10]) or pathogens causing immunodeficiency (e.g., Human herpesvirus 8 and HIV [11]). Although not as extensively studied as the human Epstein-Barr Virus (EBV) and Kaposi's sarcoma-associated herpesvirus (KSHV), gammaherpesvirus reactivation has long been linked to sometimes severe diseases in domestic animals [12–14], and also increasingly to illness in wildlife species [15–18].

European badgers, *Meles meles* (henceforth ‘badgers’), are seasonally breeding mustelids that are commonly infected with the *Mustelid gammaherpesvirus 1* (MusGHV-1, genus *Percavirus*), for which prevalence of viral DNA in blood samples can reach up to 100% in the UK and in Ireland across all age groups, with primary infection likely occurring in early life [19,20]. Furthermore, MusGHV-1 DNA has also been detected in 55–82.5% badger genital swab samples [21,22], suggesting local genital tract reactivation is common, and thus MusGHV-1 is a sexually transmittable pathogen. Sexually transmitted infections (STIs) are often associated with infertility in animals [23,24], and in badgers previous research has linked otherwise asymptomatic MusGHV-1 shedding in genital tracts to poorer female reproductive performance [22]. Considering the potentially considerable impact of STI-induced sterility on population dynamics and potential to drive extinction [25,26], understanding the patterns of MusGHV-1 genital tract reactivation is important in identifying the risk factors. Two previous studies on genital tract MusGHV-1 screening have indicated that cubs, old individuals and adult males during the mating season are at particular risk of genital MusGHV-1 shedding [21,22]. Nevertheless, the impact of social, physiological and seasonal effects, which may contribute to immune dysregulation due to stress, on MusGHV-1 reactivation has not been investigated.

Since the local extinction of brown bears (*Ursus arctos*) and wolves (*Canis lupus*) [27], badgers do not have natural predators in the UK, and their survival is mainly dependent on weather conditions [28,29], which are closely linked to food availability. Faecal corticoid levels indicate that badgers experience seasonal variation in stress levels [30,31]. There are several potential underlying causes: seasonal stress may be linked to seasonally varying food availability [30,31] that can result in mortality [28,29]. Higher cortisol levels are also associated with poor body condition and coinfection with *Mycobacterium bovis*, the pathogen causing bovine tuberculosis [32]. Therefore, metrics of body condition, especially reduced body-condition as a consequence of recent lactation [30,33], can indicate that individuals experience physiological stress [30]. Sociologically, higher social group density is associated with female reproductive suppression [34,35], reduced body condition and fecundity [36] and increased bite wounding among male badgers [37]. Furthermore, aging reduces tolerance to stress [38], specifically altering the balance between innate and acquired immunity in badgers [39], and increasing their risk of herpesvirus reactivation [22], as also observed in other carnivora species [40,41], sometimes resulting in chronic or continuous herpesvirus reactivation [42].

Here we conducted population-wide molecular screening using genital swabs taken from a free-ranging badger population in the south of England across three seasons (spring: May, summer: September, autumn: November). We investigated whether environmental factors (i.e., season, social group size and structure), host demographic parameters (i.e., sex, age and lactation) and host body conditions are risk factors for genital MusGHV-1 shedding. Furthermore, we investigated if risk factors differ between juveniles and adults, given that virus shedding in juveniles is more likely associated with primary infection, whereas in adults it is likely due to reactivation from latency.

## 2. Materials and Methods

Samples were collected from live-trapped badgers in Wytham Woods, Oxfordshire, UK (51°46'26" N, 1°19'19" W; caught in May, September and November 2018 following the methodology described in Macdonald et al. [43]; for details see Table 1). For each
capture, we recorded sex, sett (i.e., communal den used by a social group) of capture, body condition score (BCS, categorised as 1 = very thin to 5 = very fat) and lactational status (determined by teat measurements of females in spring: [44]). Because each badger in Wytham is given an individual tattoo at first capture (usually as a cub [45]), exact age (in years) was known for most (243 of 251) animals in the dataset. For the remaining eight badgers first caught as adults, age was inferred by toothwear according to the method described by Bright Ross et al. [46]. We defined four age classes based on sex-steroid levels [47] and age: (i) juveniles < 2 years old (cubs and yearlings were combined to increase sample sizes, as there was no difference between MusGHV-1 prevalence in cubs and yearlings: Fisher's exact test: *p*-value = 0.7449); (ii) young adults: 2 ≤ *x* < 5 years old; (iii) old adults: 5 ≤ *x* < 8 years old; (iv) very old adults: ≥8 years old. The number of cubs and adults resident in each sett was estimated using minimum number alive (MNA) estimates [43,46].

<table><caption>**Table 1.** Summary of sampling effort (badgers with unknown age were assigned to age group by tooth wear according to the method described by Bright Ross et al. [46]).</caption><thead><tr><th rowspan="2">Age Group</th><th rowspan="2">Age</th><th colspan="2">Spring</th><th colspan="2">Summer</th><th colspan="2">Autumn</th><th rowspan="2">Total</th></tr><tr><th>Female</th><th>Male</th><th>Female</th><th>Male</th><th>Female</th><th>Male</th></tr></thead><tbody><tr><td>Cub</td><td>0</td><td>15</td><td>15</td><td>10</td><td>4</td><td>9</td><td>9</td><td>62</td></tr><tr><td>Yearling</td><td>1</td><td>1</td><td>3</td><td>2</td><td>2</td><td></td><td>3</td><td>11</td></tr><tr><td rowspan="3">Young</td><td>2</td><td>7</td><td>7</td><td>7</td><td>9</td><td>5</td><td>6</td><td>41</td></tr><tr><td>3</td><td>4</td><td>8</td><td>3</td><td>5</td><td>2</td><td>8</td><td>30</td></tr><tr><td>4</td><td>2</td><td>5</td><td>4</td><td>5</td><td>2</td><td>8</td><td>26</td></tr><tr><td rowspan="4">Old</td><td>5</td><td>2</td><td>1</td><td>6</td><td>1</td><td>3</td><td>1</td><td>14</td></tr><tr><td>6</td><td>3</td><td>1</td><td>3</td><td>1</td><td>3</td><td></td><td>11</td></tr><tr><td>7</td><td></td><td>2</td><td></td><td>2</td><td></td><td>2</td><td>6</td></tr><tr><td>Unknown</td><td>1</td><td>1</td><td>1</td><td>1</td><td></td><td>2</td><td>6</td></tr><tr><td rowspan="6">Very old</td><td>8</td><td>4</td><td>4</td><td>4</td><td>3</td><td>1</td><td>3</td><td>19</td></tr><tr><td>9</td><td>2</td><td>1</td><td>2</td><td>1</td><td>1</td><td>1</td><td>8</td></tr><tr><td>10</td><td>3</td><td>1</td><td>3</td><td></td><td>4</td><td></td><td>11</td></tr><tr><td>11</td><td></td><td></td><td></td><td></td><td></td><td></td><td>0</td></tr><tr><td>12</td><td></td><td></td><td></td><td></td><td></td><td></td><td>0</td></tr><tr><td>13</td><td>1</td><td>1</td><td></td><td>1</td><td>1</td><td></td><td>4</td></tr><tr><td></td><td>Unknown</td><td></td><td>1</td><td></td><td></td><td>1</td><td></td><td>2</td></tr></tbody><tfoot><tr><td colspan="2">Total</td><td>45</td><td>51</td><td>45</td><td>35</td><td>31</td><td>44</td><td>251</td></tr></tfoot></table>

Sterile cotton tops with wooden shafts were used to swab the genital tracts of all females (cubs and adults) and all males (except for very small male cubs in spring for animal welfare reasons), and stored in 2 mL sterile microcentrifuge tubes. All samples were frozen and stored at -20 °C immediately after sampling. Badgers were released at their site of capture on the same day, after full recovery from anaesthesia.

Each swab was reconstituted with 400 µL sterile double distilled water and vortexed gently at room temperature for 10 min. A 200 µL aliquot was taken from the reconstituted swab fluids, and viral DNA was extracted and purified using a commercial kit (DNeasy Blood and Tissue Kit, Qiagen, Hilden, Germany) following manufacturer's instructions. Purified DNA was then eluted in 100 µL of the provided buffer.

The purified DNA was screened using a MusGHV-1-specific primer pair designed by King et al. 2004 [20], targeting 281 base pairs of the partial DNA polymerase gene. For each reaction, a total of 20 µL PCR solution was mixed with 10 µL HotStartTaq Master Mix (Qiagen, containing 1 unit of HotStartTaq DNA Polymerase, 12 µM of MgCl₂ and 1.6 µM of each dNTP), 0.5 µM of each primer and 2 µL CoralLoad gel loading dye and 5 µL DNA template. Amplification conditions were kept at 95 °C for 5 min to activate DNA polymerase, followed by 45 cycles of denaturation at 95 °C for 45 s, primer annealing
at 60 °C for 45 s, and chain elongation at 72 °C for 1 min, followed by a final extension at 72 °C for 10 min. Finally, the PCR products were loaded in 2% agarose gel to check the amplification results under UV light. Five randomly selected samples with positive results were then amplified again with substituted front primer (5' CCA AGC AGT GCA TAG GAG GT 3') to generate longer sequences (771 base pairs). PCR products were then purified and sent for genotyping using Sanger sequencing to confirm the identity of produced amplicon. Sequences returned were then aligned by Clustal W method [48] and analyzed for variation using MEGA X (10.1.7) [49]. Representative sequences were selected and deposited on GenBank under accession numbers MT332100 and MT332101.

Statistical analyses were performed with the R and R Studio software (version 1.21335) [50,51]. Detection rate of genital MusGHV-1 DNA was calculated by dividing the number of PCR positive cases by the total number of tested cases, and 5% upper and lower confidence intervals were calculated using the Wilson method [52]. Logistic regression (glmer function, R package lme4) with badger identity (tattoo) number as a random effect was used to measure univariate effects of MusGHV-1 reactivation in genital tracts with season, sex, age class, BCS, number of residents per social group (total, adult and cub MNA), and percentage of cubs per social group (calculated by cub MNA divided by total MNA) as fixed effect variables. We categorised the percentage of cubs per social group into low and high using 30% as the dividing point according to the distribution of these data (Figure S1). Effects of lactation were analyzed using Fisher's exact tests, due to low sample sizes. Because juveniles (especially cubs) are generally thinner than adults, and thus have a different body condition distribution [36], in the univariate analysis we analyzed the association of MusGHV-1 detection and individual BCS for juveniles and adults separately. The effect size of each variable is presented using odds ratios. Variables with *p*-values < 0.1 were retained in the multivariable analysis. We constructed two multivariable models based on data from juveniles and adults as we expected patterns to be different between these two groups. Manual backwards selection was used to select the final model with best fit (determined by lowest AIC value). Model residual diagnostics were conducted using the R package DHARMa (version 0.3.3.0). Model fit was established using area-under-receiver-operating characteristics (AUC) with an acceptance threshold of 0.7 [53]. Collinearity of variables was checked using variance inflation factor (VIF), where variables with a VIF value of more than 5 were excluded from the model. We defined variables in the final models with *p*-values of less than 0.1 as risk factors for genital tract MusGHV-1 shedding.

# 3. Results

## 3.1. *Univariate Analysis*

The overall detection rate of genital tract MusGHV-1 DNA was 35.5% (89/251, 95% CI: 30–41.6%), and was generally higher in summer (43.8%, 35/80) than in spring (34.4%, 33/96) and autumn (28%, 21/75) (Table 2). There was strong evidence for an effect of age on genital MusGHV-1 shedding, where genital tract MusGHV-1 DNA detection rate in juveniles (45.2%, 33/73) and very old badgers (47.7%, 21/44) was higher than in middle-aged adults (i.e., young (30.9%, 30/97) and old (13.5%, 5/37) individuals; Table 2). However, no effect of sex was observed (logistic regression analysis, *p* = 0.679). Adults with lower BCS had a higher probability of genital MusGHV-1 DNA detection rate (*p*-value = 0.012, Table 2). There was no evidence that the number of adults, or the total number of badgers resident in each social group, affected MusGHV-1 DNA detection rates, except for the number of cubs per social group (Table 2): the prevalence of genital tract MusGHV-1 was significantly higher in badgers living in social groups comprising >30% cubs (logistic regression analysis, *p*-value = 0.004). Finally, there was no evidence for recent lactation affecting the MusGHV-1 DNA detection rate in sexually mature females (Fisher's exact test, *p*-value = 1; Table 2).
<table><caption>**Table 2.** Overview of MusGHV-1 DNA detection rate in badger genital tract and univariate logistic regression analysis.
Formula: MusGHV-1 ~ Variate + (1 | Tattoo); number of observations: 251; groups by tattoo number: 150.</caption><thead><tr><th>Variable</th><th>Positive</th><th>Total</th><th>%</th><th>95% CI</th><th>Odds Ratio (OR)</th><th>OR 95% CI</th><th>p-Value</th></tr></thead><tbody><tr><td>Sex</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Male</td><td>45</td><td>130</td><td>34.62%</td><td>27-43.1%</td><td></td><td></td><td></td></tr><tr><td>Female</td><td>44</td><td>121</td><td>36.36%</td><td>28.3-45.2%</td><td>0.93</td><td>0.55-1.56</td><td>0.774</td></tr><tr><td>Season</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Spring</td><td>33</td><td>96</td><td>34.38%</td><td>25.6-44.3%</td><td>1.36</td><td>0.7-2.65</td><td>0.368</td></tr><tr><td>Summer</td><td>35</td><td>80</td><td>43.75%</td><td>33.4-54.7%</td><td>2.02</td><td>1.02-4.03</td><td>0.044</td></tr><tr><td>Autumn</td><td>21</td><td>75</td><td>28.00%</td><td>19.1-39%</td><td></td><td></td><td></td></tr><tr><td>Age group</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Juvenile (&lt;2 years old)</td><td>33</td><td>73</td><td>45.21%</td><td>34.3-56.6%</td><td>5.28</td><td>1.85-15.08</td><td>0.001</td></tr><tr><td>Young (2-4 years old)</td><td>30</td><td>97</td><td>30.93%</td><td>22.6-40.7%</td><td>2.87</td><td>1.02-8.08</td><td>0.046</td></tr><tr><td>Old (5-7 years old)</td><td>5</td><td>37</td><td>13.51%</td><td>5.9-27.8%</td><td></td><td></td><td></td></tr><tr><td>Very old (&gt;7 years old)</td><td>21</td><td>44</td><td>47.73%</td><td>33.8-62%</td><td>5.9</td><td>1.92-17.78</td><td>0.002</td></tr><tr><td>Body condition <sup>a</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Body condition score (1-5)</td><td></td><td>182</td><td></td><td></td><td>0.64</td><td>0.45-0.9</td><td>0.012</td></tr><tr><td>Social group size</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Total</td><td></td><td>251</td><td></td><td></td><td>1.02</td><td>0.94-1.1</td><td>0.653</td></tr><tr><td>Adult</td><td></td><td>251</td><td></td><td></td><td>0.96</td><td>0.88-1.05</td><td>0.405</td></tr><tr><td>Cub</td><td></td><td>251</td><td></td><td></td><td>1.18</td><td>1.01-1.38</td><td>0.035</td></tr><tr><td>Cub percentage per sett</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Low (&lt;30%)</td><td>58</td><td>191</td><td>30.37%</td><td>24.3-37.2%</td><td></td><td></td><td></td></tr><tr><td>High (&gt;30%)</td><td>31</td><td>60</td><td>51.67%</td><td>39.3-63.8%</td><td>2.35</td><td>1.03-4.24</td><td>0.004</td></tr><tr><td>Lactational status <sup>bc</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Not Lactated</td><td>3</td><td>10</td><td>30.00%</td><td>10.8-60.3%</td><td></td><td></td><td></td></tr><tr><td>Lactated</td><td>6</td><td>18</td><td>33.33%</td><td>16.3-56.3%</td><td>0.9</td><td>0.16-4.92</td><td>0.856</td></tr></tbody></table>

## 3.2. Multivariable Analysis

We generated multivariate models of MusGHV-1 detection for juveniles (n = 72) and adults (n = 171). The most parsimonious models for juveniles (Table 3; Figure S2) and adults (Table 4; Figure S3) had both diagnostically acceptable [53] AUC area of 0.727 and 0.774, respectively. In juveniles, multivariate analysis showed that being a female, being in better body condition and living in a social group with a high proportion of cubs (over 30%) were risk factors for genital MusGHV-1 shedding (Table 3). In contrast, in adults, young (2–5 years old) or very old (≥8 years) badgers in poorer body condition, and living in social groups with a high proportion of cubs, are at particular risk of shedding MusGHV-1 in their genital tract (Table 4).

<footnote><sup>a</sup>: Only adults were included in this analysis; <sup>b</sup>: only females captured in spring were included in this analysis; <sup>c</sup>: Fisher exact test.</footnote>
<table><caption>**Table 3.** Final general mixed effect model of multivariable logistic regression analysis for juveniles. Formula: MusGHV-1 ~ Sex + Body condition + (1 | Tattoo); number of observations: 72; groups by tattoo number: 48.</caption><thead><tr><th>Group</th><th>Estimate</th><th>Standard Error</th><th>z Value</th><th>Adjusted OR</th><th>95% CI</th><th>p Value</th></tr></thead><tbody><tr><td>(Intercept)</td><td>-1.204</td><td>0.813</td><td>-1.482</td><td></td><td></td><td></td></tr><tr><td>Sex</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Female</td><td>-0.985</td><td>0.528</td><td>-1.863</td><td>0.37</td><td>0.13-1.05</td><td>0.062</td></tr><tr><td>Male</td><td>-0.985</td><td>0.528</td><td>-1.863</td><td>0.37</td><td>0.13-1.05</td><td>0.062</td></tr><tr><td>Body Condition</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>BCS</td><td>0.785</td><td>0.306</td><td>2.565</td><td>2.19</td><td>1.2-3.99</td><td>0.01</td></tr><tr><td>Cub percentage</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Low (&lt;30%)</td><td>-0.977</td><td>0.539</td><td>-1.812</td><td>0.38</td><td>0.13-1.08</td><td>0.07</td></tr><tr><td>High (&gt;30%)</td><td>-0.977</td><td>0.539</td><td>-1.812</td><td>0.38</td><td>0.13-1.08</td><td>0.07</td></tr></tbody></table>

<table><caption>**Table 4.** Final general mixed effect model of multivariable logistic regression analysis for adults. Formula: MusGHV-1 ~ Sex + Season + AgeGroup + Cub percentage + Sex * Cub percentage + (1 | Tattoo); number of observations: 171; groups by tattoo number: 99.</caption><thead><tr><th>Group</th><th>Estimate</th><th>Standard Error</th><th>z Value</th><th>Adjusted OR</th><th>95% CI</th><th>p Value</th></tr></thead><tbody><tr><td>(Intercept)</td><td>0.059</td><td>1.075</td><td>0.055</td><td></td><td></td><td></td></tr><tr><td>Season</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Spring</td><td>-0.289</td><td>0.584</td><td>-0.495</td><td>0.74</td><td>0.24-2.35</td><td>0.62</td></tr><tr><td>Summer</td><td>0.684</td><td>0.520</td><td>1.318</td><td>1.98</td><td>0.72-5.48</td><td>0.188</td></tr><tr><td>Autumn</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Age Group</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Young</td><td>1.767</td><td>0.673</td><td>2.627</td><td>5.85</td><td>1.57-21.87</td><td>0.009</td></tr><tr><td>Old</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Very old</td><td>1.75</td><td>0.695</td><td>2.522</td><td>5.77</td><td>1.48-22.55</td><td>0.012</td></tr><tr><td>Body condition</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>BCS</td><td>-0.60</td><td>0.238</td><td>-2.513</td><td>0.54</td><td>0.34-0.88</td><td>0.012</td></tr><tr><td>Cub percentage</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Low (&lt;30%)</td><td>-1.15</td><td>0.544</td><td>-2.115</td><td>0.31</td><td>0.11-0.92</td><td>0.034</td></tr><tr><td>High (&gt;30%)</td><td>-1.15</td><td>0.544</td><td>-2.115</td><td>0.31</td><td>0.11-0.92</td><td>0.034</td></tr></tbody></table>

### 3.3. Genetic Diversity of MusGHV-1 in the Wytham Badger Population

We sequenced five randomly selected MusGHV-1 positive PCR products of the partial DNA polymerase gene. All sequences were trimmed to 694 base pairs and confirmed to be MusGHV-1 according to the NCBI BLAST results, returning 98.7% (n = 3) and 100% (n = 2) nucleotide identity to the published MusGHV-1 sequence isolated from a badger in Cornwall, England (accession number: AF275657).

## 4. Discussion

In this study, we investigated the patterns of MusGHV-1 shedding in badger genital tracts, and identified associated risk factors comprising age, sex, season and social group composition. We confirmed that juveniles and adults are subject to different sets of risk factors, putatively reflecting primary infection and reactivation from latency, respectively.
In terms of age class effects, the high genital reactivation rate detected in cubs and yearlings suggests that badgers contract MusGHV-1 early in life, before reaching sexual maturity. Although MusGHV-1 reactivates repeatedly throughout life, reactivation tends to be less frequent in young and old adults, although rates increase in very old individuals. This matches patterns in humans where most people become infected with herpes during their childhood/adolescence (e.g., 100% and 70% seroprevalence of EBV before age 14 in Hong Kong and the United Kingdom [54]), respectively, then typically experience viral latency during their prime, but can suffer from increasingly longer and more frequent herpesvirus reactivation that can cause mild disease (e.g., shingles [55] in old age due to lowered immune response [56,57]).

Since vertical transmission of MusGHV-1 through the placenta is unlikely [22], and the potential for infection from the vaginal tract during parturition is equally low, due to low genital MusGHV-1 detection rate in pregnant females [22], we hypothesise that cubs contract MusGHV-1 through close contact with virus-shedding conspecifics [22]. Thereafter genital virus shedding in cubs may arise after primary acute infection through non-sexual routes and subsequent latency, as observed in the murine model where *Murine herpesvirus 4* (MuHV-4, also a gammaherpesvirus), inoculation in the nasal cavity results in acute infection in the respiratory tract and lungs, and establishes latency in the spleen, but then reactivates in the vaginal tract 17-21 days after inoculation [58]. MuHV-4 nasal cavity inoculation, however, does not result in reactivation in male genital tracts, and transmission is only possible from females to males. After sexual intercourse with virus-shedding female mice, the virus replicates in the male penis for 3 weeks. Interestingly, also in badgers, female juveniles are at higher risk of MusGHV-1 reactivation in the genital tract than males. Furthermore, juveniles in better body condition, regardless of sex, exhibit higher virus detection rate. Indeed, juvenile males in better body condition enter puberty earlier than thinner males (11 months compared to 22-28 months: [59]). Once juveniles enter puberty they will experience an increased risk of contracting MusGHV-1 through sexual contact and/or that their latent infection is reactivated through mating resulting in viral shedding in the genital tract.

Faecal corticosteroid measurements from badgers in Ireland [30] evidence higher stress levels in summer likely associated with dry environmental conditions that result in lower earthworm availability (i.e., the badgers' main food type [45]); similarly, in our own study population, summer drought is an established mortality factor due to malnutrition/starvation [28,29]. This corresponds to our finding that, in all adults, seasonal MusGHV-1 reactivation rates were highest in summer, but females tended to have higher viral reactivation levels than males in spring—possibly due to reproductive stresses, while the reverse was true in autumn (Table 2) [30]. Badgers captured in summer generally have lower BCS than in other seasons. This is reflected in our multivariate model (Table 4) such that when season and BCS are both included in the model, the effect of seasons found in univariate model (Table 1) is eliminated, while BCS maintains a strong negative effect on MusGHV-1 DNA presence in the genital tract. This implies that BCS is the main driver causing the seasonal effect of MusGHV-1 DNA detection rate in the genital tract. However, we did not detect any collinearity between seasons and BCS, and removing season from the model decreased the model fit (i.e., an increase in the AIC value), suggesting season contributed to the response value through a different explanatory variable that was not included in our study.

Chronic stress has proven a significant risk factor of herpesvirus reactivation, causing immune system dysregulation, where corticosteroids inhibit the pro-inflammatory cytokine responses, allowing viruses to (re-)activate and undergo lytic proliferation unchecked [60]. This link between elevated corticosteroid levels and herpesvirus reactivation has been proven experimentally (horses: [61]; captive reindeer: [62]) and through observation (e.g., humans: [4,63]; captive Grévy's zebras (*Equus grevyi*): [3]), but has not been investigated in free-living wildlife populations. Herpesvirus reactivation triggered by stress has been widely confirmed naturally and experimentally by corticosteroid injection in humans
and domestic animals [3,61]. Linking stress and viral reactivation in wildlife, however, is particularly challenging due to the difficulties of monitoring individual stress levels in the field, and typically this relationship can only be confirmed experimentally by taking subjects into captivity [62]. Using indicators that have been linked to stress hormone levels in previous studies can thus provide an informative way to study the relationship between stress and herpesvirus reactivation in free-ranging wildlife.

Our finding of no sex-related effect on herpesvirus reactivation in adults is at odds with our recent survey of Irish badger populations [22], where the rate (over 82.5%) of genital MusGHV-1 DNA presence in adult males was significantly higher than in females (47.5%), during the peak (postpartum) mating season (i.e., from mid-January to mid-February). This implies not only a link to mating activity, but also a mechanism enhancing sexual transmission where sex-biases in viral shedding in badger genital tracts are most pronounced during the mating season. This corroborates our finding in the same study where males with more spermatozoa have a higher detection rate of genital MusGHV-1 DNA [22]; linking higher sexual activity to higher STI prevalence as predicted by theory [64] and reported also in many other studies [65–68]. We were unable to test effects relating to badger mating behaviour explicitly because we did not trap individuals during the mating season, which coincides with late pregnancy and neonatal cub care, where we sought to avoid stressing mothers and to avoid depriving cubs of maternal care. Nevertheless, reactivation rate in autumn, when food sources are most abundant, and badgers undergo a period of reproductive quiescence [69], and thus experience less implicit stress, was significantly lower compared only to summer, but not to spring. This suggests that other factor(s) (e.g., sex hormone cycles [58], oxidative stress [70], and/or genital microbiome [71] and linked coinfections [72]) might also be affecting reactivation rates, beyond the scope of our current study.

Our results also show that social group structure can affect the risk of genital her-
pesvirus shedding, particularly the higher proportion of cubs within a co-resident social
group. This may be because badger cubs generally carry higher pathogen burdens than
adults (particularly lice *Trichodectes melis* and coccidia *Eimeria melis*) [73,74], and thus
increase per capita immunity burden among all badgers resident in the respective sett.

**5. Conclusions**

Our study demonstrates, for the first time in the wild, the link between host stressors and herpesvirus shedding, which likely reflects both primary infection and reactivation from latency according to host age. Amplified stress levels induced by human disturbance as well as food insecurity and more frequent severe weather events arising from human-induced rapid environmental change could therefore not only increase the risk of disease development and promote pathogen transmission within a population, but also negatively impact host reproductive fitness through localised herpesvirus reactivation in the reproductive tract. Careful monitoring of endemic latent virus infection as well as surveillance for possible newly emerging strains should therefore be included when planning in situ and ex situ conservation programmes for endangered species.

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/biom11050716/s1. Figure S1: Data distribution of cub percentage (n = 251) according to group density. Figure S2: Residual diagnostic plot of the final multivariable logistic regression model of juvenile MusGHV-1 shedding in genital tract. Figure S3: Residual diagnostic plot of the final multivariable logistic regression model of adult MusGHV-1 shedding in genital tract. Data S1: A spreadsheet of raw data used for analysis in the present study.

**Author Contributions:** Conceptualization, M.-s.T., C.D.B. and D.W.M.; methodology, M.-s.T. and S.F.; formal analysis, M.-s.T.; investigation, M.-s.T., C.D.B. and C.N.; resources, D.W.M.; data curation, M.-s.T., C.D.B. and C.N.; writing—original draft preparation, M.-s.T.; writing—review and editing, M.-s.T., C.D.B., C.N., S.F. and D.W.M.; visualization, M.-s.T.; supervision, C.D.B., S.F. and D.W.M.; project administration, M.-s.T. All authors have read and agreed to the published version of the manuscript.
**Funding:** This research received no external funding.

**Institutional Review Board Statement:** All trapping and animal sampling protocols were approved by the University of Oxford' Animal Welfare and Ethical Review Board. Trapping was conducted under Natural England license (currently 2019-38863, Badger Act 1992) and all animal handling procedures were carried out by qualified Personal Individual License (PIL) holders under Home Office license (current PPL 30/3379, Animals (Scientific Procedures) Act 1986).

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Data are available from Supplementary Materials Data S1.

**Acknowledgments:** The authors would like to thank Nadine Sugianto, Sil van Lieshout, Tanesha Allen and Julius Bright Ross for assistance with sample collection, fieldwork and MNA data provision. MST would like to thank the Ministry of Education in Taiwan and Lady Margaret Hall, University of Oxford, for scholarship support. The authors also thank Paul Johnson and Ta-Chun Liu for providing advice on statistics. The authors claim no conflict of interests in the present work.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Roizman, B.; Baines, J. The diversity and unity of herpesviridae. *Comp. Immunol. Microbiol. Infect. Dis.* **1991**, *14*, 63–79. [CrossRef]

2. Pellett, P.E. Trunkloads of Viruses. *J. Virol.* **2014**, *88*, 13520–13522. [CrossRef] [PubMed]

3. Seeber, P.A.; Quintard, B.; Sicks, F.; Dehnhard, M.; Greenwood, A.D.; Franz, M. Environmental stressors may cause equine herpesvirus reactivation in captive Grévy's zebras (*Equus grevyi*). *PeerJ* **2018**, *2018*, 1–18. [CrossRef] [PubMed]

4. Stowe, R.P.; Pierson, D.L.; Barrett, A.D.T. Elevated stress hormone levels relate to Epstein-Barr virus reactivation in astronauts. *Psychosom. Med.* **2001**, *63*, 891–895. [CrossRef]

5. Grinde, B. Herpesviruses: Latency and reactivation-viral strategies and host response. *J. Oral Microbiol.* **2013**, *5*, 22766. [CrossRef]

6. Donofrio, G.; Cavirani, S.; Van Santen, V.; Flammini, C.F. Potential secondary pathogenic role for bovine herpesvirus 4. *J. Clin. Microbiol.* **2005**, *43*, 3421–3426. [CrossRef]

7. Shannon-Lowe, C.; Adland, E.; Bell, A.I.; Delecluse, H.-J.; Rickinson, A.B.; Rowe, M. Features Distinguishing Epstein-Barr Virus Infections of Epithelial Cells and B Cells: Viral Genome Expression, Genome Maintenance, and Genome Amplification. *J. Virol.* **2009**, *83*, 7749–7760. [CrossRef]

8. Goodman, L.B.; Loregian, A.; Perkins, G.A.; Nugent, J.; Buckles, E.L.; Mercorelli, B.; Kydd, J.H.; Palù, G.; Smith, K.C.; Osterrieder, N.; et al. A point mutation in a herpesvirus polymerase determines neuropathogenicity. *PLoS Pathog.* **2007**, *3*, 1583–1592. [CrossRef]

9. Barton, E.; Mandal, P.; Speck, S.H. Pathogenesis and Host Control of Gammaherpesviruses: Lessons from the Mouse. *Annu. Rev. Immunol.* **2011**, *29*, 351–397. [CrossRef]

10. Stalder, K.; Vaz, P.K.; Gilkerson, J.R.; Baker, R.; Whiteley, P.; Ficorilli, N.; Tatarczuch, L.; Portas, T.; Skogvold, K.; Anderson, G.A.; et al. Prevalence and clinical significance of herpesvirus infection in populations of Australian marsupials. *PLoS ONE* **2015**, *10*, e0133807. [CrossRef]

11. Zhang, P.; Wang, J.; Zhang, X.; Wang, X.; Jiang, L.; Gu, X. Identification of AIDS-Associated Kaposi Sarcoma: A Functional Genomics Approach. *Front. Genet.* **2020**, *10*, 1–12. [CrossRef]

12. Chastant-Maillard, S. Impact of bovine herpesvirus 4 (BoHV-4) on reproduction. *Transbound. Emerg. Dis.* **2015**, *62*, 245–251. [CrossRef]

13. McLuckie, A.J.; Barrs, V.R.; Lindsay, S.; Aghazadeh, M.; Sangster, C.; Beatty, J.A. Molecular diagnosis of Felis catus gamma-herpesvirus 1 (FcaGHV1) infection in cats of known retrovirus status with and without lymphoma. *Viruses* **2018**, *10*, 128. [CrossRef]

14. Marenzoni, M.L.; Sforna, M.; Stefanetti, V.; Casagrande Proietti, P.; Brignone, L.; Del Sero, A.; Falcioni, F.; Orvieto, S.; Tamantini, C.; Tiburzi, A.; et al. Detection of Equid herpesvirus type 2 and 5 DNA in uterine flushings of mares with reproductive disorders. *Vet. Microbiol.* **2014**, *174*, 570–576. [CrossRef]

15. Gagnon, C.A.; Tremblay, J.; Larochelle, D.; Music, N.; Tremblay, D. Identification of a novel herpesvirus associated with cutaneous ulcers in a fisher (*Martes pennanti*). *J. Vet. Diagn. Investig.* **2011**, *23*, 986–990. [CrossRef]

16. Tseng, M.; Fleetwood, M.; Reed, A.; Gill, V.A.; Harris, R.K.; Moeller, R.B.; Lipscomb, T.P.; Mazet, J.A.K.; Goldstein, T. Mustelid herpesvirus-2, a novel herpes infection in northern sea otters (*Enhydra lutris kenyoni*). *J. Wildl. Dis.* **2012**, *48*, 181–185. [CrossRef]

17. Nicolas de Francisco, O.; Esperón, F.; Juan-Sallés, C.; Ewbank, A.C.; das Neves, C.G.; Marco, A.; Neves, E.; Anderson, N.; Sacristán, C. Neoplasms and novel gammaherpesviruses in critically endangered captive European minks (*Mustela lutreola*). *Transbound. Emerg. Dis.* **2020**, 1–13. [CrossRef]

18. Abade dos Santos, F.A.; Monteiro, M.; Pinto, A.; Carvalho, C.L.; Peleteiro, M.C.; Carvalho, P.; Mendonça, P.; Carvalho, T.; Duarte, M.D. First description of a herpesvirus infection in genus Lepus. *BioRxiv* **2020**. [CrossRef]
19. Sin, Y.W.; Annavi, G.; Dugdale, H.L.; Newman, C.; Burke, T.; MacDonald, D.W. Pathogen burden, co-infection and major histocompatibility complex variability in the European badger (*Meles meles*). *Mol. Ecol.* **2014**, *23*, 5072–5088. [CrossRef]

20. King, D.P.; Mutukwa, N.; Lesellier, S.; Cheeseman, C.; Chambers, M.A.; Banks, M. Detection of Mustelid Herpesvirus-1 Infected European Badgers (*Meles meles*) in the British Isles. *J. Wildl. Dis.* **2004**, *40*, 99–102. [CrossRef]

21. Kent, A.; Ehlers, B.; Mendum, T.; Newman, C.; Macdonald, D.W.; Chambers, M.; Buesching, C.D. Genital tract screening finds widespread infection with mustelid gammaherpesvirus 1 in the European badger (*Meles meles*). *J. Wildl. Dis.* **2017**, *54*, 133–137. [CrossRef] [PubMed]

22. Tsai, M.S.; Fogarty, U.; Byrne, A.W.; O'keeffe, J.; Newman, C.; Macdonald, D.W.; Buesching, C.D. Effects of mustelid gammaherpesvirus 1 (Musghv-1) reactivation in european badger (*Meles meles*) genital tracts on reproductive fitness. *Pathogens* **2020**, *9*, 769. [CrossRef] [PubMed]

23. McLeod, D.V.; Day, T. Why is sterility virulence most common in sexually transmitted infections? Examining the role of epidemiology. *Evolution* **2019**, *73*, 872–882. [CrossRef] [PubMed]

24. Lockhart, A.B.; Thrall, P.H.; Antonovics, J. Sexually transmitted diseases in animals: Ecological and evolutionary implications. *Biol. Rev.* **1996**, *71*, 415–471. [CrossRef]

25. Ryder, J.J.; Miller, M.R.; White, A.; Knell, R.J.; Boots, M. Host-parasite population dynamics under combined frequency- and density-dependent transmission. *Oikos* **2007**, *116*, 2017–2026. [CrossRef]

26. Best, A.; Webb, S.; Antonovics, J.; Boots, M. Local transmission processes and disease-driven host extinctions. *Theor. Ecol.* **2012**, *5*, 211–217. [CrossRef]

27. Wilson, C.J. Could we live with reintroduced large carnivores in the UK? *Mamm. Rev.* **2004**, *34*, 211–232. [CrossRef]

28. Macdonald, D.W.; Newman, C.; Buesching, C.D.; Nouvellet, P. Are badgers "Under The Weather"? Direct and indirect impacts of climate variation on European badger (*Meles meles*) population dynamics. *Glob. Chang. Biol.* **2010**, *16*, 2913–2922. [CrossRef]

29. Nouvellet, P.; Newman, C.; Buesching, C.D.; Macdonald, D.W. A Multi-Metric Approach to Investigate the Effects of Weather Conditions on the Demographic of a Terrestrial Mammal, the European Badger (*Meles meles*). *PLoS ONE* **2013**, *8*. [CrossRef]

30. George, S.C.; Smith, T.E.; Mac Cana, P.S.S.; Coleman, R.; Montgomery, W.I. Physiological stress in the Eurasian badger (*Meles meles*): Effects of host, disease and environment. *Gen. Comp. Endocrinol.* **2014**, *200*, 54–60. [CrossRef]

31. Virgós, E.; Mangas, J.G.; Blanco-Aguiar, J.A.; Garrote, G.; Almagro, N.; Viso, R.P. Food habits of European badgers (*Meles meles*) along an altitude gradient of Mediterranean environments: A field test of the earthworm specialization hypothesis. *Can. J. Zool.* **2004**, *82*, 41–51. [CrossRef]

32. Cheeseman, C.L.; Little, T.W.A.; Mallington, P.J.; Page, R.J.C.; Wilesmith, J.W.; Pritchard, D.G. Population ecology and prevalence of tuberculosis in Badgers in an area of Staffordshire. *Mamm. Rev.* **1985**, *15*, 125–135. [CrossRef]

33. Woodroffe, R.; Macdonald, D.W. Costs of breeding status in the European badger, *Meles meles*. *J. Zool.* **1995**, *235*, 237–245. [CrossRef]

34. Woodroffe, R.; MacDonald, D.W. Female/Female Competition in European Badgers *Meles meles*: Effects on Breeding Success. *J. Anim. Ecol.* **1995**, *64*, 12. [CrossRef]

35. Creel, S.R.; Creel, N.M. Energetics, reproductive suppression and obligate communal breeding in carnivores. *Behav. Ecol. Sociobiol.* **1991**, *28*, 263–270. [CrossRef]

36. Macdonald, D.W.; Newman, C.; Stewart, P.D.; Domingo-Roura, X.; Johnson, P.J. Density-dependent regulation of body mass and condition in badgers (*Meles meles*) from Wytham Woods. *Ecology* **2002**, *83*, 2056–2061. [CrossRef]

37. Macdonald, D.W.; Harmsen, B.J.; Johnson, P.J.; Newman, C. Increasing frequency of bite wounds with increasing population density in Eurasian badgers, *Meles meles*. *Anim. Behav.* **2004**, *67*, 745–751. [CrossRef]

38. Hunter, R.G.; McEwen, B.S. Stress and anxiety across the lifespan: Structural plasticity and epigenetic regulation. *Epigenomics* **2013**, *5*, 177–194. [CrossRef]

39. Van Lieshout, S.H.J.; Badás, E.P.; Mason, M.W.T.; Newman, C.; Buesching, C.D.; Macdonald, D.W.; Dugdale, H.L. Social effects on age-related and sex-specific immune cell profiles in a wild mammal: Immune cell profiles in a wild mammal. *Biol. Lett.* **2020**, *16*, 20200234. [CrossRef]

40. Dall'Ara, P.; Labriola, C.; Sala, E.; Spada, E.; Magistrelli, S.; Lauzi, S. Prevalence of serum antibody titres against feline panleukopenia, herpesvirus and calicivirus infections in stray cats of Milan, Italy. *Prev. Vet. Med.* **2019**, *167*, 32–38. [CrossRef]

41. Troyer, R.M.; Beatty, J.A.; Stutzman-Rodriguez, K.R.; Carver, S.; Lozano, C.C.; Lee, J.S.; Lappin, M.R.; Riley, S.P.D.; Serieys, L.E.K.; Logan, K.A.; et al. Novel Gammaherpesviruses in North American Domestic Cats, Bobcats, and Pumas: Identification, Prevalence, and Risk Factors. *J. Virol.* **2014**, *88*, 3914–3924. [CrossRef]

42. Sebastiano, M.; Chastel, O.; De Thoisy, B.; Eens, M.; Costantini, D. Oxidative stress favours herpes virus infection in vertebrates: A meta-analysis. *Curr. Zool.* **2016**, *62*, 325–332. [CrossRef]

43. Macdonald, D.W.; Newman, C.; Nouvellet, P.M.; Buesching, C.D. An Analysis of Eurasian Badger (*Meles meles*) Population Dynamics: Implications for Regulatory Mechanisms. *J. Mammal.* **2009**, *90*, 1392–1403. [CrossRef]

44. Dugdale, H.L.; Davison, D.; Baker, S.E.; Ellwood, S.A.; Newman, C.; Buesching, C.D.; Macdonald, D.W. Female teat size is a reliable indicator of annual breeding success in European badgers: Genetic validation. *Mamm. Biol.* **2011**, *76*, 716–721. [CrossRef]

45. Macdonald, D.W.; Feber, R.E.; Newman, C.; Buesching, C.D. Badgers in the rural landscape-conservation paragon or farmland pariah? Lessons from the Wytham Badger Project. *Wildl. Conserv. Farml.* **2015**, *2*, 65–95. [CrossRef]
46. Bright Ross, J.G.B.; Newman, C.; Buesching, C.D.; Macdonald, D.W. What lies beneath? Population dynamics conceal pace-of-life and sex ratio variation, with implications for resilience to environmental change. *Glob. Chang. Biol.* **2020**, [CrossRef]

47. Sugianto, N.A.; Newman, C.; Macdonald, D.W.; Buesching, C.D. Reproductive and Somatic Senescence in the European Badger (*Meles meles*): Evidence from Lifetime Sex-Steroid Profiles. *Zoology* **2020**, *141*, 125803. [CrossRef]

48. Thompson, J.D.; Higgins, D.G.; Gibson, T.J. CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. *Nucleic Acids Res.* **1994**, *22*, 4673–4680. [CrossRef]

49. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. *Mol. Biol. Evol.* **2018**, *35*, 1547–1549. [CrossRef]

50. R-Development-Core-Team R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria. Available online: https://www.R-project.org/ (accessed on 1 October 2019).

51. RStudio Team RStudio: *Integrated Development for R*. RStudio Team RStudio: PBC: Boston, MA, USA, 2015. Available online: http://www.rstudio.com/ (accessed on 1 October 2019).

52. Wilson, E.B. Probable Inference, the Law of Succession, and Statistical Inference. *J. Am. Stat. Assoc.* **1927**, *22*, 209. [CrossRef]

53. Mandrekar, J.N. Receiver operating characteristic curve in diagnostic test assessment. *J. Thorac. Oncol.* **2010**, *5*, 1315–1316. [CrossRef] [PubMed]

54. Kangro, H.O.; Osman, H.K.; Lau, Y.L.; Heath, R.B.; Yeung, C.Y.; Ng, M.H. Seroprevalence of antibodies to human herpesviruses in England and Hong Kong. *J. Med. Virol.* **1994**, *43*, 91–96. [CrossRef] [PubMed]

55. Donahue, J.G. The incidence of herpes zoster. *Arch. Intern. Med.* **1995**, *155*, 1605–1609. [CrossRef] [PubMed]

56. Stowe, R.P.; Kozlova, E.V.; Yetman, D.L.; Walling, D.M.; Goodwin, J.S.; Glaser, R. Chronic herpesvirus reactivation occurs in aging. *Exp. Gerontol.* **2007**, *42*, 563–570. [CrossRef]

57. Gouin, J.P.; Hantsoo, L.; Kiecolt-Glaser, J.K. Immune dysregulation and chronic stress among older adults: A review. *Neuroimmunomodulation* **2008**, *15*, 251–259. [CrossRef]

58. François, S.; Vidick, S.; Sarlet, M.; Desmecht, D.; Drion, P.; Stevenson, P.G.; Vanderplasschen, A.; Gillet, L. Illumination of Murine Gammaherpesvirus-68 Cycle Reveals a Sexual Transmission Route from Females to Males in Laboratory Mice. *PLoS Pathog.* **2013**, *9*, e1003292. [CrossRef]

59. Sugianto, N.A.; Newman, C.; Macdonald, D.W.; Buesching, C.D. Heterochrony of puberty in the European badger (*Meles meles*) can be explained by growth rate and group-size: Evidence for two endocrinological phenotypes. *PLoS ONE* **2019**, *14*, e0203910. [CrossRef]

60. Morey, J.N.; Boggero, I.A.; Scott, A.B.; Segerstrom, S.C. Current directions in stress and human immune function. *Curr. Opin. Psychol.* **2015**, *5*, 13–17. [CrossRef]

61. Barrandeguy, M.; Vissani, A.; Olguin, C.; Becerra, L.; Miño, S.; Pereda, A.; Oriol, J.; Thiry, E. Experimental reactivation of equine herpesvirus-3 following corticosteroid treatment. *Equine Vet. J.* **2008**, *40*, 593–595. [CrossRef]

62. Das Neves, C.G.; Mørk, T.; Thiry, J.; Godfroid, J.; Rimstad, E.; Thiry, E.; Tryland, M. Cervid herpesvirus 2 experimentally reactivated in reindeer can produce generalized viremia and abortion. *Virus Res.* **2009**, *145*, 321–328. [CrossRef]

63. Glaser, R.; Pearl, D.K.; Kiecolt-Glaser, J.K.; Malarkey, W.B. Plasma cortisol levels and reactivation of latent Epstein-Barr virus in response to examination stress. *Psychoneuroendocrinology* **1994**, *19*, 765–772. [CrossRef]

64. Kokko, H.; Ranta, E.; Ruxton, G.; Lundberg, P. Sexually transmitted disease and the evolution of mating systems. *Evolution* **2002**, *56*, 1091–1100. [CrossRef]

65. Ryder, J.J.; Pastok, D.; Hoare, M.J.; Bottery, M.J.; Boots, M.; Knell, R.K.; Atkinson, D.; Hurst, G.D.D. Spatial variation in food supply, mating behavior, and sexually transmitted disease epidemics. *Behav. Ecol.* **2013**, *24*, 723–729. [CrossRef]

66. Munday, P.E.; Pritchard, G.; Harris, J.R.W.; Taylor-Robinson, D. Prevalence of chlamydial infection in promiscuous women. *Sex. Transm. Infect.* **1983**, *59*, 103–104. [CrossRef]

67. White, J.; Richard, M.; Massot, M.; Meylan, S. Cloacal bacterial diversity increases with multiple mates: Evidence of sexual transmission in female common lizards. *PLoS ONE* **2011**, *6*, e22339. [CrossRef]

68. Poiani, A. Sexually Transmitted Diseases: A Possible Cost of Promiscuity in Birds? *Auk* **2000**, *117*, 1061–1065. [CrossRef]

69. Sugianto, N.A.; Dehnhard, M.; Newman, C.; Macdonald, D.W.; Buesching, C.D. A non-invasive method to assess the reproductive status of the European badger (*Meles meles*) from urinary sex-steroid metabolites. *Gen. Comp. Endocrinol.* **2021**, *301*, 113655. [CrossRef]

70. Costantini, D.; Seeber, P.A.; Soilemetzidou, S.E.; Azab, W.; Bohner, J.; Buuveibaatar, B.; Czirják, G.; East, M.L.; Greunz, E.M.; Kaczensky, P.; et al. Physiological costs of infection: Herpesvirus replication is linked to blood oxidative stress in equids. *Sci. Rep.* **2018**, *8*, 1–10. [CrossRef]

71. Urbaniak, C.; Lorenzi, H.; Thissen, J.; Jaing, C.; Crucian, B.; Sams, C.; Pierson, D.; Venkateswaran, K.; Mehta, S. The influence of spaceflight on the astronaut salivary microbiome and the search for a microbiome biomarker for viral reactivation. *Microbiome* **2020**, *8*, 56. [CrossRef]

72. Nordengrahn, A.; Rusvai, M.; Merza, M.; Ekström, J.; Morein, B.; Belák, S. Equine herpesvirus type 2 (EHV-2) as a predisposing factor for Rhodococcus equi pneumonia in foals: Prevention of the bifactorial disease with EHV-2 immunostimulating complexes. *Vet. Microbiol.* **1996**, *51*, 55–68. [CrossRef]
73. Newman, C.; Macdonald, D.W.; Anwar, M.A. Coccidiosis in the European badger, *Meles meles* in Wytham Woods: Infection and consequences for growth and survival. *Parasitology* **2001**, *123*, 133–142. [CrossRef]

74. Albery, G.; Newman, C.; Bright Ross, J.B.; Bansal, S.; Buesching, C. Negative density-dependent parasitism in a group-living carnivore. *Proc. R. Soc. B* **2020**, *287*, 20202655. [CrossRef]