key: cord-0983945-7vxif3y5
authors: Golender, Natalia; Bumbarov, Velizar; Eldar, Avi; Lorusso, Alessio; Kenigswald, Gabriel; Varsano, Joseph Seffi; David, Dan; Schainin, Shani; Dagoni, Ilan; Gur, Iosef; Kaplan, Alon; Gorohov, Anna; Koren, Ori; Oron, Eldad; Khinich, Yevgeny; Sclamovich, Ily; Meir, Abraham; Savini, Giovanni
title: Bluetongue Serotype 3 in Israel 2013–2018: Clinical Manifestations of the Disease and Molecular Characterization of Israeli Strains
date: 2020-03-06
journal: Front Vet Sci
DOI: 10.3389/fvets.2020.00112
sha: 90637e4f46c7466511a622cb1cce5ed3c8eb4f3d
doc_id: 983945
cord_uid: 7vxif3y5

In this paper, the results of the diagnostic activities on Bluetongue virus serotype 3 (BTV-3) conducted at Kimron Veterinary Institute (Beit Dagan, Israel) between 2013 and 2018 are reported. Bluetongue virus is the causative agent of bluetongue (BT), a disease of ruminants, mostly transmitted by competent Culicoides species. In Israel, BTV-3 circulation was first detected in 2013 from a sheep showing classical BT clinical signs. It was also evidenced in 2016, and, since then, it has been regularly detected in Israeli livestock. Between 2013 and 2017, BTV-3 outbreaks were limited in sheep flocks located in the southern area only. In 2018, BTV-3 was instead found in the Israeli coastal area being one of the dominant BTV serotypes isolated from symptomatic sheep, cattle and goats. In Israeli sheep, BTV-3 was able to cause BT classical clinical manifestations and fatalities, while in cattle and goats infection ranged from asymptomatic forms to death cases, depending on either general welfare of the herds or on the occurrence of viral and bacterial co-infections. Three different BTV-3 strains were identified in Israel between 2013 and 2018: ISR-2019/13 isolated in 2013, ISR-2153/16 and ISR-2262/2/16 isolated in 2016. Sequencing and phylogenetic analysis of these strains showed more than 99% identity by segment (Seg) 2, 5, 6, 7, and 8 sequences. In contrast, a wide range of diversity among these strains was exhibited in other viral gene segments, implying the occurrence of genome reassortment between these local circulating strains and those originating from Africa. The genome sequences of the BTV-3 isolated in 2017 and 2018 were most closely related to those of the ISR-2153/16 strain suggesting their common ancestor. Comparison of BTV-3 Israeli strains with those recently detected in the Mediterranean region uncovered high percentage identity (98.19–98.28%) only between Seg-2 of all Israeli strains and the BTV-3 Zarzis/TUN2016 strain. A 98.93% identity was also observed between Seg-4 sequences of ISR-2019/13 and the BTV-3 Zarzis/TUN2016 strain. This study demonstrated that BTV-3 has been circulating in the Mediterranean region at least since 2013, but, unlike the other Mediterranean strains, Israeli BTV-3 were able to cause clinical signs also in cattle.

Bluetongue (BT) is a non-contagious, arthropod-borne viral disease of domestic and wild ruminants, listed as a notifiable disease by the World Organization for Animal Health (OIE). The Bluetongue virus (BTV) is the prototype member of the Orbivirus genus within the Reoviridae family (1, 2).

BTV is a double-stranded (ds) RNA virus (3), its genome consists of ten segments (Seg-1 to Seg-10) of linear dsRNA coding 7 structural (VP1-VP7) and 5 non-structural (NS1, NS2, NS3/NS3a, NS4, and NS5) proteins (4) At present twenty-eighth distinct BTV serotypes have been officially recognized based on Seg-2 gene sequence (5, 6) . Other putative novel BTV serotypes have also been described (7) (8) (9) (10) (11) .

Transmission between mammalian hosts and spread of the infection rely mostly on competent Culicoides species (12) (13) (14) , so the presence of the disease is then strictly related to the distribution of competent vectors (2) . Even though they don't seem to be epidemiologically important, vertical and horizontal transmissions have also been described (15) (16) (17) (18) . For BTV-26, BTV-27 v02 and, probably, BTV-X ITL2015 and BTV-28 transmission by direct contact has been demonstrated or hypothesized (6, 10, 19, 20) .

Clinical signs of BT are more severe and most commonly observed in sheep or in white-tailed deer, often leading to animal fatality especially in naïve animals (1, 2). In cattle, BTV infection is usually asymptomatic, although symptoms were reported after infection with some strains (21) (22) (23) (24) .

As a RNA virus with a segmented genome, BTV can undergo reassortment which can occur when a cell is simultaneously infected with more than one BTV strain and involves the packaging, into a single virion, of full length of genomic segments of different ancestry. Reassortment in BTV is very flexible, and can involve any genome segment (25) (26) (27) . However, the genome sequence of BTV isolates generally reflects their geographic origins (28) (29) (30) .

Seg-2 and Seg-6 encoded the BTV outer capsid proteins VP2 and VP5. They represent the primary target for neutralizing antibodies generated during infection of the mammalian host (31) (32) (33) . Their highly variable sequences are associated with virus serotype (particularly in VP2/segment) (34) (35) (36) and, within each serotype, with the geographic origin of the virus strain (28, (34) (35) (36) (37) .

Among BTV serotypes, Seg-2/VP2 sequences of BTV-3, BTV-16, and BTV-13 are closely related. All of them are included within the nucleotype B reflecting a serological relationship (3). BTV-3 strains can be subdivided into at least two main clusters. 1st cluster includes strains originated from Africa, Mediterranean Basin and North America (western topotypes, w); 2nd-includes strains originated from Japan, India, and Australia (eastern topotypes, e). The nucleotide (nt) identity between Western and Eastern topotype Seg-2 can be as small as 71.5% (7, (38) (39) (40) (41) (42) (43) .

Between 2016 and 2018, two novel BTV-3 western strains have been identified in two different geographical areas of Tunisiaone in the north-eastern part of the country (Peninsula of Cap Bon, prototype BTV-3 TUN2016) and the other in the South-East near by the border with Libya (prototype strain BTV-3 TUN2016/Zarzis). The BTV-3 TUN2016 spread in 2017 to Italy infecting a single 3-year-old female crossbred sheep belonging to a flock located in the municipality of Trapani, Sicily, which are 150 km distant from Peninsula of Cap Bon (7, 39, 40) and in 2018 in the Southern area of Sardinia causing numerous outbreaks (38) . Clinical signs in infected sheep included depression, fever, nasal discharge, submandibular edema, and crusted discharge around the nostrils. Four animals died because of the severity of infection (38) . In 2016, another BTV-3 strain closely related to TUN2016/Zarzis strain was detected in Egypt (38) .

The present paper reports the results of the diagnostic activities on BTV conducted at Kimron Veterinary Institute between 2013 and 2018 and the evidence of BTV-3 circulation in Israel. Clinical signs of infected sheep, goats and cattle along with the genetic characterization and phylogenetic analysis of the BTV-3 strains are also described.

During years 2013, 2016, 2017, and 2018, 3,149 samples (714 in 2013, 669 in 2016, 744 in 2017, and 1,022 in 2018) from domestic and wild ruminants were collected and examined at the Kimron Institute, Beit Dagan, Israel (KVI). Samples included whole blood from symptomatic animals, spleen or/and lung samples from dead animals and spleen, lung, placenta and brain samples from aborted fetuses. Details on number, samples, species from which samples were collected are shown in Table 1 .

Farms where domestic clinically ill ruminants were confirmed by laboratory tests as BTV-3 infected were recorded and numbered in Table 3 . Information in Table 3 included location of the farm or animal grazing place, type of farming, breed, and number of animals in the farm/group. All clinical events observed in the farms were followed and reported ( Table 4) .

Internal organs were examined for aerobic and anaerobic bacteria and Salmonella spp. growth according to standard procedures (45) .

From all field samples and virus isolates (chicken embryo homogenate and tissue culture supernatant) viral RNA was extracted using Invisorb Spin Virus RNA Mini Kit (STRATEC Molecular, Berlin, Germany), and by MagMAX TM CORE Nucleic Acid Purification Kit (Thermo Fisher Scientific), according to recommendations of manufacturers. To detect bovine ephemeral fever virus (BEF), the RNA extraction and RT-qPCR were performed on white blood cells from whole EDTA blood samples of cattle origin according to Erster et al. (46) .

Presence of Malignant Catarrhal Fever virus (MCFV) DNA was determined according to the method described by Cunha (47) . Most samples from cattle were also tested for the presence of Epizootic Hemorrhagic disease virus (EHDV) RNA by EHDV Real-Time PCR Kit (Applied Biosystems TM , Thermo Fisher Scientific Inc., Lissieu, France) according to manufacturer's instructions. The method described by Boxus (48) was used to No. of isolated BTV-6 2 2

No. of isolated BTV- 15 9 9 No., number; w.b, whole blood samples; s/l, spleen or lung samples; s., spleen samples; a.f., aborted fetus; VI, virus isolation. Data from 2017 was also used in publication (44) .

detect Bovine Respiratory Syncytial virus (BRSV) RNA from lung cattle samples, while an in-house RT-qPCR (unpublished) and Kishimoto et al. (49) methods were used for detecting bovine Parainfluenza 3 (BPI3) and bovine Corona virus (BCoV) RNAs, respectively. A specific RT-qPCR was used for detection and identification Simbu serogroup RNA viruses in plasma samples (50) and, finally, the VetMAX TM BTV NS3 All Genotypes Kit (Applied Biosystems TM , Thermo Fisher Scientific Inc., Lissieu, France) RT-qPCR kit targeting Seg-10 of the BTV genome was used for detecting BTV RNA. Those samples which were positive to the BTV RT-qPCR with Cycle Threshold (Ct) values ≤33, were further tested for determining the serotype. Virotype R BTV pan/4 RT-PCR and Virotype R BTV pan/8 RT-PCR kits (QIAGEN, Leipzig, Germany) were applied directly to RNA extracted from diagnostic samples to detect and identify BTV-4 and BTV-8 serotypes, respectively. The same samples were also retroactively tested for the presence of BTV-3 by an in-house specific RT-qPCR according to the method described by Lorusso et al. (7) .

In an attempt to isolate BTV, all BTV RT-qPCR positive samples were inoculated in 9 to 11 day old embryonated chicken eggs (ECE) according to the method described by Komarov and Goldsmith (51) and adopted by Golender et al. (52) . In the homogenate of the ECE organs, the presence of BTV was confirmed by BTV RT-qPCR. BTV-positive ECE homogenates were subsequently tested for the presence of the RNA of all BTV serotypes that have circulated in Israel and neighboring countries during the recent decades (BTV-2, −3, −4, −5, −8, −12, −15, −16, and−24) by using a conventional specific in-house RT-PCR (primers are listed in Table S1 ) using a One-Step RT-PCR kit (Qiagen, Hilden, Germany).

Primers 3VP2-1F (5 ′ -GTT AAA AAC GCT GTC CCG AGA-3 ′ ) and 3VP2-637R (5 ′ -GAG CGC CCA CTC TAA ATT CCT C-3 ′ ) targeting a portion of 658 bp of the Seg-2 of BTV-3 by RT-PCR were used for BTV-3 identification. All positive in in-house conventional serotype specific RT-PCR tested BT virus isolates from all serotypes were consequently sequenced.

All viral-segment amplicons of Israeli BTV-3 were sequenced by standard Sanger methods with an ABI 3730xl DNA Analyzer (Hylabs, Rehovot, Israel). The cDNA fragments were purified with a MEGAquick-spin Total Fragment DNA Purification Kit (iNtRON Biotechnology, Gyeonggi-do, South Korea). The resulting nucleotide sequences were assembled and nucleotide (nt) and amino acid (aa) sequences were aligned and pairwise compared by using Geneious version 9.0.5 (Biomatters, Auckland, New Zealand). Phylogenetic trees were constructed with the Mega 7.1 software (53) .

Data regarding the presence of BTV in the samples examined by BTV RT-qPCR between 2013 and 2018 are shown in Table 1 . Table 2) .

Of the 792 BTV RT-PCR samples inoculated into ECE, 136 BTV strains were isolated, 15 were isolated from the batch of samples 

In this study, BTV-3 was detected in 20 farms numbered from 1 to 21 (Tables 3, 4 The main clinical signs in sheep included pyrexia, nasal discharge, hyperemia and ulceration of the oral and/or nasal mucosa, facial and/or thorax and hindlimbs edema, fatigue, apnea, recumbency, and few abortions. Morbidity in flocks was seen both in young and adult animals and ranged from 5 to 33.3%, with case mortality ranging from 0 to 30% (Table 3) . In farms number 1, 2, 6, and 9, no additional pathogen was identified, while in farms number 3, 4, 5, 10, and 11 additional BTV serotypes were contemporaneously identified or/and isolated in the farm or/and in the same animals (BTV-4, 6, 8, Tables 2, 3 were seen in this farm during 2 month period and confirmed by successful virus isolation ( Table 4 ).

In 2018, BTV-3 was also detected in two goat farms (Tables 3, 4) . One was detected in the spleen of a sudden dead goat and the other in a blood sample collected from a doe right after abortion. In the first case, the low BTV-3 load (Ct 32) found in the spleen was the only laboratory finding as neither bacteriological nor toxicological investigation was conducted. Regarding the farm with abortion cases, BTV-4 was also detected in the blood of another doe after abortion. In both farms the pathogen was detected in adult does.

BTV-3 was also identified in field samples collected from 10 cattle farms situated in southern and central areas of the country (Figure 1; Tables 3, 4) . Except for the cattle farm number 14, where clinically healthy animals were tested for commercial purposes, BTV-3 was detected in samples of sick and dead cattle received from the beginning of September till the middle of December 2018. In two cattle farms (farms 8 and 15) with unusual mortality rates, Clostridium perfringens and Salmonella spp. were also found in calves, respectively (Tables 3, 4) . In farm number 18, a post parturient BTV-3 positive adult cow was euthanized 2 days after showing hypersalivation, neural ketosis, and neck tense muscle. In addition, a one-year-old calf died after suffering from fever, fatigue, and conjunctival hyperemia. From blood sample of this calf, BEFV and BTV-3 were detected and BTV-3 was successfully isolated. In farm number 20, the clinical cases associated with BTV-3 detection included a calf which died after showing foamy salivation, nystagmus, hypothermia, and mucosal cyanosis and an adult cow showing post parturient ketosis, sharp blindness, endometritis, and tachycardia, which recovered after treatment with antibiotics and anti-inflammatory drugs. BTV-4 was also identified by RT-qPCR in tested samples. Similarly, BTV-3 and BTV-4 mixed infection was detected in a cow with bloody-purulent nasal discharge in farm number 21. Oppositely to previously described cases, no additional pathogens other than BTV-3 were found in farms number 13 and 19. In farm number 13, two adult pregnant cows showed recumbency before dying several days after appearance of clinical signs. No additional laboratory investigations were done for identifying the reason of the death. In farm 19, BTV-3 was detected in an adult cow which recovered after showing fever, indifference, inappetence, and sharp decrease in milk production. Field samples from ill cattle were not investigated for BTV-15 by RT-qPCR, due to absence of in-house or commercial RT-qPCR validated with currently circulating Israeli strains. Table S2) were almost full sequenced, when most segments of ISR-2396/2/17 strain were sequenced partially. Partial sequences of segment 2 (Seg-2) of all other Israeli BTV-3 isolates were also submitted to NCBI GenBank (accession numbers MH107823, MH107824, MN398282-MN398287) (https://www.ncbi.nlm. nih.gov/nucleotide/). Israeli BTV-3 strain sequences were also compared and data on nt and aa substitutions were summarized in the Table 5 . Due to very close relationship (99.45-99. 92% of nt identity in all genome segments) between BTV-3 strains isolated during 2017-2018 and ISR-2153/16 strain isolated in 2016, only sequences of the ISR-2153/16 strain were considered in the genetic analyses when nt sequences were compared to those of the global BTV strains ( Table 6) . First Israeli BTV-3 isolate (ISR-2019/13) was used as a prototype Israeli BTV-3 strain. Based on the data on number nt and aa substitutions, as well as on the identity with local or global BTV sequences, genome sequences were considered as homologous, or reassorted ( Table 5) .

According to BLAST and pairwise analyses, the ISR-2019/13 and the ISR-2262/2/16 strains were closely related (99.59% nt and 99.61 aa identity, Table 5 ). They also showed high degree of nt sequence identity with Tunisian (TUN2000/01) and French (FRA2001/01) BTV-2 strains, with which they shared 97.63-97.71% identity, respectively. The ISR-2153/16 strain was instead most closely related to the BTV-22 84/184 strain, with which it shared 96.66% nt identity (Table 6, Figure S1A ). 

The pairwise analysis showed that all isolated BTV-3 Israeli strains were closely related, sharing 99.39-99.52% nt and 99.47-100% aa identity ( Table 6) .

The phylogenetic analysis showed that Israeli and Tunisian BTV-3 belonged to the same western topotype/lineage; based on Seg-2 analyses, two different subclusters of closely related BTV-3 isolates circulating in Africa and the Mediterranean region were revealed. One included the closely related BTV-3 Israeli and the BTV-3 TUN-Zarzis/2016 isolates, the second subcluster included BTV-3TUN2016 and the BTV-3 ZIM2002/01 strain from Zimbabwe (Figure 2A) .

The BLAST and the pairwise analyses showed that the BTV-3 ISR-2019/13 and the ISR-2262/2/16 strains had almost identical nt sequences (99.75% nt and 100% aa identity, Table 5 ). They showed also high identity (97.73-97.83%) with the South African BTV-3 strain 3/O. aries-tc/ZAF/2017/Smithfield_VR33 ( Table 6 ). The ISR-2153/16 and South African BTV-5 O. ariestc/ZAF/2011/Benoni_01012015 strains were most closely related to each other, sharing 98.41% of nt identity ( Table 6) .

The phylogenetic analysis of Seg-3 nt sequences also revealed two sub-clusters. Figure S1B ).

The BLAST and pairwise analyses showed a very high nt sequence identity (98.93%) between the BTV-3 ISR-2019/13 and TUN2016/Zarzis strains. The ISR-2153/16 strain was most closely related to Chinese, South African, and Nigerian BTV-1 and BTV-16 (strains SZ97/1, NIG1982/10, and 5012, respectively), with which it shared 95.39-95.74% nt identity. The ISR-2262/2/16 strain was most closely related to Israeli BTV-24 (the representative strain ISR2008/02), with which it shared 99.39 % nt identity (Table 6, Figure S1C) .

Except for the ISR-2153/16 -like isolates, which formed a separate branch, the phylogenetic analysis showed that Tunisian and the other Israeli BTV-3 strains were closely related and formed a single cluster with some BTV-2,−4, and−24 strains that circulated in the Mediterranean Basin between 2000 and 2010. The phylogenetic analysis also confirmed the very high nucleotide sequence identity between the ISR-2019/13 and the TUN-Zarsis/2016 strains, which clustered together, and between the ISR-2262/2/16 and Israeli BTV-24 strains (Figure S1D) .

The BLAST and the pairwise analyses confirmed that all Israeli BTV-3 strains had high nt and aa sequences (99.64-99.71% and 99.47-100% respectively, Table 5 ). They also showed high nt (98.19-98.31%) identity with BTV-1 and BTV-4 strains from the Mediterranean Basin ( Table 6 ). The phylogenetic analysis of Seg-5 resulted in a single cluster grouping Israeli BTV-3 strains and BTV-1 and BTV-4 strains circulating in the Mediterranean Basin between 2006 and 2010 ( Figure S1E) .

The pairwise analysis showed that the three Israeli BTV-3 strains were closely related, sharing 99.82-99.94% nt and 99.61-100% aa identity ( Table 5) . According to the BLAST and the pairwise analyses, sequences of the Israeli strains showed up to 97% identity with the Nigerian BTV-16 NIG1982/10 strain ( Table 6 , Figure 2B ). Phylogenetic analysis of Seg-6 showed that the Israeli BTV-3 strains formed a separate branch which clustered with Nigerian BTV-16 (NIG1982/10), ZIM2002/01 BTV-3 strain from Zimbabwe, German and Israeli BTV-6 strains (GER2008/05 and ISR-2095/3/17, respectively) and Tunisian BTV-3 TUN2016 and TUN2016/Zarzis strains (Figure 2B) .

The BLAST and the pairwise analyses showed that Israeli BTV-3 strains had almost identical nt and aa sequences (99.65-99.73% and 99.41-100%, respectively, Table 5 ). They also showed highest identity with the South African BTV-3 3/O. aries-tc/ZAF/2017/Smithfield_VR33 strain (97.93-98.17%) ( Table 6 ). Phylogenetic analysis confirmed the close relationship between Israeli BTV-3 and South African BTV-3 3/O. ariestc/ZAF/2017/Smithfield_VR33 strains ( Figure S1E ).

The BLAST and the pairwise analyses also showed a very high nt and aa sequence identity between Israeli BTV-3 strains (99.14-99.82% and 98.85-100%, respectively, Table 5 ). A high nt identity (97.30-97.61%) was also observed with the Libyan BTV-1 LIB2007/06 strain ( Table 6 ). Phylogenetic analysis showed that Israeli BTV-3 strains formed a separate branch which clustered with Libyan BTV-1 strain LIB2007/06, BTV-2 French and Tunisian strains (TUN2000/01 and FRA2001/03, respectively), BTV-4 SUD1983/01 strain from Sudan and BTV-5 South African strain 5/ZAF/O.aries-tc/ZAF/2011/Benoni_01012015 ( Figure S1F ).

The pairwise and the BLAST analyses showed a very high nt and aa sequence identity (99.33-99.39%, respectively) between BTV-3 ISR-2019/13 and ISR-2153/16 strains ( Table 5) . They also showed high nt sequence identity with the Israeli BTV-24 ISR2008/02 strain (98.76-98.95%, Table 6 ). Conversely, the BTV-3 ISR-2262/2/16 strain was most closely related to the BTV-2 FRA2001/03 French strain with which it shared 97.13% nt identity (Tables 5, 6) .

Phylogenetic analysis showed that Israeli BTV-3 ISR-2019/13 and ISR-2153/16 strains formed a separate sub-cluster with the Israeli BTV-24 2008/02 strain. The ISR-2262/2/16 strain instead clustered with the British BTV-8 8UKG2007/06 strain ( Figure S1G ).

The BLAST and pairwise analysis of Seg-10 nt sequences showed that the BTV-3 ISR-2019/13 strain was closely related to Israeli BTV-5, 16 and 24 (strains ISR-1405/11; ISR-1794/4/12; ISR-3027/6/10, respectively), and to Cypriot BTV-16 CYP2006/01 strain, sharing 99.62-99.87% nt identity ( Table 6 , the only CYP2006/01 strain is shown). The BTV-3 ISR-2153/16 strain was closely related to some European BTV-8 (UKR2007/06 and NET2006/04 strains), BTV-12 O.ariestc/ZAF/2017/Queenstown_VR55and South African BTV-18 (BT32/76 and 600578 strains), sharing 97.13% nt identity (only BTV-12 12/ZAF/O.aries-tc/ZAF/2017/Queenstown_VR55 strain is shown in Table 6 ). The BTV-3 ISR-2262/2/16 strain had the closest relationship with BTV-19 (600579 and RA 19 OP strains), sharing 91.88-92.47% nt identity ( Table 6) .

The phylogenetic analysis showed that the ISR-2019/13 strain clustered with local (Israeli and Cypriot) BTV 

The first recorded cases of BTV-3 in the Mediterranean Basin date back in the middle of 20th century, when BTV-3 was isolated from sheep samples in 1943 and 1958 (www.ncbi.nlm.nih.gov/ pubmed). For the next several decades, no evidence of BTV-3 circulation was reported. BTV-3 circulation has been recently reported in Tunisia, Egypt and Italy. Even though detected either in goats or cattle (7, 38, 54) in these countries clinical signs associated to BTV-3 infection were observed and described in sheep only (7, 38, 39) . According to the retrospective findings of this study, during 2013-2017 BTV-3 outbreaks were restricted to Negev Desert area, in the southern district, while in 2018 it spread among most coastal and some central areas of Israel (Figure 1, Table 3 ). An additional difference between 2018 and 2013-2017 BTV-3 outbreaks was the exposure of Israeli goats and cattle in 2018 outbreaks. In this year, in fact, BTV-3 infection was occasionally associated with BT typical and atypical clinical signs, occurrence which was never observed before either in Israel or in the Mediterranean region (7, 38, 54) . These events may imply an increased infectivity and pathogenicity of the BTV-3 Israeli strains circulating in 2018 or/and their adaptation to the local vectors. However, because of the exiguous number of samples from diseased and clinically healthy animals received from the BTV-3 affected farms, it was hard to estimate the real exposure of ruminants and real effect BTV-3 infection on these animals. In Israeli goats and cattle, making any conclusion regarding the role of BTV-3 in causing illness and death was even harder due to involvement of different bacterial and viral agents in the clinical cases or to the absence of bacterial and toxicological investigations in some fatal cases. According to what observed in this survey, all Israeli BTV-3 strains were definitely capable of inducing classical manifestations of BT disease in sheep. These symptoms were more severe in case of BTV mixed infections (BTV-3 and BTV-4, BTV-6, or BTV-8). Notably, cases of acutely affected animal showing classical BT clinical signs were seen in one sheep farm during two-month period and confirmed by successful virus isolation. This fact illustrated presence of newly infected animals and probably infected vector during prolonged period during two seasons 2017 and 2018 (farms number 4 and 9, Table 4 ).

Interestingly, in some cattle and sheep farms, BTV-3 as well as BTV-4 were detected in fetal tissues, placenta and in newborn animals. Although these findings were not sufficient to establish the definitive involvement of BTV-3 in determining reproductive failures, its presence in fetal tissues and/or newborn animals, provided evidence of its capability to infect placenta and probably also cross the placental barrier. As far as we know, apart from BTV-8, causing numerous abortions and malformations, and lastly BTV-1, this capability is proper of vaccine or lab derived strains (18, 55) .

Higher susceptibility observed in sheep reflected the number of successful virus isolations achieved in this species (33%) compared to cattle (7.7%), which, in turn, may imply a higher number of sheep acutely infected by BTV-3 than cattle.

Two different strains named TUN2016 and TUN2016/Zarzis have been identified as responsible of the Northern African and Italian BTV-3 outbreaks (7, 38, 39, 54) . This study revealed that three additional BTV-3 strains have been circulating in the Mediterranean Basin and in Israel, in particular, at least since 2013. Genome comparison allowed a tentative reconstruction of the ancestral viral genome of these strains. When the Israeli BTV-3 strain sequences were compared, it was evident that Seg-2, 5, 6, 7 and 8 of the Israeli strains derived from a common ancestor. According to Seg-2 phylogenetic analyses, the BTV-3 Israeli strains also shared common ancestors with the Tunisian TUN2016/Zarzis strain. For the ISR-2019/13 strain, this was evident for Seg-4 too.

In all these strains, reassortment phenomena were also found. When compared with the prototype ISR-2019/13 strain, the ISR-2262/2/16 and the ISR-2153/16-like strains have 3 and 4 reassorted segments, respectively. Moreover, their Seg-4, Seg-9, and Seg-10 sequences clearly evidenced a different origin. "Traces" of untyped South African strains and BTV serotypes (BTV-18, BTV-19, BTV-22) exotic for the region in the Israeli BTV-3 genomes were clear, at least in their last segments. These results may indicate the circulation of exotic strain/serotypes in the region.

Between 2013 and 2017, only a low proportion of the BTV strains isolated from sick domestic ruminants was identified as BTV-3. Thus, in 2013 only one out of 15, in 2016, 3 out of 58, in 2017 4 out of 40 virus isolates were BTV-3. In these years, the BTV-3 circulation was limited to southern Israel only. In 2018, the number of BTV-3 isolates among the total number of BTV isolated, sharply increased (8 out of 25). The BTV-3 circulation also spread along all coastal area of Israel, suggesting an increased infectivity of the BTV-3 ISR-2153/16 strain among susceptible Israeli domestic ruminants which can be explained by point mutations both in coding and non-coding regions (not shown in this work) or possible introduction closely related viruses to ISR-2153/16 strain into Israel.

In spite of a long history of BTV-3 infections in ruminants in South Africa, India, Caribbean and Northern America (41, 42, 56, 57) , epidemiological data as well as information on pathogenicity and infectivity of BTV-3 infection from these regions are absent or scarce. Even if it was not possible to evaluate their pathogenicity in cattle and goats, this paper still gives important information on the possible origin of the BTV-3 strains circulating in the Mediterranean basin, elucidating their possible routes of incursion. However, further investigation is needed to improve our understanding on their biological properties and their impact on livestock.

The data analyzed in this study can be found in the article/Supplementary Material.

This study did not involve any human participants and animals. Ethical approval was not required for this study according to national legislation.

GK, JV, IS, SS, AM, ID, IG, AK, OK, and EO collected field samples and epidemiological data. NG, DD, and VB performed diagnostic tests of field samples. NG developed in-house PCR systems, performed PCR tests for sequencing, analyzed and summarized data of sequencing, and phylogenetic analyses. VB, AG, and NG performed isolations of the viruses in eggs and tissue cultures. AE and AL reviewed the manuscript. NG and GS drafted the manuscript, while VB, AE, and YK discussed the results and commented on the manuscript.

This study was supported by H2020 PALE BLU project, Proposal ID 727393-2. 

Seg-4

Seg-5 JX861492

Seg-5 JX861492/BTV-1/FRA2007/18*-98

Seg-4

Seg-5 JX861492

*Strain sharing the same nt identity with several other BTV strains belonging to BTV-1 and BTV-4

The pathology and pathogenesis of bluetongue

Bluetongue virus, other orbiviruses and other reoviruses: their relationships and taxonomy

Characterization of a second open reading frame in genome segment 10 of bluetongue virus

Bluetongue virus: from BTV-1 to BTV-27

Characterization of bluetongue virus serotype 28

Analysis of bluetongue serotype 3 spread in Tunisia and discovery of a novel strain related to the bluetongue virus isolated from a commercial sheep pox vaccine

Serological and Genetic Characterisation of Putative New Serotypes of Bluetongue Virus and Epizootic Haemorrhagic Disease Virus Isolated from an Alpaca

Emergence of a novel Bluetongue virus serotype

Novel putative Bluetongue virus in healthy goats from

One after the other: a novel Bluetongue virus strain related to Toggenburg virus detected in the Piedmont region (North-western Italy), extends the panel of novel atypical BTV strains

Culicoides and the emergence of bluetongue virus in northern Europe

Climate change and the recent emergence of bluetongue in Europe

Bluetongue: history, global epidemiology, and pathogenesis

Experimental infection of European red deer (Cervus elaphus) with bluetongue virus serotypes 1 and 8

Transplacental and oral transmission of wild-type bluetongue virus serotype 8 in cattle after experimental infection

Using shared needles for subcutaneous inoculation can transmit bluetongue virus mechanically between ruminant hosts

Bluetongue serotype 2 and 9 modified live vaccine viruses as causative agents of abortion in livestock: a retrospective analysis in Italy

Evidence for transmission of bluetongue virus serotype 26 through direct contact

Bluetongue virus serotype 27: experimental infection of goats, sheep and cattle with three BTV-27 variants reveal atypical characteristics and likely direct contact transmission BTV-27 between goats

Epidemiologic characteristics of bluetongue virus serotype 8 laboratory-confirmed outbreaks in the Netherlands in 2007 and a comparison with the situation in 2006

Mortality and case fatality during the recurrence of BTV-8 in northern Europe in 2007

Bluetongue virus serotype 24 (BTV-24) in Israel: phylogenetic characterization and clinical manifestation of the disease

Unusual clinical manifestations in Israeli ruminant populations infected with Orbiviruses

Viral replication kinetics and in vitro cytopathogenicity of parental and reassortant strains of bluetongue virus serotype 1, 6 and 8

Reassortment between two serologically unrelated bluetongue virus strains is flexible and can involve any genome segment

Widespread reassortment shapes the evolution and epidemiology of bluetongue virus following european invasion

Genomic sequences amongst bluetongue viruses revealed by comparisons of capsid and outer coat protein nucleotide sequences

Relationships amongst bluetongue viruses revealed by comparisons of capsid and outer coat protein nucleotide sequences

Genetic diversity of bluetongue viruses in south east Asia

Cross-neutralization of genetic reassortants of bluetongue virus serotypes 20 and 21

The use of recombinant DNA technology for the development of a bluetongue virus subunit vaccine

Analysis of the roles of bluetongue virus outer capsid proteins VP2 and VP5 in determination of virus serotype

Analysis and phylogenetic comparisons of full-length VP2 genes of the 24 bluetongue virus serotypes

Design of primers and use of RT-PCR assays for typing European bluetongue virus isolates: differentiation of field and vaccine strains

Phylogenetic analysis of bluetongue virus genome segment 6 (encoding VP5) from different serotypes

Sequence comparison of the L2 and S10 genes of bluetongue viruses from the United States and the People's Republic of China

Western Bluetongue virus serotype 3 in Sardinia, diagnosis and characterization

A novel Bluetongue virus serotype 3 strain in Tunisia

Bluetongue virus serotype 3 in Western Sicily

A South African overview of the virus, vectors, surveillance and unique features of bluetongue

Recent US bluetongue virus serotype 3 isolates found outside of Florida indicate evidence of reassortment with co-circulating endemic serotypes

Abstracts of the 11th EPIZONE Annual Meeting

Emergence of a Novel Reassortant Strain of Bluetongue Serotype 6 in Israel, 2017: Clinical Manifestations of the Disease and Molecular Characterization

Clinical Veterinary Microbiology

High-resolution melting (hrm) for genotyping bovine ephemeral fever virus (befv)

Development of a multiplex real-time PCR for detection and differentiation of malignant catarrhal fever viruses in clinical samples

Real Time RT-PCR for the detection and quatitation of bovine respiratory syncytial virus

Development of one-run real-time PCR detection system for pathogens associated with bovine respiratory disease complex

Development and validation of a universal S-segment-based real-time RT-PCR assay for the detection of Simbu serogroup viruses

A disease similar to Blue Tongue in cattle and sheep in Israel. Ref Vet

Epizootic hemorrhagic disease virus serotype 6 outbreak in Israeli cattle in 2015

MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets

Presence of bluetongue and epizootic hemorrhagic disease viruses in Egypt in 2016 and 2017

Isolation of bluetongue virus serotype 1 from aborted goat fetuses

Reverse transcription loop-mediated isothermal amplification assays for rapid identification of eastern and western strains of bluetongue virus in India

Molecular epidemiology studies of bluetongue virus