High throughput, high resolution selection of polymorphic microsatellite loci for multiplex analysis
Plant Methodsvolume 1, Article number: 3 (2005)
Large-scale genetic profiling, mapping and genetic association studies require access to a series of well-characterised and polymorphic microsatellite markers with distinct and broad allele ranges. Selection of complementary microsatellite markers with non-overlapping allele ranges has historically proved to be a bottleneck in the development of multiplex microsatellite assays. The characterisation process for each microsatellite locus can be laborious and costly given the need for numerous, locus-specific fluorescent primers.
Here, we describe a simple and inexpensive approach to select useful microsatellite markers. The system is based on the pooling of multiple unlabelled PCR amplicons and their subsequent ligation into a standard cloning vector. A second round of amplification utilising generic labelled primers targeting the vector and unlabelled locus-specific primers targeting the microsatellite flanking region yield allelic profiles that are representative of all individuals contained within the pool. Suitability of various DNA pool sizes was then tested for this purpose. DNA template pools containing between 8 and 96 individuals were assessed for the determination of allele ranges of individual microsatellite markers across a broad population. This helped resolve the balance between using pools that are large enough to allow the detection of many alleles against the risk of including too many individuals in a pool such that rare alleles are over-diluted and so do not appear in the pooled microsatellite profile. Pools of DNA from 12 individuals allowed the reliable detection of all alleles present in the pool.
The use of generic vector-specific fluorescent primers and unlabelled locus-specific primers provides a high resolution, rapid and inexpensive approach for the selection of highly polymorphic microsatellite loci that possess non-overlapping allele ranges for use in large-scale multiplex assays.
Microsatellite analysis using fluorescently labelled primers and capillary fractionation is the pre-eminent method for the genetic analysis of eukaryotic organisms. The approach is routinely used for many applications including forensic analysis , linkage mapping and association genetics , population genetics [3, 4] and genetic analysis of diversity . The need to screen microsatellite loci for polymorphism between genotypes within the target organism, and for their suitability in multiplex analysis, is an inevitable part of such efforts. The high cost of fluorescently labelled primers has meant that selection of microsatellite markers has typically relied on initial, low-resolution screens of unlabelled primers prior to high-resolution marker selection using labelled primers. However, the preliminary screen is inevitably crude and inefficient, making it either prone to error or reliant upon the more expensive high-resolution selection. There is therefore a need for a high throughput and high-resolution single-step method of selecting appropriate microsatellite markers for genetic studies .
For multiplex analysis, greatest efficiency is achieved when utilising many polymorphic loci possessing closely spaced, non-overlapping allelic ranges. Unexpected allelic range overlap between multiplexed microsatellite loci yields ambiguous alleles that may be misassigned to an inappropriate locus, compromising the integrity of the data set. One inevitable problem lies in the possibility that the screen does not encompass all alleles present in the population under study. Confidence in the definition of allelic ranges is invariably a function of the number and diversity of genotypes screened. There is therefore a balance between the desire to examine many individuals and the cost of doing so using fluorescently labelled primers. Thus, screening invariably becomes expensive as the number of genotypes tested grows, and as the number of discarded markers increases. Common approaches to selecting microsatellite markers for multiplex use include assembling panels from previously fluorescently characterised individual markers [7, 8], and pre-screening markers on polyacrylamide gels utilising radioactivite labelling of PCR products . Several authors have proposed low cost alternatives for preliminary screens using direct DNA staining following polyacrylamide gel electrophoresis [10–13]. Such strategies have merit, but are labour-intensive, cannot assign actual size ranges and generally lack the resolution required to accurately predict polymorphism in dinucleotide markers [10–13]. One methodology that is able to generate high resolution allelic ladders in a similar fashion to the method reported here is that of Oetting . This method employs the use of locus specific primers tailed with generic sequence allowing a second round of labelled PCR and subsequent capillary fractionation. This method however suffers from a number of potential disadvantages relative to our method. The use of long oligonucleotides for PCR of genomic templates at below optimal annealing temperature allows for an increased frequency in the production non specific amplification products. The PCR amplification conditions required for locus specific amplification using tailed oligonucleotides are often different to those conditions optimal for amplification with equal length 20 mer oligonucleotides. The method of Oetting is not suited to the genotypic analysis of dinucleotide repeat markers due the possibility of extensive stutter profiles generated by the second round of PCR complicating the allelic profiles and so is only considered of merit for marker selection.
Here, we propose a simple but novel approach in which microsatellite amplicons generated from pooled genomic DNA templates are ligated into a standard cloning vector, re-amplified using a labelled universal primer targeting the plasmid insert flanking region, and an unlabelled locus-specific primer. The resultant profiles represents allelic ladders derived from the component alleles contained by the pooled DNA. This procedure thereby offers a single assay, high resolution and inexpensive means of screening microsatellite loci for polymorphism and allelic size range. The profile also offers a qualitative indication of the locus with regard to stutter, a problem often associated with the use of dinucleotide repeat markers for genetic analysis, but also of interest to laboratories utilising tri- and tetranucletide repeat markers.
When employing a pooling strategy, there is a balance between sampling extensively to encompass the full range of variation, and dilution of individuals within the pool such that rare alleles are not detected. DNA pools were created by combining equal amounts of individual DNA samples before dilution with nano-pure water to 5 ng·μL-1. To select the most appropriate pool size, while allowing the detection of rare alleles, DNA pools of 8, 12, 16, 24, 32, and 48 individuals were compared. PCR amplification of pooled DNA utilising unlabelled primer pairs specific to single microsatellite loci  were performed incorporating 5 ng template DNA with AccuPrime Taq DNA Polymerase in supermix I, using half recommended volumes (Invitrogen Ltd). The thermal cycling protocol was 96°C for 2 min; 35 cycles of 96°C for 30 s, 51°C or 46°C for 30 s dependent on primer annealing characteristics, 72°C for 2 min; followed by 72°C for 10 min, in a MJ Research PTC-100 thermal cycler (Genetic Research Instrumentation Ltd). Successful PCR was confirmed by 1.5% (w/v) agarose gel electrophoresis . PCR products for multiple individual microsatellite loci, amplified from aliquots of the same template DNA pool, were combined then purified using NucleoFast 96 PCR cleanup plates (Macherey-Nagel GmbH & Co. KG), before ligation into pDrive vector (Qiagen Ltd). Ligation products were diluted 1/10 with HPLC grade water and used as template for a second round of PCR using the 'reverse' microsatellite specific primer and a generic fluorescently labelled primer, M13 (-40), targeting the plasmid. Labelled amplicons were diluted 1/100 in HPLC grade water and fractionated by capillary electrophoresis on an ABI 3100 and viewed with genotyper 3.7 software (Applied Biosystems UK Ltd). The allele size reported by this method is that expected from the microsatellite primers plus an additional 150 bases of vector sequence. Comparison of profiles from pooled amplifications and those of constituent members of the pools demonstrated that homozygous individuals possessing a rare allele could be detected reliably in pools of 12 individuals or less. Thus, one strategy would be to assemble several small pools, allowing variance in allelic limits to be described, and continue screening until the addition of more pools no longer increases the allele range. In practice, however, it may be preferable to use much larger pools and accommodate for uncertainty over rare alleles by imposing buffer zones around detected allelic ranges prior to multiplexing. We empirically tested this approach. Template DNA from 96 diverse genotypes of Theobroma cacao was adjusted to 5 ng·uL-1, pooled and individually amplified by PCR for 84 dinucleotide cocoa microsatellite markers described by Pugh et al . The complex profiles generated (Figure 1) broadly represent the array of alleles present when genotypes were assayed individually. We therefore selected 36 markers generating the widest range of homogeneous peaks for further study. The allelic range of peak sizes is taken to be indicative of the allelic range in the unsampled gene pool. In general, loci generating large numbers of peaks with approximately even height (Figure 1B,C) were highly informative whereas those producing few peaks (Figure 1A) or profiles dominated by one peak (Figure 1D) have less utility for genetic analysis.
We then examined the relationship between predicted allele range in pooled profiles and that observed after wider genotype sampling. For this, we employed multiplex PCR microsatellite analyses performed on 672 individual cocoa genotypes (Table 1). Two loci predicted to yield few alleles on the basis of the sample pool (mTcCIR080 and mTcCIR155) produced the same number of alleles when the sample range was expanded to include 672 individual cocoa genotypes. However, the number of alleles in more variable loci (mTcCIR131 and mTcCIR190) increased from 9 to 12 when the sample range was expanded, with the size range increasing by 67% and 62% respectively. Given a modest increase in allelic range when sample size was increased six fold, one approach would be to accommodate undetected alleles by imposing a buffer between the ranges of neighbouring loci prior to multiplexing. In this case, a spacing of 1× predicted range either side of the mean allele size would appear adequate. Overall, adoption of this protocol allows for improved selection of compatible polymorphic microsatellite markers, with reduced likelihood of producing overlapping profiles in multiplexed microsatellite reactions.
Adoption of this methodology allows for both a qualitative and semi quantitative characterisation of polymorphism at individual microsatellite loci. When using pooled samples, combining DNA from up to 12 individuals allowed for the reliable detection of single copy alleles within that sample. If characterising microsatellite loci using DNA pools of greater than 12 individuals the incorporation of a buffer zone into the final genotyping assay, based on the observed range of allele sizes, can allow for the variability likely to be encountered in a larger sample size. The methodology is suitable for high throughput applications by the combination of differing fluorescent dyes in association with convenient liquid handling formats. This protocol benefits from initially utilising unlabelled primers identical to those used in the final genotyping assay, reducing the possibilities of unexpected banding patterns due to changes in primer sequence or assay conditions. The high resolution DNA size measurement makes this protocol suitable for characterising dinucleotide microsatellite loci.
Polymerase chain reaction
Carey L, Mitnik L: Trends in DNA forensic analysis. Electrophoresis. 2002, 23: 1386-1397. 10.1002/1522-2683(200205)23:10<1386::AID-ELPS1386>3.0.CO;2-M.
Laan M, Paabo S: Demographic history and linkage disequilibrium in human populations. Nat Genet. 1997, 17: 435-438. 10.1038/ng1297-435.
PerezLezaun A, Calafell F, Seielstad M, Mateu E, Comas D, Bosch E, Bertranpetit J: Population genetics of y-chromosome short tandem repeats in humans. J Mol Evol. 1997, 45: 265-270.
PerezLezaun A, Calafell F, Mateu E, Comas D, Ruiz-Pacheco R, Bertranpetit J: Microsatellite variation and the differentiation of modern humans. Hum Genet. 1997, 1: 1-7.
Tsutsui ND, Suarez AV, Holway DA, Case TJ: Reduced genetic variation and the success of an invasive species. Proc Natl Acad Sci USA. 2000, 97: 5948-5953. 10.1073/pnas.100110397.
Frasier TR, Wilson PJ, White BN: Rapid screening of microsatellite markers for polymorphisms using SYBR® green 1 and a DNA sequencer. BioTechniques. 2004, 36: 408-409.
Narvel JM, Chu WC, Fehr WR, Cregan PB, Shoemaker RC: Development of multiplex sets of simple sequence repeat DNA markers covering the soybean genome. Molecular Breeding. 2000, 6: 175-183. 10.1023/A:1009637119947.
Tang S, Kishore VK, Knapp SJ: PCR-multiplexes for a genome-wide framework of simple sequence repeat marker loci in cultivated sunflower. Theor Appl Genet. 2003, 107: 6-19.
Tommasini L, Batley J, Arnold GM, Cooke RJ, Donini P, Law JR, Lowe C, Moule C, Trick M, Edwards KJ: The development of multiplex simple sequence repeat (SSR) markers to complement distinctness, uniformity and stability testing of rape (Brassica napus L.) varieties. Theor Apl Genet. 2003, 106: 1091-1101.
Morin PA, Smith DG: Non-radioactive detection of hypervariable simple sequence repeats in short polyacrylamide gels. BioTechniques. 1995, 19: 223-228.
Scrimshaw BJ: Non-radioactive detection of hypervariable simple sequence repeats in short polyacrylamide gels. BioTechniques. 1992, 13: 188-
White HW, Kusukawa N: Agarose-based system for separation of short tandem repeat loci. BioTechniques. 1997, 22: 976-980.
Houriham RN, O'Sullivan GC, Morgan JG: High-resolution detection of loss of heterozygosity of dinucleotide microsatellite markers. BioTechniques. 2001, 30: 342-346.
Oetting WS, Lee HK, Flanders DJ, Wiesner GL, Sellers TA, King RA: Linkage analysis with multiplexed short tandem repeat polymorphisms using infrared fluorescence and M13 tailed primers. Genomics. 1995, 30: 450-458. 10.1006/geno.1995.1264.
Pugh T, Fouet O, Risterucci AM, Brottier P, Abouladze M, Deletrez C, Courtois B, Clement D, Larmande P, N'Goran JAK, Lanaud C: A new cacao linkage map based on codominant markers: development and integration of 201 new microsatellite markers. Theor App Genet. 2004, 108: 1151-1161. 10.1007/s00122-003-1533-4.
Sambrook J, Fitsch EF, Maniatis T: Molecular Cloning: A Laboratory Manual. 1989, Cold Spring Harbor, Cold Spring Harbor Press
We thank Olivier Sounigo, and Claire Lanaud for providing DNA, Didier Clement for providing as then unpublished microsatellite information (Pugh et al., 2004) and Steve Brown for providing a linkage map of cocoa microsatellite markers. This work forms part of a collaboration between The University of Reading and the Cocoa Research Unit, The University of West Indies, Trinidad and Tobago, funded by The Biscuit, Cake, Chocolate and Confectionery Association, United Kingdom.
The author(s) declare that they have no competing interests.
NCC conceived and developed the microsatellite selection protocol and drafted the manuscript. The overall project was conceived by DRB. MJW aided the experimental design and played a major role in developing the manuscript. All authors have read and approved the final manuscript.