A large-scale genomic approach affords unprecedented resolution for the molecular epidemiology and evolutionary history of contagious caprine pleuropneumonia
© Dupuy et al. 2015
Received: 10 February 2015
Accepted: 18 May 2015
Published: 6 July 2015
Contagious caprine pleuropneumonia (CCPP), caused by Mycoplasma capricolum subsp. capripneumoniae (Mccp), is a devastating disease of domestic goats and of some wild ungulate species. The disease is currently spreading in Africa and Asia and poses a serious threat to disease-free areas. A comprehensive view of the evolutionary history and dynamics of Mccp is essential to understand the epidemiology of CCPP. Yet, analysing the diversity of genetically monomorphic pathogens, such as Mccp, is complicated due to their low variability. In this study, the molecular epidemiology and evolution of CCPP was investigated using a large-scale genomic approach based on next-generation sequencing technologies, applied to a sample of strains representing the global distribution of this disease. A highly discriminatory multigene typing system was developed, allowing the differentiation of 24 haplotypes among 25 Mccp strains distributed in six genotyping groups, which showed some correlation with geographic origin. A Bayesian approach was used to infer the first robust phylogeny of the species and to date the principal events of its evolutionary history. The emergence of Mccp was estimated only at about 270 years ago, which explains the low genetic diversity of this species despite its high mutation rate, evaluated at 1.3 × 10−6 substitutions per site per year. Finally, plausible scenarios were proposed to elucidate the evolution and dynamics of CCPP in Asia and Africa, though limited by the paucity of Mccp strains, particularly in Asia. This study shows how combining large-scale genomic data with spatial and temporal data makes it possible to obtain a comprehensive view of the epidemiology of CCPP, a precondition for the development of improved disease surveillance and control measures.
Contagious caprine pleuropneumonia (CCPP) is a severe respiratory disease affecting goats and some wild ruminant species. The disease, listed by the World Organisation for Animal Health (OIE)  has a great economic impact on livestock production in fragile rural economies and poses a serious threat to disease-free areas. CCPP is caused by Mycoplasma capricolum subsp. capripneumoniae (Mccp), a member of the Mycoplasma mycoides cluster . This cluster comprises five mycoplasmas which are pathogenic for ruminants, including Mycoplasma capricolum subsp. capricolum (Mcc), the closest relative of Mccp, and Mycoplasma mycoides subsp. mycoides (Mmm), the agent of contagious bovine pleuropneumonia (CBPP). Since it was first isolated in 1976 , Mccp has only been isolated in 17 countries, mainly because of its fastidiousness in culture. However, clinical descriptions have been published in nearly 40 countries in Africa and Asia, suggesting a much wider distribution . The disease is present in the Arabian Peninsula, North, Central and East Africa and Asia, but its boundaries are still uncertain, particularly in western and southern Africa and in Asia. In the last decade, an increasing number of outbreaks have been reported in both domestic and wild ruminants . New detections were often the result of improved diagnosis, confirming the presence of the disease in suspected regions [6–8], but in certain cases, they indicated that the disease had spread to new territories [9,10].
A comprehensive view of the evolutionary history and the dynamics of Mccp strains is essential to understand the epidemiology of CCPP. Molecular epidemiology investigations require robust genotyping tools with sufficient resolution, while evolutionary analyses are also constrained by the need for a reliable molecular clock to infer evolutionary timescales. Various molecular methods have been developed for the analysis of Mccp strains. The first study on the molecular evolution of Mccp was based on the 16S rRNA gene , which provided a basis for the phylogeny and systematics of bacteria since the evolutionary studies of Woese . Mccp strains showed a surprisingly high degree of polymorphism between their two 16S rRNA gene copies , which allowed the use of this housekeeping gene to combine epidemiological and evolutionary analyses. Still, the study was limited by the low discriminatory power of this molecular marker. Despite the heterogeneity observed in the sequence of its two rRNA operons, the Mccp genome has been shown to be rather monomorphic. Thus, a phylogenic study of the M. mycoides cluster based on five partial housekeeping gene sequences showed a very low distance between Mccp strains . The low diversity of such a monomorphic bacterium precluded the use of standard multilocus sequence typing (MLST), which relies on housekeeping genes . Approaches based on alternative sequences, independently of their coding capacity, were thus preferred. Analysis of the H2 locus allowed the discrimination of four groups showing a good correlation with geographic origin . This system was improved by the addition of seven loci, increasing the discriminatory power and leading to the description of five groups . Although this multilocus sequence analysis (MLSA) scheme targets polymorphic loci by focusing on non-housekeeping genes (e.g., intergene regions, pseudogenes) it interrogates less than 1% of the genome and therefore has limited discriminatory power for epidemiological investigations. Moreover, it is not reliable for phylogenetic analyses because of differences in the molecular evolutionary clock among target sequences.
These limitations can now be overcome by the recent development of high-throughput methods . Even though single nucleotide polymorphism (SNP) frequencies are low in monomorphic pathogens, their numbers can be dramatically increased by enlarging the scale of the analysis. Using high-throughput data can enhance genetic investigations by providing a way to disclose genome-wide variations. Combining large-scale genomic data with spatial and temporal data already enabled a comprehensive view of the molecular epidemiology and evolution of bacterial pathogens like Salmonella typhi , Yersinia pestis , Mmm  and Mycobacterium tuberculosis . Following the increasing use and affordability of DNA sequencing with next-generation sequencing (NGS) technologies, complete annotated genomes of several Mccp strains have become available [22–24], thereby making large-scale genomic investigations possible.
In this study, we investigated the molecular epidemiology and evolution of CCPP using a large-scale genomic approach based on NGS data, on a sample of strains representing the global distribution of this disease. Our main objective was to develop a discriminatory genotyping method to investigate the genetic diversity and population structure of Mccp. A robust phylogeny was also inferred from a large phylogenomic data set and divergence time of Mccp strains was estimated to reconstruct the evolutionary history of CCPP.
Materials and methods
List of Mycoplasma capricolum subsp. capripneumoniae strains analysed and corresponding MLSA types
Doha, Al Wabra
Sample preparation and sequencing
Twenty-one strains for which the genome sequence was not available were cultured in modified Hayflick’s medium  at 37 °C, 5% CO2. Culture purity was ensured by phenotypic control on solid medium and specific Mccp QPCR amplification . DNA was extracted using a standard phenol/chloroform method . DNA purity, quality and quantity were checked using NanoDrop™ ND-1000 Spectrophotometer (Thermo Fisher Scientific, MA, USA), gel electrophoresis and Qubit® 2.0 fluorometer (Invitrogen, USA), respectively. Then, 21 tagged standard genomic libraries were constructed and pooled to be sequenced in 100 bp single reads on an Illumina HiSeq 2000 (GATC, Constanz, Germany).
The choice of Mccp genes was based on the choice previously made for the analysis of the evolutionary history of Mmm , consisting in 62 rigorously selected genes. Pseudogenes and duplicated genes had been excluded from this set, as well as genes coding for membrane proteins or restriction enzymes, and those known to be involved in horizontal transfer. The annotated, circularised genome sequence of strain 9231-Abomsa  (Table 1) was used as reference. Four genes did not exist in the Mccp genome (guaC, gntR, suk, bgl), while dnaC was duplicated and was therefore excluded. As a result, a subset of 57 genes, comprising 47 coding sequences and 10 pseudogenes of the core genome, were used (Additional file 1). The 57 genes are evenly distributed along the chromosome of strain 9231-Abomsa (Additional file 2).
Data set collection
The sequences of the 57 selected genes from strain 9231-Abomsa, including flanking regions (up to 350 bp), were concatenated and annotated, resulting in an “enlarged” sequence of 107 050 bp (Additional file 3). This sequence was used as reference to automatically retrieve the entire gene set from the whole genome sequence data of 21 Mccp strains by mapping raw data using Seqman NGen (2.0) software (DNASTAR, Madison WI, USA). First, this procedure allowed the correct mapping of reads on the entire sequence of each corresponding gene, thanks to the presence of flanking regions. Second, it allowed the visual verification of sequencing depth and the identification of any incongruities on all coding sequences in the Seqman genome browser. On average, a read depth of 500X was obtained. SNPs and indels were called when more than 85% of the reads supported the change.
An in-house software was developed at CIRAD Montpellier to extract and concatenate the gene sequences corresponding to each strain from the enlarged reference sequence. Sequence searches are based on the Needleman-Wunsch algorithm, using tags to frame sequences of interest. These tags are short sequences homologous to the extremities of each target gene. The algorithm described here is implemented in C++ as a stand-alone program, “SelectRegion”, and the source is freely available from the authors on request. The process can also be automated using graphical-interfaces within the web-based Galaxy .
Otherwise, sequence data were retrieved from the published genome sequences of strains M1601, 99108 and ILRI181 (Table 1).
For diversity analyses, the 57 genes (comprising both coding sequences and pseudogenes) corresponding to each of the 25 strains analysed were selected and concatenated as described above. Sequences were aligned using Clustal W with default parameters (Additional file 4). Haplotypes were estimated taking into account all sites, including gaps, removing invariant sites, with DnaSP . A median-joining network was reconstructed using NETWORK V4.6 .
The discriminatory power of the genotyping system was calculated using Simpson’s index of diversity , which expresses the probability of two unrelated strains sampled from the test population being placed into different genotyping groups. A total of 95% confidence intervals (CI) were determined as previously described .
Phylogeny and molecular dating analyses
From the 57 genes initially selected, a subset of 47 coding sequences was retained for evolutionary analyses. The 47 genes of 25 Mccp strains were extracted from each corresponding enlarged reference sequence using the internal software SelectRegion. The type strain of the subspecies Mcc (California KidT), which is the closest relative of Mccp, was chosen as outgroup and corresponding sequence data were obtained from the published genome sequence [Genbank:CP000123].
Sequences were aligned using Clustal W with default parameters (Additional file 5). The Tamura-Nei, 93 (TN93) model was selected as the best fitting model by Modeltest V3.7 . A maximum-likelihood phylogenetic tree was inferred using PhyML V3.0  on the Galaxy web platform. A bootstrap resampling procedure with 1000 replicates was used to assess the reliability of key tree nodes. To infer a temporal framework from dated sequences, a Bayesian approach, implemented in the flexible Bayesian phylogenetic analysis package BEAST V1.6.2 , was used. This allowed the simultaneous estimation of the tree structure, the time of the most recent common ancestor (MRCA), the divergence time of nodes, and the mutation rate. The TN93 evolution model and partition codon (1 + 2, 3) were selected. Two molecular clocks (strict-clock and uncorrelated lognormal-clock) and various demographic models (constant, expansion, exponential and extended Bayesian skyline plot) were tested. Convergence was evaluated using Tracer V1.6. To choose the best fitting clock and demographic model, the Bayes Factors [35,36] were calculated from the marginal likelihoods on the Akaike’s information criterion through MCMC (AICM) using Tracer V1.6. The Maximum Clade Credibility tree was constructed using TreeAnnotator V1.6 and visualised using FigTree V1.3.
Genetic diversity and molecular genotyping of Mccp strains
A discriminatory genotyping system based on a large-scale genomic approach was developed to characterise the diversity of Mccp by analysing 25 strains representing the known global distribution of this species (Table 1). A set of 57 genes, comprising 47 coding sequences and 10 pseudogenes (Additional file 1), was analysed, covering 77 898 base pairs (7.7% of the genome). Fifty-two of these genes were polymorphic as a result of 239 polymorphic positions consisting in 212 SNPs and 17 indels (Additional file 6). Nine events corresponded to indels of either one/two bases or one/two codons and eight events to variations in homopolymer size. The average frequency of polymorphic sites in the gene set was of one event per 577 pb for coding sequences, and one event per 210 pb for pseudogenes. The 239 polymorphic positions made it possible to define 24 haplotypes among the 25 strains analysed, resulting in a Simpson’s diversity index of 0.997 (0.988-1.000). Polymorphic positions consisted in 115 informative sites and 124 sites specific to single haplotypes.
Three groups (A, B and F) presented at least four isolates originating from the same geographic region, whereas too few strains from groups D and E were available to confirm a correlation between genotyping group and geographic origin. In group C no clear correlation could be found with geographic origin, also arguably due to insufficient sampling. In the Arabian Peninsula many different groups (all except B and D) were found and in Turkey two groups (C and E) were present. Also, two distinct groups (A and F) were found in East Africa.
Evolutionary history of Mccp
Increasing the scale of analysis to tackle the genetic population structure and evolution of the monomorphic pathogen Mccp
Analysing the diversity in genetically monomorphic pathogens such as Mccp is a challenge because of their low variability. Classical molecular approaches based on sequencing a few gene fragments provide so little appreciable diversity that, to date, few questions regarding the epidemiology and evolution of CCPP have been answered. Technological advances in DNA sequencing now allow the rapid sequencing of complete bacterial genomes and offer the possibility to multiplex many different isolates, thereby dramatically increasing the discovery of polymorphic positions to analyse species diversity . In this study, using data from NGS technologies made it possible to develop new approaches based on the analysis of a large number of genes.
The choice of genes for Mccp was based on a recent study of the evolutionary history of Mmm . This pathogen, responsible for CBPP, is also a very monomorphic species for which classical methods like MLST do not afford sufficient discriminatory power. In this previous study, the rigorous selection of 62 genes from the core genome allowed the differentiation of 19 genotypes out of 20 strains analysed and made it possible to elucidate the evolutionary history of CBPP. Therefore, the same approach was chosen to investigate the molecular epidemiology and evolution of CCPP. In the present study, the initial set of 62 genes was reduced to 57 genes after exclusion of missing and duplicated genes in Mccp. The selection finally comprised 47 coding sequences and 10 pseudogenes of the core genome. However, it should be noted that the complete set of 57 genes was only used to conduct genotyping analysis. Keeping the pseudogenes in this gene set gave the genotyping system increased discriminatory power. Indeed, because of their neutral evolution, pseudogenes benefit from a much higher rate of nucleotide substitution than coding sequences . As a result, in our study, the inclusion of pseudogenes made it possible to distinguish strains that could not be distinguished in the phylogenetic analysis based only on 47 coding sequences. Furthermore, almost 20% of the selected Mccp genes were pseudogenes, implying that the analysis of new strains in future genotyping studies may possibly disclose additional pseudogenisation events. On the other hand, for evolutionary analysis, all pseudogenes were discarded. Precise phylogenetic inference and reliable molecular dating were thus based on 47 coding gene sequences presenting similar rates of evolution. Furthermore, detailed sequence analysis did not reveal any mutation hotspots that may have suggested inter-specific recombination. Therefore, the set of genes was considered appropriate for reliable evolutionary analysis.
This multigene approach was applied to a sample of 25 epidemiologically unrelated strains, representing the known global diversity of this species. A large number of stable molecular markers were then selected from this sample by fixing the threshold of polymorphic positions to 85%, therefore excluding potential hypervariable positions. The strategy used here to identify polymorphic positions, based on the comparison of multiple independent strains representing the global species diversity, made it possible to minimise the phylogenetic discovery bias .
In the near future, “Whole genome sequencing” (WGS) is expected to become customary for molecular genotyping, since it provides access to information on the entire genome . However, to enable the use of WGS in routine diagnostic laboratories, the main challenge remains the provision of accessible bioinformatics tools to enable the extraction of relevant information out of the enormous amount of data generated, for a fast and reliable analysis [40,41]. The use of WGS data from diverse strains would also facilitate new approaches to elucidate the evolution of monomorphic species. However, for evolutionary analyses, the problem is further complicated by the need for reliable molecular clock rates .
A new discriminatory system for CCPP epidemiological investigations
The genotyping system developed here improves robustness and resolution for CCPP epidemiological analysis compared to classical typing methods. The MLST method could not be applied due to insufficient discriminatory power. As an illustration, an MLST scheme based on seven partial housekeeping gene sequences used for the study of the evolution of the M. mycoides cluster revealed only six polymorphic positions between 14 Mccp strains . Seven of these strains were included in our study in which they were distinguished by 135 polymorphisms. As an alternative to MLST, an MLSA system had been developed , which remained limited in terms of variability. Our new system, which includes most of the strains previously analysed by MLSA (Table 1), is based on a higher number of polymorphic positions (239 polymorphisms identified in 77 912 bp, versus 53 polymorphisms found in 6747 bp for the MLSA system), resulting in increased discriminatory power (Simpson’s index of diversity of 0.997 (0.988-1.000) versus 0.964 (0.937-0.991) for the MLSA system).
The six genotyping groups identified in this study (groups A to F; Figure 1) are basically in accordance with MLSA clustering (groups 1 to 5; Table 1) except that diversity analyses now provide evidence for an additional group represented by a single Chinese strain. This strain (M1601) was previously included in MLSA group 3, together with two other strains from the Asian continent, although it branched separately in the MLSA tree, rendering this group polyphyletic. This anomaly was previously accepted because strain M1601 was the only one in the group that did not present a large deletion of around 1000 bp, which may have biased distance analysis. However, the present analysis, which is more robust and discriminatory, called for a re-classification of this genotype in a new group (D), the distribution and variability of which remain to be characterised. Most importantly, the intra-group resolution was dramatically increased in the present system, allowing all but two of the strains analysed to be differentiated, while MLSA only afforded further distribution of strains into sub-groups, resulting in a biased presentation of the evolutionary distance between strains (Table 1).
A highly discriminatory genotyping tool is now available for the molecular epidemiology of CCPP, allowing precise epidemiological investigations. Recent works have reported the spread of CCPP to new territories, such as the Indian Ocean  or the Thrace region of Turkey , which constituted a serious threat to Europe. Also, Mccp has been isolated from an increasing variety of wild ungulate species [43–45]. This striking discovery has changed our view of the strict host specificity of Mccp and poses new challenges for CCPP surveillance and control. Tracing new outbreaks requires accurate knowledge of the genetic population structure and the geographic distribution of strains, while the main limitation to such a study remains the paucity of Mccp isolations due to the fastidiousness of this bacterium. However, the use of high throughput sequencing, applied directly on Mccp-rich pleural fluid samples, could circumvent this problem and ease the establishment of a representative Mccp genome database.
The monomorphic bacterium Mccp: a recently emerged pathogen undergoing rapid evolution
A large-scale genomic approach made it possible to infer the first robust and discriminatory phylogeny of Mccp, to estimate its evolutionary rate, and to date its emergence. A previous study on the molecular evolution of Mccp, based on polymorphisms in the two 16S rRNA gene copies , only provided evidence for two major lineages for the species. Our study confirmed the presence of these two lineages but, most importantly, it enabled a dramatic increase in the resolution and robustness of the phylogenetic analysis, resulting in the identification of six well-supported clusters. In addition, this multigene approach made it possible to overcome the bias caused by single gene specificities, providing a more reliable inference of the phylogeny of Mccp .
The mutation rate of Mccp was estimated by BEAST at 1.3 × 10−6 substitutions per site per year, based exclusively on coding gene sequences. This estimate is in agreement with previous reports for other mycoplasma species [20,47], and exceeds the mutation rate estimated for Mmm at 5.0 × 10−7 based on practically the same gene set. Mycoplasmas are among the fastest evolving bacteria, which may be explained by the loss of DNA repair genes during reductive genome evolution . These high mutation rates may facilitate the accumulation of nonsense and frame-shift mutations resulting in gene decay, as evidenced in the analysis of the complete genome of strain 9231-Abomsa, which revealed the presence of 248 pseudogenes representing 25% of the genome . Increased fixation of detrimental mutations by Mccp may be associated with the severe population bottleneck which characterised the emergence of this species. Indeed, evolutionary and dating analysis made it possible to recreate a plausible scenario for the emergence of the species: i.e., the specialisation of Mcc, a goat pathogen presenting a wider tissue tropism, in a restricted ecological niche (the lung) resulted in the emergence of Mccp. This emergence is indeed a recent event, estimated by BEAST analysis at about 270 years ago (between 120 and 736), in agreement with the estimate of Fischer et al. , dating the origin of the species between 56 and 490 years ago. Mccp, like other young pathogens, has evolved too recently to allow the accumulation of many mutations [49,50], thus explaining its low genetic diversity despite a high mutation rate.
Combining large-scale genomic data with spatial and temporal data provides an unprecedented resolution to elaborate plausible scenarios regarding CCPP evolution and epidemiology
BEAST analysis showed that the ancestor of all Mccp strains emerged about 270 years ago, around 1740. The first unequivocal accounts of CCPP were reported more than a century after this predicted emergence in Algeria  and then in South Africa, involving goats imported from Turkey . It is therefore difficult to establish where this emergence took place simply by examining the literature, particularly since, at that time, CCPP was confused with diseases caused by other mycoplasmas of the M. mycoides cluster which frequently infect goats and induce similar clinical signs and lesions. Historical data must therefore be analysed very carefully to ensure that Mccp was indeed the etiologic agent.
In East Africa, two distinct groups (A and F) were found, each belonging to a different lineage (Figure 3), with MRCAs estimated around 35 and 90 years ago, respectively. Intra-group variability was extremely low for group A strains, while Group F strains showed higher genetic diversity, in agreement with longer evolution times allowing further clonal expansion. The presence of two extremely divergent groups in East Africa suggests that Mccp emerged in this region on at least two different occasions. Strains from both groups then spread thanks to trade and other animal movements between neighbouring countries. However, the origin and dynamics of co-evolution of these two distinctive groups in East Africa remain to be elucidated.
In Central Africa, a strict correlation between genotype (group B) and geographic origin was identified. In addition, the four strains belonging to this group showed very little genetic diversity and the date of their MRCA was estimated at around 40 years ago, suggesting that Mccp must have emerged very recently in Central Africa, as stated by Lefevre et al. when the first strain was isolated in Chad in 1987 . In fact, CCPP was not previously suspected by veterinarians established in this region during the colonial times, although they were familiar with the typical clinical picture of this disease, which was described in scientific journals and reference works . If we consider the distribution of other contagious diseases of goats, such as the “peste des petits ruminants”, the presence of CCPP may also be suspected in West Africa. An active search involving improved epidemiological surveillance networks will be required to determine the western limits of CCPP distribution in Africa.
In conclusion, this study illustrates how using a large scale genomic approach can provide an unprecedented resolution to analyse the dynamics and evolution of Mccp. Combining high-throughput genomic data with spatial and temporal data enabled a comprehensive view of the epidemiology of CCPP, which will facilitate the development of improved disease surveillance and control measures. Further efforts are now needed to better define the genetic diversity and distribution of Mccp, particularly in North and Central Africa and in Asia, from where very few strains have been isolated to date.
The authors wish to thank Caroline COSTE for supervision of bioinformatics development, Dr Amirbekov and Dr Ostrowski for providing the sample from which strain 12002 was isolated. This work was supported by the high performance cluster of the UMR AGAP CIRAD.
- OIE, Contagious caprine pleuropneumonia, Chapter 2.7.6. in: OIE Terrestrial Manual [on line] (2014). http://www.oie.int/fileadmin/Home/eng/Health_standards/tahm/2.07.06_CCPP.pdf
- Manso-Silvan L, Vilei EM, Sachse K, Djordjevic SP, Thiaucourt F, Frey J (2009) Mycoplasma leachii sp. nov. as a new species designation for Mycoplasma sp. bovine group 7 of Leach, and reclassification of Mycoplasma mycoides subsp. mycoides LC as a serovar of Mycoplasma mycoides subsp. capri. Int J Syst Evol Microbiol 59:1353–1358PubMedView ArticleGoogle Scholar
- MacOwan KJ, Minette JE (1976) A mycoplasma from acute contagious caprine pleuropneumonia in Kenya. Trop Anim Health Prod 8:91–95PubMedView ArticleGoogle Scholar
- Manso-Silvan 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:86PubMed CentralPubMedView ArticleGoogle Scholar
- Nicholas R, Churchward C (2012) Contagious caprine pleuropneumonia: new aspects of an old disease. Transbound Emerg Dis 59:189–196PubMedView ArticleGoogle Scholar
- Amirbekov M, Murvatulloev S, Ferrari G (2010) Contagious caprine pleuropneumonia detected for the first time in Tajikistan. EMPRES Transboundary Anim Dis Bull 35:20–22Google Scholar
- Awan MA, Abbas F, Yasinzai M, Nicholas RA, Babar S, Ayling RD, Attique MA, Ahmed Z, Wadood A, Khan FA (2010) 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–3406PubMedView ArticleGoogle Scholar
- 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 (in Chinese)PubMedGoogle Scholar
- Srivastava AK, Meenowa D, Barden G, Salguero FJ, Churchward C, Nicholas RA (2010) Contagious caprine pleuropneumonia in Mauritius. Vet Rec 167:304–305PubMedView ArticleGoogle Scholar
- Ozdemir U, Ozdemir E, March JB, Churchward C, Nicholas RA (2005) Contagious caprine pleuropneumonia in the Thrace region of Turkey. Vet Rec 156:286–287PubMedView ArticleGoogle Scholar
- Pettersson B, Bolske G, Thiaucourt F, Uhlen M, Johansson KE (1998) Molecular evolution of Mycoplasma capricolum subsp. capripneumoniae strains, based on polymorphisms in the 16S rRNA genes. J Bacteriol 180:2350–2358PubMed CentralPubMedGoogle Scholar
- Woese CR (1987) Bacterial evolution. Microbiol Rev 51:221–271PubMed CentralPubMedGoogle Scholar
- Pettersson B, Leitner T, Ronaghi M, Bolske G, Uhlen M, Johansson KE (1996) Phylogeny of the Mycoplasma mycoides cluster as determined by sequence analysis of the 16S rRNA genes from the two rRNA operons. J Bacteriol 178:4131–4142PubMed CentralPubMedGoogle Scholar
- Manso-Silvan L, Perrier X, Thiaucourt F (2007) Phylogeny of the Mycoplasma mycoides cluster based on analysis of five conserved protein-coding sequences and possible implications for the taxonomy of the group. Int J Syst Evol Microbiol 57:2247–2258PubMedView ArticleGoogle Scholar
- Maiden MC, Bygraves JA, Feil E, Morelli G, Russell JE, Urwin R, Zhang Q, Zhou J, Zurth K, Caugant DA, Feavers IM, Achtman M, Spratt BG (1998) Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms. Proc Natl Acad Sci U S A 95:3140–3145PubMed CentralPubMedView ArticleGoogle Scholar
- Lorenzon S, Wesonga H, Ygesu L, Tekleghiorgis T, Maikano Y, Angaya M, Hendrikx P, Thiaucourt P (2002) Evolution of M.capricolum subsp. capripneumoniae strains and molecular epidemiology of contagious caprine pleuropneumonia. Vet Microbiol 85:111–123PubMedView ArticleGoogle Scholar
- Achtman M (2008) Evolution, population structure, and phylogeography of genetically monomorphic bacterial pathogens. Annu Rev Microbiol 62:53–70PubMedView ArticleGoogle Scholar
- Roumagnac P, Weill FX, Dolecek C, Baker S, Brisse S, Chinh NT, Le TA, Acosta CJ, Farrar J, Dougan G, Achtman M (2006) Evolutionary history of Salmonella typhi. Science 314:1301–1304PubMed CentralPubMedView ArticleGoogle Scholar
- Morelli G, Song Y, Mazzoni CJ, Eppinger M, Roumagnac P, Wagner DM, Feldkamp M, Kusecek B, Vogler AJ, Li Y, Cui Y, Thomson NR, Jombart T, Leblois R, Lichtner P, Rahalison L, Petersen JM, Balloux F, Keim P, Wirth T, Ravel J, Yang R, Carniel E, Achtman M (2010) Yersinia pestis genome sequencing identifies patterns of global phylogenetic diversity. Nat Genet 42:1140–1143PubMed CentralPubMedView ArticleGoogle Scholar
- Dupuy V, Manso-Silvan L, Barbe V, Thebault P, Dordet-Frisoni E, Citti C, Poumarat F, Blanchard A, Breton M, Sirand-Pugnet P, Thiaucourt F (2012) Evolutionary history of contagious bovine pleuropneumonia using next generation sequencing of Mycoplasma mycoides Subsp. mycoides “Small Colony”. PLoS One 7:e46821PubMed CentralPubMedView ArticleGoogle Scholar
- Roetzer A, Diel R, Kohl TA, Ruckert C, Nubel U, Blom J, Wirth T, Jaenicke S, Schuback S, Rusch-Gerdes S, Supply P, Kalinowski J, Niemann S (2013) Whole genome sequencing versus traditional genotyping for investigation of a Mycobacterium tuberculosis outbreak: a longitudinal molecular epidemiological study. PLoS Med 10:e1001387PubMed CentralPubMedView ArticleGoogle Scholar
- Chu Y, Gao P, Zhao P, He Y, Liao N, Jackman S, Zhao Y, Birol I, Duan X, Lu Z (2011) Genome sequence of Mycoplasma capricolum subsp. capripneumoniae strain M1601. J Bacteriol 193:6098–6099PubMed CentralPubMedView ArticleGoogle Scholar
- Dupuy V, Thiaucourt F (2014) Complete Genome Sequence of Mycoplasma capricolum subsp. capripneumoniae Strain 9231-Abomsa. Genome Announc 2:e01067–14PubMed CentralPubMedView ArticleGoogle Scholar
- 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–14PubMed CentralPubMedView ArticleGoogle Scholar
- Lorenzon S, Manso-Silvan L, Thiaucourt F (2008) Specific real-time PCR assays for the detection and quantification of Mycoplasma mycoides subsp. mycoides SC and Mycoplasma capricolum subsp. capripneumoniae. Molecul Cell Probes 22:324–328View ArticleGoogle Scholar
- Sambrook J, Fritsch EF, Maniatis T (1989) Molecular Cloning: A Laboratory Manual. Cold Spring Harbor, New YorkGoogle Scholar
- Galaxy - Data intensive biology for everyone [on line] http://galaxyproject.org/
- Librado P, Rozas J (2009) DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25:1451–1452PubMedView ArticleGoogle Scholar
- Bandelt HJ, Forster P, Rohl A (1999) Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol 16:37–48PubMedView ArticleGoogle Scholar
- 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–2466PubMed CentralPubMedGoogle Scholar
- 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–4192PubMed CentralPubMedView ArticleGoogle Scholar
- Posada D, Crandall KA (1998) MODELTEST: testing the model of DNA substitution. Bioinformatics 14:817–818PubMedView ArticleGoogle Scholar
- Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O (2010) New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol 59:307–321PubMedView ArticleGoogle Scholar
- Drummond AJ, Rambaut A (2007) BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol 7:214PubMed CentralPubMedView ArticleGoogle Scholar
- Kass RE, Raftery AE (1995) Bayes Factors. J Am Stat Assoc 90:773–795View ArticleGoogle Scholar
- Suchard MA, Weiss RE, Sinsheimer JS (2001) Bayesian selection of continuous-time Markov chain evolutionary models. Mol Biol Evol 18:1001–1013PubMedView ArticleGoogle Scholar
- Baker S, Hanage WP, Holt KE (2010) Navigating the future of bacterial molecular epidemiology. Curr Opin Microbiol 13:640–645PubMed CentralPubMedView ArticleGoogle Scholar
- Li WH, Gojobori T, Nei M (1981) Pseudogenes as a paradigm of neutral evolution. Nature 292:237–239PubMedView ArticleGoogle Scholar
- Eyre DW, Golubchik T, Gordon NC, Bowden R, Piazza P, Batty EM, Ip CL, Wilson DJ, Didelot X, O’Connor L, Lay R, Buck D, Kearns AM, Shaw A, Paul J, Wilcox MH, Donnelly PJ, Peto TE, Walker AS, Crook DW (2012) A pilot study of rapid benchtop sequencing of Staphylococcus aureus and Clostridium difficile for outbreak detection and surveillance. BMJ Open 2:e001124PubMed CentralPubMedView ArticleGoogle Scholar
- Carrico JA, Sabat AJ, Friedrich AW, Ramirez M, ESCMID Study Group for Epidemiological Markers (ESGEM) (2013) Bioinformatics in bacterial molecular epidemiology and public health: databases, tools and the next-generation sequencing revolution. Euro Surveill 18:20382PubMedGoogle Scholar
- Sabat AJ, Budimir A, Nashev D, Sa-Leao R, van Dijl J, Laurent F, Grundmann H, Friedrich AW, ESCMID Study Group for Epidemiological Markers (ESGEM) (2013) Overview of molecular typing methods for outbreak detection and epidemiological surveillance. Euro Surveill 18:20380PubMedGoogle Scholar
- Fischer A, Shapiro B, Muriuki C, Heller M, Schnee C, Bongcam-Rudloff E, Vilei EM, Frey J, Jores J (2012) The origin of the ‘Mycoplasma mycoides cluster’ coincides with domestication of ruminants. PLoS One 7:e36150PubMed CentralPubMedView ArticleGoogle Scholar
- 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 Wild Med 38:93–96View ArticleGoogle Scholar
- Chaber A, Lignereux L, Qassimi MA, Saegerman C, Manso-Silvan L, Dupuy V, Thiaucourt F (2014) Fatal transmission of contagious caprine pleuropneumonia to an Arabian oryx (Oryx leucoryx). Vet Microbiol 173:156–159PubMedView ArticleGoogle Scholar
- 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–2053PubMed CentralPubMedView ArticleGoogle Scholar
- Coenye T, Gevers D, Van de Peer Y, Vandamme P, Swings J (2005) Towards a prokaryotic genomic taxonomy. FEMS Microbiol Rev 29:147–167PubMedView ArticleGoogle Scholar
- Delaney NF, Balenger S, Bonneaud C, Marx CJ, Hill GE, Ferguson-Noel N, Tsai P, Rodrigo A, Edwards SV (2012) Ultrafast Evolution and Loss of CRISPRs Following a Host Shift in a Novel Wildlife Pathogen, Mycoplasma gallisepticum. PLoS Genet 8:e100251View ArticleGoogle Scholar
- Carvalho FM, Fonseca MM, Batistuzzo De Medeiros S, Scortecci KC, Blaha CA, Agnez-Lima LF (2005) DNA repair in reduced genome: the Mycoplasma model. Gene 360:111–119PubMedView ArticleGoogle Scholar
- Kidgell C, Reichard U, Wain J, Linz B, Torpdahl M, Dougan G, Achtman M (2002) Salmonella typhi, the causative agent of typhoid fever, is approximately 50,000 years old. Infect Genet Evol 2:39–45PubMedView ArticleGoogle Scholar
- Sreevatsan S, Pan X, Stockbauer KE, Connell ND, Kreiswirth BN, Whittam TS, Musser JM (1997) Restricted structural gene polymorphism in the Mycobacterium tuberculosis complex indicates evolutionarily recent global dissemination. Proc Natl Acad Sci U S A 94:9869–9874PubMed CentralPubMedView ArticleGoogle Scholar
- Thomas P. (1873) Rapport médical sur le Bou Frida. In: l’Algérie G.g.c.d. (Ed.), Jourdan A, Alger, pp. 35Google Scholar
- Hutcheon D (1881) Contagious pleuro-pneumonia in angora goats. Vet J 13:171–180Google Scholar
- Walker GK (1914) Pleuro-pneumonia of goats in the Kangra district, Punjab, India. J Comp Pathol Ther 27:68–71View ArticleGoogle Scholar
- Lefevre PC, Breard A, Alfarouk I, Buron S (1987) Mycoplasma species F 38 isolated in Chad. Vet Rec 121:575–576PubMedGoogle Scholar
- Curasson G (1936) Maladies à ultravirus. In: Curasson G (ed) Traité de pathologie exotique vétérinaire et comparée. Vigot frères, Paris, pp 590–606Google Scholar
- Jones GE, Wood AR (1988) Microbiological and serological studies on caprine pneumonias in Oman. Res Vet Sci 44:125–131PubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.