# Infection with a recently discovered gammaherpesvirus variant in european badgers, Meles meles, is associated with higher relative viral loads in blood

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. Infection with a recently discovered gammaherpesvirus variant in european badgers, Meles meles, is associated with higher relative viral loads in blood. Pathogens, 2022, 11 (10), pp.1154-1173. 10.3390/pathogens11101154. hal-04338414

HAL Id: hal-04338414

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

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
<figure><img src="image_2.png" /></figure>

Article

# Infection with a Recently Discovered Gammaherpesvirus Variant in European Badgers, *Meles meles*, Is Associated with Higher Relative Viral Loads in Blood

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

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

check for
updates

**Citation:** Tsai, M.-s.; François, S.; Newman, C.; Macdonald, D.W.; Buesching, C.D. Infection with a Recently Discovered Gammaherpesvirus Variant in European Badgers, *Meles meles*, Is Associated with Higher Relative Viral Loads in Blood. *Pathogens* **2022**, *11*, 1154. https://doi.org/10.3390/pathogens11101154

Academic Editor: Lawrence S. Young

Received: 29 August 2022

Accepted: 3 October 2022

Published: 6 October 2022

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

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

**Copyright:** © 2022 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:** Herpesviruses are ubiquitous pathogens infecting most animals. Although host immunity continually coevolves to combat virulence, viral variants with enhanced transmissibility or virulence occasionally emerge, resulting in disease burdens in host populations. Mustelid gammaherpesvirus 1 (MusGHV-1) is the only herpesvirus species identified thus far in European badgers, *Meles meles*. No MusGHV-1 associated pathomorbidity has been reported, but reactivation of MusGHV-1 in genital tracts is linked to impaired female reproductive success. An analysis of a short sequence from the highly conserved DNA polymerase (DNApol) gene previously identified two variants in a single host population. Here we compared genetic variance in blood samples from 66 known individuals of this same free-ranging badger population using a partial sequence comprising 2874 nucleotides of the DNApol gene, among which we identified 15 nucleotide differences resulting in 5 amino acid differences. Prevalence was 86% (59/66) for the common and 17% (11/66) for the novel variant, with 6% (4/66) of badgers presenting with coinfection. MusGHV-1 variants were distributed unevenly across the population, with individuals infected with the novel genotype clustered in 3 of 25 contiguous social groups. Individuals infected with the novel variant had significantly higher MusGHV-1 viral loads in their blood (p = 0.002) after adjusting for age (juveniles > adults, p < 0.001) and season (summer > spring and autumn, p = 0.005; mixed-effect linear regression), likely indicating higher virulence of the novel variant. Further genome-wide analyses of MusGHV-1 host resistance genes and host phenotypic variations are required to clarify the drivers and sequelae of this new MusGHV-1 variant.

**Keywords:** red queen hypothesis; coevolution; sexually transmitted infection; one health; genetic epidemiology

## 1. Introduction

Under the concept of 'One Health' (a framework recognising the interconnection between human, animal and environmental health), emerging pathogen strain monitoring and pathotype identification are considered essential for disease surveillance and prevention, not only in humans and domestic animals but also in wildlife [1,2]. The majority (60%) of emerging infectious diseases are of animal origin, predominantly (72%) from wildlife [1]. Due to a lack of fidelity during replication, viruses have a higher mutation rate than cellular organisms. This results in higher evolutionary plasticity, which favours viral adaptation to changing environments (i.e., novel hosts, new host immune strategies, different medical treatments, etc.) and the emergence of new strains [3,4]. The resulting arms race between pathogen virulence and host resistance continually shapes the evolution of competitive coevolving genes or organisms through repeated cycles of adaptation

<footnote n="1">¹ Wildlife Conservation Research Unit, Recanati-Kaplan Centre, Department of Biology, University of Oxford, Tubney, Oxfordshire OX13 5QL, UK</footnote>

<footnote n="2">² Department of Biology, University of Oxford, Peter Medawar Building, South Park Road, Oxford OX1 3SY, UK</footnote>

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

<footnote>* Correspondence: cindy150051@gmail.com</footnote>
and counter-adaptation, a process formalised under the Red Queen hypothesis (RQH) [5]. When applied to host-pathogen coevolution, the RQH predicts that, under natural conditions, changes in host allele frequencies within a population cause pathogens to adjust allele frequencies to match those of the host and vice versa [6–8]. This can drive intrahost speciation (i.e., several closely related virus strains or species can evolve independently within the same host species) [9,10]. A recent study shows that this so-called ‘duplication’ has occurred more frequently than expected in the evolution of herpesviruses (previously believed to be dominated by cospeciation with their host species), especially among viruses belonging to the *Gammaherpesvirinae* subfamily [11].

European badgers (*Meles meles*) are promiscuous, polygynandrous mustelids [12]. Their social system is highly flexible [13], ranging from solitary or pair-living [14] to the formation of extended family groups comprising upwards of 30 individuals [15] Mustelid gammaherpesvirus 1 (MusGHV-1) is the only gammaherpesvirus known to infect badgers [16]. Like other gammaherpesviruses, after primary infection, MusGHV-1 may establish latency in white blood cells [17]. Screening badger blood samples revealed a 98–100% prevalence of MusGHV-1 in wild British and Irish populations [18,19]. Reactivation from latency causes viral replication and shedding in the mucosa in oral, nasal, and/or genital tracts. Depending on the season and population [20], the prevalence of genital shedding of MusGHV-1 ranges from 55% to 62% [21,22]. No specific signs or pathomorbidity have been linked to MusGHV-1 reactivation in badgers. However, viral shedding is associated with lower female reproductive success and lesions in the genital tract of both sexes [21]. In our previous analysis of partial DNA polymerase (DNApol) gene sequences [20], we identified two distinct variants through the genotyping of MusGHV-1 in a single closed wild badger population located in Wytham Woods, Oxfordshire (for details of the study population, sees Buesching et al., 2010 [23]).

Since the herpesviruses DNApol gene encodes exonuclease and polymerase proteins essential for virus replication and viral genome proofreading, its quantified expression could then be used as a proxy for determining virulence [24]. The DNApol gene is thus one of the most conserved genes in herpesviral genomes and likely to be under strong purifying selection with a low substitution rate. Nevertheless, any mutations found in these domains are typically associated with lower replication fidelity (see review: [25]), contributing to higher mutation rates. If some of these mutations are beneficial to virus replication efficiency, they can, in turn, result in faster evolution rates [26]. Consequently, herpesvirus strains belonging to the same species that show differences in their DNApol genes can result in different disease expressions and outcomes in their hosts. For example, a single nucleotide difference in the equine herpesvirus 1 (EHV-1) polymerase gene on amino acid position 752 (N to D) is associated with neurological disease outbreaks [27–29] and higher viraemia levels in infected horses compared with those infected with the EHV-1 N<sub>752</sub> strain, suggesting a higher virus replication efficiency [30]. Differences in host susceptibility to the two variants have not yet been explored, but simultaneous infections with both strains have been reported [31,32].

Here we investigated genetic, epidemiological, and phenotypic differences between the two previously reported MusGHV-1 variants in a natural host population [20]. Specifically, we aimed to:

(1) acquire longer sequences of the MusGHV-1 DNA polymerase gene to investigate for additional mutations/genetic differences between the two variants;

(2) establish the relative prevalence of the two MusGHV-1 variants within the population, investigate the rate of coinfection, and whether potential differences in variant prevalence are associated with demographic (i.e., host sex and age) and/or socio-geographic (i.e., badger sett and social group) patterns;

(3) analyse potential differences in viral loads in blood samples between individuals infected with each or both MusGHV-1 variants.
## 2. Materials and Methods

### 2.1. Animal Sampling

We captured free-ranging badgers in Wytham Woods, Oxfordshire, United Kingdom, (51°46'26" N, 1°19'19" W; for details of the study site, see [33]; for details of badger trapping and handling protocols see [20,34]) in spring (21st May to 2nd June), summer (3rd to 15th September) and autumn (13th to 24th November) 2018. In summary, individual badgers were identified by their serial tattoo numbers, given at first capture as part of a long-term badger research programme [15,35]. After sedation with ketamine hydrochloride [36], oral, nasal, genital, and rectal swab samples were collected using wooden shafted cotton tips (Dynarex). Blood samples (not exceeding 2 mL) were collected by jugular venipuncture in heparin-coated tubes and stored immediately at -20 °C. All protocols and procedures were approved by the Zoology Department's (University of Oxford) Animal Welfare and Ethical Review Body and were conducted under the Animals (Scientific Procedures) Act, 1986 (PPL: 30/3379).

### 2.2. Genotyping of MusGHV-1 DNA

DNA was extracted and purified from swab samples using DNeasy Blood and Tissue Kits (Qiagen), and from 200 µL whole blood samples using Monarch Genomic DNA Purification Kits (New England Biolabs Inc., Ipswich, MA, USA) following manufacturers' instructions (for full methodological details see Tsai et al. [20]). MusGHV-1 DNA was amplified from each sample by PCR using primers (pol3-F and pol3-R, Table 1), specifically targeting the MusGHV-1 DNA polymerase gene. After checking amplification results on 1.5% agarose gels following electrophoresis, we sent 66 amplified PCR products for genotyping using Sanger sequencing (Zoology Sequencing, University of Oxford, Oxford, UK). Variant allocation was determined using ClustalW alignment (with the respective default parameters) in MEGA11 (Molecular Evolutionary Genetics Analysis version 11.0.8, Tamura, Stecher, and Kumar) according to nine substitutions across 691 nucleotides identified between the two previously reported MusGHV-1 sequences (accession numbers MT332100 and MT332101) in the same study population: at position 363 (C to T), 378 (T to C), 383 (A to G), 480 (T to C), 486 (T to C), 537 (A to G), 540 (G to A), 546 (G to A), and 558 (C to A) [20]. Coinfection was reported if the aforementioned signals for both nucleotides were present in a given sequencing product.

<table><caption>**Table 1.** Primers used for MusGHV-1 DNA polymerase gene PCR amplification.</caption><thead><tr><th></th><th>Start Position</th><th>End Position</th><th>Primer Sequence (5'-3')</th><th>Product Length</th></tr></thead><tbody><tr><td>pol5-F</td><td>3</td><td>22</td><td>GGCAGGGAATTTTTATAACC</td><td></td></tr><tr><td>pol5-R</td><td>681</td><td>701</td><td>CCACCCAAAAGTAGAAAATCC</td><td>699</td></tr><tr><td>pol6-F</td><td>345</td><td>364</td><td>CCCCTCTGGAACTGTGCTAA</td><td></td></tr><tr><td>pol6-R</td><td>962</td><td>981</td><td>AGGGTCACATGTCCCCAAAT</td><td>637</td></tr><tr><td>pol7-F</td><td>863</td><td>846</td><td>ATGTCTGGGGGAAAATGG</td><td></td></tr><tr><td>pol7-R</td><td>1329</td><td>1348</td><td>GACCTCCTATGCACTGCTTG</td><td>486</td></tr><tr><td>pol10-F</td><td>1185</td><td>1204</td><td>TGAAGTTCACACACCCCAGA</td><td></td></tr><tr><td>pol10-R</td><td>1524</td><td>1545</td><td>TCCATCGGTCAGCACTCTC</td><td>361</td></tr><tr><td>pol3-F</td><td>1328</td><td>1347</td><td>CCAAGCAGTGCATAGGAGGT</td><td></td></tr><tr><td>pol3-R</td><td>2072</td><td>2098</td><td>TGGACTTCTCCAACATGCGTCGCCCTT</td><td>771</td></tr><tr><td>pol11-F</td><td>1684</td><td>1703</td><td>CCGATCTTGGTGGTTGATTT</td><td></td></tr><tr><td>pol11-R</td><td>2291</td><td>2310</td><td>CTTAATTGGCTCCTCGAACA</td><td>627</td></tr><tr><td>pol12-F</td><td>2015</td><td>2034</td><td>CAGGTGTGTCCTCGGGTATT</td><td></td></tr><tr><td>pol12-R</td><td>2556</td><td>2575</td><td>TCACTTTGAAAAGTGGAAGTGG</td><td>561</td></tr><tr><td>pol9-F</td><td>2393</td><td>2412</td><td>TGATGAAGGGAGTGGATCTC</td><td></td></tr><tr><td>pol9-R</td><td>2936</td><td>2989</td><td>TCACAGCTTTGTCTGCACTG</td><td>597</td></tr></tbody></table>

### 2.3. MusGHV-1 DNA Polymerase Gene Mapping and Phylogenetic Analysis

One sample of each variant was chosen for extended DNAPol gene genotyping using primer pairs designed specifically for this study (Table 1). The same PCR procedures
and conditions were used as described above, and successfully amplified PCR products were sent for Sanger sequencing (Zoology Sequencing, University of Oxford). Sequences were then mapped to the reference gene (DNA polymerase gene of MusGHV-1, Accession number: AF275657), assembled, translated to amino acids, and aligned using the ClustalW method (with set default parameters) [37] to 7 representatives of published gammaherpesvirus sequences in GenBank (having either complete or > 1000 bp of the DNA polymerase gene), using the software Unipro UGENE (version 37, Okonechnikov, Golosova, Fursov, the UGENE team, Novosibirsk, Russia) [38]. For each isolate, a purifying selection test using the Nei-Gojobori method was employed to determine whether the number of nonsynonymous substitutions per nonsynonymous site (dNS) was fewer than the number of synonymous substitutions per synonymous site (dS) [39]. Nucleotide sequences were aligned using ClustalW with default parameters. A maximum likelihood (ML) phylogenetic tree was constructed under the GTR + I + gamma nucleotide model using 100 bootstrap replicates using MAFFT version 7.450 [40]. The ML tree was mid-point rooted.

## 2.4. Geographic Distribution of MusGHV-1 Strains and Social Group-Specific Prevalence

The social group affiliation of each individual was determined using previously established residency rules [41,42] (Figure S1), with social group ranges established through bait-marking surveys [34,42]. Social group-specific MusGHV-1 variant prevalence was calculated by dividing the number of total positive badgers by the total of tested badgers resident in each social group. All maps were drawn using QGIS (version 3.14.15, QGIS Development Team).

## 2.5. Blood Viral Load Quantification

We selected blood samples from individuals with known MusGHV-1 variant infections. Some individuals (see Table 2) were sampled in more than one trapping season. We analysed each sample's relative MusGHV-1 genome copy numbers using StepOne-Plus PCR Systems (Applied Biosystems, ABI) to obtain Ct values. We used primers designed by Sin et al. [18] targeting the MusGHV-1 DNA polymerase gene (forward: 5'-GGAGAGTGCTGACCGATGGA; reverse: 5'-AAAAGCCTGGAATTGGATCAATAA, product length: 150 bp). We prepared 20 µL of reaction mix containing 0.5 mM of each primer, 10 µL of Luna Universal qPCR Master Mix, 3 µL of RNase-free water and 5 µL of DNA template. Amplification conditions were as follows: 95 °C for 60 s (initial denaturation), 45 cycles of 95 °C (denaturation) for 15 s, and 60 °C for 30 s (extension). StepOne Software (version 2.3) obtained Ct values of each sample automatically.

<table><caption>**Table 2.** Blood samples used for relative viral load analysis.</caption><thead><tr><th>Tattoo Number</th><th>Sex</th><th>Age</th><th>Spring</th><th>Summer</th><th>Autumn</th></tr></thead><tbody><tr><td>1012</td><td>Male</td><td>13</td><td></td><td>v</td><td></td></tr><tr><td>1045x</td><td>Female</td><td>13</td><td>v</td><td></td><td>v</td></tr><tr><td>1232</td><td>Female</td><td>10</td><td>v</td><td></td><td>v</td></tr><tr><td>1234x</td><td>Female</td><td>10</td><td></td><td>v</td><td>v</td></tr><tr><td>1283</td><td>Male</td><td>9</td><td>v</td><td>v</td><td></td></tr><tr><td>1330</td><td>Male</td><td>8</td><td></td><td>v</td><td></td></tr><tr><td>1379</td><td>Male</td><td>&gt;8*</td><td>v</td><td></td><td>v</td></tr><tr><td>1435</td><td>Female</td><td>6</td><td></td><td>v</td><td>v</td></tr><tr><td>1478</td><td>Female</td><td>5</td><td></td><td>v</td><td></td></tr><tr><td>1487</td><td>Male</td><td>7</td><td>v</td><td></td><td>v</td></tr><tr><td>1520</td><td>Female</td><td>4</td><td>v</td><td>v</td><td></td></tr><tr><td>1622</td><td>Male</td><td>3</td><td>v</td><td>v</td><td></td></tr><tr><td>1694</td><td>Female</td><td>2</td><td></td><td>v</td><td>v</td></tr><tr><td>1699</td><td>Male</td><td>1</td><td></td><td>v</td><td>v</td></tr><tr><td>1735</td><td>Male</td><td>0</td><td>v</td><td></td><td></td></tr><tr><td>1746</td><td>Male</td><td>0</td><td></td><td>v</td><td></td></tr><tr><td>1749</td><td>Female</td><td>0</td><td>v</td><td></td><td></td></tr></tbody></table>
<table><caption>**Table 2.** *Cont.*</caption><thead><tr><th>Tattoo Number</th><th>Sex</th><th>Age</th><th>Spring</th><th>Summer</th><th>Autumn</th></tr></thead><tbody><tr><td>1750</td><td>Female</td><td>0</td><td></td><td>v</td><td>v</td></tr><tr><td>1751</td><td>Female</td><td>0</td><td>v</td><td>v</td><td>v</td></tr><tr><td>1753</td><td>Female</td><td>0</td><td>v</td><td></td><td></td></tr><tr><td>1754</td><td>Male</td><td>0</td><td>v</td><td></td><td></td></tr><tr><td>1755</td><td>Male</td><td>0</td><td>v</td><td>v</td><td>v</td></tr><tr><td>1778</td><td>Male</td><td>0</td><td></td><td>v</td><td></td></tr><tr><td>1785</td><td>Female</td><td>0</td><td></td><td></td><td>v</td></tr></tbody></table>

## 2.6. Statistical Analyses

All statistical analyses were performed using R studio (version 1.2.1335, RStudio Team, Boston, MA, USA) [44,45]. Fisher's exact tests were used to test if infections with the MusGHV-1 variant differed with sex, age group (adult: ≥2 years old, and juvenile: <2 years old) or social group (n = 25). Differences in viral loads between variants were analysed with a mixed-effect linear regression analysis using the R package lmerTest. We used Ct values as the response variable and variants (common, coinfection and novel), age group and sampling seasons (spring, summer and autumn) as fixed effects. Tattoo ID was included as a random effect, as some individuals were sampled more than once.

# 3. Results

## 3.1. Substitutions in the DNApol Gene of the Two MusGHV-1 Variants

We sequenced 66 MusGHV-1 PCR products amplified from the partial DNApol gene. All sequences were trimmed to 693 base pairs and were confirmed to be MusGHV-1, showing between 98.7% and 100% identity with the published MusGHV-1 sequence isolated from a badger in Cornwall, England (Accession number: AF275657).

The extended genotyping of the DNApol gene resulted in the acquisition of a total sequence length of 2874 nucleotides for each of the two MusGHV-1 variants. This enabled us to test for additional variations between the two variants. In total, we found 15 nucleotides and five amino acid differences in the DNA polymerase domain between the two variants, which we classified here as the 'common variant' and the 'novel variant', as the former was almost identical (99.9% nucleotide identity) to the reference sequence. In contrast, the novel variant only shared 99.4% identity with the reference sequence. Both MusGHV-1 variants had three nucleotides and two amino acid substitutions in common with the reference sequence, all located at the exonuclease domain (details of these substitutions are listed in Table 3 and Figure 1). There was support for strong purifying selection between the two variants in the Wytham badger population (dNS < dS, *p* value = 0.004, dS/dNS = 2.66).

<table><caption>**Table 3.** Amino acid substitutions (nonsynonymous differences) between the two variants in the MusGHV-1 DNA polymerase gene.</caption><thead><tr><th>Position</th><th>Reference (AF275657 <sup>a</sup>)</th><th>Common Variant</th><th>Novel Variant</th><th>Conserved Protein Domain Family <sup>b</sup></th><th>In Conserved Region <sup>c</sup></th><th>In Conserved Region within Host Order <sup>d</sup></th></tr></thead><tbody><tr><td>253</td><td>P</td><td>L</td><td>L</td><td>3'-5' exonuclease Exo I</td><td>Yes</td><td>Yes</td></tr><tr><td>274</td><td>N</td><td>Y</td><td>Y</td><td>3'-5' exonuclease Exo I</td><td>Yes</td><td>Yes</td></tr><tr><td>591</td><td>K</td><td>K</td><td>R</td><td>Polymerase</td><td>No</td><td>No</td></tr><tr><td>649</td><td>N</td><td>N</td><td>K</td><td>Polymerase</td><td>No</td><td>Yes</td></tr><tr><td>819</td><td>E</td><td>E</td><td>G</td><td>Polymerase</td><td>No</td><td>No</td></tr><tr><td>872</td><td>T</td><td>T</td><td>N</td><td>Polymerase</td><td>No</td><td>No</td></tr><tr><td>889</td><td>S</td><td>S</td><td>L</td><td>Polymerase</td><td>No</td><td>Yes</td></tr></tbody></table>

<footnote>* Age not determined as this individual was first captured as an adult; according to its tooth wear [43], it was estimated as older than 8 years of age.</footnote>

<footnote>a: GenBank accession number. b: Data from the NCBI Conserved Domain Database (CDD) [46]. c: Determined by whether the substitution is located in previously identified conserved motif of viral DNA polymerase genes [47]. d: Determined by whether the substitution is located in a position with a higher conservation degree within herpesvirus species confirmed in the order Carnivora (Figure 1).</footnote>
<figure><img src="image_5.png" /><figcaption>**Figure 1.** DNA polymerase gene alignment containing three MusGHV-1 variants and seven other gammaherpesviruses belonging to the genus *Percavirus*. Red rectangles show amino acid substitutions between the three MusGHV-1 variants (numbers 1 to 3 on the list).</figcaption></figure>

Among these 66 individuals, the common variant was detected in 59 (89.4%), including genital swabs (n = 49), oral swabs (n = 5), rectal swabs (n = 2), and blood samples (n = 5). The novel variant was detected in 11 (16.7%) individuals, including in genital swabs (n = 9) and blood samples (n = 5). Coinfection was confirmed in four individuals (6.1%) in three blood samples and one rectal swab (Supplementary Material Table S1).

## 3.2. Phylogenetic Relationships between the Two MusGHV-1 Variants and Other Gammaherpesvirus Species

All Carnivora-infecting gammaherpesviruses are grouped in a monophyletic clade, including the MusGHV-1 virus, which likely belongs to the genus *Percavirus* [48]. A phylogenic analysis using the DNApol gene sequence showed that the novel and the common variant share a common ancestor (Figure 2). The Cornish variant (NC_038266) was derived from the common variant more recently, making the more recently discovered variant the actual ancestral variant.

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Maximum likelihood phylogenetic tree based on a portion of the polymerase gene of 10 gammaherpesvirus sequences. The alignment of sequences trimmed at 693 nucleotides was produced using MAFFT v7.450 using the G-INS-i algorithm. The tree is mid-point rooted. Bootstrap</figcaption></figure>
values (100 replicates) in excess of 30% are indicated at each node. The scale bar corresponds to nucleotide substitutions per site. The MusGHV-1 variants discussed in this paper (i.e., common and novel variants) are coloured in red.

### 3.3. *Novel MusGHV-1 Variant Infections Were Clustered within Just Three Social Groups*

MusGHV-1 variant frequencies did not differ significantly between males and females (Fisher's exact test, *p* = 0.19) or between age groups (*p* = 0.19) but varied somewhat between social groups (*p* = 0.052) (Table S3). The novel variant (and coinfection with both variants) was confined to badgers from three neighbouring social groups, while the common strain was evenly distributed across the Woods (Figure 3). To investigate further, we screened 14 DNA samples from blood collected from badgers resident in these three social groups during 2009 and 2010. Two individuals were infected with the novel variant, one was infected with both variants and the rest were infected with the common variant. This retrospective scrutiny indicated that the novel variant originated in the CHO social group before spreading to the CH and RC groups (Table S4).

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Map of genital tract MusGHV-1 occurrence and MusGHV-1 variants distribution in badger social groups in Wytham Woods in 2018. The name of each social group is highlighted in yellow. The areas with grey colour indicate sample sizes of less than 3.</figcaption></figure>

### 3.4. *The MusGHV-1 Novel Variant Was Associated with Higher Viral Load in the Blood Stream*

The Ct values from real-time PCR on 40 blood samples collected from 24 individuals in spring, summer and autumn 2018 (Table 2) showed that the MusGHV-1 relative viral load was highest (i.e., low Ct value) among badger cubs (age class) then decreased thereafter in yearlings (age class 1) (Figure 4). Furthermore, juveniles had higher relative viral loads than adults (*p* < 0.001), and samples taken in summer had higher loads than those in spring or autumn (*p* = 0.006) (Table 4). The mixed-effect linear regression analysis results suggest
that badgers infected with the novel variant generally had a higher MusGHV-1 viral load than those infected with the common variant (*p* = 0.004).

<figure><img src="image_8.png" /><figcaption>**Figure 4.** Ct values representing relative blood MusGHV-1 viral load of 39 blood samples collected from 23 individuals in different seasons (see Table S2 for details). The relative viral load was higher in juveniles (less than 2 years old) than in any other age class (linear mixed-effect model, *t*-test, *p* < 0.001, Table 4).</figcaption></figure>

<table><caption>**Table 4.** Results of linear mixed-effect model fit by REML (restricted maximum likelihood). *t*-tests using Satterthwaite's method (formula= Ct ~ Variant + Age Group + Season + (1 | Tattoo)). Number of observations: 39, groups (by tattoo ID: 23).</caption><thead><tr><th>Variable</th><th>Estimate</th><th>Std. Error</th><th>df</th><th>t Value</th><th>p Value</th></tr></thead><tbody><tr><td>(Intercept)</td><td>34.855</td><td>0.464</td><td>28.794</td><td>75.072</td><td></td></tr><tr><td>Genotype</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Common</td><td>(Reference)</td><td></td><td></td><td></td><td></td></tr><tr><td>Novel</td><td>-1.815</td><td>0.547</td><td>19.447</td><td>-3.316</td><td>0.004</td></tr><tr><td>Coinfection</td><td>-0.346</td><td>0.834</td><td>32.885</td><td>-0.415</td><td>0.681</td></tr><tr><td>Age Group</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Adult</td><td>(Reference)</td><td></td><td></td><td></td><td></td></tr><tr><td>Juvenile</td><td>-3.030</td><td>0.543</td><td>18.615</td><td>-5.583</td><td>&lt;0.001</td></tr><tr><td>Season</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Spring</td><td>(Reference)</td><td></td><td></td><td></td><td></td></tr><tr><td>Summer</td><td>-1.379</td><td>0.453</td><td>23.179</td><td>-3.045</td><td>0.006</td></tr><tr><td>Autumn</td><td>0.012</td><td>0.416</td><td>17.381</td><td>0.028</td><td>0.978</td></tr></tbody></table>

## 4. Discussion

Understanding epidemiological traits can aid the interpretation of viral load data across wildlife species. Despite this study's relatively small sample size, epidemiological variations were apparent in viral loads. Higher relative viral loads of MusGHV-1 in juveniles corroborated a previous study using animals from the same population [18], where primary herpes infection with a higher viral load typically occurs early in life [49]. Higher viral loads in summer were correlated with a higher prevalence of genital tract viral DNA reactivation [20]. Most intriguingly, we discovered that the relative blood viral load was also dependent on the MusGHV-1 variant. This could indicate that the two MusGHV-1 variants found in Wytham have not only genotypic but also phenotypic differences. A higher viral load could cause a higher pathogen burden and additional stress on the host immune
system, which can lead to higher susceptibility to other pathogens [18]. Accordingly, we previously found that badgers infected with the novel variant, which had higher relative viral loads in this study, had a higher risk of infection with *Clostridium perfringens*, a zoonotic bacterial pathogen that can cause severe enteritis [50].

The genetic diversity of herpesviruses strains found within a single host species is typically very low, even between host populations separated by geographical barriers. For example, the gammaherpesvirus (LcaGHV-1) variants found in two geographically distinct Canadian lynx (*Lynx canadensis*) populations (from Maine, USA, and Newfoundland, Canada) are identical apart from just two synonymous substitutions across the 3.4 kb nucleotide sequences of its partial gB protein and DNApol gene [51]. Considering the generally low mutation rate of herpesviruses and their strong host cospeciation traits [50], the coexistence of two closely related gammaherpesvirus variants in the same population is unusual. A molecular analysis of the DNApol gene of twelve EHV-2 strains circulating in horse herds in Iceland that had been isolated for more than 1000 years, still showed 99% to 100% similarity across a region of 483 base pairs to seven EHV-2 strains from other European countries [52]. Furthermore, country-wide screening for MusGHV-1 in Irish badgers in 2019 showed that 691 base pair sequences from 23 individuals were 100% identical to the reference MusGHV-1 DNApol sequence [21]; all of which demonstrates the stability of species-specific herpesvirus DNApol genes.

Establishing whether the MusGHV-1 variant in the Wytham Woods population originated from introduction or evolution is challenging without a genome-wide analysis of mutations and sequencing results from other badger populations across the UK and nearby countries. Potentially, MusGHV-1 may have coevolved with badgers from a different population (likely subject to greater isolation by distance than the Cornish and Irish populations, e.g., continental Europe [53]) and then have been introduced into the Wytham population more recently, but at least before 2009 [54]. However, there are no historical records of badgers either colonising or being introduced into Wytham from other subpopulations, and badger genetic diversity in this population is low [55]; although any such introduced badger could, plausibly, have caused contagion with the novel variant and died without leaving any descendants. Ultimately, it seems unlikely that the novel variant originated from an ex situ source.

A phylogenetic analysis revealed that the node from/at the Cornwall variant branched from the ancestral MusGHV-1 root occurred earlier than the node at which either Wytham population branched which separated and evolved independently from an earlier time point. Although the common variant is dominant in the Wytham and the Irish populations, this does not infer that the novel variant emerged (branched) later than the common variant. It is also possible that the common variant became dominant by selectively outcompeting the novel variant, even if that novel strain branched earlier. For instance, if the higher virulence of the novel variant (rising from replicative advantage) caused higher host mortality prior to some immune adaptation, or triggered a more robust host immune response to eliminate the pathogen [56], then a less virulent variant (common) would spread among surviving badgers and become dominant. Some explanation based on balancing or fluctuating selection between specific variant advantages [57,58], or among host immune alleles [59], could account for occasional coinfection and why infections of the novel variants were limited to only a few social groups due to incomplete over-dominance.

The hypothesis that virulence may be a trade-off against transmission is debated [56,60]. Although a higher viral load can indicate more efficient within-host replication or higher virulence [60], various vertebrate and invertebrate models show that selection for an intermediate replication level can maximise pathogen fitness [61,62]. This may prevent the novel MusGHV-1 variant found in this study from being completely wiped out from the badger population. Although the high replication efficiency of the novel variant may result in it being outcompeted by the common variant due to host immune activities, this advantage may aid the novel variant to transmit and seed in the population through virus-shedding individuals. The lifelong infection trait of the herpesvirus may allow the novel variant to
persist in a latency stage in the host at subclinical levels, possibly beyond the detection
limit of PCR.

An alternative but not mutually exclusive explanation is that differences in host-pathogen resistance determined by immune-related genes, such as MHC genes, may be a causal factor for new strain emergence and the coexistence of multiple strains within a population. In badgers, Sin et al. (2014) found that individuals with specific MHC variants had different blood viral loads of MusGHV-1. Similarly, in horses, the EHV-1-induced downregulation of cell surface molecules expressed by MHC I is allele-dependent [63], and the extent to which MHC I is downregulated is EHV-1 strain-dependent [64]. In humans, herpesvirus infection has been linked to dementia [65], as some human leucocyte antigen (HLA; equivalent to MHC) alleles that protect against dementia have a high binding ability to a variety of human herpesvirus strains, suggesting that this link may be dependent on host variant [66]. Moreover, the subtypes of human herpesvirus 8 (HHV-8, or Kaposi's sarcoma-associated herpesvirus, KSHV) are clustered by ethnic group, and infection with different subtypes affects the risk of Kaposi's sarcoma occurring [67]. However, Wytham's badgers engage extensively in extra-group mating. Therefore, if a specific MHC variant is more susceptible to the novel variant reactivation, it should spread quickly across social groups; however, this is contrary to our observations where the novel variant expanded from one to three social groups over nine badger cohorts (from 2009 to 2018). This would require a recessive MHC gene circulating within these social groups only expressed through inbreeding [42]. Currently, we lack the data to corroborate or refute this hypothesis.

**5. Conclusions**

In this study, we report a newly discovered variant of MusGHV-1 in badgers, indicating that a highly host-specialised pathogen maintains a certain degree of genetic diversity at the host-population level. This could illustrate an ongoing coevolution between badgers and MusGHV-1. Furthermore, we also found that variant variation in MusGHV-1 results in different epidemiological and spatial patterns of infection. By better understanding viral genetics and pathotypes and how these interact with coevolving host immunity, we can accurately inform conservation management and initiatives, such as One Health, tasked with disease surveillance and prevention in wildlife.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/pathogens11101154/s1, Figure S1: Location of badger setts and range of social groups in Wytham Wood; Table S1: List of samples with MusGHV-1 DNA variants and host details of tattoo ID, sex, age, natal group, social group and sett visited; Table S2: Ct values of MusGHV-1 genome by variant in 40 blood samples from badgers of previously known identity (Tattoo ID) and age class; Table S3: Counts of individuals infected with the common or novel MusGHV-1 variant categorised by sex, age class and social groups; Table S4: List of samples tested for MusGHV-1 in blood and genotyping.

**Author Contributions:** Conceptualisation, M.-s.T. and C.D.B.; methodology, M.-s.T., S.F. and C.D.B.; formal analysis, M.-s.T. and S.F.; investigation, M.-s.T., C.D.B. and C.N.; resources, D.W.M.; data curation, M.-s.T.; writing—original draft preparation, M.-s.T.; writing—review and editing, S.F., C.D.B., C.N. and D.W.M.; visualisation, M.-s.T.; supervision, S.F. and D.W.M.; project administration, C.N.; funding acquisition, D.W.M. 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's 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 Licence (PIL) holders under Home Office licence (current PPL 30/3379, Animals (Scientific Procedures) Act 1986).

**Informed Consent Statement:** Not applicable.
**Data Availability Statement:** Not applicable.

**Acknowledgments:** The authors would like to thank Nadine Sugianto, Sil van Lieshout, Tanesha Allen, and Julius Bright-Ross for assistance with sample collection and fieldwork. M.T. would like to thank the Ministry of Education in Taiwan and Lady Margaret Hall, University of Oxford, for doctoral scholarship support.

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

## References

1. Xiao, X.; Newman, C.; Buesching, C.D.; Macdonald, D.W.; Zhou, Z.M. Animal sales from Wuhan wet markets immediately prior to the COVID-19 pandemic. *Sci. Rep.* **2021**, *11*, 11898. [CrossRef] [PubMed]

2. Worobey, M.; Levy, J.I.; Serrano, L.M.; Crits-Christoph, A.; Pekar, J.E.; Goldstein, S.A.; Rasmussen, A.L.; Kraemer, M.U.G.; Newman, C.; Koopmans, M.P.G.; et al. The Huanan Seafood Wholesale Market in Wuhan was the early epicenter of the COVID-19 pandemic. *Science* **2022**, *959*, abp8715. [CrossRef] [PubMed]

3. Biek, R.; Pybus, O.G.; Lloyd-Smith, J.O.; Didelot, X. Measurably evolving pathogens in the genomic era. *Trends Ecol. Evol.* **2015**, *30*, 306–313. [CrossRef] [PubMed]

4. Duffy, S.; Shackelton, L.A.; Holmes, E.C. Rates of evolutionary change in viruses: Patterns and determinants. *Nat. Rev. Genet.* **2008**, *9*, 267–276. [CrossRef] [PubMed]

5. Scheinberg, L. The Red Queen. *Arch. Neurol.* **1983**, *40*, 189. [CrossRef]

6. Carius, H.J.; Little, T.J.; Ebert, D. Genetic variation in a host-parasite association: Potential for coevolution and frequency-dependent selection. *Evolution* **2001**, *55*, 1136–1145. [CrossRef]

7. Lively, C.M.; Dybdahl, M.F. Parasite adaptation to locally common host genotypes. *Nature* **2000**, *405*, 679–681. [CrossRef]

8. Salathé, M.; Kouyos, R.D.; Bonhoeffer, S. The state of affairs in the kingdom of the Red Queen. *Trends Ecol. Evol.* **2008**, *23*, 439–445. [CrossRef]

9. Escalera-Zamudio, M.; Rojas-Anaya, E.; Kolokotronis, S.-O.; Taboada, B.; Loza-Rubio, E.; Méndez-Ojeda, M.L.; Arias, C.F.; Osterrieder, N.; Greenwood, A.D. Bats, Primates, and the Evolutionary Origins and Diversification of Mammalian Gammaherpesviruses. *mBio* **2016**, *7*, e01425-16. [CrossRef]

10. Davison, A.J. Evolution of the herpesviruses. *Vet. Microbiol.* **2002**, *86*, 69–88. [CrossRef]

11. Brito, A.F.; Baele, G.; Nahata, K.D.; Grubaugh, N.D.; Pinney, J.W. Intrahost speciations and host switches played an important role in the evolution of herpesviruses. *Virus Evol.* **2021**, *7*, veab025. [CrossRef] [PubMed]

12. Carpenter, P.J.; Pope, L.C.; Greig, C.; Dawson, D.A.; Rogers, L.M.; Erven, K.; Wilson, G.J.; Delahay, R.J.; Cheeseman, C.L.; Burke, T. Mating system of the Eurasian badger, Meles meles, in a high density population. *Mol. Ecol.* **2005**, *14*, 273–284. [CrossRef] [PubMed]

13. Johnson, D.D.P.; Macdonald, D.W.; Dickman, A.J. An analysis and review of models of the sociobiology of the Mustelidae. *Mamm. Rev.* **2000**, *30*, 171–196. [CrossRef]

14. Woodroffe, R.; Macdonald, D.W. Badger sociality—Models of spatial grouping. *Mamm. Predators* **1993**, *65*, 145–169.

15. Macdonald, D.W.; Newman, C. Population dynamics of badgers (*Meles meles*) in Oxfordshire, UK: Numbers, density and cohort life histories, and a possible role of climate change in population growth. *J. Zool.* **2002**, *256*, 121–138. [CrossRef]

16. Banks, M.; King, D.P.; Daniells, C.; Stagg, D.A.; Gavier-Widen, D. Partial characterization of a novel gammaherpesvirus isolated from a European badger (*Meles meles*). *J. Gen. Virol.* **2002**, *83*, 1325–1330. [CrossRef]

17. 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]

18. 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]

19. 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]

20. 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. [CrossRef]

21. Tsai, M.; 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]

22. 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.* **2018**, *54*, 133–137. [CrossRef] [PubMed]

23. Buesching, C.D.; Clarke, J.R.; Ellwood, S.A.; King, C.; Newman, C.; Macdonald, D.W. The Mammals of Wytham Woods. In *Wytham Woods: Oxford's Ecological Laboratory*; Oxford University Press: Oxford, UK, 2011; ISBN 9780191810039.

24. Dourmishev, L.A.; Dourmishev, A.L.; Palmeri, D.; Schwartz, R.A.; Lukac, D.M. Molecular Genetics of Kaposi's Sarcoma-Associated Herpesvirus (Human Herpesvirus 8) Epidemiology and Pathogenesis. *Microbiol. Mol. Biol. Rev.* **2003**, *67*, 175–212. [CrossRef] [PubMed]
25. Trimpert, J.; Osterrieder, N. Herpesvirus DNA Polymerase Mutants—How Important Is Faithful Genome Replication? *Curr. Clin. Microbiol. Rep.* **2019**, *6*, 240–248. [CrossRef]

26. Tian, W.; Hwang, Y.T.; Lu, Q.; Hwang, C.B.C. Finger Domain Mutation Affects Enzyme Activity, DNA Replication Efficiency, and Fidelity of an Exonuclease-Deficient DNA Polymerase of Herpes Simplex Virus Type 1. *J. Virol.* **2009**, *83*, 7194–7201. [CrossRef] [PubMed]

27. Lunn, D.P.; Davis-Poynter, N.; Flaminio, M.J.B.F.; Horohov, D.W.; Osterrieder, K.; Pusterla, N.; Townsend, H.G.G. Equine herpesvirus-1 consensus statement. *J. Vet. Intern. Med.* **2009**, *23*, 450–461. [CrossRef]

28. Nugent, J.; Smith, K.C.; Mumford, J.A.; Newton, J.R.; Bowden, R.J.; Allen, G.P. Analysis of Equid Herpesvirus 1 Strain Variation Reveals a Point Mutation of the DNA Polymerase Strongly Associated with Neuropathogenic versus Nonneuropathogenic Disease Outbreaks. *J. Virol.* **2006**, *80*, 4047–4060. [CrossRef]

29. Stasiak, K.; Rola, J.; Ploszay, G.; Socha, W.; Zmudzinski, J.F. Detection of the neuropathogenic variant of equine herpesvirus 1 associated with abortions in mares in Poland. *BMC Vet. Res.* **2015**, *11*, 102. [CrossRef]

30. 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]

31. Pusterla, N.; Wilson, W.D.; Mapes, S.; Finno, C.; Isbell, D.; Arthur, R.M.; Ferraro, G.L. Characterization of viral loads, strain and state of equine herpesvirus-1 using real-time PCR in horses following natural exposure at a racetrack in California. *Vet. J.* **2009**, *179*, 230–239. [CrossRef]

32. Pusterla, N.; Mapes, S.; Wilson, W.D. Prevalence of equine herpesvirus type 1 in trigeminal ganglia and submandibular lymph nodes of equids examined postmortem. *Vet. Rec.* **2010**, *167*, 376–379. [CrossRef] [PubMed]

33. Savill, P.; Perrins, C.; Kirby, K.; Fisher, N. *Wytham Woods: Oxford's Ecological Laboratory*; Oxford University Press: Oxford, UK, 2010; Volume 48.

34. 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]

35. MacDonald, D.W.; 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]

36. Sugianto, N.A.; Buesching, C.D.; MacDonald, D.W.; Newman, C. The Importance of Refining Anaesthetic Regimes to Mitigate Adverse Effects in Very Young and Very Old Wild Animals: The European Badger (*Meles meles*). *J. Zool. Res.* **2019**, *3*, 10–17.

37. 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] [PubMed]

38. Okonechnikov, K.; Golosova, O.; Fursov, M.; Varlamov, A.; Vaskin, Y.; Efremov, I.; German Grehov, O.G.; Kandrov, D.; Rasputin, K.; Syabro, M.; et al. Unipro UGENE: A unified bioinformatics toolkit. *Bioinformatics* **2012**, *28*, 1166–1167. [CrossRef]

39. Nei, M.; Gojobori, T. Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. *Mol. Biol. Evol.* **1986**, *3*, 418–426. [CrossRef]

40. Miyata, T.; Katoh, K.; Misawa, K.; Kuma, K. MAFFT: A novel method for rapid multiple sequence alignment based on fast Fourier transform. *Nucleic Acids Res.* **2002**, *30*, 3059–3066.

41. 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]

42. Annavi, G.; Newman, C.; Dugdale, H.L.; Buesching, C.D.; Sin, Y.W.; Burke, T.; Macdonald, D.W. Neighbouring-group composition and within-group relatedness drive extra-group paternity rate in the European badger (*Meles meles*). *J. Evol. Biol.* **2014**, *27*, 2191–2203. [CrossRef]

43. Hancox, M. Field age determination in the european badger. *Rev. D'écologie* **1988**, *43*, 399–404.

44. R-Development-Core-Team R: A Language and Environment for Statistical Computing. Available online: https://www.r-project.org/ (accessed on 28 August 2022).

45. RStudio Team. RStudio: Integrated Development for R. RStudio. Available online: https://rstudio.com (accessed on 28 August 2022).

46. Lu, S.; Wang, J.; Chitsaz, F.; Derbyshire, M.K.; Geer, R.C.; Gonzales, N.R.; Gwadz, M.; Hurwitz, D.I.; Marchler, G.H.; Song, J.S.; et al. CDD/SPARCLE: The conserved domain database in 2020. *Nucleic Acids Res.* **2020**, *48*, D265–D268. [CrossRef] [PubMed]

47. Knopf, C.W. Evolution of viral DNA-dependent DNA polymerases. *Virus Genes* **1998**, *16*, 47–58. [CrossRef] [PubMed]

48. Ehlers, B.; Dural, G.; Yasmum, N.; Lembo, T.; de Thoisy, B.; Ryser-Degiorgis, M.-P.; Ulrich, R.G.; McGeoch, D.J. Novel Mammalian Herpesviruses and Lineages within the Gammaherpesvirinae: Cospeciation and Interspecies Transfer. *J. Virol.* **2008**, *82*, 3509–3516. [CrossRef]

49. Bauer, C.C.; Aberle, S.W.; Popow-Kraupp, T.; Kapitan, M.; Hofmann, H.; Puchhammer-Stöckl, E. Serum Epstein-Barr virus DNA load in primary Epstein-Barr virus infection. *J. Med. Virol.* **2005**, *75*, 54–58. [CrossRef]

50. Tsai, M.; Newman, C.; Macdonald, D.W.; Buesching, C.D. Stress-Related Herpesvirus Reactivation in Badgers Can Result in Clostridium Proliferation. *Ecohealth* **2021**, *18*, 440–450. [CrossRef] [PubMed]
51. Hendrikse, L.D.; Kambli, A.; Kayko, C.; Canuti, M.; Rodrigues, B.; Stevens, B.; Vashon, J.; Lang, A.S.; Needle, D.B.; Troyer, R.M. Identification of a novel gammaherpesvirus in Canada lynx (*Lynx canadensis*). *Viruses* **2019**, *11*, 363. [CrossRef] [PubMed]

52. Thorsteinsdóttir, L.; Torfason, E.G.; Torsteinsdóttir, S.; Svansson, V. Genetic diversity of equine gammaherpesviruses (γ-EHV) and isolation of a syncytium forming EHV-2 strain from a horse in Iceland. *Res. Vet. Sci.* **2013**, *94*, 170–177. [CrossRef] [PubMed]

53. Marmi, J.; López-Giráldez, F.; Macdonald, D.W.; Calafell, F.; Zholnerovskaya, E.; Domingo-Roura, X. Mitochondrial DNA reveals a strong phylogeographic structure in the badger across Eurasia. *Mol. Ecol.* **2006**, *15*, 1007–1020. [CrossRef] [PubMed]

54. Allen, A.; Guerrero, J.; Byrne, A.; Lavery, J.; Presho, E.; Courcier, E.; O’Keeffe, J.; Fogarty, U.; Delahay, R.; Wilson, G.; et al. Genetic evidence further elucidates the history and extent of badger introductions from Great Britain into Ireland. *R. Soc. Open Sci.* **2020**, *7*, 200288. [CrossRef]

55. Domingo-Roura, X.; Macdonald, D.W.; Roy, M.S.; Marmi, J.; Terradas, J.; Woodroffe, R.; Burke, T.; Wayne, R.K. Confirmation of low genetic diversity and multiple breeding females in a social group of Eurasian badgers from microsatellite and field data. *Mol. Ecol.* **2003**, *12*, 533–539. [CrossRef] [PubMed]

56. Alizon, S.; Hurford, A.; Mideo, N.; Van Baalen, M. Virulence evolution and the trade-off hypothesis: History, current state of affairs and the future. *J. Evol. Biol.* **2009**, *22*, 245–259. [CrossRef] [PubMed]

57. Abdul-Rahman, F.; Tranchina, D.; Gresham, D. Fluctuating Environments Maintain Genetic Diversity through Neutral Fitness Effects and Balancing Selection. *Mol. Biol. Evol.* **2021**, *38*, 4362–4375. [CrossRef] [PubMed]

58. Gibson, A.K.; Xu, J.Y.; Lively, C.M. Within-population covariation between sexual reproduction and susceptibility to local parasites. *Evolution* **2016**, *70*, 2049–2060. [CrossRef] [PubMed]

59. Tollenaere, C.; Bryja, J.; Galan, M.; Cadet, P.; Deter, J.; Chaval, Y.; Berthier, K.; Ribas Salvador, A.; Voutilainen, L.; Laakkonen, J.; et al. Multiple parasites mediate balancing selection at two MHC class II genes in the fossorial water vole: Insights from multivariate analyses and population genetics. *J. Evol. Biol.* **2008**, *21*, 1307–1320. [CrossRef] [PubMed]

60. Acevedo, M.A.; Dillemuth, F.P.; Flick, A.J.; Faldyn, M.J.; Elderd, B.D. Virulence-driven trade-offs in disease transmission: A meta-analysis. *Evolution* **2019**, *73*, 636–647. [CrossRef]

61. De Roode, J.C.; Yates, A.J.; Altizer, S. Virulence-transmission trade-offs and population divergence in virulence in a naturally occurring butterfly parasite. *Proc. Natl. Acad. Sci. USA* **2008**, *105*, 7489–7494. [CrossRef]

62. Ben-Shachar, R.; Koelle, K. Transmission-clearance trade-offs indicate that dengue virulence evolution depends on epidemiological context. *Nat. Commun.* **2018**, *9*, 2355. [CrossRef]

63. Rappocciolo, G.; Birch, J.; Ellis, S.A. Down-regulation of MHC class I expression by equine herpesvirus-1. *J. Gen. Virol.* **2003**, *84*, 293–300. [CrossRef]

64. Ma, G.; Feineis, S.; Osterrieder, N.; Van de Walle, G.R. Identification and Characterization of Equine Herpesvirus Type 1 pUL56 and Its Role in Virus-Induced Downregulation of Major Histocompatibility Complex Class I. *J. Virol.* **2012**, *86*, 3554–3563. [CrossRef]

65. Itzhaki, R.F. Corroboration of a major role for herpes simplex virus type 1 in Alzheimer's disease. *Front. Aging Neurosci.* **2018**, *10*, 1–11. [CrossRef] [PubMed]

66. Charonis, S.; James, L.M.; Georgopoulos, A.P. In silico assessment of binding affinities of three dementia-protective Human Leukocyte Antigen (HLA) alleles to nine human herpes virus antigens. *Curr. Res. Transl. Med.* **2020**, *68*, 211–216. [CrossRef] [PubMed]

67. Hulaniuk, M.L.; Mojsiejczuk, L.; Jauk, F.; Remondegui, C.; Mammana, L.; Bouzas, M.B.; Zapiola, I.; Ferro, M.V.; Ajalla, C.; Blejer, J.; et al. Genetic diversity and phylogeographic analysis of human herpesvirus type 8 (HHV-8) in two distant regions of Argentina: Association with the genetic ancestry of the population. *Infect. Genet. Evol.* **2020**, *85*, 104523. [CrossRef] [PubMed]