Skip to main content
  • Research article
  • Open access
  • Published:

Near-saturated and complete genetic linkage map of black spruce (Picea mariana)

Abstract

Background

Genetic maps provide an important genomic resource for understanding genome organization and evolution, comparative genomics, mapping genes and quantitative trait loci, and associating genomic segments with phenotypic traits. Spruce (Picea) genomics work is quite challenging, mainly because of extremely large size and highly repetitive nature of its genome, unsequenced and poorly understood genome, and the general lack of advanced-generation pedigrees. Our goal was to construct a high-density genetic linkage map of black spruce (Picea mariana, 2n = 24), which is a predominant, transcontinental species of the North American boreal and temperate forests, with high ecological and economic importance.

Results

We have developed a near-saturated and complete genetic linkage map of black spruce using a three-generation outbred pedigree and amplified fragment length polymorphism (AFLP), selectively amplified microsatellite polymorphic loci (SAMPL), expressed sequence tag polymorphism (ESTP), and microsatellite (mostly cDNA based) markers. Maternal, paternal, and consensus genetic linkage maps were constructed. The maternal, paternal, and consensus maps in our study consistently coalesced into 12 linkage groups, corresponding to the haploid chromosome number (1n = 1x = 12) of 12 in the genus Picea. The maternal map had 816 and the paternal map 743 markers distributed over 12 linkage groups each. The consensus map consisted of 1,111 markers distributed over 12 linkage groups, and covered almost the entire (> 97%) black spruce genome. The mapped markers included 809 AFLPs, 255 SAMPL, 42 microsatellites, and 5 ESTPs. Total estimated length of the genetic map was 1,770 cM, with an average of one marker every 1.6 cM. The maternal, paternal and consensus genetic maps aligned almost perfectly.

Conclusion

We have constructed the first high density to near-saturated genetic linkage map of black spruce, with greater than 97% genome coverage. Also, this is the first genetic map based on a three-generation outbred pedigree in the genus Picea. The genome length in P. mariana is likely to be about 1,800 cM. The genetic maps developed in our study can serve as a reference map for various genomics studies and applications in Picea a nd Pinaceae.

Background

Genetic maps provide an important genomic resource for understanding genome organization and evolution, comparative genomics, mapping genes and quantitative trait loci, and associating genes and genomic segments with phenotypic traits, especially in those species whose genomes are not yet completely sequenced. For understanding the genetic architecture of species, genetic maps with high levels of genome coverage and confidence in the marker order are required. High-density genetic maps and identification of genes or genetic factors controlling traits related to productivity, health, and adaptation to climatic change could accelerate forest tree improvement programs. Conifers are economically and ecologically important, and are the dominant tree species of the boreal and temperate forests. Genetic mapping and other genomics research is challenging in conifers, mainly because of their very large genome size (~25-30 Gbp) [1], the long time required to reach sexual maturity, inbreeding depression, and a general lack of advanced-generation pedigrees. The genome of any conifer species has yet to be completely sequenced.

Picea (spruce) is the second largest genus after Pinus (pine) in the family Pinaceae of conifers. Black spruce (Picea mariana (Mill.) B.S.P.) is a widespread transcontinental species of the North American boreal and temperate forests [2], and has great ecological and economic importance. It is one of the most important trees in Canada for the production of pulp and paper, and is one of the most important reforestation species in Canada [3]. Black spruce is a diploid species with haploid chromosome number (n) of 12 (2n = 2x = 24), like most other Pinaceae members. The estimated genome size of black spruce is large (1C = 15.8 pg; [4]; http://www.rbgkew.org.uk/cval/homepage.html) with an approximate 2C genome length of about 31,000 mpb http://www.rbgkew.org.uk/cval/homepage.html.

Although the first genetic linkage map in conifers was constructed for a single white spruce (Picea glauca (Moench) Voss) tree in 1992 from the analysis of haploid megagametophytes [5], the progress in the spruce genome mapping has been rather slow, particularly compared with the genus Pinus. Genetic linkage maps have been constructed for Norway spruce (Picea abies L.) [6–9], white spruce [5, 10, 11], and a black × red spruce (Picea rubens Sarg.) hybrid complex with an unknown proportion of the black spruce and red spruce genetic contribution to this hybrid [12, 13]. A parentage test with species-specific DNA markers revealed that the crosses used in [12] and [13] harbored a substantial amount of the red spruce genetic background. The markers used in the above-reviewed genome mapping studies were random amplified polymorphic DNA (RAPD) or a combination of RAPD, amplified fragment length polymorphism (AFLP), microsatellite/simple sequence repeat (SSR), expressed sequence tag polymorphism (ESTP), selectively amplified microsatellite polymorphic loci (SAMPL), single nucleotide polymorphism (SNPs), and/or 5 S rDNA. With the exception of the maps constructed for Norway spruce and white spruce from F1 mapping populations [9, 11] and the map constructed for black × red spruce hybrids from F1 and BC1 mapping populations [12, 13], all other maps were constructed for single trees from the segregation of a small number of markers in haploid megagametophytes. Single-tree genetic maps are of limited value. In predominantly outcrossing plants, such as conifers, a three-generation outbred pedigree (TGOP) is considered to be more informative than F1 or BC1 pedigree [14]. However, there is no published genetic linkage map in the genus Picea based on a TGOP. Moreover, most of the published spruce genetic linkage maps have not coalesced into 12 linkage groups corresponding to the haploid chromosome number of Picea. The linkage groups have ranged from 12 to 29. The first single tree genetic linkage map of white spruce developed from 47 RAPD markers coalesced into 12 linkage groups. Although the consensus map of Picea abies [9] and the composite map of P. mariana X P. rubens hybrids [12, 13] coalesced into 12 linkage groups, the maternal and/or paternal maps in these species coalesced into 13-23 linkage groups. There is no information published on a genetic linkage map in pure black spruce.

The objective of this study was to construct a high-density genetic linkage map of black spruce. We have used a three-generation outbred pedigree (TGOP) to develop a high density to near-saturated genetic map of black spruce. Here, we report maternal, paternal, and near-saturated consensus genetic linkage maps developed for black spruce using AFLP, SAMPL, ESTP, and SSR markers.

Results

AFLP markers

Forty AFLP primer combinations generated 809 markers segregating according to the expected Mendelian ratios. The number of polymorphic fragments ranged from 2 to 52, with an average of 20 polymorphic fragments per primer combination (Table 1). The average number of polymorphic fragments obtained per primer combination was 27, 18, 16 with the use of three, four, and five selective nucleotides (at the selective amplification step), respectively. The size of the segregating polymorphic fragments without adapters ranged from 40 to 660 bp, with only one mapped fragment below 50 bp and seven between 50 and 60 bp. Of the 809 markers, 485 segregated in ratio of 1:1 and 324 in ratio of 3:1. The number of markers segregating in the 1:1 ratio was 253 in the maternal parent and 232 in the paternal parent.

Table 1 AFLP primer combinations used, and the number and size of polymorphic fragments, and their segregation ratios.

SAMPL markers

A total of 255 SAMPL markers, segregating according to the expected Mendelian ratios, was obtained from the 12 SAMPL-Mse I primer combinations (Table 2; Table 3). The fragment size without adapters ranged from 32 to 661 bp (Table 3). Only six mapped SAMPL markers were of < 50 bp and 6 from 50 to 60 bp without adapters. The number of polymorphic fragments ranged from 2 to 44, with an average of 21 polymorphic fragments per SAMPL primer combination. Of the 255 SAMPL markers, 149 segregated in the ratio of 1:1, whereas 106 SAMPL markers segregated in the ratio of 3:1. The number of SAMPL markers segregating in the 1:1 ratio was 94 in the maternal parent and 54 in the paternal parent.

Table 2 SAMPL primers developed from Lactuca (from Witsenboer et al. [41]) compound microsatellite repeats and used for SAMPL marker mapping in the black spruce mapping population.
Table 3 SAMPL primer and Mse I primer extension combinations used, and the number and size of polymorphic fragments, and their segregation ratios.

Microsatellite markers

Twenty white spruce EST-based SSR loci were polymorphic between the parents of the mapping population and were mapped on the consensus map (Table 4). Six of the 20 polymorphic loci were heterozygous in the female parent as dominant makers and fourteen in both parents as co-dominant makers. Two alleles at a SSR locus heterozygous in the male or female parent, segregated in a 1:1 ratio in the progeny. Where both the parents were heterozygous, the progeny segregated either in a 1:2:1 or 1:1:1:1 ratio for their parental alleles (2-4). The known genes mapped using the microsatellites from white spruce EST sequences were as follows: Linkage Group (LG) I - Cytochrome B561 (RPGSE40); LG III - TIR/P-loop/LRR (RPGSE37); LG VI - Cytochrome P450 (RPGSE10) and putative UNC-50 (RPGSE46); LG VII - Early light inducible protein (RPGSE11) and Chloroplast nucleoid DNA binding protein (RPGSE34); LG VIII - Homeotic protein BEL1 (RPGSE02); LG IX - Auxin-induced protein 1 (RPGSE13); and LG XI- Putative Beta glycosidase (RPGSE03) and Gibberellin 12-oxidase (RPGSE04).

Table 4 Microsatellite DNA markers used for genetic mapping in black spruce (Picea mariana).

Six black spruce EST-based SSR loci that were polymorphic between the parents as dominant markers were mapped on the consensus map; four were heterozygous in the female and two in the male parent for a null allele (Table 4). These markers segregated in a 1:1 ratio in the progeny. Sixteen black spruce genomic sequence-based SSR loci were polymorphic between the parents and were mapped on the consensus map (Table 4). Eight microsatellite markers were heterozygous for a null allele in the female parent and five were heterozygous in the male parent. These markers segregated in a 1:1 ratio as dominant markers. Three microsatellite loci were heterozygous in both parents and the progeny segregated in a 1:2:1 or 1:1:1:1 ratio as co-dominant markers.

ESTP markers

Five ESTP markers showed Mendelian segregation in the progeny. Three ESTP loci were heterozygous in the female parent and two were heterozygous in the male parent. These markers segregated in a 1:1 ratio. Two ESTP markers were mapped on LG XI as well as one each on three linkage groups I, V, and IX (Figures 1, 2, 3).

Figure 1
figure 1

Linkage groups I to IV of the consensus genetic linkage map of black spruce ( Picea mariana ) constructed using 809 AFLPs (*A), 255 SAMPL (*S), 5 ESTPs (italicized and in red) and 42 SSR (italicized and in red color) markers. Names of the markers are provided on the right side of the linkage groups, with the DNA fragment size in bp. Genetic map distances, in cM, are provided on the left side of the linkage groups. AFLP and SAMPL markers in bold and italics are intercross markers.

Figure 2
figure 2

Linkage groups V to VIII of the consensus genetic linkage map of black spruce ( Picea mariana ) constructed using 809 AFLPs (*A), 255 SAMPL (*S), 5 ESTPs (italicized and in red) and 42 SSR (italicized and in red color) markers. Names of the markers are provided on the right side of the linkage groups, with the DNA fragment size in bp. Genetic map distances, in cM, are provided on the left side of the linkage groups. AFLP and SAMPL markers in bold and italics are intercross markers.

Figure 3
figure 3

Linkage groups IX to XII of the consensus genetic linkage map of black spruce ( Picea mariana ) constructed using 809 AFLPs (*A), 255 SAMPL (*S), 5 ESTPs (italicized and in red) and 42 SSR (italicized and in red color) markers. Names of the markers are provided on the right side of the linkage groups, with the DNA fragment size in bp. Genetic map distances, in cM, are provided on the left side of the linkage groups. AFLP and SAMPL markers in bold and italics are intercross markers.

Genetic linkage maps

The maternal map consisted of 816 markers distributed on 12 linkage groups covering 1,597 cM (Tables 5, 6). The number of mapped markers ranged from 35 to 93, with an average of 68 markers per linkage group. The length of the linkage groups ranged from 112 to 158 cM, with an average of 133 cM per linkage group (Table 6). The paternal map consisted of 743 markers assigned on 12 linkage groups, which covered 1,636 cM (Tables 5, 6). The number of mapped markers ranged from 25 to 89, with an average of 62 markers per linkage group. The length of the linkage groups ranged from 110 to 181 cM, with an average of 136 cM per linkage group (Table 6).

Table 5 Marker systems used for the construction of genetic linkage map and the number of markers mapped.
Table 6 Linkage groups, markers mapped, and marker density for the maternal, paternal, and consensus linkage maps in black spruce.

The homologous linkage groups between the parents were identified on the basis of segregating intercross AFLP, SAMPL, and SSR markers in the maternal and paternal maps. At least three intercross markers per linkage group were used. The integrated data set from the maternal and paternal maps allowed construction of a consensus linkage map. The consensus linkage map composed of 1,111 markers (Tables 5, 6) mapped to 12 linkage groups (Figures1, 2, 3). The linkage groups correspond to the haploid chromosome number (n = 12) of black spruce. It is worth noting that we have consistently obtained 12 linkage groups for the maternal, paternal, and consensus linkage maps, unlike other mapping studies where parental and/or consensus maps did not coalesce into 12 linkage groups. The consensus map covered 1,770 cM, with an average of 93 markers per linkage group and an average of one marker every 1.6 cM. The size of linkage groups varied from 131 to 172 cM, with an average of 148 cM (Table 6; Figures 1, 2, 3). The maternal, paternal and consensus genetic maps for each of the 12 linkage groups aligned very well and showed almost perfect colinearity of the marker order (Figure S1- Additional File 1).

Genome length and map coverage

The estimated length of the black spruce genome was 1,786 cM based on the method of Chakravarti et al. [15], and 1,819 cM according to the method of Fishmann et al. [16]. The observed length of the black spruce genome obtained from the consensus genetic linkage map was 1,770 cM. Thus, the consensus genetic map constructed in our study covered more than 97% of the estimated genome length of black spruce.

Distribution of markers along linkage groups

Significant deviations from the Poisson distribution of markers were observed for marker intervals of 2.5 cM, 5 cM, 10 cM, 20 cM, and 40 cM. For a 10 cM interval, the significant deviation (P < 0.001) is shown in Figure 4, indicating that the markers were not randomly distributed in the black spruce linkage groups. Marker distribution for other intervals (2.5 cM, 5 cM, 20 cM, and 40 cM) also showed clustering of markers (P < 0.05) along linkage groups. The independent analysis for testing the random distribution of AFLP (P < 0.001) and SAMPL (P < 0.001) markers indicated deviations from the random distribution. No correlation was observed between the number of mapped markers and the size of linkage groups. These results further support the clustering of makers on the linkage map.

Figure 4
figure 4

Poisson distribution function for the observed and expected frequencies of the markers distributed at 10 cM interval.

The distance between two adjacent markers on the linkage groups varied from 0 to 23.7 cM, with an average distance of 1.6 cM between any two adjacent markers (Figure 5; Table 6). This distance distribution reveals a strong skewness (P < 0.05), further indicating the non-random distribution of the markers along the linkage groups (Figure 5). Among the 1,099 intervals on 12 different linkage groups, 820 intervals were smaller than 2 cM (74.6%), and 89 intervals were larger than 6 cM (8.1%).

Figure 5
figure 5

Distribution of the map distance between two adjacent mapped markers.

Discussion

Genetic linkage map

We have developed a high density to near-saturated and complete genetic linkage map of black spruce. This is the first genetic map for black spruce, although genetic maps have recently been reported for a black spruce × red spruce hybrid complex [12, 13]. Except in the southern part of its range, red spruce is largely sympatric with black spruce. These two species hybridize in nature, although interspecific crossability represents a substantial but imperfect reproductive barrier for maintaining the separation of the species [17]. The differentiation of black and red spruce and their interspecific hybrids based on DNA markers, as used in Pelgas et al. [12] and Pavy et al. [13], is quite tenuous. The parents of the mapping pedigree in our study represent pure black spruce [18, 19].

Our black spruce genetic map is the first map in the genus Picea based on a three-generation outbred pedigree (TGOP). There are only three other published genetic maps in the genus Picea that are based on pedigreed material: one each in Norway spruce and white spruce based on F1 mapping populations [9, 11] and one in the black × red spruce hybrids based on F1 and BC1 mapping populations [12, 13]. Almost all other reported maps are for single trees (Table 7). The single-tree genetic maps are based on segregation of markers in haploid megagametophytes of maternal trees and do not take into account the segregation or recombination of markers in the paternal trees. Also, in conifers, the recombination rates were reported to be lower in female gametes than those in male gametes [20, 21]. Thus, single-tree linkage maps are not as informative as genetic maps developed from diploid segregating pedigreed populations, and are specific to a single tree (genotype), and thus, having limited or often no use for QTL mapping. Black spruce, like most other spruce or pine species, is highly outbred [22, 23]. For genome and QTL mapping in outbred plants, TGOP is more informative than any other pedigree used so far in the genus Picea [14]. Indeed, TGOP not only allows differentiation of up to four segregating alleles at a locus but also establishment of linkage phase among alleles in the mapping population [14]. This information is required to use genetic linkage maps for QTL mapping. Furthermore, the outbred pedigrees are more representative of natural populations in an outbred plant. In Pinus, where the genome and QTL mapping work has been more advanced than in its sister genus Picea, high-density genetic and QTL maps have been prepared using TGOP in loblolly (Pinus taeda L.) and maritime (Pinus pinaster Aiton) pine [24–28].

Table 7 Comparison of the genetic linkage maps constructed in black spruce (Picea mariana) with those constructed (published) for other species in the genus Picea. Mega = megagametophytes.

Genome length and map coverage

The black spruce genome length estimated in our study is 1,770 cM (Kosambi). This is comparable with the genome length of 1845.5 or 1849.8 cM (Kosambi) reported for black × red spruce hybrids based on a composite map of F1 and/or BC1 mapping populations (Table 7) [12, 13] and that of 1,865 cM estimated for black × red spruce controlled cross hybrids from a BC1 mapping population (Kang et al. unpublished data). These results are consistent with the length of the genetic map for black spruce that is likely to be between 1800 and 1900 cM. The length of the black spruce genetic map observed in our study is about 15% shorter than that reported for Norway spruce (2,035 cM) [9] and white spruce (2,007, 2,059, or 1933.5 cM) [10, 11] maps. The genome size (1C nuclear DNA contents) of Norway spruce (18.6 pg) [29] and white spruce (20.2 pg) [4] is 17.8% and 27.8%, respectively, higher than that of black spruce (15.8 pg) [4]. Although no direct relationships between the nuclear DNA contents and genetic map lengths apparently exist, the shorter genetic map length in black spruce than in Norway spruce or white spruce is apparently consistent with its comparatively smaller genome size. Our results suggest that the genome length observed in our study covers more than 97% of the estimated black spruce genome length. Thus, the black spruce genetic map reported here can very well be considered as almost complete. This is the highest map coverage so far for any Picea species. The estimated length and the extent of coverage of genetic maps in different species could vary owing to differences in the mapping populations used, variation in recombination rates of the parents, and the number and types of markers used in linkage map construction [14]. The mapping pedigree and markers used in our study are different from those used in Norway spruce and white spruce (Table 7).

Linkage groups and marker density

In our study, the paternal, maternal, and consensus linkage maps consistently coalesced into 12 linkage groups, corresponding to the haploid chromosome number (n = 12) in black spruce. By contrast, in all other studies reported on genome mapping in the genus Picea (Table 7), with one exception, maternal, paternal, and/or consensus map did not coalesce into 12 linkage groups (see Table 7). The consensus map reported here has 1,111 markers distributed over 12 linkage groups, which represents an almost complete coverage of the black spruce genome, as the number of linkage groups corresponds to the haploid chromosome number. The average distance observed among adjacent markers mapped for the genetic map of black spruce in our study (1.6 cM) is comparable with or lower than that reported for the composite map of black spruce × red spruce complex [12, 13], as well as lower than that reported for Norway spruce (2.6 cM) (Table 7). However, the marker density of the black spruce genetic map reported here is the highest for any genetic map based on a single cross in the genus Picea, as well for the maternal and paternal genetic maps (see Table 7 for comparisons).

Marker systems

The genetic map of black spruce was constructed using AFLP, SAMPL, SSR, and ESTP marker systems. The AFLP and SAMPL systems provided a sufficient number of anonymous polymorphic and informative markers to construct a high density to near-saturated genetic map, whereas the SSR and ESTP systems provided highly informative and co-dominant markers. Although SSR and ESTP markers, due primarily to their codominant nature, would be preferred for genome mapping, availability of limited numbers of these markers precludes constructing a high density to saturated genetic map in conifers using only these markers. We mapped 809 AFLP markers that were resolved by 40 primer combination, showing an average multiplex ratio of 20 markers per primer combination. This multiplex ratio is comparable to that observed in Norway spruce (14) [9] and loblolly pine (21) [30]. We mapped 255 SAMPL markers on 12 different linkage groups. The multiplex ratio for SAMPL markers (21 polymorphic mapped markers per primer combinations) was comparable to that observed for AFLP markers. The only other report where SAMPL markers have been used for genetic linkage mapping in conifers is for Norway spruce [8], where 20 SAMPL markers were mapped using two primer combinations. A large number of AFLP and SAMPL markers segregated in the 3:1 ratio, which suggests that the parental genomes are highly heterozygous. As AFLP and SAMPL markers were dominant, the 3:1 segregating (intercross) markers were useful for aligning the parental maps to construct the consensus map, which cannot be established directly. Also the intercross markers can help to identify additional linkage groups that were not represented in the parental maps [31].

SSR loci provided highly informative markers. Forty-two SSR loci were mapped onto 11 linkage groups; 17 of these were highly informative for integrating the maternal and paternal maps to construct the consensus map. The mapped cDNA-based SSR markers are excellent candidates for comparative and composite mapping because these markers are expected to show high intraspecific homology and high interspecific orthology [e.g., [11]]. Also these markers allowed mapping of 10 known genes on seven different linkage groups.

Most EST primer pairs resolved multilocus patterns, which is not surprising given the occurrence of multigene families in conifers [32]. Also, the rate of polymorphism observed in ESTPs was very low, which suggests that more powerful methods, such as detection of single nucleotide polymorphisms (SNPs) can be used to increase the resolution of polymorphism [13]. Nevertheless, the mapped ESTP loci along with cDNA-based SSR markers provide good candidates for comparative mapping in Picea, Pinaceae, or other conifers [12, 33, 34].

The genetic map presented here is the first-generation genetic map for black spruce, which provides a framework to map SNP and other markers in the second-generation genetic map. We are planning to map SNP markers from candidate genes and genes expressed differentially in response to climate change and SNPs mapped by Pavy et al. [13] onto this map in future. This is further discussed below under the Future Perspective section. The mapped SNP markers from expressed genes are quite useful for QTL and association mapping of relevant traits because allelic variation in the genes could be linked or associated with trait phenotypes. However, these markers represent a very small proportion of the spruce genome. Assuming 50,000 expressed genes of average 1 kb size and genome size of 30 Gbp in black spruce, the transcribed genome represents less than 0.2% of the whole genome. SNPs in about 500 candidate genes are normally mapped, further reducing the proportion of the genome sampled to less than 0.002%, although the mapped genes may have genome-wide distribution. The anonymous markers, such as AFLP and SAMPL markers generally provide a whole-genome scan; thus a genetic map using these markers may cover a large proportion of the genome, if not the whole genome. However, anonymous markers such as AFLP and SAMPL, have limitations in tagging genes controlling traits of interest via QTL and association mapping because variation in these markers may not represent functional genetic variation. Nevertheless, the black spruce genetic map reported here provides a valuable genomic resource in Picea and Pinaceae.

Clustering of markers

Even though only those markers segregating in Mendelian ratios and not those showing a distorted segregation were used for the linkage analysis, clustering of AFLP and SAMPL markers was detected in the linkage groups of black spruce. These results agree with the clustering of AFLP markers reported for genetic maps of Picea abies [35], Pinus taeda [30], and Pinus sylvestris L[36], but in contrast to random distribution of AFLP markers in the genetic maps reported for Norway spruce [9] and black × red spruce hybrids [12]. It should be noted that the study of Scotti et al. [35] was specifically performed to examine the distribution of marker classes in a genetic linkage map of Norway spruce.

The non-random distribution of markers may be caused by non-random and unequal crossing over and recombination along the chromosome length. The recombination is suppressed in the centromeric and heterochromatic pericentromeric regions [37], and the presence of heterochromatin in pericentromeric regions is a general feature of plant chromosomes. Assuming a random distribution of markers, low levels of meiotic recombination may well cause markers that are physically well separated, to cluster on a linkage map.

Future Perspectives

The genetic map reported here is suitable for constructing a composite map of two TGOP in black spruce and for the envisaged comparative mapping with red spruce and black spruce × red spruce hybrid. The mapped EST-based microsatellites and ESTPs will provide very useful markers for this work. We are also mapping additional EST-based SSR markers. Black spruce and red spruce are highly genetically related species but their evolutionary relationships are controversial. Comparative genome mapping may shed some light on comparative genome organization (orthology, synteny and order of the markers) and evolution of these species. The map could also be used for comparative mapping with other genera of Pinaceae.

We are genotyping and phenotyping the larger mapping populations of the TGOP used in this study and another TGOP (> 300) of black spruce for QTL mapping of traits related to growth and adaptation to climate change. We plan to use about 1,500 SNPs from candidate genes and genes expressed differentially in response to climate change conditions as well as SNPs mapped by Pavy et al. [13] for our QTL mapping work, using a high throughput SNP genotyping platform such as Illumina's Golden Gate Genotyping Assay http://www.illumina.com/technology/goldengate_genotyping_assay.ilmn. Thus, the technologies, markers and genetic map developed in the present study provide an invaluable genomic resource for basic and applied genomics studies.

Conclusions

We have constructed the first high density to near-saturated genetic linkage map of black spruce, with above 97% genome coverage. The maternal, paternal, and consensus maps in our study consistently coalesced into 12 linkage groups, corresponding to the haploid chromosome number (1n = 1x = 12) of 12 in the genus Picea. Also, this is the first genetic map based on a three-generation outbred pedigree in the genus Picea. The genome length in P. mariana is likely to be about 1,800 cM. The genetic maps developed in our study can serve as a reference map for various genomics studies and applications in Picea a nd Pinaceae. It will provide a foundation and a valuable resource for comparative mapping, constructing composite maps, and mapping quantitative trait loci and determining the genetic basis of complex quantitative traits of interest, such as growth and adaptation to climate change.

Methods

Mapping population

A three generation outbred pedigree (TGOP), including the grandparents, parents (F1), and F2 progeny, was used to construct the black spruce genetic linkage map. The grandparents of this pedigree were part of a 7 × 7 diallel F1 controlled-cross experiment, performed by Dr. E.K. Morgenstern in the early 1970s at the Petawawa National Forestry Institute (PNFI), in Chalk River, Ontario, Canada (46° N, 77° 30' W) [18]. The seven parental trees used for the diallel cross were from a pure black spruce plantation established at the Petawawa Research Forest (PRF), but the exact origin of the trees is unknown, other than that they were grown from seeds collected from the Lake Simcoe-Rideau region in Ontario [18]. The F1 seedlings from the full-sib families of this diallel were planted in genetic tests at three sites at PRF in 1973 [18]. The parents of the mapping pedigree were crossed from the F1 genetic tests in 1987 and 1988 by Dr. Tim Boyle at PNFI to produce F2 controlled crosses [19]. The F2 family 643 (32 × 40) was selected for genetic mapping purposes based on near-top and bottom ranking of its parents for growth and 13C discrimination rate [38] and availability of sufficient number of F2 seeds. Grandparents, parents, and 90 F2 individuals from this family were used as the mapping population. The F2 progeny were raised and grown at the Canadian Forest Service-Atlantic Forestry Centre, Fredericton, New Brunswick, Canada (45° 52' N, 66° 31' W).

DNA extraction

Genomic DNA was extracted from needle tissues of the female grandparent and megagametophtyes of the male grandparent, and needle tissues of the parents and their F2 progenies, using the Qiagen DNeasy Plant® Mini Kit, following the manufacturer's protocol (Qiagen Inc. Mississauga, ON, Canada). Needle tissues from the paternal grandparent were not available as the tree was harvested from the plantation, but its open-pollinated seeds were stored at the Atlantic Forestry Centre. Therefore, to genotype this grandparent, we used DNA extracted from pooled megagametophyte tissues from 20 to 30 seeds. In spruce, the genetic constitution of haploid megagametophytes is the same as that of female gametes of the mother tree. The quality and quantity of DNA preparations were determined by electrophoresing the DNA samples along with a standard of Lambda DNA on 0.8% agarose gels followed by staining with ethidium bromide.

Marker systems

Four different marker systems were used to genotype the grandparents, parents, and F2 progeny of the TGOP: AFLPs, SAMPL, microsatellites or SSRs, and ESTPs. We used four marker systems in order to achieve better genome coverage because different marker types target different genomic regions.

AFLP markers

Because of the extremely large genome size of spruce, a standard AFLP protocol, based on Eco RI-Mse I digestion, Eco RI, and Mse I primer extension by 1 nucleotide in the preamplification step and 3 base extension to the Eco RI and Mse I primers in the selective amplification step [39], produced complex AFLP fragment patterns. We developed methods for high throughput resolution of high-quality and clearly scorable AFLP markers for black spruce, using LI-COR 4200L® (LI-COR Biosciences, Lincoln, NE, USA) or Beckmann Coulter CEQ 8000 GENETIC ANALYSIS SYSTEM® (Beckmann Coulter, Fullerton, CA, USA), by evaluating a variety of conditions, including Eco RI and Mse I restriction digestion time of black spruce genomic DNA, and the number of selective nucleotides used in the preamplification and selective amplification steps. The primer combinations producing consistent, clear, and easily scorable polymorphic AFLP markers were identified and used for genotyping the mapping population.

The AFLP method was essentially as described in Vos et al. [39], with some modifications. Black spruce genomic DNA (500 ng) was digested with 2U each of Eco RI and Mse I (New England Biolab Inc. Ipswich, MA, USA) for 3 h at 37°C, followed by incubation at 70°C for 20 min. The digested DNA was ligated overnight with the Eco RI and Mse I adapters in a total volume of 20 μl at 25°C, followed by incubation at 70°C for 20 min. This restriction-ligation (RL) mixture was diluted 1:5 with autoclaved deionized distilled water before using it in the preamplification step.

A 3 μl aliquot of the RL mixture was preamplified using Eco RI (E) and Mse I (M) preamplification primers with an extension of one or two selective nucleotides at the 3' end. Eco RI preamp primers (+1/+2):

(+1) 5'- GAC TGC GTA CCA ATT CA - 3'

(+2) 5'- GAC TGC GTA CCA ATT CAC -3'

Mse I preamp primers (+1/+2):

(+1) 5'-GAT GAG TCC TGA GTA AC - 3'

(+2) 5'-GAT GAG TCC TGA GTA ACC - 3'

The PCR profile for the preamplification step consisted of 20 cycles each of denaturation at 94°C for 30 sec, annealing at 56°C for 1 min, and extension at 72°C for 1 min, followed by a final soak at 10°C using a PTC-200 thermal cycler (MJ Research, Reno, NV, USA). After the preamplification step, the reaction mixture was diluted 1:50 with sterile deionized distilled water. A total of 54 different Eco RI and Mse I primer pairs were tested with one or two selective nucleotides at the preamplification step and three to five selective nucleotides at the selective amplification step. From these, 40 AFLP primer combinations were selected for further use in mapping. Selective amplifications were performed using these primer combinations with various selective nucleotide extensions (E+3/M+3, E+3/M+4, E+3/M+5) (Table 1). The reaction mixture for the selective amplification consisted of 2 μl of diluted preamplified template DNA, 1 U Taq polymerase, 2.5 ng of Eco RI labeled primer (IRDye 700 label for LI-COR and D2 or D3 label for Beckmann Coulter CEQ 8000 GENETIC ANALYSIS SYSTEM), 12.5 ng Mse I primer, 10× PCR buffer (MBI Fermentas Inc, Burlington, ON, Canada), 1.5 mM MgCl2, 0.2 mM each of all four dNTPs (MBI Fermentas Inc, Burlington, ON, Canada), and BSA (1 μg/μl) (Sigma-Aldrich, Oakville, ON, Canada). PCR amplification profile consisted of 12 cycles each of denaturation at 94°C for 30 sec, annealing at 65°C for 30 sec (with lowering of 0.7°C per cycle) and extension at 72°C for 1 min, followed by 23 cycles each of denaturation at 94°C for 30 sec, annealing at 56°C for 60 sec and extension at 72°C for 1 min, followed by a final soak at 10°C.

Reaction products following selective amplification were resolved either on LI-COR 4200L or Beckmann Coulter CEQ 8000 GENETIC ANALYSIS SYSTEM. For LI-COR, selective amplification products were resolved on 6.5% denaturing Long Ranger polyacrylamide gels (LI-COR Biosciences, Lincoln, NE, USA). Approximately 0.5 μl of each sample (10 μl of PCR product and 15 ul of loading dye) was loaded on the gel. IRD-labeled molecular-weight markers were loaded in three lanes as a size-standard. Electrophoresis was carried out using 1× TBE running buffer, with run parameters of 1500 V, 35 mA, 70 W, signal channel 3, motor speed 3, 50°C plate temperature and 16-bit pixel depth for collection of TIFF image files. Polymorphic fragments were visually scored in the TIFF image files. Only those markers that were segregating in a Mendelian ratio (χ2- test, P < 0.05) were scored. For Beckmann Coulter CEQ 8000 GENETIC ANALYSIS SYSTEM, 2 μl of the selective amplification product was added to 27.5 μl of sample loading solution and 0.5 μl of CEQ DNA size standard-600 (Beckmann Coulter, Fullerton, CA, USA), followed by overlaying a drop of mineral oil. Samples were injected into a 33 cm capillary at 2.0 KV for 90 sec and electrophoresed at 7.5 KV for 70 min at 35°C. The AFLP fragments data were exported to an Excel® file using fragment analysis software for further analysis of genetic linkage parameters.

SAMPL markers

Selectively amplified microsatellite polymorphic loci markers, based on a combination of AFLP and microsatellite technology, can combine good features of both AFLP and microsatellite markers, and can reduce the marker complexity of AFLPs in spruce. The SAMPL technology is a modified AFLP technique, in which a compound microsatellite sequence is used as one of the two AFLP primers in selective amplification, generally in place of Eco RI primers [40]. We developed SAMPL markers using the compound microsatellite repeats from Lactuca species [41] as SAMPL primers (Table 2) in place of the Eco RI primer in the selective amplification step. The SAMPL markers were analyzed on the LI-COR and Beckman CEQ 8000 systems, using the protocol described above for AFLP analysis as well as in Gupta et al. [40].

Sixteen combinations of four SAMPL and four Mse I primers (with an extension of three selective nucleotides) were tested to screen SAMPL marker polymorphisms between the parents of the mapping population. Of these, 12 primer combinations were selected for genotyping of the mapping population based on the quality and polymorphism of the markers resolved (Table 3). The SAMPL marker data were scored as described above for AFLP markers.

Microsatellite/SSR markers

Seventy-eight microsatellites developed from black spruce cDNA (EST) or genomic DNA sequences and white spruce ESTs in our lab were used to screen polymorphisms between the parents of the mapping population. Forty-two microsatellites showed inter-parental polymorphisms, and were used to genotype the mapping population (Table 4). Out of these, 20 were from the white spruce ESTs, 6 from black spruce ESTs, and 16 from black spruce genomic sequences (SSR-enriched and AFLP-SSR libraries). White spruce EST sequences were obtained from the publicly available NCBI GenBank EST database. EST sequences containing microsatellites were identified. Primers were designed and microsatellite markers were optimized. The mapped six black spruce EST-SSRs were developed from the EST sequences obtained from a cambium-transcript enriched cDNA library constructed from the male parent (40) of the mapping population. The details of the study on cDNA library construction and EST sequencing, analysis and annotation will be published elsewhere. Here, we provide the primer sequences and annealing temperatures for six mapped EST-SSR markers in the Additional File 2. The results on the development and characterization of microsatellite DNA markers from the white spruce EST sequences (RPGSE series) and from the black spruce genomic sequences (RPMSA and RPMSG series) will be published elsewhere because this work was performed by other researchers in the Rajora lab. However, pending publication, information on these mapped markers, including primer sequences, is available from the Principal Investigator O.P. Rajora of the Spruce Genomics Program. The microsatellite markers were resolved on the LI-COR system and data scored as described in Rajora et al. [42, 43].

ESTP markers

Primers for 198 ESTs, obtained from sequencing of a black spruce standard cDNA library prepared from needle tissue (to be published separately), were designed using Primer 3.0 software [44]. The parents of the mapping population were screened for ESTPs (length polymorphism). Most EST primers resolved multilocus patterns. Only five ESTP markers, showing inter-parental polymorphism, displayed unambiguous single-locus patterns and therefore, were used to genotype the mapping population. The primer sequences and annealing temperatures for these mapped ESTP markers are provided in the Additional File 2. The PCR amplification profile consisted of initial denaturation at 94°C for 5 min, 40 cycles each of denaturation at 94°C for 1 min, annealing at 55°C for 1 min, and extension at 72°C for 1.3 min, followed by final extension at 72°C for 10 min. The ESTP markers were resolved by electrophoresis on either 2% agarose or 6% polyacrylamide gels.

Nomenclature or labeling of markers on the genetic maps

The AFLP and SAMPL markers were named, starting with letters A, and S, respectively, followed by the primer number, and then the size of the fragment. The AFLP and SAMPL fragment sizes reported in this manuscript are without adapter sequences. The SSR markers were named with a prefix of five letters. The first letter represents the Principal Investigator/Lab (R = Rajora), the next two letters the species name (PG = Picea glauca, PM = Picea mariana) from which the markers were developed, the next letter S representing SSR, and the last letter denoting the source of sequences or library type (E = EST; G = Genomic; A = AFLP-SSR genomic). These prefix letters were followed by the marker number. Thus, SSR markers developed from the white spruce EST database have a prefix of RPGSE, SSR markers developed from black spruce EST sequences a prefix of RPMSE, SSR markers developed from the genomic library a prefix of RPMSG, and SSR markers developed from the black spruce SSR-enriched AFLP sequences a prefix of RPMSA. The ESTP markers developed from the EST sequences from the black spruce cDNA library were named starting with RPMEP, followed by the marker number.

Statistical analysis

Segregation analysis and map construction

Individual paternal and maternal maps were constructed according to two-way pseudo-testcross mapping strategy [45]. All linkage analysis and genetic map construction, including marker order and map length estimations, were performed using JOINMAP® 3.0 software [46] with the maximal threshold jump value of 5 and ripple value of 1. The Kosambi [47] mapping function was used for map length estimations. Our marker data set for genetic linkage mapping included three different segregation patterns: 1:1 for markers heterozygous in one parent and homozygous or null in the other, 3:1 for dominant markers heterozygous in both parents, and 1:2:1 or 1:1:1:1 for co-dominant markers heterozygous in both parents. The JOINMAP command "similarity of loci" was used to identify the similar loci. Only one of the markers was kept from the similar loci for linkage mapping analysis.

For framework, preliminary linkage grouping of AFLPs, SAMPL, SSRs, and ESTPs markers were ordered using the "Group" command, with LOD (log of odds) threshold maximum 5.0, minimum 4.0, recombination ratio 0.35. The map orders were found by calculating pairwise recombination frequencies, and map distances were estimated by a least-squares procedure. The two parental maps based on segregating markers were grouped and ordered using a minimum LOD score of 3.0 and recombination fraction of 0.4 as the grouping criterion. The marker order obtained from the third round of analysis was retained with the JOINMAP command "calculate map". This order was fixed to allow positioning of additional markers. To construct consensus maps, the maternal maps and paternal were aligned based on the co-dominant markers (1:2:1, 1:1:1:1 and 1:1) and intercross markers (3:1), and then consensus maps were constructed by using JoinMap function "Combine Groups for Map Integration" command. The maternal, paternal and consensus genetic maps were aligned using the JOINMAP. Graphic presentation of the individual linkage groups was drawn using Mapchart® version 2.0 software [48].

Estimation of genome length and map coverage

The length of the black spruce genome (G) was estimated using the Method 4 of Chakravarti et al. [15] after each length had been adjusted by the factor m+1/m-1, where m is the number of framework markers on the linkage group, as well by the method described by Fishmann et al. [16] that twice the value of map density (d) was added to the length of each linkage group to account for chromosome ends beyond the terminal markers, and G was calculated by summing up the resulting lengths of 12 linkage groups. The observed genome length was obtained by summing up the map lengths of the 12 individual linkage groups. The map coverage was calculated as the ratio of the observed to the estimated genome length. The number of markers required to cover the whole genome of black spruce was calculated according to Lange and Boehnke [49].

Marker distribution analysis

To evaluate whether the mapped markers were randomly distributed on the linkage map, the linkage groups were divided into 2.5, 5, 10, 20, and 40 cM blocks, and the number of markers per block was counted. Observed frequencies of the number of markers per block were compared with the expected ones by performing a Chi-square test [30, 36, 50], using a Poisson distribution function, P(x) = e-μμx/x!, where x is the number of markers per block and μ is the average marker density in the consensus map. Average marker density (μ) was used to calculate the expected binomial frequencies for each marker class per block interval for all the linkage groups. The distribution of markers on the linkage groups was also evaluated separately for the AFLP and SAMPL markers. The SSR and ESTP markers could not be considered independently for this analysis because of their small numbers or low frequencies on each linkage group.

References

  1. Neale DB, Williams CG: Restriction fragment length polymorphism mapping in conifers and applications to forest genetics and tree improvement. Can J For Res. 1991, 21: 545-554. 10.1139/x91-076.

    Article  CAS  Google Scholar 

  2. Viereck LA, Johnston WF: Picea mariana (Mill.) B.S.P. - Black spruce. Conifers. Silvics of North America. Edited by: Burns RM, Honkala BH. 1990, US Dept Agric, Forest Service, Agriculture Handbook 654, Washington DC, 1: 227-237.

    Google Scholar 

  3. Morgenstern EK, Wang BSP: Trends in forest depletion, seed supply, and reforestation in Canada during the past four decades. For Chron. 2001, 6: 1014-1021.

    Article  Google Scholar 

  4. Ohri D, Khoshoo TN: Genome size in gymnosperms. Plant Sys Evol. 1986, 153: 119-132. 10.1007/BF00989421.

    Article  Google Scholar 

  5. Tulsieram LK, Glaubitz JC, Kiss G, Carlson JE: Single tree genetic linkage analysis in conifers using haploid DNA from megagametophytes. BioTechnology. 1992, 10: 686-690. 10.1038/nbt0692-686.

    Article  CAS  PubMed  Google Scholar 

  6. Binelli G, Bucci G: A genetic linkage map of Picea abies Karst., based on RAPD markers, as a tool in population genetics. Theor Appl Genet. 1994, 88: 283-288. 10.1007/BF00223633.

    CAS  PubMed  Google Scholar 

  7. Bucci G, Kubisiak TL, Nance WL, Menozzi P: A population consensus partial linkage map of Picea abies Karst. based on RAPD markers. Theor Appl Genet. 1997, 95: 643-654. 10.1007/s001220050607.

    Article  CAS  Google Scholar 

  8. Paglia GP, Olivieri AM, Morgante M: Towards second generation STS (sequence-tagged sites) linkage maps in conifers: a genetic map of Norway spruce (Picea abies K.). Mol Gen Genet. 1998, 258: 466-478. 10.1007/s004380050757.

    Article  CAS  PubMed  Google Scholar 

  9. Acheré V, Faivre-Rampant P, Jeandroz S, Besnard G, Markussen T, Aragones A, Fladung M, Ritter E, Favre JM: A full saturated linkage map of Picea abies including AFLP, SSR, ESTP, 5 S rDNA and morphological markers. Theor Appl Genet. 2004, 108: 1602-1613. 10.1007/s00122-004-1595-y.

    Article  PubMed  Google Scholar 

  10. Gosselin I, Zhou Y, Bousquet J, Isabel N: Megagametophyte-derived linkage maps of white spruce (Picea glauca) based on RAPD, SCAR and ESTP markers. Theor Appl Genet. 2002, 104: 987-997. 10.1007/s00122-001-0823-y.

    Article  CAS  PubMed  Google Scholar 

  11. Pelgas B, Beauseigle S, Acheré V, Jeandroz S, Bousquet J, Isabel N: Comparative genome mapping among Picea glauca, P. mariana x P. rubens and P. abies, and correspondence with other Pinaceae. Theor Appl Genet. 2006, 113: 1371-1393. 10.1007/s00122-006-0354-7.

    Article  CAS  PubMed  Google Scholar 

  12. Pelgas B, Bousquet J, Beauseigle S, Isabel N: A composite linkage map from two crosses for the species complex Picea mariana X Picea rubens and analysis of synteny with other Pinaceae. Theor Appl Genet. 2005, 111: 1466-1488. 10.1007/s00122-005-0068-2.

    Article  CAS  PubMed  Google Scholar 

  13. Pavy N, Pelgas B, Beauseigle S, Blais S, Gagnon F, Gosselin I, Lamothe M, Isabel N, Bousquet J: Enhancing genetic mapping of complex genomes through the design of highly-multiplexed SNP arrays: application to the large and unsequenced genomes of white spruce and black spruce. BMC Genomics. 2008, 9: 21-10.1186/1471-2164-9-21.

    Article  PubMed Central  PubMed  Google Scholar 

  14. Liu BH: Statistical genomics: linkage mapping and QTL analysis. 1998, CRC Press, New York, NY

    Google Scholar 

  15. Chakravarti A, Lasher LK, Reefer JE: A maximum-likelihood method for estimating genome length using genetic linkage data. Genetics. 1991, 128: 175-182.

    CAS  PubMed Central  PubMed  Google Scholar 

  16. Fishman L, Kelly AJ, Morgan E, Willis JH: A genetic map in the Mimulus guttatus species complex reveals transmission ratio distortion due to heterospecific interactions. Genetics. 2001, 159: 1701-1716.

    CAS  PubMed Central  PubMed  Google Scholar 

  17. Major JE, Mosseler A, Johnsen K, Rajora OP, Barsi DC, Kim K-H, Park J-M, Campbell M: Reproductive barriers and hybridity in two spruces, Picea rubens and Picea mariana, sympatric in eastern North America. Can J Bot. 2005, 83: 163-175. 10.1139/b04-161.

    Article  Google Scholar 

  18. Morgenstern EK: A diallel cross in black spruce Picea mariana (Mill.) B.S.P. Silvae Genet. 1974, 23: 67-70.

    Google Scholar 

  19. Boyle TJB: A diallel cross in black spruce. Genome. 1987, 29: 180-186.

    Article  Google Scholar 

  20. Groover A, Williams CG, Devey ME, Lee JM, Neale DB: Sex-related differences in meiotic recombination frequency in Pinus taeda. J Hered. 1995, 86: 157-158.

    Google Scholar 

  21. Plomion C, O'Malley DM: Recombination rate differences for pollen parents and seed parents in pine. Heredity. 1996, 77: 341-350. 10.1038/hdy.1996.152.

    Article  CAS  Google Scholar 

  22. Boyle TJB, Morgenstern EK: Estimates of outcrossing rates in six populations of black spruce in central New Brunswick. Silvae Genet. 1986, 35: 102-106.

    Google Scholar 

  23. Sproulè AT, Dancik BP: The mating system of black spruce in north-central Alberta, Canada. Silvae Genet. 1996, 45: 159-164.

    Google Scholar 

  24. Devey ME, Fiddler TA, Liu BH, Knapp SJ, Neale BD: An RFLP linkage map for loblolly pine based on a three-generation outbred pedigree. Theor Appl Genet. 1994, 88: 273-278. 10.1007/BF00223631.

    CAS  PubMed  Google Scholar 

  25. Groover A, Devey M, Fiddler T, Lee J, Megraw R, Sherman B, Vujcic S, Williams C, Neale D: Identification of quantitative trait loci iflencing wood specific gravity in an outbred pedigree of loblolly pine. Genetics. 1994, 138: 1293-1300.

    CAS  PubMed Central  PubMed  Google Scholar 

  26. Sewell MM, Sherman BK, Neale DB: A consensus map for loblolly pine (Pinus taeda L.). I. Construction and integration of individual linkage maps from two outbred three-generation pedigrees. Genetics. 1999, 151: 321-330.

    CAS  PubMed Central  PubMed  Google Scholar 

  27. Sewell MM, Davis MF, Tuskan GA, Wheeler NC, Elam CC, Bassoni DI, Neale DB: Identification of QTLs influencing wood property traits in loblolly pine (Pinus taeda L.). II. Chemical wood properties. Theor Appl Genet. 2002, 104: 214-222. 10.1007/s001220100697.

    Article  CAS  PubMed  Google Scholar 

  28. Chagne D, Brown GR, Lalanne C, Madur D, Pot D, Neale DB, Plomion C: Comparative genome and QTL mapping between maritime and loblolly pines. Mol Breed. 2003, 12: 185-195. 10.1023/A:1026318327911.

    Article  CAS  Google Scholar 

  29. Siljak-Yakovlev S, Cerbah M, Coulaud J, Stoian V, Brown SC, Zoldos V, Jelenic S, Papes D: Nuclear DNA content, base composition, heterochromatin and rDNA in Picea omorika and Picea abies. Theor Appl Genet. 2002, 104: 505-512. 10.1007/s001220100755.

    Article  CAS  PubMed  Google Scholar 

  30. Remington DL, Whetten RW, Liu BH, O'Malley DM: Construction of an AFLP genetic map with nearly complete genome coverage in Pinus taeda. Theor Appl Genet. 1999, 98: 1279-1292. 10.1007/s001220051194.

    Article  CAS  PubMed  Google Scholar 

  31. Hanley S, Barker A, Van Ooijen W, Aldam C, Harris L, Ahman I, Larsson S, Karp A: A genetic linkage map of willow (Salix viminalis) based on AFLP and microsatellite markers. Theor Appl Genet. 2002, 105: 1087-1096. 10.1007/s00122-002-0979-0.

    Article  CAS  PubMed  Google Scholar 

  32. Kinlaw CS, Neale DB: Complex gene families in pine genomes. Trends Plant Sci. 1997, 2: 356-359. 10.1016/S1360-1385(97)84624-9.

    Article  Google Scholar 

  33. Komulainen P, Brown GR, Mikkonen M, Karhu A, Garcia-Gil MR, O'Malley D, Lee B, Neale DB, Savolainen O: Comparing EST-based genetic maps between Pinus sylvestris and P. taeda. Theor Appl Genet. 2003, 107: 667-678. 10.1007/s00122-003-1312-2.

    Article  CAS  PubMed  Google Scholar 

  34. Krutovsky KV, Troggio M, Brown GR, Jermstad KD, Neale DB: Comparative mapping in the Pinaceae. Genetics. 2004, 168: 447-461. 10.1534/genetics.104.028381.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  35. Scotti I, Burelli A, Cattonaro F, Chagné D, Fuller J, Hedley PE, Jansson G, Lalanne C, Madur D, Neale D, Plomion C, Powell W, Troggio M, Morgante M: Analysis of the distribution of marker classes in a genetic linkage map: a case study in Norway spruce (Picea abies Karst). Tree Genet Genom. 2005, 1: 93-102. 10.1007/s11295-005-0012-2.

    Article  Google Scholar 

  36. Yin TM, Wang XR, Andersson B, Lerceteau-Kohler E: Nearly complete genetic maps of Pinus sylvestris L. (Scots pine) constructed by AFLP marker analysis in a full-sib family. Theor Appl Genet. 2003, 106: 1075-1083.

    CAS  PubMed  Google Scholar 

  37. Tanksley SD, Ganal MW, Prince JP, de Vicente MC, Bonierbale MW, Broun P, Fulton TM, Giovannoni JJ, Grandillo S, Martin GB, Messeguer R, Miller JC, Miller L, Paterson AH, Pineda O, Riider MS, Wing RA, Wu W, Young ND: High-density molecular linkage maps of the tomato and potato genomes. Genetics. 1992, 132: 1141-1160.

    CAS  PubMed Central  PubMed  Google Scholar 

  38. Johnsen KH, Flanagan LB, Huber DA, Major JE: Genetic variation in growth and carbon isotope discrimination in Picea mariana: analyses from a half-diallel mating design using field grown trees. Can J For Res. 1999, 29: 1727-1735. 10.1139/cjfr-29-11-1727.

    Article  Google Scholar 

  39. Vos P, Hogers R, Bleeker M, Reijans M, Van de Lee T, Hornes M, Frijters A, Pot J, Peleman J, Kuiper M, Zabeau M: AFLP: a new technique for DNA fingerprints. Nucl Acids Res. 1995, 23: 4407-4414. 10.1093/nar/23.21.4407.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  40. Gupta AK, Kang BY, Roy JK, Rajora OP: Large-scale development of selectively amplified polymorphic microsatellite loci (SAMPL) markers in spruce (Picea). Mol Eco Not. 2005, 5: 481-483. 10.1111/j.1471-8286.2005.00963.x.

    Article  CAS  Google Scholar 

  41. Witsenboer H, Vogel J, Michelmore RW: Identification, genetic localization, and allelic diversity of selectively amplified microsatellite polymorphic loci in lettuce and wild relatives (Lactuca spp.). Genome. 1997, 40: 923-936. 10.1139/g97-119.

    Article  CAS  PubMed  Google Scholar 

  42. Rajora OP, Rahman MH, Dayanandan S, Mosseler A: Isolation, characterization, inheritance and linkage of microsatellite DNA markers in white spruce (Picea glauca) and their usefulness in other spruce species. Mol Gen Genet. 2001, 264: 871-882. 10.1007/s004380000377.

    Article  CAS  PubMed  Google Scholar 

  43. Rajora OP, Mann IK, Shi Y-Z: Genetic diversity and population structure of boreal white spruce (Picea glauca) in pristine conifer-dominated and mixedwood forest stands. Can J Bot. 2005, 83: 1096-1105. 10.1139/b05-083.

    Article  CAS  Google Scholar 

  44. Rozen S, Skaletsky HJ: Primer 3 on the WWW for general users and for biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz S, Misener S. 2000, Humana Press, Totowa, NJ, USA, 365-386.

    Google Scholar 

  45. Grattapaglia D, Sederoff R: Genetic linkage maps of Eucalyptus grandis and Eucalyptus urophylla using a pseudo testcross mapping strategy and RAPD markers. Genetics. 1994, 137: 1121-1137.

    CAS  PubMed Central  PubMed  Google Scholar 

  46. Van Ooijen JW, Voorrips RE: JOINMAP 3.0, Software for the calculation of genetic linkage maps. 2001, Plant Research International, Wageningen, the Netherlands, [http://www.plant.wageningen-ur.nl]

    Google Scholar 

  47. Kosambi DD: The estimation of map distances from recombination values. Ann Eugen. 1944, 12: 172-175.

    Article  Google Scholar 

  48. Voorrips RE: MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered. 2002, 93: 77-78. 10.1093/jhered/93.1.77.

    Article  CAS  PubMed  Google Scholar 

  49. Lange K, Boehnke M: How many polymorphic genes will it take to span the human genome?. Am J Hum Genet. 1982, 34: 842-845.

    CAS  PubMed Central  PubMed  Google Scholar 

  50. Cervera MT, Storme V, Ivens B, Gusmao J, Liu BH, Hostyn V, Slycken JV, Montagu MV, Boerjan W: Dense genetic linkage maps of three Populus species (Populus deltoides, P. nigra and P. trichocarpa) based on AFLP and microsatellite markers. Genetics. 2001, 158: 787-809.

    CAS  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgements

We thank (late) Dr. Kris Morgenstern for making F1 crosses; Dr. Tim Boyle for making F2 controlled crosses; Debby Barsi, Moira Campbell, and Stephanie West for raising, cloning, and maintaining the mapping population and plant material collection; Joy Roy, Nicholas Buckley and Wendy Wiesser for their laboratory assistance in microsatellite genotyping. I.K. Mann was financially supported by a Dalhousie University Graduate Scholarship and through an NSERC (Natural Sciences and Engineering Research Council of Canada) Strategic Project grant (STPGP 234783-00) to O.P. Rajora. B.Y. Kang was supported from the Genome Canada and NSERC (STPGP 234783-00) funds to O.P. Rajora. O.P. Rajora held the Stora Enso Senior Chair in Forest Genetics and Biotechnology at Dalhousie University, which was supported by Stora Enso Port Hawkesbury Ltd., and currently holds the Senior Canada Research Chair in Forest and Conservation Genomics and Biotechnology at UNB, which is supported by the Canada Research Chair Program (CRC950-201869). The research was funded by an NSERC Strategic Project grant (STPGP 234783-00) and the Canada Research Chair Program (CRC950-201869) to O.P. Rajora and partially by the Canadian Forest Service operating funds to J.E. Major. We thank two anonymous reviewers for their thorough review of the manuscript and providing useful comments and suggestions.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Om P Rajora.

Additional information

Authors' contributions

BYK: SAMPL genotyping, data analysis and construction and alignment of the genetic maps; IKM: AFLP and ESTP genotyping, early data analysis and genetic map construction, preparation of preliminary drafts of the parts of the manuscript; JEM: mapping population and its selection, propagation and sampling, and manuscript revision; OPR: project conception, design, funding, guidance and overall supervision, manuscript preparation and revision. All authors have read and approved the manuscript.

Bum-Yong Kang, Ishminder K Mann contributed equally to this work.

Electronic supplementary material

12864_2010_3109_MOESM1_ESM.PDF

Additional file 1: Figure S1 Alignment of the maternal (Maternal), paternal (Paternal) and consensus (PM643) linkage maps of black spruce. Names of the markers are provided on the right side of the linkage groups, with the DNA fragment size in bp. Genetic map distances, in cM, are provided on the left side of the linkage groups. AFLP markers start with A, SAMPL markers with S. The SSR and ESTP markers are italicized and are in red color. (PDF 127 KB)

12864_2010_3109_MOESM2_ESM.DOC

Additional file 2: Table S1 Expressed sequence tag polymorphism (ESTP) and EST-based microsatellite (SSR) markers developed from black spruce ESTs and mapped on the black spruce genetic map in this study. (DOC 37 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Kang, BY., Mann, I.K., Major, J.E. et al. Near-saturated and complete genetic linkage map of black spruce (Picea mariana). BMC Genomics 11, 515 (2010). https://doi.org/10.1186/1471-2164-11-515

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-11-515

Keywords