The aim of the present study is to characterize the dipteran larvae infesting sheep at SRRC, Mannavanur by means of molecular phylogenetic analysis.
The area of the present study falls under subtemperate zone, which is located about 2030 metres above mean sea level and is receiving an annual rainfall of 1055 mm throughout the year. The grazing area of sheep is invaded by grass species including Kikuyu ( Pennisetum clandestinum) and speargrass (Heteropogon contortus) (Swain et al.2004).
At SRRC, Mannavanur, the cases of sheep infested with Oestrus ovis larvae are mild and asymptomatic and the larvae in the nasal sinuses of sheep are noticed only upon post-mortem inspection. So far, the typical systemic cases of Oesrosis have not been observed among sheep at SRRC, Mannavanur. Numerous factors including the host susceptibility, chrobinoology of O. ovis in a particular geographical region and system of animal management would be influencing the responses of the host to oestrosis and larval burden (Sotiraki and Hall 2012). Very recently, Bello et al. (2022) concluded that there is no correlation between the larvae infestation intensity & the clinical signs of Oestrosis and ELISA is better when compared to PCR technique for the diagnosis of O.ovis in live animals.
In addition, the oral administration of commercial preparation of Ivermectin (HITEKR : Each ml contains 0.800 mg of Ivermectin : Dose – 25 ml/100 Kg body weight) for the treatment of gastrointestinal parasitism among sheep at SRRC, Mannavanur, is a periodical practice and this could be the reason for the absence of typical clinical cases of Oestrosis among sheep in the present study. Bello et al. 2022) reported that in sheep, oral drenching of with Closantel at a dose rate 10 mg/kg and subcutaneous injection with ivermectin at a dose rate of 0.2 mg/kg could be effective in the control of oestrosis.
As per the published report (Soulsby 1987), the dipteran larvae in the nasal sinuses of sheep at SRRC, Mannavanur was identified as the maggot of O.ovis based on the possession of 10 segments with a dark transverse band on the dorsal surface of each segment and triangular posterior spiracle with radiating slits (Fig. 1). Additionally, the maggots were also onfirmed as the third stage larvae of O.ovis as there were spines on the ventral aspect (Allaie et al.2016). As per the findings of Metwally et al. (2021), it was found that third instar larvae (L3) of O.ovis in Sheep and Goats from Riyadh, Saudi Arabia would be yellow in colour while immature stage and light brown in colour while mature stage. The larvae were consisting of broad transverse blackish bands on the dorsal side. There was a presence of variable number of small denticles on the dorsal side of the second segment .The following segments were plain with a rough, leather like skin pattern, which was clearly observed only on the darkened parts. The ventral side of the larvae was characterized by the presence of series of strong spines on the segments and the spines were haphazardly located on the segment but the spines were following a regular pattern in the succeeding segments. The preanal bulge was naked where as the postanal bulge were having fewer spines. Further, the posterior peritremes were circular in shape having a central button and devoid of a clear-cut suture. There was a channel in the posterior peritremes of O.ovis larvae. It is therefore recommended that it is essential to carry out a detailed morphological examination of the larvae of Oestrus ovis in sheep from SRRC, Mannavanur in the future.
Further, there are three stages in the development of Oestridae larvae. The most frequently observed larvae by veterinary inspectors at slaughterhouses are L3 larvae. It is very easy to separate the internal tissues of L3 from the chitin rich external cuticle and therefore, there would be a huge quantity of larval tissues available for genomic DNA isolation. Additionally, when compared to first (L1) and second (L2) stage larvae, fully developed L3 would be having low amount of enzymes which are having the inhibitory activities in the downstream processing of the genomic DNA ( Otranto et al.2003). In the present study also, the yield and quantity of the genomic DNA extracted from the sheep maggots were found to be good and eventually, the amplicons obtained out of PCR experiments were also having a single intense band ( around 1600 bp size) while visualizing the agarose gel under UV light after running electrophoresis. It is therefore concluded that in the present study, the stage of dipteran larvae utilized for genomic DNA extraction and PCR experiments could undoubtedly be the third stage larvae of O.ovis.
Insect mitochondrial cytochrome oxidase I (COI) genes are known for their huge size, carriers of both conserved and variable regions and the hot spots for mutations. Owing to the said genetic features, COI genes would be beneficial in the accurate calculation of the rate of evolution among the insect species (Lunt et al.1996). Therefore, COI gene based PCR was used for the molecular data analysis of the dipteran larvae collected from the nasal sinuses of sheep in the present study.
Since there is no published report about the molecular characterization of O.ovis from Indian sheep, the present study is the first report of its own kind.
As per the GenBank Accession No. NC_059851.1 representing the complete mitochondrial genome of O.ovis, in the present study, the size of the PCR amplified fragment specifying cytochrome c oxidase subunit I (COI) gene of O.ovis from Sheep at SRRC, Mannavanur, Tamil Nadu, India is 1534 bp and it is nothing but the full length nucleotide sequence of COI gene for O.ovis (covering the full length coding sequence of COI gene of O.ovis – from 1499 bp to 3032 of GenBank Accession No. NC_059851.1).
As reported earlier [16], in the present study, the COI gene of O.ovis was also having high AT content (69%) (A: 30% 475; T: 39% 568; G: 13% 212; C: 18% 279). In a recent study on molecular analysis of O.ovis larvae in Small ruminants from Saudi Arabia, it was found that the overall GC content of 23 mtCOI gene sequences was 35.6% (i.e., 64.4% AT content) (Metwally et al.2021).
While looking for the type II restriction enzyme sites within the nucleotide sequences of both COI gene of O.ovis from Mannavanur ( South Indian isolate ; GenBank Accession No. ON000070.1) & Jammu (North Indian isolate; GenBank Accession No. ON912055.1) using the GenScript Restriction Enzyme Map Analysis Tools (free software available in the public domain), it was found that the common hexabase cutter (Sac I – GAGCTC) was having only one cutting site at the position 857/862 in the former sequence, where as the latter sequence was having two cutting sites for Sac I at positions 257/262 & 857/862. Similarly, the COI gene sequences of Mannavanur isolate was having one cutting site for the restriction enzymes, EcoRV (GATATC) and Xho I (CTCGAG) at positions 1181/1184 and 896/897, respectively but no recognition site for the said two six base cutters was present in COI gene sequences of Jammu isolate from India. This interesting feature of restriction enzyme analysis could be explored for the differentiation of large number of specimens of O.ovis isolates from different geographical areas of India.
The COI gene sequences of Oestrus ovis from Sheep at SRRC, Mannavanur, India, obtained out of sequencing experiments by M/s. Eurofins Genomics India Pvt. Ltd., Bengaluru-560048, Karnataka, India, were compared to that of 29 different dipteran species available in the NCBI database. In the present study, the alignment of last 90 amino acids at carboxyl terminal ranging from the amino acid position 421 to amino acid position 510 of COI protein of Oestrus ovis from Mannavanur with that of other dipteran species is shown in Fig. 2. The above said amino acid alignment region is falling within the hypervariable part of the COI protein corresponding to the region from external loop 4 (E4) to the carboxyl terminal (–COOH) of myiasis-causing Oestridae (Otranto et al.2003).
Similar to the findings discussed in earlier report (Otranto et al.2003), in the present study also, the amino acid alignment revealed that the presence of five most conserved residues in COI protein within the subfamily of Oestrinae (i.e. K-405, L-450, I-464, F-493 and N-510) was noticed. But in case of Rhinoestrus usbekistanicus (GenBank Accession No. NC_045882.1), the amino acid Lysine (K) was replaced by the amino acid Serine (S) at position 405 (Supplementary Fig. 1).
The sequence analysis of the current study revealed that the recently published Ostrus ovis sequence from Jammu has 99% sequence identity with our sequence. At position 118, amino acid Thr appears in Oestrinae subfamily compared to Lys in Calliphora, Haematobia and Stomoxys family, whereas it was replaced by Asn in other members included in this study. Similarly, at position 253, Ser was replaced by Pro in Oestrus ovis sequence published from Jammu compared to other sequences that were included in this study, which were conserved with Ser. At position 268, Ser was found to be conserved in all reported Oestrus species, however, it was substituted with Ala in a few species like Rhinoestrus, Cephenemyia, Muscina, Gasterophilus and Culicoides, whereas it was replaced with Thr in Phlebotomus. Similarly, at position 332, amino acid Tyr/Ser was replaced by Leu in Culicoides species and Ostrus ovis sequence of Jammu, whereas it was replaced by Phe in Rhinoestrus and Cephenemyia species. At positions 335 and 339, either Ala/Thr was replaced by Ser or vice versa. At position 483, the residue Pro was replaced by Thr in all the Ostrus ovis sequences, whereas it was replaced by Ser in case of Rhinoestrus sps. At position 499, Ala was replaced by Ser in sequence of Oestrinae subfamily members, Glossina, Hhypoderma and Gasterophilus species, whereas in case of Muscina by Thr (Supplementary Fig. 1).
As per the published report (Schwartz et al.2019), the amino acid Ser substituted with either Pro or Ala occur only in the conserved position, whereas Ser substituted with either Asn or Gly occur at less conserved position. According to the above said report, the occurrence of Ser with Ala substitution (in most cases) at position 44, 253, 268, 335, 339 and 499 could highly be conserved, where as substitution of Ser with Asn or vice versa at position 487 might be less conserved. However, substitution of Pro to Thr/Ser or vice versa at position 483 could be important for functional phosphorylation (Supplementary Fig. 1).
While searching for the functional motifs within the protein encoded by COI gene of O.ovis using the Scanprosite software, the presence of cytochrome oxidase subunit I signature and profile was observed in the region ranging from 1–510 (Fig. 3). Apart from the family signature, the protein sequence also possesses Heme-copper oxidase catalytic subunit, copper B binding region signature (WFFGHPEVYILILPGFGMISHIISQESGKKETFGSLGMIYAMLAIGLLGFIVWAHH ) in the region ranging from 234–289 (Fig. 4) .
As far as the sequence analysis of 30 COI gene sequences pertaining to different dipteran spcies in the present study is concerned, it was observed that COI gene of O.ovis from sheep at Mannavanur, Tamil Nadu, India (South India), exhibited 100% sequence identity both at nucleotide and amino acid levels with that of O.ovis from Spain (GenBank Accession No. NC_059851.1). O.ovis from Jammu, India (North India) (Gen Bank Accession No. ON912055.1) shared 92 and 99 percent identity with O.ovis from Mannavanur at nucleotide and amino acid level, respectively. With other members of Oestrinae, the range of nucleotide and amino acid identity of O.ovis from South India was 85–86% and 95–96%, respectively. Hypoderma lineatum (GenBank Accession No. NC_013932.1) was having 81% identity at nucleotide level and 90% at amino acid level with O. ovis from South Indian. On the other hand, with Gasterophilus intestinalis (GenBank Accession No. NC_029834.1), the South Indian isolate of O.ovis was sharing 80% and 90% at nucleotide and amino acid levels, respectively. The COI gene of the members belonging to the subfamilies Caalliphorinae, Chrysomyinae and Lucilinae members were displaying 85 & 92% identity at nucleotide and amino acid levels with that of South Indian isolate of O.ovis (Table 1).
Phylogenetic relationship was constructed employing MEGA-X on the multiple sequence alignment of COI amino acid sequences of 30 dipteran flies. In the phylogenetic analysis, the nucleotide sequences having the size of more than 1000 bp of COI gene of the dipternan flies (including faculatative myiasis causing flies as well as biting flies) were only retrieved from the NCBI database and some of the nucleotide sequences were full length COI gene. There were only three number of nucleotide sequences (more than 1000 bp) representing the COI gene of O.ovis in the GenBank database.
Other than Cuterebrinae, the COI gene sequences (more than 1000 bp) belonging to the members of the other three subfamilies Oesstrinae, Gasterophilinae and Hypodermatinae were included in the molecular analyses of the present study. Since the nucleotide sequences representing the COI gene of members of the subfamily Cuterebrinae available in the NCBI database were less than 1000 bp, no COI gene sequences of the said subfamily within Oestridae family was included in the molecular phylogenetics of the present study.
Phylogenetic analysis based on amino acid sequences ruled out the close relationship between O.ovis from South India and other oestrid flies within the subfamily Oestrinae from different geographical area across the world (Fig. 5).
Based on the phylogenetic tree, it was found that COI gene could clearly differentiate the species in all three subfamilies (Oestrinae, Hypodermatinaedae and Gasterophilinae of estridae (except Cuterebrinae) (Fig. 5 ). Further, the branches of the subfamilies Gasterophilinae and Hypordematinae were forming the clusters and divided from the subfamily Oestrinae (Fig. 5). This finding is also consistent with the earlier reports on Oestrinae species from Turkey (Karademir et al.2020), Spain (Moreno et al.2015) and Italy ( Otranto et al.2003), where in, partial COI nucleotide sequences were used in the phylogenetic analysis and there was a clear cut cluster formed by the members of the subfamily Oestrinae & deviating from the members belonging to the other three subfamilies within the family Oestridae. Interestingly, the members of the subfamilies, Calliphorinae, Chrysomyinae and Lucilinae within the family Calliphoridae were grouped and clustered separately in the present phylogentic tree analysis (Fig. 5) and is agreeing with the recent report (Nasser et al. 2021).
In the present study, the positions of Oestrus ovis and Rhinoestrus usbekistanicus in the phylogenetic tree (Fig. 5) are supporting the findings of the earlier report ( Aleix-Mata et al. 2021), where in the said two Oestrinae species could be sister groups as per the phylogenetic tree analysis involving the whole mitochondrial genome of the 13 Oestridae species.
Even though O. ovis is the usual parasite of small ruminants across the continents, it is also known to infest other ruminants including Old World camels (Fekry et al.1997), New World camels ( Gomez-Puerta et al.2013), the Siberian ibex (Capra ibex sibirica), argali (Ovis ammon), bighorn sheep (O. canadensis), Barbary sheep (Ammotragus lervia) and the European mouflon (O. orientalis musimon) as accidental hosts (Moreno et al.1999). The grazing area of sheep SRRC, Mannavanur, is co-grazed by cattle, buffaloes and horses, which are belonging to the villagers of the local area. In addition to the domestic animal species, the area of the present study is also a habitat for wild animals such as Indian Guar, Deer, wild pig and red dogs (Nagarajan 2020). Owing to the co-sharing of both domestic and wild animals, the area of the present study is speculated to have the richness of dipteran flies (belonging to the superfamilies, Muscoidea and Oestroidea), which are responsible for facultative myiasis in the total fauna population.