Rajesh Man Rajbhandari1,2,3,4
Rajindra Napit4
Prajwol Manandhar4
Roji Raut4Anupama Gurung4
Ajit Poudel4Nisha Shrestha4Amir Sadaula5
Dibesh Karmacharya4,6
Christian Gortázar3
Paulo Célio Alves1,2,7
José de la Fuente3,8
João Queirós1,2,7*- 1Departamento de Biologia, Faculdade de Ciencias da Universidade do Porto, Porto, Portugal
- 2CIBIO, Centro de Investigação em Biodiversidade e Recursos Genéticos, InBIO Laboratório Associado—BIOPOLIS Program in Genomics, Biodiversity and Land Planning, Campus de Vairão, Universidade do Porto, Vairão, Portugal
- 3SaBio, Instituto de Investigación en Recursos Cinegéticos, IREC (UCLM and CSIC), Ciudad Real, Spain
- 4Center for Molecular Dynamics Nepal, Kathmandu, Nepal
- 5National Trust for Nature Conservation, Biodiversity Conservation Center, Chitwan, Nepal
- 6School of Biological Sciences, Faculty of Science, The University of Queensland, Brisbane, QLD, Australia
- 7EBM, Estação Biológica de Mértola, Praça Luís de Camões, Mértola, Portugal
- 8Department of Veterinary Pathobiology, Center for Veterinary Health Sciences, Oklahoma State University, Stillwater, OK, United States
Introduction: Tuberculosis is an infectious disease caused by a group of acid-fast bacilli known as Mycobacterium tuberculosis complex (MTC), which has a major impact on humans. Transmission of MTC across the human-animal interface has been demonstrated by several studies. However, the reverse zoonotic transmission from humans to animals (zooanthroponosis) has often been neglected.
Methods: In this study, we used Nanopore MinION and Illumina MiSeq approaches to sequence the whole genome of M. tuberculosis strains isolated from two deceased Asian elephants (Elephas maximus) and one human in Chitwan, Nepal. The evolutionary relationships and drug resistance capacity of these strains were assessed using the whole genome data generated by the stand-alone tool Tb-Profiler. Phylogenomic trees were also constructed using a non-synonymous SNP alignment of 2,596 bp, including 94 whole genome sequences representative of the previously described M. tuberculosis lineages from elephants worldwide (lineages 1 and 4) and from humans in Nepal (lineages 1, 2 and 3).
Results and Discussion: The new genomes achieved an average coverage of 99.6%, with an average depth of 55.67x. These M. tuberculosis strains belong to lineage 1 (elephant DG), lineage 2 (elephant PK) and lineage 4 (human), and none of them were found to have drug-resistant variants. The elephant-derived isolates were evolutionarily closely related to human-derived isolates previously described in Nepal, both in lineages 1 and 2, providing additional support for zooanthroponosis or bidirectional transmission between humans and elephants. The human-derived isolate clustered together with other published human isolates from Argentina, Russia and the United Kingdom in the lineage 4 clade. This complex multi-pathogen, multi-host system is challenging and highlights the need for a One Health approach to tuberculosis prevention and control at human-animal interface, particularly in regions where human tuberculosis is highly endemic.
1. Introduction
Tuberculosis (TB) is a significant global burden and is widely reported to be a major public health and economic problem, costing the world $617 billion between 2000 and 2015 and projected to cost $1 trillion between 2015 and 2030 (1). It is the second leading cause of death after COVID-19, with an estimated 10 million cases worldwide and 1.5 million deaths in 2020 (2). Drug-resistant TB is also a major threat to global disease control (3). It is a major contributor to global antimicrobial resistance (4). In 2021, 450,000 cases of rifampicin-resistant and multidrug-resistant (MDR/MDR-TB) TB were forecasted worldwide, of which 191,000 died (5). In Nepal, the number of cases of drug-resistant TB is estimated to be around 1,500 per year (6).
Human and animal tuberculosis is caused by a group of closely related acid-fast bacilli known as the Mycobacterium tuberculosis complex (MTC) (7, 8). The MTC includes several species of mycobacteria, namely M. tuberculosis, Mycobacterium bovis (Bacillus Calmette–Guerin), Mycobacterium africanum, Mycobacterium caprae, Mycobacterium microti, Mycobacterium canettii and Mycobacterium pinnipedii (9). M. tuberculosis is the human-adapted variant and is the major cause of human TB (9). However, most publications on emerging infectious agents such as MTC often focus on their zoonotic origin, but comparatively fewer reports are published on possible zooanthroponosis and the human origin of infectious diseases affecting domestic or wild animals (10). Zooanthroponosis, is an important ongoing debate regarding the transmission of pathogens from humans to animals. Evidence from a global survey showed that humans are capable of transmitting at least 21 bacterial, 12 viral and 7 fungal pathogens to animals (10). Although there is evidence that pathogens can be transmitted from humans to animals, most reported cases involve captive or domestic animals (11). Influenza is one of the best known examples of bidirectional transmission between humans and domestic pigs (12). The emergence and re-emergence of these pathogens has important implications for human and animal health (13). This risk appears to be even greater when bidirectional transmission between humans and animals is frequent, as in the case of TB in elephants (14).
Tuberculosis in elephants has been reported worldwide, although most reported cases involve captive elephants with M. tuberculosis, the ethological agent of human TB (15–19). In countries such as Nepal, where there is a high prevalence of active M. tuberculosis infection in humans (117,000 people living with TB; 20) and a large elephant population (more than 200 captive and 200–250 wild elephants; 21), this issue is particularly relevant. Elephant-human interactions are particularly high in regions where captive elephants are used for ecotourism and patrolling protected areas, such as Chitwan National Park. This poses a serious risk of zooanthroponosis and disease to both captive and wild elephant populations, as well as a risk to humans given the possibility of bidirectional transmission (15, 22). A recent study in four national parks in Nepal showed a high seroprevalence (21.56%) of TB in captive elephants, with most cases detected in Chitwan National Park (23). In addition, some evidence of zooanthroponosis has been reported from previous studies in Nepal that also isolated M. tuberculosis from diseased elephants using multi-locus variable number of tandem repeats (MLVA) (15, 24), large sequence polymorphism (24) or whole genome data (25). However, the latter study lacks a phylogenetic analysis of the two M. tuberculosis genomes (26) and their integration with previously published data (25).
In this study, a whole genome sequencing approach was used to genotype M. tuberculosis isolates from a human and two deceased elephants in Chitwan, Nepal. These new data, combined with genomic data from elephants and human isolates available in the literature, allowed a comprehensive evolutionary study of M. tuberculosis in Nepal. These results contribute to the understanding of the possible bidirectional transmission of this pathogen at human-elephant interface.
2. Methods and materials
2.1. Study area
Chitwan district is located in the southwestern part of Bagmati Province in Nepal. It covers an area of 2,238.39 km2 and had a human population of 719,859 in 2022. Nepal’s first national park, Chitwan National Park (CNP), is located in this district (Figure 1) and was established in 1973 (27). The park covers an area of 932 km2 and is located in the subtropical Terai lowlands of south-central Nepal. It lies in a river valley basin or dun, along the floodplains of the Rapti, Reu and Narayani rivers. The Chitwan valley consists of tropical and subtropical forests. Sal forest covers 70% of the park. Sal leaves are used locally for plates in festivals and religious offerings, and are also known for their Ayurvedic use in the treatment of wounds, coughs and other ailments (28). Grasslands cover 20 per cent of the park. There are more than 50 different species of grasses, including elephant grass (Saccharum spp), known for its immense height. The climate is mainly dominated by the summer monsoon, and the valley experiences three distinct seasons each year: winter, summer and monsoon.
Figure 1. Geographical map of Nepal showing the location of Chitwan National Park (light green) and the origin of the human and elephant samples used in this study (icons). Forest communities around Chitwan National Park are shown in dark green.
Chitwan National Park is regarded as one of the best national parks for wildlife viewing in Asia (29). It is home to 68 species of mammals, more than 576 species of birds, 49 species of reptiles and amphibians, 120 species of fish and several species of invertebrates that contribute significantly to ecosystem processes in the region. The reserve is home to many of Nepal’s charismatic wildlife species, including the one-horned rhinoceros (Rhinoceros unicornis), Bengal tiger (Panthera tigris), leopard (Panthera pardus) and Asian elephant (Elephas maximus), among others (29, 30). Of the various recreational activities available in the park, elephant safaris are the most popular, especially in and around the buffer zones. Elephants are also used by the government for transport and patrolling the park. Buffer zones in the CNP were established in 1996 with the aim of involving people in the management of park resources for the conservation of biodiversity and the livelihoods of buffer zone communities (31). The buffer zone covers an area of 750 km2 and is home to approximately 45,516 houses and 2,60,352 people (32). Buffer zones are primarily designed to create human-wildlife coexistence by providing an ecological and socio-economic buffer for communities, and in these areas, where wildlife and humans share the same landscape in close proximity, their interaction is very obvious (33).
2.2. Sample collection
Samples were collected from two deceased adult (>30 years) elephants in Sauraha (Elephant PK and Elephant DG), Chitwan, owned by the government, and used for tourism, grass collection and patrolling. Sauraha is a village at the eastern entrance to Chitwan National Park, near the East Rapti River. These animals had died in captivity and at necropsy showed granulomatous lesions in the lungs compatible with tuberculosis. Tissue samples from lung biopsies (with lesions) were taken from both elephants. The human sputum sample (No. 52G111) collected in Chitwan as part of a previous study (20) was also genotyped, although it was not closely related to the diseased elephants. It was included to understand the evolutionary relationship with the M. tuberculosis strains isolated from elephants, as it was collected from the same region.
2.3. Molecular confirmation of Mycobacterium tuberculosis infection
The elephant lung tissue samples (PK and DG) and human sputum samples (52G111) were sent to GENETUP (German Nepal Tuberculosis Project microbiology lab) for culture confirmation of M. tuberculosis as described in (34). The heat-killed cultures were then resuspended in 100 μL of PBS for DNA extraction. The DNA extraction was performed using a DNA-SorB kit according to the manufacturer’s instructions (Sacace Biotechnologies, Italy). Briefly, 300 μL of the lysis solution was mixed with 100 μL of the heat lysed M. tuberculosis PBS solution and incubated at 60°C for 5 min. The samples were centrifuged at 15000 g for 5 min to separate the cell debris and DNA in the supernatant. The supernatant was then mixed with 20 μL of the sorbent to bind the DNA to the sorbent. The sorbent was washed twice, dried at 65°C for 5 min and eluted with 25 μL of DNA eluent. Elution was further enhanced by incubating the tubes with DNA eluent at 65°C for 5 min and periodically vortexing. The extracted DNA was then quantified using Qubit (Thermo Fisher Scientific, USA) and M. tuberculosis infection was diagnosed using the specific real-time PCR (35).
2.4. Whole-genome sequencing of elephant and human Mycobacterium tuberculosis isolates
Two complementary whole genome sequencing approaches were used in this study: Illumina short reads and MinION long reads. In the case of Illumina short reads, genomic DNA from the heat-killed cultures of M. tuberculosis was purified using AMPure beads (Beckman Coulter, United States) (1:1), quantified using Qubit (Thermo Fisher Scientific, United States) and normalised to 0.2 ng/μL. The normalised samples were tagged and indexed using the Nextera Indexing Kit (Illumina, United States). The indexed library was then purified using (0.8X) AMPure beads and quantified using HS kit (Thermo Fisher Scientific, United States). Each library was analysed using a bioanalyzer and normalised to a final concentration of 2 nM. The three normalised libraries were then pooled in equal volumes of 5 μL to give a final pooled library of 2 nM. The pooled library was spiked with 5% of 2 nM PhiX and denatured with NaOH. The denatured library was diluted to 20 pM, then 10 pM using hybridisation buffer, and the final 10 pM library was loaded onto the Illumina MiSeq platform using the MiSeq Reagent Kit V2 300-cycles (Illumina, United States), with an estimated coverage of approximately100× per sample.
For long-read sequencing in Nanopore MinION, normalised genomic DNA (0.2 ng/μL) from the heat-killed cultures of M. tuberculosis isolated from the two elephants was also used for long-read sequencing in Nanopore MinION, using the Nextera XT Library Preparation Kit (Illumina, USA). The 20 μL of tagged genomic DNA was end-repaired using Ultra II End Prep Buffer and Ultra II End Prep Enzyme (NEB, United Kingdom) according to the manufacturer’s instructions. The end-repaired DNA was then subjected to Nanopore DNA preparation using Ligation sequencing kit (Oxford Nanopore, United Kingdom) with PCR barcoding expansion kit (Oxford Nanopore, United Kingdom). The final library was purified using 0.4X AMPure beads and Long Fragment Buffer (LFB) to obtain the eluate in 15 μL elution buffer. The adapter-ligated library was quantified using the Qubit Fluorometer and approximately 15 fm of the total library was finally loaded onto the MinION for sequencing using 15 μL of sequencing buffer and 10 μL loading beads.
2.5. Lineage and multi drug resistance status determination
Tb-Profiler version 4.1.0, a stand-alone tool (36), was used for lineage identification and drug resistance assessment of M. tuberculosis isolates from whole genome sequences. The tool accepts raw data from both MinION Nanopore and Illumina MiSeq sequences as well as BAM and FASTA files. It is a command-line tool that can use whole-genome data to predict lineages and identify 21 types of drug resistance (based on small variants and large deletions associated with drug resistance).
2.6. Processing and analysis of Illumina and Nanopore sequence data
Quality control and filtering of Illumina reads was performed using Fastp v0.20 (37), while NanoFilt v2.6.0 (38) was used to filter Nanopore reads. After quality filtering, the Illumina and Nanopore reads were mapped separately to the M. tuberculosis reference genome (GenBank accession NC_000962.3) using BWA v0.7.17 (39). SAMtools v1.15 (40) was used to convert the SAM files of the mapped reads from both Illumina and Nanopore to BAM files, and to sort and index the BAM file reads using the view, sort and index commands. SAMtools merge was then used to merge the index and sorted BAM files into a single BAM file containing both long and short reads for each sample (40).
2.7. Phylogenomic analysis
To performe a comparative phylogenomic analysis of M. tuberculosis isolates, we selected a sequence dataset comprising 94 full-length M. tuberculosis whole genomes, including data from human isolates in Nepal (n = 8), and a genome of M. africanum to use as an outgroup. These genomes were from 33 different countries and the four different geographical lineages previously described from elephants worldwide (lineages 1 and 4) and from humans in Nepal (lineages 1, 2 and 3), and they were retrieved from the NCBI Sequence Read Archive (SRA) database (Supplementary Table S1). The prefetch command in the SRA toolkit v2.11.0 (41) was used to download all whole-genome sequence files of our dataset in SRA format, which were later converted to FASTQ format using the fastq-dump command. Fastp v0.20 was used for quality control and read filtering (37). Read mapping to the reference genome (GenBank accession: NC_000962.3) was performed using BWA v0.7.17 (39). SAMtools v1.15 was used to convert the SAM files of the mapped reads to BAM files, and then to sort and index the BAM file reads (40). The SAMtools mpileup command generated a text pileup output, summarizing the base calls of the aligned reads to the reference sequence. The consensus call was performed using the call command in BCFtools v1.10 (40), which produced a VCF file that was converted to a FASTQ file using the vcf2fq function in the vcfutils.pl script (42). Seqtk v1.3 (42) was used to convert the FASTQ files to FASTA files using the seq command with the masking of bases with quality less than 20.
Variant calling for all the downloaded M. tuberculosis sequences and three M. tuberculosis isolates from this study was performed using PhaME v1.0.4 (43), with the H37Rv genome (GenBank accession: NC_000962.3) as the reference. The PhaME analysis tool was used to remove repeats from each genome in the dataset, and genome alignment to the reference genome was performed using the nucmer2 tool from MUMmer v3.0, as described in the pipeline (43). The output genome alignments were filtered, and SNP tables were generated for each alignment using the delta-filter and show-snps utilities, respectively, both in MUMmer v3.0. The PhaME pipeline then collated the alignment files to produce a core genome alignment and output a multiple alignment FASTA file (~200 kbp) with the variant sites. This was further filtered with SNP-sites to remove non-informative synonymous sites, resulting in a final non-synonymous SNP alignment dataset of 2,596 base pairs.
Using the non-synonymous SNP alignment data, phylogenomic analysis was performed using maximum likelihood and Bayesian inference. For maximum likelihood, the best substitution models and likelihood trees were evaluated using IQTREE v1.6.11 (44) with 500 bootstrap replicates and default parameters. Similarly, for Bayesian inference, the best nucleotide substitution model, TVM, was selected based on the Bayesian Information Criterion score in jModeltest2 v2.1.8 (45, 46), followed by tree reconstruction in MrBayes v3.2.7 (47) with 2,000,000 iterations, discarding 25% as burn-in.
3. Results
3.1. Molecular screening results
The two heat-killed culture samples from elephants yielded 1.07 ng/μL and 0.35 ng/μL DNA for DG and PK, respectively, allowing confirmation of M. tuberculosis by real-time PCR. The M. tuberculosis isolated from the human sample 52G11 was already a confirmed cultured M. tuberculosis with a DNA concentration of 2.34 ng/μL.
3.2. TB-profiler results
Lineage analysis using TB-Profiler showed that the three isolates belonged to three different M. tuberculosis lineages (Table 1). The elephant isolates DG and PK belong to lineage/sub-lineage 1.2.2.2 and 2.2, respectively, whereas the human isolate 52G111 belongs to lineage/sub-lineage 4.1.2.1. None of the three isolates had drug-resistant variants.
Table 1. Mycobacterium tuberculosis lineages, sub-lineages, and drug resistance type of the two elephant and one human isolates collected and sequenced in this study and identified using TB-Profiler v4.1.0.
3.2.1. Whole genome and phylogenomic analysis
The mean coverage of the three M. tuberculosis genomes was 99.6%, while the mean depth was 55.67x, with the isolate DG having the lowest mean depth of 39x (Table 2). The PhaME pipeline generated a core genome alignment of ~200 kbp and a non-synonymous SNP alignment of 2,596 bp size, which included the elephant and human isolates and the genomes downloaded from the NCBI Sequence Read Archive database. The SNP alignment was used for phylogenomic analysis (Table 2).
Table 2. Size, coverage, and mean depth of the analysed M. tuberculosis genomes after merging of the Illumina short reads and Nanopore long reads.
The topology of both maximum-likelihood and Bayesian inference trees were similar, as shown in Figure 2. The phylogenomic analysis of M. tuberculosis sequences revealed that the two elephant-derived isolates, DG and PK, clustered into clades of two different lineages, lineage 1 and lineage 2, respectively. Similarly, the human-derived isolate, 52G111, clustered within into lineage 4 clade. This result is consistent with the lineages/sub-lineages defined by TB-Profiler. Sample DG is evolutionary closely related to two other published genomes of elephant-derived isolates sampled in Nepal (SRR14514371 and SRR14514372). These three elephant-derived sequences clustered together, in the same sub-clade, with published human-derived isolates sampled in Nepal (ERR1213884) and India (SRR5341274). Similarly, sample PK clustered closely in the same clade with published human-derived isolates sampled in Nepal (ERR553171 and ERR551520). The human sample 52G111 clustered in a clade together with other published human-derived samples from Argentina (ERR757160), Russia (ERR067670) and the United Kingdom (ERR047885). The topology of the tree also showed that previously published M. tuberculosis isolates derived from captive Asian elephants in the USA and Switzerland belonged to lineage 4 and clustered in a clade with human-derived isolates from nearby geographical regions. Similarly, an isolate originating from wild African elephants (SRR6487127) clustered very closely with a human-derived isolate (SRR1140950), both sampled in South Africa.
Figure 2. Phylogenomic trees constructed using the 2,596 bp non-synonymous SNP alignment comprising the three new M. tuberculosis genomes described in this study and the 94 genome sequences retrieved from GenBank: (A) IQTREE maximum likelihood tree, the values in the nodes are 500 bootstrap scores; (B) MrBayes Bayesian inference tree, the values in the nodes are posterior probability scores. The elephant (PK, DG) and human (52G111) samples generated in this study are shown in red and orange font, respectively, while the published elephant and human derived samples are shown in blue and black font, respectively.
4. Discussion
In this study, we characterised elephant-derived M. tuberculosis isolates from two diseased animals (PK and DG), and a human-derived M. tuberculosis isolate (52G111) using a whole-genome sequencing approach combining Illumina short reads and Nanopore long reads, resulting in an average genome coverage of 99.6% with an average depth of 55.67x. Phylogenomic analyses were used to explore possible spillover of M. tuberculosis between humans and elephants. We also used an online-based TB profiling tool (TB-Profiler) to further identify lineages/sublineages of the sequenced isolates and to assess drug-resistant variants.
The two elephant-derived M. tuberculosis isolates appear to be of human origin based on both the TB-Profiler result and the phylogenomic trees. These strains are evolutionarily closely related to previous human isolates sampled in Nepal (Figure 1), suggesting a bidirectional transmission from humans to elephants. This may be due to the high prevalence of TB in humans and the continuous exposure of elephants to their handlers, the mahouts, rather than to the sporadic contact with tourists (48). Despite living in the same locality, the two elephant-derived M. tuberculosis isolates clustered in different clades of TB lineages, suggesting a complex and diverse M. tuberculosis transmission dynamics in the region. There must be multiple sources of infection and cross-species (human-mediated) transmission other than direct elephant-to-elephant transmission. This conclusion could be drawn from the fact that lineage 2 is the second most prevalent lineage in Nepal according to human samples from 2009 to 2010, while lineage 1 is the least prevalent according to a study conducted in 2018 (20). The possibility of bidirectional interspecies transmission of M. tuberculosis between humans and elephants has also been supported by previous reports (49–51). Our results reinforce previous evidence suggesting a close association between M. tuberculosis isolates from elephants and humans in Nepal using M. tuberculosis spoligotyping (15, 24), large sequence polymorphism (24) or whole genome data (25).
Tuberculosis is known as a re-emerging disease, but recently it has been considered a zooanthroponosis (52). The observed cases of zooanthroponosis have implications for global and national elephant conservation. Human tuberculosis already contributes to a huge annual burden of morbidity and mortality, and cases are concentrated in countries such as Nepal where poverty and high population density overlap. Multi-drug resistant M. tuberculosis is an additional concern, and the global burden of zoonotic TB is increasing due to the uncontrolled use of antibiotics in animals (53). Zoonotic TB is severely under-reported due to diagnostic challenges and inadequate public health surveillance (54). However, this study attempts to fill this gap in Nepal. Furthermore, the risk factors prevalent in South Asian countries, including Nepal, such as high human-animal density, close and frequent contact with infected animals, inadequate disease control measures, consumption of unpasteurised milk and milk products (55), and use of elephants in tourism (56), explain the possibility of zoonotic TB cases in humans as well as zooanthroponosis in animals. Nevertheless, a serious risk factor analysis should be undertaken to better understand the complex and dynamic transmission of TB in Nepal.
In conclusion, the prospect of an elephant potentially acting as a carrier/reservoir for the transmission of drug-resistant M. tuberculosis as a result of human-animal interactions is a serious concern that may pose future risks not only to humans but also to TB control in livestock and wildlife (57). A One Health approach is fundamental to understanding the spillover or transmission dynamics of infectious diseases such as TB. Therefore, regular screening and detection of these infectious diseases in elephants, other wild and domestic animals, human populations and the environment is essential. The Food and Agriculture Organization of the United Nations (FAO) has already identified in its ‘Roadmap for Zoonotic Tuberculosis’ (58) a wide range of actions involving a wide range of stakeholders to reduce the risk of transmission, thereby bringing economic benefits and improvements in animal welfare. The One Health Strategic Framework 2019 has already been approved in Nepal (59). In Nepal, implementation of the One Health framework appears to be possible only with a focus on zoonoses and zooanthroponosis. Therefore, it should be monitored through surveillance systems and better screening and diagnostic strategies using information on circulating M. tuberculosis genotypes as identified in this study. However, to limit transmission of the pathogen, oral BCG administration provides significant protection against human (60) and animal (61) TB in addition to surveillance prevention strategies. In the context of Nepal, regular M. tuberculosis screening of all captive elephants and people in close contact with them, especially mahouts, should be promoted to control and prevent the spread of TB. This will also help prevent the spread of TB to other endangered species in the forest (62). In addition, good practices and guidelines for tourism-related activities involving direct human-animal interactions, such as mahout-elephant interaction, should be implemented as a preventive measure.
Data availability statement
The datasets presented in this study can be found in online repositories. The raw FASTQ genome data is available in NCBI SRA database under BioProject ID: PRJNA884899.
Ethics statement
The studies involving human participants were reviewed and approved by Ethical approval for the study was obtained from Nepal Ethical Review Board, Nepal Health Research Council (IRC number 312/2018). The patients/participants provided their written informed consent to participate in this study. Ethical review and approval was not required for the animal study because we took the samples from natural deceased animals.
Author contributions
RMR and JQ prepared the overall concept and design. RN, AG, RR, and PM performed the bioinformatics analysis, lab processing, and analysis. AP and NS contributed to the literature review. AS and RMR were involved in sample collection. RMR, NS, CG, DK, PA, JF, and JQ prepared the main writing part of the manuscript and methodology design. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by National Funds through FCT-Fundação para a Ciência e a Tecnologia in the scope of the project UIDP/50027/2020.
Acknowledgments
The authors thank the Department of National Parks and Wildlife Conservation, Government of Nepal for providing research permits to conduct this study. The authors all the staff members of the National Trust for Nature Conservation—Biodiversity Conservation Centre (NTNC-BCC) for their assistance in handling and transportating the samples. The authors thank GENETUP lab for assistance in preparation and confirmation of bacterial cultures. The authors also thank all the staff members of Intrepid Nepal lab especially Bishwo Shrestha, Nisha Shrestha, Dhiraj Puri and Ashok Chaudhary for their assistance during field and laboratory analysis.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2023.1133823/full#supplementary-material
References
1. Burki, T. The global cost of tuberculosis. Lancet Respir Med. (2018) 6:13. doi: 10.1016/S2213-2600(17)30468-X
2. World Health Organization (WHO). Tuberculosis. (2021). Available at: https://www.who.int/news-room/fact-sheets/detail/tuberculosis
3. Chiang, SS, and Starke, JR. Mycobacterium tuberculosis. Princ Pract Pediatr Infect Dis. (2018):790–806.e5. doi: 10.1016/B978-0-323-40181-4.00134-1
4. World Health Organization (WHO). Tackling the drug-resistant TB crisis. (2022). Available at: https://www.who.int/activities/tackling-the-drug-resistant-tb-crisis
5. Wandwalo, E. Why drug-resistant tuberculosis poses a major risk to global health security. The global fund. (2020). Available at: https://www.theglobalfund.org/en/blog/2020-11-19-why-drug-resistant-tuberculosis-poses-a-major-risk-to-global-health-security
6. Nepal National Tuberculosis Centre. Drug resistant tuberculosis management: national guideline. (2019). Available at: https://nepalntp.gov.np/wp-content/uploads/2019/10/DR-TB-Guidelines-2019_Nepal
7. Brites, D, and Gagneux, S. Co-evolution of mycobacterium tuberculosis and Homo sapiens. Immunol Rev. (2015) 264:6–24. doi: 10.1111/imr.12264
8. Malone, KM, and Gordon, SV. Mycobacterium tuberculosis complex members adapted to wild and domestic animals. Adv Exp Med Biol. (2017) 1019:135–54. doi: 10.1007/978-3-319-64371-7_7
9. Talip, BA, Sleator, RD, Lowery, CJ, Dooley, JSG, and Snelling, WJ. An update on global tuberculosis (TB). Infect Dis Res Treat. (2013) 6:39–50. doi: 10.4137/IDRT.S11263
10. Messenger, AM, Barnes, AN, and Gray, GC. Reverse zoonotic disease transmission (zooanthroponosis): a systematic review of seldom-documented human biological threats to animals. PLoS One. (2014) 9:e89055. doi: 10.1371/journal.pone.0089055
11. Miller, M, and Olea-Popelka, F. One health in the shrinking world: experiences with tuberculosis at the human-livestock-wildlife interface. Comp Immunol Microbiol Infect Dis. (2013) 36:263–8. doi: 10.1016/j.cimid.2012.07.005
12. Chastagner, A, Enouf, V, Peroz, D, Hervé, S, Lucas, P, Quéguiner, S, et al. Bidirectional human–swine transmission of seasonal influenza A(H1N1) pdm09 virus in pig herd, France, 2018. Emerg Infect Dis. (2019) 25:1940–3. doi: 10.3201/eid2510.190068
13. Epstein, JH, and Price, JT. The significant but understudied impact of pathogen transmission from humans to animals. Mt Sinai J Med. (2009) 76:448–55. doi: 10.1002/msj.20140
14. Michalak, K, Austin, C, Diesel, S, Bacon, JM, Zimmerman, P, and Maslow, JN. Mycobacterium tuberculosis infection as a zoonotic disease: transmission between humans and elephants. Emerg Infect Dis. (1998) 4:283–7. doi: 10.3201/eid0402.980217
15. Paudel, S, Mikota, SK, Nakajima, C, Gairhe, KP, Maharjan, B, Thapa, J, et al. Molecular characterization of Mycobacterium tuberculosis isolates from elephants of Nepal. Tuberculosis. (2014) 94:287–92. doi: 10.1016/j.tube.2013.12.008
16. Mikota, SK, Peddie, L, Peddie, J, Lindsay, W, Larsen, RS, Iyer, A, et al. Epidemiology and diagnosis of Mycobacterium tuberculosis in captive Asian elephants (Elephas maximus). J Zoo Wildl Med. (2001) 32:1–16. doi: 10.1638/1042-7260(2001)032
17. Zachariah, A, Pandiyan, J, Madhavilatha, GK, Mundayoor, S, Chandramohan, B, Sajesh, PK, et al. Mycobacterium tuberculosis in wild Asian elephants, southern India. Emerg Infect Dis. (2017) 23:504–6. doi: 10.3201/eid2303.161741
18. Angkawanish, T, Wajjwalku, W, Sirimalaisuwan, A, Mahasawangkul, S, Kaewsakhorn, T, Boonsri, K, et al. Mycobacterium tuberculosis infection of domesticated Asian elephants, Thailand. Emerg Infect Dis. (2010) 16:1949–51. doi: 10.3201/eid1612.100862
19. Payeur, JB, Jarnagin, JL, Marquardt, JG, and Whipple, DL. Mycobacterial isolations in captive elephants in the United States. Ann N Y Acad Sci. (2002) 969:256–8. doi: 10.1111/j.1749-6632.2002.tb04388.x
20. Nepal Ministry of Health and Population. National Tuberculosis prevalence survey report 2020. (2020). Available at: https://nepalntp.gov.np/wp-content/uploads/2021/03/NTPS-Report-Bodypages
21. Shrestha, S, and Shrestha, J. Asian elephants and their status in Nepal: a review. J Agric Nat Resour. (2021) 4:227–37. doi: 10.3126/janr.v4i2.33828
22. Pradhan, NMB, Williams, AC, and Dhakal, M. Current status of Asian elephants in Nepal. Gajah. (2011) 35:87–92.
23. Thapa, J, Mikota, SK, Gairhe, KP, Paudel, S, Singh, DK, Dhakal, IP, et al. Tuberculosis seroprevalence and comparison of hematology and biochemistryparameters between seropositive and seronegative captive Asian elephants of Nepal. J Vet Med Sci. (2021) 83:1278–83. doi: 10.1292/JVMS.21-0113
24. Paudel, S, Nakajima, C, Mikota, SK, Gairhe, KP, Maharjan, B, Subedi, S, et al. Mixed Mycobacterium tuberculosis lineage infection in 2 elephants, Nepal. Emerg Infect Dis. (2019) 25:1031–2. doi: 10.3201/eid2505.181898
25. Paudel, S, Brenner, EP, Hadi, SA, Suzuki, Y, Nakajima, C, Tsubota, T, et al. Genome sequences of two Mycobacterium tuberculosis isolates from Asian elephants in Nepal. Microbiol Resour Announc. (2021) 10:e00614–21. doi: 10.1128/MRA.00614-2
26. Martinez, E, Bustamante, A, Menon, R, Wang, Q, Jelfs, P, Marais, B, et al. Whole-genome sequencing of Mycobacterium tuberculosis for rapid diagnostics: feasibility of a decentralised model. Lancet Respir Med. (2016) 4:e13–4. doi: 10.1016/S2213-2600(16)00092-8
27. Chitwan District Information: Chitwan Tourism (2013). Nepal. (2013). Available at: https://www.chitwantourism.com/articles-and-information/chitwan-district-information
28. Lonely Planet. Chitwan National Park Travel—Lonely Planet. Nepal, Asia. (2013). Available at: https://www.lonelyplanet.com/nepal/the-terai-and-mahabharat-range/royal-chitwan-national-park
29. CNP. Managment plan for Chitwan National Park and it’s buffer zone 2013–2017. Chitwan Natl Park Off. (2015) 53:1689–99.
30. CNP. Annual report Chitwan National Park. Government of Nepal. (2022). Available at: https://chitwannationalpark.gov.np
31. Bhusal, NP. Buffer zone management system in protected areas of Nepal. Third Pole J Geogr Educ. (2014) 11:34–44. doi: 10.3126/ttp.v11i0.11558
32. Budhathoki, P. Linking communities with conservation in developing countries: buffer zone management initiatives in Nepal. Oryx. (2004) 38:334–41. doi: 10.1017/S0030605304000584
33. Heinen, JT, and Mehta, JN. Emerging issues in legal and procedural aspects of buffer zone management with case studies from Nepal. J Env Dev. (2000) 9:45–67. doi: 10.1177/107049650000900103
34. Nazarova, EV, and Russell, DG. Growing and handling of Mycobacterium tuberculosis for macrophage infection assays. Methods Mol Biol. (2017) 1519:325–31. doi: 10.1007/978-1-4939-6581-6_22
35. Sacace Biotechnologies S.r.l. Via Scalabrini. Sacace Biotechnologies. (2022). Available at: https://www.sacace.com
36. Phelan, JE, O’Sullivan, DM, Machado, D, Ramos, J, Oppong, YEA, Campino, S, et al. Integrating informatics tools and portable sequencing technology for rapid detection of resistance to anti-tuberculous drugs. Genome Med. (2019) 11:41–7. doi: 10.1186/s13073-019-0650-x
37. Chen, S, Zhou, Y, Chen, Y, and Gu, J. fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. (2018) 34:i884–90. doi: 10.1093/bioinformatics/bty560
38. De Coster, W, D’Hert, S, Schultz, DT, Cruts, M, and Van Broeckhoven, C. NanoPack: visualizing and processing long-read sequencing data. Bioinformatics. (2018) 34:2666–9. doi: 10.1093/bioinformatics/bty149
39. Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. (2013). arXiv. [Epub ahed of preprint], doi: 10.48550/arxiv.1303.3997
40. Danecek, P, Bonfield, JK, Liddle, J, Marshall, J, Ohan, V, Pollard, MO, et al. Twelve years of SAMtools and BCFtools. Gigascience. (2021) 10:1–4. doi: 10.1093/gigascience/giab008
41. NCBI/NLM/NIH. Sequence read archive—download software. (2022). Available at: https://trace.ncbi.nlm.nih.gov/Traces/sra/sra.cgi?view=software
42. GitHub. lh3/seqtk: Toolkit for processing sequences in FASTA/Q formats. (2022). Available at: https://github.com/lh3/seqtk
43. Shakya, M, Ahmed, SA, Davenport, KW, Flynn, MC, Lo, CC, and Chain, PSG. Standardized phylogenetic and molecular evolutionary analysis applied to species across the microbial tree of life. Sci Rep. (2020) 10:1723–15. doi: 10.1038/s41598-020-58356-1
44. Minh, BQ, Schmidt, HA, Chernomor, O, Schrempf, D, Woodhams, MD, Von Haeseler, A, et al. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol. (2020) 37:1530–4. doi: 10.1093/molbev/msaa015
45. Darriba, D, Taboada, GL, Doallo, R, and Posada, D. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. (2012) 9:772. doi: 10.1038/nmeth.2109
46. Guindon, S, and Gascuel, O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. (2003) 52:696–704. doi: 10.1080/10635150390235520
47. Ronquist, F, Teslenko, M, van der Mark, P, Ayres, DL, Darling, A, Höhna, S, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. (2012) 61:539–42. doi: 10.1093/sysbio/sys029
48. Szydlowski, M. Elephants in Nepal: correlating disease, tourism, and welfare. J Appl Anim Welf Sci. (2022) 25:126–38. doi: 10.1080/10888705.2022.2028628
49. Malla, B, Stucki, D, Borrell, S, Feldmann, J, Maharjan, B, Shrestha, B, et al. First insights into the phylogenetic diversity of Mycobacterium tuberculosis in Nepal. PLoS One. (2012) 7:e52297. doi: 10.1371/journal.pone.0052297
50. Abraham, D, and Pillai, VKP. Cross-species transmission of Mycobacterium tuberculosis in mahouts and captive elephants: implications to health policy. Int J Infect Dis. (2016) 45:463. doi: 10.1016/j.ijid.2016.02.981
51. Ong, BL, Ngeow, YF, Razak, MFAA, Yakubu, Y, Zakaria, Z, Mutalib, AR, et al. Tuberculosis in captive Asian elephants (Elephas maximus) in peninsular Malaysia. Epidemiol Infect. (2013) 141:1481–7. doi: 10.1017/S0950268813000265
52. Shimshony, A. Tuberculosis in elephants: a reverse zoonosis. Infect Dis News. (2008) 21:8. Available at: https://www.healio.com/news/infectious-disease/20120225/tuberculosis-in-elephants-a-reverse-zoonosis
53. Acharya, KP, and Wilson, RT. Antimicrobial Resistance in Nepal. Front Med. (2019) 6:105. doi: 10.3389/fmed.2019.00105
54. Kock, R, Michel, AL, Yeboah-Manu, D, Azhar, EI, Torrelles, JB, Cadmus, SI, et al. Zoonotic tuberculosis—the changing landscape. Int J Infect Dis. (2021) 113:S68–72. doi: 10.1016/j.ijid.2021.02.091
55. Bapat, PR, Dodkey, RS, Shekhawat, SD, Husain, AA, Nayak, AR, Kawle, AP, et al. Prevalence of zoonotic tuberculosis and associated risk factors in central Indian populations. J Epidemiol Glob Health. (2017) 7:277–83. doi: 10.1016/j.jegh.2017.08.007
56. Szydlowski, M. Framing conservation, colonialism and care: Captive endangered Asian elephants (Elephas maximus) in Nepal. (2021). Doctoral Theses. Available at: https://ore.exeter.ac.uk/repository/handle/10871/127765
57. Botelho, A, Perdigão, J, Canto, A, Albuquerque, T, Leal, N, Macedo, R, et al. Pre-multidrug-resistant Mycobacterium tuberculosis Beijing strain associated with disseminated tuberculosis in a pet dog. J Clin Microbiol. (2014) 52:354–6. doi: 10.1128/JCM.02800-13
58. Food and Agriculture Organization of the United Nations (FAO). Roadmap for zoonotic tuberculosis. (2017). Available at: https://www.fao.org/3/i7807e/i7807e
59. One Health Strategy, 2076. Ministry of Agriculture and Livestoke Development, Nepal. (2022). Available at: https://publichealthupdate.com/one-health-strategy-2076
60. Brewer, TF. Preventing tuberculosis with bacillus Calmette–Guérin vaccine: a meta-analysis of the literature. Clin Infect Dis. (2000) 31 Suppl 3:S64–7. doi: 10.1086/314072
61. Buddle, BM, Vordermeier, HM, Chambers, MA, and de Klerk-Lorist, LM. Efficacy and safety of BCG vaccine for control of tuberculosis in domestic livestock and wildlife. Front Vet Sci. (2018) 5:259. doi: 10.3389/fvets.2018.00259
Keywords: Elephas maximus, whole-genome sequencing, one health, tuberculosis, zooanthroponosis
Citation: Rajbhandari RM, Napit R, Manandhar P, Raut R, Gurung A, Poudel A, Shrestha N, Sadaula A, Karmacharya D, Gortázar C, Alves PC, de la Fuente J and Queirós J (2023) Phylogenomic analysis supports Mycobacterium tuberculosis transmission between humans and elephants. Front. Vet. Sci. 10:1133823. doi: 10.3389/fvets.2023.1133823
Edited by:
Jose Angel Barasona, Complutense University of Madrid, SpainReviewed by:
Azlan Che-Amat, Putra Malaysia University, MalaysiaDavid González-Barrio, Carlos III Health Institute (ISCIII), Spain
Copyright © 2023 Rajbhandari, Napit, Manandhar, Raut, Gurung, Poudel, Shrestha, Sadaula, Karmacharya, Gortázar, Alves, de la Fuente and Queirós. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: João Queirós, joao.queiros@cibio.up.pt