# Low cost, low tech SNP genotyping tools for resource-limited areas: Plague in Madagascar as a model

Cedar L Mitchell, Voahangy Andrianaivoarimanana, Rebecca E Colman,

Joseph Busch, Heidie Hornstra-O'Neill, Paul Keim, David Wagner,

Minoarisoa Rajerison, Dawn N Birdsell

► To cite this version:

Cedar L Mitchell, Voahangy Andrianaivoarimanana, Rebecca E Colman, Joseph Busch, Heidie Hornstra-O'Neill, et al.. Low cost, low tech SNP genotyping tools for resource-limited areas: Plague in Madagascar as a model. PLoS Neglected Tropical Diseases, 2017, 11 (12), pp.e0006077. 10.1371/journal.pntd.0006077. pasteur-01930482

HAL Id: pasteur-01930482

https://riip.hal.science/pasteur-01930482

Submitted on 22 Nov 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Low cost, low tech SNP genotyping tools for resource-limited areas: Plague in Madagascar as a model

Cedar L. Mitchell¹ᵃ, Voahangy Andrianaivoarimanana², Rebecca E. Colman¹ᵇ, Joseph Busch¹, Heidie Hornstra-O'Neill¹, Paul S. Keim¹,³, David M. Wagner¹, Minoarisoa Rajerison²*, Dawn N. Birdsell¹*

¹ The Pathogen and Microbiome Institute, Northern Arizona University, Flagstaff, Arizona, United States of America, ² Plague Unit, Institut Pasteur de Madagascar, Antananarivo, Madagascar, ³ Translational Genomics Research Institute, Flagstaff, Arizona, United States of America

ᵃ Current address: Department of Epidemiology, Gillings School of Global Public Health, University of North Carolina, Chapel Hill, North Carolina, United States of America

ᵇ Current address: Department of Medicine, University of California, San Diego, La Jolla, California, United States of America

* mino@pasteur.mg (MR); Dawn.Birdsell@nau.edu (DNB)

<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

## OPEN ACCESS

**Citation:** Mitchell CL, Andrianaivoarimanana V, Colman RE, Busch J, Hornstra-O'Neill H, Keim PS, et al. (2017) Low cost, low tech SNP genotyping tools for resource-limited areas: Plague in Madagascar as a model. PLoS Negl Trop Dis 11 (12): e0006077. https://doi.org/10.1371/journal.pntd.0006077

**Editor:** Cheryl Andam, University of New Hampshire, UNITED STATES

**Received:** April 4, 2017

**Accepted:** October 27, 2017

**Published:** December 11, 2017

**Copyright:** © 2017 Mitchell et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This work was supported by financial aid from Experiment.com, and the Center for International Education, the University College, and the Hooper Undergraduate Research Award program at Northern Arizona University. Cedar Mitchell was employed by Northern Arizona University at the time of this work. CM currently is

## Abstract

### Background

Genetic analysis of pathogenic organisms is a useful tool for linking human cases together and/or to potential environmental sources. The resulting data can also provide information on evolutionary patterns within a targeted species and phenotypic traits. However, the instruments often used to generate genotyping data, such as single nucleotide polymorphisms (SNPs), can be expensive and sometimes require advanced technologies to implement. This places many genotyping tools out of reach for laboratories that do not specialize in genetic studies and/or lack the requisite financial and technological resources. To address this issue, we developed a low cost and low tech genotyping system, termed agarose-MAMA, which combines traditional PCR and agarose gel electrophoresis to target phylogenetically informative SNPs.

### Methodology/Principal findings

To demonstrate the utility of this approach for generating genotype data in a resource-constrained area (Madagascar), we designed an agarose-MAMA system targeting previously characterized SNPs within *Yersinia pestis*, the causative agent of plague. We then used this system to genetically type pathogenic strains of *Y. pestis* in a Malagasy laboratory not specialized in genetic studies, the Institut Pasteur de Madagascar (IPM). We conducted rigorous assay performance validations to assess potential variation introduced by differing research facilities, reagents, and personnel and found no difference in SNP genotyping results. These agarose-MAMA PCR assays are currently employed as an investigative tool at IPM, providing Malagasy researchers a means to improve the value of their plague epidemiological investigations by linking outbreaks to potential sources through genetic characterization of isolates and to improve understanding of disease ecology that may contribute to a long-term control effort.
a Ph.D. student in the department of epidemiology,
Gillings School of Global Public Health, University
of North Carolina Chapel Hill. The funders had no
role in study design, data collection and analysis,
decision to publish, or preparation of the
manuscript.

**Competing interests:** The authors have declared
that no competing interests exist.

Conclusions

The success of our study demonstrates that the SNP-based genotyping capacity of labora-
tories in developing countries can be expanded with manageable financial cost for resource
constraint laboratories. This is a practical formula that reduces resource-driven limitations to
genetic research and promises to advance global collective knowledge of infectious dis-
eases emanating from resource limited regions of the world.

Author summary

Although genetic characterization of pathogenic organisms is a powerful tool for investi-
gating outbreak origins and transmission, associated high upfront costs and demanding
technological maintenance exclude this tool for many under-resourced laboratories. Para-
doxically, resource constrained regions commonly suffer from high rates of infectious
diseases and could benefit most from genetic tracking tools. One such country is Mada-
gascar, which lacks resources to acquire high tech genetic typing equipment, yet suffers
from seasonal human plague outbreaks. A serious disease, plague is caused by the clonal
bacterium, *Yersinia pestis*, and is capable of causing human outbreaks. Using plague as a
model organism, we developed a genetic typing method that requires only basic, widely
used molecular machinery. Our tools target unique single mutations in the *Y. pestis*
genome to assign isolates to distinct phylogenetic groups with known geographical distri-
butions. Transfer of this technology to Madagascar permits genetic characterization of
strains from current outbreaks. This eliminates the need for external genetic analysis and
expands the research capacity of this resource-constrained laboratory by allowing rapid,
in-house strain typing. Ultimately, our goal is to help improve the ability of local institutes
to genetically characterize circulating strains, link outbreaks to originating sources, and
improve our understanding of the ecology of tropical diseases in resource-limited regions
of the world.

**Introduction**

Single nucleotide polymorphisms (SNPs) are highly valuable genetic markers that have advanced our knowledge of diverse biological fields such as human health [1,2], infectious disease epidemiology [3-5], agriculture [6], and ecology [7], among others. In the study of infectious diseases SNPs can be informative of bacterial phenotype, such as antibiotic susceptibility [8,9], and also can be used to classify unknown strains. For non-recombining bacterial pathogens, most of their SNPs become fixed in the genome and are faithfully replicated throughout future generations [3,10]. These stable signatures can be used to classify unknown strains into known phylogenetic groups according to SNP profiles [3,11,12]. Within the context of epidemiological investigations, these SNP profiles can link isolates from active outbreak sites to possible sources and help track disease transmission patterns [5,13,14].

Genotyping assays that use real-time PCR to identify single SNPs remain in demand despite
the wide-scale availability of whole genome sequence (WGS) data and continued reductions in
WGS costs. For many research facilities that are interested in small-scale studies or face
resource limitations, a WGS-based approach to SNP typing is not a feasible nor a desirable
option. A variety of other technological platforms have been employed for SNP typing and
have been extensively described in several publications [15–18]. Popular platforms for SNP typing use real-time PCR instruments that employ Dual Probe TaqMan assays or melt-MAMA SNP assays [19–21]. But real-time platforms are not commonly available in resource constrained laboratories, due to their high upfront costs and the need for ongoing highly technical instrument maintenance. However, a more simplified method for SNP genotyping that employs conventional PCR coupled with standard agarose gel electrophoresis (agarose-MAMA) is a viable alternative in these settings. The advantage of this alternative method is that it utilizes relatively inexpensive instruments that are almost universally available even in developing nations where it is used for a variety of molecular applications. Much of this is due to the simplicity of the agarose gel electrophoresis platform, in contrast to the complex instrumentation of the real-time platform [20,22].

To illustrate the effectiveness of agarose-MAMA as a SNP genotyping tool in resource constrained laboratories, we developed *Y. pestis* assays for use at the Institut Pasteur de Madagascar (IPM). *Y. pestis* is the bacterium infamously known as the causal agent of the disease plague. *Y. pestis* is ecologically established on nearly every inhabited continent [12,23] and remains a particularly significant threat to human health in developing nations in Africa and especially the island country of Madagascar [24,25]. Primarily a zoonotic agent, *Y. pestis* has a complex ecological cycle involving rodent-host populations and flea vectors. In unfortunate circumstances, humans are incidental hosts [26]. Without prompt antibiotic treatment, human death rates vary from 30–60% to 100% depending on the route of exposure, with pneumonic plague being the most deadly [24,25]. Within the last two decades, Madagascar has reported some of the highest incidences of human plague infections throughout the world [24]. Given that *Y. pestis* is a re-emerging public health threat in Madagascar [27–29] and other developing African nations, performing SNP typing studies through the use of agarose-MAMA tools could lead to greatly improved epidemiological investigations and to more effective disease management.

Here, we describe how we resolved technological limitations that prevented *Y. pestis* SNP genotyping studies at IPM by: 1) re-designing a SNP-based *Y. pestis* genotyping system (from melt-MAMA real-time platform to agarose-MAMA) to be compatible with existing resources in Madagascar, and 2) validating these genetic tools at IPM for adoption. Collectively, these achievements have removed a research barrier at IPM once imposed by technological disparities and serve as a promising model for building similar research capacities in other developing nations affected by plague and/or other pathogens common to low-resource settings. In addition, these tools generate data that can be compared to other sequence-based methods and can feed into existing global databases. This ability can strengthen scientific exchange between affected countries and the international scientific community, thereby advancing the global collective knowledge of dangerous infectious diseases.

## Results and discussion

### Removal of technological barriers: Agarose MAMA—a highly accurate genetic typing tool

We present a SNP genotyping tool that uses standard genetic equipment and reagents that are accessible to nearly any laboratory. The successful development of agarose-MAMA tools removes the dependence on real-time instruments for conducting SNP-based epidemiological studies. Due to emerging infectious diseases being most globally prevalent in resource challenged countries [30,31], building research capacities in these regions has been a goal for the World Health Organization and other agencies charged with global biosafety, response, and biosecurity efforts [32,33]. It is now possible for researchers in resource constrained locations
<figure><img src="image_4.png" /><figcaption>**Fig 1.** Simplified SNP phylogeny of *Y. pestis* depicting subgroups identified in Madagascar. Colored circles indicate phylogenetic groups as previously described [35][11,34]. Group names are assigned as letters and sometimes followed by a number (e.g. s4). Black bars indicate the phylogenetic positions of the 18 SNPs targeted for MAMA PCR design in this study.</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0006077.g001

to conduct SNP studies to obtain bacterial phenotype information such as antibiotic resistance or phylogenetic information to more robustly understand the dynamics of disease transmission and potentially identify sources of human infections.

The phylogeny of *Y. pestis* in Madagascar has recently expanded [34] to include more phylogenetic groups (Fig 1) and the 18 assays developed for this study define a subset of these phylogenetic groups or lineages as illustrated on a simplified phylogeny (Fig 1). The genotyping results generated from each agarose-MAMA was identical to independent SNP genotyping technologies [34] when tested across the same diverse panel of 16 *Y. pestis* DNA strains. Assay specificity remained intact when these agarose-MAMAs were tested on diverse types of negative controls (high levels of human, *Leptospira* spp., and *Bacillus anthracis* DNA, or a no template water control). When the assays were employed in a stepwise, hierarchical order (Fig 2), each isolate could be assigned to a single lineage or one subgroup as depicted in the simplified Malagasy *Y. pestis* phylogeny (Fig 1). The equivalent performance of our agarose-MAMA tools to melt-MAMAs [11] and other independent SNP technologies [34] demonstrate that MAMA tools preserve their genotyping accuracy independent of the technology platform used.
<figure><img src="image_5.png" /><figcaption>**Fig 2.** Flow chart providing the hierarchical organization of 18 SNP-genotyping assays. The hierarchy provides the sequential order to allow for stepwise identification of the genotype of an unknown isolate. Genetic subgroups are represented by colored circles labeled with each genetic subgroup. Individual assays are represented by the black bar labeled with each assay ID. The derived and ancestral allele states for each assay are indicated immediately above and below the black bar.</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0006077.g002

## Optimization of assays to obtain accurate genotypes

For the 18 selected SNPs, we were able to successfully design agarose-MAMAs to achieve 100% genotyping accuracy through a few key optimization steps. Only one of our assays failed initial optimization efforts; however, when we redesigned the assay to target the reverse complement of the reference template DNA (strain CO92), the assay was rescued to full functionality. Essential validation steps included the identification of the optimal ratios for the concentration of the two forward primers (ancestral:derived) as described [20], ideal annealing temperature for each assay, appropriate number of PCR cycles, and an occasional need for MgCl₂ concentration alteration in the reaction mix. These optimization strategies worked well
<figure><img src="image_6.png" /><figcaption>**Fig 3. Size differences in allele-specific PCR products.** A) Schematic of PCR amplicons originating from ancestral and derived genomic templates. For both PCR amplicons, gray represents sequence originating from the primers (forward and reverse) and synthesized internal sequence from PCR extension. The primer sequence and synthesized internal sequence within the amplicon are indicated by the bracket. For the derived amplicon, the blue represents the incorporated GC-clamp which originates from the 5' end of the derived MAMA forward primer. The SNP region for each amplicon is represented in green (ancestral) and yellow (derived) and the deliberate antepenultimate mutation is represented as red (ancestral) and light blue (derived). B) Allele-specific PCR products migrate at different rates on a 2% agarose gel due to their size difference as conferred by the GC-clamp.</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0006077.g003

on 8 other bacterial pathogens [20] and should be applicable to many other pathogenic organisms as well.

The GC-clamp (21-28 oligo base), added to derived MAMA forward primers but excluded from ancestral MAMA forward primers, was sufficient to provide visible size differences between the two allelic states when visualized on an agarose gel (Fig 3). A repeated pattern of 5'-cggttcgggttcgggttcgggttcggg-3' did not appear to induce non-specific binding nor interact with template DNA in an inhibitory manner. However, we observed that the GC-clamp on the derived MAMA primer did frequently confer a competitive advantage over the ancestral MAMA, as previously described for MAMA tools [20]. This competitive advantage resulted in cross-hybridization of primers for a subset of assays. This phenomenon may be based on the ability for GC rich DNA regions to anneal at lower temperatures compared to GC poor regions. The GC-clamp of the derived MAMA primer likely anneals to its target amplicon at lower temperatures compared to the no-clamp ancestral primer. As a consequence, the GC-clamp primer anneals to the template at an earlier time point per PCR cycle than the no-clamp primer, resulting in a competitive kinetic advantage. To correct for this competitive advantage, we followed published guidelines [20] by increasing the ancestral primer concentration relative to the derived primer concentration of the affected assays. In these skewed reaction mixes, the ancestral primer had >1x-4x concentration compared to its paired derived primer. We have detected best results with 4:1 and 2:1 (ancestral:derived) primer ratios, depending on the assay.

We obtained maximal assay specificity by a combination of customizing the annealing temperature per assay and/or adjusting the MgCl₂ concentrations. Most of our assays accurately genotyped SNP alleles at annealing temperatures around 60°C; however, some assays had non-specific product at this condition. For these assays, non-specific PCR products were present in our negative controls and/or in our sample PCR product, visualized as extra banding on
a gel. A reduction of MgCl₂ concentration by 0.5 mM was sufficient to greatly increase the specificity of most assays, as evidenced by the elimination of non-specific amplification. But this reduction in MgCl₂ decreased the PCR robustness in a subset of assays. To address this loss, the annealing temperatures were reduced in affected assays while still preserving SNP specificity. For the few assays that did not respond to lowered MgCl₂ strategy, the raising of annealing temperatures while still maintaining the normal 2.0 mM concentration of MgCl₂ yielded improvements in reducing non-specific products.

## Assay function on clinical samples

The assays described here are capable of genotyping directly from complex clinical samples if pathogen DNA levels are sufficient (Fig 4). This capability is highly important as nearly 43.32% of 775 F1 RDT [36] positive human plague cases between 2011–2015 in Madagascar and many environmental samples (rodents and fleas) do not yield live isolate culture (IPM records). Using two agarose MAMA tools (Mad-05 and Mad-43) to demonstrate proof of principal, we were able to determine that the three complex clinical samples belonged to group I lineage and not group II. This genetic assignment for these samples matched the results of a recent publication [34]. These samples were not further tested on additional agarose-MAMAs. Out of five complex clinical samples positive for the high copy plasmid *pla* gene [37] only four were positive for *Y. pestis* chromosomal DNA, assessed by a TaqMan assay targeting the 3a

<figure><img src="image_7.png" /><figcaption>**Fig 4.** Agarose-MAMA is capable of genotyping *Y. pestis* directly from complex clinical samples if pathogen targets are at sufficient levels. (A) Agarose-MAMA (Mad-43) gel showing the PCR products for two template controls at the expected size for each respective ancestral and derived allele state (#1 and #2, respectively). Three complex clinical samples (#3-#5) also yielded PCR products with the size expected of an ancestral genotype. Two other complex clinical samples (#6 & #7) showed no PCR products but displayed a banding pattern consistent with NTC negative controls. (B) To assess the relative quantity of the *Y. pestis* target in the five clinical samples (#3, #4, #5, #6, #7), we generated amplification plots of these clinical samples on a TaqMan 3a assay used to target *Y. pestis* chromosomal DNA. Three clinical samples (#3-#5) showed amplified at a mid-range cycle-time (Ct) value consistent with high copy numbers of template DNA [21]. These same samples showed a robust signal of a PCR product on the MAMA gel. The two samples (#6 & #7) that failed on the MAMA gel showed a high Ct value and failed amplification, respectively, when tested by real-time PCR. The pairing of TaqMan 3a assay with MAMA gel results on the same templates provided insight to the genotyping capability of MAMA tools on complex clinical samples with low-level target template.</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0006077.g004
gene (Fig 4 and S3 Appendix). Of the four 3a-positive samples, one isolate (Yp3182) gave a late amplification with the 3a assay resulting in a Ct value of 36 using real-time PCR. This indicates a very low concentration amount for this clinical extract, near a single copy of pathogen DNA [20,21]. The other three samples (Yp2483, Yp2486, Yp2485) (Fig 4) generated Ct values ranging from 24–27, which indicates higher concentrations of *Y. pestis* chromosomal DNA in these clinical samples. These more concentrated *Y. pestis* clinical samples were successfully genotyped by agarose-MAMA tools following an increase in the number of PCR cycles. The low-level sample (Yp3182) failed to amplify PCR product using agarose-MAMA tools. Together, these results indicate that agarose-MAMAs can genotype directly from complex clinical samples if the pathogen target (*Y. pestis* DNA) is of sufficient concentration. Published studies show that TaqMan assays are highly sensitive [20,21] and can readily detect minute amount of template not detectable using melt-MAMAs [20]. Our results suggest that the same is true for agarose-MAMAs (Fig 4).

## Transfer of technology—the minor effects of different research facilities and reagents

Following the transfer of our agarose-MAMA tools to the IPM facility, a subset of our agarose-MAMAs was validated for genotyping accuracy. The differing instruments and reagents used at IPM introduced very little variability on assay performance when tested on the same DNA panel previously used at NAU (Fig 5). The optimized PCR conditions identified at NAU for the assays were suitable for most of the assays when conducted at IPM.

<figure><img src="image_8.png" /><figcaption>**Fig 5.** Mad-05 MAMA tool showed congruent genotyping results when generated across two research institutes using identical DNA samples. Side by side comparison of Agarose-MAMA gel images generated at NAU (top) and IPM (bottom) using different reagents and instruments. *Y. pestis* sample names follow the designation scheme practiced at IPM. Individual samples are designated an ID # according to consecutive order of collection per a given year. Strain 87/11 was the 87th isolate collected in 2011.</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0006077.g005
The most apparent difference in assay performance between the two laboratories was an increase in non-specific amplification at high fragment size (Fig 5) and, for some assays, also in our negative controls at IPM, which was not observed at NAU. These non-specific fragments did not affect the MAMA's capability to accurately genotype isolates. Band profiles were not distinctly in line with either the ancestral or derived product but rather appeared either as a smear across both profiles or as a much fainter band falling between the ancestral and derived fragment sizes. We suspect that the basis of this performance difference is due to different Taq polymerases used between the two institutions. NAU employed an antibody-immobilized Taq polymerase (Invitrogen, Carlsbad, CA), which has been characterized as having no polymerase activity prior to a hot-start step in the thermal cycle protocol [38]. However, IPM utilizes regular Taq polymerase, which may begin product synthesis with primers and template in the master mix prior to PCR thermal cycling [38]. Although master mix preparation was done on ice to suppress premature Taq polymerase activity, additional bands above the target PCR product and faint bands in the NTC amplification suggest that pre-PCR Taq activity was not completely suppressed (Fig 5). We therefore reduced the concentration of MgCl₂ in five of our assays (Mad-43, Mad-46, Mad-12, Mad-36 and Mad-58) to 1.5mM and observed the elimination of amplification in the negative controls in most of the affected assays. Reduction of MgCl₂ also imparted a general improvement of assay specificity for positive control and test isolates. Surprisingly, the reduction of MgCl₂ did not necessitate a corresponding increase in the number of PCR cycles for most assays, contrary to the results observed in our laboratory at NAU. Once again, this may be a result of differences in Taq polymerase activity between the two institutions although we have not found any published evidence of this occurring elsewhere. As a way to rule out possible contamination as a source of non-specific banding, all surfaces used for the preparation of PCR reaction mixtures were sterilized with UV light for 15 minutes and decontaminated with 70% ethanol. Following slight modifications, we found that results of the majority of our assays at IPM aligned very well with results produced at NAU (Fig 5 and S5 Appendix).

## Advancing global collective knowledge by building research capacity

Our success in developing agarose-MAMA tools and transferring them to IPM facility in Madagascar demonstrates that this SNP genotyping strategy can be achieved with existing technologies routinely used in developing nations (Fig 5). The transferred agarose-MAMA technology is now in regular use at IPM and was recently used to infer the source of a 2015 pneumonic outbreak [39]. There is ample evidence that this same assay design strategy would transfer well to many other pathogenic organisms [20]. This would allow institutions in developing countries to perform molecular studies in-house and on local infectious organisms. The current methods largely used for short-term control efforts are presence/absence assays (PCR based [40] or protein based [36]) that diagnose the causative agents of disease outbreaks but provide no genetic resolution. Having in-house capabilities to genetically discriminate goes beyond what presence/absence assays can provide, therefore, resource-constrained laboratories will be able to advance their epidemiological capabilities over the status quo. The SNP data they generate locally can be shared among research institutes and compared to existing global databases. These advanced capabilities would accelerate the understanding of plague ecology, persistence, and evolution; which in turn could beneficially inform strategies for disease control.

Building research capacity in low-resource settings is an important endeavor for global bio-safety, response, and biosecurity preparedness [32,33]. Our study is a successful model of achieving this goal pragmatically. We successfully developed and transferred genetic tools to a
developing nation. The design principles for this technology, which we detail above, can be applied to diverse pathogen species [20]. The use of MAMA technology will give the scientific community the means to gain insight into the genetic patterns and population structure of many neglected diseases. This is a practical formula that will advance global collective knowledge of infectious diseases emanating from more impoverished regions of the world.

# Materials and methods

## Ethics statement

The 16 archival strains and five clinical samples were not subject to IRB regulations because they did not meet the federal definition of human subjects research according to 45 CFR 46.102 (f). All samples underwent de-identification of patient information prior to Northern Arizona University (NAU) transfer. They were collected as part of the medical workup mandated by the Ministry of Health in Madagascar and not for the purpose of this study. For this reason the strains used in this study does not meet the federal definition of human subjects research according to 45 CFR 46.102 (f) and therefore are not subject to review from NAU Institutional Review Board.

## DNA

DNA samples utilized in this study were obtained from 16 *Y. pestis* isolates and 5 clinical samples (bubo aspirates and sputum) collected from suspected human plague cases (Table 1). The 16 archival strains and 5 clinical samples (Table 1) originated from diverse geographic locations in Madagascar and were collected as described in a recently published study [34]. Our assays worked well on DNA concentrations that ranged from 1 ng to100 pg. Molecular confirmation of *Y. pestis* in five clinical samples was based on PCR detection of the *pla* gene located in a high copy number plasmid PCP1 in addition to positive F1 RDT [28,36,37,41,42].

## SNP groups

Eighteen previously published SNPs [11,34,35] specific, or canonical [3], for a subset of distinct phylogenetic groups within the *Y. pestis* Malagasy phylogeny (Fig 1) were selected as the targets of agarose MAMA genotyping assays following published guidelines [20]. These selected SNPs can be used in a hierarchical way (Fig 2; Table 2) to assign an unknown strain to one of the most common lineages or phylogenetic subgroups in Madagascar [11,34,35,43].

## Assay design

The common reverse primer and two forward allele-specific MAMA primers for each assay were designed using NetPrimer analysis software (Premier Biosoft, Palo Alto, CA). One forward primer represents the original SNP allele, referred to as “ancestral” and the other represents the mutated SNP allele, referred to as “derived”. The MAMA primers for each assay were designed to compete for the same SNP locus on the template and the resulting amplicon product is generated by the allele-specific MAMA primer that most closely matches the template.

To differentiate between the amplicon products of the derived and ancestral genotypes, additional length was added on the derived amplicon product but not the ancestral product. This was achieved by adding 21–28 oligonucleotides rich in GC content at the 5’end of the derived MAMA forward primer (GC-clamp) (Fig 3). Since the primers are incorporated in the final amplicon product, the addition of the GC-clamp on the derived MAMA forward primer resulted in derived amplicon products that were 21–28 bp longer than their ancestral amplicon counterparts.
<table><caption>Table 1. Y. pestis isolates from Madagascar used in this study.</caption><thead><tr><th>ID</th><th>Original ID</th><th>DNA Type</th><th>pla status</th><th>3a status</th><th>Tested at IPM?</th><th>Group<sup>a</sup></th><th>Vogler et al. 2011 SNP Node<sup>a</sup></th><th>Vogler et al. 2017 SNP Node<sup>a</sup></th><th>Originating Lab<sup>b</sup></th><th>Province</th><th>District</th><th>Commune</th><th>Village</th><th>Source</th><th>Year</th><th>Previously published<sup>c</sup></th></tr></thead><tbody><tr><td>Yp2554</td><td>194/95</td><td>bacterial culture</td><td>+</td><td>+</td><td>yes</td><td>I</td><td>q05</td><td>q05</td><td>IPM</td><td>Antananarivo</td><td>Ajizorozobe</td><td>Ambongamarina</td><td>Ambohimiaramanana</td><td>human</td><td>1995</td><td>Vogler et al. 2017</td></tr><tr><td>Yp2564</td><td>67/96</td><td>bacterial culture</td><td>+</td><td>+</td><td>yes</td><td>I</td><td>r</td><td>r</td><td>IPM</td><td>Antananarivo</td><td>Antsirabe II</td><td>Vinaninkarena</td><td>Fiakarandava</td><td>human</td><td>1996</td><td>Vogler et al. 2017</td></tr><tr><td>Yp2669</td><td>172/98</td><td>bacterial culture</td><td>+</td><td>+</td><td>yes</td><td>I</td><td>s9</td><td>s09</td><td>IPM</td><td>Mahajanga</td><td>Mahajanga I</td><td>Mahabibo</td><td>Abattoir</td><td>Suncus murinus</td><td>1998</td><td>Vogler et al. 2017</td></tr><tr><td>Yp2975</td><td>145/06</td><td>bacterial culture</td><td>+</td><td>+</td><td>yes</td><td>I</td><td>q2</td><td>q10</td><td>IPM</td><td>Toamasina</td><td>Moramanga</td><td>Andasibe</td><td>Ampangalantsary</td><td>human</td><td>2006</td><td>Vogler et al. 2017</td></tr><tr><td>Yp2610</td><td>402/97</td><td>bacterial culture</td><td>+</td><td>+</td><td>yes</td><td>I</td><td>q06</td><td>q06</td><td>IPM</td><td>Antananarivo</td><td>Antsirabe II</td><td>lbity</td><td>Ampopoka</td><td>human</td><td>1997</td><td>Vogler et al. 2017</td></tr><tr><td>Yp2924</td><td>26/05</td><td>bacterial culture</td><td>+</td><td>+</td><td>yes</td><td>I</td><td>q06</td><td>q06</td><td>IPM</td><td>Antananarivo</td><td>Manjakandriana</td><td>Ambatomena</td><td>Antsahakely</td><td>human</td><td>2005</td><td>Vogler et al. 2017</td></tr><tr><td>Yp2559</td><td>17/96</td><td>bacterial culture</td><td>+</td><td>+</td><td>yes</td><td>II</td><td>d</td><td>v01</td><td>IPM</td><td>Antananarivo</td><td>Betato</td><td>Mandritsara</td><td>Miadana</td><td>human</td><td>1996</td><td>Vogler et al. 2017</td></tr><tr><td>Yp2560</td><td>19/96</td><td>bacterial culture</td><td>+</td><td>+</td><td>yes</td><td>II</td><td>d</td><td>v01</td><td>IPM</td><td>Antananarivo</td><td>Betato</td><td>Madana Maranasoand</td><td>Ambohijato</td><td>human</td><td>1996</td><td>Vogler et al. 2017</td></tr><tr><td>Yp2561</td><td>26/96</td><td>bacterial culture</td><td>+</td><td>+</td><td>yes</td><td>II</td><td>d</td><td>v02</td><td>IPM</td><td>Antananarivo</td><td>Betato</td><td>Tritriva</td><td>Ivonarivo</td><td>human</td><td>1996</td><td>Vogler et al. 2017</td></tr><tr><td>Yp2568</td><td>119/96</td><td>bacterial culture</td><td>+</td><td>+</td><td>yes</td><td>II</td><td>h1</td><td>h11</td><td>IPM</td><td>Antananarivo</td><td>Antsirabe II</td><td>Margarano</td><td>Maninarivo</td><td>human</td><td>1996</td><td>Vogler et al. 2017</td></tr><tr><td>Yp2572</td><td>154/96</td><td>bacterial culture</td><td>+</td><td>+</td><td>yes</td><td>II</td><td>h1</td><td>h11</td><td>IPM</td><td>Antananarivo</td><td>Antsirabe I</td><td>Antsenakely Andraikba I</td><td>Andrany</td><td>human</td><td>1996</td><td>Vogler et al. 2017</td></tr><tr><td>Yp2672</td><td>32/99</td><td>bacterial culture</td><td>+</td><td>+</td><td>yes</td><td>I</td><td>11</td><td>101</td><td>IPM</td><td>Antananarivo</td><td>Antsirabe II</td><td>Alatsininy Ibity</td><td>Sahamalaza</td><td>human</td><td>1999</td><td>Vogler et al. 2017</td></tr><tr><td>Yp2675</td><td>63/99</td><td>bacterial culture</td><td>+</td><td>+</td><td>yes</td><td>I</td><td>s03</td><td>s03</td><td>IPM</td><td>Antananarivo</td><td>Betato</td><td>Amb himanambola</td><td>Vinanisoa</td><td>human</td><td>1999</td><td>Vogler et al. 2017</td></tr><tr><td>Yp2933</td><td>71/05</td><td>bacterial culture</td><td>+</td><td>+</td><td>yes</td><td>I</td><td>s05</td><td>s05</td><td>IPM</td><td>Antananarivo</td><td>Manjakandriana</td><td>Marjakandriana</td><td>Anosimanarivo</td><td>human</td><td>2005</td><td>Vogler et al. 2017</td></tr><tr><td>Yp3020</td><td>02/08</td><td>bacterial culture</td><td>+</td><td>+</td><td>yes</td><td>I</td><td>s05</td><td>s05</td><td>IPM</td><td>Antananarivo</td><td>Manjakandriana</td><td>Manjakandriana</td><td>Manakasikely</td><td>human</td><td>2008</td><td>Vogler et al. 2017</td></tr><tr><td>Yp3070</td><td>87/11</td><td>bacterial culture</td><td>+</td><td>+</td><td>yes</td><td>I</td><td>j</td><td>j10</td><td>IPM</td><td>Antananarivo</td><td>Betato</td><td>Mandoto</td><td>Analavory</td><td>human</td><td>2011</td><td>Vogler et al. 2017</td></tr><tr><td>Yp2483<sup>a</sup></td><td>644/M 004</td><td>complex clinical sample</td><td>+</td><td>+</td><td>no</td><td>I</td><td>k</td><td>k</td><td>IPM/BIM<sup>b</sup></td><td>Fianaranisca</td><td>Amthalavao</td><td>Sendrisoa</td><td>Manakony</td><td>human</td><td>2007</td><td>Riehm et al. 2015</td></tr><tr><td>Yp2486<sup>a</sup></td><td>571/M 044</td><td>complex clinical sample</td><td>+</td><td>+</td><td>no</td><td>I</td><td>q3</td><td>q12</td><td>IPM/BIM<sup>b</sup></td><td>Mahajanga</td><td>Tsaratanana</td><td>Keliloha</td><td>Ambatomitsangana</td><td>human</td><td>2007</td><td>Riehm et al. 2015</td></tr><tr><td>Yp3181<sup>a</sup></td><td>TLO432/13</td><td>complex clinical sample</td><td>+</td><td>-</td><td>no</td><td>I</td><td>s3</td><td>s13</td><td>IPM/BIM<sup>b</sup></td><td>Antananarivo</td><td>Faratisiho</td><td>Ramainandro</td><td>Alatsininy Bevohoka</td><td>human</td><td>2007</td><td>Riehm et al. 2015</td></tr><tr><td>Yp3182<sup>a</sup></td><td>TLO464/13</td><td>complex clinical sample</td><td>+</td><td>+</td><td>no</td><td>unknown</td><td>unknown</td><td>unknown</td><td>IPM</td><td>Mahajanga</td><td>Mandritsara</td><td>Amptatakamaroreny</td><td>Beranimbo</td><td>human</td><td>2013</td><td>-</td></tr></tbody></table>

<sup>a</sup>Strain ID in the Northern Arizona University DNA collection.

<sup>b</sup>Strain ID from the originating laboratory.

<sup>c</sup>Indicates Group I or II, as described (Vogler et al. 2011).

<sup>d</sup>Indicates SNP determined node (Vogler et al. 2011).

<sup>e</sup>Indicates SNP determined node (Vogler et al. 2017). For samples where a specific node could not be determined, the lineage letter is indicated without a specific node number.

<sup>f</sup>IPM: Institut Pasteur de Madagascar, BIM: Institut für Mikrobiologie der Bundeswehr.

<sup>g</sup>These samples were analyzed previously in the listed references.
<table><caption>Table 2. The SNP allele state profiles of genetic subgroups targeted for agarose MAMA design where the Derived SNP state is shaded and Ancestral state is unshaded.</caption><thead><tr><th rowspan="2">Major Group<sup>a</sup></th><th colspan="2">Subgroup<sup>a</sup></th><th>Tier 1<sup>b</sup></th><th>Tier 2<sup>b</sup></th><th>Tier 2<sup>b</sup></th><th>Tier 3<sup>b</sup></th><th>Tier 4<sup>b</sup></th><th>Tier 5<sup>b</sup></th><th>Tier 6<sup>b</sup></th><th>Tier 7<sup>b</sup></th><th>Tier 2<sup>b</sup></th><th>Tier 4<sup>b</sup></th><th>Tier 5<sup>b</sup></th><th>Tier 6<sup>b</sup></th><th>Tier 7<sup>b</sup></th><th>Tier 2<sup>b</sup></th><th>Tier 3<sup>b</sup></th></tr><tr><th></th><th></th><th>*Mad-05<sup>c</sup></th><th>Tier 12<sup>c</sup></th><th>Tier 2<sup>b</sup></th><th>Tier 2<sup>b</sup></th><th>Tier 3<sup>b</sup></th><th>Tier 4<sup>b</sup></th><th>Tier 5<sup>b</sup></th><th>Tier 6<sup>b</sup></th><th>Tier 7<sup>b</sup></th><th>Tier 2<sup>b</sup></th><th>Tier 3<sup>b</sup></th><th>Tier 5<sup>b</sup></th><th>Tier 6<sup>b</sup></th><th>Tier 7<sup>b</sup></th><th>Tier 2<sup>b</sup></th><th>Tier 3<sup>b</sup></th></tr><tr><th></th><th></th><th></th><th>C-T</th><th>G-A</th><th>C</th><th>A-G</th><th>G-A</th><th>T</th><th>C</th><th>Mad-64<sup>c</sup></th><th>Mad-65<sup>c</sup></th><th>C-T</th><th>G-A</th><th>G-A</th><th>T</th><th>C</th><th>T</th></tr></thead><tbody><tr><td>I</td><td>k<sup>2</sup></td><td>T</td><td></td><td>G</td><td>C</td><td>A</td><td>G</td><td>T</td><td>C</td><td>G</td><td>G</td><td>C</td><td>G</td><td>G</td><td>T</td><td>G</td><td>T</td><td>T</td></tr><tr><td>I</td><td>11</td><td>T</td><td>A</td><td>C</td><td>C</td><td>A</td><td>G</td><td>T</td><td>C</td><td>G</td><td>G</td><td>C</td><td>G</td><td>G</td><td>T</td><td>G</td><td>T</td><td>T</td></tr><tr><td>I</td><td>j</td><td>T</td><td>G</td><td>T</td><td>A</td><td>A</td><td>G</td><td>T</td><td>C</td><td>G</td><td>C</td><td>G</td><td>G</td><td>G</td><td>T</td><td>G</td><td>T</td><td>T</td></tr><tr><td>I</td><td>s2</td><td>T</td><td>G</td><td>C</td><td>G</td><td>A</td><td>G</td><td>T</td><td>C</td><td>G</td><td>C</td><td>G</td><td>G</td><td>G</td><td>T</td><td>G</td><td>T</td><td>T</td></tr><tr><td>I</td><td>s3</td><td>T</td><td>G</td><td>C</td><td>G</td><td>A</td><td>G</td><td>T</td><td>C</td><td>G</td><td>C</td><td>G</td><td>G</td><td>G</td><td>T</td><td>G</td><td>T</td><td>T</td></tr><tr><td>I</td><td>s4</td><td>T</td><td>G</td><td>C</td><td>G</td><td>A</td><td>G</td><td>T</td><td>C</td><td>G</td><td>C</td><td>G</td><td>G</td><td>G</td><td>T</td><td>G</td><td>T</td><td>T</td></tr><tr><td>I</td><td>s5</td><td>T</td><td>G</td><td>C</td><td>G</td><td>A</td><td>C</td><td>T</td><td>C</td><td>G</td><td>C</td><td>G</td><td>G</td><td>G</td><td>T</td><td>G</td><td>T</td><td>T</td></tr><tr><td>I</td><td>s8</td><td>T</td><td>G</td><td>C</td><td>G</td><td>A</td><td>G</td><td>T</td><td>C</td><td>G</td><td>C</td><td>G</td><td>G</td><td>G</td><td>T</td><td>G</td><td>T</td><td>T</td></tr><tr><td>I</td><td>s9</td><td>T</td><td>G</td><td>C</td><td>G</td><td>A</td><td>G</td><td>T</td><td>C</td><td>G</td><td>A</td><td>T</td><td>G</td><td>T</td><td>G</td><td>T</td><td>T</td><td>T</td></tr><tr><td>I</td><td>q1</td><td>T</td><td>G</td><td>C</td><td>A</td><td>A</td><td>G</td><td>T</td><td>C</td><td>G</td><td>A</td><td>T</td><td>G</td><td>T</td><td>G</td><td>T</td><td>T</td><td>T</td></tr><tr><td>I</td><td>q2</td><td>T</td><td>G</td><td>C</td><td>A</td><td>A</td><td>G</td><td>T</td><td>C</td><td>G</td><td>A</td><td>A</td><td>G</td><td>T</td><td>G</td><td>T</td><td>T</td><td>T</td></tr><tr><td>I</td><td>q3</td><td>T</td><td>G</td><td>C</td><td>A</td><td>A</td><td>G</td><td>T</td><td>C</td><td>G</td><td>A</td><td>A</td><td>A</td><td>A</td><td>T</td><td>G</td><td>T</td><td>T</td></tr><tr><td>I</td><td>q4</td><td>T</td><td>G</td><td>C</td><td>A</td><td>A</td><td>G</td><td>T</td><td>C</td><td>G</td><td>A</td><td>A</td><td>A</td><td>A</td><td>T</td><td>G</td><td>T</td><td>T</td></tr><tr><td>I</td><td>q5</td><td>T</td><td>G</td><td>C</td><td>A</td><td>A</td><td>G</td><td>T</td><td>C</td><td>G</td><td>A</td><td>A</td><td>A</td><td>A</td><td>T</td><td>G</td><td>T</td><td>T</td></tr><tr><td>I</td><td>q6</td><td>T</td><td>G</td><td>C</td><td>A</td><td>A</td><td>G</td><td>T</td><td>C</td><td>G</td><td>A</td><td>A</td><td>A</td><td>A</td><td>T</td><td>G</td><td>T</td><td>T</td></tr><tr><td>I</td><td>r</td><td>T</td><td>G</td><td>C</td><td>A</td><td>A</td><td>G</td><td>T</td><td>C</td><td>G</td><td>A</td><td>A</td><td>A</td><td>A</td><td>T</td><td>G</td><td>T</td><td>T</td></tr><tr><td>II</td><td>d</td><td>T</td><td>G</td><td>C</td><td>A</td><td>A</td><td>G</td><td>T</td><td>C</td><td>G</td><td>A</td><td>A</td><td>A</td><td>A</td><td>T</td><td>G</td><td>T</td><td>T</td></tr><tr><td>II</td><td>h1</td><td>T</td><td>G</td><td>C</td><td>A</td><td>A</td><td>G</td><td>T</td><td>C</td><td>G</td><td>A</td><td>A</td><td>A</td><td>A</td><td>T</td><td>G</td><td>T</td><td>T</td></tr></tbody></table>

<sup>a</sup>Described in Vogler *et al* 2011 and Vogler *et al* 2017

<sup>b</sup>Corresponds to levels indicated in Fig 2 in this study.

<sup>c</sup>canSNP that define Malagasy subgroups.

<sup>d</sup> Collapsed node containing many subgroups.

https://doi.org/10.1371/journal.pntd.0006077.t002
To maximize the visible size differences between the two allelic-specific amplicons when viewed on an agarose gel, the size of the PCR amplicon was restricted to ≤ 80 bases total length. Amplicons within ~80 bases show the greatest migration difference on a gel when small size differences exist. This is the case with our derived and ancestral allele-specific PCR products, which differ between 21–28 bases (Fig 3, Table 3, S2 Appendix). Additionally, through the use of the GC-clamp, our assays retained the capability of SNP genotype discrimination on a real-time PCR platform, which is based on differential melt-curve properties of each SNP-specific PCR product [20].

## PCR conditions and size based genotype discrimination

Initial PCR conditions were identified at NAU. PCR conditions for different assays varied and are described in Table 3. PCR amplification per assay was carried out in 20 µL volume with the following reagents (see S1 Appendix for volume of each reagent): for one reaction, 1x PCR buffer without MgCl₂, MgCl₂ range of 1.5–2.5 mM, 0.30 mM deoxynucleoside triphosphate, 1.6 units of platinum *Taq* DNA polymerase (Invitrogen, Carlsbad, CA), both sets of forward MAMA primers (derived and ancestral allele-specific) with one common reverse primer at 0.40 µL each (for a 1:1 ratio), molecular grade water to achieve 18 µL total volume and 2.0 µL of diluted DNA template at ~1ng/µL per reaction. We did not directly test our assays on genomic DNA concentrations below 100 pg but previously published work suggests that the MAMA approach is sensitive to DNA amounts below 100 pg [20]. For each set of reactions, at least one of each ancestral and derived allele templates were used as positive controls as well as at least two no-template controls (NTC). Thermal cycling parameters for the eighteen assays are as follows: initial denaturation at 94°C for 5 min followed by 30–40 cycles of 94°C for 30 s, 51°C -67.3°C (Table 3) for 30 s, and 72°C for 30 s, with a final extension at 72°C for 5 min. All PCR amplifications were performed with a MJ Research PTC 200 thermal cycler (BioRad, Hercules, CA).

Conditions of agarose gel electrophoresis for PCR amplicons included adding 4 µL of 6x loading dye 0.25 w/v xylene cyanol FF and 30% v/v glycerol, water (Thermo Fisher scientific, Waltham, MA) to individual PCR products to achieve a 1x final dye concentration. Individual reactions (20 µL) mixed with loading dye were loaded onto a 2% agarose gel matrix; 100 bp DNA ladder (Invitrogen, Carlsbad, CA) was used for size referencing. Gels were prepared in 1x lithium borate buffer (S2 Appendix) and stained with SybrSafe dye (Life Technologies, Carlsbad, CA). Electrophoresis was conducted at 300V for 25–30 minutes and viewed under UV transillumination.

The genotyping accuracy of our SNP assays was validated using positive DNA controls that represented known ancestral and derived allele states for each SNP target. Assay accuracy was further assessed by testing them on 16 genomic DNA extracts of Malagasy *Y. pestis* strains belonging to known diverse phylogenetic groups (Table 1) based on amplicon sequencing [34]. To assess assay capability to genotype *Y. pestis* directly from complex DNA samples (containing high levels of host DNA), we tested the performance of each assay on a positive control sample comprised of high concentrations of human DNA with our positive control *Y. pestis* DNA extract (strain A1122). To further assess this capability, we tested two assays (Mad-05 and Mad-43) on five human complex samples confirmed to be plague positive. This confirmation was based on *Y. pestis*-specific TaqMan assay targeting a high copy plasmid *pla* gene [37]. We assessed the limit of detection of our MAMA tools on these five human clinical samples by testing them on a new TaqMan assay designed for the chromosomal *3a* gene (Fig 4 and S3 Appendix).

To confirm the specificity of the assays to the *Y. pestis* genome we tested assay performance on DNA extracts of *Bacillus anthracis*, human DNA background, and no template water
<table><caption>Table 3. Agarose-MAMA primers targeting 18 previously published SNP positions.</caption><thead><tr><th>SNP<sup>a</sup></th><th>CO-92 Position<sup>a</sup></th><th>SNP State (A→D)<sup>b</sup></th><th>Vogler et al. 2011 and 2013 SNP group</th><th>Vogler et al. 2017 SNP group</th><th>Primers<sup>c,d</sup></th><th>eConc. (μM)</th><th>NAU-Agarose Annealing temperature (°C) and Cycle Repeats</th><th>Facility site of validation</th></tr></thead><tbody><tr><td>Mad-05</td><td>3525070</td><td>C→T</td><td>All</td><td>NA</td><td>A: GTGGCTGGCAGCGGTCIC<br>D: cgggttcgggttcgggttcggggttcgggGTGGCTGGCAGCGGTCaT<br>C: GAAGTGAACAAAATGCGACTAATA</td><td>0.200<br>0.200<br>0.200</td><td>60°C<br>33 cycles</td><td>NAU and<br>IPM</td></tr><tr><td>Mad-43</td><td>1348724</td><td>C→T</td><td>d/e/i/h1/h2</td><td>k-d</td><td>A: GGTTGCTGATGAACACGGIC<br>D: cgggttcgggttcgggttcgggGTTGCTGATGAACACGGCT<br>C: GCAACCAAATCAGCAAAATAGAGA</td><td>0.200<br>0.200<br>0.200</td><td>60°C<br>33 cycles</td><td>NAU and<br>IPM<sup>g</sup></td></tr><tr><td>Mad-46</td><td>1791464</td><td>T→C</td><td>h1/h2</td><td>hu-h</td><td>A: ATACAGGATTATTTTAAATGGcA<br>D: cgggttcgggttcgggttcgggATACAGGATTATTTTAAATGGg<br>C: GCATAAGTATTGCAAAATTTAATTTC</td><td>0.200<br>0.200<br>0.200</td><td>55°C<br>33 cycles</td><td>NAU and<br>IPM<sup>g</sup></td></tr><tr><td>Mad-08</td><td>3210210</td><td>C→T</td><td>j</td><td>k-j01</td><td>A: CAGCTTCACGCGACGACIAc<br>D: cgggttcgggttcggggttcgggAGCTTCACGCGACGACaAT<br>C: GGAGAGGAAGACCATCAACC</td><td>0.200<br>0.200<br>0.200</td><td>67.3°C<br>33 cycles</td><td>NAU and<br>IPM</td></tr><tr><td>Mad-12</td><td>805205</td><td>G→A</td><td>I1/I2</td><td>k-I01</td><td>A: GGCTGTGGATGCGGGGTaG<br>D: cgggttcgggttcgggttcggggttcgggGGCTGTGGATGCGGGGTa<br>C: GTTGTCCATCGGTAGCATCTTG</td><td>0.200<br>0.200<br>0.200</td><td>60°C<br>35 cycles</td><td>NAU and<br>IPM<sup>g</sup></td></tr><tr><td>Mad-42</td><td>3926933</td><td>T→C</td><td>r</td><td>k-r</td><td>A: ATAGTACAATACAGTAAAGTGACAATAAgT<br>D: cgggttcgggttcgggttcgggGTAACATACAGTAAAGTGACAATAAC<br>C: CTAAATAGAGTGTACGCTTACAAC</td><td>0.800<br>0.100<br>0.200</td><td>54°C<br>35 cycles</td><td>NAU and<br>IPM<sup>g</sup></td></tr><tr><td>Mad-30</td><td>2899068</td><td>G→A</td><td>q1-q7</td><td>k-q1</td><td>A: GCACGAAACGCCTCATGCGICC<br>D: cgggttcgggttcgggttcgggACGAAACGCCTCATGCGaCT<br>C: TGTTGTGGCGGCGCGGTG</td><td>0.800<br>0.200<br>0.200</td><td>60°C<br>31 cycles</td><td>NAU and<br>IPM<sup>g</sup><br>MgCl<sub>2</sub> =<br>1.5mM</td></tr><tr><td>Mad-32</td><td>1632446</td><td>T→A</td><td>q2-q7</td><td>k-q1</td><td>A: AAACATAACCACCAGCCAAATAgT<br>D: cgggttcgggttcgggttcgggACATAACCACCAGCCAAATcAA<br>C: TGCGGGAGGGGCTTTACT</td><td>0.200<br>0.200<br>0.200</td><td>60°C<br>33 cycles</td><td>NAU and<br>IPM<sup>g</sup></td></tr><tr><td>Mad-35</td><td>3897676</td><td>G→A</td><td>q3-q7</td><td>qI-qII</td><td>A: GAGATGCTCAGCGAGCGAIg<br>D: cgggttcgggttcgggttcgggGATGCTCAGCGAGCGgAA<br>C: AAGGCAGATGCAACGGATAAC</td><td>0.200<br>0.200<br>0.200</td><td><sup>1</sup>MgCl<sub>2</sub> [1.5mM]<br>60°C<br>40 cycles</td><td>NAU and<br>IPM<sup>g</sup><br>MgCl<sub>2</sub> =<br>1.5mM</td></tr><tr><td>Mad-36</td><td>3533308</td><td>C→T</td><td>q4-q7</td><td>qII-q04</td><td>A: ACCCGTGAGCAAAACCcC<br>D: cgggttcgggttcgggttcgggATTTAACCCGTGAGCAAAACaGT<br>C: ATGGCCATAGCAAAAGGTGACA</td><td>0.400<br>0.200<br>0.200</td><td>60°C<br>33 cycles</td><td>NAU and<br>IPM<sup>g</sup><br>MgCl<sub>2</sub> =<br>1.5mM</td></tr><tr><td>Mad-37</td><td>606374</td><td>G→A</td><td>q5</td><td>q04-q05</td><td>A: ATCACCATCCCGCAACGATTAG<br>D: cgggttcgggttcgggttcgggATCACCATCCCGCAACGATAcA<br>C: ACGCCCCAGAACTTTCAATAAG</td><td>0.200<br>0.200<br>0.200</td><td>60°C<br>33 cycles</td><td>NAU and<br>IPM<sup>g</sup></td></tr><tr><td>Mad-38</td><td>1950363</td><td>A→G</td><td>q6/q7</td><td>q04-q06</td><td>A: GAATATCCAAGCGTTGCTGAgT<br>D: cgggttcgggttcgggttcgggAATATCCAAGCGTTGCTGAcC<br>C: TTGCTGGAAGGTGAAATG</td><td>0.800<br>0.200<br>0.200</td><td>MgCl<sub>2</sub> [2.5mM]<br>60°C<br>30 cycles</td><td>NAU and<br>IPM<sup>g</sup></td></tr><tr><td>Mad-58</td><td>2681067</td><td>A→G</td><td>s2-s9</td><td>k-s</td><td>A: AGATGTGGCCAAACACGgA<br>D: cgggttcgggttcgggttcgggAAAGATGTGGCCAAACACGgA<br>C: CGACGAAACGAGTCTTGTTGATAA</td><td>0.400<br>0.200<br>0.200</td><td>60°C<br>35 cycles</td><td>NAU and<br>IPM<sup>g</sup></td></tr><tr><td>Mad-62</td><td>1902978</td><td>G→A</td><td>s3-s9</td><td>s-s03</td><td>A<sup>h</sup>: CAGCATGCACCAGCCCGaGG<br>D: cgggttcgggttcgggttcgggCAGCATGCACCAGCCCGIGA<br>C<sup>h</sup>: GGATGATGCCGTCGGTATTTTC</td><td>0.800<br>0.200<br>0.200</td><td>60°C<br>31 cycles</td><td>NAU</td></tr></tbody></table>

(Continued)
<table><caption>Table 3. (Continued)</caption><thead><tr><th>SNP<sup>a</sup></th><th>CO-92 Position<sup>a</sup></th><th>SNP State (A→D)<sup>b</sup></th><th>Vogler et al. 2011 and 2013 SNP group</th><th>Vogler et al. 2017 SNP group</th><th>Primers<sup>c,d</sup></th><th><sup>e</sup>Conc. (µM)</th><th>NAU-Agarose Annealing temperature (°C) and Cycle Repeats</th><th>Facility site of validation</th></tr></thead><tbody><tr><td>Mad-64</td><td>2393962</td><td>T→C</td><td>s4-s9</td><td>s03-s04</td><td>A<sup>h</sup>: GACGCACGGTGTGACAATAAAAgTT<br/>D: cgggttcggggttcgggttcggggGACGCACGGTGTGACAATAAAcTC<br/>C<sup>h</sup>: GCTTTTTGCATTTGGATTCTCCTT</td><td>0.800<br/>0.200<br/>0.200</td><td>64°C<br/>35 cycles<br/>Note: double band on derived genotype</td><td>NAU</td></tr><tr><td>Mad-65</td><td>435838</td><td>C→T</td><td>s5-s9</td><td>s04-s05</td><td>A<sup>h</sup>: GCTGAGAAAAATTATGACAACTTTTATCgCC<br/>D: cgggttcgggttcgggttcgggCTGAGAAAAATTATGACAACTTTATCcCT<br/>C<sup>h</sup>: TGCTGCACAGCATTTCACTACGG</td><td>0.300<br/>0.200<br/>0.200</td><td>60°C<br/>30 cycles</td><td>NAU</td></tr><tr><td>Mad-74</td><td>1355526</td><td>G→T</td><td>s8/s9</td><td>s18-s08</td><td>A: GAACCATAATTCACAATAAAAgG<br/>D: cgggttcggggttcggggttcggggAAGAACCATAATTCACAATAAAAcT<br/>C: TATCAGTTCTAATTTGAATGTAAGG</td><td>0.800<br/>0.200<br/>0.200</td><td>51°C<br/>33 cycles</td><td>NAU</td></tr><tr><td>Mad-78</td><td>1780902</td><td>C→T</td><td>s9</td><td>s08-s09</td><td>A: ATTGCCAAGCTATTTGTCCG<br/>D: cgggttcggggttcggggttcgggATAACATTGCCAAGCTATTTGTgTA<br/>C: GTGTCG ATGCGATAGCACTAATA</td><td>0.200<br/>0.100<br/>0.200</td><td>60°C<br/>35 cycles</td><td>NAU</td></tr></tbody></table>

<footnote><sup>a</sup>SNP position and naming as published in Vogler et al. 2011, 2013, 2017.</footnote>

<footnote><sup>b</sup>SNP states are presented according to their orientation in the Yp CO-92 reference genome; D: Derived SNP state; A: Ancestral SNP state.</footnote>

<footnote><sup>c</sup>D: Derived; A: Ancestral; C: Common Reverse.</footnote>

<footnote><sup>d</sup>Primer tails and antepenultimate/and or penultimate mismatch bases are in lower case.</footnote>

<footnote><sup>e</sup>Ratios are skewed in some assays.</footnote>

<footnote><sup>f</sup>Standard MgCl<sub>2</sub> is 2 mM.</footnote>

<footnote><sup>g</sup>MgCl<sub>2</sub> reduced to 1.5 mM at IPM from the standard 2 mM used at NAU.</footnote>

<footnote><sup>h</sup>Primer originally published in Vogler et al 2013.</footnote>

<footnote>https://doi.org/10.1371/journal.pntd.0006077.t003</footnote>
control. To compare assay performance at different research facilities, NAU and IPM jointly conducted a second validation study on the IPM laboratory premises using IPM PCR reagents and instrument (Fig 5). Thirteen of the eighteen assays were selected for validation (Table 3). At IPM, specificity was confirmed by testing on *Leptospira interrogans* serovar *Canicola* and no template controls. PCR was prepared for each assay using the *Taq* Core Kits 10, Cat# EPTQK300 PCR reagents (MP Biomedicals, Santa Ana, CA) and amplification was conducted using AB Applied Biosystems Veriti 96 Well Thermal Cycler thermal cycle (ThermoFisher Scientific, Waltham, MA). Electrophoresis was performed using 2% agarose gels visualized on a Gelscan (Bio-Rad, Hercules, CA).

**Supporting information**

**S1 Appendix. PCR master mix calculation sheet indicating final concentration and volume addition of each reagent.**
(DOCX)

**S2 Appendix. Side by side comparison of 2% agarose gel electrophoresis in a 1x Lithium Borate and 1x TAE matrix.** Electrophoresis was conducted with identical PCR products, equipment and agarose gel conditions. The two runs differ only in buffer solution used.
(DOCX)

**S3 Appendix. *Y. pestis* 3a TaqMan assay design.**
(DOCX)

**S4 Appendix. SNP allele state for the 16 isolates shown in Table 1.**
(DOCX)

**S5 Appendix. Agarose gel pictures for each assays generated at both NAU and IPM or NAU only.**
(DOCX)

**S6 Appendix. Agarose gel pictures of negative controls (human background DNA and *Bacillus anthracis* A0635 and A0643) on a subset of MAMA tools.**
(DOCX)

**Acknowledgments**

We thank Dr. Amy J. Vogler for sharing her considerable expertise to understand phylogenetic groups in Madagascar and the review of this manuscript. We thank Isaiah K. Self, Fehivola Andriamiarimanana, and Soloandry Rahajandraibe for technical assistance.

**Author Contributions**

**Conceptualization:** Minoarisoa Rajerison, Dawn N. Birdsell.

**Formal analysis:** Cedar L. Mitchell, Voahangy Andrianaivoarimanana, Dawn N. Birdsell.

**Funding acquisition:** Cedar L. Mitchell, Paul S. Keim, David M. Wagner.

**Methodology:** Cedar L. Mitchell, Voahangy Andrianaivoarimanana, Rebecca E. Colman, Dawn N. Birdsell.

**Supervision:** Voahangy Andrianaivoarimanana, Minoarisoa Rajerison, Dawn N. Birdsell.

**Validation:** Cedar L. Mitchell, Voahangy Andrianaivoarimanana, Rebecca E. Colman, Dawn N. Birdsell.
**Writing – original draft:** Cedar L. Mitchell, Joseph Busch, Heidie Hornstra-O'Neill, David M. Wagner, Dawn N. Birdsell.

**Writing – review & editing:** Cedar L. Mitchell, Voahangy Andrianavoarimanana, Rebecca E. Colman, Joseph Busch, Heidie Hornstra-O'Neill, Paul S. Keim, David M. Wagner, Minoarisoa Rajerison, Dawn N. Birdsell.

## References

1. Belmont JW, Hardenbol P, Willis T, Yu F, Yang H, et al. (2003) The international HapMap project. Nature 426: 789–796. https://doi.org/10.1038/nature02168 PMID: 14685227

2. Abecasis G, Auton A, Brooks L, DePristo M, Durbin R, et al. (2012) 1000 Genomes Project Consortium: An integrated map of genetic variation from 1,092 human genomes. Nature 491: 56–65. https://doi.org/10.1038/nature11632 PMID: 23128226

3. Keim P, Van Ert MN, Pearson T, Vogler AJ, Huynh LY, et al. (2004) Anthrax molecular epidemiology and forensics: using the appropriate marker for different evolutionary scales. Infect Genet Evol 4: 205–213. https://doi.org/10.1016/j.meegid.2004.02.005 PMID: 15450200

4. Ozsurekci Y, Birdsell DN, Celik M, Karadag-Oncel E, Johansson A, et al. (2015) Diverse Francisella tularensis strains and oropharyngeal tularemia, Turkey. Emerg Infect Dis 21: 173–175. https://doi.org/10.3201/eid2101.141087 PMID: 25531237

5. Wahab T, Birdsell DN, Hjertqvist M, Mitchell CL, Wagner DM, et al. (2014) Insights to genetic characterization tools for epidemiological tracking of Francisella tularensis in Sweden. PLoS One 9: e112167. https://doi.org/10.1371/journal.pone.0112167 PMID: 25401326

6. Morris GP, Ramu P, Deshpande SP, Hash CT, Shah T, et al. (2013) Population genomic and genome-wide association studies of agroclimatic traits in sorghum. Proc Natl Acad Sci U S A 110: 453–458. https://doi.org/10.1073/pnas.1215985110 PMID: 23267105

7. Seddon JM, Parker HG, Ostrander EA, Ellegren H (2005) SNPs in ecological and conservation studies: a test in the Scandinavian wolf population. Mol Ecol 14: 503–511. https://doi.org/10.1111/j.1365-294X.2005.02435.x PMID: 15660941

8. Farhat MR, Shapiro BJ, Kieser KJ, Sultana R, Jacobson KR, et al. (2013) Genomic analysis identifies targets of convergent positive selection in drug-resistant Mycobacterium tuberculosis. Nat Genet 45: 1183–1189. https://doi.org/10.1038/ng.2747 PMID: 23995135

9. Yamaguchi T, Suzuki S, Okamura S, Miura Y, Tsukimori A, et al. (2015) Evolution and single-nucleotide polymorphisms in methicillin-resistant Staphylococcus aureus strains with reduced susceptibility to vancomycin and daptomycin, based on determination of the complete genome. Antimicrob Agents Chemother 59: 3585–3587. https://doi.org/10.1128/AAC.05159-14 PMID: 25779584

10. Keim PS, Wagner DM (2009) Humans and evolutionary and ecological forces shaped the phylogeography of recently emerged diseases. Nat Rev Microbiol 7: 813–821. https://doi.org/10.1038/nrmicro2219 PMID: 19820723

11. Vogler AJ, Chan F, Wagner DM, Roumagnac P, Lee J, et al. (2011) Phylogeography and molecular epidemiology of Yersinia pestis in Madagascar. PLoS Negl Trop Dis 5: e1319. https://doi.org/10.1371/journal.pntd.0001319 PMID: 21931876

12. Morelli G, Song Y, Mazzoni CJ, Eppinger M, Roumagnac P, et al. (2010) Yersinia pestis genome sequencing identifies patterns of global phylogenetic diversity. Nat Genet 42: 1140–1143. https://doi.org/10.1038/ng.705 PMID: 21037571

13. Kilic S, Birdsell DN, Karagoz A, Celebi B, Bakkaloglu Z, et al. (2015) Water as Source of Francisella tularensis Infection in Humans, Turkey. Emerg Infect Dis 21: 2213–2216. https://doi.org/10.3201/eid2112.150634 PMID: 26583383

14. Hendriksen RS, Price LB, Schupp JM, Gillece JD, Kaas RS, et al. (2011) Population genetics of Vibrio cholerae from Nepal in 2010: evidence on the origin of the Haitian outbreak. MBio 2: e00157–00111. https://doi.org/10.1128/mBio.00157-11 PMID: 21862630

15. Black WCt, Vontas JG (2007) Affordable assays for genotyping single nucleotide polymorphisms in insects. Insect Mol Biol 16: 377–387. https://doi.org/10.1111/j.1365-2583.2007.00736.x PMID: 17488301

16. Kim S, Misra A (2007) SNP genotyping: technologies and biomedical applications. Annu Rev Biomed Eng 9: 289–320. https://doi.org/10.1146/annurev.bioeng.9.060906.152037 PMID: 17391067

17. Shi MM (2001) Enabling large-scale pharmacogenetic studies by high-throughput mutation detection and genotyping technologies. Clin Chem 47: 164–172. PMID: 11159763
18. Ding C, Jin S (2009) High-throughput methods for SNP genotyping. Methods Mol Biol 578: 245–254.
https://doi.org/10.1007/978-1-60327-411-1_16 PMID: 19768599

19. Birdsell DN, Vogler AJ, Buchhagen J, Clare A, Kaufman E, et al. (2014) TaqMan real-time PCR assays for single-nucleotide polymorphisms which identify Francisella tularensis and its subspecies and sub-populations. PLoS One 9: e107964. https://doi.org/10.1371/journal.pone.0107964 PMID: 25238067

20. Birdsell DN, Pearson T, Price EP, Hornstra HM, Nera RD, et al. (2012) Melt analysis of mismatch amplification mutation assays (Melt-MAMA): a functional study of a cost-effective SNP genotyping assay in bacterial models. PLoS One 7: e32866. https://doi.org/10.1371/journal.pone.0032866 PMID: 22438886

21. Foster JT, Okinaka RT, Svensson R, Shaw K, De BK, et al. (2008) Real-time PCR assays of single-nucleotide polymorphisms defining the major Brucella clades. J Clin Microbiol 46: 296–301. https://doi.org/10.1128/JCM.01496-07 PMID: 18032628

22. Gaudet M, Fara AG, Beritognolo I, Sabatti M (2009) Allele-specific PCR in SNP genotyping. Methods Mol Biol 578: 415–424. https://doi.org/10.1007/978-1-60327-411-1_26 PMID: 19768609

23. Anisimov AP, Lindler LE, Pier GB (2004) Intraspecific diversity of Yersinia pestis. Clin Microbiol Rev 17: 434–464. https://doi.org/10.1128/CMR.17.2.434-464.2004 PMID: 15084509

24. WHO (2010) Human plague: review of regional morbidity and mortality, 2004–2009. World Health Organization. 37–48 p.

25. Dennis DT, Gage KL, Gratz N, Poland JD, Tikhomirov E (2010) Plague manual: epidemiology, distribution, surveillance and control. World Health Organization.

26. Perry RD, Fetherston JD (1997) Yersinia pestis—etiologic agent of plague. Clin Microbiol Rev 10: 35–66. PMID: 8993858

27. Andrianaitoarimanana V, Kreppel K, Elissa N, Duplantier JM, Carniel E, et al. (2013) Understanding the persistence of plague foci in Madagascar. PLoS Negl Trop Dis 7: e2382. https://doi.org/10.1371/journal.pntd.0002382 PMID: 24244760

28. Richard V, Riehm JM, Herindrainy P, Soanandrasana R, Ratsitoharina M, et al. (2015) Pneumonic plague outbreak, Northern Madagascar, 2011. Emerg Infect Dis 21: 8–15. https://doi.org/10.3201/eid2101.131828 PMID: 25530466

29. World Health Organization (2014) Plague—Madagascar. Emergencies preparedness, response.

30. Heymann DL, Rodier GR (2001) Hot spots in a wired world: WHO surveillance of emerging and re-emerging infectious diseases. Lancet Infect Dis 1: 345–353. https://doi.org/10.1016/S1473-3099(01)00148-7 PMID: 11871807

31. Stoler J, Awandare GA (2016) Febrile illness diagnostics and the malaria-industrial complex: a socio-environmental perspective. BMC Infect Dis 16: 683. https://doi.org/10.1186/s12879-016-2025-x PMID: 27855644

32. Dickmann P, Sheeley H, Lightfoot N (2015) Biosafety and Biosecurity: A Relative Risk-Based Framework for Safer, More Secure, and Sustainable Laboratory Capacity Building. Front Public Health 3: 241. https://doi.org/10.3389/fpubh.2015.00241 PMID: 26539427

33. Paramasivan CN, Lee E, Kao K, Mareka M, Kubendiran G, et al. (2010) Experience establishing tuberculosis laboratory capacity in a developing country setting. Int J Tuberc Lung Dis 14: 59–64. PMID: 20003696

34. Vogler AJ, Andrianaitoarimanana V, Telfer S, Hall CM, Sahl JW, et al. (2017) Temporal phylogeography of Yersinia pestis in Madagascar: Insights into the long-term maintenance of plague. PLoS Negl Trop Dis 11: e0005887. https://doi.org/10.1371/journal.pntd.0005887 PMID: 28873412

35. Vogler AJ, Chan F, Nottingham R, Andersen G, Drees K, et al. (2013) A decade of plague in Maha-janga, Madagascar: insights into the global maritime spread of pandemic plague. MBio 4: e00623–00612. https://doi.org/10.1128/mBio.00623-12 PMID: 23404402

36. Chanteau S, Rahalison L, Ralafiarisoa L, Foulon J, Ratsitorahina M, et al. (2003) Development and testing of a rapid diagnostic test for bubonic and pneumonic plague. Lancet 361: 211–216. https://doi.org/10.1016/S0140-6736(03)12270-2 PMID: 12547544

37. Iqbal SS, Chambers JP, Goode MT, Valdes JJ, Brubaker RR (2000) Detection of Yersinia pestis by pesticide fluorogenic probe-coupled PCR. Mol Cell Probes 14: 109–114. https://doi.org/10.1006/mcpr.2000.0295 PMID: 10799272

38. Nilsson J, Bosnes M, Larsen F, Nygren PA, Uhlen M, et al. (1997) Heat-mediated activation of affinity-immobilized Taq DNA polymerase. Biotechniques 22: 744–751. PMID: 9105627

39. Ramasindrazana B, Andrianaitoarimanana V, Rakotondramanga JM, Birdsell DN, Ratsitorahina M, et al. (2017) Pneumonic Plague Transmission, Moramanga, Madagascar, 2015. Emerg Infect Dis 23: 521–524. https://doi.org/10.3201/eid2303.161406 PMID: 28221119
40. Ratovonjato J, Rajerison M, Rahelinirina S, Boyer S (2014) Yersinia pestis in Pulex irritans fleas during plague outbreak, Madagascar. Emerg Infect Dis 20: 1414–1415. https://doi.org/10.3201/eid2008.130629 PMID: 25061697

41. Rahalison L, Vololonirina E, Ratsitorahina M, Chanteau S (2000) Diagnosis of bubonic plague by PCR in Madagascar under field conditions. J Clin Microbiol 38: 260–263. PMID: 10618097

42. Riehm JM, Rahalison L, Scholz HC, Thoma B, Pfeffer M, et al. (2011) Detection of Yersinia pestis using real-time PCR in patients with suspected bubonic plague. Mol Cell Probes 25: 8–12. https://doi.org/10.1016/j.mcp.2010.09.002 PMID: 20933595

43. Riehm JM, Projahn M, Vogler AJ, Rajerison M, Andersen G, et al. (2015) Diverse Genotypes of Yersinia pestis Caused Plague in Madagascar in 2007. PLoS Negl Trop Dis.