Mycobacterium bovis (M. bovis) is the main causative agent for bovine tuberculosis (BTB) and can also be the cause of zoonotic tuberculosis in humans. In view of its zoonotic nature, slaughterhouse surveillance, potentially resulting in total or partial condemnation of the carcasses and organs, is conducted routinely. Spoligotyping, VNTR profiling, and whole genome sequencing (WGS) of M. bovis isolated from tissues with tuberculosis-like lesions collected from 14 cattle at Eritrea’s largest slaughterhouse in the capital Asmara, were conducted. The 14 M. bovis isolates were classified into three different spoligotype patterns (SB0120, SB0134 and SB0948) and six VNTR profiles. WGS results matched those of the conventional genotyping methods and further discriminated the six VNTR profiles into 14 strains. Furthermore, phylogenetic analysis of the M. bovis isolates suggests two independent introductions of BTB into Eritrea possibly evolving from a common ancestral strain in Europe. This molecular study revealed the most important strains of M. bovis in Eritrea and their (dis) similarities with the strains generally present in East Africa and Europe, as well as potential routes of introduction of M. bovis. Though the sample size is small, the current study provides important information as well as platform for future in-depth molecular studies on isolates from both the dairy and the traditional livestock sectors in Eritrea and the region. This study provides information onthe origin of some of the M. bovis strains in Eritrea, its genetic diversity, evolution and patterns of spread between dairy herds. Such information is essential in the development and implementation of future BTB control strategy for Eritrea.
Mycobacterium bovis (M. bovis) is the causative agent of bovine tuberculosis (BTB), achronic, infectious and contagious disease that also affects other domestic animals as well as humans [1,2]. Although BTB is prevalent in dairy cattle in Eritrea as shown by Omer et al. (2001) [3] and Ghebremariam et al. (2016) [4] by skin-test based survey, detection and isolation of the causative agent has never been done. Routine meat inspection at municipal slaughterhouses is performed for identifying tuberculosis-like lesions (TBL) that usually results in either total or partial condemnation of carcasses depending on the level of TBL dissemination, however, confirmatory testing or trace back epidemiological investigations are not conducted in Eritrea. Genotyping is a vital tool for trace back in epidemiological investigations, and according to Biek et al. (2012) [5] results from WGS alone can provide insight into TB epidemiology even in the absence of detailed contact data. Despite the usefulness of genotyping, it is rarely used in developing countries, i. e. , in Africa, Asia, and South America [6–9]. The routine use of such tool in these countries could be instrumental in complementing BTB control strategies. Spoligotyping and variable number of tandem repeat (VNTR) profiling have been used extensively in many countries to document the molecular epidemiology of Mycobacterium tuberculosis complex (MTBC) species [7,10–14]. For this reason, the digital MTBC molecular genotypes are predominantly stored in these two forms globally[15–18]. The recent technological advancements in molecular genetics imply that we can now more than ever understand the molecular epidemiology of MTBC at amore granular level. In the last few years, whole genome sequencing (WGS) for typing of pathogens has been explored and yielded important additional information on strain diversity in comparison to the classical DNA typing methods. Analysis of data from WGS also allows detection of minute differences in genetic diversity and this has contributed retrospectively to outbreak investigations [19–23]. Significantly, WGS allows for better genomic coverage withsingle nucleotide polymorphisms (SNP) profiling than the two classical typing methods [24,25]. WGS has also led to a significant growth in quantitative methodology that allows for a robust estimation of phylogenetic and temporal relationships between samples[26]. All these aspects are essential in enhancing our understanding of local and distant, recent and historical dynamics of BTB [5,24]. Although several reports predict that the use of WGS for genotyping will eclipse the classical MTBC typing tools [27], this will likely take longer to occurin Africa. It is therefore important to compare their utility in resource limited settings. Although such tools have never been used in Eritrea, their use would greatly enhance our understanding of: a) the genetic diversity of M. bovis, b) its evolution and c) the patterns of spread (spatial and temporal) between dairy herds, in the country and region. Such data (information) would be critical for safeguarding and further development of the dairy industry of Eritrea. In the present study, the classical MTBC typing tools (Spoligotyping and MIRU-VNTR) as well as WGS were used to gain insight into the spatial and temporal dynamics, genetic diversity and evolution of M. bovis strains circulating in Eritrean dairy cattle. Furthermore, to infer local and international historical phylogenetic relationships.
Pooled tissue samples (lungs and pleura, mediastinal, bronchial, deep inguinal and lung lymph nodes), were collected from 15 animals that showed TBL in gross pathology, at the Asmara municipal slaughterhouse from March 2014 to May 2015. These 15 animals were all those with TBL during the study period. The animals were slaughtered for meat purpose and processed as part of the normal work of the abattoir. Approximately 5–10 grams of pooled tissues from each sampled animal were collected in sterile specimen containers, and immediately transported on icepacks to the National Animal and Plant Health Laboratory (NAPHL), Asmara, and stored at -20°C until processing for culture. Data collected from individual animals (Table 1) included: source of the animal slaughtered, date of slaughter, species, breed, sex, age, pregnancy status (pregnant/non-pregnant), ante mortem clinical signs, post mortem lesions, and type of the tissue samples collected. In addition, retrospective meat inspection data for the period 2010 to 2015 were retrieved from the logbook of the slaughterhouse. Samples were processed for M. bovis culture as follows: approximately 5 g of each pooled tissue sample with TBL per animal was cut into small pieces and covered with 100 ml of sterile distilled water in a biohazard cabinet (Esco Class II BSC; Labotec, SA). The samples were homogenized using an Ultra-Turrax homogenizer at 17500 rpm (Separation Scientific, SA). Seven millilitres of the homogenate was poured into each of two separate 15 ml falcon tubes, and the remaining homogenate was poured into individual 50 ml centrifuge tubes and stored at -20°C as reference samples. The samples were decontaminated with 7 ml of 2% HCL (final concentration of 1%) and 7 ml of 4% NaOH (final concentration of 2%), respectively, and incubated at room temperature (18–25°C) for 10 min. After subsequent centrifugation (HeraeusLabofuge 400) of the samples at 3500 rpm for 10 min. , supernatants were poured off and 7 ml of sterile distilled water was added. After washing, the centrifugation step was repeated and most of the supernatant was poured off. The pellets were re-suspended in a volume of approximately 1ml using a sterile inoculation loop. Two loops of each of the pellets were spread evenly onto two Löwenstein-Jensen (L-J) media slants supplemented with pyruvate (National Health Laboratory Service, SA) and onto one L-J medium slant supplemented with glycerol (BD Diagnostics), and incubated at 37°C for up to ten weeks. The slants were monitored weekly for mycobacterial growth. Ziehl-Neelsen staining was conducted andlysate (DNA) of acid fast bacteria was subjected to polymerase chain reaction (PCR) testing to identify bacteria as MTBCas previously described [28,29]. Subsequently, deletion analysis was performed on the isolates using PCR primers targeting the RD4 (region of difference-4) as previously described for M. bovis identification [30]. The three features used to distinguish M. bovis clonal complexes were: a) they are a derivative of most recent clonal ancestors (MRCA) spoligotype b) region of difference deletion and c) geographic restriction (Example: African 1 is localized in West Africa) To obtain the whole genome sequences, DNA of the 14 Eritrea M. bovis isolates was extracted (dx. doi. org/10. 17504/protocols. io. nsgdebw) and sequenced on a MiSeq instrument (Illumina, San Diego, CA) using 2x250 paired-end chemistry and the Nextera XT library preparation kit (Illumina, San Diego, CA). FASTQ files from the instrument were put through the National Veterinary Services Laboratories (NVSL) in-house pipeline (see https: //github. com/USDA-VS). Briefly, reads were aligned to the reference genome AF2122/97, NCBI accession number NC_0002945, using BWA and Samtools[32,33]. A depth of coverage of 80X was targeted. BAM files were processed using Genome Analysis Toolkit (GATK) ’s best practice workflow. SNPs were called using GATK’s HaplotypeCaller outputting them to variant call files (VCF) [34–36]. Results were filtered using a minimum QUAL score of 150 and AC = 2. From the VCFs, SNPs gathered were outputted to three formats: an aligned FASTA file; tab-delimited files sorted by position location and by SNP groups; and a maximum likelihood phylogenetic tree created with RAxML[37]. The tree was built using a GTR-CAT model with input taken as an alignment file containing only informative and validated SNPs. SNPs were visually validated using Integrative Genomics Viewer (IGV) [38]. Because WGS isolates from this region of the globe are not readily available, databases from three laboratories (United States Department of Agriculture, Centre de Recercaen Sanitat Animal (CReSA) —Institute de Recerca i Technologia Agroalimentáries (IRTA), Spain, and Tuberculosis Research Laboratory, Department of Veterinary Public Health and Preventive Medicine, University of Ibadan, Nigeria) that are actively sequencing M. bovis isolates were queried and field isolates that were within 150 SNPs of the Eritrea isolates were included in our analysis. Also included for perspective were widely available reference strains, AN5, Ravenel, 95–1315, AF2122/97, BCG, and BZ-31150. “FASTQ files from the isolates sequenced were uploaded into NCBI short read archive. Accession numbers Bioproject and sample numbers are listed in supplemental S1 Table.
During the period 2010 to 2015,78,820 cattle were slaughtered and 38 carcasses, originating from Maekel and Debub regions, were totally condemned due to generalized TBLs showing caseous necrosis identified in gross pathology in the lungs, livers, pleura (chest), peritoneum, and lymph nodes. Besides, fore quarters of three animals, plucks, shoulders, chests and heads of six cattle were partially condemned due to the presence of TBL (Table 2). All, except one (local breed), of the condemned carcasses were of the exotic HF breed or their crosses. Out of the 15 animals sampled from March 2014 to May 2015, nine originated from Maekel and one from Debub, while the origin of the other five slaughtered animals was unknown due to lack of records. Detailed gross pathology information on the tissues collected is presented in Table 1. During this period 26,603 cattle were slaughtered and nine out of the 15 carcasses sampled, were totally condemned due to generalized TBL. In addition, the entire plucks and shoulders of three animals were partially condemned (Table 2), and from three other animals, tissues with TBL were collected and the carcasses passed for consumption. Out of the 15 pooled tissue samples cultured on L-J media slants supplemented with pyruvate, 14 yielded smooth dysgonic growth, suggestive of M. bovis presence. All the 14 isolates were identified as MTBC. Subsequent examination by M. bovis specific PCR targeting the RD4, yielded banding patterns typical of M. bovis with a 268 bp product indicating RD4 deletion.
The dominant spoligotype identified in our study was SB0120, named BCG-like by Haddad et al. [13] and considered as parental strain for the M. bovis vaccine strain. It accounted for 64% of the isolates, whereas the other two spoligotypes SB0134 and SB0948 did so for 29% and 7%, respectively. The first two strains are widely distributed in a number of African countries, namely; Ethiopia, Algeria, Zambia, South Africa [6,10,16,40–43] as well as in Italy, Spain, other European countries and Mexico [13,44–49]. In addition to cattle, SB0120 affects wildlife and humans in Africa and Europe [50–53]. The third spoligotype (SB0948) has been reported in France, Italy, and in Zambia [13,41,44]. The relatively high frequency of the spoligotype SB0120 found in the present study may indicate its predominance in Eritrean dairy cattle, though difficult to conclude with such small sample size. The second predominant spoligotype (SB0134) appears to have evolved from spoligotype SB0120 by the loss of spacers 4 and 5 in addition to spacers 3,9, 16, and 39–43 that classify spoligotype SB0120. This finding might not be surprising, in view of the past trade relations between Eritrea and Ethiopia, as both SB0120 and SB0134 spoligotypes are also present in Ethiopia. Besides, these two countries share open borders that consequently allowed the uncontrolled movement of animals as obtained in most African countries. Therefore, it might be plausible to speculate that these strains of M. bovis are shared between Eritrea and Ethiopia. On the other hand, it might also be plausible to suggest Italy as a possible source of these strains, on the following grounds: a) long historical ties (1900 to 1970s) between Eritrea and Italy existed, b) Italian settlers initiated the establishment of dairying in Eritrea in the 19th century by importing exotic breeds (Holstein–Friesian), c) the M. bovis spoligotypes detected in our study are also wide spread in Italy. Although the reason for the apparent predominance of the two spoligotypes (SB0120 and SB0134) needs further study, as this may indicate an epidemiological link between different dairy farms/regions in Eritrea, as buying and selling of cows between dairy farms is common in the country[54] without following strict sanitary rules. Since not all the slaughtered cattle with TBL had records of their farm of origin, it may also be possible to suspect that some of the slaughteredanimals might have originated from the same farm. It is noteworthy, however, that based on the WGS data there appears to be at least two introductions of M. bovis into Eritrean dairy cattle, an SB0120 strain and SB0134 strain. The SB1517 (Ethiopian strain; Fig 2) is an offspring of SB0134 suggesting that the common ancestor of the cluster was SB0134. Spoligotype SB0948 was found in only one animal. It is a descendant of spoligotype SB0120 as it differs by the absence of spacer 1 only, and deviates only by the Mtub21 locus in its VNTR profile from the other members of the SB0120 group (Fig 1). Further, the WGS data confirmed that SB0948 is a recent descendant of a sub-cluster of SB0120 isolates. Though unclear what its relevance is in neighboring Ethiopia, this spoligotype was reported in several countries in Africa and Europe [13,41,44,48,55]. The African and global comparisons of spoligotype profiles (Fig 3and S1 Fig) demonstrated the regional and global distribution of the spoligotype and VNTR profiles and their similarities with the Eritrean ones. These similarities could be attributed to the following two plausible reasons: a) inter-regional and global livestock trade, b) colonial livestock and livestock product trade within their then colonies and outside. Variable number tandem repeat (VNTR) profiles are considered appropriate to complement spoligotyping due to their ability to discriminate between M. bovis strains as defined by spoligotyping[15,55,56]. The three spoligotypes were clustered into six VNTR profiles (Fig 1). The diversity seen in the VNTR profiles may suggest that M. bovis has been circulating in the dairy herds of the country for quite a long time with only minor mutations as the BCG-like spoligotype (SB0120) is the predominant one. Four of the VNTR profiles (ER-2, -3, -4 and 6) may have derived from the predominant VNTR profile ER-1, that corresponds to spoligotype SB0120. According to Smith et al. [49], strains bearing the same spoligotype pattern are assumed to be a set of individuals derived relatively recently by clonal replication from a single ancestral cell. On the basis of the VNTR profile, both strains, SB0948 and SB0134, are clustered within the SB0120 group with a loss of only one locus (Mtub12) in the former and two loci (ETR-B and ETR-E), in the latter strain, respectively. One of the VNTR profiles within the SB0134 strain exhibited two different VNTR alleles (3 and 4 tandem repeats) for locus ETR-E (Fig 1), suggesting either a mixed infection with two distinct strains or a microevolution in this strain. The VNTR profiles found in our study showed clonal variants differing at their loci as compared to what was reported in other parts of Africa (i. e. , Zambia) (Fig 1 & Fig 3), though they were all M. bovis strains belonging to SB0120 spoligotype. This clonal difference (Fig 1 & Fig 3) seen in our study may have been attributed to the absence of active livestock (dairy cattle) trade between Eritrea and other parts of Africa (Zambia) or due to the different geographical locations and livestock management systems between the countries, that might have dictated the microevolutions (mutations) differently. The possible reason for having the same spoligotype (SB0120) in Eritrea and other African countries (S1 Fig), might be that the source of the cattle for Eritrea and the other countries was Europe, as Europe is the source for the high yielding dairy cows, like the Holstein Friesians, that are imported by most African countries with the aim of improving milk production in their countries in order to realize food security. The investigation of the 14 M. bovis isolates for clonal complex differentiation revealed that they belonged to none of the complexes identified so far i. e. , African 1 (RDAf1), African 2 (RDAf2), European 1 (RDEu1) and European 2 (RDEu2) [15–18]. The absence of members that belong to clonal complex African 1& 2 in our samples could suggest limited introduction of such strains from the neighboring Eastern and Western African livestock movement routes. It is noteworthy, that these two strains (SB0120 and SB0134) are present in Ethiopia [16], although most of the other strains in this country belong to Af2. In the current study, little strain diversity is recorded (Fig 1) as compared to studies conducted in other countries with similar agricultural setting like Eritrea[6,42]. The WGS results matched the conventional laboratory methods with better resolution. These data support two separate introductions of M. bovis into Eritrea, with subsequent localized spread. The common ancestor of these two groups is shared widely with isolates in the USA and Spain, with greater diversity found in Spain suggesting an introduction from Europe. The presence of a common ancestor in these distantly located countries may be due to the international livestock trade between these countries, geographical proximity and similar livestock production systems. Example: the origin of the high yielding dairy breed (Holstein Friesian) is Europe. As indicated in the spoligotyping section above, the spoligotype SB0120, predominant in our study, is also ubiquitous in Europe, especially in France[13], Italy [44], Portugal[45], and Spain[48], most likely as a result of geographical closeness and trade relations between these countries. Therefore, our finding may not be a surprise, given the historical establishment of ‘intensive’ dairy farming by the Italian settlers in Eritrea through the importation of high yielding dairy breeds (Holstein Friesians) to meet the high demand for milk and dairy products. The fact that the Eritrean strains are between (close to) Spain samples (Fig 2) may suggest two introductions or may be just one introduction; i. e. , from Europe (Italy). Since we do not have information that shows historical, political or trade ties between Spain and Eritrea, we can speculate that either the strains are circulating in Italy and Spain. Or that, the Italian settlers during the establishment of dairy farms may have imported the cattle from Spain or other European countries where the same strains of M. bovis might have been circulating. A classical analogy for this speculation may be rinderpest that was brought to Sub-Saharan Africa by Italian forces in 1889, with infected cattle they had imported from India, Aden, South Russia to feed their army that had then occupied Massawa (Eritrea) [57]. However, although phylogenetic comparison with Italian M. bovis isolates could not be done in our study, we cannot refute the possibility that these strains originate from Italy or via the above indicated routes from other countries. The second probable route of introduction for one of the groups of the Eritrean strains, but not for the other, may be Ethiopia considering the long and close historical relationship and uncontrolled livestock movement between these two countries. The Ethiopian and Eritrean samples have accumulated 8–16 additional SNPs since sharing a common ancestor suggesting a recent common source and regional spread. But the four Eritrean samples (strains) are within 5–6 SNPs from sharing a common ancestor suggesting these isolates have established and spread within Eritrea, though it might be premature to reach into conclusion with such small sample size. Eritrea, on the other hand, might have introduced this strain to Ethiopia. This is plausible because both intensive dairy farming, established 100 years ago in Eritrea and the first report of BTB (Pirani, 1929), cited by Omer et. al. [3], occurred earlier than in Ethiopia where ‘intensive’ dairy farming started in the 1950s (1947) by importing Friesians and Brown Swiss [58]. This was followed by the detection of acid fast bacilli in a cow’s milk, in one study, and detection of what was called ‘Mycobacteria tuberculosis bovine type’ seemingly, M. bovis from 18 cattle, in another study, in Eritrea, by Sfroza in 1944[3]. The samples collected in this study are not considered representative of all strains possibly circulating in Eritrea. However, Asmara slaughterhouse, as the country’s biggest facility mostly slaughters exotic cattle breeds from various regions in Eritrea in which previously a high BTB prevalence was reported [4]. Therefore, the panel of samples still provides a valuable insight in the genetic strain composition from mostly dairy producing regions in Eritrea and a valuable basis for future investigations. The current study characterized strains of M. bovis in Eritrea and revealed their (dis) similarities with the strains generally present in Africa and Europe, as well as potential routes of introduction of M. bovis. Though the sample size is small, our study provides important information as well as availability of technology for future in-depth molecular studies including more samples from dairy cattle as well as cattle and goats from the traditional livestock sector. This study provides information on the origin of the M. bovis strain in Eritrea, its genetic diversity, evolution and patterns of spread (spatial and temporal) between dairy herds. The information obtained will be instrumental in making informed decisions in future BTB control strategy for Eritrea. Our study has some limitations. The samples were collected from one slaughterhouse and were few due to the absence of tissues with TBL during the study period. The low prevalence of BTB in the traditional livestock raising system [59] where majority of slaughtered animals come from, has limited the possibility of detecting more M. bovis strains from different geographical regions of Eritrea. Genetic profiling of M. bovis strains is a highly useful approach which can aid in the study and control of the temporal and geographical disease spread in the country and the African continent where BTB is largely uncontrolled. We recommend future studies in Eritrea to include genetic profiling of Italian isolates so as to support or negate our hypothesis with certainty than just live with speculation that the origin of the Eritrean M. bovis strains was Italy. In future studies in Eritrea, inclusion of more regional slaughterhouses including animal traceability will enable us gain greater insight into the epidemiology of BTB in the country which will allow the M. bovis genotype to be linked to the population from which it was obtained. We also recommend that simultaneous detection and strain differentiation of M. bovis isolates should become a reality in the routine of human tuberculosis reference laboratories, as well as in the routine meat inspection at municipal slaughterhouses. Therefore, using the One Health paradigm (i. e. interdependence between the medical and veterinary fields), greater integration between agriculture and health sectors could be an important strategy to control M. bovis in several places in the world where the agent is disseminated between animals and humans.
The livestock sector plays a major role in poverty and hunger reduction in the vast majority of Africa, as a source of food, cash income, manure, draught power, transportation, savings, insurance and social status. However, for livestock to play this vital role, the impact of diseases of economic and zoonotic importance need to be reduced. Bovine tuberculosis, mainly caused by Mycobacterium bovis, is such an infectious disease. Slaughterhouse (gross pathology) surveillance, followed by bacterial culture and genotyping, are options to identify the disease-causing agents, their distribution, and enabling trace back of the sources of infections, in order to prevent their re-introduction and spread. Unfortunately, genotyping is by far not generally introduced in the continent. In the present study, tissues with tuberculosis-like lesions were collected from the Asmara municipal slaughterhouse, the largest slaughterhouse in Eritrea, and bacterial culture, classical Mycobacterium tuberculosis complex typing (Spoligotyping and VNTR profiling), as well as whole genome sequencing (WGS) were used to gain insight into the spatial and temporal distribution, genetic diversity and evolution of M. bovis strains circulating in Eritrean dairy cattle. The results revealed (dis) similarities of the Eritrean M. bovis strains with the strains generally present in Africa and Europe, potential routes of introduction to Eritrea and genetic diversity of the M. bovis strains. Future in-depth molecular studies including more samples from dairy cattle as well as cattle and goats from the traditional livestock sector are recommended.
livestock
evolutionary biology
geographical locations
eritrea
molecular genetics
molecular biology techniques
genotyping
bacteria
africa
research and analysis methods
artificial gene amplification and extension
actinobacteria
molecular biology
evolutionary genetics
agriculture
people and places
polymerase chain reaction
genetics
biology and life sciences
europe
organisms
mycobacterium bovis
Genetic profiling of Mycobacterium bovis strains from slaughtered cattle in Eritrea