Highly Informative Single-Copy Nuclear Microsatellite DNA Markers Developed Using an AFLP-SSR Approach in Black Spruce (Picea mariana) and Red Spruce (P. rubens)

Background

Microsatellites or simple sequence repeats (SSRs) are highly informative molecular markers for various biological studies in plants. In spruce (Picea) and other conifers, the development of single-copy polymorphic genomic microsatellite markers is quite difficult, owing primarily to the large genome size and predominance of repetitive DNA sequences throughout the genome. We have developed highly informative single-locus genomic microsatellite markers in black spruce (Picea mariana) and red spruce (Picea rubens) using a simple but efficient method based on a combination of AFLP and microsatellite technologies.

Principal Findings

A microsatellite-enriched library was constructed from genomic AFLP DNA fragments of black spruce. Sequencing of the 108 putative SSR-containing clones provided 94 unique sequences with microsatellites. Twenty-two of the designed 34 primer pairs yielded scorable amplicons, with single-locus patterns. Fourteen of these microsatellite markers were characterized in 30 black spruce and 30 red spruce individuals drawn from many populations. The number of alleles at a polymorphic locus ranged from 2 to 18, with a mean of 9.3 in black spruce, and from 3 to 15, with a mean of 6.2 alleles in red spruce. The polymorphic information content or expected heterozygosity ranged from 0.340 to 0.909 (mean = 0.67) in black spruce and from 0.161 to 0.851 (mean = 0.62) in red spruce. Ten SSR markers showing inter-parental polymorphism inherited in a single-locus Mendelian mode, with two cases of distorted segregation. Primer pairs for almost all polymorphic SSR loci resolved microsatellites of comparable size in Picea glauca, P. engelmannii, P. sitchensis, and P. abies.

Significance

The AFLP-based microsatellite-enriched library appears to be a rapid, cost-effective approach for isolating and developing single-locus informative genomic microsatellite markers in black spruce. The markers developed should be useful in black spruce, red spruce and other Picea species for various genetics, genomics, breeding, forensics, conservation studies and applications.     

Development of a set of highly polymorphic genomic microsatellites (gSSRs) in Sitka spruce (Picea sitchensis (Bong.) Carr.)

Robust, polymorphic microsatellite DNA markers (simple sequence repeats—SSRs) are valuable tools for a range of tree conservation and breeding applications. SSRs are routinely used in the study of population genetic structure and diversity, pedigree reconstruction and genetic linkage mapping. Their abundance in the genome, co-dominant inheritance and potential for cross-species amplification make microsatellites highly prized markers. This paper characterises 22 novel genomic polymorphic microsatellite loci for Sitka spruce (Picea sitchensis (Bong.) Carr.). Amplification of DNA from Sitka spruce material was carried out both with a set of unrelated trees to obtain diversity statistics for each locus, and with the progeny of a full-sib family to test simple Mendelian inheritance. Observed heterozygosity ranged from 0.38 to 0.91 and allele number per locus ranged from 6 to 21, with a mean of 12.2. In addition, the primer pairs were tested with DNA from Norway spruce (P. abies) and white spruce (P. glauca) to investigate their potential for cross-species amplification and ten loci amplified in all three species. The results from these genomic microsatellites are compared to data generated from microsatellites derived from Picea EST libraries. In summary, this novel, highly polymorphic markers represent a significant addition to the rapidly expanding Picea genomics tool-box. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A set of microsatellite markers for use in Sitka spruce (Picea sitchensis) developed from Picea glauca ESTs

Few microsatellite markers have been specifically developed for Picea sitchensis. In January 2004 the appearance of over 10 000 sequences of expressed regions of DNA from P. glauca in GenBank presented an opportunity for the development of additional microsatellite markers in Sitka spruce. Mono‐ and dinucleotide repeat sequences were located in these sequences and primers were designed around these regions. Amplification was attempted in Sitka material from a broad geographical range and the level of polymorphism was assessed. Primers were also tested in progeny of a controlled cross. Nine polymorphic loci that demonstrated Mendelian inheritance in Sitka were discovered in this study.     

Pcr-based multiplex DNA fingerprinting techniques for the analysis of conifer genomes

We utilised fingerprinting techniques such as amplified restriction fragment polymorphism (AFLP) and selective amplification of microsatellite polymorphic loci (SAMPL) in order to identify molecular markers in the conifer species Norway spruce (Picea abies K). A total of 366 polymorphic AFLP and 20 polymorphic SAMPL bands were scored, with an average of 12.6 polymorphic bands obtained per primer combination. Of the bands, 94% segregated in a Mendelian fashion, and 14% of the AFLP bands behaved as codominant markers. Norway spruce, as well as other conifers, has a very large genome (1.5–2.5 × 10¹⁰ bp) which is thought to be largely composed of repetitive sequences of different origin, either organised in tandem or interspersed in the genome. The high complexity of these genomes has so far hindered the efficient use of multiplexed fingerprinting techniques for their analysis. The use of the methylation-sensitive restriction enzyme PstI as 6 bp cutter allows a reduction in pattern complexity as compared to EcoRI. We thus demonstrated that high-throughput PCR-based multiplex fingerprinting techniques such as AFLP and SAMPL, when opportunely modified, can be effectively applied also to these species and be useful for mapping purposes.     

Development of microsatellite markers for white spruce (Picea glauca) and related species

We report the development of 13 primer pairs that allow the unambiguous amplification of 15 microsatellite (SSR) loci in white spruce (Picea glauca). Fourteen of these loci were polymorphic in trees sampled at three geographically separated regions of western Canada. Segregation analysis carried out on these loci confirmed a Mendelian inheritance pattern for all except two, which showed significant segregation distortion. All of these primer pairs amplified SSR loci in at least one of the other Picea species tested [black spruce (P. mariana), red spruce (P. rubens), Norway spruce (P. abies), Colorado spruce (P. pungens), sitka spruce (P. sitchensis) and Engelmann spruce (P. engelmannii)]. Given the important commercial and ecological roles of these species, this set of markers will be invaluable for their management, the improvement of commercially important traits, and the study of their ecology and genetics. Received: 18 August 2000 / Accepted: 28 September 2000     

A high-density genetic linkage map of a black spruce (Picea mariana) × red spruce (Picea rubens) interspecific hybrid

Genetic maps provide an important genomic resource of basic and applied significance. Spruce (Picea) has a very large genome size (between 0.85 × 1010 and 2.4 × 1010 bp; 8.5-24.0 pg/1C, a mean of 17.7 pg/1C ). We have constructed a near-saturated genetic linkage map for an interspecific backcross (BC1) hybrid of black spruce (BS; Picea mariana (Mill.) B.S.P.) and red spruce (RS; Picea rubens Sarg.), using selectively amplified microsatellite polymorphic loci (SAMPL) markers. A total of 2284 SAMPL markers were resolved using 31 SAMPL-MseI selective nucleotide primer combinations. Of these, 1216 SAMPL markers showing Mendelian segregation were mapped, whereas 1068 (46.8%) SAMPL fragments showed segregation distortion at α = 0.05. Maternal, paternal, and consensus maps consistently coalesced into 12 linkage groups, corresponding to the haploid chromosome number (1n = 1x = 12) of 12 in the genus Picea. The maternal BS map consisted of 814 markers distributed over 12 linkage groups, covering 1670 cM, with a mean map distance of 2.1 cM between adjacent markers. The paternal BS × RS map consisted of 773 markers distributed over 12 linkage groups, covering 1563 cM, with a mean map distance of 2.0 cM between adjacent markers. The consensus interspecific hybrid BC1 map consisted of 1216 markers distributed over 12 linkage groups, covering 1865 cM (98% genome coverage), with a mean map distance of 1.5 cM between adjacent markers. The genetic map reported here provides an important genomic resource in Picea, Pinaceae, and conifers.     

DNA fingerprinting based on microsatellite-anchored fragment length polymorphisms,and isolation of sequence-specific PCR markers in lupin (Lupinus angustifolius L.)

We report a method of microsatellite-anchored fragment length polymorphisms for DNA fingerprinting. The method combines the concept of AFLP and the microsatellite-anchor primer technique. Genomic DNA was digested by one restriction enzyme MseI. One AFLP adaptor (MseI adaptor) was ligated onto the restriction fragments. DNA fingerprints were produced by PCR using one microsatellite-anchor primer in combination with one MseI-primer. The method allows co-amplification of over 100 DNA fragments containing microsatellite motifs per PCR. Polymorphisms detected from lupin by this method included those arising from variation in the number of microsatellite repeat units targeted by the microsatellite-anchor primers, from variation on the annealing sites for the SSR-anchor primers, from insertions/deletions outside the SSR region, and from variation in restriction sites. The first three types of polymorphisms were readily converted into sequence-specific PCR markers suitable for marker-assisted breeding.     

Testing of microsatellite primers with different populations of Eurasian spruces Picea abies (L.) Karst. and Picea obovata Ledeb.

From a clone library containing microsatellite DNA fragments of Norway spruce, seven pairs of primers were selected. These primers were tested as markers in the genetic structure analysis of nine populations of Eurasian spruce species Picea abies (L.) Karst. and Picea obovata Ledeb. Five pairs of these primers identified polymorphic loci with allele numbers ranging from 6 to 15. In the populations examined, the observed and expected heterozygosity values assessed at five loci varied from 0.1778 to 0.6556 and from 0.7800 to 0.900, respectively. In the populations examined, the values of F _st index varied from 0.0691 to 0.2551 with the mean value of 0.1318. On the dendrogram based on Nei genetic distances, the populations formed three groups: Pskov-Ciscarpathya, Komi-Tatarstan-Arkhangelsk, Kazakhstan-Karelia(natural)-Karelia(plantation)-Krasnoyarsk. Five of the primer pairs tested proved to be useful for analysis of the population genetic structure in Eurasian spruce species.     

Rapid identification of white-Engelmann spruce species by RAPD markers

Fragments of random amplified polymorphic DNA (RAPDs) were used as markers to distinguish Picea glauca (Moench) Voss (white spruce) and Picea engelmannii Parry (Engelmann spruce). These species and their putative hybrids are difficult to differentiate morphologically and are collectively known as interior spruce. Four oligodeoxynucleotide decamer primers showed species-specific amplification products between white spruce and Engelmann spruce. These fragments are highly conserved among seed lots and individual trees of each species from diverse geographic origins. The consistency and reproducibility of these species-specific amplification products were tested in more than two amplification reactions. Therefore, RAPD markers can provide genetic markers for easy and rapid identification of the specific genetic entry of these spruce species and their reported putative hybrids. According to the frequencies of the species-specific RAPD markers, it is possible to estimate the hybrid fraction, indicative of true introgression between the two species. These results are useful for quick identification of both species and their hybrid swarms at any stage in the sporophyte phase of the life cycle, for determining the occurrence and the magnitude of introgressive hybridization in an overlap zone between the two species, and for certification purposes in operational re-forestation and tree-improvement programs.     

A genetic map of an interspecific cross in Allium based on amplified fragment length polymorphism (AFLPTM) markers

Segregation of 692 polymorphic AFLP^TM (amplified fragment length polymorphism) fragments was determined in an F₂ of the interspecific cross A. roylei x A. cepa. Two different enzyme combinations were used, PstI/MseIand EcoRI/MseI; in the latter one extra selective nucleotide was added to the MseI primer. The map based on A. cepa markers consisted of eight linkage groups with 262 markers covering 694 cM of the expected 800 cM. The map based on A. roylei markers comprised 15 linkage groups with 243 markers and had a length of 626 cM. The two maps were not integrated, and 25% of the markers remained unlinked. One of the alliinase genes and a SCAR marker linked to the disease resistance gene to downy mildew are present on this map. Of the AFLP markers, 50—80% were polymorphic between A. cepa and A. roylei; the level of polymorphic markers between different A. cepa accessions was4-8%. Received: 28 August 1998 / Accepted: 31 March 1999     

Identification, genetic localization, and allelic diversity of selectively amplified microsatellite polymorphic loci in lettuce and wild relatives (Lactuca spp.).

Selectively amplified microsatellite polymorphic locus (SAMPL) analysis is a method of amplifying microsatellite loci using generic PCR primers. SAMPL analysis uses one AFLP primer in combination with a primer complementary to microsatellite sequences. SAMPL primers based on compound microsatellite sequences provided the clearest amplification patterns. We explored the potential of SAMPL analysis in lettuce to detect PCR-based codominant microsatellite markers. Fifty-eight SAMPLs were identified and placed on the genetic map. Seventeen were codominant. SAMPLs were dispersed with RFLP markers on 11 of the 12 main linkage groups in lettuce, indicating that they have a similar genomic distribution. Some but not all fragments amplified by SAMPL analysis were confirmed to contain microsatellite sequences by Southern hybridization. Forty-five cultivars of lettuce and five wild species of Lactuca were analyzed to determine the allelic diversity for codominant SAMPLs. From 3 to 11 putative alleles were found for each SAMPL; 2-6 alleles were found within Lactuca sativa and 1-3 alleles were found among the crisphead genotypes, the most genetically homogeneous plant type of L. sativa. This allelic diversity is greater than that found for RFLP markers. Numerous new alleles were observed in the wild species; however, there were frequent null alleles. Therefore, SAMPL analysis is more applicable to intraspecific than to interspecific comparisons. A phenetic analysis based on SAMPLs resulted in a dendrogram similar to those based on RFLP and AFLP markers.     

Genetic mapping of a fusarium wilt resistance gene (Fom-2) in melon (Cucumis melo L.)

Fusarium wilt caused by Fusarium oxysporum f.sp. melonis is one of the most devastating diseases in melon production worldwide. The most effective control measure available is the use of resistant varieties. Identifying molecular markers linked to resistance genes can serve as a valuable tool for the selection of resistant genotypes. Bulked segregant analysis was used to identify markers linked to the Fom-2 genes, which confers resistance to races 0 and 1 of the fungal pathogen. Pooled DNA from homozygous resistant or homozygous susceptible progeny of F₂ cross between MR-1 and AY was screened using 240 PstI/MseI and 200 EcoRI/MseI primer combinations to identify AFLP markers linked to Fom-2. Fifteen markers potentially linked to Fom-2 were identified, all with EcoRI/MseI primer pairs. These were mapped relative to Fom-2 in a backcross (BC) population of 60 progeny derived from MR-1 × AY with AY as recurrent parent. Two AFLP markers (ACT/CAT1 and AAC/CAT1) flanked the gene at 1.7 and 3.3 cM, respectively. Moreover, AFLP marker AGG/CCC and the previously identified RAPD marker 596-1 cosegregated with Fom-2. These two dominant markers were converted to co-dominant markers by designing specific PCR primers that produced product length polymorphisms between the parents. A survey of 45 melon genotypes from diverse geographic origins with the co-dominant markers demonstrated a high correlation between fragment size and the resistance phenotype. These markers may therefore be useful in marker-assisted breeding programs.     

PCR-based molecular markers for assessment of somaclonal variation in <Emphasis Type="Italic">Pinus pinea</Emphasis> clones micropropagated <Emphasis Type="Italic">in vitro</Emphasis>

Four different markers [random amplified polymorphic DNA (RAPD), inter simple sequence repeat (ISSR), amplified fragment length polymorphism (AFLP), and selective amplified microsatellite polymorphism length (SAMPL)] were applied for evaluating somaclonal variation of micropropagated genotypes of stone pine (Pinus pinea L.). The total number of primers tested was 130, with 223 combinations assayed. A high number of them amplified successfully (178), representing 79.82 % of the total, and the average number of amplified fragments ranged from 2.47 (ISSR) to 65.76 (SAMPL). Based on internal controls, no problem of reproducibility was detected. Almost no somaclonal variation was detected within the clones. Of the tested markers, ISSR, AFLP, and SAMPL showed monomorphic amplification profiles, with only RAPD markers showing some interclonal variation.     

Testing of microsatellite primers with different populations of Eurasian spruces Picea abies (L.) Karst. and Picea obovata Ledeb

From a clone library containing microsatellite DNA fragments of Norwegian spruce, seven pairs of primers were selected. These primers were tested to be the markers in the genetic structure analysis of nine populations of Eurasian spruce species Picea abies (L.) Karst. and Picea obovata Ledeb. Five pairs of these primers identified polymorphic loci with the allele numbers from 6 to 15. In the populations examined, the observed and expected heterozygosity values assessed at five loci varied from 0.1778 to 0.6556 and from 0.7800 to 0.900, respectively. In the populations examined, the values of F(st) index varied from 0.0691 to 0.2551 with the mean value of 0.1318. On the dendrogram based on Nei genetic distances, the populations formed three groups: Pskov-Ciscarpathia, Komi-Tatarstan-Arkhangelsk, Kazakhstan-Karelia(natural)-Karelia(culture)-Krasnoyarsk. Five of the primer pairs tested proved useful for analysis of the population genetic structure in Eurasian spruce species.     

Dinucleotide microsatellite primers designed for a critically endangered primate,the black lion tamarin (Leontopithecus chrysopygus)

Black lion tamarin (BLT) monkeys (Leontopithecus chrysopygus) have suffered a severe reduction in their natural range and are consequently critically endangered. Because allozyme data showed very low levels of variation, it was not clear if these monkeys had much genetic diversity. We designed microsatellite primers for BLTs, and from them we identified nine polymorphic loci, seven of which were tested on golden lion tamarins (GLTs) (Leontopithecus rosalia). All of the seven polymorphic loci and two other monomorphic BLT loci were polymorphic in GLTs. The microsatellite markers identified here are directly applicable to ongoing lion tamarin population and conservation genetics studies.     

Segregation of random amplified DNA markers in F1 progeny of conifers

Summary The recently developed approach to deriving genetic markers via amplification of random DNA segments with single primers of arbitrary nucleotide sequence was tested for its utility in genetic linkage mapping studies with conifers. Reaction conditions were optimized to reproducibly yield clean and specific amplification products. Template DNA from several genotypes of Douglas-fir (Pseudotsuga menziesii) and white spruce (Picea glauca) were tested against eight ten-base oligonucleotide primers. Most of the tested primer/parent tree combinations yielded polymorphic PCR products (RAPD markers). Selected primers were then used in PCR reactions with template DNA isolated from offspring in Douglas-fir and black spruce diallel crosses among the same parental lines. The diallel study confirmed the appropriate inheritance of RAPD markers in the F₁ generation. The value of these dominant RAPD markers for genetic linkage mapping in trees was established from both theoretical and applied perspectives.     

Analysis of genetic diversity through AFLP,SAMPL, ISSR and RAPD markers in <Emphasis Type="Italic">Tribulus terrestris</Emphasis>, a medicinal herb

Tribulus terrestris is well known for its medicinal importance in curing urino-genital disorders. Amplified fragment length polymorphism (AFLP), selective amplification of microsatellite polymorphic loci (SAMPL), inter-simple sequence repeat (ISSR) and randomly amplified polymorphic DNA (RAPD) markers were used for the first time for the detection of genetic polymorphism in this medicinal herb from samples collected from various geographical regions of India. Six assays each of AFLP and SAMPL markers and 21 each of ISSR and RAPD markers were utilized. AFLP yielded 500 scorable amplified products, of which 82.9% were polymorphic. SAMPL primers amplified 488 bands, 462 being polymorphic (94.7%). The range of amplified bands was 66 [(TC)₈G + M-CAG] to 98 [(CA)₆AG + M-CAC] and the percentage polymorphism, 89.9 [from (CT)₄C (AC)₄A + M-CTG] to 100 [from (GACA)₄ + M-CTA]. The ISSR primers amplified 239 bands of 0.4–2.5 kb, 73.6% showed polymorphism. The amplified products ranged from 5 to 16 and the percentage polymorphism 40–100. RAPD assays produced 276 bands, of which 163 were polymorphic (59%). Mantel test employed for detection of goodness of fit established cophenetic correlation values above 0.9 for all the four marker systems. The dendrograms and PCA plots derived from the binary data matrices of the four marker systems are highly concordant. High bootstrap values were obtained at major nodes of phenograms through WINBOOT software. The relative efficiency of the four molecular marker systems calculated on the basis of multiplex ratio, marker index and average heterozygosity revealed SAMPL to be the best. Distinct DNA fingerprinting profile, unique to every geographical region could be obtained with all the four molecular marker systems. Clustering can be a good indicator for clear separation of genotypes from different regions in well-defined groups that are supported by high bootstrap values.     

Development of polymorphic microsatellite markers in barfin flounder (<Emphasis Type="Italic">Verasper moseri</Emphasis>) and spotted halibut (<Emphasis Type="Italic">Verasper variegatus</Emphasis>) by the cross-species amplification

Barfin flounder (Verasper moseri) and spotted halibut (Verasper variegatus) are two commercially important flatfish species in the Northeast Asia. In the present study, we reported polymorphic microsatellite markers in V. moseri and V. variegatus by the cross-species amplification of microsatellite primers developed previously in two other related marine fish species. A total of 244 polymorphic microsatellite markers were selected for cross-species amplification in V. moseri and V. variegatus, of which 182 markers deriving from Atlantic halibut (Hippoglossus hippoglossus) and 62 markers deriving from Japanese flounder (Paralichthys olivaceus). A sample of 10 individuals were detected. As a result, a total of 67 loci showed polymorphisms in V. moseri, and 62 loci showed polymorphisms in V. variegatus, with the observed number of alleles per locus ranging from two to five in V. moseri, and from two to seven in V. variegatus, respectively. This paper provided more candidate microsatellite markers which could be useful for construction of genetic linkage maps, evaluation of population genetic structure and stock management of V. moseri and V. variegatus.     

RAPD variations among pure and hybrid populations ofPicea mariana, P. rubens, andP. glauca (Pinaceae), and cytogenetic stability ofPicea hybrids: Identification of species-specific RAPD markers

Random amplified polymorphic DNA (RAPD) analysis was used to determine genetic relationships amongP. mariana (black spruce),P. rubens (red spruce), andP. glauca (white spruce) and to assess the degree of polymorphism within populations from different provenances and among spruce hybrids. Eleven arbitrary decamer primers were used to amplify genomic DNAs extracted from embryogenic cultures and seedlings. Species-specific RAPD markers were identified.Picea mariana andP. rubens showed similar RAPD profiles confirming their close genetic relationship. Species-specific RAPD markers were identified and were useful in distinguishing white spruce from black and red spruces. RAPD differentiation between populations within each species was small. The level of polymorphism was much higher in spruce hybrid populations than in the pure species. Cytological analysis ofP. mariana ×P. rubens hybrids showed normal mitotic behaviour at prophase, metaphase, anaphase, and telophase. All the hybrids analyzed from different cross combinations were euploids.     

Characterization of microsatellite loci in the spruce bark beetle Ips typographus (Coleoptera: Scolytinae)

Ips typographus is an economically important pest of Norway spruce stands. We developed five polymorphic microsatellite markers using a biotin enrichment protocol. The number of alleles ranged from three to 11 per locus and no strong evidence for null alleles was found. Heterozygosity ranged from 0.3 to 0.97. These markers could be useful tools to study the population structure and genetic consequences of I. typographus population outbreaks.