The flies were mostly observed on the soft part of the body in the dewlap, perpetual flap, udder and peri-anal regions. Hard cutaneous nodules formation in the various body parts of 2-3 cm size were noted during the study. Similar findings of extreme annoyance and blood loss (up to 0.5 ml per fly) and dermal nodules (local reactions to bites) were recorded by
Foil and Hogsette (1994). The Tabanus flies (Fig 1) recorded from body surface of mithun were prominent size (17-18 mm), proboscis was stout, pointed downwardly projecting with large eyes and hyline wings. Thorax black with yellowish tomentum, black pubescence. Abdomen reddish yellow on first 3 segments with large black median spot on 2
nd and 3
rd segment, 4
th abdominal tergite blackish at its posterior border (Fig 2). Legs are black and fore pair of tibiae with black apex. Wing venation was prominent and R5 cell was open. Wing vein R4 and R5 were widely splayed either side of the wing tip with R5 converging on vein M1(Fig 3).
Tabanus flies’ infestation increased during May to August (premonsoon and monsoon) and moderately declines thereafter which may be due to the breeding season of the flies (Table 1). These periods are also characterized by higher humidity and rainfall, which create favorable breeding conditions for Tabanus flies. As the seasons progress and weather conditions change, the incidence of infestations usually declines. In a study in Himachal Pradesh on cattle animals,
Maity et al., (2015) also recorded abundant Tabanidae mainly throughout the year with a little decline in mid to higher elevations during winter.
Barros (2001) noted the tabanid peak, in October, While
Lucus et al., (2020) observed horse fly season in Tacuarembó, Uruguay started in September and ended in May.
Based on morphological findings followed by DNA barcoding with the PCR amplification of marker genes: ITS-2 and mitochondrial enzyme COX1(Fig 4), species were identified as
Tabanus rufofrater. The sequence analysis confirmed 99% similarity with
Tabanus rufofrater (DQ631993.1) and 98% similiarity with
Tabanus atratus (KU568511).
Phylogenetic analysis
The evolutionary history was inferred by using the Maximum Likelihood method based on the Kimura 2-parameter model as per
Kimura M. (1980). The tree with the highest log likelihood (-2316.98) is shown (Fig 5). The percentage of trees in which the associated taxa clustered together is shown next to the branches. Initial tree(s) for the heuristic search were obtained automatically by applying Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances estimated using the Maximum Composite Likelihood (MCL) approach and then selecting the topology with superior log likelihood value. The analysis involved 11 nucleotide sequences.Codon position sinc luded wer 1
st+2
nd+3
rd+Noncoding. All positions containing gaps and missing data were eliminated. There were a total of 648 positions in the final dataset. Evolutionary analyses were conducted in MEGA7 as described by
Kumar et al., (2015).
Distance matrix
The number of base substitutions per site from between sequences are shown (Table 2). Standard error estimate(s) are shown above the diagonal and were obtained by a bootstrap procedure (
Felsenstein, 1985) (1000 replicates). Analyses were conducted using the Kimura 2-parameter model. The analysis involved 11 nucleotide sequences. Codon positions included were 1
st + 2
nd + 3
rd + Noncoding. All positions containing gaps and missing data were eliminated. There was a total of 648 positions in the final dataset. Evolutionary analyses were conducted in MEGA7 as per
Kumar et al., (2016). There were a total of 648 positions in the final dataset. Evolutionary analyses were conducted in MEGA7.