Skip to main content

Method: a single nucleotide polymorphism genotyping method for Wheat streak mosaic virus



The September 11, 2001 attacks on the World Trade Center and the Pentagon increased the concern about the potential for terrorist attacks on many vulnerable sectors of the US, including agriculture. The concentrated nature of crops, easily obtainable biological agents, and highly detrimental impacts make agroterrorism a potential threat. Although procedures for an effective criminal investigation and attribution following such an attack are available, important enhancements are still needed, one of which is the capability for fine discrimination among pathogen strains. The purpose of this study was to develop a molecular typing assay for use in a forensic investigation, using Wheat streak mosaic virus (WSMV) as a model plant virus.


This genotyping technique utilizes single base primer extension to generate a genetic fingerprint. Fifteen single nucleotide polymorphisms (SNPs) within the coat protein and helper component-protease genes were selected as the genetic markers for this assay. Assay optimization and sensitivity testing was conducted using synthetic targets. WSMV strains and field isolates were collected from regions around the world and used to evaluate the assay for discrimination. The assay specificity was tested against a panel of near-neighbors consisting of genetic and environmental near-neighbors.


Each WSMV strain or field isolate tested produced a unique SNP fingerprint, with the exception of three isolates collected within the same geographic location that produced indistinguishable fingerprints. The results were consistent among replicates, demonstrating the reproducibility of the assay. No SNP fingerprints were generated from organisms included in the near-neighbor panel, suggesting the assay is specific for WSMV. Using synthetic targets, a complete profile could be generated from as low as 7.15 fmoles of cDNA.


The molecular typing method presented is one tool that could be incorporated into the forensic science tool box after a thorough validation study. This method incorporates molecular biology techniques that are already well established in research and diagnostic laboratories, allowing for an easy introduction of this method into existing laboratories. Keywords: single nucleotide polymorphisms, genotyping, plant pathology, viruses, microbial forensics, Single base primer extension, SNaPshot Multiplex Kit


The September 11, 2001 attacks on the World Trade Center and the Pentagon increased concern about terrorist attacks on many vulnerable sectors of the United States (US), including agriculture [1]. The planting density of crops and the facts that biological agents are easily obtainable and can have highly detrimental impacts make agroterrorism a potential threat.

A developing area of forensic science is microbial forensics, ‘a scientific discipline dedicated to analyzing evidence from a bioterrorism act, biocrime, or inadvertent microorganism/toxin release for attribution purposes’ [2]. As the threat of the use of microorganisms as weapons has become more apparent, there is a national effort to increase capability in this new field [2, 3]. Enhancements of the investigative process are still being developed to close gaps that remain in key areas, such as sample collection, packaging, shipping and storage, as well as pathogen identification and discrimination [2, 4].

In addition to identifying the perpetrator of a biological crime, it is also important to identify the source of the pathogen and when and how it was introduced [4]. The stringent demands of a criminal investigation and the pressures of the courtroom require that existing methods of plant pathogen detection and disease diagnosis [5, 6] be revised and validated for forensic use [4].

Molecular comparisons using single nucleotide polymorphisms (SNPs) are often used for evolutionary, ecological and forensic studies [715], and can be especially helpful in typing organisms having high mutation rates, such as RNA viruses, in which mutations are caused primarily by an error-prone RNA-dependent RNA polymerase [16]. Specific regions within the genome can be targeted based on the level of conservation of the nucleic acid sequence. For example, the coat protein (CP) gene of Wheat streak mosaic virus (WSMV) is the most variable region of the genome, with several conserved regions corresponding to the functionally important residues [17]. Other genes with highly variable regions include helper-component protease (HC-Pro) protein 1 (P1), and protein 3 (P3) [17, 18]. The high mutation rate characteristic of RNA viruses contributes to the evolution of this virus and to the presence of multiple genetic variants in a single population. Analyzing the SNPs of these variants can facilitate discrimination among populations, making this procedure a powerful tool that provides benefits that techniques such as sequencing, which focus on identifying the consensus of the population, do not offer.

WSMV, an important pathogen of wheat in the US, was selected as the model pathogen for the development of a molecular typing method. WSMV belongs to the genus Tritimovirus within the family Potyviridae[19]. The 9,384 nucleotide (nt), positive sense, single-stranded RNA virus genome [20] is translated into a single polyprotein that is then processed by three viral proteinases [19]. Because most plant viruses possess an RNA genome, the method described here can be applied almost universally, with only minor modifications, to most plant viruses. Wheat streak mosaic (WSM), the disease caused by WSMV, is found in most of the wheat growing regions of the US nearly every year. Therefore, the high prevalence of WSMV facilitates collection of field isolates from across the nation, as well as internationally. Previous studies have focused on WSMV’s genetic diversity and evolution, including mutation rates, phylogenetic positions and polymorphic sites [16, 18, 19, 2124]. Excellent PCR detection methods have been developed and provide a basis for generating the cDNA that is subjected to the SNP typing method described here [2527]. Extensive knowledge of the viral genome and the availability of robust detection assays facilitate the development of a more stringent, discriminatory assay for variant characterization.

Even though viruses do not represent the majority of plant pathogens, the small size and simplicity of their genomes make them great laboratory models. Bacteria and fungi, the most predominant plant pathogens, present their own unique challenges for genetic analysis, such as genome size, multiple life stages, diploid genomes, and presence of plasmids [28]. However, a method like the one presented here could be adapted for such organisms by applying what is known about the organism’s lifecycle, genomic features, and mutation rates to the selection of genetic markers and reaction conditions.

In this work, a SNP method for molecular typing of WSMV was developed by adapting protocols designed originally to test human DNA in forensic applications [7, 9, 12, 13, 29]. Fifteen SNPs, located in the CP and HC-Pro genes, were included to generate a genetic fingerprint using the ABI PRISM® SNaPshot™ Multiplex Kit. This procedure provides a mechanism of SNP identification by utilizing only fluorescently labeled dideoxynucleotides (ddNTPs) during the elongation of template-specific primers during thermocycling. The chemical composition of ddNTPs prevents any further elongation, resulting in a product that is one nucleotide longer than the primer itself. This nucleotide represents the SNP element and can be detected by fluorescence scanning of a capillary electrophoretic separation of products. We evaluated the method by determining the assay sensitivity, specificity and reproducibility, as well as its ability to discriminate among known strains and several field isolates of WSMV.


Generation of SNP-specific primers

One hundred and twenty-two SNPs within the CP of WSMV were previously identified [19]. Ten of these SNPs were selected for inclusion in this assay because of their high variability among strains and the presence of a semi-conserved region adjacent to the SNP, providing a suitable template for primer binding. The remaining CP SNPs for this assay were identified by multiple sequence alignment of 85 different WSMV sequences obtained from the National Center for Biotechnology Information (NCBI), including the strains Sidney 81, Type, and El Batan 3 using the same criteria. The SNPs in the HC-Pro region of the genome were identified by performing a multiple sequence alignment of 13 WSMV HC-Pro sequences obtained from the NCBI database. The same criteria were used for selection of the SNPs. Due to the limited number of SNPs having an adjacent template sufficient for primer binding, both the positive and negative sense strands of the viral template were used, an approach previously used [30, 31].

A 20 to 25 nt hybridizing region was designed for each SNP with the addition of a [GACT]n tail to the 5′-end to provide a unique length for each primer, ranging from 30 to 65 nt. Using Oligocalc software [32], the thermodynamics and folding structures were checked for each primer. Each primer was checked for specificity to WSMV using the Basic Local Alignment Search Tool for nucleotides (BLASTn) with the available NCBI database. Any primer that was not specific to WSMV or possessed folding structures with a ∆G greater than 3.0 kcal/mol was eliminated from the study.

The remaining SNP-specific primers were screened for inclusion in the assay by performing monoplex and multiplex in vitro tests using synthetic targets that were specific to each primer. Primers generating a single fluorescent peak in the SNaPshot electropherogram and no peaks in the negative control were selected for inclusion in the WSMV SNaPshot assay (Figure 1). These final fifteen primers are listed in Table 1.

Figure 1

WSMV CP multiplex electropherogram. The height of each peak indicates the concentration of fluorescent units and the location of the peaks along the x-axis indicates the size of the fragment in nucleotides; left to right showing smallest to largest. A) Amplification of synthetic targets using eight SNP-specific primers for the WSMV CP, showing a single amplified target for each primer. B) Negative control (using sterile nanopure water instead of the synthetic templates) for the CP amplification demonstrating the absence of non-target amplification. The peaks present (labeled according to fragment size) are size standards to help with the assignment of the fragment sizes. CP, coat protein; SNP, single nucleotide polymorphism; WSMV, Wheat streak mosaic virus.

Table 1 SNaPshot primers designed for WSMV

Sensitivity and limit of detection

The sensitivities of the CP and HC-Pro assays were determined using the plasmids containing the synthetic target regions. A 1:3 dilution series was created from a starting solution of 0.075 pmole/μl. The SNaPshot reaction was performed and the electropherogram generated by capillary electrophoresis. There was a linear relationship between the fluorescent units and the amount of DNA (pmoles) in each sample (Figure 2).

Figure 2

Sensitivity of SNaPshot assay for the (A) CP and the (B) HC-Pro primers. Both graphs demonstrate a linear relationship for each SNP (individual lines; see Figure 1 for product names). The sensitivity of the assay was determined by identifying the pmoles that correspond to 800 FU using the linear equation of the SNP with the steepest slope (C55 and H45). Linear equations: y = 336285x-475.19 and y = 285041x-1239, CP and HC-Pro, respectively. CP, coat protein; FU, fluorescent units, HC-Pro, helper-component protease, SNP, single nucleotide polymorphism.

The background noise in the electropherograms of the negative and positive controls appeared consistently between 300 and 750 fluorescent units (FUs) (data not shown). To reduce the possibilities of false positives from background peaks, the limit of detection was set at 800 FUs, which was subsequently subtracted from all fluorescent values. The sensitivity of detection was determined by calculating the number of pmoles that correspond to 800 FUs; 3.79 fmoles and 7.15 fmoles for CP and HC-Pro, respectively.

Specificity to WSMV

A panel of near-neighbors to WSMV was tested for SNaPshot assay detection using the same protocols as outlined for WSMV samples. This panel included the phylogenetic near neighbors Oat necrotic mottle virus (ONMV) and Wheat spindle streak mosaic virus (WSSMV), and the environmental near-neighbors Maize dwarf mosaic virus (MDMV), High plains virus (HPV), Barley yellow dwarf virus PAV (BYDV-PAV) and Cereal yellow dwarf virus RPV (CYDV-RPV), as well as the healthy host Triticum aestivum cv. Chisholm. Gel electrophoresis results demonstrated the lack of the WSMV characteristic approximately 1,300 bp amplicon in the near-neighbor samples and the negative controls and the presence of this amplicon in the WSMV positive controls (Figure 3). Amplicons were present in three samples from the CP assay: HPV, CYDV-RPV, and BYDV-PAV. In silico analysis using NCBI BLASTn demonstrated that portions of the C1 and XC1 primers are complementary to regions within each of their genomes. However, the primer sequences are not 100% complementary to any of the three, which may explain the presence of the faint bands. When these three samples were subjected to the SNaPshot assay no fluorescent peaks were present in electropherograms, demonstrating that the SNaPshot assay is specific to WSMV.

Figure 3

Gel electrophoresis of CP and HCPro amplification products from WSMV near-neighbors. Lanes 2–8 and 11–17 are ONMV, WSSMV, MDMV, HPV, CYDV-RPV, BYDV-PAV, and healthy wheat PCR products for CP and HC-Pro assays, respectively. Lanes 1 and 20 are a 100 bp DNA ladder. Lanes 9 and 18 are the positive control CP and HC-PRO WSMV plasmid, respectively. Lanes 10 and 19 are negative controls. Only the positive controls produced approximately 1,300 bp amplicons, corresponding to WSMV. BYDV-PAV, Barley yellow dwarf virus PAV; CP, coat protein; CYDV-RPV, Cereal yellow dwarf virus RPV; HC-Pro, helper-component; HPV, High plains virus; MDMV, Maize dwarf mosaic virus; ONMV, Oat necrotic mottle virus; PCR, polymerase chain reaction; WSMV, Wheat streak mosaic virus; WSSMV, Wheat spindle streak mosaic virus.

Characterization of WSMV strains and field isolates

The discriminatory capability of the SNaPshot assay was tested using known strains of WSMV (Sidney 81, Type, OSU, Merredin 1, Merredin 2, and Ginnindera) and 16 field isolates. A PCR amplicon of the appropriate size for the CP and HC-Pro genes was obtained for each sample (Figure 3). Each strain and most field isolates produced a unique SNaPshot fingerprint representing the genetic variants within each sample (Table 2). However, some field isolates collected from a common geographical region had nearly identical fingerprints, such as three samples from Montana (Conrad 1, Kalispell, and Huntley B).

Table 2 SNaPshot fingerprints of WSMV strains and field isolates using the IUB nucleotide code a

The relationships among these samples were depicted using principal component analysis (PCA) with the inclusion of the in silico SNaPshot results of El Batan 3 and Czech WSMV strains (Figure 4). The strains Czech, El Batan 3, Sidney 81, and Type provide a dimension of diversity as these strains have been thoroughly tested in previous phylogenetic studies [16, 19, 23, 24, 34]. The genetic variability among these samples as depicted in the PCA is consistent with the current understanding of WSMV genetic diversity. The known strains and the field isolates show a high level of diversity with no apparent segregation of field isolates based on geographic location. The results indicate that these samples, even those obtained from Australia, are more closely related to Sidney 81 and Type than to El Batan 3 or Czech, a conclusion that is consistent with previous analyses of the US WSMV population [16, 18, 21, 24].

Figure 4

Principle component analysis of WSMV strains and field isolates. The majority of the samples form a weak cluster between the Sidney 81 and Type strains, with El Batan 3 and Czech strains appearing as outliers. * = profiles generated in silico using NCBI sequence. WSMV, Wheat streak mosaic virus.

Reproducibility of assay

To test for reproducibility of the assay, each SNaPshot profile was generated three times from cDNA of the six WSMV strains. Each replicate produced identical SNaPshot fingerprints for each sample based on qualitative results (data not shown). To determine if the results (fluorescent units) were quantitatively reproducible, analysis of variance (ANOVA) was performed (see materials and methods for complete description of analysis). During the analysis, each nucleotide was compared to the respective nucleotide within each SNP position for each sample replicate. A subsequent analysis was performed using the individual P-values for each replicate group to calculate an overall P-value for each nucleotide. The level of variance for each nucleotide was determined based on these P-values. Based on the results of the peak height-absolute fluorescence dataset analysis, the fluorescence values for all four nucleotides were statistically significant (Table 3). However, only the fluorescence values of nucleotides G and T were found to be statistically significant when the analysis was performed on the peak height-nucleotide percent dataset. The same results were found for both of the area under the curve datasets with the exception of nucleotide G in the area under the curve – nucleotide percent data. The statistical significance indicates that the fluorescence values are different among each strain replicate; however, the difference is minimal when using the area under the curve dataset.

Table 3 ANOVA results of WSMV strain replicates


The ability to discriminate among viral isolates can be useful in evolutionary and ecological studies and for forensic investigations in which determining the origin of a biological agent is critical to criminal attribution. For example, identifying the source of the pathogen became a crucial aspect of the investigation of the anthrax letter attacks in 2001, but because appropriate scientific tools for strain discrimination were not available at that time it was not until seven years later that a suitable test was developed and validated, ultimately leading to attribution [35]. The development and validation of universal, robust, sensitive, and specific detection and discriminatory assays is important for a rapid and efficient response to biological attacks that may occur against our agricultural industry.

The method described here is one example of an assay that addresses this need. During a criminal investigation of a WSMV outbreak, this SNP typing method can be used to generate a genetic profile of the outbreak strain. This profile can then be compared to profiles generated from reference WSMV samples or samples found in a suspect’s possession. The results of the comparison may provide guidance to investigators in determining a likely geographic location of the origin of the virus. The unique patterns in the profile that are known to exist in a specific WSMV strain could help pinpoint the location because certain strains are known to persist in specific areas of the world. Ultimately, the strength in this assay is in comparing one sample to another to determine suspect exclusion for the crime. Determining the threshold of exclusion in terms of number of similar SNPs needs to be determined after further validation testing has occurred.

The SNP assay described here is intended for use in a forensic investigation, but it can be useful also for other molecular studies in which viral variant discrimination may be necessary. Due to the inclusion of 15 SNPs in the assay, we were able to discriminate between multiple strains and isolates of WSMV. The relationship between SNP number and discrimination ability is comparable to that of human DNA typing, in which the use of as many as thirteen short tandem repeat (STR) loci results in an extremely low probability that a genetic fingerprint could be found in any other individual based on allele frequencies in the population [36]. Because none of the probes failed to reveal a variation, each SNP is informative and should be included in the assay. Although we did not test strain El Batan 3 or any of the European strains in vitro, their theoretical profiles were included in the data analysis to demonstrate the estimated discriminatory ability of the method. The fifteen SNPs we used for WSMV provide a high level of confidence that the fingerprint is unique to the variant population without testing a prohibitive number of SNPs.

In our laboratory, the SNP method discriminated among several known strains of the virus and also among field isolates collected throughout the US and in Australia. Several of the profiles are similar due to most of the strains and isolates belonging to the same phylogenetic clade of WSMV [19], but minor differences exist in the secondary peaks that appear for some of the SNPs. These secondary peaks are due to the presence of multiple variants in the population. The ability for this SNP typing assay to detect multiple variants in the population is a feature that adds a higher level of discrimination than can be achieved through consensus sequence analysis. However, to fully understand if all variants or just the dominant variants are represented in the SNP profile, further testing needs to be conducted. In our hands, the limit of detection for the CP assay and HC-Pro assay, respectively, were 3.79 fmoles and 7.15 fmoles of WSMV cDNA, both of which are lower than the detection limit of 10 fmoles reported in the SNaPshot Kit manual. The level of sensitivity of this assay is comparable to that of other SNP-based techniques, reverse-transcription PCR (RT-PCR), and real-time PCR [5].

Including an internal control is important for assay validation. One of the SNaPshot primers (C55) was initially included in the assay to serve as the internal control. This primer is an excerpt of the WSMV diagnostic primer used by National Plant Diagnostic Network laboratories [3739]. Unfortunately, this SNaPshot primer did not appear to be universal for all WSMV isolates tested. Therefore, it was no longer considered the internal control. To validate this assay for forensic investigations, an alternative internal control would need to be identified and included in the assay. At that point, a thorough evaluation of error rate would also be conducted for the assay.

In applications of the SNaPshot assay to field samples it is important to consider that many organisms and viruses other than the target may be encountered, and that some of these, such as Triticum mosaic virus (TriMV) and HPV, which are commonly found together in the same infected tissue [4042], may be genetically related to WSMV. We tested many of the viruses known to co-infect wheat with WSMV. However, since TriMV is almost always found with WSMV [43], we were unable to obtain TriMV infected wheat tissue that was not also infected with WSMV in the necessary time frame; therefore, TriMV was not included in the specificity study. However, the primers used in this assay were confirmed to be WSMV specific through in silico analyses. Our in vitro specificity testing demonstrated that the genetically similar viruses and ecologically similar organisms included in this study would neither interfere with the assay nor provide false positives. This conclusion was confirmed by the lack of fluorescent signal during capillary electrophoresis of the SNaPshot products for HPV, CYDV-RPV and BYDV-PAV.

Implementation of this assay in a forensic investigation requires that the SNP profile of a single sample be reproducible, and this requirement was fulfilled through qualitative comparisons in this study. The inability to apply this assay in a quantitative capacity does not hamper its utility in forensic investigations, as the presence/absence of SNPs in a profile should provide sufficient information. However, quantitative analyses may be beneficial for epidemiological or ecological studies in which the amount of genetic change may provide insight into organismal interactions or disease development.

Even though a wheat virus was used as the model for the development of this assay, SNP typing may prove to be a successful molecular comparison tool for other organisms as well. Ultimately, the strength and reliability of the assay relies on appropriate marker selection which is facilitated by understanding the biology of the organism, the pressures that result in genetic evolution, and the functionally important regions of the genome. SNP analysis may provide additional benefits that other genetic analysis techniques do not. Specifically, DNA does not need to be intact to perform SNP analysis, which would be beneficial for analysis of fungi where the large genomic DNA may be easily sheared during processing.

This study demonstrates a new method for viral isolate comparison that could potentially be utilized in a forensic investigation. Before this method can be incorporated into the forensic science toolbox, a more thorough and exhaustive validation procedure must be conducted. Due to resource constraints, we were unable to perform this level of validation, which would include larger inclusivity and exclusivity panels and standardization of procedures to understand fully the variant representation and error rates. However, this preliminary evaluation of the technology supports the utility of this assay in a forensic setting and suggests a successful validation could be achieved.


This specific, sensitive, and discriminatory SNaPshot assay has been developed using WSMV as a model system. Similar assays could be developed for application to many other plant, animal, or human pathogens by designing organism-specific primers, a task that is becoming ever more feasible with the increasing availability of microbial sequence data and primer designing software. This method already has been shown to be useful for other pathogens, such as Potato virus Y [15], and for genetic comparisons of humans and other organisms [79, 1113, 4447]. The assay uses molecular biology techniques that are well established already in many research and diagnostic laboratories, allowing for a seamless introduction of the technology into existing laboratories.


Selection of virus strains

WSMV strains were received from virus collections at Oklahoma State University in Stillwater, OK, the University of Nebraska in Lincoln, NE, and from AGWEST Plant Laboratories in South Perth, West Australia. Other infected wheat samples were collected in Oklahoma wheat fields and from the Great Plains Diagnostic Network (GPDN) wheat virus survey of 2008. Strains of ONMV, WSSMV, and MDMV strain A were received from a virus collection at Oklahoma State University. Strains of HPV, BYDV-PAV, and CYDV-RPV were received from the GPDN wheat virus survey. All samples were maintained as dry, frozen leaf material at −80°C, with the exception of the Australian WSMV isolate Ginnindera, which was stored at −80°C as an RNA pellet.

SNP identification

SNPs were identified by comparing the sequences of 85 different WSMV sequences obtained from the NCBI database, including the strains Sidney 81 [GenBank: AF057533], Type [GenBank: AF285169.1], and El Batan 3 [GenBank: AF285170.1]. The CP and HC-Pro segments of the genome were targeted for the SNP-typing method because of the high variability among strains [17, 18]. The polymorphic sites selected contained high variability among the strains and were surrounded by semi-conserved regions that would provide a suitable template for primer binding.

Nucleic acid extraction

Total RNA was extracted from healthy or infected wheat (Triticum aestivum) leaves using the RNeasy Plant Mini Kit (Qiagen, Germantown, MD, USA) according to the manufacturer’s instructions. A total of 100 mg of tissue was homogenized at full speed (3,450 oscillations/minute) for 90 seconds using a mini bead-beater (Biospec Products, Bartlesville, OK, USA) and 550 μl of the RLT lysis buffer. The RNA was eluted using 50 μl of RNase-free water supplied with the kit and the concentration determined using a Nanodrop v.2000 spectrophotometer (Thermo Scientific, Waltham, MA, USA).

Reaction controls

A synthetic plasmid carrying the target gene sequences derived from the cloning vector pUC57 was created and produced (GenScript, Piscataway, NJ, USA) for each of the target regions of WSMV to serve as a positive control. The CP synthetic plasmid contains nt 8061–9428 (1,367 bp) and the HC-Pro synthetic plasmid contains nt 1081–2507 (1,426 bp) (based on Sidney 81 genome [GenBank: AF057533]). The plasmids were used during the one-step RT-PCR at a concentration of 2 ng/μl and were processed in the same fashion as the samples.

In addition to the positive control, a negative control was included with each sample set in every process. The negative control consisted of all reaction reagents with sterile nanopure water instead of the viral sample. These controls were used to verify the performance of the assay and confirm the lack of contamination. Additionally, Sidney 81 derived synthetic plasmids confirm the generated profile remains authentic after the sample is subjected to PCR.

cDNA synthesis and amplification

cDNA was synthesized and amplified from total RNA in a single step using the SuperScript™ III One-Step RT-PCR kit with Platinum© Taq (Invitrogen, Carlsbad, CA, USA), following the protocol supplied by the manufacturer. Briefly, in the first part of the RT-PCR, primer C1 (5′-TACTTGACTGGGACCCGAA-3′; Sidney 81 nt 8117 to 8134) was used for the CP reverse-transcription and primer HCR (5′- CATGCTTGTATACTGAGAACAGTCTCTTG-3′; Sidney 81 nt 2373 to 2345) was used for the HC-Pro reverse-transcription, in separate reactions, to synthesize the DNA strand complementary to the RNA sequence [16, 48, 49]. In the second part, the CP region of the cDNA was amplified using the forward primer C1 and the reverse primer XC1 (5′-AACCCACACATAGCTACCAAG-3′; Sidney 81 nt 9371 to 9351) [16, 29, 48, 49]. The HC-Pro region of the cDNA was amplified using the forward primer HCF (5′- GAAATGCACACATGGACTTAGATGGTAT-3′; Sidney 81 nt 1159 to 1186) and the reverse primer HCR [50]. These primers produced approximately 1.3 kb amplicons for both CP and HC-Pro which contained the gene and flanking regions (based on the Sidney 81 genome [GenBank: AF057533]).

Each RT-PCR reaction contained the following components: 25 μl 2X reaction mix (supplied with kit), 1 μg template RNA (6 ng of synthetic plasmid), 0.2 μM of each primer (C1 and XC1, or HCF and HCR), 2 mM MgSO4, 1 μl of RT/Platinum Taq mix (supplied with kit), and sterile nanopure water (17 Ω) up to 50 μl. The kit protocol consisted of cDNA synthesis (one cycle of 50°C for 60 minutes and 94°C for two minutes), amplification (35 cycles of 94°C for 30 seconds, 55°C for 30 seconds, and 68°C for one minute) and elongation (68°C for five minutes).

The amplification of the CP and HC-Pro genes was verified by electrophoresis using a 1.5% agarose gel in 1X TAE buffer containing SYBR® Safe DNA gel stain (Invitrogen). The samples and a 100 bp DNA ladder (Invitrogen) were loaded with 1X loading dye containing bromophenol blue and xylene cyanol and were electrophoresed for 30 minutes at 100 mV. The gel was visualized with a UV-transilluminator and photographed using a camera. The RT-PCR amplification products were treated with a mixture containing one unit of shrimp alkaline phosphatase (SAP), two units of exonuclease 1 (Exo1) (Affymetrix, Inc., Cleveland, OH, USA), and sterile water up to 16 μl for 30 minutes at 37°C followed by 75°C for 15 minutes to inactivate the enzymes. Samples were used immediately in SNaPshot assays or stored at −20°C.

SNaPshot technology

The PRISM® SNaPshot™ Multiplex Kit (Applied Biosystems, Foster City, CA, USA) was used to generate genetic fingerprints from each PCR product. The methods below briefly describe the manufacturer’s process.

Specific WSMV primers were designed for each targeted SNP (Table 1). Each primer was of a different length, ranging from 30 to 65 nt, to facilitate product recognition. The primers were complementary to the DNA template beginning one nucleotide 3′ of the polymorphic site. Controls that were supplied with the SNaPshot kit were included in each assay to verify the performance of the PCR components.

Each reaction contained 5 μl SNaPshot Multiplex Ready Reaction Mix (provided with kit), 3 μl purified PCR product, 1 μl of pooled SNaPshot primers (final concentration of each ranging from 0.1 to 0.4 μM each; see Table 1), and 1 μl of sterile nanopure water (17 Ω). The products were amplified using the following two-step program: 25 cycles of 96°C for 10 seconds and 60°C for 35 seconds. The amplified products were treated with SAP to prevent further binding of ddNTPs. One μl (1 unit) of SAP, 2 μl of 10X SAP buffer, and 7 μl of nuclease-free water were added to each SNaPshot reaction. The mixtures were incubated for 30 minutes at 37°C followed by 15 minutes incubation at 75°C.

The SNaPshot products were separated by capillary electrophoresis on an ABI 3730 Genetic Analyzer (Applied Biosystems) after combining 0.5 μl purified SNaPshot product, 0.5 μl GeneScan 120-LIZ size standard, and 9 μl Hi-Di Formamide (Applied Biosystems). The parameters suggested in the SNaPshot protocol were used. The separated products were analyzed using Peak Scanner Software v1.1 (Applied Biosystems).

Data analysis

Multiple datasets were created from the SNP profiles for quantitative analyses, including data from each peak having a peak height greater than 800 FUs (this limit was set based on background noise, explained in results section). Of the first two datasets created using peak height values, one was composed of the absolute fluorescence values after normalization based on positive control values; the second was composed of percentages of each nucleotide at each SNP position. Two additional datasets were created similarly, but based on the area under the curve instead of peak height. Both peak height and area under the curve were included in the analyses to determine which set of fluorescence values was more reliable for profile comparisons, determined through statistical analyses. Each dataset was analyzed using ANOVA, assuming a randomized complete block model and a P-value <0.05. The strain and the SNP position were the random blocking effects and the replicate was the factor of interest to assess reproducibility.

Using the percent nucleotide dataset from the area under the fluorescent peak values, a PCA was performed using Canoco software (Biometris-Plant Research International, Wageningen, The Netherlands) with square-root transformation [51]. The WSMV strains and field isolates were included in the analysis, in addition to the in silico SNP profile data for strains El Batan 3 and Czech to provide a known dimension of diversity.



Applied Biosystems Inc.


Analysis of variance


Basic local alignment search tool for nucleotides


Base pair


Barley yellow dwarf virus PAV


Coat protein


Cereal yellow dwarf virus RPV


Dideoxynucleoside triphosphates


Exonuclease 1


Fluorescent units


Great Plains Diagnostic Network


Helper component – protease


High plains virus


Maize dwarf mosaic virus


National Center for Biotechnology Information




Oat necrotic mottle virus


Protein 1


Protein 3


Principal component analysis


Polymerase chain reaction


Reverse transcription polymerase chain reaction


Shrimp alkaline phosphatase


Single nucleotide polymorphism


Short tandem repeat




Triticum mosaic virus


United States




Wheat streak mosaic


Wheat streak mosaic virus


Wheat spindle streak mosaic virus.


  1. 1.

    Cupp OS, Walker DE, Hillison J: Agroterrorism in the US: Key security challenge for the 21st century. Biosecur Bioterror. 2004, 2: 97-105. 10.1089/153871304323146397.

    Article  PubMed  Google Scholar 

  2. 2.

    Budowle B, Schutzer SE, Einseln A, Kelley L, Walsh AC, Smith J, Marrone BL, Robertson J, Campos J: Public health: Building microbial forensics as a response to bioterrorism. Science. 2003, 301: 1852-1853. 10.1126/science.1090083.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Budowle B, Burans JP, Breeze R, Wilson M, Chakraborty R: Microbial Forensics. Edited by: Breeze R, Budowle B, Schutzer SE. 2005, Elsevier Academic Press, San Diego, CA, 1-26.

    Google Scholar 

  4. 4.

    Fletcher J, Bender C, Budowle B, Cobb WT, Gold S, Ishimaru C, Luster D, Melcher U, Murch R, Scherm H, Seem RC, Sherwood JL, Sobral BW, Tolin S: Plant pathogen forensics: capabilities, needs, and recommendations. Microbiol Mol Biol Rev. 2006, 70: 450-471. 10.1128/MMBR.00022-05.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  5. 5.

    Henson JM, French R: The polymerase chain reaction and plant disease diagnosis. Ann Rev Phytopath. 1993, 31: 81-109. 10.1146/

    CAS  Article  Google Scholar 

  6. 6.

    Vincelli P, Tisserat N: Nucleic acid-based pathogen detection in applied plant pathology. Plant Dis. 2008, 92: 660-669. 10.1094/PDIS-92-5-0660.

    CAS  Article  Google Scholar 

  7. 7.

    Bouakaze C, Keyser C, Amory S, Crubezy E, Ludes B: First successful assay of Y-SNP typing by SNaPshot minisequencing on ancient DNA. Int J Legal Med. 2007, 121: 493-499. 10.1007/s00414-007-0177-3.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Bujalkova M, Zavodna K, Krivulcik T, Ilencikova D, Wolf B, Kovac M, Karner-Hanusch J, Heinimann K, Marra G, Jiricny J, Bartosova Z: Multiplex SNaPshot genotyping for detecting loss of heterozygosity in the mismatch-repair genesMLH1andMSH2in microsatellite-unstable tumors. Clin Chem. 2008, 54: 1844-1854. 10.1373/clinchem.2008.108902.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Grignani P, Peloso G, Achilli A, Turchi C, Tagliabracci A, Alu M, Beduschi G, Ricci U, Giunti L, Robino C, Gino S, Previdere C: Subtyping mtDNA haplogroup H by SNaPshot minisequencing and its application in forensic individual identification. Int J Legal Med. 2006, 120: 151-156. 10.1007/s00414-005-0059-5.

    CAS  Article  PubMed  Google Scholar 

  10. 10.

    Jungerius BJ, Veenendaal AT, van Oost BA, te Pas MF, Groenen MA: Typing single-nucleotide polymorphisms using a gel-based sequencer. Mol Biotechnol. 2003, 25: 283-287. 10.1385/MB:25:3:283.

    CAS  Article  PubMed  Google Scholar 

  11. 11.

    Kawuki RS, Ferguson M, Labuschagne M, Herselman L, Kim DJ: Identification, characterisation and application of single nucleotide polymorphisms for diversity assessment in cassava (Manihot esculentaCrantz). Mol Breeding. 2009, 23: 669-684. 10.1007/s11032-009-9264-0.

    CAS  Article  Google Scholar 

  12. 12.

    Kohnemann S, Sibbing U, Pfeiffer H, Hohoff C: A rapid mtDNA assay of 22 SNPs in one multiplex reaction increases the power of forensic testing in European Caucasians. Int J Legal Med. 2008, 122: 517-523. 10.1007/s00414-008-0267-x.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Quintans B, Alvarez-Iglesias V, Salas A, Phillips C, Lareu M, Carracedo A: Typing of mitochondrial DNA coding region SNPs of forensic and anthropological interest using SNaPshot minisequencing. Forensic Sci Int. 2004, 140: 251-257. 10.1016/j.forsciint.2003.12.005.

    CAS  Article  PubMed  Google Scholar 

  14. 14.

    Ranade K, Chang MS, Ting CT, Pei D, Hsiao CF, Olivier M, Pesich R, Hebert J, Chen YDI, Dzau VJ, Curb D, Olshen R, Risch N, Cox DR, Botstein D: High-throughput genotyping with single nucleotide polymorphisms. Genome Res. 2001, 11: 1262-1268.

    PubMed Central  CAS  PubMed  Google Scholar 

  15. 15.

    Rolland M, Glais L, Kerlan C, Jacquot E: A multiple single nucleotide polymorphisms interrogation assay for reliablePotato virusY group and variant characterization. J Virol Methods. 2008, 147: 108-117. 10.1016/j.jviromet.2007.08.022.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    French R, Stenger DC: Population structure within lineages ofWheat streak mosaic virusderived from a common founding event exhibits stachastic variation inconsistent with the deterministic quasi-species model. Virology. 2005, 343: 179-189. 10.1016/j.virol.2005.08.037.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Shukla D, Frenkel M, Ward C: Structure and function of the potyvirus genome with special reference to the coat protein coding region. Can J Plant Pathol. 1991, 13: 178-191. 10.1080/07060669109500953.

    CAS  Article  Google Scholar 

  18. 18.

    Choi IR, Hall JS, Henry M, Zhang L, Hein GL, French R, Stenger DC: Contributions of genetic drift and negative selection on the evolution of three strains ofWheat streak mosaic tritimovirus. Arch Virol. 2001, 146: 619-628. 10.1007/s007050170167.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Stenger DC, Seifers DL, French R: Patterns of polymorphism inWheat streak mosaic virus: sequence space explored by a clade of closely related viral genotypes rivals that between the most divergent strains. Virology. 2002, 302: 58-70. 10.1006/viro.2001.1569.

    CAS  Article  PubMed  Google Scholar 

  20. 20.

    Niblett CL, Zagula KR, Calvert LA, Kendall TL, Stark DM, Smith CE, Beachy RN, Lommel SA: cDNA cloning and nucleotide sequence of theWheat streak mosaic viruscapsid protein gene. J Gen Virol. 1991, 72: 499-504. 10.1099/0022-1317-72-3-499.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    French R, Stenger DC: Evolution ofWheat streak mosaic virus: Dynamics of population growth within plants may explain limited variation. Annu Rev Phytopathol. 2003, 41: 199-214. 10.1146/annurev.phyto.41.052002.095559.

    CAS  Article  PubMed  Google Scholar 

  22. 22.

    Hall JS, French R, Hein GL, Morris TJ, Stenger DC: Three distinct mechanisms facilitate genetic isolation of sympatricWheat streak mosaic viruslineages. Virology. 2001, 282: 230-236. 10.1006/viro.2001.0841.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Hall JS, French R, Morris TJ, Stenger DC: Structure and temporal dynamics of populations withinWheat streak mosaic virusisolates. J Virol. 2001, 75: 10231-10243. 10.1128/JVI.75.21.10231-10243.2001.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  24. 24.

    Rabenstein F, Seifers DL, Schubert J, French R, Stenger DC: Phylogenetic relationships, strain diversity, and biogeography of tritimoviruses. J Gen Virol. 2002, 83: 895-906.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Berger PH, Shiel PJ: Potyvirus isolation and RNA extraction. Method Mol Biol. 1998, 81: 151-160.

    CAS  Google Scholar 

  26. 26.

    Chen J, Chen J, Adams M: A universal PCR primer to detect members of thePotyviridaeand its use to examine the taxonomic status of several members of the family. Arch Virol. 2001, 146: 757-766. 10.1007/s007050170144.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    French R, Robertson NL: Simplified sample preparation for detection ofWheat streak mosaic virusandBarley yellow dwarf virusby PCR. J Virol Methods. 1994, 49: 93-100. 10.1016/0166-0934(94)90059-0.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Agrios GE: Plant pathology. 2005, Elsevier Academic Press, Oxford, 5

    Google Scholar 

  29. 29.

    Carver J: Molecular typing of Wheat streak mosaic virus for forensic applications. MS thesis. 2007, Oklahoma State University-Center for Health Sciences, Department of Forensic Science,

    Google Scholar 

  30. 30.

    Filippini S, Blanco A, Fernández-Marmiesse A, Àlvarez-Iglesias V, RuÍz-Ponte C, Carracedo A, Vega A: Multiplex SNaPshot for detection ofBRCA 1/2common mutations in Spanish and Spanish related breast/ovarian cancer families. BMC Med Genet. 2007, 8: 40-

    PubMed Central  Article  PubMed  Google Scholar 

  31. 31.

    Lurkin R, Stoehr R, Hurst CD, van Tilborg AAG, Knowles MA, Hartmann A, Zwarthoff EC: Two multiplex assays that simultaneously identify 22 possible mutation sites in theKRAS, BRAF, NRASandPIK3CAgenes. PLoS One. 2010, 5: e8802-10.1371/journal.pone.0008802.

    PubMed Central  Article  PubMed  Google Scholar 

  32. 32.

    Kibbe WA: OligoCalc: an online oligonucleotide properties calculator. Nucleic Acids Res. 2007, 35: W43-W46. 10.1093/nar/gkm234.

    PubMed Central  Article  PubMed  Google Scholar 

  33. 33.

    Nomenclature Committee of the International Union of Biochemistry (NC-IUB): Nomenclature for incompletely specified bases in nucleic acid sequences: recommendations 1984. Proc Natl Acad Sci U S A. 1986, 83: 4-8. 10.1073/pnas.83.1.4.

    Article  Google Scholar 

  34. 34.

    Stenger DC, Hall JS, Choi IR, French R: Phylogenetic relationships within the familyPotyviridae:Wheat streak mosaic virusandBrome streak mosaic virusare not members of the genusRymovirus. Phytopathology. 1998, 88: 782-787. 10.1094/PHYTO.1998.88.8.782.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    Budowle B, Johnson MD, Fraser CM, Leighton TJ, Murch RS, Chakraborty R: Genetic analysis and attribution of microbial forensics evidence. Crit Rev Microbiol. 2005, 31: 233-254. 10.1080/10408410500304082.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Budowle B, Shea B, Niezgoda S, Chakraborty R: CODIS STR loci data from 41 sample populations. J Forensic Sci. 2001, 46: 453-489.

    CAS  PubMed  Google Scholar 

  37. 37.

    Stack J, Hammerschmidt R, Hudler G, Luke E, Bostock R: National Plant Diagnostic Network: a record of accomplishment. CSREES five-year review. 2007, Retrieved from

    Google Scholar 

  38. 38.

    Stack J, Cardwell K, Hammerschmidt R, Byrne J, Loria R, Snover-Clift K, Baldwin W, Wisler G, Beck H, Bostock RM, Thomas CS, Luke E: The National Plant Diagnostic Network. Plant Dis. 2006, 90: 128-136. 10.1094/PD-90-0128.

    Article  Google Scholar 

  39. 39.

    Myslik JT, Nassuth A: Rapid detection of viruses, transgenes, and mRNAs in small plant leaf samples. Plant Mol Biol Reporter. 2001, 19: 329-340. 10.1007/BF02772831.

    CAS  Article  Google Scholar 

  40. 40.

    Brakke M: Virus disease in wheat. Wheat and Wheat Improvement, second edition. Edited by: Heyne E. 1987, American Society of Agronomy, Crop Science Society of America, Soil Science Society of America, Madison, WI, 585-603.

    Google Scholar 

  41. 41.

    Jensen SG, Lane LC, Seifers DL: A new disease of maize and wheat in the high plains. Plant Dis. 1996, 80: 1387-1390. 10.1094/PD-80-1387.

    Article  Google Scholar 

  42. 42.

    Tatineni S, Ziems AD, Wegulo SN, French R: Triticum mosaic virus: A distinct member of the FamilyPotyviridaewith an unusually long leader sequence. Phytopathology. 2009, 99: 943-950. 10.1094/PHYTO-99-8-0943.

    CAS  Article  PubMed  Google Scholar 

  43. 43.

    Burrows ME, Franc G, Rush CM, Blunt T, Ito D, Olson J, Kinzer K, O’Mara J, Price JA, Tande C, Ziems A, Stack J: Occurence of viruses in wheat in the Great Plains region. Plant Health Progress. 2009, 10:1094/PHP-1009-0706-01-RS.

    Google Scholar 

  44. 44.

    Ceriotti G, Chessa S, Bolla P, Budelli E, Bianchi L, Duranti E, Caroli A: Single nucleotide polymorphisms in the ovine casein genes detected by polymerase chain reaction-single stranded conformation polymorphism. J Dairy Sci. 2004, 87: 2606-2613. 10.3168/jds.S0022-0302(04)73386-X.

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Civaova K, Knoll A: SNaPshot minisequencing and a panel of candidate genes for complex routine testing of meat performance traits in pigs. Anim Biotechnol. 2007, 18: 109-115. 10.1080/10495390601088657.

    Article  Google Scholar 

  46. 46.

    Murphy SJ, Hughes AE, Patterson CC, Anderson LA, Watson RGP, Johnston BT, Comber H, McGuigan J, Reynolds JV, Murray LJ: A population-based association study of SNPs ofGSTP1, MnSOD, GPX2and Barrett’s esophagus and esophageal adenocarcinoma. Carcinogenesis. 2007, 28: 1323-1328. 10.1093/carcin/bgm007.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Ovesna J, Polakova K, Kucera L, Rulcova J: SNP typing in cereals: comparison of SSCP and SNaPshot markers using BarleyMlolocus as a model. Czech J Genet Plant Breed. 2003, 39: 109-112.

    Google Scholar 

  48. 48.

    Chenault KD, Hunger RM, Sherwood JL: Comparison of the nucleotide sequences of the coat protein open reading frame of nine isolates ofWheat streak mosaic rymovirus. Virus Genes. 1996, 13: 187-188. 10.1007/BF00366978.

    CAS  Article  PubMed  Google Scholar 

  49. 49.

    McNeil JE, French R, Hein GL, Baenziger PS, Eskridge KM: Characterization of genetic variability among natural populations ofWheat streak mosaic virus. Phytopathology. 1996, 86: 1222-1227. 10.1094/Phyto-86-1222.

    CAS  Article  Google Scholar 

  50. 50.

    Stenger DC, French R: Functional replacement ofWheat streak mosaic virusHC-Pro with the corresponding cistron from a diverse array of viruses in the family Potyviridae. Virology. 2004, 323: 257-267. 10.1016/j.virol.2004.03.014.

    CAS  Article  PubMed  Google Scholar 

  51. 51.

    Lepš J, Šmilauer P: Multivariate Analysis of Ecological Data using Canoco. 2003, Cambridge University Press, Cambridge, UK

    Google Scholar 

Download references


The authors acknowledge Jacob Price from Texas AgriLife Research Center, Mary Burrows from the Montana State University, Roy French from USDA-ARS, Geoffrey Dwyer from AGWEST plant laboratories in South Perth, Australia, TeeCie Brown from Oklahoma State University, and Trenna Bladgen from Oklahoma State University for supplying virus samples; Mike Palmer from Oklahoma State University and Andrew Doust from Oklahoma State University for their assistance with data analysis; and Mohammed Arif and Peter Hoyt from Oklahoma State University for proof reading the article. This work was approved for publication by the Director of the Oklahoma Agricultural Experiment Station and supported in part under project OKL 2052, and in part by USDA NIFA.

Author information



Corresponding author

Correspondence to Stephanie M Rogers.

Additional information

Competing interests

The authors declare they have no competing interests.

Authors’ contributions

SMR participated in the design of the study and carried out the experimental design, genotyping, data analysis, and drafted the manuscript. MP participated in the design of the study and performed the statistical analysis. RWA conceived of the study and participated in the design of study and data analysis. UM participated in design of the study and data analysis. JC participated in the concept and design of the study and genotyping. JF participated in conception and design of study, data analysis, and drafting of manuscript. All authors read and reviewed the final manuscript.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. 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 work is properly cited.

Reprints and Permissions

About this article

Cite this article

Rogers, S.M., Payton, M., Allen, R.W. et al. Method: a single nucleotide polymorphism genotyping method for Wheat streak mosaic virus. Investig Genet 3, 10 (2012).

Download citation


  • Coat Protein
  • Field Isolate
  • Forensic Investigation
  • Wheat Streak Mosaic Virus
  • Maize Dwarf Mosaic Virus