Characterization of variable EST SSR markers for Norway spruce (Picea abies L.)

Natural populations of Picea abies L. (Norway spruce) are found from north-western Europe outside permafrost areas down to northern Greece, westwards to the Massif Central (France) and east to the Ural Mountains. Picea abies is growing above 400-500 m and ascends close to 2000 m in the Alps. Studies on genetic variation based on allozymes have shown that Picea abies genetic differentiation among populations is rather low over its whole distribution range [1, 2]. Previous studies on the genetic structure of P. abies using organelle markers showed pronounced differentiation between north-east boreal origins and areas in the central European mountains [3, 4], supporting the hypothesis of two distinct main glacial refugia as postulated from pollen data [5].

Initial reports on the occurrence of SSRs in conifers such as Pinus radiata[6], Pinus sylvestris[7] or Picea abies[8] have shown that marker development for such complex genomes is difficult. Frequently several DNA fragments in addition to the expected ones are amplified when using primers flanking putative SSR regions. This can be attributed to the very large size of the conifer genomes, with Norway spruce having 39.140 Mb in the diploid genome (2n = 24) as well as the high proportion of repetitive DNA in the conifer genomes [9]. These repetitive elements as well as pseudogenes frequently produce complex amplification products from multiple loci [10, 11]. To overcome this problem, Scotti et al. [12] showed that by isolating di-nucleotide SSRs from cDNA libraries of conifer species this problem can be overcome relying on the fact that expressed genes are less likely to be of repetitive nature in the genome.

Following this line of exploiting the advantages of publicly available data on expressed sequences (ESTs) as source for marker development [13] as described also by Rungis et al. [14] for conifers, we used Norway spruce EST sequences from NCBI for in silico identification of SSR regions which could serve for marker development.

Several studies describe the development of SSR markers for Norway spruce [8, 12]. In the present work we used EST mining of sequences available in databases instead of classical approaches like screening of genomic libraries or fragment enrichment and subsequent cloning. The newly developed markers can be used to complete unsaturated maps, and might be useful in marker assisted breeding and population genetics.

The high number of loci deviating from Hardy-Weinberg equilibrium could partly be explained by the presence of null alleles (Table 3), partly be a result of sampling or selective pressure on the coding regions which is supported by the fact that loci deviating from HWE are not consistent between the two populations. Null alleles can occur due to mutations in primer binding sites and lead to the overestimation of homozygosity as shown by Callen [16]. The presence of such null alleles as estimated in this study may further be confirmed by synthesis of alternative oligonucleotide primer pairs. Increased appearance of homozygotes in some loci supports the hypothesis of selective pressure which might be caused by advantages of recessive or dominant homozygotes known as directional selection. Positive selection was confirmed by outlier detection for two loci which show a high frequency of homozygotes. Both of them are located in coding regions - Pa_42 shows homology to a cysteine protease and Pa_51 is related to a cadmium transporting ATPase. These genes play important roles in various physiological processes, including response to biotic and abiotic stress [17, 18]. Therefore, additional data of both populations (e.g. environmental growth) conditions would be of interest. Especially locus Pa_51 may be an interesting target for diversity studies because departure from HWE was only detected in the Mayrhofen population. But although these two loci show considerable degrees of population differentiation, the low F_st values present in the other loci support previous findings [3, 4] that there is no major degree of differentiation among Alpine Norway spruce populations.

The new SSR markers characterized in this study are a good resource for assessing diversity in Alpine spruce populations and might be of further use for applications in forestry. Additionally, 31 monomorphic SSR loci were identified in this study. It is still possible to derive polymorphic data from them when screening Norway spruce populations from different origin, as it has already been shown that flanking regions of monomorphic SSRs show a moderate to high degree of variability [19]. Possible variation within these flanking regions in Norway spruce could serve as alternative source for diversity studies applying different techniques for the identification of variability within these regions.

All variable SSRs are located within ORF sequences and 60% of them gave significant hits to functional genes present in GenBank. A frameshift caused by length variation within an ORF might lead to distortion in the coding region and thus to the abortion of the coded protein. In all of the newly identified EST-SSRs no frameshift is occurring because either, as in 50% of the cases, the SSRs are trinucleotide repeats or the repeat number of polymorphic di- and pentanucleotide repeats alters in such a way that the sum of the basepair variation is a multiple of three.

We have identified 92 new SSR regions with 48 different repeat motifs in the expressed part of the Norway spruce genome. In two screening populations of 48 individuals each, 23 SSR regions showed a moderate to high degree of polymorphism. Although a high number of SSR markers were already developed for Norway spruce, most of the available SSRs consist of dinucleotide repeats and do not include our tri- and pentanucleotide repeats. We have demonstrated that on the one hand genetic differentiation between distinct stands is low but on the other hand selective forces are likely to have been acting on several genes because deviations from HWE could not only be explained due to null alleles and test for outliers revealed two loci under positive selection. The results gave an insight into the abundant repeat classes and will be of use in analysis of variation linked to expressed genes in Alpine Norway spruce populations as well as evolutionary forces acting at these loci.

14,022 Norway Spruce EST sequences were extracted from the NCBI dbEST and screened for SSR sequences with the SciRoKo software program [20]. Only perfect SSRs were searched in two iterations, one for mono- and dinucleotides (at least 4 repeats), and another one for tri-, tetra-, penta- and hexanucleotides (at least 3 repeats). SSR motif containing sequences were subjected to clustering and annotation using the EST2uni pipeline [21]. For clustering, the default settings (30 bp minimum overlap with at least 94% identity for pairwise alignment using BLAST, 93% overlap identity cutoff for assembly using CAP3) were applied. The resulting Unigene Set (contigs and singletons) was compared against three protein databases with BLASTX, in particular NCBI NR, Arabidopsis TAIR7, and Uniprot/Swissprot, using an e-value of 1e-10 as cutoff. The descriptions of the BLAST hits obtained with the different BLAST runs were parsed and merged to yield a descriptive annotation for each unigene. The annotations were attributed with modifiers like "Similar to" or "Highly similar to", depending on the e-value of the alignment with the corresponding BLAST hit.

Plant material used: 16 samples from all over Austria (testing population) and two times 48 individuals from two distinct populations (Gusswerk, Styria [N 47° 44|, E 15° 18|]; Mayrhofen, Tyrol [N 47° 11', E 11° 52']; screening population). DNA was kindly provided by the Federal Research and Training Centre for Forests, Natural Hazards and Landscape (BFW).

60 primer pairs were designed using Primer3 web interface [22] with default settings. PCR reactions were performed in 25 μl volumes containing 50 ng genomic DNA, 1X PCR Buffer (QIAGEN), 1 mM MgCl₂, 0.2 mM each dNTP, 0.3 μM FAM-labelled forward primer, 0.3 μM reverse primer and 0.625 U HotStar Taq polymerase (QIAGEN). The PCR thermal profile consisted of 15 min. initial denaturation at 95 °C followed by 35 cycles denaturation at 95 °C for 50 sec., a primer-specific annealing-temperature (Table 1) for 50 sec., extension at 72 °C for 105 sec., and a final extension at 72 °C for 10 minutes. Fragments amplified in 16 individuals from all over Austria (testing population) and 96 individuals from 2 distinct locations in Austria (screening population) were analyzed on an ABI 3100 Genetic Analyzer. Fragment sizes were extracted using PeakScanner 1.0 (Applied Biosystems, Foster City, California, USA) and allele calls were manually performed in Excel. The testing population was only used to test for polymorphic loci and was not included in further calculations. Number of alleles, effective number of alleles, observed and expected heterozygosities and deviation from Hardy-Weinberg equilibrium were calculated using Genepop 4.0 [23]. Fixation indices and fixation indices excluding null alleles (ENA), as well as estimates of null alleles were calculated with FreeNA [24] and SSR presence in coding regions was analyzed with the GetORF software using default settings. Test for outlier loci were performed with LOSITAN [25] using an infinite allele model with 100.000 simulations, a confidence interval of 0.95 and a subsample size of 75.