Figures
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.
Conclusions
The success of our study demonstrates that the SNP-based genotyping capacity of laboratories 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 diseases emanating from resource limited regions of the world.
Author summary
Although genetic characterization of pathogenic organisms is a powerful tool for investigating outbreak origins and transmission, associated high upfront costs and demanding technological maintenance exclude this tool for many under-resourced laboratories. Paradoxically, resource constrained regions commonly suffer from high rates of infectious diseases and could benefit most from genetic tracking tools. One such country is Madagascar, 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 distributions. 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.
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: 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 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.
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 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.
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.
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.
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 MgCl2 concentration alteration in the reaction mix. These optimization strategies worked well 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’-cgggttcgggttcgggttcgggttcggg-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.
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.
We obtained maximal assay specificity by a combination of customizing the annealing temperature per assay and/or adjusting the MgCl2 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 MgCl2 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 MgCl2 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 MgCl2 strategy, the raising of annealing temperatures while still maintaining the normal 2.0 mM concentration of MgCl2 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 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).
(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.
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.
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.
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 MgCl2 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 MgCl2 also imparted a general improvement of assay specificity for positive control and test isolates. Surprisingly, the reduction of MgCl2 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 biosafety, 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.
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 MgCl2, MgCl2 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 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.
https://doi.org/10.1371/journal.pntd.0006077.s001
(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.
https://doi.org/10.1371/journal.pntd.0006077.s002
(DOCX)
S3 Appendix. Y. pestis 3a TaqMan assay design.
https://doi.org/10.1371/journal.pntd.0006077.s003
(DOCX)
S4 Appendix. SNP allele state for the 16 isolates shown in Table 1.
https://doi.org/10.1371/journal.pntd.0006077.s004
(DOCX)
S5 Appendix. Agarose gel pictures for each assays generated at both NAU and IPM or NAU only.
https://doi.org/10.1371/journal.pntd.0006077.s005
(DOCX)
S6 Appendix. Agarose gel pictures of negative controls (human background DNA and Bacillus anthracis–A0635 and A0643) on a subset of MAMA tools.
https://doi.org/10.1371/journal.pntd.0006077.s006
(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.
References
- 1. Belmont JW, Hardenbol P, Willis T, Yu F, Yang H, et al. (2003) The international HapMap project. Nature 426: 789–796. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. 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. pmid:21862630
- 15. Black WCt, Vontas JG (2007) Affordable assays for genotyping single nucleotide polymorphisms in insects. Insect Mol Biol 16: 377–387. pmid:17488301
- 16. Kim S, Misra A (2007) SNP genotyping: technologies and biomedical applications. Annu Rev Biomed Eng 9: 289–320. 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. 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 subpopulations. PLoS One 9: e107964. 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. 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. pmid:18032628
- 22. Gaudet M, Fara AG, Beritognolo I, Sabatti M (2009) Allele-specific PCR in SNP genotyping. Methods Mol Biol 578: 415–424. pmid:19768609
- 23. Anisimov AP, Lindler LE, Pier GB (2004) Intraspecific diversity of Yersinia pestis. Clin Microbiol Rev 17: 434–464. 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. Andrianaivoarimanana 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. 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. 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. 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. 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. 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, Andrianaivoarimanana 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. pmid:28873412
- 35. Vogler AJ, Chan F, Nottingham R, Andersen G, Drees K, et al. (2013) A decade of plague in Mahajanga, Madagascar: insights into the global maritime spread of pandemic plague. MBio 4: e00623–00612. 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. pmid:12547544
- 37. Iqbal SS, Chambers JP, Goode MT, Valdes JJ, Brubaker RR (2000) Detection of Yersinia pestis by pesticin fluorogenic probe-coupled PCR. Mol Cell Probes 14: 109–114. 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, Andrianaivoarimanana V, Rakotondramanga JM, Birdsell DN, Ratsitorahina M, et al. (2017) Pneumonic Plague Transmission, Moramanga, Madagascar, 2015. Emerg Infect Dis 23: 521–524. 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. 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. 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.