# Circulation of Bluetongue Virus Serotypes 1, 4, 8, 10 and 16 and Epizootic Hemorrhagic Disease Virus in the Sultanate of Oman in 2020–2021

Emmanuel Bréard, Lydie Postic, Mathilde Gondard, Cindy Bernelin-Cottet, Aurélie Le Roux, Mathilde Turpaud, Pierrick Lucas, Yannick Blanchard, Damien Vitour, Labib Bakkali-Kassimi, et al.

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

## To cite this version:

Emmanuel Bréard, Lydie Postic, Mathilde Gondard, Cindy Bernelin-Cottet, Aurélie Le Roux, et al. "Circulation of Bluetongue Virus Serotypes 1, 4, 8, 10 and 16 and Epizootic Hemorrhagic Disease Virus in the Sultanate of Oman in 2020–2021." Viruses, 2023, 15 (6), pp.1259. 10.3390/v15061259. anses-04145201

HAL Id: anses-04145201

https://anses.hal.science/anses-04145201

Submitted on 7 Mar 2024

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_2.png" /></figure>
# Circulation of Bluetongue Virus Serotypes 1, 4, 8, 10 and 16 and Epizootic Hemorrhagic Disease Virus in the Sultanate of Oman in 2020–2021

Emmanuel Bréard, Lydie Postic, Mathilde Gondard, Cindy Bernelin-Cottet, Aurélie Le Roux, Mathilde Turpaud, Pierrick Lucas, Yannick Blanchard, Damien Vitour, Labib Bakkali-Kassimi, et al.

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

## ► To cite this version:

Emmanuel Bréard, Lydie Postic, Mathilde Gondard, Cindy Bernelin-Cottet, Aurélie Le Roux, et al. "Circulation of Bluetongue Virus Serotypes 1, 4, 8, 10 and 16 and Epizootic Hemorrhagic Disease Virus in the Sultanate of Oman in 2020–2021." Viruses, 2023, 15 (6), pp.1259. 10.3390/v15061259. anses-04145201

HAL Id: anses-04145201

https://anses.hal.science/anses-04145201v1

Submitted on 7 Mar 2024

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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Circulation of Bluetongue Virus Serotypes 1, 4, 8, 10 and 16 and Epizootic Hemorrhagic Disease Virus in the Sultanate of Oman in 2020–2021

Emmanuel Bréard <fnref n="1" />,* ¹D, Lydie Postic ¹D, Mathilde Gondard ¹D, Cindy Bernelin-Cottet <fnref n="1" />, Aurélie Le Roux <fnref n="2" />, Mathilde Turpaud <fnref n="1" />, Pierrick Lucas <fnref n="2" />, Yannick Blanchard <fnref n="2" />, Damien Vitour <fnref n="1" />, Labib Bakkali-Kassimi <fnref n="1" />, Stéphan Zientara ¹D, Wafaa Al Rawahi <fnref n="3" />,<fnref n="4" /> and Corinne Sailleau <fnref n="1" />

**Abstract:** The circulation of Bluetongue (BT) and Epizootic Hemorrhagic Disease (EHD) in the Middle East has already been reported following serological analyses carried out since the 1980s, mostly on wild ruminants. Thus, an EHD virus (EHDV) strain was isolated in Bahrain in 1983 (serotype 6), and more recently, BT virus (BTV) serotypes 1, 4, 8 and 16 have been isolated in Oman. To our knowledge, no genomic sequence of these different BTV strains have been published. These same BTV or EHDV serotypes have circulated and, for some of them, are still circulating in the Mediterranean basin and/or in Europe. In this study, we used samples from domestic ruminant herds collected in Oman in 2020 and 2021 for suspected foot-and-mouth disease (FMD) to investigate the presence of BTV and EHDV in these herds. Sera and whole blood from goats, sheep and cattle were tested for the presence of viral genomes (by PCR) and antibodies (by ELISA). We were able to confirm the presence of 5 BTV serotypes (1, 4, 8, 10 and 16) and the circulation of EHDV in this territory in 2020 and 2021. The isolation of a BTV-8 strain allowed us to sequence its entire genome and to compare it with another BTV-8 strain isolated in Mayotte and with homologous BTV sequences available on GenBank.

**Keywords:** bluetongue; epizootic hemorrhagic disease; the Sultanate of Oman; BTV serotype 8

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

check for
updates

**Citation:** Bréard, E.; Postic, L.; Gondard, M.; Bernelin-Cottet, C.; Le Roux, A.; Turpaud, M.; Lucas, P.; Blanchard, Y.; Vitour, D.; Bakkali-Kassimi, L.; et al. Circulation of Bluetongue Virus Serotypes 1, 4, 8, 10 and 16 and Epizootic Hemorrhagic Disease Virus in the Sultanate of Oman in 2020–2021. *Viruses* **2023**, *15*, 1259. https://doi.org/10.3390/v15061259

Academic Editor: Yannick Simonin

Received: 15 May 2023

Revised: 24 May 2023

Accepted: 25 May 2023

Published: 27 May 2023

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

**Copyright:** © 2023 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/).

1. Introduction

Bluetongue (BT) and Epizootic Hemorrhagic Disease (EHD) are vector-borne diseases
affecting domestic and wild ruminants caused by two orbiviruses which can induce signifi-
cant production losses in livestock industries [1,2]. BT virus (BTV) and EHD virus (EHDV)
are arboviruses mainly transmitted by hematophagous insects (Culicoides midges) [3]. It is
assumed that the epidemiological characteristics of EHDV reflect those of BTV. Before the
2000s, BTV and EHDV were restricted to between 40° and 50° North latitude and 35° South
latitude, these latter being confined to the tropical and subtropical regions [1,2]. Since the
beginning of the 21st century, the situation has changed, especially in Europe where BTV
has been detected up to 60° north latitude [4].

BTV has a genome of 10 double-stranded RNA segments (S) encoding 7 structural (viral protein 1 (VP1) to VP7) and 6 non-structural proteins (NS1-NS4, NS3a and NS5) [5,6].

<footnote n="1">¹ UMR 1161 VIROLOGIE, INRA, ENVA, ANSES, Laboratory for Animal Health, Paris Est University, 94701 Maisons-Alfort, France; lydie.postic@anses.fr (L.P.); mathilde.gondard@anses.fr (M.G.); cindy.bernelin-cottet@anses.fr (C.B.-C.); mathilde.turpaud@anses.fr (M.T.); damien.vitour@anses.fr (D.V.); labib.bakkali-kassimi@anses.fr (L.B.-K.); stephan.zientara@anses.fr (S.Z.); corinne.sailleau@anses.fr (C.S.)</footnote>

<footnote n="2">² Laboratory of Ploufragan, ANSES, Unit of Viral Genetics and Biosafety, 22440 Ploufragan, France; aurelie.leroux@anses.fr (A.L.R.); pierrick.lucas@anses.fr (P.L.); yannick.blanchard@anses.fr (Y.B.)</footnote>

<footnote n="3">³ Department of Biology, College of Science, Sultan Qaboos University, Muscat P.C. 123, Oman; s69556@student.squ.edu.om</footnote>

<footnote n="4">⁴ Central Laboratory of Animal Health, Ministry of Agriculture, Fisheries Wealth and Water Resources, Muscat P.C. 100, Oman</footnote>

<footnote>* Correspondence: emmanuel.breard@anses.fr</footnote>
Article

# Circulation of Bluetongue Virus Serotypes 1, 4, 8, 10 and 16 and Epizootic Hemorrhagic Disease Virus in the Sultanate of Oman in 2020–2021

Emmanuel Bréard <fnref n="1" />,* ¹D, Lydie Postic ¹D, Mathilde Gondard ¹D, Cindy Bernelin-Cottet <fnref n="1" />, Aurélie Le Roux <fnref n="2" />, Mathilde Turpaud <fnref n="1" />, Pierrick Lucas <fnref n="2" />, Yannick Blanchard <fnref n="2" />, Damien Vitour <fnref n="1" />, Labib Bakkali-Kassimi <fnref n="1" />, Stéphan Zientara ¹D, Wafaa Al Rawahi <fnref n="3" />,<fnref n="4" /> and Corinne Sailleau <fnref n="1" />

**Abstract:** The circulation of Bluetongue (BT) and Epizootic Hemorrhagic Disease (EHD) in the Middle East has already been reported following serological analyses carried out since the 1980s, mostly on wild ruminants. Thus, an EHD virus (EHDV) strain was isolated in Bahrain in 1983 (serotype 6), and more recently, BT virus (BTV) serotypes 1, 4, 8 and 16 have been isolated in Oman. To our knowledge, no genomic sequence of these different BTV strains have been published. These same BTV or EHDV serotypes have circulated and, for some of them, are still circulating in the Mediterranean basin and/or in Europe. In this study, we used samples from domestic ruminant herds collected in Oman in 2020 and 2021 for suspected foot-and-mouth disease (FMD) to investigate the presence of BTV and EHDV in these herds. Sera and whole blood from goats, sheep and cattle were tested for the presence of viral genomes (by PCR) and antibodies (by ELISA). We were able to confirm the presence of 5 BTV serotypes (1, 4, 8, 10 and 16) and the circulation of EHDV in this territory in 2020 and 2021. The isolation of a BTV-8 strain allowed us to sequence its entire genome and to compare it with another BTV-8 strain isolated in Mayotte and with homologous BTV sequences available on GenBank.

**Keywords:** bluetongue; epizootic hemorrhagic disease; the Sultanate of Oman; BTV serotype 8

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

check for
updates

**Citation:** Bréard, E.; Postic, L.; Gondard, M.; Bernelin-Cottet, C.; Le Roux, A.; Turpaud, M.; Lucas, P.; Blanchard, Y.; Vitour, D.; Bakkali-Kassimi, L.; et al. Circulation of Bluetongue Virus Serotypes 1, 4, 8, 10 and 16 and Epizootic Hemorrhagic Disease Virus in the Sultanate of Oman in 2020–2021. *Viruses* **2023**, *15*, 1259. https://doi.org/10.3390/v15061259

Academic Editor: Yannick Simonin

Received: 15 May 2023

Revised: 24 May 2023

Accepted: 25 May 2023

Published: 27 May 2023

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

**Copyright:** © 2023 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/).

## 1. Introduction

Bluetongue (BT) and Epizootic Hemorrhagic Disease (EHD) are vector-borne diseases affecting domestic and wild ruminants caused by two orbiviruses which can induce significant production losses in livestock industries [1,2]. BT virus (BTV) and EHD virus (EHDV) are arboviruses mainly transmitted by hematophagous insects (Culicoides midges) [3]. It is assumed that the epidemiological characteristics of EHDV reflect those of BTV. Before the 2000s, BTV and EHDV were restricted to between 40° and 50° North latitude and 35° South latitude, these latter being confined to the tropical and subtropical regions [1,2]. Since the beginning of the 21st century, the situation has changed, especially in Europe where BTV has been detected up to 60° north latitude [4].

BTV has a genome of 10 double-stranded RNA segments (S) encoding 7 structural (viral protein 1 (VP1) to VP7) and 6 non-structural proteins (NS1-NS4, NS3a and NS5) [5,6].

<footnote n="1">¹ UMR 1161 VIROLOGIE, INRA, ENVA, ANSES, Laboratory for Animal Health, Paris Est University, 94701 Maisons-Alfort, France; lydie.postic@anses.fr (L.P.); mathilde.gondard@anses.fr (M.G.); cindy.bernelin-cottet@anses.fr (C.B.-C.); mathilde.turpaud@anses.fr (M.T.); damien.vitour@anses.fr (D.V.); labib.bakkali-kassimi@anses.fr (L.B.-K.); stephan.zientara@anses.fr (S.Z.); corinne.sailleau@anses.fr (C.S.)</footnote>

<footnote n="2">² Laboratory of Ploufragan, ANSES, Unit of Viral Genetics and Biosafety, 22440 Ploufragan, France; aurelie.leroux@anses.fr (A.L.R.); pierrick.lucas@anses.fr (P.L.); yannick.blanchard@anses.fr (Y.B.)</footnote>

<footnote n="3">³ Department of Biology, College of Science, Sultan Qaboos University, Muscat P.C. 123, Oman; s69556@student.squ.edu.om</footnote>

<footnote n="4">⁴ Central Laboratory of Animal Health, Ministry of Agriculture, Fisheries Wealth and Water Resources, Muscat P.C. 100, Oman</footnote>

<footnote>* Correspondence: emmanuel.breard@anses.fr</footnote>
BTV serotypes are recognized on the basis of specific interactions between neutralizing
antibodies and two structural proteins (VP2 and VP5) forming an outer-capsid layer.

The occurrence of BTV in the Middle East region is already known [7–9]. In Saudi Arabia and the United Arab Emirates [9], it has been shown that these viruses were enzootic following serological analyses carried out mostly on wild ruminants. In the Sultanate of Oman, ELISA and serum neutralization test identified at least three BTV serotypes (BTV-3, 4 and 22) present in 1987–1988 [7,8]. In 2009, BTV-1, 4, 8 and 16 were isolated in Oman (www.reoviridae.org/dsRNA_virus_proteins/outbreaks.htm#AHS-2010/, accessed on 14 May 2023), and in 2010, BTV-1 and 16 were detected in four Oryx antelopes imported from Oman to Croatia [10]. To our knowledge, no genomic sequence of these different BTV strains have been published to date.

In 1983, an EHDV was isolated in Bahrain. This virus was first classified as serotype 3, and since 2009, it is recognized that this strain is serotype 6 [11]. In Israel, several serotypes of this virus have been detected (serotypes 1, 6 and 7). In the Maghreb countries, the same serotypes have been detected with, since 2021 (in Tunisia), serotype 8 [12]. This strain was detected in Spain, Sardinia and Sicily in 2022, ending the free status of Europe with respect to this disease [13].

In Europe and in the Mediterranean Basin, continuous BTV incursions are observed involving different serotypes: BTV-1, 2, 3, 4, 6, 8, 9, 11 and 16 [14]. Remarkably, the BTV-8 strain present in Europe since 2006 is one of the rare BTV strains capable of inducing clinical signs in cattle due to its pathogenic and transplacental transmission properties [14].

In this study, we investigated BTV and EHDV genomes from 61 EDTA blood samples from domestic goats and cattle from different farms in Oman. These samples were collected from animals with clinical signs of FMD virus in 2020 and 2021. BTV antibodies were tested in 145 sera from sheep (n = 23), goats (n = 26) and cattle (n = 96). Anti-EHDV were also tested in the 96 cattle sera. The results obtained show a circulation of these two viruses in the Sultanate of Oman, with notably the detection of BTV-1, 4, 8, 10 and 16. Partial sequences of segment 2 of BTV-1, 4 and 10 were obtained by Sanger sequencing. A BTV-8 strain was isolated, and the whole genome sequences, determined by NGS, were notably compared with the BTV-8 genomes from strains present in Europe, in the Mediterranean Basin or in Africa and another BTV-8 isolated in Mayotte in 2016.

**2. Materials and Methods**

**2.1. Samples**

Sera (n = 145) and whole blood (n = 61) from goats, sheep and cattle were collected from domestic ruminant herds showing clinical signs of FMD in 2020 and 2021. The whole blood samples taken on EDTA were frozen at -80 °C and sent to the ANSES laboratory with dry ice. The sera were stored and sent at -20 °C. Due to the potential presence of FMD virus in the samples, all analyses performed for this study were performed in an A3 containment laboratory.

**2.2. Serological Analyses**

BTV and EHDV antibodies were tested using two VP7 competition enzyme-linked immunosorbent assays (ELISA) (ID SCREEN Bluetongue or EHDV competition kits, Innovative Diagnostics, France), according to the manufacturer's instructions. On the basis of the ELISAs cut-off values, the tested samples with an inhibition percentage <35% for bluetongue and <30% for EHD were considered positive.

**2.3. RT-qPCR Analyses**

Total RNA from EDTA blood samples was extracted from 100 µL of blood using a Kingfisher 96 robot and the ID extraction kit (Innovative Diagnostics, France) according to the manufacturer's instructions. Finally, the RNAs were eluted in 80 µL of water, and 5 µL of the eluted RNA previously heated to 95 °C for 3 min was tested with EHDV and BTV group-specific real time RT-PCRs (rt-qPCR). BTV rt-qPCRs were performed using a
BTV serotypes are recognized on the basis of specific interactions between neutralizing
antibodies and two structural proteins (VP2 and VP5) forming an outer-capsid layer.

The occurrence of BTV in the Middle East region is already known [7–9]. In Saudi Arabia and the United Arab Emirates [9], it has been shown that these viruses were enzootic following serological analyses carried out mostly on wild ruminants. In the Sultanate of Oman, ELISA and serum neutralization test identified at least three BTV serotypes (BTV-3, 4 and 22) present in 1987–1988 [7,8]. In 2009, BTV-1, 4, 8 and 16 were isolated in Oman (www.reoviridae.org/dsRNA_virus_proteins/outbreaks.htm#AHS-2010/, accessed on 14 May 2023), and in 2010, BTV-1 and 16 were detected in four Oryx antelopes imported from Oman to Croatia [10]. To our knowledge, no genomic sequence of these different BTV strains have been published to date.

In 1983, an EHDV was isolated in Bahrain. This virus was first classified as serotype 3, and since 2009, it is recognized that this strain is serotype 6 [11]. In Israel, several serotypes of this virus have been detected (serotypes 1, 6 and 7). In the Maghreb countries, the same serotypes have been detected with, since 2021 (in Tunisia), serotype 8 [12]. This strain was detected in Spain, Sardinia and Sicily in 2022, ending the free status of Europe with respect to this disease [13].

In Europe and in the Mediterranean Basin, continuous BTV incursions are observed involving different serotypes: BTV-1, 2, 3, 4, 6, 8, 9, 11 and 16 [14]. Remarkably, the BTV-8 strain present in Europe since 2006 is one of the rare BTV strains capable of inducing clinical signs in cattle due to its pathogenic and transplacental transmission properties [14].

In this study, we investigated BTV and EHDV genomes from 61 EDTA blood samples from domestic goats and cattle from different farms in Oman. These samples were collected from animals with clinical signs of FMD virus in 2020 and 2021. BTV antibodies were tested in 145 sera from sheep (n = 23), goats (n = 26) and cattle (n = 96). Anti-EHDV were also tested in the 96 cattle sera. The results obtained show a circulation of these two viruses in the Sultanate of Oman, with notably the detection of BTV-1, 4, 8, 10 and 16. Partial sequences of segment 2 of BTV-1, 4 and 10 were obtained by Sanger sequencing. A BTV-8 strain was isolated, and the whole genome sequences, determined by NGS, were notably compared with the BTV-8 genomes from strains present in Europe, in the Mediterranean Basin or in Africa and another BTV-8 isolated in Mayotte in 2016.

**2. Materials and Methods**

2.1. *Samples*

Sera (n = 145) and whole blood (n = 61) from goats, sheep and cattle were collected from domestic ruminant herds showing clinical signs of FMD in 2020 and 2021. The whole blood samples taken on EDTA were frozen at -80 °C and sent to the ANSES laboratory with dry ice. The sera were stored and sent at -20 °C. Due to the potential presence of FMD virus in the samples, all analyses performed for this study were performed in an A3 containment laboratory.

## 2.2. *Serological Analyses*

BTV and EHDV antibodies were tested using two VP7 competition enzyme-linked immunosorbent assays (ELISA) (ID SCREEN Bluetongue or EHDV competition kits, Innovative Diagnostics, France), according to the manufacturer's instructions. On the basis of the ELISAs cut-off values, the tested samples with an inhibition percentage <35% for bluetongue and <30% for EHD were considered positive.

2.3. *RT-qPCR Analyses*

Total RNA from EDTA blood samples was extracted from 100 µL of blood using a Kingfisher 96 robot and the ID extraction kit (Innovative Diagnostics, France) according to the manufacturer's instructions. Finally, the RNAs were eluted in 80 µL of water, and 5 µL of the eluted RNA previously heated to 95 °C for 3 min was tested with EHDV and BTV group-specific real time RT-PCRs (rt-qPCR). BTV rt-qPCRs were performed using a
commercial kit (ADI-352, Bio-X Diagnostics, Rochefort, Belgium). This pan-rt-qPCR BTV
kit allows all 36 BTV serotypes to be detected by amplification of BTV S10. EHDV RNA
was detected using a pan-EHDV rt-qPCR targeting the S9 [15]. These two pan rt-qPCRs
allow specific detection of BTV or EHDV regardless of serotype.

## 2.4. Virus Isolation

The rt-qPCR EHDV and BTV positive blood samples were inoculated to KC cells (Culi-
coides sonorensis cell line) or to embryonated chicken eggs. Briefly, a confluent monolayer of
KC cells was inoculated with EDTA-blood samples diluted 1 to 10 in sterile PBS. In paral-
lel, groups of three embryonated chicken eggs were each inoculated intravenously with
0.1-0.2 mL of this diluted blood. For KC cells, the inoculum was removed 24 h after inocu-
lation, and the cells were incubated 7 days with cell culture medium. The dead embryos or
KC culture cells were tested by pan EHDV or BTV rt-qPCR and then used to inoculate BSR
as already described [16].

## 2.5. Serotype Determination

End-point RT-PCRs were performed in the presence of RNA extracted from EHDV or BTV pan rt-qPCR positive bloods or cell cultures in order to amplify part of S2 of EHDV and BTV. Gel-based subgroup-specific RT-PCRs were described previously [15-17]. Commercial type-specific rt-qPCRs kits as well as rt-qPCRs methods published by Maan [18,19], Lorusso [13] and Viarouge [15] were also used for serotyping. These different PCRs were performed in presence of 5 µL of total RNA extracted and denatured at 95 °C for 3 min. After migration of the amplification products on agarose gel, the amplified S2 portions were then sequenced using Sanger technology (Eurofins, France).

## 2.6. Full Genome Sequencing

Next-generation sequencing was performed on total RNA extracts from infected BSR
cells using the QiaAmp Viral RNA mini extraction kit (Qiagen, Hilden, Germany) and the
Qiacube extractor. Then, 140 µL of infected culture media was extracted without carrier
in a final elution volume of 60 µL. One µL of the eluted RNA was tested by BTV pan
rt-qPCRs. After a rRNA depletion with the Low Input Ribominus Kit (Ambion, Hongkong,
China), the RNA library was performed with the Ion Total RNA-Seq Kit v2 (Life Tech-
nologies, Carlsbad, CA, USA), using manufacturer recommendations, then sequenced
with the Ion Torrent Proton Life Technologies, Shanghai, China. The resulting reads were
cleaned with the Trimmomatic 0.32 software, then a Bowtie 2 alignment was performed
on BTV genome references. The reads were down-sampled to fit a global coverage esti-
mation of 80 x and submitted to the SPAdes 3.1.1 de novo assembler. The de novo contigs
were then submitted to BLAST (Basic Local Alignment Search Tool, Bethesda, MD, USA)
on a local nt database. For each segment, the best matches were selected for a Bowtie
2 alignment producing clean and robust 5' and 3' ends. Finally, the de novo assemblies and
the alignment on the references were compared, and the strict identities of the de novo and
aligned sequences were assessed [20].

## 2.7. Genome Annotation and Phylogenetic Analyses

Each of the 10 segments of the Omani BTV-8 strain described in this study were
annotated using Geneious Prime (version 2022.0.2) and deposited in GenBank (acces-
sion number (AN) from S1 to S10: OQ860824 to OQ860833). An online BLAST search
was used to compare the ten segments of the Omani BTV-8 strain to reference sequences
listed in GenBank sequence databases (NCBI) and also to the 10 segments of a Mayotte
Island BTV-8 (Indian Ocean) isolated in 2016 in the laboratory and a South African (SA)
BTV-8 isolated in 1937. Among the 1000 first hits, we then selected hit sequences accord-
ing to genomic diversity and a geographical diversity for alignment and phylogenetic
analysis using MEGA X (version 10.2.0) [21]. Alignments were first performed using MUS-
CLE [22], and then phylogenetic trees were reconstructed using the Maximum Likelihood
commercial kit (ADI-352, Bio-X Diagnostics, Rochefort, Belgium). This pan-rt-qPCR BTV
kit allows all 36 BTV serotypes to be detected by amplification of BTV S10. EHDV RNA
was detected using a pan-EHDV rt-qPCR targeting the S9 [15]. These two pan rt-qPCRs
allow specific detection of BTV or EHDV regardless of serotype.

## 2.4. Virus Isolation

The rt-qPCR EHDV and BTV positive blood samples were inoculated to KC cells (Culi-
coides sonorensis cell line) or to embryonated chicken eggs. Briefly, a confluent monolayer of
KC cells was inoculated with EDTA-blood samples diluted 1 to 10 in sterile PBS. In paral-
lel, groups of three embryonated chicken eggs were each inoculated intravenously with
0.1-0.2 mL of this diluted blood. For KC cells, the inoculum was removed 24 h after inocu-
lation, and the cells were incubated 7 days with cell culture medium. The dead embryos or
KC culture cells were tested by pan EHDV or BTV rt-qPCR and then used to inoculate BSR
as already described [16].

## 2.5. Serotype Determination

End-point RT-PCRs were performed in the presence of RNA extracted from EHDV or BTV pan rt-qPCR positive bloods or cell cultures in order to amplify part of S2 of EHDV and BTV. Gel-based subgroup-specific RT-PCRs were described previously [15-17]. Commercial type-specific rt-qPCRs kits as well as rt-qPCRs methods published by Maan [18,19], Lorusso [13] and Viarouge [15] were also used for serotyping. These different PCRs were performed in presence of 5 µL of total RNA extracted and denatured at 95 °C for 3 min. After migration of the amplification products on agarose gel, the amplified S2 portions were then sequenced using Sanger technology (Eurofins, France).

## 2.6. Full Genome Sequencing

Next-generation sequencing was performed on total RNA extracts from infected BSR
cells using the QiaAmp Viral RNA mini extraction kit (Qiagen, Hilden, Germany) and the
Qiacube extractor. Then, 140 µL of infected culture media was extracted without carrier
in a final elution volume of 60 µL. One µL of the eluted RNA was tested by BTV pan
rt-qPCRs. After a rRNA depletion with the Low Input Ribominus Kit (Ambion, Hongkong,
China), the RNA library was performed with the Ion Total RNA-Seq Kit v2 (Life Tech-
nologies, Carlsbad, CA, USA), using manufacturer recommendations, then sequenced
with the Ion Torrent Proton Life Technologies, Shanghai, China. The resulting reads were
cleaned with the Trimmomatic 0.32 software, then a Bowtie 2 alignment was performed
on BTV genome references. The reads were down-sampled to fit a global coverage esti-
mation of 80 x and submitted to the SPAdes 3.1.1 de novo assembler. The de novo contigs
were then submitted to BLAST (Basic Local Alignment Search Tool, Bethesda, MD, USA)
on a local nt database. For each segment, the best matches were selected for a Bowtie
2 alignment producing clean and robust 5' and 3' ends. Finally, the de novo assemblies and
the alignment on the references were compared, and the strict identities of the de novo and
aligned sequences were assessed [20].

## 2.7. Genome Annotation and Phylogenetic Analyses

Each of the 10 segments of the Omani BTV-8 strain described in this study were
annotated using Geneious Prime (version 2022.0.2) and deposited in GenBank (acces-
sion number (AN) from S1 to S10: OQ860824 to OQ860833). An online BLAST search
was used to compare the ten segments of the Omani BTV-8 strain to reference sequences
listed in GenBank sequence databases (NCBI) and also to the 10 segments of a Mayotte
Island BTV-8 (Indian Ocean) isolated in 2016 in the laboratory and a South African (SA)
BTV-8 isolated in 1937. Among the 1000 first hits, we then selected hit sequences accord-
ing to genomic diversity and a geographical diversity for alignment and phylogenetic
analysis using MEGA X (version 10.2.0) [21]. Alignments were first performed using MUS-
CLE [22], and then phylogenetic trees were reconstructed using the Maximum Likelihood
method and Tamura-Nei model, with a bootstrap of 1000 [23]. The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. All positions containing gaps and missing data were eliminated. Further information is provided in the figure legends.

## 3. Results

### 3.1. *Serological Results*

* The seroprevalence for BT in Oman is high (74%) in cattle and illustrates a significant circulation of BTV in the country (Table 1). It seems less important in small ruminants with 50 and 30.4% of goats and sheep, respectively, showing antibodies against BTV (Table 2). Seroprevalence is also high for EHDV: more than one out of two cattle has antibodies against this virus (Table 1). In total, 46 (90.2%) of the 51 cattle with antibodies to EHDV were also positive in BTV ELISA.

<table><caption>**Table 1.** EHD and BTV ELISA results in cattle.</caption><thead><tr><th></th><th colspan="3">BTV ELISA</th><th></th></tr><tr><th></th><th>Negative</th><th>Doubtful</th><th>Positive</th><th>Number Tested</th></tr></thead><tbody><tr><th>number</th><td>23</td><td>2</td><td>71</td><td rowspan="2">96</td></tr><tr><th>%</th><td>24.0</td><td>2.1</td><td>74.0</td></tr><tr><th></th><th colspan="3">EHDV ELISA</th><th></th></tr><tr><th>number</th><td>44</td><td>1</td><td>51</td><td rowspan="2">96</td></tr><tr><th>%</th><td>45.8</td><td>1</td><td>53.1</td></tr></tbody></table>

<table><caption>**Table 2.** BTV ELISA results in goats and sheep.</caption><thead><tr><th></th><th></th><th colspan="3">BTV ELISA</th><th></th></tr><tr><th></th><th></th><th>Negative</th><th>Doubtful</th><th>Positive</th><th>Number tested</th></tr></thead><tbody><tr><th rowspan="2">Goats</th><th>number</th><td>13</td><td>0</td><td>13</td><td rowspan="2">26</td></tr><tr><th>%</th><td>50</td><td>0</td><td>50</td></tr><tr><th rowspan="2">Sheep</th><th>number</th><td>15</td><td>1</td><td>7</td><td rowspan="2">23</td></tr><tr><th>%</th><td>65.2</td><td>4.3</td><td>30.4</td></tr></tbody></table>

### 3.2. *Pan rt-qPCR Results and Virus Isolation*

* In total, 15 of the 61 animals tested were positive for BTV rt-qPCR (Ct value range: 22 to 34), and 11 of the 56 cattle tested were positive for pan EHDV rt-qPCR (Ct value range: 30 to 35). Three cattle were positive for both viruses. Viral isolation assays on embryonated eggs and KC cells were tested several times on each rt-qPCR-positive blood sample. Only one isolate of BTV was obtained, and no EHDV isolate was recovered.

### 3.3. *Typing and Sequencing*

* The Figure 1 illustrates the location of the herds where BTV and EHDV genome detection occurred. The rt-qPCRs from Maan et al. [18,19], performed on the RNA extracts from the blood, have shown the presence of the BTV-1, 4, 8, 10 and 16. The endpoint RT-PCRs [16] performed on the RNAs from the bloods allowed the determination of partial sequences of S2 of BTV-1, 4 and 10. No S2 sequence could be obtained for BTV-16-positive blood by classical RT-PCRs. The complete genome sequences of the BTV-8 strain were obtained by NGS.
method and Tamura-Nei model, with a bootstrap of 1000 [23]. The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. All positions containing gaps and missing data were eliminated. Further information is provided in the figure legends.

### 3. Results

#### 3.1. *Serological Results*

* The seroprevalence for BT in Oman is high (74%) in cattle and illustrates a significant circulation of BTV in the country (Table 1). It seems less important in small ruminants with 50 and 30.4% of goats and sheep, respectively, showing antibodies against BTV (Table 2). Seroprevalence is also high for EHDV: more than one out of two cattle has antibodies against this virus (Table 1). In total, 46 (90.2%) of the 51 cattle with antibodies to EHDV were also positive in BTV ELISA.

<table><caption>**Table 1.** EHD and BTV ELISA results in cattle.</caption><thead><tr><th></th><th colspan="3">BTV ELISA</th><th></th></tr><tr><th></th><th>Negative</th><th>Doubtful</th><th>Positive</th><th>Number Tested</th></tr></thead><tbody><tr><th>number</th><td>23</td><td>2</td><td>71</td><td rowspan="2">96</td></tr><tr><th>%</th><td>24.0</td><td>2.1</td><td>74.0</td></tr><tr><th></th><th colspan="3">EHDV ELISA</th><th></th></tr><tr><th>number</th><td>44</td><td>1</td><td>51</td><td rowspan="2">96</td></tr><tr><th>%</th><td>45.8</td><td>1</td><td>53.1</td></tr></tbody></table>

<table><caption>**Table 2.** BTV ELISA results in goats and sheep.</caption><thead><tr><th></th><th></th><th colspan="3">BTV ELISA</th><th></th></tr><tr><th></th><th></th><th>Negative</th><th>Doubtful</th><th>Positive</th><th>Number tested</th></tr></thead><tbody><tr><th rowspan="2">Goats</th><th>number</th><td>13</td><td>0</td><td>13</td><td rowspan="2">26</td></tr><tr><th>%</th><td>50</td><td>0</td><td>50</td></tr><tr><th rowspan="2">Sheep</th><th>number</th><td>15</td><td>1</td><td>7</td><td rowspan="2">23</td></tr><tr><th>%</th><td>65.2</td><td>4.3</td><td>30.4</td></tr></tbody></table>

#### 3.2. *Pan rt-qPCR Results and Virus Isolation*

* In total, 15 of the 61 animals tested were positive for BTV rt-qPCR (Ct value range: 22 to 34), and 11 of the 56 cattle tested were positive for pan EHDV rt-qPCR (Ct value range: 30 to 35). Three cattle were positive for both viruses. Viral isolation assays on embryonated eggs and KC cells were tested several times on each rt-qPCR-positive blood sample. Only one isolate of BTV was obtained, and no EHDV isolate was recovered.

#### 3.3. *Typing and Sequencing*

* The Figure 1 illustrates the location of the herds where BTV and EHDV genome detection occurred. The rt-qPCRs from Maan et al. [18,19], performed on the RNA extracts from the blood, have shown the presence of the BTV-1, 4, 8, 10 and 16. The endpoint RT-PCRs [16] performed on the RNAs from the bloods allowed the determination of partial sequences of S2 of BTV-1, 4 and 10. No S2 sequence could be obtained for BTV-16-positive blood by classical RT-PCRs. The complete genome sequences of the BTV-8 strain were obtained by NGS.
<figure><img src="image_5.png" /><figcaption>**Figure 1.** Location of the BTV serotypes and EHDV positive samples.</figcaption></figure>

* The EHDV type specific rt-qPCRs as well as the EHDV end-point RT-PCRs did not yield positive result.

### 3.4. BTV Genome Annotation and Phylogenetic Analyses

S2 partial sequences were obtained for BTV-1, 10 and 4 from pan-BTV rt-qPCR-positive samples tested with end-point RT-PCRs. These sequences were compared with homologous segments in GenBank (Table 3).

<table><caption>**Table 3.** Annotation of the BTV-1, 4 and 10 OMAN S2 and homology comparison with other BTV S2 deposited in GenBank.</caption><thead><tr><th>AN</th><th>Length<br>bp</th><th>CDS (nt<br>Position/Length AA)</th><th>Product</th><th>(% nt/AA)—Better Homology in<br>GenBank (AN)</th></tr></thead><tbody><tr><td>OQ860822</td><td>342</td><td>20–335/114</td><td>Partial VP2 BTV-1</td><td>(98.0/98.2)—BTV-1 Spa 2008 (KP821020)</td></tr><tr><td>OQ860823</td><td>318</td><td>1182–1523/105</td><td>Partial VP2 BTV-10</td><td>(93.4/99.0)—BTV-10 SA 2017 (MT078290)</td></tr><tr><td>OQ860904</td><td>2874</td><td>20–2890/956</td><td>Partial VP2 BTV-4</td><td>(97.5/98.6)—BTV-4 Spa 2003 (KP821067)</td></tr></tbody></table>
<figure><img src="image_5.png" /><figcaption>**Figure 1.** Location of the BTV serotypes and EHDV positive samples.</figcaption></figure>

* The EHDV type specific rt-qPCRs as well as the EHDV end-point RT-PCRs did not yield positive result.

### 3.4. BTV Genome Annotation and Phylogenetic Analyses

S2 partial sequences were obtained for BTV-1, 10 and 4 from pan-BTV rt-qPCR-positive samples tested with end-point RT-PCRs. These sequences were compared with homologous segments in GenBank (Table 3).

<table><caption>**Table 3.** Annotation of the BTV-1, 4 and 10 OMAN S2 and homology comparison with other BTV S2 deposited in GenBank.</caption><thead><tr><th>AN</th><th>Length<br>bp</th><th>CDS (nt<br>Position/Length AA)</th><th>Product</th><th>(% nt/AA)—Better Homology in<br>GenBank (AN)</th></tr></thead><tbody><tr><td>OQ860822</td><td>342</td><td>20–335/114</td><td>Partial VP2 BTV-1</td><td>(98.0/98.2)—BTV-1 Spa 2008 (KP821020)</td></tr><tr><td>OQ860823</td><td>318</td><td>1182–1523/105</td><td>Partial VP2 BTV-10</td><td>(93.4/99.0)—BTV-10 SA 2017 (MT078290)</td></tr><tr><td>OQ860904</td><td>2874</td><td>20–2890/956</td><td>Partial VP2 BTV-4</td><td>(97.5/98.6)—BTV-4 Spa 2003 (KP821067)</td></tr></tbody></table>
* The complete length of the S1, 2, 3, 4, 6, 7, 8, 9 and 10 were obtained for the BTV-8 genome, with the usual 5' and 3' extremities sequences (5'-GTTAAA... ACTTAC-3'). Regarding S5, only three nucleotides were missing at the 3' extremity. Thus, all expected CDS were recovered (Table 4).

* Table 4 shows the sequence homologies obtained between the 10 segments of the BTV-8 Oman strain and 3 other BTV-8 strains: one isolated in France in 2015 (AN: MN495893 to MN495902), the second isolated in Mayotte (Indian Ocean) in 2016 (AN: OQ860834 to OQ860843) and a 3rd isolated in SA in 1937 (AN: MT078269 to MT078278). Table 4 also shows the best homology found between BTV-8 Oman sequences and homologous sequences deposited on GenBank. Sequence alignments (in AA) of the VP3, NS1, VP5, VP7 and NS3 of the four BTV-8 strains (Table 4) show more than 99% for all these VPs.

<table><caption>**Table 4.** Annotation of the BTV-8 OMAN RNA segments/VP and homology comparison with other BTV segments deposited in GenBank.</caption><thead><tr><th rowspan="2">Seg.</th><th rowspan="2">AN</th><th rowspan="2">Length<br>bp</th><th rowspan="2">CDS (nt<br>Position/Length AA)</th><th rowspan="2">Product</th><th colspan="3">% nt and (AA) Homology between BTV8 Oman and</th><th rowspan="2">(%)—Better Homology in<br>GenBank (AN)</th></tr><tr><th>BTV-8<br>France</th><th>BTV-8<br>Mayotte</th><th>BTV-8<br>SA</th></tr></thead><tbody><tr><td>1</td><td>OQ860824</td><td>3944</td><td>12–3920/1302</td><td>VP1</td><td>93.6<br>(98.8)</td><td>92.5 (98.6)</td><td>96<br>(99.1)</td><td>97.8—BTV-1 Isr 2019<br>(OM502362)</td></tr><tr><td>2</td><td>OQ860825</td><td>2939</td><td>18–2903/961</td><td>VP2</td><td>94.6<br>(97.2)</td><td>98.7 (99)</td><td>92.5<br>(97.1)</td><td>95.9—BTV-8 Nig 1982<br>(AJ585184)</td></tr><tr><td>3</td><td>OQ860826</td><td>2772</td><td>18–2723/901</td><td>VP3</td><td>94.5<br>(99.7)</td><td>94.3 (99.7)</td><td>95.5<br>(99.8)</td><td>97.1—BTV-24 SA 2017<br>(MG255591)</td></tr><tr><td>4</td><td>OQ860827</td><td>1982</td><td>9–1943/644</td><td>VP4</td><td>94.2<br>(97.8)</td><td>94 (98.4)</td><td>95.2<br>(98.8)</td><td>97.5—BTV-3 Isr 2016<br>(MG344993)</td></tr><tr><td>5</td><td>OQ860828</td><td>1773</td><td>35–1693/552</td><td>NS1</td><td>93.7<br>(99.3)</td><td>96.9 (99.3)</td><td>93<br>(99.1)</td><td>97.7—BTV-4 Fr 2020<br>(OK018214)</td></tr><tr><td>6</td><td>OQ860829</td><td>1637</td><td>28–1608/526</td><td>VP5</td><td>96.8<br>(99.8)</td><td>96.6 (99.8)</td><td>96.2<br>(99.6)</td><td>96.9—BTV-8 Nig 1982<br>(AJ586706)</td></tr><tr><td>7</td><td>OQ860830</td><td>1156</td><td>18–1067/349</td><td>VP7</td><td>97.3<br>(100)</td><td>94.7 (99.4)</td><td>95.1<br>(100)</td><td>97.7—BTV-1 Sud 1987<br>(KP821626)</td></tr><tr><td>8</td><td>OQ860831</td><td>1125</td><td>20–1084/354</td><td>NS2</td><td>92.9<br>(96.9)</td><td>92.4 (97.7)</td><td>94.7<br>(97.5)</td><td>97.5—BTV-1 Isr 2019<br>(OM502359)</td></tr><tr><td>9</td><td>OQ860832</td><td>1049</td><td>16–1005/329</td><td>VP6</td><td>95.9<br>(94.8)</td><td>95 (93.0)</td><td>95.3<br>(95.5)</td><td>96.8—BTV-24 SA 2020<br>(MT090653)</td></tr><tr><td rowspan="3">10</td><td rowspan="3">OQ860833</td><td rowspan="3">822</td><td>182–415/77</td><td>NS4</td><td>98.7<br>(98.7)</td><td>98.7 (100)</td><td>98.7<br>(100)</td><td rowspan="3">98.8—BTV-15 Isr 2006<br>(JX272377)</td></tr><tr><td>20–709/229</td><td>NS3</td><td>95.7<br>(100)</td><td>98.6 (99.6)</td><td>92.3<br>(99.1)</td></tr><tr><td>108–287/59</td><td>NS5</td><td>93.9<br>(100)</td><td>93.3 (99.6)</td><td>93.3<br>(99.1)</td></tr></tbody></table>

* The phylogenetic tree of S2 (Figure 2) was performed with selected sequences of BTV-8, 18 and 23 that constitute the nucleotype D of the BTV S2 [24]. The tree shows two clusters of S2 of BTV-8: the Oman and Mayotte S2, which are very similar, form one of these clusters also consisting of S2s of European, Mediterranean and Nigerian BTV-8 strains. The second cluster includes the S2 of the SA and Kenia strains.
* The complete length of the S1, 2, 3, 4, 6, 7, 8, 9 and 10 were obtained for the BTV-8 genome, with the usual 5' and 3' extremities sequences (5'-GTTAAA... ACTTAC-3'). Regarding S5, only three nucleotides were missing at the 3' extremity. Thus, all expected CDS were recovered (Table 4).

* Table 4 shows the sequence homologies obtained between the 10 segments of the BTV-8 Oman strain and 3 other BTV-8 strains: one isolated in France in 2015 (AN: MN495893 to MN495902), the second isolated in Mayotte (Indian Ocean) in 2016 (AN: OQ860834 to OQ860843) and a 3rd isolated in SA in 1937 (AN: MT078269 to MT078278). Table 4 also shows the best homology found between BTV-8 Oman sequences and homologous sequences deposited on GenBank. Sequence alignments (in AA) of the VP3, NS1, VP5, VP7 and NS3 of the four BTV-8 strains (Table 4) show more than 99% for all these VPs.

<table><caption>**Table 4.** Annotation of the BTV-8 OMAN RNA segments/VP and homology comparison with other BTV segments deposited in GenBank.</caption><thead><tr><th rowspan="2">Seg.</th><th rowspan="2">AN</th><th rowspan="2">Length<br>bp</th><th rowspan="2">CDS (nt<br>Position/Length AA)</th><th rowspan="2">Product</th><th colspan="3">% nt and (AA) Homology between BTV8 Oman and</th><th rowspan="2">(%)—Better Homology in<br>GenBank (AN)</th></tr><tr><th>BTV-8<br>France</th><th>BTV-8<br>Mayotte</th><th>BTV-8<br>SA</th></tr></thead><tbody><tr><td>1</td><td>OQ860824</td><td>3944</td><td>12–3920/1302</td><td>VP1</td><td>93.6<br>(98.8)</td><td>92.5 (98.6)</td><td>96<br>(99.1)</td><td>97.8—BTV-1 Isr 2019<br>(OM502362)</td></tr><tr><td>2</td><td>OQ860825</td><td>2939</td><td>18–2903/961</td><td>VP2</td><td>94.6<br>(97.2)</td><td>98.7 (99)</td><td>92.5<br>(97.1)</td><td>95.9—BTV-8 Nig 1982<br>(AJ585184)</td></tr><tr><td>3</td><td>OQ860826</td><td>2772</td><td>18–2723/901</td><td>VP3</td><td>94.5<br>(99.7)</td><td>94.3 (99.7)</td><td>95.5<br>(99.8)</td><td>97.1—BTV-24 SA 2017<br>(MG255591)</td></tr><tr><td>4</td><td>OQ860827</td><td>1982</td><td>9–1943/644</td><td>VP4</td><td>94.2<br>(97.8)</td><td>94 (98.4)</td><td>95.2<br>(98.8)</td><td>97.5—BTV-3 Isr 2016<br>(MG344993)</td></tr><tr><td>5</td><td>OQ860828</td><td>1773</td><td>35–1693/552</td><td>NS1</td><td>93.7<br>(99.3)</td><td>96.9 (99.3)</td><td>93<br>(99.1)</td><td>97.7—BTV-4 Fr 2020<br>(OK018214)</td></tr><tr><td>6</td><td>OQ860829</td><td>1637</td><td>28–1608/526</td><td>VP5</td><td>96.8<br>(99.8)</td><td>96.6 (99.8)</td><td>96.2<br>(99.6)</td><td>96.9—BTV-8 Nig 1982<br>(AJ586706)</td></tr><tr><td>7</td><td>OQ860830</td><td>1156</td><td>18–1067/349</td><td>VP7</td><td>97.3<br>(100)</td><td>94.7 (99.4)</td><td>95.1<br>(100)</td><td>97.7—BTV-1 Sud 1987<br>(KP821626)</td></tr><tr><td>8</td><td>OQ860831</td><td>1125</td><td>20–1084/354</td><td>NS2</td><td>92.9<br>(96.9)</td><td>92.4 (97.7)</td><td>94.7<br>(97.5)</td><td>97.5—BTV-1 Isr 2019<br>(OM502359)</td></tr><tr><td>9</td><td>OQ860832</td><td>1049</td><td>16–1005/329</td><td>VP6</td><td>95.9<br>(94.8)</td><td>95 (93.0)</td><td>95.3<br>(95.5)</td><td>96.8—BTV-24 SA 2020<br>(MT090653)</td></tr><tr><td></td><td></td><td></td><td>182–415/77</td><td>NS4</td><td>98.7<br>(98.7)</td><td>98.7 (100)</td><td>98.7<br>(100)</td><td></td></tr><tr><td>10</td><td>OQ860833</td><td>822</td><td>20–709/229</td><td>NS3</td><td>95.7<br>(100)</td><td>98.6 (99.6)</td><td>92.3<br>(99.1)</td><td>98.8—BTV-15 Isr 2006<br>(JX272377)</td></tr><tr><td></td><td></td><td></td><td>108–287/59</td><td>NS5</td><td>93.9<br>(100)</td><td>93.3 (99.6)</td><td>93.3<br>(99.1)</td><td></td></tr></tbody></table>

* The phylogenetic tree of S2 (Figure 2) was performed with selected sequences of BTV-8, 18 and 23 that constitute the nucleotype D of the BTV S2 [24]. The tree shows two clusters of S2 of BTV-8: the Oman and Mayotte S2, which are very similar, form one of these clusters also consisting of S2s of European, Mediterranean and Nigerian BTV-8 strains. The second cluster includes the S2 of the SA and Kenia strains.
<figure><img src="image_6.png" /><figcaption>**Figure 2. Phylogenetic analysis of BTV S2 sequences.** Phylogenetic analysis of S2 sequences of BTV-8, 18 and 23 strains using the maximum likelihood method and Tamura-Nei model with 1000 bootstrap replicates in MEGA X. This analysis involved 27 nucleotide sequences, and there were a total of 2943 positions in the final dataset. Bootstrap values appeared at the corresponding nodes. In the phylogenetic tree, GenBank sequences, bluetongue serotype, country and year of sample collection are given. The sequences investigated in the present study are marked in blue, red, orange and purple.</figcaption></figure>

* The S6 phylogenetic tree (Figure 3) illustrates the homologies between S6 of BTV-8 and BTV-18, constituting the G nucleotype of BTV VP5 [25]. As observed with the S2 phylogenetic tree, the S6 of the Oman and Mayotte strains are closely related and form a cluster with BTV strains S6 isolated in Europe, in the Mediterranean basin and in Nigeria. In the cluster of the BTV-8 S6, two BTV-18 S6 from SA were found having a strong homology with BTV-8 S6 from the same area. These two SA BTV-18 S6 are clearly distinct from those isolated from American strains (Figure 3).
<figure><img src="image_6.png" /><figcaption>**Figure 2. Phylogenetic analysis of BTV S2 sequences.** Phylogenetic analysis of S2 sequences of BTV-8, 18 and 23 strains using the maximum likelihood method and Tamura-Nei model with 1000 bootstrap replicates in MEGA X. This analysis involved 27 nucleotide sequences, and there were a total of 2943 positions in the final dataset. Bootstrap values appeared at the corresponding nodes. In the phylogenetic tree, GenBank sequences, bluetongue serotype, country and year of sample collection are given. The sequences investigated in the present study are marked in blue, red, orange and purple.</figcaption></figure>

* The S6 phylogenetic tree (Figure 3) illustrates the homologies between S6 of BTV-8 and BTV-18, constituting the G nucleotype of BTV VP5 [25]. As observed with the S2 phylogenetic tree, the S6 of the Oman and Mayotte strains are closely related and form a cluster with BTV strains S6 isolated in Europe, in the Mediterranean basin and in Nigeria. In the cluster of the BTV-8 S6, two BTV-18 S6 from SA were found having a strong homology with BTV-8 S6 from the same area. These two SA BTV-18 S6 are clearly distinct from those isolated from American strains (Figure 3).
<figure><img src="image_7.png" /><figcaption>**Figure 3. Phylogenetic analysis of BTV S6 sequences.** Phylogenetic analysis of S6 sequences of BTV-8 and 18 strains using the Maximum Likelihood method and Tamura-Nei model (1000 replicates). This analysis involved 27 nucleotide sequences, and there were a total of 1591 positions in the final dataset. In the phylogenetic tree, GenBank sequences, species designations and strain names are given. The sequences investigated in the present study are marked in blue, red, orange and purple.</figcaption></figure>

* S1, 3 to 5 and 7 to 10 of the Oman strain have the strongest homologies with segments all originating from BTV of many serotypes strains circulating further west: in the Mediterranean basin (Israel, Lebanon, Cyprus, Maghreb...), in the European countries, in Central Africa (Sudan, Nigeria, Cameroon) or SA (See Supplementary Figures S1–S8 for phylogenetic analysis of each BTV segment). Note the presence of a S4 of a Chinese BTV-7 strain in the cluster otherwise presenting only BTV S4 from strains isolated in Africa, Europe or the Arabian Peninsula (Supplementary Figure S4), suggesting for this Chinese BTV-7 segment 4 a "western" origin.

**4. Discussion**

In this study, serological results on bovine samples show that 74% of the animals are seropositive against BTV and 53.1% against EHDV, suggesting a significant circulation of these two viruses in the south-eastern Arabian Peninsula. In comparison, BTV ELISA results obtained on small ruminants seem to show a lower seroprevalence (50% and 30.4% of positive results for goat and sheep, respectively). This result, to be taken with caution because few small ruminants were tested in this study, seems, however, to illustrate that cattle are better sentinels for BTV and EHDV circulation. BTV surveillance can be focused on cattle as their larger body size results in a greater range of attraction to Culicoides than
<figure><img src="image_7.png" /><figcaption>**Figure 3. Phylogenetic analysis of BTV S6 sequences.** Phylogenetic analysis of S6 sequences of BTV-8 and 18 strains using the Maximum Likelihood method and Tamura-Nei model (1000 replicates). This analysis involved 27 nucleotide sequences, and there were a total of 1591 positions in the final dataset. In the phylogenetic tree, GenBank sequences, species designations and strain names are given. The sequences investigated in the present study are marked in blue, red, orange and purple.</figcaption></figure>

* S1, 3 to 5 and 7 to 10 of the Oman strain have the strongest homologies with segments all originating from BTV of many serotypes strains circulating further west: in the Mediterranean basin (Israel, Lebanon, Cyprus, Maghreb...), in the European countries, in Central Africa (Sudan, Nigeria, Cameroon) or SA (See Supplementary Figures S1–S8 for phylogenetic analysis of each BTV segment). Note the presence of a S4 of a Chinese BTV-7 strain in the cluster otherwise presenting only BTV S4 from strains isolated in Africa, Europe or the Arabian Peninsula (Supplementary Figure S4), suggesting for this Chinese BTV-7 segment 4 a "western" origin.

**4. Discussion**

In this study, serological results on bovine samples show that 74% of the animals are seropositive against BTV and 53.1% against EHDV, suggesting a significant circulation of these two viruses in the south-eastern Arabian Peninsula. In comparison, BTV ELISA results obtained on small ruminants seem to show a lower seroprevalence (50% and 30.4% of positive results for goat and sheep, respectively). This result, to be taken with caution because few small ruminants were tested in this study, seems, however, to illustrate that cattle are better sentinels for BTV and EHDV circulation. BTV surveillance can be focused on cattle as their larger body size results in a greater range of attraction to Culicoides than
small ruminants and are therefore more likely to be involved in virus transmission [26,27]. The low seroprevalence rates of EHDV reported in small ruminants in previous studies, compared to those observed in cattle, explain why small ruminant samples were not tested in our study for either EHDV antibodies or genomes [28].

Despite our efforts, the EHDV serotype(s) present in the Sultanate of Oman samples remained undetermined. EHDV typing assays based on specific rt-qPCRs were performed on pan-EHDV-positive samples, and none of them gave a positive result. The circulation of an EHDV-8 could be ruled out as the rt-qPCR used in this study was developed to specifically detect the EHDV-8 strain recently described in the Maghreb and in southern Europe [12,13]. Concerning the other rt-qPCRs used for EHDV typing assays, these tools were developed 10 years ago from EHDV S2 sequences deposited in GenBank at this time. Thus, primers and probes used in these assays might not be able to detect divergent EHDV S2 sequences. Finally, end-point PCR, allowing amplification of S2 portions of EHDV belonging to the same nucleotype, also failed [11]. Here, again, the lack of sensitivity if sequence mismatches are present, as well as the low amount of EHDV viral genome detected in the blood (Ct values in pan-rt-qPCR between 30 and 34), could explain the absence of amplification.

In total, 15 of 61 animals tested were pan-BTV rt-qPCR-positive. Serotype 4 is the most detected serotype, with 7 cattle detected positive (out of 15) in BTV-4 specific rt-qPCR, located in different farms (Figure 1). Two animals from the same farm located in northern Oman were positive for BTV-8. Despite repeated viral isolation tests on these 15 BTV-positive samples, only one BTV-8 strain could be obtained and its genome fully sequenced. The presence of BTV-1, 4, 8, 10 and 16 was determined using the rt-qPCRs developed by Maan et al. [18], demonstrating the robustness of this method and its efficiency. Partial S2 sequences of BTV-1, 4 and 10 were obtained in classical RT-PCRs from the blood samples and sequenced by the Sanger method.

The genomic segments of the Mayotte and Oman BTV-8 strains have their strongest homologies with the multiple BTV serotypes that have circulated or are still circulating in Africa, in the Mediterranean basin and in Europe (Tables 3 and 4, Figures 2 and 3, and the Supplementary Figures S1–S8). None of the 10 segments of the BTV-8 Oman strain appear to have an eastern, Indian origin. This is also the case with the amplified partial S2 of BTV-1, 4 and 10. The highest homologies of homologous sequences from GenBank with the amplified BTV-1 and 4 partial S2 are with BTV-1 and 4 strains that circulated in the Mediterranean basin (Table 2). The BTV-10 S2 circulating in Oman has the highest nt homology with the SA vaccine strain (93.6% nt).

The S6 phylogenic tree shows two BTV-18 S6 having 98% of nt homology with the S6 of SA BTV-8. The VP5 of these BTV-18 have a homology of more than 99% with VP5 of the different BTV-8 strains, showing a common origin for these S6. The conformation of the VP2s of these serotypes 8 and 18, belonging to the same nucleotype, allows these VP2s to be associated with the same VP5.

Many BTV-8 strains circulating since 2006 in Europe and the Mediterranean basin, including Cyprus and Israel, have been isolated, sequenced and had their genomes deposited in GenBank. However, all of these BTV-8 strains have a common origin whose first isolate was detected in the Netherlands in 2006 [29]. They are very similar, and the strains circulating in Europe today are identical to the 2006 prototype. In Europe, few serotypes co-circulate at the same time and place. However, in 2008, it has been described that a wild BTV-8 strain isolated in France had reassorted with genomic segments of BTV-1 circulating in the same region [30], but this reassorted strain does not appear to have circulated in the field. The known reassortant BTV-8 strains that have circulated have been observed in Israel, where multiple BTV serotypes co-circulate, increasing the likelihood of reassortant strains [31]. The other genomic sequences of BTV-8 strains found in GenBank are African strains, mostly from SA, and they have between them an important genomic homology showing a common origin as illustrated by the phylogenic trees of S2 and 6. For these reasons, in our study, the genomic sequences obtained from the Oman and the Mayotte
small ruminants and are therefore more likely to be involved in virus transmission [26,27]. The low seroprevalence rates of EHDV reported in small ruminants in previous studies, compared to those observed in cattle, explain why small ruminant samples were not tested in our study for either EHDV antibodies or genomes [28].

Despite our efforts, the EHDV serotype(s) present in the Sultanate of Oman samples remained undetermined. EHDV typing assays based on specific rt-qPCRs were performed on pan-EHDV-positive samples, and none of them gave a positive result. The circulation of an EHDV-8 could be ruled out as the rt-qPCR used in this study was developed to specifically detect the EHDV-8 strain recently described in the Maghreb and in southern Europe [12,13]. Concerning the other rt-qPCRs used for EHDV typing assays, these tools were developed 10 years ago from EHDV S2 sequences deposited in GenBank at this time. Thus, primers and probes used in these assays might not be able to detect divergent EHDV S2 sequences. Finally, end-point PCR, allowing amplification of S2 portions of EHDV belonging to the same nucleotype, also failed [11]. Here, again, the lack of sensitivity if sequence mismatches are present, as well as the low amount of EHDV viral genome detected in the blood (Ct values in pan-rt-qPCR between 30 and 34), could explain the absence of amplification.

In total, 15 of 61 animals tested were pan-BTV rt-qPCR-positive. Serotype 4 is the most detected serotype, with 7 cattle detected positive (out of 15) in BTV-4 specific rt-qPCR, located in different farms (Figure 1). Two animals from the same farm located in northern Oman were positive for BTV-8. Despite repeated viral isolation tests on these 15 BTV-positive samples, only one BTV-8 strain could be obtained and its genome fully sequenced. The presence of BTV-1, 4, 8, 10 and 16 was determined using the rt-qPCRs developed by Maan et al. [18], demonstrating the robustness of this method and its efficiency. Partial S2 sequences of BTV-1, 4 and 10 were obtained in classical RT-PCRs from the blood samples and sequenced by the Sanger method.

The genomic segments of the Mayotte and Oman BTV-8 strains have their strongest homologies with the multiple BTV serotypes that have circulated or are still circulating in Africa, in the Mediterranean basin and in Europe (Tables 3 and 4, Figures 2 and 3, and the Supplementary Figures S1–S8). None of the 10 segments of the BTV-8 Oman strain appear to have an eastern, Indian origin. This is also the case with the amplified partial S2 of BTV-1, 4 and 10. The highest homologies of homologous sequences from GenBank with the amplified BTV-1 and 4 partial S2 are with BTV-1 and 4 strains that circulated in the Mediterranean basin (Table 2). The BTV-10 S2 circulating in Oman has the highest nt homology with the SA vaccine strain (93.6% nt).

The S6 phylogenic tree shows two BTV-18 S6 having 98% of nt homology with the S6 of SA BTV-8. The VP5 of these BTV-18 have a homology of more than 99% with VP5 of the different BTV-8 strains, showing a common origin for these S6. The conformation of the VP2s of these serotypes 8 and 18, belonging to the same nucleotype, allows these VP2s to be associated with the same VP5.

Many BTV-8 strains circulating since 2006 in Europe and the Mediterranean basin, including Cyprus and Israel, have been isolated, sequenced and had their genomes deposited in GenBank. However, all of these BTV-8 strains have a common origin whose first isolate was detected in the Netherlands in 2006 [29]. They are very similar, and the strains circulating in Europe today are identical to the 2006 prototype. In Europe, few serotypes co-circulate at the same time and place. However, in 2008, it has been described that a wild BTV-8 strain isolated in France had reassorted with genomic segments of BTV-1 circulating in the same region [30], but this reassorted strain does not appear to have circulated in the field. The known reassortant BTV-8 strains that have circulated have been observed in Israel, where multiple BTV serotypes co-circulate, increasing the likelihood of reassortant strains [31]. The other genomic sequences of BTV-8 strains found in GenBank are African strains, mostly from SA, and they have between them an important genomic homology showing a common origin as illustrated by the phylogenic trees of S2 and 6. For these reasons, in our study, the genomic sequences obtained from the Oman and the Mayotte
BTV-8 strains were compared with sequences of one BTV-8 prototype that circulated in Europe/Mediterranean Basin (i.e., BTV-8 France isolated in 2015) and one in SA isolated in 1937 (Tables 3 and 4).

The sequences of the BTV-8 strains from Oman and Mayotte add to the data on this serotype. This study shows that this serotype also circulates in the Indian Ocean, that their genomes are close and suggest that all these BTV-8 strains (from the Middle East region, Mayotte, Africa, Europe and the Mediterranean basin) have a common African origin.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v15061259/s1, Figure S1: Phylogenetic analysis of BTV S1 sequences. Figure S2: Phylogenetic analysis of BTV S3 sequences. Figure S3: Phylogenetic analysis of BTV S4 sequences. Figure S4: Phylogenetic analysis of BTV S5 sequences. Figure S5: Phylogenetic analysis of BTV S7 sequences. Figure S6: Phylogenetic analysis of BTV S8 sequences. Figure S7: Phylogenetic analysis of BTV S9 sequences. Figure S8: Phylogenetic analysis of BTV S10 sequences.

**Author Contributions:** Conceptualization, L.P., C.S. and E.B.; methodology, L.P., C.S., A.L.R., P.L., M.G. and E.B.; software, M.G., A.L.R., P.L. and Y.B.; validation, W.A.R., D.V., Y.B., E.B. and S.Z.; formal analysis, W.A.R., D.V., Y.B., E.B. and S.Z.; investigation, C.B.-C.; resources, C.B.-C., L.B.-K. and W.A.R.; data curation, L.P., M.G. and Y.B.; writing—original draft preparation, E.B.; writing—review and editing, all co-authors; visualization, E.B. and M.G.; supervision, L.B.-K. and D.V.; project administration, C.B.-C. and L.B.-K.; funding acquisition, S.Z. and E.B. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was partially funded by PREPMEDVET, grant number ANR-20-SEBM-0004.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data presented in this study are available on request from the corresponding author.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

## References

1. Gibbs, E.P.; Greiner, E.C. The epidemiology of bluetongue. *Comp. Immunol. Microbiol. Infect. Dis.* **1994**, *17*, 207–220. [CrossRef] [PubMed]

2. Maclachlan, N.J.; Zientara, S.; Savini, G.; Daniels, P.W. Epizootic haemorrhagic disease. *Rev. Sci. Tech.* **2015**, *34*, 341–351. [CrossRef] [PubMed]

3. Savini, G.; Afonso, A.; Mellor, P.; Aradaib, I.; Yadin, H.; Sanaa, M.; Wilson, W.; Monaco, F.; Domingo, M. Epizootic haemorragic disease. *Res. Vet. Sci.* **2011**, *91*, 1–17. [PubMed]

4. Hulten, C.; Frossling, J.; Chenais, E.; Sternberg Lewerin, S. Seroprevalence after vaccination of cattle and sheep against Bluetongue Virus (BTV) serotype 8 in Sweden. *Transbound. Emerg. Dis.* **2013**, *60*, 438–447. [CrossRef]

5. Mertens, P.P.; Brown, F.; Sangar, D.V. Assignment of the genome segments of bluetongue virus type 1 to the proteins which they encode. *Virology* **1984**, *135*, 207–217. [CrossRef]

6. Mohd Jaafar, F.; Monsion, B.; Mertens, P.P.C.; Attoui, H. Identification of Orbivirus Non-Structural Protein 5 (NS5), Its Role and Interaction with RNA/DNA in Infected Cells. *Int. J. Mol. Sci.* **2023**, *24*, 6845. [CrossRef]

7. al-Busaidy, S.M.; Mellor, P.S. Isolation and identification of arboviruses from the Sultanate of Oman. *Epidemiol. Infect.* **1991**, *106*, 403–413.

8. al-Busaidy, S.M.; Mellor, P.S. Epidemiology of bluetongue and related orbiviruses in the Sultanate of Oman. *Epidemiol. Infect.* **1991**, *106*, 167–178. [CrossRef]

9. Frolich, K.; Hamblin, C.; Jung, S.; Ostrowski, S.; Mwanzia, J.; Streich, W.J.; Anderson, J.; Armstrong, R.M.; Anajariyah, S. Serologic surveillance for selected viral agents in captive and free-ranging populations of Arabian oryx (Oryx leucoryx) from Saudi Arabia and the United Arab Emirates. *J. Wildl. Dis.* **2005**, *41*, 67–79. [CrossRef]

10. Bosnic, S.; Beck, R.; Listes, E.; Lojkic, I.; Savini, G.; Roic, B. Bluetongue virus in Oryx antelope (Oryx leucoryx) during the quarantine period in 2010 in Croatia. *Vet. Ital.* **2015**, *51*, 139–143.

11. Anthony, S.J.; Maan, S.; Maan, N.; Kgosana, L.; Bachanek-Bankowska, K.; Batten, C.; Darpel, K.E.; Sutton, G.; Attoui, H.; Mertens, P.P. Genetic and phylogenetic analysis of the outer-coat proteins VP2 and VP5 of epizootic haemorrhagic disease virus (EHDV):
BTV-8 strains were compared with sequences of one BTV-8 prototype that circulated in Europe/Mediterranean Basin (i.e., BTV-8 France isolated in 2015) and one in SA isolated in 1937 (Tables 3 and 4).

The sequences of the BTV-8 strains from Oman and Mayotte add to the data on this serotype. This study shows that this serotype also circulates in the Indian Ocean, that their genomes are close and suggest that all these BTV-8 strains (from the Middle East region, Mayotte, Africa, Europe and the Mediterranean basin) have a common African origin.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v15061259/s1, Figure S1: Phylogenetic analysis of BTV S1 sequences. Figure S2: Phylogenetic analysis of BTV S3 sequences. Figure S3: Phylogenetic analysis of BTV S4 sequences. Figure S4: Phylogenetic analysis of BTV S5 sequences. Figure S5: Phylogenetic analysis of BTV S7 sequences. Figure S6: Phylogenetic analysis of BTV S8 sequences. Figure S7: Phylogenetic analysis of BTV S9 sequences. Figure S8: Phylogenetic analysis of BTV S10 sequences.

**Author Contributions:** Conceptualization, L.P., C.S. and E.B.; methodology, L.P., C.S., A.L.R., P.L., M.G. and E.B.; software, M.G., A.L.R., P.L. and Y.B.; validation, W.A.R., D.V., Y.B., E.B. and S.Z.; formal analysis, W.A.R., D.V., Y.B., E.B. and S.Z.; investigation, C.B.-C.; resources, C.B.-C., L.B.-K. and W.A.R.; data curation, L.P., M.G. and Y.B.; writing—original draft preparation, E.B.; writing—review and editing, all co-authors; visualization, E.B. and M.G.; supervision, L.B.-K. and D.V.; project administration, C.B.-C. and L.B.-K.; funding acquisition, S.Z. and E.B. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was partially funded by PREPMEDVET, grant number ANR-20-SEBM-0004.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data presented in this study are available on request from the corresponding author.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

## References

1. Gibbs, E.P.; Greiner, E.C. The epidemiology of bluetongue. *Comp. Immunol. Microbiol. Infect. Dis.* **1994**, *17*, 207–220. [CrossRef] [PubMed]

2. Maclachlan, N.J.; Zientara, S.; Savini, G.; Daniels, P.W. Epizootic haemorrhagic disease. *Rev. Sci. Tech.* **2015**, *34*, 341–351. [CrossRef] [PubMed]

3. Savini, G.; Afonso, A.; Mellor, P.; Aradaib, I.; Yadin, H.; Sanaa, M.; Wilson, W.; Monaco, F.; Domingo, M. Epizootic haemorragic disease. *Res. Vet. Sci.* **2011**, *91*, 1–17. [PubMed]

4. Hulten, C.; Frossling, J.; Chenais, E.; Sternberg Lewerin, S. Seroprevalence after vaccination of cattle and sheep against Bluetongue Virus (BTV) serotype 8 in Sweden. *Transbound. Emerg. Dis.* **2013**, *60*, 438–447. [CrossRef]

5. Mertens, P.P.; Brown, F.; Sangar, D.V. Assignment of the genome segments of bluetongue virus type 1 to the proteins which they encode. *Virology* **1984**, *135*, 207–217. [CrossRef]

6. Mohd Jaafar, F.; Monsion, B.; Mertens, P.P.C.; Attoui, H. Identification of Orbivirus Non-Structural Protein 5 (NS5), Its Role and Interaction with RNA/DNA in Infected Cells. *Int. J. Mol. Sci.* **2023**, *24*, 6845. [CrossRef]

7. al-Busaidy, S.M.; Mellor, P.S. Isolation and identification of arboviruses from the Sultanate of Oman. *Epidemiol. Infect.* **1991**, *106*, 403–413.

8. al-Busaidy, S.M.; Mellor, P.S. Epidemiology of bluetongue and related orbiviruses in the Sultanate of Oman. *Epidemiol. Infect.* **1991**, *106*, 167–178. [CrossRef]

9. Frolich, K.; Hamblin, C.; Jung, S.; Ostrowski, S.; Mwanzia, J.; Streich, W.J.; Anderson, J.; Armstrong, R.M.; Anajariyah, S. Serologic surveillance for selected viral agents in captive and free-ranging populations of Arabian oryx (Oryx leucoryx) from Saudi Arabia and the United Arab Emirates. *J. Wildl. Dis.* **2005**, *41*, 67–79. [CrossRef]

10. Bosnic, S.; Beck, R.; Listes, E.; Lojkic, I.; Savini, G.; Roic, B. Bluetongue virus in Oryx antelope (Oryx leucoryx) during the quarantine period in 2010 in Croatia. *Vet. Ital.* **2015**, *51*, 139–143.

11. Anthony, S.J.; Maan, S.; Maan, N.; Kgosana, L.; Bachanek-Bankowska, K.; Batten, C.; Darpel, K.E.; Sutton, G.; Attoui, H.; Mertens, P.P. Genetic and phylogenetic analysis of the outer-coat proteins VP2 and VP5 of epizootic haemorrhagic disease virus (EHDV):
Comparison of genetic and serological data to characterise the EHDV serogroup. *Virus Res.* **2009**, *145*, 200–210. [CrossRef] [PubMed]

12. Sghaier, S.; Sailleau, C.; Marcacci, M.; Thabet, S.; Curini, V.; Ben Hassine, T.; Teodori, L.; Portanti, O.; Hammami, S.; Jurisic, L.; et al. Epizootic Haemorrhagic Disease Virus Serotype 8 in Tunisia, 2021. *Viruses* **2022**, *15*, 16. [CrossRef] [PubMed]

13. Lorusso, A.; Cappai, S.; Loi, F.; Pinna, L.; Ruiu, A.; Puggioni, G.; Guercio, A.; Purpari, G.; Vicari, D.; Sghaier, S.; et al. Epizootic Hemorrhagic Disease Virus Serotype 8, Italy, 2022. *Emerg. Infect. Dis.* **2023**, *29*, 1063–1065. [CrossRef]

14. Kundlacz, C.; Caignard, G.; Sailleau, C.; Viarouge, C.; Postic, L.; Vitour, D.; Zientara, S.; Breard, E. Bluetongue Virus in France: An Illustration of the European and Mediterranean Context since the 2000s. *Viruses* **2019**, *11*, 672. [CrossRef] [PubMed]

15. Viarouge, C.; Breard, E.; Zientara, S.; Vitour, D.; Sailleau, C. Duplex Real-Time RT-PCR Assays for the Detection and Typing of Epizootic Haemorrhagic Disease Virus. *PLoS ONE* **2015**, *10*, e0132540. [CrossRef]

16. Viarouge, C.; Lancelot, R.; Rives, G.; Breard, E.; Miller, M.; Baudrimont, X.; Doceul, V.; Vitour, D.; Zientara, S.; Sailleau, C. Identification of bluetongue virus and epizootic hemorrhagic disease virus serotypes in French Guiana in 2011 and 2012. *Vet. Microbiol.* **2014**, *174*, 78–85. [CrossRef] [PubMed]

17. Sailleau, C.; Zanella, G.; Breard, E.; Viarouge, C.; Desprat, A.; Vitour, D.; Adam, M.; Lasne, L.; Martrenchar, A.; Bakkali-Kassimi, L.; et al. Co-circulation of bluetongue and epizootic haemorrhagic disease viruses in cattle in Reunion Island. *Vet. Microbiol.* **2012**, *155*, 191–197. [CrossRef] [PubMed]

18. Maan, N.S.; Maan, S.; Belaganahalli, M.N.; Ostlund, E.N.; Johnson, D.J.; Nomikou, K.; Mertens, P.P. Identification and differentiation of the twenty six bluetongue virus serotypes by RT-PCR amplification of the serotype-specific genome segment 2. *PLoS ONE* **2012**, *7*, e32601. [CrossRef]

19. Maan, S.; Maan, N.S.; Belaganahalli, M.N.; Potgieter, A.C.; Kumar, V.; Batra, K.; Wright, I.M.; Kirkland, P.D.; Mertens, P.P. Development and Evaluation of Real Time RT-PCR Assays for Detection and Typing of Bluetongue Virus. *PLoS ONE* **2016**, *11*, e0163014. [CrossRef]

20. Briand, F.X.; Schmitz, A.; Ogor, K.; Le Prioux, A.; Guillou-Cloarec, C.; Guillemoto, C.; Allee, C.; Le Bras, M.O.; Hirchaud, E.; Quenault, H.; et al. Emerging highly pathogenic H5 avian influenza viruses in France during winter 2015/16: Phylogenetic analyses and markers for zoonotic potential. *Euro Surveill.* **2017**, *22*, 30473. [CrossRef]

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

22. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. *Nucleic Acids Res.* **2004**, *32*, 1792–1797. [CrossRef] [PubMed]

23. Tamura, K.; Nei, M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. *Mol. Biol. Evol.* **1993**, *10*, 512–526.

24. Maan, S.; Maan, N.S.; Samuel, A.R.; O'Hara, R.; Meyer, A.J.; Rao, S.; Mertens, P.P. Completion of the sequence analysis and comparisons of genome segment 2 (encoding outer capsid protein VP2) from representative isolates of the 24 bluetongue virus serotypes. *Vet. Ital.* **2004**, *40*, 484–488. [PubMed]

25. Maan, S.; Maan, N.S.; Nomikou, K.; Veronesi, E.; Bachanek-Bankowska, K.; Belaganahalli, M.N.; Attoui, H.; Mertens, P.P. Complete genome characterisation of a novel 26th bluetongue virus serotype from Kuwait. *PLoS ONE* **2011**, *6*, e26147. [CrossRef]

26. Elbers, A.R.; Meiswinkel, R. Culicoides (Diptera: Ceratopogonidae) host preferences and biting rates in the Netherlands: Comparing cattle, sheep and the black-light suction trap. *Vet. Parasitol.* **2014**, *205*, 330–337. [CrossRef]

27. Ayllon, T.; Nijhof, A.M.; Weiher, W.; Bauer, B.; Allene, X.; Clausen, P.H. Feeding behaviour of Culicoides spp. (Diptera: Ceratopogonidae) on cattle and sheep in northeast Germany. *Parasit. Vectors* **2014**, *7*, 34. [CrossRef]

28. Duan, Y.; Yang, Z.; Zhu, P.; Xiao, L.; Li, Z.; Li, Z.; Li, L.; Zhu, J. A serologic investigation of epizootic hemorrhagic disease virus in China between 2014 and 2019. *Virol. Sin.* **2022**, *37*, 513–520. [CrossRef]

29. Maan, S.; Maan, N.S.; Ross-smith, N.; Batten, C.A.; Shaw, A.E.; Anthony, S.J.; Samuel, A.R.; Darpel, K.E.; Veronesi, E.; Oura, C.A.; et al. Sequence analysis of bluetongue virus serotype 8 from the Netherlands 2006 and comparison to other European strains. *Virology* **2008**, *377*, 308–318. [CrossRef]

30. Nomikou, K.; Hughes, J.; Wash, R.; Kellam, P.; Breard, E.; Zientara, S.; Palmarini, M.; Biek, R.; Mertens, P. Widespread Reassortment Shapes the Evolution and Epidemiology of Bluetongue Virus following European Invasion. *PLoS Pathog.* **2015**, *11*, e1005056. [CrossRef]

31. Golender, N.; Eldar, A.; Ehrlich, M.; Kenigswald, G.; Shlamovitz, I.; Even-Tov, B.; Zamir, L.; Klement, E.; Bumbarov, V. Genomic Analysis Illustrated a Single Introduction and Evolution of Israeli Bluetongue Serotype 8 Virus Population 2008–2019. *Microorganisms* **2021**, *9*, 1955. [CrossRef] [PubMed]
Comparison of genetic and serological data to characterise the EHDV serogroup. *Virus Res.* **2009**, *145*, 200–210. [CrossRef] [PubMed]

12. Sghaier, S.; Sailleau, C.; Marcacci, M.; Thabet, S.; Curini, V.; Ben Hassine, T.; Teodori, L.; Portanti, O.; Hammami, S.; Jurisic, L.; et al. Epizootic Haemorrhagic Disease Virus Serotype 8 in Tunisia, 2021. *Viruses* **2022**, *15*, 16. [CrossRef] [PubMed]

13. Lorusso, A.; Cappai, S.; Loi, F.; Pinna, L.; Ruiu, A.; Puggioni, G.; Guercio, A.; Purpari, G.; Vicari, D.; Sghaier, S.; et al. Epizootic Hemorrhagic Disease Virus Serotype 8, Italy, 2022. *Emerg. Infect. Dis.* **2023**, *29*, 1063–1065. [CrossRef]

14. Kundlacz, C.; Caignard, G.; Sailleau, C.; Viarouge, C.; Postic, L.; Vitour, D.; Zientara, S.; Breard, E. Bluetongue Virus in France: An Illustration of the European and Mediterranean Context since the 2000s. *Viruses* **2019**, *11*, 672. [CrossRef] [PubMed]

15. Viarouge, C.; Breard, E.; Zientara, S.; Vitour, D.; Sailleau, C. Duplex Real-Time RT-PCR Assays for the Detection and Typing of Epizootic Haemorrhagic Disease Virus. *PLoS ONE* **2015**, *10*, e0132540. [CrossRef]

16. Viarouge, C.; Lancelot, R.; Rives, G.; Breard, E.; Miller, M.; Baudrimont, X.; Doceul, V.; Vitour, D.; Zientara, S.; Sailleau, C. Identification of bluetongue virus and epizootic hemorrhagic disease virus serotypes in French Guiana in 2011 and 2012. *Vet. Microbiol.* **2014**, *174*, 78–85. [CrossRef] [PubMed]

17. Sailleau, C.; Zanella, G.; Breard, E.; Viarouge, C.; Desprat, A.; Vitour, D.; Adam, M.; Lasne, L.; Martrenchar, A.; Bakkali-Kassimi, L.; et al. Co-circulation of bluetongue and epizootic haemorrhagic disease viruses in cattle in Reunion Island. *Vet. Microbiol.* **2012**, *155*, 191–197. [CrossRef] [PubMed]

18. Maan, N.S.; Maan, S.; Belaganahalli, M.N.; Ostlund, E.N.; Johnson, D.J.; Nomikou, K.; Mertens, P.P. Identification and differentiation of the twenty six bluetongue virus serotypes by RT-PCR amplification of the serotype-specific genome segment 2. *PLoS ONE* **2012**, *7*, e32601. [CrossRef]

19. Maan, S.; Maan, N.S.; Belaganahalli, M.N.; Potgieter, A.C.; Kumar, V.; Batra, K.; Wright, I.M.; Kirkland, P.D.; Mertens, P.P. Development and Evaluation of Real Time RT-PCR Assays for Detection and Typing of Bluetongue Virus. *PLoS ONE* **2016**, *11*, e0163014. [CrossRef]

20. Briand, F.X.; Schmitz, A.; Ogor, K.; Le Prioux, A.; Guillou-Cloarec, C.; Guillemoto, C.; Allee, C.; Le Bras, M.O.; Hirchaud, E.; Quenault, H.; et al. Emerging highly pathogenic H5 avian influenza viruses in France during winter 2015/16: Phylogenetic analyses and markers for zoonotic potential. *Euro Surveill.* **2017**, *22*, 30473. [CrossRef]

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

22. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. *Nucleic Acids Res.* **2004**, *32*, 1792–1797. [CrossRef] [PubMed]

23. Tamura, K.; Nei, M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. *Mol. Biol. Evol.* **1993**, *10*, 512–526.

24. Maan, S.; Maan, N.S.; Samuel, A.R.; O'Hara, R.; Meyer, A.J.; Rao, S.; Mertens, P.P. Completion of the sequence analysis and comparisons of genome segment 2 (encoding outer capsid protein VP2) from representative isolates of the 24 bluetongue virus serotypes. *Vet. Ital.* **2004**, *40*, 484–488. [PubMed]

25. Maan, S.; Maan, N.S.; Nomikou, K.; Veronesi, E.; Bachanek-Bankowska, K.; Belaganahalli, M.N.; Attoui, H.; Mertens, P.P. Complete genome characterisation of a novel 26th bluetongue virus serotype from Kuwait. *PLoS ONE* **2011**, *6*, e26147. [CrossRef]

26. Elbers, A.R.; Meiswinkel, R. Culicoides (Diptera: Ceratopogonidae) host preferences and biting rates in the Netherlands: Comparing cattle, sheep and the black-light suction trap. *Vet. Parasitol.* **2014**, *205*, 330–337. [CrossRef]

27. Ayllon, T.; Nijhof, A.M.; Weiher, W.; Bauer, B.; Allene, X.; Clausen, P.H. Feeding behaviour of Culicoides spp. (Diptera: Ceratopogonidae) on cattle and sheep in northeast Germany. *Parasit. Vectors* **2014**, *7*, 34. [CrossRef]

28. Duan, Y.; Yang, Z.; Zhu, P.; Xiao, L.; Li, Z.; Li, Z.; Li, L.; Zhu, J. A serologic investigation of epizootic hemorrhagic disease virus in China between 2014 and 2019. *Virol. Sin.* **2022**, *37*, 513–520. [CrossRef]

29. Maan, S.; Maan, N.S.; Ross-smith, N.; Batten, C.A.; Shaw, A.E.; Anthony, S.J.; Samuel, A.R.; Darpel, K.E.; Veronesi, E.; Oura, C.A.; et al. Sequence analysis of bluetongue virus serotype 8 from the Netherlands 2006 and comparison to other European strains. *Virology* **2008**, *377*, 308–318. [CrossRef]

30. Nomikou, K.; Hughes, J.; Wash, R.; Kellam, P.; Breard, E.; Zientara, S.; Palmarini, M.; Biek, R.; Mertens, P. Widespread Reassortment Shapes the Evolution and Epidemiology of Bluetongue Virus following European Invasion. *PLoS Pathog.* **2015**, *11*, e1005056. [CrossRef]

31. Golender, N.; Eldar, A.; Ehrlich, M.; Kenigswald, G.; Shlamovitz, I.; Even-Tov, B.; Zamir, L.; Klement, E.; Bumbarov, V. Genomic Analysis Illustrated a Single Introduction and Evolution of Israeli Bluetongue Serotype 8 Virus Population 2008–2019. *Microorganisms* **2021**, *9*, 1955. [CrossRef] [PubMed]