Development and characterization of 96 microsatellite markers suitable for QTL mapping and accession control in an Arabidopsis core collection
Plant Methodsvolume 10, Article number: 2 (2014)
To identify plant genes involved in various key traits, QTL mapping is a powerful approach. This approach is based on the use of mapped molecular markers to identify genomic regions controlling quantitative traits followed by a fine mapping and eventually positional cloning of candidate genes. Mapping technologies using SNP markers are still rather expensive and not feasible in every laboratory. In contrast, microsatellite (also called SSR for Simple Sequence Repeat) markers are technologically less demanding and less costly for any laboratory interested in genetic mapping.
In this study, we present the development and the characterization of a panel of 96 highly polymorphic SSR markers along the Arabidopsis thaliana genome allowing QTL mapping among accessions of the Versailles 24 core collection that covers a high percentage of the A. thaliana genetic diversity. These markers can be used for any QTL mapping analysis involving any of these accessions. We optimized the use of these markers in order to reveal polymorphism using standard PCR conditions and agarose gel electrophoresis. In addition, we showed that the use of only three of these markers allows differentiating all 24 accessions which makes this set of markers a powerful tool to control accession identity or any cross between any of these accessions.
The set of SSR markers developed in this study provides a simple and efficient tool for any laboratory focusing on QTL mapping in A. thaliana and a simple means to control seed stock or crosses between accessions.
Identifying plant genes that control key traits is central to understanding the mechanism of their gene action and to target them in breeding programs. Few genes controlling quantitative traits have been identified in plants and most have been cloned by linkage mapping. Mapping of Quantitative Trait Loci (QTL) is a standard procedure that requires both genotypic (i.e. based on molecular markers) and phenotypic data generated from a segregating population[2, 3]. In a first step, mapped QTL intervals of several tens of centimorgans are identified for which a portion of the phenotypic variation of trait is explained. A fine mapping approach that aims to narrow down these QTL to candidate genes is then implemented.
This initial step necessitates having a set of molecular markers distributed along the plant genome and polymorphic between the two parents of the progeny. In Arabidopsis thaliana, major inroads have been made for gene mapping, through the identification of thousands of Single Nucleotide Polymorphism (SNP) markers[4–11] or through extensive genome re-sequencing of individual accessions[12–14]. However, mapping technologies using SNPs markers are still rather expensive and not applicable in every laboratory. On the other hand, due to their relative abundance in the genome and simple cost effective detection microsatellite markers (also named simple sequence repeat or SSR markers) remain for many laboratories the markers of choice in a primary QTL mapping screen. Microsatellites are tandem repeat motifs of 1 to 6 bp that have frequent occurrence in all genomes analyzed to date[15, 16]. The high variability of microsatellites is mainly due to different number of repeats in the region of the repeat motif but also to short insertion/deletion events. Thus, SSR markers usually exhibit many alleles in comparison to SNP or other polymerase chain reaction (PCR)-based markers. Moreover, they have the advantage of being codominant. Consequently, only a limited number of markers are sufficient to detect polymorphism among parental lines. As the flanking regions of the repeat motif are in many cases highly conserved (particularly in the case of EST derived SSRs), microsatellite markers are easily and reproducibly amplified by PCR in most accessions and visualized by gel electrophoresis (agarose, polyacrylamide) or automated sequencing. In addition, their identification is relatively simple when plant genome sequences are available. In addition, in the case of A. thaliana, several studies showed that microsatellites are abundant and highly variable[18–20].
In this paper, we describe the development and the characterization of a panel of 96 SSR markers distributed along the A. thaliana genome suitable for QTL mapping in an Arabidopsis core collection.
Results and discussion
The primary aim of this study was to produce molecular markers suitable for QTL mapping in A. thaliana without the use of expensive and specific technical platforms. Therefore, we decided to develop a panel of polymorphic SSR markers among accessions of the Versailles core collection 24 which maximizes genetic diversity in a limited number of individuals and has been used to generate F2 and Recombinant Inbred Lines (RIL)[21, 22]. We showed previously in several studies that these accessions display a diversity in behavior toward plant viruses allowing development of mapping projects to identify plant genes involved in plant/virus interactions[23, 24]. Our strategy was to identify markers distributed along the Arabidopsis genome, approximately every 4 Mb (approximately 20 cM), which are as polymorphic as possible among the Versailles core-24 accessions. Those criteria optimize their use for mapping analysis within progenies obtained from crosses between any of the 24 Versailles accessions. Six other accessions (Col-0, C24, Ws-2, Nd-1, Bay-0 and Ler-2) were added to the study as they are frequently used as parents in our studies[23–27]. The criteria of 4-Mb intervals along each chromosome generated 8 regions (1.1 to 1.8) in chromosome 1 (chr1), 5 regions (2.1 to 2.5) in chr2, 6 regions (3.1 to 3.6) in chr3, 5 regions (4.1 to 4.5) in chr4 and 7 regions (5.1 to 5.7) in chr5 (Figure 1) in which we searched for tandem repeat motifs. Thirty-one SSR markers distributed throughout these regions were found in the Versailles MSAT database or from The Arabidopsis Information Resource (TAIR). They displayed high genetic diversity (Additional file1: Table S1) but they were not sufficient to cover the 31 genomic regions of interest (Figure 1). Moreover, they are not sufficiently polymorphic among our set of accessions (Additional file1: Table S2). To complete this marker set, new loci containing tandem repeat motifs were identified within the Col-0 sequence as described in the Methods section. We preferentially selected perfect dinucleotide repeat motifs with a minimum of ten repeats except for a few regions where we had to select mono or trinucleotide repeat motifs or imperfect motifs (Additional file1: Table S1). Primer sets were developed to amplify these repeat containing loci and tested among the 30 accessions. Our aim was to select SSRs in each region that displayed at least one polymorphism among all accessions and also for which polymorphism can be preferentially detected on agarose rather than polyacrylamide gels. Following this procedure, we generated 65 new SSRs, named BSATX.Y, X being the chromosome number and Y, the development rank number, with two to four markers generated in each selected chromosomal region. They are described in Additional file1: Table S1 and their relative position on the five Arabidopsis chromosomes is presented in Figure 1. For some regions, especially those located around the centromere, we extended the size of the interval in order to identify suitable SSRs (Figure 1). Indeed, it has been previously shown that the occurrence of microsatellites is consistently lower in the centromeric regions than in the rest of the chromosomes. To facilitate the use of these markers, primers were designed such that almost all SSRs can be amplified using the same PCR conditions. An example of SSR polymorphism displayed on an agarose gel is presented in Figure 2. The size of each SSR allele for each accession is presented in Additional file1: Table S2. A few SSRs in some accessions could not be amplified at all, most likely due to polymorphisms or deletions in primer binding sites compared with the Col-0 sequence used for primer design. In other cases two amplified products were produced, most likely due to residual heterozygosity (Additional file1: Table S2). From the size polymorphism of each SSR, several features were deduced. The average number of alleles for each marker is 14, ranging from 2 to 27 alleles per marker (Additional file1: Table S1). The average polymorphism information content (PIC, which corresponds to the average allelic diversity of all the markers, see Methods) for the markers reaches 0.87 with values ranging from 0.50 to 0.96 (Additional file1: Table S1) and the percentage of polymorphic markers between each pair of Arabidopsis accessions ranges from 64% to 99% with an average value of 90% (Additional file1: Table S3). Thus, except for only a very few chromosomal regions, the SSR set distinguishes all the accessions (Additional file1: Table S2 and Additional file1: Table S4).
With this panel of SSR markers, any primary QTL mapping can be easily achieved with progenies (F2 or RIL populations) issued from any cross between the 24 Versailles accessions or between these accessions and Col-0, C24, Ws-2, Nd-1, Bay-0 and Ler-2. As soon as two accessions diverge for a trait of interest, the genetic factors associated with this phenotype can be located on the five Arabidopsis chromosomes. In addition, to simply revealing SSR polymorphisms, these SSRs are amplified using similar PCR conditions and subsequently displayed on agarose gels.
The high level of genetic diversity of these SSRs is in agreement with the high genetic diversity of the Versailles core collection[22, 31]. Within the collection these 96 markers display more than 1,300 different alleles. This high level of polymorphism also makes them markers of choice for other Arabidopsis segregating populations which do not belong to the Versailles collection as illustrated in our recent genetic mapping study. In addition, as only PCR and electrophoresis are needed to reveal polymorphism, this set of SSR markers has significant advantages and is complementary to the use of SNP markers[31, 32] or Simple Sequence Length Polymorphism (SSLP) markers which are based on indels.
Accession identification using a set of three highly polymorphic SSR markers
Due to the high level of genetic diversity displayed by the developed markers, we examined the possibility of differentiating each of the accessions with a minimum number of SSR markers. In this regard, we determined that using a select combination of only 3 SSR markers we could discriminate each accession from the 29 others. These markers are BSAT1.024, MSAT3.1 and MSAT 5.22 (Figure 2). Thus, these markers constitute a powerful tool to control seed stocks and crosses for the accessions studied in this work and likely for many other worldwide accessions as well.
Mining of microsatellites
The Arabidopsis BAC clones file containing the GenBank versions (fasta format) of sequences that make up the TAIR AGI map was subsequently searched for the presence of one to three perfect repeat motifs using the Sputnik software. The minimum length for each type of SSR is set as follows: mononucleotide repeats ≥30 nucleotides; dinucleotide repeats ≥30 nucleotides and trinucleotide repeats ≥24 nucleotides.
Primer design and PCR conditions
Sequences of the SSRs that were identified in the genomic areas of interest were then subjected to similarity searches in the TAIR Whole Genome (BAC clones) dataset through the BLASTN algorithm. SSRs with several matches along the genome were eliminated. Oligonucleotide primers were designed for selected SSR loci using the PRIMER3 software and were also subjected to BLASTN analysis. The parameters for primer design were as follows: amplicon size of 100–350 bp, primer size of 18–25 bp and primer melting temperature of 55–60°C with an optimum at 57°C.
Genomic DNA used for PCR was extracted from Arabidopsis young leaves using the NucleoSpin® Plant kit (Macherey-Nagel, Düren, Germany). PCRs were performed in a final volume of 15 μl consisting of 15 ng of DNA template, 1 U of Taq polymerase (Fermentas, Glen Burnie, MD, USA), 1x Taq reaction buffer, 0.2 mM dNTPs, 1.5 mM MgCl2, 200 nM of each primer. Conditions used for amplification were as follows: preincubation at 94°C for 2 min; followed by 35 cycles of denaturation at 94°C for 30 s; annealing at 57°C for the SSRs developed in this study, 50°C for the MSAT SSRs and 55°C for the other SSRs for 30 s; elongation at 72°C for 30s and a final extension step at 72°C for 10 min.
To reveal polymorphism, the PCR products were subsequently separated in 3% agarose gels. To determine the precise PCR product size, 6% acrylamide gels were used.
Polymorphic SSR markers were scored for the presence or absence of the corresponding bands among the tested accessions. Stutter and background bands were excluded. The scores ‘1’ and ‘0’ indicated presence and absence of bands, respectively. The polymorphism information content (PIC) of each marker was calculated according to as follows: PIC = 1-∑Pi2, where Pi is the band frequency of the ith allele. PIC provides an estimate of the discriminatory power of a locus by taking into account, not only the number of alleles but also the relative frequencies of those alleles. PIC values vary from 0 (monomorphic) to 1 (very highly discriminative, with many alleles in equal frequencies).
Nordborg M, Weigel D: Next-generation genetics in plants. Nature. 2008, 456: 720-723. 10.1038/nature07629.
Borevitz JO, Chory J: Genomics tools for QTL analysis and gene discovery. Curr Opin Plant Biol. 2004, 7: 132-136. 10.1016/j.pbi.2004.01.011.
Mauricio R: Mapping quantitative trait loci in plants: uses and caveats for evolutionary biology. Nat Rev Genet. 2001, 2: 370-381. 10.1038/35072085.
Clark RM, Schweikert G, Toomajian C, Ossowski S, Zeller G, Shinn P, Warthmann N, Hu TT, Fu G, Hinds DA: Common sequence polymorphisms shaping genetic diversity in arabidopsis thaliana. Science. 2007, 317: 338-342. 10.1126/science.1138632.
Bergelson J, Roux F: Towards identifying genes underlying ecologically relevant traits in Arabidopsis thaliana. Nat Rev Genet. 2010, 11: 867-879. 10.1038/nrg2896.
Brachi B, Faure N, Horton M, Flahauw E, Vazquez A, Nordborg M, Bergelson J, Cuguen J, Roux F: Linkage and association mapping of Arabidopsis thaliana flowering time in nature. PLoS Genet. 2010, 6: e1000940-10.1371/journal.pgen.1000940.
El-Lithy ME, Bentsink L, Hanhart CJ, Ruys GJ, Rovito D, Broekhof JLM, van der Poel HJA, van Eijk MJT, Vreugdenhil D, Koornneef M: New Arabidopsis recombinant inbred line populations genotyped using SNPWave and their use for mapping flowering-time quantitative trait loci. Genetics. 2006, 172: 1867-1876.
Filiault DL, Maloof JN: A genome-wide association study identifies variants underlying the Arabidopsis thaliana shade avoidance response. PLoS Genet. 2012, 8: e1002589-10.1371/journal.pgen.1002589.
Fournier-Level A, Korte A, Cooper MD, Nordborg M, Schmitt J, Wilczek AM: A map of local adaptation in Arabidopsis thaliana. Science. 2011, 334: 86-89. 10.1126/science.1209271.
Seymour DK, Filiault DL, Henry IM, Monson-Miller J, Ravi M, Pang A, Comai L, Chan SWL, Maloof JN: Rapid creation of Arabidopsis doubled haploid lines for quantitative trait locus mapping. Proc Natl Acad Sci. 2012, 109: 4227-4232. 10.1073/pnas.1117277109.
Simon M, Loudet O, Durand S, Bérard A, Brunel D, Sennesal F-X, Durand-Tardif M, Pelletier G, Camilleri C: Quantitative trait loci mapping in five new large recombinant inbred line populations of Arabidopsis thaliana genotyped with consensus single-nucleotide polymorphism markers. Genetics. 2008, 178: 2253-2264. 10.1534/genetics.107.083899.
Cao J, Schneeberger K, Ossowski S, Gunther T, Bender S, Fitz J, Koenig D, Lanz C, Stegle O, Lippert C: Whole-genome sequencing of multiple Arabidopsis thaliana populations. Nat Genet. 2011, 43: 956-963. 10.1038/ng.911.
Gan X, Stegle O, Behr J, Steffen JG, Drewe P, Hildebrand KL, Lyngsoe R, Schultheiss SJ, Osborne EJ, Sreedharan VT: Multiple reference genomes and transcriptomes for Arabidopsis thaliana. Nature. 2011, 477: 419-423. 10.1038/nature10414.
1001 Genomes. Adv Physiol Educ. [http://1001genomes.org/],
Jones N, Ougham H, Thomas H, Pašakinskienė I: Markers and mapping revisited: finding your gene. New Phytol. 2009, 183: 935-966. 10.1111/j.1469-8137.2009.02933.x.
Kalia R, Rai M, Kalia S, Singh R, Dhawan AK: Microsatellite markers: an overview of the recent progress in plants. Euphytica. 2011, 177: 309-334. 10.1007/s10681-010-0286-9.
Decroocq V, Favé MG, Hagen L, Bordenave L, Decroocq S: Development and transferability of apricot and grape EST microsatellite markers across taxa. Theor Appl Genet. 2003, 106: 912-922.
Casacuberta E, Puigdomènech P, Monfort A: Distribution of microsatellites in relation to coding sequences within the Arabidopsis thaliana genome. Plant Sci. 2000, 157: 97-104. 10.1016/S0168-9452(00)00271-5.
Symonds VV, Lloyd AM: An analysis of microsatellite loci in Arabidopsis thaliana: mutational dynamics and application. Genetics. 2003, 165: 1475-1488.
Zhang L, Yuan D, Yu S, Li Z, Cao Y, Miao Z, Qian H, Tang K: Preference of simple sequence repeats in coding and non-coding regions of Arabidopsis thaliana. Bioinformatics. 2004, 20: 1081-1086. 10.1093/bioinformatics/bth043.
Versailles Arabidopsis stock center.http://publiclines.versailles.inra.fr/,
McKhann HI, Camilleri C, Berard A, Bataillon T, David JL, Reboud X, Le Corre V, Caloustian C, Gut IG, Brunel D: Nested core collections maximizing genetic diversity in Arabidopsis thaliana. Plant J. 2004, 38: 193-202. 10.1111/j.1365-313X.2004.02034.x.
Cosson P, Schurdi-Levraud V, Le QH, Sicard O, Caballero M, Roux F, Le Gall O, Candresse T, Revers F: The RTM resistance to potyviruses in Arabidopsis thaliana: natural variation of the RTM genes and evidence for the implication of additional genes. PLoS One. 2012, 7: e39169-10.1371/journal.pone.0039169.
Pagny G, Paulstephenraj PS, Poque S, Sicard O, Cosson P, Eyquard J-P, Caballero M, Chague A, Gourdon G, Negrel L: Family-based linkage and association mapping reveals novel genes affecting Plum pox virus infection in Arabidopsis thaliana. New Phytol. 2012, 196: 873-886. 10.1111/j.1469-8137.2012.04289.x.
Decroocq V, Sicard O, Alamillo JM, Lansac M, Eyquard JP, Garcia JA, Candresse T, Le Gall O, Revers F: Multiple resistance traits control Plum pox virus infection in Arabidopsis thaliana. Mol Plant Microbe Interact. 2006, 19: 541-549. 10.1094/MPMI-19-0541.
Revers F, Guiraud T, Houvenaghel M-C, Mauduit T, Le Gall O, Candresse T: Multiple resistance phenotypes to lettuce mosaic virus among Arabidopsis thaliana accessions. Mol Plant Microbe Interact. 2003, 16: 608-616. 10.1094/MPMI.2003.16.7.608.
Sicard O, Loudet O, Keurentjes JJB, Candresse T, Le Gall O, Revers F, Decroocq V: Identification of quantitative trait loci controlling symptom development during viral infection in Arabidopsis thaliana. Mol Plant Microbe Interact. 2008, 21: 198-207. 10.1094/MPMI-21-2-0198.
The V.A.S.T lab.http://www7.inra.fr/vast/msat.php,
The Arabidopsis information resource.http://www.arabidopsis.org/index.jsp,
Guo W-J, Ling J, Li P: Consensus features of microsatellite distribution: microsatellite contents are universally correlated with recombination rates and are preferentially depressed by centromeres in multicellular eukaryotic genomes. Genomics. 2009, 93: 323-331. 10.1016/j.ygeno.2008.12.009.
Simon M, Simon A, Martins F, Botran L, Tisné S, Granier F, Loudet O, Camilleri C: DNA fingerprinting and new tools for fine-scale discrimination of Arabidopsis thaliana accessions. Plant J. 2012, 69: 1094-1101. 10.1111/j.1365-313X.2011.04852.x.
Anastasio AE, Platt A, Horton M, Grotewold E, Scholl R, Borevitz JO, Nordborg M, Bergelson J: Source verification of mis-identified Arabidopsis thaliana accessions. Plant J. 2011, 67: 554-566. 10.1111/j.1365-313X.2011.04606.x.
Păcurar DI, Păcurar ML, Street N, Bussell JD, Pop TI, Gutierrez L, Bellini C: A collection of INDEL markers for map-based cloning in seven Arabidopsis accessions. J Exp Bot. 2012, 63: 2491-2501. 10.1093/jxb/err422.
The European Arabidopsis stock centre.http://arabidopsis.info/,
The Arabidopsis BAC clones file.ftp://ftp.arabidopsis.org/home/tair/Sequences/clones/AGI_BAC,
Sputnik: DNA microsatellite repeat search utility.http://services.cbib.u-bordeaux2.fr/pise/sputnik.html,
Smith JSC, Chin ECL, Shu H, Smith OS, Wall SJ, Senior ML, Mitchell SE, Kresovich S, Ziegle J: An evaluation of the utility of SSR loci as molecular markers in maize (Zea mays L.): comparisons with data from RFLPS and pedigree. Theor Appl Genet. 1997, 95: 163-173. 10.1007/s001220050544.
The authors thank Thierry Mauduit, Marylin Roncoroni and Aurélie Bailly for the production and maintenance of the Arabidopsis plants and Albert G. Abbott for his critical reading of the manuscript.