Exploring genetic diversity and variation of Ovar-DRB1 gene in Sudan Desert Sheep using targeted next-generation sequencing.

INTRODUCTION: The Ovar-DRB1 gene, a crucial element of the Major Histocompatibility Complex (MHC) Class II region, initiates adaptive immunity by presenting antigens to T-cells. Genetic diversity in sheep, particularly in MHC Class II genes like Ovar-DRB1, directly influences the specturm of presented antigens impacting immune responses and disease susceptability. Understanding the allelic diversity of Ovar-DRB1 gene in Sudan Desert Sheep (SDS) is essential for uncovering the genetic basis of immune responses and disease resistance, given the the breeds significance in Sudan's unique environment. METHODS: Utilizing Targeted Next-Generation Sequencing (NGS) we explore allelic diversity in Ovar-DRB1 gene within SDS. Successfully ampliying and and sequencing the second exon of this gene in 288 SDS samples representing six breeds provided a comprehensive allelic profile, enabling a detalied examination of the gene's genetic makeup. RESULTS: We identifed forty-six alleles, including four previously unreported, enrichness the genetic diversity of SDS breeds. These alleles exhibiting non-uniform distribution, varying frequencies across breeds, indicating a breed-specific genetic landscape. Certain alleles, known and novel, show higher frequencies in specific populations, suggesting potential associations with adaptive immune responses. Identifying these alleles sets the stage for investigating their functional roles and implications for disease resistance. Genetic differentiation among SDS breeds, as indicated by FST values and clustering analyses, highlights a unique genetic makeup shaped by geographic and historical factors. These differentiation patterns among SDS breeds have broader implications for breed conservation and targeted breeding to enhance disease resistance in specific populations. CONCLUSION: This study unveils Ovar-DRB1 gene allelic diversity in SDS breeds through targeted NGS and genetic analyses, revealing new alleles that underscore the breeds' unique genetic profile. Insights into the genetic factors governing immune responses and disease resistance emerge, promising for optimization of breeding strategies for enhanced livestock health in Sudan's unique environment.

Intra- and interspecies variation and population dynamics of Fasciola gigantica among ruminants in Sudan.

Fasciola gigantica is a widespread parasite that causes neglected disease in livestock worldwide. Its high transmissibility and dispersion are attributed to its ability to infect intermediate snail hosts and adapt to various mammalian definitive hosts. This study investigated the variation and population dynamics of F. gigantica in cattle, sheep, and goats from three states in Sudan. Mitochondrial cytochrome c oxidase subunit I (COI) and NADH dehydrogenase subunit 1 (ND1) genes were sequenced successfully to examine intra and interspecific differences. ND1 exhibited higher diversity than COI, with 15 haplotypes and 10 haplotypes, respectively. Both genes had high haplotype diversity but low nucleotide diversity, with 21 and 11 polymorphic sites for ND1 and COI, respectively. Mismatch distribution analysis and neutrality tests revealed that F. gigantica from different host species was in a state of population expansion. Maximum likelihood phylogenetic trees and median networks revealed that F. gigantica in Sudan and other African countries had host-specific and country-specific lineages for both genes. The study also indicated that F. gigantica-infected small ruminants were evolutionarily distant, suggesting deep and historical interspecies adaptation.

Genetic variation and demographic history of Sudan desert sheep reveal two diversified lineages.

More than 400 million sheep are raised on the African continent, the majority of which are indigenous and are primarily reared for sustenance. They have effectively adapted to various climatic and production environments, surviving and flourishing. The genetic relationships among these sheep populations remain understudied. Herein, we sequenced the entire mitochondrial DNA control region of 120 animals from Hamary and Kabashi and their crossbreed (Hamary x Kabashi) of Sudan desert sheep (SDS) to understand their maternal-inherited genetic variation and demographic history profiles and relate those to the history of sheep pastoralism on the African continent. The results show a diversified and predominant D- loop haplogroup B (n = 102, 85%), with all other sequences belonging to haplogroup A. Most of the maternal genetic variation was partitioned between haplogroup (76.3%) while within haplogroup accounted for 23.7% of the variation. However, little genetic differentiation was observed among the two breeds and their crosses, with our results supporting a Hamari maternal origin for the crossbreed. Bayesian coalescent-based analysis reveals distinct demographic history between the two haplogroups, two breeds and their crosses. Comparison of the two haplogroup showed that haplogroup B experienced an earlier expansion than haplogroup A. Unlike the breed-based comparison, the expansion of the two breeds started roughly at the same time, around 6500 years ago, with Kabashi having a slightly greater effective population size. The maternal ancestors of SDS may have diverged before their introduction to the African continent. This study provides novel insights into the early history of these two main breeds of Sudan desert sheep and their crosses.

Genomic signatures for drylands adaptation at gene-rich regions in African zebu cattle.

BACKGROUND: Indigenous Sudanese cattle are mainly indicine/zebu (humped) type. They thrive in the harshest dryland environments characterised by high temperatures, long seasonal dry periods, nutritional shortages, and vector disease challenges. Here, we sequenced 60 indigenous Sudanese cattle from six indigenous breeds and analysed the data using three genomic scan approaches to unravel cattle adaptation to the African dryland region. RESULTS: We identified a set of gene-rich selective sweep regions, detected mostly on chromosomes 5, 7 and 19, shared across African and Gir zebu. These include genes involved in immune response, body size and conformation, and heat stress response. We also identified selective sweep regions unique to Sudanese zebu. Of these, a 250 kb selective sweep on chromosome 16 spans seven genes, including PLCH2, PEX10, PRKCZ, and SKI, which are involved in alternative adaptive metabolic strategies of insulin signalling, glucose homeostasis, and fat metabolism. CONCLUSIONS: Our results suggest that environmental adaptation may involve recent and ancient selection at gene-rich regions, which might be under a common regulatory genetic control, in zebu cattle.

Demographic expansion and high level of matrilineal diversity in two populations of East African Baggara cattle.

Western Baggara cattle breed (WBCB) is an East African zebu inhabiting Sudan, well-known as beef-producing cattle. We investigated herein two phenotypically and geographically distinct populations of this breed, namely Nyalawi and Daeinawi, which are renowned for their unique meat production capabilities and adaptation attributes, with the aim to contribute to our understanding of their maternal genetic diversity and demography dynamics. Genetic polymorphism analysis of the full-length D-loop mtDNA region revealed 44 and 35 polymorphic sites defining 28 and 24 distinct haplotypes in the Nyalawi and the Daeinawi, respectively. Observed genetic diversity is high within the population with a low level of genetic differentiation between populations. Approximate Bayesian computation via the calculation of Bayesian skyline plots and neutrality tests support past expansion with a higher maternal effective population size (Ne ) in Nyalawi compared with the Daeinawi population and a population expansion beginning around 4,500 YBP and 3,500 YBP, respectively, before the arrival of zebu into the continent.

Ancient and modern DNA reveal dynamics of domestication and cross-continental dispersal of the dromedary.

Dromedaries have been fundamental to the development of human societies in arid landscapes and for long-distance trade across hostile hot terrains for 3,000 y. Today they continue to be an important livestock resource in marginal agro-ecological zones. However, the history of dromedary domestication and the influence of ancient trading networks on their genetic structure have remained elusive. We combined ancient DNA sequences of wild and early-domesticated dromedary samples from arid regions with nuclear microsatellite and mitochondrial genotype information from 1,083 extant animals collected across the species' range. We observe little phylogeographic signal in the modern population, indicative of extensive gene flow and virtually affecting all regions except East Africa, where dromedary populations have remained relatively isolated. In agreement with archaeological findings, we identify wild dromedaries from the southeast Arabian Peninsula among the founders of the domestic dromedary gene pool. Approximate Bayesian computations further support the "restocking from the wild" hypothesis, with an initial domestication followed by introgression from individuals from wild, now-extinct populations. Compared with other livestock, which show a long history of gene flow with their wild ancestors, we find a high initial diversity relative to the native distribution of the wild ancestor on the Arabian Peninsula and to the brief coexistence of early-domesticated and wild individuals. This study also demonstrates the potential to retrieve ancient DNA sequences from osseous remains excavated in hot and dry desert environments.

First molecular characterization of Echinococcus granulosus (sensu stricto) genotype 1 among cattle in Sudan.

BACKGROUND: Echinococcus granulosus sensu lato (s.l.) is the causative agent of cystic echinococcosis (CE), which is a cosmopolitan zoonotic parasitic disease infecting humans and a wide range of mammalian species including cattle. Currently, little information is available on the genetic diversity of Echinococcus species among livestock in Sudan. In the present study, fifty (n = 50) hydatid cysts were collected from cattle carcasses (one cyst sample per animal) at Al-kadarou slaughterhouse, Khartoum North, Sudan. DNA was extracted from protoscolices and the germinal layer of each cyst and subsequently amplified by PCR targeting the mitochondrial NADH dehydrogenase subunit 1 (NADH-1) gene. The amplified PCR products were purified and subjected to direct sequencing for subsequent construction of phylogenetic tree and net work analysis. RESULTS: The phylogenetic tree revealed the presence of Echinococcus canadenesis genotype 6 (G6) in 44 cysts (88.0%), Echinococcus ortleppi genotype 5 (G5) in 4 cysts (8.0%) and Echinococcus granulosus sensu stricto (s.s) genotype 1 (G1) in 2 cysts (4.0%). The phylogenetic network analysis revealed genetic variation among the different haplotypes/genotypes. This report has provided, for the first time, an insight of the role of cattle in the transmission of the zoonotic G1 echinococosis. CONCLUSIONS: The results of the study illustrate that Sudanese breeds of cattle may play an important role in the transmission dynamics and the epidemiology of cystic echinococcosis in Sudan. This study reports the first molecular identification of E. granulosus s.s. in cattle in Central Sudan.

The incrimination of three trypanosome species in clinically affected German shepherd dogs in Sudan.

Canine trypanosomosisis (CT) is a common disease caused by tsetse- and non-tsetse-transmitted trypanosomes worldwide. The severity of the disease varies from acute, sub-acute to chronic with non-specific clinical signs. Here, we attempt in a cross-sectional study to assess the current situation of CT and the role of dogs in transmitting trypanosomes to other domesticated animals. The study was carried out during July 2016 on 50 caged German shepherd dogs in Khartoum State to investigate the prevalence of dog trypanosomosis using both serological (CATT/Trypanosoma evansi) and molecular (KIN-PCR, RoTat1.2 VSG-PCR and TviCatL-PCR) tests to detect possible trypanosome infections. CATT/T. evansi detected antibodies against T. evansi in 15 (30%) dogs, while parasite DNA was detected in 17 (34%) dogs by RoTat1.2 PCR. In contrast, a KIN-PCR detected the subgenus Trypanozoon, Trypanosoma congolense savannah, T. congolense Kenya and T. vivax in 36 (72%), 3 (6%), 1 (2%), and 2 (4%) dogs, respectively. However, a species-specific PCR for Trypanosoma vivax was detected 7 (14%) positive cases. We concluded that CT was caused by at least three species of trypanosomes, namely T. evansi, T. vivax and T. congolense. Trypanozoon other than T. evansi could not be ruled out since other tsetse-transmitted trypanosomes have also been detected and species-specific PCRs were not used. This study illustrates that dogs play an important role in the transmission dynamic and the epidemiology of the abovementioned trypanosome species.

Historical demographic profiles and genetic variation of the East African Butana and Kenana indigenous dairy zebu cattle.

Butana and Kenana breeds from Sudan are part of the East African zebu Bos indicus type of cattle. Unlike other indigenous zebu cattle in Africa, they are unique due to their reputation for high milk production and are regarded as dairy cattle, the only ones of their kind on the African continent. In this study, we sequenced the complete mitochondrial DNA (mtDNA) D-loop of 70 animals to understand the maternal genetic variation, demographic profiles and history of the two breeds in relation to the history of cattle pastoralism on the African continent. Only taurine mtDNA sequences were identified. We found very high mtDNA diversity but low level of maternal genetic structure within and between the two breeds. Bayesian coalescent-based analysis revealed different historical and demographic profiles for the two breeds, with an earlier population expansion in the Butana vis a vis the Kenana. The maternal ancestral populations of the two breeds may have diverged prior to their introduction into the African continent, with first the arrival of the ancestral Butana population. We also reveal distinct demographic history between the two breeds with the Butana showing a decline in its effective population size (Ne ) in the recent past ~590 years. Our results provide new insights on the early history of cattle pastoralism in Sudan indicative of a large ancient effective population size.

Robustness in population-structure and demographic-inference results derived from the Aedes aegypti genotyping chip and whole-genome sequencing data.

The mosquito Aedes aegypti is the primary vector of many human arboviruses such as dengue, yellow fever, chikungunya, and Zika, which affect millions of people worldwide. Population genetic studies on this mosquito have been important in understanding its invasion pathways and success as a vector of human disease. The Axiom aegypti1 SNP chip was developed from a sample of geographically diverse A. aegypti populations to facilitate genomic studies on this species. We evaluate the utility of the Axiom aegypti1 SNP chip for population genetics and compare it with a low-depth shotgun sequencing approach using mosquitoes from the native (Africa) and invasive ranges (outside Africa). These analyses indicate that results from the SNP chip are highly reproducible and have a higher sensitivity to capture alternative alleles than a low-coverage whole-genome sequencing approach. Although the SNP chip suffers from ascertainment bias, results from population structure, ancestry, demographic, and phylogenetic analyses using the SNP chip were congruent with those derived from low-coverage whole-genome sequencing, and consistent with previous reports on Africa and outside Africa populations using microsatellites. More importantly, we identified a subset of SNPs that can be reliably used to generate merged databases, opening the door to combined analyses. We conclude that the Axiom aegypti1 SNP chip is a convenient, more accurate, low-cost alternative to low-depth whole-genome sequencing for population genetic studies of A. aegypti that do not rely on full allelic frequency spectra. Whole-genome sequencing and SNP chip data can be easily merged, extending the usefulness of both approaches.

Mitochondrial DNA control-region and coding-region data highlight geographically structured diversity and post-domestication population dynamics in worldwide donkeys.

Donkeys (Equus asinus) have been used extensively in agriculture and transportations since their domestication, ca. 5000-7000 years ago, but the increased mechanization of the last century has largely spoiled their role as burden animals, particularly in developed countries. Consequently, donkey breeds and population sizes have been declining for decades, and the diversity contributed by autochthonous gene pools has been eroded. Here, we examined coding-region data extracted from 164 complete mitogenomes and 1392 donkey mitochondrial DNA (mtDNA) control-region sequences to (i) assess worldwide diversity, (ii) evaluate geographical patterns of variation, and (iii) provide a new nomenclature of mtDNA haplogroups. The topology of the Maximum Parsimony tree confirmed the two previously identified major clades, i.e. Clades 1 and 2, but also highlighted the occurrence of a deep-diverging lineage within Clade 2 that left a marginal trace in modern donkeys. Thanks to the identification of stable and highly diagnostic coding-region mutational motifs, the two lineages were renamed as haplogroup A and haplogroup B, respectively, to harmonize clade nomenclature with the standard currently adopted for other livestock species. Control-region diversity and population expansion metrics varied considerably between geographical areas but confirmed North-eastern Africa as the likely domestication center. The patterns of geographical distribution of variation analyzed through phylogenetic networks and AMOVA confirmed the co-occurrence of both haplogroups in all sampled populations, while differences at the regional level point to the joint effects of demography, past human migrations and trade following the spread of donkeys out of the domestication center. Despite the strong decline that donkey populations have undergone for decades in many areas of the world, the sizeable mtDNA variability we scored, and the possible identification of a new early radiating lineage further stress the need for an extensive and large-scale characterization of donkey nuclear genome diversity to identify hotspots of variation and aid the conservation of local breeds worldwide.

Genomic Reference Resource for African Cattle: Genome Sequences and High-Density Array Variants.

The diversity in genome resources is fundamental to designing genomic strategies for local breed improvement and utilisation. These resources also support gene discovery and enhance our understanding of the mechanisms of resilience with applications beyond local breeds. Here, we report the genome sequences of 555 cattle (208 of which comprise new data) and high-density (HD) array genotyping of 1,082 samples (537 new samples) from indigenous African cattle populations. The new sequences have an average genome coverage of ~30X, three times higher than the average (~10X) of the over 300 sequences already in the public domain. Following variant quality checks, we identified approximately 32.3 million sequence variants and 661,943 HD autosomal variants mapped to the Bos taurus reference genome (ARS-UCD1.2). The new datasets were generated as part of the Centre for Tropical Livestock Genetics and Health (CTLGH) Genomic Reference Resource for African Cattle (GRRFAC) initiative, which aspires to facilitate the generation of this livestock resource and hopes for its utilisation for complete indigenous breed characterisation and sustainable global livestock improvement.

Rapid detection and identification of Theileria equi and Babesia caballi by high-resolution melting (HRM) analysis.

The application of high-resolution melting (HRM) analysis in the differentiation between Theileria equi and Babesia caballi was evaluated using control samples from the United States Department of Agriculture and field samples collected from horses in Sudan and China. A region of the 18S rRNA gene, with four known nucleotide differences between the two parasites, was selected for primer design. HRM analysis successfully allowed the detection and differentiation of T. equi and B. caballi without the necessity of performing time-consuming and expensive post-PCR procedures such as sequencing or restriction digestion. Our results suggest that HRM could be an ideal method for rapid genotyping, which is required to determine a drug of choice or to administer an appropriate vaccine during an outbreak.

Epidemics of Crimean-Congo Hemorrhagic Fever (CCHF) in Sudan between 2010 and 2020.

Crimean-Congo hemorrhagic fever (CCHF) is a zoonotic arboviral disease that poses a great threat to global health in the Old World, and it is endemic in Europe, Asia, and Africa, including Sudan. In this retrospective study, we reviewed previous epidemiological reports about the major epidemics of CCHF throughout Sudan between 2010 and 2020. During these epidemics, the infection of humans with Crimean-Congo hemorrhagic fever virus (CCHFV), the causative agent of CCHF, was diagnosed using qRT-PCR. We have identified 88 cases of CCHF, including 13 fatalities reported during five epidemics that occurred in 2010, 2011, 2015, 2019, and 2020. The two epidemics in 2010 and 2011 were by far the largest, with 51 and 27 cases reported, respectively. The majority of cases (78%) were reported in the endemic region of Kordofan. Here, we document that the first emergence of CCHFV in the Darfur region, West Sudan, occurred in 2010. We were not able to investigate outbreak dynamics through phylogenetic analysis due to the limited diagnostic capacity and the lack of sequencing services in the country. These findings call for establishing a genomic-based integrated One Health surveillance and response system for the early preparedness, prevention, and control of CCHF in the country.

Population Genetic Analysis of Aedes aegypti Mosquitoes From Sudan Revealed Recent Independent Colonization Events by the Two Subspecies.

Increases in arbovirus outbreaks in Sudan are vectored by Aedes aegypti, raising the medical importance of this mosquito. We genotyped 12 microsatellite loci in four populations of Ae. aegypti from Sudan, two from the East and two from the West, and analyzed them together with a previously published database of 31 worldwide populations to infer population structure and investigate the demographic history of this species in Sudan. Our results revealed the presence of two genetically distinct subspecies of Ae. aegypti in Sudan. These are Ae. aegypti aegypti in Eastern Sudan and Ae. aegypti formosus in Western Sudan. Clustering analysis showed that mosquitoes from East Sudan are genetically homogeneous, while we found population substructure in West Sudan. In the global context our results indicate that Eastern Sudan populations are genetically closer to Asian and American populations, while Western Sudan populations are related to East and West African populations. Approximate Bayesian Computation Analysis supports a scenario in which Ae. aegypti entered Sudan in at least two independent occasions nearly 70-80 years ago. This study provides a baseline database that can be used to determine the likely origin of new introductions for this invasive species into Sudan. The presence of the two subspecies in the country should be consider when designing interventions, since they display different behaviors regarding epidemiologically relevant parameters, such as blood feeding preferences and ability to transmit disease.

Emergence of the invasive malaria vector Anopheles stephensi in Khartoum State, Central Sudan.

The emergence of the Asian invasive malaria vector, Anopheles stephensi, has been identified in Khartoum, the capital city of Sudan. This is the first report that confirms the geographical expansion of this urban mosquito into Central Sudan. We urgently recommend the launch of a national entomological survey to determine the distribution of this invasive disease vector and to generate essential information about its bionomics and susceptibility to available malaria control measures.

Cattle anaplasmosis and babesiosis: Major tick-borne diseases affecting the cattle industry in Khartoum State, Sudan.

Tick-borne hemoparasitic (TBH) infections are a major problem affecting livestock industries worldwide, particularly in tropical and subtropical regions. This study was carried out in response to repeated reports from local veterinarians in Khartoum State, Sudan, where TBH infections are prevalent in dairy farms. This cross-sectional study was undertaken from October 2017 to April 2018 with the objective of assessing the prevalence and potential risk factors associated with cattle anaplasmosis and babesiosis in the localities of Omdurman, Khartoum, and Khartoum North, Khartoum State. A total of 292 cattle blood samples collected from apparently healthy animals were examined for the presence of A. marginale, Babesia bigemina, and B. bovis using PCR. The overall prevalence of A. marginale and B. bigemina was found to be 40.41% and 3.42%, respectively, while B. bovis was not detected. Mixed infections with A. marginale and B. bigemina were detected in four (1.37%) cattle. The prevalence of the two pathogens was found to be significantly higher in Khartoum and Omdurman than in Khartoum North. However, no significant difference was observed for the prevalence based on sex, age, breed, and mean packed cell volume values. Our findings indicated that A. marginale is a highly prevalent parasite in Khartoum State, which may be a primary constraint to the cattle industry. Inclusion of this pathogen in the diagnostic protocols, and consequent treatment and tick control are necessary. Moreover, the role of B. bigemina infection may exacerbate the situation to some extent in this region.

Molecular Characterization of Ticks and Tick-Borne Pathogens in Cattle from Khartoum State and East Darfur State, Sudan.

Ticks transmit many pathogens with public health and veterinary importance. Despite the wide distribution of tick-borne pathogens in Sudan, the information on the tick-pathogen relationship needs to be updated, particularly using modern molecular techniques. This cross-sectional study, conducted between September and November 2019, used morphology, PCR, and sequencing to confirm the identity of adult cattle ticks (male and female; n = 536) from Khartoum State (n = 417) and East Darfur State (n = 119). Moreover, the presence of Theileria annulata, Babesia bigemina, B. bovis, Anaplasma marginale, and Ehrlichia ruminantium was detected and confirmed in each tick using species-specific PCR or nested PCR and sequencing. The most economically important tick genera, Rhipicephalus, Hyalomma, and Amblyomma, were prevalent in the study area, and 13 different tick species were identified. The most prevalent tick species were Rhipicephalusevertsi evertsi (34.3%) and Hyalomma anatolicum (57.3%) in Khartoum State, and Rhipicephalus annulatus (27%), Rhipicephalus decoloratus (25%), and Hyalomma rufipes (29%) in East Darfur State. We detected all five pathogens in both states. To the best of our knowledge, this is the first study to report the presence of E. ruminantium, its vector Amblyomma variegatum, and B. bovis in Khartoum State. Further, this is the first report on most tick and pathogen species identified in East Darfur State. Our findings indicate the migration of some tick and pathogen species beyond their distribution areas in the country, and this consideration is necessary to develop future control strategies.

BoLA-DRB3 gene haplotypes show divergence in native Sudanese cattle from taurine and indicine breeds.

Autochthonous Sudanese cattle breeds, namely Baggara for beef and Butana and Kenana for dairy, are characterized by their adaptive characteristics and high performance in hot and dry agro-ecosystems. They are thus used largely by nomadic and semi-nomadic pastoralists. We analyzed the diversity and genetic structure of the BoLA-DRB3 gene, a genetic locus linked to the immune response, for the indigenous cattle of Sudan and in the context of the global cattle repository. Blood samples (n = 225) were taken from three indigenous breeds (Baggara; n = 113, Butana; n = 60 and Kenana; n = 52) distributed across six regions of Sudan. Nucleotide sequences were genotyped using the sequence-based typing method. We describe 53 alleles, including seven novel alleles. Principal component analysis (PCA) of the protein pockets implicated in the antigen-binding function of the MHC complex revealed that pockets 4 and 9 (respectively) differentiate Kenana-Baggara and Kenana-Butana breeds from other breeds. Venn analysis of Sudanese, Southeast Asian, European and American cattle breeds with 115 alleles showed 14 were unique to Sudanese breeds. Gene frequency distributions of Baggara cattle showed an even distribution suggesting balancing selection, while the selection index (ω) revealed the presence of diversifying selection in several amino acid sites along the BoLA-DRB3 exon 2 of these native breeds. The results of several PCA were in agreement with clustering patterns observed on the neighbor joining (NJ) trees. These results provide insight into their high survival rate for different tropical diseases and their reproductive capacity in Sudan's harsh environment.

Exploring Prokaryotic and Eukaryotic Microbiomes Helps in Detecting Tick-Borne Infectious Agents in the Blood of Camels.

Dromedary camels (Camelus dromedarius) are widely distributed in Africa, the Middle East and northern India. In this study, we aimed to detect tick-borne pathogens through investigating prokaryotic and eukaryotic microorganisms in camel blood based on a metagenomic approach and then to characterize potentially pathogenic organisms using traditional molecular techniques. We showed that the bacteria circulating in the blood of camels is dominated by Proteobacteria, Bacteroidetes, Firmicutes and Actinobacteria. At the genus level, Sediminibacterium, Hydrotalea, Bradyrhizobium and Anaplasma were the most abundant taxa. Eukaryotic profile was dominated by Fungi, Charophyta and Apicomplexa. At the genus level, Theileria was detected in 10 out of 18 samples, while Sarcocystis, Hoplorhynchus and Stylocephalus were detected in one sample each. Our metagenomic approach was successful in the detection of several pathogens or potential pathogens including Anaplasma sp., Theileria ovis, Th. separata, Th. annulate, Th. mutans-like and uncharacterized Theileria sp. For further characterization, we provided the partial sequences of citrate synthase (gltA) and heat-shock protein (groEL) genes of Candidatus Anaplasma camelii. We also detected Trypanosoma evansi type A using polymerase chain reaction (PCR) targeting the internal transcribed spacer 1 (ITS1) region. This combined metagenomic and traditional approach will contribute to a better understanding of the epidemiology of pathogens including tick-borne bacteria and protozoa in animals.