BackgroundPlague was introduced to Madagascar in 1898 and continues to be a significant human health problem. It exists mainly in the central highlands, but in the 1990s was reintroduced to the port city of Mahajanga, where it caused extensive human outbreaks. Despite its prevalence, the phylogeography and molecular epidemiology of Y. pestis in Madagascar has been difficult to study due to the great genetic similarity among isolates. We examine island-wide geographic-genetic patterns based upon whole-genome discovery of SNPs, SNP genotyping and hypervariable variable-number tandem repeat (VNTR) loci to gain insight into the maintenance and spread of Y. pestis in Madagascar.Methodology/Principal FindingsWe analyzed a set of 262 Malagasy isolates using a set of 56 SNPs and a 43-locus multi-locus VNTR analysis (MLVA) system. We then analyzed the geographic distribution of the subclades and identified patterns related to the maintenance and spread of plague in Madagascar. We find relatively high levels of VNTR diversity in addition to several SNP differences. We identify two major groups, Groups I and II, which are subsequently divided into 11 and 4 subclades, respectively. Y. pestis appears to be maintained in several geographically separate subpopulations. There is also evidence for multiple long distance transfers of Y. pestis, likely human mediated. Such transfers have resulted in the reintroduction and establishment of plague in the port city of Mahajanga, where there is evidence for multiple transfers both from and to the central highlands.Conclusions/SignificanceThe maintenance and spread of Y. pestis in Madagascar is a dynamic and highly active process that relies on the natural cycle between the primary host, the black rat, and its flea vectors as well as human activity.

A cluster of human plague cases occurred in the seaport city of Mahajanga, Madagascar, from 1991 to 1999 following 62 years with no evidence of plague, which offered insights into plague pathogen dynamics in an urban environment. We analyzed a set of 44 Mahajanga isolates from this 9-year outbreak, as well as an additional 218 Malagasy isolates from the highland foci. We sequenced the genomes of four Mahajanga strains, performed whole-genome sequence single-nucleotide polymorphism (SNP) discovery on those strains, screened the discovered SNPs, and performed a high-resolution 43-locus multilocus variable-number tandem-repeat analysis of the isolate panel. Twenty-two new SNPs were identified and defined a new phylogenetic lineage among the Malagasy isolates. Phylogeographic analysis suggests that the Mahajanga lineage likely originated in the Ambositra district in the highlands, spread throughout the northern central highlands, and was then introduced into and became transiently established in Mahajanga. Although multiple transfers between the central highlands and Mahajanga occurred, there was a locally differentiating and dominant subpopulation that was primarily responsible for the 1991-to-1999 Mahajanga outbreaks. Phylotemporal analysis of this Mahajanga subpopulation revealed a cycling pattern of diversity generation and loss that occurred during and after each outbreak. This pattern is consistent with severe interseasonal genetic bottlenecks along with large seasonal population expansions. The ultimate extinction of plague pathogens in Mahajanga suggests that, in this environment, the plague pathogen niche is tenuous at best. However, the temporary large pathogen population expansion provides the means for plague pathogens to disperse and become ecologically established in more suitable nonurban environments.

No abstract