- Short report
- Open Access
Multi-locus sequence analysis reveals great genetic diversity among Mycoplasma capricolum subsp. capripneumoniae strains in Asia
Veterinary Research volume 53, Article number: 92 (2022)
Multi-Locus Sequence Analysis (MLSA) of Mycoplasma capricolum subsp. capripneumoniae (Mccp) strains from Asia revealed unforeseen diversity and a central position for genotyping groups representing strains from Central/East Asia, suggesting a possible origin of contagious caprine pleuropneumonia in this continent. A better assessment of the emergence, diversity and distribution of Mccp in Asia and Africa calls for renewed efforts to dramatically enlarge the sample of strains. Availability and affordability in the field, added to superior typeability (directly from poor samples) and high stability, discriminatory power and concordance with epidemiological and phylogenetic analyses, make MLSA an excellent tool for such investigations.
Introduction, methods and results
Contagious caprine pleuropneumonia (CCPP) is a devastating disease affecting domestic goats and several wild ungulate species in arid and semiarid regions of Africa, Middle East and Asia, where goat rearing plays an essential role in food security and poverty alleviation . Owing to its high contagiousness, morbidity and mortality, CCPP is included in the list of notifiable diseases of the World Organisation for Animal Health (WOAH, founded as OIE; ). Its etiologic agent, a fastidious bacterium known as Mycoplasma capricolum subsp. capripneumoniae (Mccp), is very rarely isolated and CCPP is hardly ever reported. As a consequence, the distribution, prevalence and impact of CCPP are not well established .
To improve our understanding on the epidemiology of CCPP, a molecular typing scheme based on the analysis of eight genetic markers, known as Multi-Locus Sequence Analysis (MLSA), was developed in 2011 . This tool was extremely robust and allowed genotyping directly from infected tissues from which Mccp could not be isolated. The scheme was applied to 27 strains of diverse origins, resulting in the identification of two lineages and 5 groups, which were correlated to the geographic origin of the strains (with the remarkable exception of the Arabian Peninsula, where strains from 4 out of the 5 groups were found). Notably, the identification of a distinct Asian cluster represented by two recent strains from Tajikistan and China (sole representatives of Central and East Asia available at the time) indicated a local evolution of strains and excluded a recent introduction of CCPP in the continent.
Thanks to the democratisation of high throughput sequencing technologies more sophisticated Mccp strain genotyping methods have been developed, from a multi-gene scheme  to a whole-genome sequence (WGS) analysis pipeline , attaining optimum strain typing for molecular epidemiology studies and outbreak investigations. However, WGS-based genotyping is not available to diagnostic laboratories, particularly in the regions where CCPP is prevalent, and MLSA may still be a valuable alternative, especially when isolation cannot be achieved. Only a few Mccp isolates and WGS have been made available since the MLSA work of 2011 and subsequent reports relating to Mccp strains from wildlife in the United Arab Emirates [7, 8], but MLSA has been conducted following investigations of CCPP outbreaks in Tibetan wild ungulates first identified in 2012  and, more recently, in Pakistani goats in 2019 . The objective of our study was thus to explore the diversity of Mccp strains in Asia, by analysing new MLSA data from Pakistan and China, including strains originating from wildlife. This was also the opportunity to update the global Mccp MLSA, by including all the data generated since 2011, and to analyse its value and performance in comparison to subsequent typing techniques based on WGS data.
The 43 Mccp strains and/or corresponding genomic sequences analysed in this study are presented in Table 1, including 8 strains from wild ungulate species. Thirty-three of them were included in subsequent typing schemes [5, 6] and corresponding phylogenetic groups are presented when available. Sixteen new strains were added to 27 previously published . MLSA data of 6 new strains were extracted from available WGS, while the remaining 10 were obtained by PCR amplification and sequencing of the corresponding eight loci as previously described , with the exception that Sanger sequencing was performed by Macrogen (South Korea), while Geneious 10.2.6  was used for sequence assembly and alignment.
The sequences of epidemiologically-related strains collected in nearby locations during CCPP epizootics in Uganda, Tunisia and Tibet or obtained by in vitro passage (Table 1) were identical, showing that the molecular markers were stable and there were no laboratory-introduced variations. Furthermore, MLSA results obtained by locus amplification and sequencing versus extraction from WGS data (for 15 strains analysed by PCR and sequencing with WGS available in GenBank, Table 1) were also identical. The only exception was strain F38, for which a single nucleotide polymorphism (SNP) in the H2 locus differentiated MLSA sequences obtained by the two methods. However, since two different laboratory stocks of this strain were used for PCR and sequencing (CIRAD) and WGS (NCTC 10192 T), this SNP may result from divergent evolution undergone by the two laboratory stocks from the original 1974 isolate . When the scheme was applied to the remaining 39 “unrelated” strains in Table 1, 24 sequence types (ST) (9 new) were discriminated based on 68 polymorphic positions (16 new), which are shown in Table 2, with locus sequences from Mccp type strain F38 serving as reference. This resulted in a Simpson’s index of diversity of 0.970 (0.953–0.987), which expresses the probability of two unrelated strains being characterised as the same type [13, 14]. All the strains could be discriminated individually by WGS analysis  and all but two of those analysed by Dupuy et al.  provided distinct genotypes (Table 1). However, these two isolates were actually discriminated by MLSA, which allowed typing of non-viable strains (n = 6, Table 1) with no added difficulty or cost.
A robust tree (Figure 1) was obtained by distance analysis of MLSA data using DARwin 6  as previously described. Seven genotyping groups were identified, distributed in the two lineages previously described. Pre-existing MLSA groups 1–5 were unchanged, with the exception of several additional ST identified in group 1, corresponding to East African and Emirati strains originating from domestic goat and wild ungulates respectively. The remaining new ST identified in this work corresponded to Asian strains and were clustered in two additional groups, positioned within lineage II. A highly variable cluster located near the centre of the tree and represented by Chinese strains from Shandong and Tibet was designated group 6, whereas the Pakistani strain constituted the single representative of group 7. All Asian strains (disregarding those originating from the Middle East) were found spread among three clusters (groups 3, 6 and 7) within lineage II, together with group 4 (represented by strains from North Africa, the Arabian Peninsula and Turkey) and group 5 (comprising mainly East African strains). As shown in Figure 2, a generally good correlation between ST and geographic origin was retained, with the exception of the Arabian Peninsula, where animals from diverse origins are imported every year, particularly at the occasion of Muslim feasts . A similar situation was now observed in Turkey, since strains from Thrace and Elazig (East Turkey) were positioned in groups 3 and 4 respectively.
The relevance of the MLSA scheme for Mccp genotyping and epidemiology analyses is unquestionable, particularly when we consider its accessibility, affordability, ease of use and superior typeability, allowing direct genotyping from poor samples. Furthermore, its stability, regardless the method used to obtain the data, was remarkable and MLSA clustering was highly congruent with both epidemiological and phylogenetic analyses [5, 6]. Finally, its high discriminatory power was compatible with epidemiological investigations.
Analysis of new strains from Pakistan and China allowed a better representation of the spread of CCPP in Asia (Figure 2) and revealed unpredicted diversity in this continent (Figure 1). The Pakistani strain was the sole representative of a new cluster (group 7), the diversity and distribution of which remain to be disclosed. This was unfortunately the only strain available from South Asia, where the occurrence of CCPP was documented as early as 1914  and where CCPP is known to be prevalent [10, 17, 18]. The new Chinese strains constituted a distinct cluster (group 6), separate from previously described Tajik and Chinese strains (group 3). A strain from Tibetan sheep collected in the Nagqu region of Tibet (Table 1), where devastating CCPP outbreaks have been reported in both domestic goat and antelope since 2012 ,was placed at the base of this group. This strain was more closely related to strains from domestic goats collected at Shandong than to strains from Tibetan antelope collected at Nagqu. This may be explained by the wide area of distribution of domestic and wild ungulate species across the Qinghai-Tibetan plateau and its peripheral mountains. It was assumed that Tibetan antelopes were infected due to close contact with domestic goats, which are progressively invading their habitat . Furthermore, strains from CCPP outbreaks affecting domestic sheep in Uganda  and four different wild ungulate species in the Middle East [7, 8, 20] (Table 1), were placed in group 1, very distant to those from Tibetan wildlife, and shared or were closely related to ST from goat isolates, indicating that the same strains can affect a wide variety of species. Again, the assumption was that domestic goats were the source of the infection in sheep and wildlife, though direct Mccp transmission among infected wild ungulates of different species has been demonstrated, at least in captivity [8, 20].
Analysis of new Emirati strains from wildlife and additional strains from East Africa resulted in the identification of three new ST in group 1, revealing greater diversity for this cluster, which is spreading in the region. The strain introduced in Mauritius in 2009 , shared ST with a highly virulent Kenyan isolate from 2012 [22, 23] and was closely related to a strain that was responsible for CCPP outbreaks across Tanzania in 2013 . The relatively low diversity of this group, and generally of lineage I compared to lineage II, deserves further investigation. Similarly, the presence in East Africa of two distant genotyping groups (one from each lineage), suggesting two different introductions of CCPP in this region, needs to be elucidated for a better understanding of the origin and evolution of CCPP in Africa.
CCPP was suspected in India and China since the beginning of the twentieth century [16, 24], but its presence in Asia was only confirmed in 2007 . Already in 2011, MLSA genotyping suggested that CCPP was present for a long time in Asia , which was substantiated by subsequent large-scale genomic analyses [5, 6]. The great genetic diversity observed here among Asian Mccp strains in spite of the limited number of samples analysed, together with the position of MLSA groups 3 and 6 (represented by Central and East Asian strains) at the centre of the tree, point towards a possible origin of CCPP in Asia. Again, the scarcity of Mccp strains hampers a precise determination of the emergence, diversity and distribution of Mccp.
A better assessment of the molecular evolution and epidemiology of CCPP in Asia and Africa calls for renewed efforts to dramatically enlarge the sample of strains from diverse origins representing the real distribution of CCPP, which is yet to be established (Figure 2). MLSA can be an excellent tool to do this, provided CCPP cases are investigated, since these analyses can be achieved from simple samples such as dried filter paper imbedded in infected material, which can be easily stored and shipped at room temperature.
Availability of data and materials
The locus sequences corresponding to new MLSA sequence types 1–011, 1–040, 6–010, 6–020 and 7–010 obtained in this study, for which no representative sequences are available, were submitted to GenBank (accession numbers: OP076701–OP076740).
Manso-Silván L, Thiaucourt F (2019) Contagious caprine pleuropneumonia. In: Kardjadj M, Diallo A, Lancelot R (eds) Transboundary animal diseases in sahelian Africa and connected regions. Springer Nature Switzerland AG, Cham
WOAH (2022) World Organization for Animal Health https://www.woah.org Accessed 22 Sept 2022
WOAH (2021) Contagious caprine pleuropneumonia. In: Manual of diagnostic tests and vaccines for terrestrial animals of the world organization for animal health https://www.woah.org/fileadmin/Home/eng/Health_standards/tahm/3.08.04_CCPP.pdf. Accessed 22 Sept 2022
Manso-Silván L, Dupuy V, Chu Y, Thiaucourt F (2011) Multi-locus sequence analysis of Mycoplasma capricolum subsp. capripneumoniae for the molecular epidemiology of contagious caprine pleuropneumonia. Vet Res 42:86
Dupuy V, Verdier A, Thiaucourt F, Manso-Silván L (2015) A large-scale genomic approach affords unprecedented resolution for the molecular epidemiology and evolutionary history of contagious caprine pleuropneumonia. Vet Res 46:74
Loire E, Ibrahim AI, Manso-Silván L, Lignereux L, Thiaucourt F (2020) A whole-genome worldwide molecular epidemiology approach for contagious caprine pleuropneumonia. Heliyon 6:e05146
Chaber AL, Lignereux L, Al Qassimi M, Saegerman C, Manso-Silván L, Dupuy V, Thiaucourt F (2014) Fatal transmission of contagious caprine pleuropneumonia to an Arabian oryx (Oryx leucoryx). Vet Microbiol 173:156–159
Lignereux L, Chaber AL, Saegerman C, Manso-Silván L, Peyraud A, Apolloni A, Thiaucourt F (2018) Unexpected field observations and transmission dynamics of contagious caprine pleuropneumonia in a sand gazelle herd. Prev Vet Med 157:70–77
Yu Z, Wang T, Sun H, Xia Z, Zhang K, Chu D, Xu Y, Xin Y, Xu W, Cheng K, Zheng X, Huang G, Zhao Y, Yang S, Gao Y, Xia X (2013) Contagious caprine pleuropneumonia in endangered Tibetan antelope, China, 2012. Emerg Infect Dis 19:2051–2053
Ahmad F, Khan H, Khan FA, Carson BD, Sadique U, Ahmad I, Saeed M, Rehman FU, Rehman HU (2020) The first isolation and molecular characterization of Mycoplasma capricolum subsp. capripneumoniae Pakistan strain: a causative agent of contagious caprine pleuropneumonia. J Microbiol Immunol Infect 54:710–717
Geneious 10.2.6 (https://www.geneious.com) Accessed 22 Sept 2022.
MacOwan KJ, Minette JE (1976) A mycoplasma from acute contagious caprine pleuropneumonia in Kenya. Trop Anim Health Prod 8:91–95
Hunter PR, Gaston MA (1988) Numerical index of the discriminatory ability of typing systems: an application of Simpson’s index of diversity. J Clin Microbiol 26:2465–2466
Grundmann H, Hori S, Tanner G (2001) Determining confidence intervals when measuring genetic diversity and the discriminatory abilities of typing methods for microorganisms. J Clin Microbiol 39:4190–4192
DARwin (2022) Dissimilarity Analysis and Representation for Windows https://darwin.cirad.fr/ Accessed 22 Sept 2022.
Walker GK (1914) Pleuro-pneumonia of goats in the Kangra district, Punjab, India. J Comp Pathol 27:68–71
Awan MA, Abbas F, Yasinzai M, Nicholas RA, Babar S, Ayling RD, Attique MA, Ahmed Z, Wadood A, Khan FA (2009) First report on the molecular prevalence of Mycoplasma capricolum subspecies capripneumoniae (Mccp) in goats the cause of contagious caprine pleuropneumonia (CCPP) in Balochistan province of Pakistan. Mol Biol Rep 37:3401–3406
Parray OR, Yatoo MI, Muheet BRA, Malik HU, Bashir ST, Magray SN (2019) Seroepidemiology and risk factor analysis of contagious caprine pleuropneumonia in Himalayan Pashmina Goats. Small Ruminant Res 171:23–36
Bolske G, Johansson KE, Heinonen R, Panvuga PA, Twinamasiko E (1995) Contagious caprine pleuropneumonia in Uganda and isolation of Mycoplasma capricolum subspecies capripneumoniae from goats and sheep. Vet Rec 137:594–594
Arif A, Schulz J, Thiaucourt F, Taha A, Hammer S (2007) Contagious caprine pleuropneumonia outbreak in captive wild ungulates at Al Wabra Wildlife Preservation, State of Qatar. J Zoo Wildl Med 38:93–96
Srivastava AK, Meenowa D, Barden G, Salguero FJ, Churchward C, Nicholas RA (2010) Contagious caprine pleuropneumonia in Mauritius. Vet Rec 167:304–305
Falquet L, Liljander A, Schieck E, Gluecks I, Frey J, Jores J (2014) Complete genome sequences of virulent Mycoplasma capricolum subsp. capripneumoniae strains F38 and ILRI181. Genome Announc 2:e01041-e1114
Liljander A, Sacchini F, Stoffel MH, Schieck E, Stokar-Regenscheit N, Labroussaa F, Heller M, Salt J, Frey J, Falquet L, Goovaerts D, Jores J (2019) Reproduction of contagious caprine pleuropneumonia reveals the ability of convalescent sera to reduce hydrogen peroxide production in vitro. Vet Res 50:10
Animal Husbandry Department of Xinjiang Uygur Autonomous Region (1995) Investigation and control of animal epidemics in Xinjiang. Xinjiang Science and Technology Health Press, Urumqi, pp 495–503
Li Y, Zhang JH, Hu SP, Wang L, Xin JQ (2007) Reclassification of the four China isolated strains of the pathogen for contagious caprine pleuropneumonia. Wei Sheng Wu Xue Bao 47:769–773
Lorenzon S, Wesonga H, Ygesu L, Tekleghiorgis T, Maikano Y, Angaya M, Hendrikx P, Thiaucourt F (2002) Genetic evolution of Mycoplasma capricolum subsp. capripneumoniae strains and molecular epidemiology of contagious caprine pleuropneumonia by sequencing of locus H2. Vet Microbiol 85:111–123
Lefevre PC, Breard A, Alfarouk I, Buron S (1987) Mycoplasma species F 38 isolated in Chad. Vet Rec 121:575–576
Chu Y, Yan X, Gao P, Zhao P, He Y, Liu J, Lu Z (2011) Molecular detection of a mixed infection of Goatpox virus, Orf virus, and Mycoplasma capricolum subsp. capripneumoniae in goats. J Vet Diagn Invest 23:786–789
Chen S, Hao H, Zhao P, Thiaucourt F, He Y, Gao P, Guo H, Ji W, Wang Z, Chu LuZ, Y, Liu Y, (2017) Genome-wide analysis of the first sequenced Mycoplasma capricolum subsp capripneumoniae strain M1601. G3 (Bethesda) 7:2899–2906
Ozdemir U, Ozdemir E, March JB, Churchward C, Nicholas RA (2005) Contagious caprine pleuropneumonia in the Thrace region of Turkey. Vet Rec 156:286–287
Amirbekov M, Murvatulloev S, Ferrari F (2010) Contagious caprine pleuropneumonia detected for the first time in Tajikistan. EMPRESS Transbound Anim Bul 35:20–22
Perreau P, Breard A, Le Goff C (1984) Experimental infection in goats caused by mycoplasma strain F.38 (contagious caprine pleuropneumonia). Ann Microbiol (Paris) 135A:119–124
Cetinkaya B, Kalin R, Karahan M, Atil E, Manso-Silvan L, Thiaucourt F (2009) Detection of contagious caprine pleuropneumonia in East Turkey. Rev Sci Tech 28:1037–1044
Jones GE, Wood AR (1988) Microbiological and serological studies on caprine pneumonias in Oman. Res Vet Sci 44:125–131
Harbi MS, Mageed IA, El Tahir MS (1983) Experimental contact transmission of contagious pleuropneumonia of goats (Abu Nini) in the Sudan. Vet Res Commun 6:139–143
Thiaucourt F, Breard A, Lefevre PC, Mebratu GY (1992) Contagious caprine pleuropneumonia in Ethiopia. Vet Rec 131:585
Dupuy V, Thiaucourt F (2014) Complete genome sequence of Mycoplasma capricolum subsp. capripneumoniae strain 9231-Abomsa. Genome Announc 2:e01067-e1114
Xin J-q, Li Y, Zhang J-h, Hu S-p, Wang L (2007) Molecule characterization of a strain of Mycoplasma capricolum subsp. capripnemonia isolated from goat. Chinese J of Prev Vet Med 29:243–248 (in Chinese)
Li Y, Wang R, Sun W, Song Z, Bai F, Zheng H, Xin J (2020) Comparative genomics analysis of Mycoplasma capricolum subsp. capripneumoniae 87001. Genomics 112:615–620
Hao H, Chen S, Li Y, Sun H, Zhao P, Jian Y, Gao Y, Wu C, Liu Y, Chu Y (2017) Complete genome sequence of Mycoplasma capricolum subsp. capripneumoniae strain zly1309F, isolated from endangered Tibetan Antelope. Genome Announc 5:e00496-e517
The authors wish to thank the Livestock & Dairy Development (L&DD, Khyber Pakhtunkhwa) and Veterinary Research Institute, Peshawar, for provision of samples, as well as Yuwei Gao, Key laboratory of Jilin Province for Zoonosis, Prevention and Control, Changchun, Jilin, China, for facilitating sampling of Tibetan animals.
This work was financially supported by the French Agricultural Research Centre for International Development (CIRAD) in France; the LVRI Yuan-Heng Talent Programme (NKLS2020-119), the Key Talents Program of Gansu Province(2021RCXM047)and the Natural Science Foundation of Gansu Province (20JR5RA583) in China; the joint research project of the University of Agriculture, Peshawar and Sandia National Laboratories, New Mexico, USA under the Pak-US Science and Technology Cooperation Program, Phase 7, 2017. This program is supported and implemented by the National Academy of Sciences in the USA and by the Higher Education Commission in Pakistan.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Akhtar, A., Boissière, A., Hao, H. et al. Multi-locus sequence analysis reveals great genetic diversity among Mycoplasma capricolum subsp. capripneumoniae strains in Asia. Vet Res 53, 92 (2022). https://doi.org/10.1186/s13567-022-01107-z
- Mycoplasma capricolum subsp. capripneumoniae
- contagious caprine pleuropneumonia
- multi-locus sequence analysis
- molecular epidemiology