BACKGROUND
single nucleotide polymorphims  have become a genomic commodity as they are becoming indispensable in various genome scans aimed at mapping genomes  <cit> , finding associations with complex traits  <cit> , and population genomics  <cit> .
they are distributed along the various regions of the genomes  <cit>  and are frequent in coding regions of angiosperms  <cit>  and conifers  <cit> .
however, the efficiency of genome scans is not only dependant on a wide genomic distribution of snps.
indeed, it also relies on the ability to genotype large numbers of snps over large sets of individuals.
snp genome scans in non model species usually involve two steps: the discovery of snps and genotyping.
with no a priori knowledge of dna polymorphisms, snps are usually discovered through various strategies of individual or pool dna sequencing  <cit> , or by using tilling techniques, a high-throughput strategy relying on the enzymatic cleavage of mismatches  <cit> .
for a number of crop species, current resequencing efforts have led to the development of snp databases and generate a wealth of snps usable in genome scans.
in conifers, large-scale est sequencing projects have been initiated  <cit> , providing a starting point to develop snp resources in pine  <cit>  and spruce  <cit> .
several snp genotyping array approaches have been developed with variable success.
the accuracy of innovative snp genotyping technologies has been assessed mostly through the development of assays suitable for analysing variations in the human genome.
broadly speaking, four reaction principles govern snp genotyping assays: hybridization with allele-specific oligonucleotide probes, oligonucleotide ligation, single nucleotide primer extension, and enzymatic cleavage reviewed in  <cit> .
among these approaches, the goldengate assay developed by illumina and relying on the bead array technology has demonstrated high performance with high levels of call rate, reproducibility, and overall success rate for the analysis of the human genome  <cit> .
high-throughput snp assays have recently been applied to plants.
large datasets of snp-based markers are being developed in barley through the development of genotyping assays relying on illumina's technologies  <cit> , leading to the undertaking of an international snp project  <cit> .
the same genotyping approach has made it possible to map large datasets of snps even in complex and duplicated genomes such as soybean  <cit> , and projects are underway in hexaploid wheat  <cit>  and poplar  <cit> .
in the present study, we are asking whether high-throughput snp genotyping technologies developed for human population genomics applications, such as the illumina goldengate snp assay, are applicable to large and essentially unsequenced genomes as seen in conifers.
conifer genomes reach very large sizes, around  <dig> – <dig>  mb  <cit> , consisting mostly of repetitive sequences  <cit> .
for the two conifers considered herein, white spruce and black spruce, genome sizes are well in excess of 10e <dig> bp  <cit> .
moreover, the partial knowledge of the large and redundant genomes of conifers can be a limiting factor to design an efficient snp genotyping assay.
indeed, sequences located upstream and downstream the snp cannot be fully validated for locus specificity and the possible presence of repetitive elements  <cit> .
the possible effect of such a drawback remains to be verified for most crop and tree species which genomes are essentially not sequenced.
based on est sequence data available for white spruce  <cit> , we have designed primers and resequenced genomic dna for hundreds of genes in white spruce and black spruce.
the high quality snp datasets developed were used to select snps amenable to the goldengate genotyping assay and test the technology for these two species.
then, we integrated these snp data into linkage maps of expressed genes and illustrated the possibility to rapidly improve the density of existing genetic maps for spruce species.
methods
plant material for genotyping assays and linkage mapping
plant material consisted of one outbred f <dig> cross # c <dig>  for white spruce  with  <dig> progeny derived from two parents selected for their high level of heterozygosity for estp anchor markers and for their intermediate performance for a number of traits such as somatic embryogenic capacity  <cit> .
plant material for black spruce  was representative of the hybridizing species complex picea mariana × picea rubens  and consisted of  <dig> progeny derived from the backcross-like cross bc <dig> , as previously described  <cit> .
spruces all harbor 2n = 2x =  <dig> chromosomes, as for most pinaceae  <cit> .
the lineages leading to these two taxa are thought to have diverged at least  <dig> myr ago  <cit> .
gene selection and pcr amplification
for snp discovery, we considered a collection of  <dig>  unigenes derived from a white spruce est database  <cit> .
for each gene, coding regions were identified based on alignments with similar sequences from uniprot-swissprot protein databases.
pcr primers for amplification and resequencing were designed using primer <dig>  <cit> .
whenever possible, one of the primers was anchored outside of the coding regions in order to increase amplification specificity.
the same primers were also used for snp discovery in black spruce.
the lists of unigene and primer sequences for the genes represented on the two illumina snp bead arrays for white spruce and black spruce are provided .
for each of white spruce and black spruce, dna was isolated from the two parents of the mapping population and their progeny, as well as from haploid megagametophyte tissue using a dneasy plant mini kit .
about 5– <dig> ng of template dna were used for pcr amplification.
reactions were done in  <dig> μl containing  <dig> mm tris-hcl ,  <dig> mm kcl,  <dig> – <dig>  mm mgcl <dig>   <dig> μm of each dntp,  <dig> μm of both 5' and 3' primers and  <dig>  unit of platinum taq dna polymerase .
peltier thermal cycler  was used with a initial denaturation of  <dig> minutes at 94°, followed by  <dig> cycles of  <dig> seconds at 94°,  <dig> seconds at annealing temperature optimized between  <dig> and 60° for each pair of primers, and  <dig> minutes at 72°, completed with an additional  <dig> minutes at 72°.
each pcr fragment was sequenced with the amplification primers using bigdye terminator v <dig>  cycle sequencing ready reaction kit  and an automated abi prism®  <dig> genetic sequencer .
sequences were analysed and assembled with seqmerge .
snp discovery
snps were detected for each of the two parents for each species from heterozygous positions indicated by double-peak signatures in sequence chromatograms.
for each snp detected, haploid dna sequences from individual megagametophyte tissue were used as a control for paralogy.
because of their haploid nature, any double peak signature in the sequence chromatogram from a megagametophyte would indicate a snp resulting from polymorphism between gene loci, hence paralogous variation.
because these snps would result in fixed polymorphisms in the progeny, they were not considered for genotyping and mapping.
snp selection for bead array construction
gene regions were amplified by pcr and resequenced in order to identify in excess of  <dig>  candidate snps for each of white spruce and black spruce.
out of them,  <dig>  snps were used to construct two species-specific illumina bead arrays of  <dig> snps using the goldengate assay .
for each species, when more than one snp was available for a gene, snps were chosen so each snp was specific to one parental genotype.
this scheme allowed to validate the snp mapping approach, as snps of a same gene are expected to map at the same position on composite maps.
in optimizing the choice of snps for a given gene in a given species, those with maximum goldengate assay functionality score were chosen.
the functionality score is an a priori measure of snp adequation to the goldengate assay and takes into account a number of parameters, including sequence conformation around the snp, lack of repetitive elements in the surrounding sequence  and sequence redundancy against the available sequence database of the recipient species  <cit> .
for white spruce, a small subset of  <dig> snps representative of  <dig> genes was drawn from in silico identification of snps in contigs resulting from the assembly of est sequences  <cit> .
all chosen in silico snps involved est sequences from at least two different cdna clones and had a probability of occurrence of  <dig>  or more, according to the statistical assessment conducted with a bayesian method  <cit> .
snp genotyping assay
the illumina bead array technology was used to carry out all genotyping reactions in accordance with the manufacturer's protocol for the snp goldengate assay  <cit> .
highly multiplexed allele-specific extension reactions were conducted with two allele-specific primers per snp for each of the two species-specific 768-snp arrays using  <dig> ng of template dna per sample  for each progeny in each species and for positive controls consisting of five replicates of each parent of the mapping populations also used to identify snps from resequencing.
negative controls were also added to each 96-well sample plate.
ligation was completed with a third locus-specific primer.
this step was followed by pcr amplification on the extension-ligation product using primers labeled with either cy <dig> or cy <dig> dye to distinguish between alleles at each snp.
products of the pcr reaction were hybridized onto a decoded sentrix array matrix  .
bundles of the sam include beadtypes coated with oligonucleotides complementary to a primer address on the pcr product.
each beadtype is represented with an average redundancy of 30x on the array to optimize the accuracy of the final genotype signal.
following hybridization, the signal in each wavelength was determined using a bead array reader that converts images to intensity data.
the intensity data for each snp was normalized and assigned a cluster position  with the beadstudio software , and a quality score for each genotype was generated.
a gencall score cutoff of  <dig>  was used to determine valid genotypes at each snp and the snps retained had to get a minimum gentrain score of  <dig>   <cit> .
gentrain scores measure the reliability of snp detection based on the distribution of genotypic classes  <cit> .
dna reports, locus summaries, and final reports were generated for downstream analysis using the beadstudio software .
estimation of linkage maps
aflp, rapd, ssr, and estp marker data previously used for linkage mapping in the two species  <cit>  were considered together with the new snp data for constructing linkage maps.
for each cross, locus segregation was tested for goodness-of-fit to expected mendelian segregation ratios using chi-square tests .
distorted loci were excluded from further analyses.
linkage analyses were conducted with the male and female datasets independently to obtain two individual linkage maps for each species.
each snp was considered as an independent marker.
then, a composite linkage map was assembled from the two parental maps for each species, where snps from the same gene were considered simultaneously as a single haplotype.
individual and composite linkage maps were estimated using procedures described by pelgas et al.
<cit> .
both crosses were analysed using the "two-way pseudo-testcross" mapping approach  <cit> .
all linkage analyses and map estimations were performed with joinmap  <dig>   <cit> .
in addition, markers were ordered with the monte carlo maximum likelihood mapping algorithm implemented in joinmap  <dig>  and using standard parameters  <cit> .
with both versions of joinmap, the parameter cp  was used with a maximal threshold value of  <dig> for the jump, a ripple value of  <dig>  and kosambi's mapping function  <cit> .
for marker grouping and linked loci ordering, a lod score of  <dig>  to  <dig>  and a minimum recombination fraction  of  <dig>  were used.
the expected genome map length ge was estimated under the assumption of random marker distribution according to the formula of  <cit> .
an estimate of genome map coverage ce was obtained according to the formula of  <cit>  for the same lod value  as for previous coverage estimates  <cit> .
distribution of markers on linkage groups
for each species, randomness of gene locus distribution within and among linkage groups, heterogeneity of marker distribution  among linkage groups, and marker dispersion were analysed from the composite map as previously described  <cit> .
analyses were conducted by considering 1) all marker types and 2) only gene markers .
for aflp markers only, previously published analyses of randomness of distribution showed no aggregated pattern of distribution  <cit> .
as no additional aflp markers were added in the present study, no distribution test was conducted with this type of marker alone.
validation of marker orthology between species
the homoeology of linkage groups between white spruce and black spruce was determined according to the same criterion as described by pelgas et al.
<cit> .
the recognition of orthologous from paralogous loci was also carried out according to pelgas et al.
<cit> .
to validate exceptions to linkage group synteny, resequencing from haploid megagametophyte tissue was performed for presumed orthologous markers positioned on non-homoeologous linkage groups.
any sequence polymorphism detected in the chromatogram of the haploid dna sequence was considered as evidence for paralogy.
changes in synteny were validated on a second mapping population available for each species whenever necessary.
RESULTS
construction of snp-arrays
a total of  <dig> expressed sequences were amplified and resequenced in white spruce parents, of which  <dig> were found with at least one orthologous snp .
primers designed for amplification and resequencing of expressed genes in white spruce could be transferred to black spruce at a rate of  <dig> %, which is in line with results obtained previously for a more limited set of genes  <cit> .
the transfer procedure resulted in the amplification and sequencing of  <dig> genes in black spruce parents, of which  <dig> contained at least one orthologous snp .
for the  <dig> genes simultaneously amplified, resequenced, and found with at least one snp in each of white spruce and black spruce, only  <dig> % of the observed snps were shared between the two species .
this fraction is based on gene resequencing from two parents in each species.
accordingly, two illumina bead arrays of  <dig> snps for the goldengate assay were constructed, one for each species.
the array built for white spruce contained snps representative of  <dig> genes and that for black spruce contained snps for  <dig> genes.
the additional file  <dig>  provides primer sequences used for pcr amplification, unigene identifiers, and links with the foresttreedb database hosting the unigene sequences and their annotations  <cit> .
a total of  <dig> genes were represented simultaneously on both arrays .
all in all, the white spruce snp-array resulted in  <dig>  snp calls and the black spruce snp-array in  <dig>  snp calls.
<dig> in silico snps were detected in aligned ests derived from white spruce cdna libraries.
the resource is described in pavy et al.
<cit> .
for black spruce, all assayed snps were obtained after resequencing from genomic dna.
reproducibility of the snp assay and effect of template concentration
the reproducibility of the assay was evaluated using five replicates for each of the mapping parents also used to identify snps from resequencing.
when estimated over all valid snps , on average,  <dig> % of the snp calls were concordant with expectations when using the recommended amount of dna template .
the rate of concordant snp calls was lower when testing with less template dna .
overall success rate of the snp bead arrays
gentrain scores correspond to the reliability of snp detection based on the distribution of genotypic classes.
thus, it is a measure of reliability based on the total array of calls for a given snp.
according to illumina, for a snp to be retained, a minimum gentrain score of  <dig>  is advisable  <cit> .
in the present study, a snp had to get a minimum gentrain score of  <dig>  and had to be segregating in the related mapping population to be declared successful.
in white spruce,  <dig> % of snps identified by resequencing had gentrain score of  <dig>  or more , which is in the range of that obtained for human snps identified from resequencing for polymorphism discovery  <cit> .
in black spruce, the corresponding percentage was  <dig> %.
contrary to expectations, a number of these snps were monomorphic in the mapping populations .
it is likely that one of the two allele-specific primers in the goldengate assay defaulted for these snps.
when discarding these monomorphic snps, the overall rate of success for the genotyping of resequenced snps was  <dig> % in white spruce and  <dig> % in black spruce .
the genotyping success rate on the basis of the number of genes assayed was higher  because of redundancy of snp sampling for some genes.
<dig> in silico snps were detected in aligned est derived from white spruce cdna libraries.
the resource is described in pavy et al.
<cit> .
for black spruce, all snps were obtained after resequencing from genomic dna.
<dig> for snps with a gentrain score ≥  <dig> , valid gencall scores were obtained for  <dig> % of samples, on average .
for white spruce, a number of in silico snps identified from redundancy in est contigs  <cit>  were also included on the snp array, and  <dig> % of them had gentrain score of  <dig>  or more , which is comparable to the percentage obtained for resequenced snps.
the overall success rate taking into account segregation in the mapping population was lower at  <dig> %, because monomorphism for in silico snps was much more frequent than that for resequenced snps.
individuals previously used for est sequencing and in silico identification of snps did not include the parents of the present white spruce mapping population.
hence, it is likely that white spruce parents were homozygous for many of these in silico snps, which resulted in a much higher rate of monomorphism than for resequenced snps.
indeed, when scoring these in silico snps for individuals from natural populations, many of these snps exhibited the expected polymorphism and the overall success rate obtained was 60% .
rate of missing data according to snp gentrain scores
a missing data is generated when the gencall score for a particular individual and snp is below  <dig>   <cit> .
gencall is a measure of the reliability of an individual snp call relative to the distribution of genotypic classes for that snp.
the call rate, which is  <dig> minus the rate of missing data, could be estimated for all snps with acceptable gentrain scores.
in agreement with data on human snps  <cit> , our results indicated that snps with gentrain scores of  <dig>  or more were highly reliable with a low rate of missing data .
the average call rate per valid snp with gentrain score of  <dig>  or more was  <dig> % for white spruce and  <dig> % for black spruce.
for snps with gentrain scores higher than  <dig> , the rate of missing data became infinitesimal and the average rate of missing data per successful snp was very low, with an average of  <dig> % for white spruce and  <dig> % for black spruce , thus less than 1%.
<dig> average call rate is 100% minus average number of missing data in %.
snp success rate according to a priori snp functionality score
before construction of the snp bead array, a functionality score was calculated for each candidate snp using the illumina oligodesigner software  <cit> .
the functionality score relies much on the uniqueness and lack of repetitive elements in the surrounding sequence of each snp  <cit> .
the higher the score, the more likely the snp will be functional at the genotyping stage under the goldengate assay.
snps with a predicted functionality score above  <dig>  had a much higher rate of success than those below  <dig>   .
thus, in spite of the incompleteness of the spruce genome sequence used to estimate the appropriateness of candidate snps for the goldengate assay, the functionality score was still a valuable predictor of the likelihood of success of the designed oligonucleotides.
indeed, most of the snps selected for arraying had functionality scores equal or higher than  <dig>  .
1according to illumina oligodesigner software  <cit>
2excluding failed and monomorphic snps.
individual and composite linkage maps for white spruce and black spruce
depending on the cross and parent analysed, between  <dig> and  <dig> genetic markers were available to estimate each individual linkage map .
of these, between  <dig> % and 91% could be mapped.
they were distributed over  <dig> major linkage groups, except for parent  <dig>  for which one additional minor linkage group derived exclusively from aflps was obtained.
afor individual linkage maps, only markers segregating 1: <dig> or 1:1:1: <dig> were used.
bbonferroni correction: p ≤  <dig> /number of loci.
cnumber of linkage group composed of  <dig> sub-groups having more than  <dig> markers.
dcould not be calculated due to the merging of data.
for each cross, male and female datasets were integrated into one composite linkage map representative of each species.
for white spruce,  <dig> markers  could be mapped over the  <dig> .
<dig> cm, including nine accessory marker loci .
on this map, the average marker spacing was  <dig>  cm.
for black spruce, a total of  <dig> markers  could be mapped over the  <dig> .
<dig> cm, including  <dig> accessory marker loci .
the average marker spacing was  <dig>  cm.
marker distribution
to analyse the distribution of marker loci, g-tests for goodness-of-fit were used.
for each composite linkage map, markers were homogeneously distributed across linkage groups .
therefore, coefficients of dispersion could be estimated for each composite map see  <cit> .
whether we considered all markers or gene loci only , the coefficients of dispersion calculated with a sliding window of  <dig> cm ranged between  <dig>  and  <dig>  for the white spruce composite map and between  <dig>  and  <dig>  for the black spruce composite map.
such values around  <dig>  indicates the absence of significant marker islands and a random distribution of gene loci at the present mapping intensity.
interspecific comparisons
the composite linkage maps of white spruce and black spruce shared  <dig> homologous anchor markers .
over the  <dig> homoeologous linkage groups,  <dig> % of homologous markers were in synteny .
one estp locus  and three gene snps  were involved in synteny discrepancies.
the differential positioning of the estp locus was already pointed out and validated by sequencing  <cit> .
the examination of chromatograms obtained from sequencing haploid megagametophyte dna around the two snp markers  <dig> and  <dig> respectively positioned onto lgs ii/vii and lgs xii/iii of white spruce/black spruce revealed polymorphisms.
such evidence indicates that these two anchor markers differentially positioned in white spruce and black spruce correspond to two paralogous gene loci.
another discrepancy involved snp  <dig> positioned onto white spruce lg v and on black spruce lg vi.
however, linkage mapping analyses independently conducted in a second white spruce mapping population  led to localize this gene marker on lg vi in both species.
therefore, exceptions to synteny were rejected regarding the differential positioning of these three snp gene loci between the two species.
along with synteny, macrocolinearity was also well conserved among homoeologous linkage groups between the two spruce taxa.
on average, 82% of syntenic anchor markers were positioned in the same order .
minor inversions involving closely spaced markers were observed within each homoeologous linkage group, involving a total of  <dig> ssr,  <dig> estps and  <dig> snps.
these inversions were also detected when comparing individual linkage maps within species , thus lending support to statistical artifacts related to joining individual linkage maps when constructing composite maps.
discussion
applicability of the goldengate snp assay
the illumina goldengate snp assay, together with the bead array technology, has been extensively used in studies of human polymorphisms  <cit> .
its use has also been extended to animal genetics, especially with regard to efforts undertaken by the bovine hapmap consortium  <cit> .
because of its well-established reliability with human data and high level of multiplexing, there is a growing interest in using the goldengate snp assay in plants when large numbers of snps need to be surveyed.
indeed, the flexibility of its design and the large number of snps screened per assay make the technology appropriate for genome scan applications or plant molecular breeding purposes.
however, with the exception of reports from biotech companies, very few studies with detailed supportive evidence have been published on the ease of implementing the technology in plants and non model organisms with largely unsequenced genomes.
we have considered this issue in spruce by examining the genotyping success rates obtained over a set of high quality snps derived from independent resequencing.
call rate
the call rate is the fraction of genotypes called over the possible snp calls after having excluded unsuccessful assays  <cit> .
in our data, the snps retained for genetic mapping analysis had call rates greater than 96%.
their average call rate was in excess of 99%.
this value is comparable to call rates obtained with human snps using the goldengate assay, which are near 100% when following the same stringent criterion of minimum gentrain score of  <dig>  as that used in the present study for considering as valid the genotyping of a snp  <cit> .
snp conversion rate
the rate of snps successfully genotyped using the goldengate assay, or snp conversion rate, is calculated by counting only valid snps displaying gentrain score above a given cutoff.
the gentrain score of a snp reflects the degree of separation between homozygote and heterozygote clusters and the ease of placement of individual calls within a cluster, which are key measures of signal-to-noise in the assay data  <cit> .
as recommended by illumina, we followed a conservative approach and filtered out snps with a gentrain score below  <dig> , which tend to show low call rate and hence, high rates of missing data  <cit> .
based on this criterion alone and not considering monomorphic snps, the snp conversion rates were  <dig> % and  <dig> % of the resequenced snps for the white spruce and black spruce snp arrays, respectively.
these rates were down to  <dig> % and  <dig> %, respectively, when excluding monomorphic snps.
high snp conversion rates have been obtained in studies conducted in species with completely sequenced genomes such as for the human genome.
besides studies involving human snps, gentrain score cutoffs have been barely mentioned in the literature although they affect the snp conversion rates and the average call rate.
thus, the following comparisons of our results with the published literature may not be orthogonal.
in a study encompassing  <dig>  resequenced snps derived from the human genome,  <dig> % of snps were called after removal of low frequency snps  <cit> .
two other sets each encompassing  <dig> hapmap snps were successfully assayed with  <dig> % and  <dig> % conversion rates  <cit> .
the application of the goldengate snp assay in non model species resulted in slightly lower snp conversion rates than that obtained for human snps.
the conversion rate was  <dig> % over  <dig>  resequenced snps in barley  <cit>  and  <dig> % over a panel of  <dig> bovine snps  <cit> .
in boechera stricta, a species from the cruciferae family, a conversion rate of  <dig> % was obtained over  <dig> resequenced snps by avoiding highly similar sequences such as for members of complex gene families  <cit> .
for candidate genes that were members of large gene families, the authors searched for markers in flanking genes that were single copy in arabidopsis  <cit> .
under these circumstances, the design of illumina probes was efficient and specific.
however, such an approach is limited to plant species close to arabidopsis, which genome is completely sequenced.
two main factors may explain the lower snp conversion rate obtained with spruces as compared to other species analysed to date with the goldengate snp assay.
first, we have adopted the severe criterion of gentrain scores <  <dig>  to reject snps, while this factor remains unknown for most of the non human studies submentioned.
in the present study, such a conservative criterion translated in snps with a high call rate, which was necessary for accurate gene mapping.
second, the complexity of conifer genomes e.g.
<cit>  may obstruct the development of specific probes for the assay.
indeed, our snp assays incorporated mostly sequences belonging to multigene families including many transcription factors .
with the present incomplete knowledge of conifer genomes, it was not possible to take into account gene family structures to improve the oligonucleotide design for the goldengate snp assay.
the level of duplication has not been quantified yet in conifer genomes contrary to angiosperm model species such as legumes or grasses.
however, phylogenetic analysis of multigene families in conifers has revealed an organization different from that observed in angiosperms.
examples include the adh genes  <cit> , knox-i  <cit>  and myb  <cit>  regulatory genes, with many gene duplications at least recent enough not to be shared with angiosperms.
paralogous snps generate a background signal giving rise to cluster compression and therefore, decreasing gentrain scores and the snp conversion rate, unless manual editing of the clusters is used to eliminate all snps that do no cluster well  <cit> .
if the conifer genome is highly duplicated, as suspected, interpreting goldengate snp data in the context of maximizing the conversion rate could be very challenging and imply much lower call rates per snp recovered.
without an exhaustive knowledge of gene sequences within a gene family, there is no easy way to avoid such snps with potential lack of specificity of flanking sequences.
this trend is even more likely, given that all snps tested herein were in gene sequences.
indeed, the probability for these snps to represent paralogous variation across family members is likely higher than that for snps located in non coding dna.
for snps exhibiting monomorphism , it is likely that one of the two allele-specific primers in the goldengate assay defaulted for these snps, given that resequencing was performed to discover snps and that we are confident that they were not sequencing artifacts.
these failures might be related to the same factors as above.
thus, considering the fragmentary knowledge we have of conifer genomes and the emphasis put on sampling snps from regulatory genes from large gene families, the overall success rate obtained was decent.
the success rate was also repeatable between the two species-specific snp-arrays that we have independently tested.
as the genome of conifers is becoming better known at the sequence level, due to several large-scale est and bac sequencing projects, the rates of genotyping success for new snps are expected to be even higher in the near future.
transcript linkage maps
with around  <dig>  to  <dig>  mb  <cit> , spruce and conifer genomes are more than  <dig> times larger than that of arabidopsis and three times larger than the human genome.
for the largely unsequenced conifer genomes, sequencing their coding regions through est sequencing and gene resequencing currently represents the most efficient approach to comparative and structural genomics using gene linkage maps  <cit> .
the snp genotyping assays developed in this study enabled to map hundreds of expressed genes in the conifer picea, which represent a large increase over any previous gene mapping effort in gymnosperms e.g.
<cit> .
a total of  <dig> major linkage groups were delineated in each species composite map, which is in agreement with cytogenetic studies indicating similar karyotypes and same numbers of chromosomes  for both white spruce and black spruce and most other pinaceae  <cit> .
because of the agreement between number of linkage groups and number of chromosomes, the present maps could be considered as saturated .
high genome coverage values also point to this observation.
however, a direct pairing between linkage groups and chromosomes was not possible, as picea and most conifer chromosomes are difficult to differentiate based on standard cytological techniques  <cit> .
with the addition of several hundreds new gene markers, the observed total length of individual parental maps was increased by 15% to 25%, and the genome coverage values were increased by nearly 10% compared to previous maps mostly based on anonymous markers  <cit> .
with more than  <dig> gene anchor markers positioned on each parental map, marker density was more than 50% higher and gene density was an order of magnitude higher than previously obtained  <cit> .
as compared to aflp and other types of dominant markers, these codominant snp gene markers mapped at such high density will contribute towards improving qtl mapping precision and power  <cit> .
more than  <dig> gene loci were shared between the two composite maps developed herein, most of which being derived from the new gene snp developing effort.
while previous efforts to increase the number of gene or ssr anchor markers relied on using several crosses per species e.g.
<cit> , the ease of mapping large numbers of transcripts using snps relaxes the need to implicate more than one cross for increasing the number of mapped anchor markers.
however, using an additional cross may be highly useful to validate orthology of gene loci between species when breaks of synteny are observed .
interspecific comparisons
comparisons between the composite maps of white spruce and black spruce revealed high synteny and colinearity between their  <dig> homoeologous linkage groups, in spite of the divergence of their lineages more than  <dig> myr ago  <cit> .
while previous reports of genome comparisons between the two species reached similar conclusions about synteny and colinearity  <cit> , the present observations are based on a more than four-fold increase in the number of mapped anchor loci in common between the two species.
synteny between genomic regions can only be established if markers are true orthologs  <cit> .
the breakdown in synteny previously noted between white spruce and black spruce for lg iii and lg iv  <cit>  was confirmed in the present study.
in addition, three new cases of putative inter-chromosomal translocations between both species were observed, but turned out to be false positives, after checks for locus orthology using haploid megagametophyte dna sequences or validating gene marker positions in a second cross of white spruce.
thus, these three gene snps corresponded to paralogous loci between white spruce and black spruce.
such false positives are likely to be frequent in conifer mapping studies, as previously observed  <cit> .
while synteny was well conserved, exceptions to colinearity between the two genomes were observed.
small inversions between the two species composite maps involving a few closely spaced anchor loci were noted for all linkage groups, but no translocation within linkage groups was observed.
most likely, these inversions resulted from analytical artifacts since the same discrepancies were observed between individual linkage maps within species.
thus, it is safe to assume that these inversions resulted from the integration of both parental datasets for each species rather than from true chromosomal inversions.
such a trend has also been observed when integrating individual linkage maps in other species  <cit> .
CONCLUSIONS
the present report illustrates how new highly multiplexed snp genotyping approaches can be used to accelerate the structural analyses of complex and largely unknown genomes.
the present maps are the most advanced genetic maps for spruce with regard to gene density and will open up several opportunities.
these maps are currently used for qtl and eqtl detection and their next more densified versions will help target specific regions of the genome for future bac sequencing  <cit> .
the present maps will also be used to explore in more detail the distribution of multigene families over linkage groups and hence, the organization of the spruce genome and its evolution.
such an approach, though at a small scale, has been previously used for the knox i multigene family, underlining a major chromosomal duplication in the spruce genome  <cit> .
maps with increasing gene densities are also required across the conifers to extend comparative genome studies and decipher the evolution of genome structure.
such comparisons have already been undertaken based on a limited number of anchor markers e.g.
<cit> .
however, comparative genome studies in these species remain challenging since snps may have to be identified de novo in each species.
this trend is best exemplified by the small number of gene snps shared between mapping populations belonging to different congeneric taxa: only  <dig> snps were shared between white spruce and black spruce mapping populations out of  <dig>  candidate gene snps discovered in both species.
using more crosses per species could contribute towards increasing the number of shared snps between species, but it is unlikely to change the order of magnitude of this number.
on a more optimistic note, the primer transfer rate between white spruce and black spruce was high, which might facilitate resequencing efforts towards gene snp discovery in other spruce species.
however, as species from different conifer genera are targeted, the primer transfer rate is expected to diminish  <cit> .
identifying sets of homologous transcripts in other pinaceae has improved in the recent years through efficient mining of sequence databases e.g.
<cit>  and because of large sets of est sequences for several species  <cit> .
nevertheless, when changes in chromosomal structure are observed between species based on gene linkage maps, gene orthology must be verified.
distinguishing between orthology and paralogy can be precarious when the assessment is only based on homology scores, and without an adequate phylogenetic landscape for each gene family involved  <cit> .
in such conditions, a sensitive analysis at the sequence level is warranted before declaring such structural changes  <cit> .
thus, comparative mapping studies between conifer species call for further developments of gene maps and gene sequence collections, but also for genomics and bioinformatics tools enabling to assess more efficiently orthology relationships.
authors' contributions
np, jb: genotyping technology evaluation; bp, ig, ni: genetic mapping analyses; sbe, sbl, fg: resequencing, snp discovery, design of the snp arrays, and snp data analyses; np, ml, candidate gene identification; np, bp, jb: manuscript preparation; jb, ni: project design, funding, and overall supervision.
all authors have read and approved the final manuscript.
supplementary material
additional file 1
data about the snps.
this table provide references about the genes analysed: unigeneid in foresttreedb, annotations, marker names, gene names, links to annotation page in foresttreedb, primer sequences used to amplify sequences encompassing the unigenes.
click here for file
acknowledgements
we thank alexandre montpetit and his team at the genome québec innovation centre  for assistance with the acquisition of snp data on the illumina platform.
this work was funded by grants from genome canada, genome québec, and the canadian biotechnology strategy through the arborea project  <cit>  lead by j. mackay and j. bousquet.

 single nucleotide polymorphims have become a genomic commodity as they are becoming indispensable in various genome scans aimed at mapping genomes, finding associations with complex traits, and population genomics. they are distributed along the various regions of the genomes and are frequent in coding regions of angiosperms and conifers. however, the efficiency of genome scans is not only dependant on a wide genomic distribution of snps. indeed, it also relies on the ability to genotype large numbers of snps over large sets of individuals. snp genome scans in non model species usually involve two steps: the discovery of snps and genotyping. with no a priori knowledge of dna polymorphisms, snps are usually discovered through various strategies of individual or pool dna sequencing, or by using tilling techniques, a high-throughput strategy relying on the enzymatic cleavage of mismatches. for a number of crop species, current resequencing efforts have led to the development of snp databases and generate a wealth of snps usable in genome scans. in conifers, large-scale est sequencing projects have been initiated, providing a starting point to develop snp resources in pine and spruce. several snp genotyping array approaches have been developed with variable success. the accuracy of innovative snp genotyping technologies has been assessed mostly through the development of assays suitable for analysing variations in the human genome. broadly speaking, four reaction principles govern snp genotyping assays: hybridization with allele-specific oligonucleotide probes, oligonucleotide ligation, single nucleotide primer extension, and enzymatic cleavage reviewed in. among these approaches, the goldengate assay developed by illumina and relying on the bead array technology has demonstrated high performance with high levels of call rate, reproducibility, and overall success rate for the analysis of the human genome. high-throughput snp assays have recently been applied to plants. large datasets of snp-based markers are being developed in barley through the development of genotyping assays relying on illumina's technologies, leading to the undertaking of an international snp project. the same genotyping approach has made it possible to map large datasets of snps even in complex and duplicated genomes such as soybean, and projects are underway in hexaploid wheat and poplar. in the present study, we are asking whether high-throughput snp genotyping technologies developed for human population genomics applications, such as the illumina goldengate snp assay, are applicable to large and essentially unsequenced genomes as seen in conifers. conifer genomes reach very large sizes, around – mb, consisting mostly of repetitive sequences. for the two conifers considered herein, white spruce and black spruce, genome sizes are well in excess of 10e bp. moreover, the partial knowledge of the large and redundant genomes of conifers can be a limiting factor to design an efficient snp genotyping assay. indeed, sequences located upstream and downstream the snp cannot be fully validated for locus specificity and the possible presence of repetitive elements. the possible effect of such a drawback remains to be verified for most crop and tree species which genomes are essentially not sequenced. based on est sequence data available for white spruce, we have designed primers and resequenced genomic dna for hundreds of genes in white spruce and black spruce. the high quality snp datasets developed were used to select snps amenable to the goldengate genotyping assay and test the technology for these two species. then, we integrated these snp data into linkage maps of expressed genes and illustrated the possibility to rapidly improve the density of existing genetic maps for spruce species. methods plant material for genotyping assays and linkage mapping plant material consisted of one outbred f cross # c for white spruce with progeny derived from two parents selected for their high level of heterozygosity for estp anchor markers and for their intermediate performance for a number of traits such as somatic embryogenic capacity. plant material for black spruce was representative of the hybridizing species complex picea mariana × picea rubens and consisted of progeny derived from the backcross-like cross bc, as previously described. spruces all harbor 2n = 2x = chromosomes, as for most pinaceae. the lineages leading to these two taxa are thought to have diverged at least myr ago. gene selection and pcr amplification for snp discovery, we considered a collection of unigenes derived from a white spruce est database. for each gene, coding regions were identified based on alignments with similar sequences from uniprot-swissprot protein databases. pcr primers for amplification and resequencing were designed using primer. whenever possible, one of the primers was anchored outside of the coding regions in order to increase amplification specificity. the same primers were also used for snp discovery in black spruce. the lists of unigene and primer sequences for the genes represented on the two illumina snp bead arrays for white spruce and black spruce are provided. for each of white spruce and black spruce, dna was isolated from the two parents of the mapping population and their progeny, as well as from haploid megagametophyte tissue using a dneasy plant mini kit. about 5– ng of template dna were used for pcr amplification. reactions were done in μl containing mm tris-hcl, mm kcl, – mm mgcl μm of each dntp, μm of both 5' and 3' primers and unit of platinum taq dna polymerase. peltier thermal cycler was used with a initial denaturation of minutes at 94°, followed by cycles of seconds at 94°, seconds at annealing temperature optimized between and 60° for each pair of primers, and minutes at 72°, completed with an additional minutes at 72°. each pcr fragment was sequenced with the amplification primers using bigdye terminator v cycle sequencing ready reaction kit and an automated abi prism® genetic sequencer. sequences were analysed and assembled with seqmerge. snp discovery snps were detected for each of the two parents for each species from heterozygous positions indicated by double-peak signatures in sequence chromatograms. for each snp detected, haploid dna sequences from individual megagametophyte tissue were used as a control for paralogy. because of their haploid nature, any double peak signature in the sequence chromatogram from a megagametophyte would indicate a snp resulting from polymorphism between gene loci, hence paralogous variation. because these snps would result in fixed polymorphisms in the progeny, they were not considered for genotyping and mapping. snp selection for bead array construction gene regions were amplified by pcr and resequenced in order to identify in excess of candidate snps for each of white spruce and black spruce. out of them, snps were used to construct two species-specific illumina bead arrays of snps using the goldengate assay. for each species, when more than one snp was available for a gene, snps were chosen so each snp was specific to one parental genotype. this scheme allowed to validate the snp mapping approach, as snps of a same gene are expected to map at the same position on composite maps. in optimizing the choice of snps for a given gene in a given species, those with maximum goldengate assay functionality score were chosen. the functionality score is an a priori measure of snp adequation to the goldengate assay and takes into account a number of parameters, including sequence conformation around the snp, lack of repetitive elements in the surrounding sequence and sequence redundancy against the available sequence database of the recipient species. for white spruce, a small subset of snps representative of genes was drawn from in silico identification of snps in contigs resulting from the assembly of est sequences. all chosen in silico snps involved est sequences from at least two different cdna clones and had a probability of occurrence of or more, according to the statistical assessment conducted with a bayesian method. snp genotyping assay the illumina bead array technology was used to carry out all genotyping reactions in accordance with the manufacturer's protocol for the snp goldengate assay. highly multiplexed allele-specific extension reactions were conducted with two allele-specific primers per snp for each of the two species-specific 768-snp arrays using ng of template dna per sample for each progeny in each species and for positive controls consisting of five replicates of each parent of the mapping populations also used to identify snps from resequencing. negative controls were also added to each 96-well sample plate. ligation was completed with a third locus-specific primer. this step was followed by pcr amplification on the extension-ligation product using primers labeled with either cy or cy dye to distinguish between alleles at each snp. products of the pcr reaction were hybridized onto a decoded sentrix array matrix. bundles of the sam include beadtypes coated with oligonucleotides complementary to a primer address on the pcr product. each beadtype is represented with an average redundancy of 30x on the array to optimize the accuracy of the final genotype signal. following hybridization, the signal in each wavelength was determined using a bead array reader that converts images to intensity data. the intensity data for each snp was normalized and assigned a cluster position with the beadstudio software, and a quality score for each genotype was generated. a gencall score cutoff of was used to determine valid genotypes at each snp and the snps retained had to get a minimum gentrain score of. gentrain scores measure the reliability of snp detection based on the distribution of genotypic classes. dna reports, locus summaries, and final reports were generated for downstream analysis using the beadstudio software. estimation of linkage maps aflp, rapd, ssr, and estp marker data previously used for linkage mapping in the two species were considered together with the new snp data for constructing linkage maps. for each cross, locus segregation was tested for goodness-of-fit to expected mendelian segregation ratios using chi-square tests. distorted loci were excluded from further analyses. linkage analyses were conducted with the male and female datasets independently to obtain two individual linkage maps for each species. each snp was considered as an independent marker. then, a composite linkage map was assembled from the two parental maps for each species, where snps from the same gene were considered simultaneously as a single haplotype. individual and composite linkage maps were estimated using procedures described by pelgas et al.. both crosses were analysed using the "two-way pseudo-testcross" mapping approach. all linkage analyses and map estimations were performed with joinmap. in addition, markers were ordered with the monte carlo maximum likelihood mapping algorithm implemented in joinmap and using standard parameters. with both versions of joinmap, the parameter cp was used with a maximal threshold value of for the jump, a ripple value of and kosambi's mapping function. for marker grouping and linked loci ordering, a lod score of to and a minimum recombination fraction of were used. the expected genome map length ge was estimated under the assumption of random marker distribution according to the formula of. an estimate of genome map coverage ce was obtained according to the formula of for the same lod value as for previous coverage estimates. distribution of markers on linkage groups for each species, randomness of gene locus distribution within and among linkage groups, heterogeneity of marker distribution among linkage groups, and marker dispersion were analysed from the composite map as previously described. analyses were conducted by considering 1) all marker types and 2) only gene markers. for aflp markers only, previously published analyses of randomness of distribution showed no aggregated pattern of distribution. as no additional aflp markers were added in the present study, no distribution test was conducted with this type of marker alone. validation of marker orthology between species the homoeology of linkage groups between white spruce and black spruce was determined according to the same criterion as described by pelgas et al.. the recognition of orthologous from paralogous loci was also carried out according to pelgas et al.. to validate exceptions to linkage group synteny, resequencing from haploid megagametophyte tissue was performed for presumed orthologous markers positioned on non-homoeologous linkage groups. any sequence polymorphism detected in the chromatogram of the haploid dna sequence was considered as evidence for paralogy. changes in synteny were validated on a second mapping population available for each species whenever necessary.  construction of snp-arrays a total of expressed sequences were amplified and resequenced in white spruce parents, of which were found with at least one orthologous snp. primers designed for amplification and resequencing of expressed genes in white spruce could be transferred to black spruce at a rate of %, which is in line with results obtained previously for a more limited set of genes. the transfer procedure resulted in the amplification and sequencing of genes in black spruce parents, of which contained at least one orthologous snp. for the genes simultaneously amplified, resequenced, and found with at least one snp in each of white spruce and black spruce, only % of the observed snps were shared between the two species. this fraction is based on gene resequencing from two parents in each species. accordingly, two illumina bead arrays of snps for the goldengate assay were constructed, one for each species. the array built for white spruce contained snps representative of genes and that for black spruce contained snps for genes. the additional file provides primer sequences used for pcr amplification, unigene identifiers, and links with the foresttreedb database hosting the unigene sequences and their annotations. a total of genes were represented simultaneously on both arrays. all in all, the white spruce snp-array resulted in snp calls and the black spruce snp-array in snp calls. in silico snps were detected in aligned ests derived from white spruce cdna libraries. the resource is described in pavy et al.. for black spruce, all assayed snps were obtained after resequencing from genomic dna. reproducibility of the snp assay and effect of template concentration the reproducibility of the assay was evaluated using five replicates for each of the mapping parents also used to identify snps from resequencing. when estimated over all valid snps, on average, % of the snp calls were concordant with expectations when using the recommended amount of dna template. the rate of concordant snp calls was lower when testing with less template dna. overall success rate of the snp bead arrays gentrain scores correspond to the reliability of snp detection based on the distribution of genotypic classes. thus, it is a measure of reliability based on the total array of calls for a given snp. according to illumina, for a snp to be retained, a minimum gentrain score of is advisable. in the present study, a snp had to get a minimum gentrain score of and had to be segregating in the related mapping population to be declared successful. in white spruce, % of snps identified by resequencing had gentrain score of or more, which is in the range of that obtained for human snps identified from resequencing for polymorphism discovery. in black spruce, the corresponding percentage was %. contrary to expectations, a number of these snps were monomorphic in the mapping populations. it is likely that one of the two allele-specific primers in the goldengate assay defaulted for these snps. when discarding these monomorphic snps, the overall rate of success for the genotyping of resequenced snps was % in white spruce and % in black spruce. the genotyping success rate on the basis of the number of genes assayed was higher because of redundancy of snp sampling for some genes. in silico snps were detected in aligned est derived from white spruce cdna libraries. the resource is described in pavy et al.. for black spruce, all snps were obtained after resequencing from genomic dna. for snps with a gentrain score ≥, valid gencall scores were obtained for % of samples, on average. for white spruce, a number of in silico snps identified from redundancy in est contigs were also included on the snp array, and % of them had gentrain score of or more, which is comparable to the percentage obtained for resequenced snps. the overall success rate taking into account segregation in the mapping population was lower at %, because monomorphism for in silico snps was much more frequent than that for resequenced snps. individuals previously used for est sequencing and in silico identification of snps did not include the parents of the present white spruce mapping population. hence, it is likely that white spruce parents were homozygous for many of these in silico snps, which resulted in a much higher rate of monomorphism than for resequenced snps. indeed, when scoring these in silico snps for individuals from natural populations, many of these snps exhibited the expected polymorphism and the overall success rate obtained was 60%. rate of missing data according to snp gentrain scores a missing data is generated when the gencall score for a particular individual and snp is below. gencall is a measure of the reliability of an individual snp call relative to the distribution of genotypic classes for that snp. the call rate, which is minus the rate of missing data, could be estimated for all snps with acceptable gentrain scores. in agreement with data on human snps, our results indicated that snps with gentrain scores of or more were highly reliable with a low rate of missing data. the average call rate per valid snp with gentrain score of or more was % for white spruce and % for black spruce. for snps with gentrain scores higher than, the rate of missing data became infinitesimal and the average rate of missing data per successful snp was very low, with an average of % for white spruce and % for black spruce, thus less than 1%. average call rate is 100% minus average number of missing data in %. snp success rate according to a priori snp functionality score before construction of the snp bead array, a functionality score was calculated for each candidate snp using the illumina oligodesigner software. the functionality score relies much on the uniqueness and lack of repetitive elements in the surrounding sequence of each snp. the higher the score, the more likely the snp will be functional at the genotyping stage under the goldengate assay. snps with a predicted functionality score above had a much higher rate of success than those below. thus, in spite of the incompleteness of the spruce genome sequence used to estimate the appropriateness of candidate snps for the goldengate assay, the functionality score was still a valuable predictor of the likelihood of success of the designed oligonucleotides. indeed, most of the snps selected for arraying had functionality scores equal or higher than. 1according to illumina oligodesigner software 2excluding failed and monomorphic snps. individual and composite linkage maps for white spruce and black spruce depending on the cross and parent analysed, between and genetic markers were available to estimate each individual linkage map. of these, between % and 91% could be mapped. they were distributed over major linkage groups, except for parent for which one additional minor linkage group derived exclusively from aflps was obtained. afor individual linkage maps, only markers segregating 1: or 1:1:1: were used. bbonferroni correction: p ≤ /number of loci. cnumber of linkage group composed of sub-groups having more than markers. dcould not be calculated due to the merging of data. for each cross, male and female datasets were integrated into one composite linkage map representative of each species. for white spruce, markers could be mapped over the. cm, including nine accessory marker loci. on this map, the average marker spacing was cm. for black spruce, a total of markers could be mapped over the. cm, including accessory marker loci. the average marker spacing was cm. marker distribution to analyse the distribution of marker loci, g-tests for goodness-of-fit were used. for each composite linkage map, markers were homogeneously distributed across linkage groups. therefore, coefficients of dispersion could be estimated for each composite map see. whether we considered all markers or gene loci only, the coefficients of dispersion calculated with a sliding window of cm ranged between and for the white spruce composite map and between and for the black spruce composite map. such values around indicates the absence of significant marker islands and a random distribution of gene loci at the present mapping intensity. interspecific comparisons the composite linkage maps of white spruce and black spruce shared homologous anchor markers. over the homoeologous linkage groups, % of homologous markers were in synteny. one estp locus and three gene snps were involved in synteny discrepancies. the differential positioning of the estp locus was already pointed out and validated by sequencing. the examination of chromatograms obtained from sequencing haploid megagametophyte dna around the two snp markers and respectively positioned onto lgs ii/vii and lgs xii/iii of white spruce/black spruce revealed polymorphisms. such evidence indicates that these two anchor markers differentially positioned in white spruce and black spruce correspond to two paralogous gene loci. another discrepancy involved snp positioned onto white spruce lg v and on black spruce lg vi. however, linkage mapping analyses independently conducted in a second white spruce mapping population led to localize this gene marker on lg vi in both species. therefore, exceptions to synteny were rejected regarding the differential positioning of these three snp gene loci between the two species. along with synteny, macrocolinearity was also well conserved among homoeologous linkage groups between the two spruce taxa. on average, 82% of syntenic anchor markers were positioned in the same order. minor inversions involving closely spaced markers were observed within each homoeologous linkage group, involving a total of ssr, estps and snps. these inversions were also detected when comparing individual linkage maps within species, thus lending support to statistical artifacts related to joining individual linkage maps when constructing composite maps. discussion applicability of the goldengate snp assay the illumina goldengate snp assay, together with the bead array technology, has been extensively used in studies of human polymorphisms. its use has also been extended to animal genetics, especially with regard to efforts undertaken by the bovine hapmap consortium. because of its well-established reliability with human data and high level of multiplexing, there is a growing interest in using the goldengate snp assay in plants when large numbers of snps need to be surveyed. indeed, the flexibility of its design and the large number of snps screened per assay make the technology appropriate for genome scan applications or plant molecular breeding purposes. however, with the exception of reports from biotech companies, very few studies with detailed supportive evidence have been published on the ease of implementing the technology in plants and non model organisms with largely unsequenced genomes. we have considered this issue in spruce by examining the genotyping success rates obtained over a set of high quality snps derived from independent resequencing. call rate the call rate is the fraction of genotypes called over the possible snp calls after having excluded unsuccessful assays. in our data, the snps retained for genetic mapping analysis had call rates greater than 96%. their average call rate was in excess of 99%. this value is comparable to call rates obtained with human snps using the goldengate assay, which are near 100% when following the same stringent criterion of minimum gentrain score of as that used in the present study for considering as valid the genotyping of a snp. snp conversion rate the rate of snps successfully genotyped using the goldengate assay, or snp conversion rate, is calculated by counting only valid snps displaying gentrain score above a given cutoff. the gentrain score of a snp reflects the degree of separation between homozygote and heterozygote clusters and the ease of placement of individual calls within a cluster, which are key measures of signal-to-noise in the assay data. as recommended by illumina, we followed a conservative approach and filtered out snps with a gentrain score below, which tend to show low call rate and hence, high rates of missing data. based on this criterion alone and not considering monomorphic snps, the snp conversion rates were % and % of the resequenced snps for the white spruce and black spruce snp arrays, respectively. these rates were down to % and %, respectively, when excluding monomorphic snps. high snp conversion rates have been obtained in studies conducted in species with completely sequenced genomes such as for the human genome. besides studies involving human snps, gentrain score cutoffs have been barely mentioned in the literature although they affect the snp conversion rates and the average call rate. thus, the following comparisons of our results with the published literature may not be orthogonal. in a study encompassing resequenced snps derived from the human genome, % of snps were called after removal of low frequency snps. two other sets each encompassing hapmap snps were successfully assayed with % and % conversion rates. the application of the goldengate snp assay in non model species resulted in slightly lower snp conversion rates than that obtained for human snps. the conversion rate was % over resequenced snps in barley and % over a panel of bovine snps. in boechera stricta, a species from the cruciferae family, a conversion rate of % was obtained over resequenced snps by avoiding highly similar sequences such as for members of complex gene families. for candidate genes that were members of large gene families, the authors searched for markers in flanking genes that were single copy in arabidopsis. under these circumstances, the design of illumina probes was efficient and specific. however, such an approach is limited to plant species close to arabidopsis, which genome is completely sequenced. two main factors may explain the lower snp conversion rate obtained with spruces as compared to other species analysed to date with the goldengate snp assay. first, we have adopted the severe criterion of gentrain scores < to reject snps, while this factor remains unknown for most of the non human studies submentioned. in the present study, such a conservative criterion translated in snps with a high call rate, which was necessary for accurate gene mapping. second, the complexity of conifer genomes e.g. may obstruct the development of specific probes for the assay. indeed, our snp assays incorporated mostly sequences belonging to multigene families including many transcription factors. with the present incomplete knowledge of conifer genomes, it was not possible to take into account gene family structures to improve the oligonucleotide design for the goldengate snp assay. the level of duplication has not been quantified yet in conifer genomes contrary to angiosperm model species such as legumes or grasses. however, phylogenetic analysis of multigene families in conifers has revealed an organization different from that observed in angiosperms. examples include the adh genes, knox-i and myb regulatory genes, with many gene duplications at least recent enough not to be shared with angiosperms. paralogous snps generate a background signal giving rise to cluster compression and therefore, decreasing gentrain scores and the snp conversion rate, unless manual editing of the clusters is used to eliminate all snps that do no cluster well. if the conifer genome is highly duplicated, as suspected, interpreting goldengate snp data in the context of maximizing the conversion rate could be very challenging and imply much lower call rates per snp recovered. without an exhaustive knowledge of gene sequences within a gene family, there is no easy way to avoid such snps with potential lack of specificity of flanking sequences. this trend is even more likely, given that all snps tested herein were in gene sequences. indeed, the probability for these snps to represent paralogous variation across family members is likely higher than that for snps located in non coding dna. for snps exhibiting monomorphism, it is likely that one of the two allele-specific primers in the goldengate assay defaulted for these snps, given that resequencing was performed to discover snps and that we are confident that they were not sequencing artifacts. these failures might be related to the same factors as above. thus, considering the fragmentary knowledge we have of conifer genomes and the emphasis put on sampling snps from regulatory genes from large gene families, the overall success rate obtained was decent. the success rate was also repeatable between the two species-specific snp-arrays that we have independently tested. as the genome of conifers is becoming better known at the sequence level, due to several large-scale est and bac sequencing projects, the rates of genotyping success for new snps are expected to be even higher in the near future. transcript linkage maps with around to mb, spruce and conifer genomes are more than times larger than that of arabidopsis and three times larger than the human genome. for the largely unsequenced conifer genomes, sequencing their coding regions through est sequencing and gene resequencing currently represents the most efficient approach to comparative and structural genomics using gene linkage maps. the snp genotyping assays developed in this study enabled to map hundreds of expressed genes in the conifer picea, which represent a large increase over any previous gene mapping effort in gymnosperms e.g.. a total of major linkage groups were delineated in each species composite map, which is in agreement with cytogenetic studies indicating similar karyotypes and same numbers of chromosomes for both white spruce and black spruce and most other pinaceae. because of the agreement between number of linkage groups and number of chromosomes, the present maps could be considered as saturated. high genome coverage values also point to this observation. however, a direct pairing between linkage groups and chromosomes was not possible, as picea and most conifer chromosomes are difficult to differentiate based on standard cytological techniques. with the addition of several hundreds new gene markers, the observed total length of individual parental maps was increased by 15% to 25%, and the genome coverage values were increased by nearly 10% compared to previous maps mostly based on anonymous markers. with more than gene anchor markers positioned on each parental map, marker density was more than 50% higher and gene density was an order of magnitude higher than previously obtained. as compared to aflp and other types of dominant markers, these codominant snp gene markers mapped at such high density will contribute towards improving qtl mapping precision and power. more than gene loci were shared between the two composite maps developed herein, most of which being derived from the new gene snp developing effort. while previous efforts to increase the number of gene or ssr anchor markers relied on using several crosses per species e.g., the ease of mapping large numbers of transcripts using snps relaxes the need to implicate more than one cross for increasing the number of mapped anchor markers. however, using an additional cross may be highly useful to validate orthology of gene loci between species when breaks of synteny are observed. interspecific comparisons comparisons between the composite maps of white spruce and black spruce revealed high synteny and colinearity between their homoeologous linkage groups, in spite of the divergence of their lineages more than myr ago. while previous reports of genome comparisons between the two species reached similar conclusions about synteny and colinearity, the present observations are based on a more than four-fold increase in the number of mapped anchor loci in common between the two species. synteny between genomic regions can only be established if markers are true orthologs. the breakdown in synteny previously noted between white spruce and black spruce for lg iii and lg iv was confirmed in the present study. in addition, three new cases of putative inter-chromosomal translocations between both species were observed, but turned out to be false positives, after checks for locus orthology using haploid megagametophyte dna sequences or validating gene marker positions in a second cross of white spruce. thus, these three gene snps corresponded to paralogous loci between white spruce and black spruce. such false positives are likely to be frequent in conifer mapping studies, as previously observed. while synteny was well conserved, exceptions to colinearity between the two genomes were observed. small inversions between the two species composite maps involving a few closely spaced anchor loci were noted for all linkage groups, but no translocation within linkage groups was observed. most likely, these inversions resulted from analytical artifacts since the same discrepancies were observed between individual linkage maps within species. thus, it is safe to assume that these inversions resulted from the integration of both parental datasets for each species rather than from true chromosomal inversions. such a trend has also been observed when integrating individual linkage maps in other species.  the present report illustrates how new highly multiplexed snp genotyping approaches can be used to accelerate the structural analyses of complex and largely unknown genomes. the present maps are the most advanced genetic maps for spruce with regard to gene density and will open up several opportunities. these maps are currently used for qtl and eqtl detection and their next more densified versions will help target specific regions of the genome for future bac sequencing. the present maps will also be used to explore in more detail the distribution of multigene families over linkage groups and hence, the organization of the spruce genome and its evolution. such an approach, though at a small scale, has been previously used for the knox i multigene family, underlining a major chromosomal duplication in the spruce genome. maps with increasing gene densities are also required across the conifers to extend comparative genome studies and decipher the evolution of genome structure. such comparisons have already been undertaken based on a limited number of anchor markers e.g.. however, comparative genome studies in these species remain challenging since snps may have to be identified de novo in each species. this trend is best exemplified by the small number of gene snps shared between mapping populations belonging to different congeneric taxa: only snps were shared between white spruce and black spruce mapping populations out of candidate gene snps discovered in both species. using more crosses per species could contribute towards increasing the number of shared snps between species, but it is unlikely to change the order of magnitude of this number. on a more optimistic note, the primer transfer rate between white spruce and black spruce was high, which might facilitate resequencing efforts towards gene snp discovery in other spruce species. however, as species from different conifer genera are targeted, the primer transfer rate is expected to diminish. identifying sets of homologous transcripts in other pinaceae has improved in the recent years through efficient mining of sequence databases e.g. and because of large sets of est sequences for several species. nevertheless, when changes in chromosomal structure are observed between species based on gene linkage maps, gene orthology must be verified. distinguishing between orthology and paralogy can be precarious when the assessment is only based on homology scores, and without an adequate phylogenetic landscape for each gene family involved. in such conditions, a sensitive analysis at the sequence level is warranted before declaring such structural changes. thus, comparative mapping studies between conifer species call for further developments of gene maps and gene sequence collections, but also for genomics and bioinformatics tools enabling to assess more efficiently orthology relationships. authors' contributions np, jb: genotyping technology evaluation; bp, ig, ni: genetic mapping analyses; sbe, sbl, fg: resequencing, snp discovery, design of the snp arrays, and snp data analyses; np, ml, candidate gene identification; np, bp, jb: manuscript preparation; jb, ni: project design, funding, and overall supervision. all authors have read and approved the final manuscript. supplementary material additional file 1 data about the snps. this table provide references about the genes analysed: unigeneid in foresttreedb, annotations, marker names, gene names, links to annotation page in foresttreedb, primer sequences used to amplify sequences encompassing the unigenes. click here for file acknowledgements we thank alexandre montpetit and his team at the genome québec innovation centre for assistance with the acquisition of snp data on the illumina platform. this work was funded by grants from genome canada, genome québec, and the canadian biotechnology strategy through the arborea project lead by j. mackay and j. bousquet.
the white spruce composite map contained markers including gene loci. each array included snps, identified by resequencing genomic dna from parents of each mapping population. to explore the potential value of high-throughput genotyping assays in the analysis of large and complex genomes, we designed two highly multiplexed illumina bead arrays using the goldengate snp assay for gene mapping in white spruce and black spruce. for white spruce and black spruce, respectively, % and % of genotyped snps had valid goldengate assay scores and segregated in the mapping populations. snp data were integrated to pre-existing alfp, estp, and ssr markers to construct two individual linkage maps and a composite map for white spruce and black spruce genomes. considering lineage divergence at least myr ago between the two spruces, interspecific comparison of homoeologous linkage groups revealed remarkable synteny and marker colinearity. also, markers including gene loci were positioned on the black spruce composite map. the design of customized highly multiplexed illumina snp arrays appears as an efficient procedure to enhance the mapping of expressed genes and make linkage maps more informative and powerful in such species with poorly known genomes. this genotyping approach will open new avenues for co-localizing candidate genes and qtls, partial genome sequencing, and comparative mapping across conifers.
 to explore the potential value of high-throughput genotyping assays in the analysis of large and complex genomes, we designed two highly multiplexed illumina bead arrays using the goldengate snp assay for gene mapping in white spruce and black spruce.  each array included snps, identified by resequencing genomic dna from parents of each mapping population. for white spruce and black spruce, respectively, % and % of genotyped snps had valid goldengate assay scores and segregated in the mapping populations. for each of these successful snps, on average, valid genotyping scores were obtained for over 99% of progeny. snp data were integrated to pre-existing alfp, estp, and ssr markers to construct two individual linkage maps and a composite map for white spruce and black spruce genomes. the white spruce composite map contained markers including gene loci. also, markers including gene loci were positioned on the black spruce composite map. in total, anchor markers were shared between the two species. considering lineage divergence at least myr ago between the two spruces, interspecific comparison of homoeologous linkage groups revealed remarkable synteny and marker colinearity.  the design of customized highly multiplexed illumina snp arrays appears as an efficient procedure to enhance the mapping of expressed genes and make linkage maps more informative and powerful in such species with poorly known genomes. this genotyping approach will open new avenues for co-localizing candidate genes and qtls, partial genome sequencing, and comparative mapping across conifers. 