Species-specific RAPD fingerprints for the closely related Picea mariana and P. rubens

Species-specific molecular markers were designed to assist in the identification of closely related black spruce (Picea mariana [B.S.P.] Mill.) and red spruce (P. rubens Sarg.) in northeastern North America. Trees from six provenances of black spruce and three provenances of red spruce were sampled from outside the sympatric zone. They were first classified using a composite index of five qualitative morphological traits. The species-specific genetic markers were developed using random amplified polymorphic DNAs (RAPD) and a combination of bulk sample and individual tree analyses. Each species bulk sample was constructed from DNAs obtained from 12 trees that were from outside the sympatric zone and showed a morphological composite index specific of each species. A total of 161 primers were screened with the bulk samples. From these, 52 primers showing segregating fingerprints were further screened with the individual trees. Most of the markers observed were shared by the two species, and there was less diversity in P. rubens. A small number of markers were found to be monomorphic or nearly monomorphic and specific to either P. mariana or P. rubens. These markers remained species-specific when F₁ progenies derived from independent intraspecific crosses were screened, and they were subsequently found to co-segregate in hybrids derived from independent interspecific crosses here used as controls.     

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

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

Genetic validation and characterization of RAPD markers differentiating black and red spruces: molecular certification of spruce trees and hybrids

 Picea mariana (black spruce) and P. rubens (red spruce) are closely related species which are difficult to differentiate morphologically. RAPD markers differentiating black and red spruces have been previously identified. In the present study, genetic validity of these markers was determined using samples representing range–wide provenances. Their applicability for certifying genetic identity of individual black, red trees and their hybrids from several sympatric and allopatric locations was demonstrated. These diagnostic fragments of both red and black spruce were present at a frequency of over 0.95 in allopatric provenances, but at a lower frequency in some sympatric provenances (0.43–1.00). Natural populations of red spruce exhibiting typical red spruce phenotype contained black spruce diagnostic RAPD fragments and black spruces growing in bogs with typical bog black spruce morphology, contained red spruce-specific RAPD markers. Some major RAPD markers were cloned and sequenced. The results reveal an extremely high degree of identity between the random primer and the primer binding sites on the genome. Amplification of black and red spruce genomic DNA with designed primers flanking the species-diagnostic RAPD markers indicates that most of RAPD markers used to differentiate black spruce from red spruce are not species specific since these sequences were detected in several spruce species using a more sensitive detection method. Received June 17, 2002; accepted August 5, 2002 Published online: February 4, 2003     

Tracking the progression of speciation: variable patterns of introgression across the genome provide insights on the species delimitation between progenitor–derivative spruces (Picea mariana × P. rubens)

The genic species concept implies that while most of the genome can be exchanged somewhat freely between species through introgression, some genomic regions remain impermeable to interspecific gene flow. Hence, interspecific differences can be maintained despite ongoing gene exchange within contact zones. This study assessed the heterogeneous patterns of introgression at gene loci across the hybrid zone of an incipient progenitor–derivative species pair, Picea mariana (black spruce) and Picea rubens (red spruce). The spruce taxa likely diverged in geographic isolation during the Pleistocene and came into secondary contact during late Holocene. A total of 300 SNPs distributed across the 12 linkage groups (LG) of black spruce were genotyped for 385 individual trees from 33 populations distributed across the allopatric zone of each species and within the zone of sympatry. An integrative framework combining three population genomic approaches was used to scan the genomes, revealing heterogeneous patterns of introgression. A total of 23 SNPs scattered over 10 LG were considered impermeable to introgression and putatively under diverging selection. These loci revealed the existence of impermeable genomic regions forming the species boundary and are thus indicative of ongoing speciation between these two genetic lineages. Another 238 SNPs reflected selectively neutral diffusion across the porous species barrier. Finally, 39 highly permeable SNPs suggested ancestral polymorphism along with balancing selection. The heterogeneous patterns of introgression across the genome indicated that the speciation process between black spruce and red spruce is young and incomplete, albeit some interspecific differences are maintained, allowing ongoing species divergence even in sympatry. The approach developed in this study can be used to track the progression of ongoing speciation processes.     

Identification of mitochondrial plasmid-like DNAs in Picea abies (L.) Karst.

The Norway spruce (Picea abies (L.) Karst.) mitochondrial DNA has been extracted from embryonal suspensor masses. In addition to a master chromosome, a family of plasmid-like DNAs were identified. These latter shared cross homologies but had no evident sequence homology with the master chromosome. The occurrence of mitochondrial plasmid-like DNAs was investigated in trees from different provenances. A vast majority of trees displayed extrachromosomal DNA elements of variable stoechiometry. For some trees, the sequences homologous to the extrachromosomal DNA elements were found associated with high molecular weight DNA. Received: 9 July 1997 / Revision received: 29 January 1998 / Accepted: 29 November 1998     

A set of polymorphic EST‐derived markers for Picea species

A pooled DNA method was used to produce fully informative EST (expressed sequence tag)‐derived markers for the Picea genus. Nine markers were produced from 10 cDNA identified as candidates for cold tolerance or embryogenesis. Indels and SNPs (single nucleotide polymorphisms) were characterized from sequences obtained from pools of 10 individuals for each of the three species: Picea glauca (white spruce), Picea mariana (black spruce) and Picea abies (Norway spruce). Indels were present in 28% of the sequences and SNPs with a frequency greater than 10% were present on average in 1.2% of the positions.     

Asymmetry matters: A genomic assessment of directional biases in gene flow between hybridizing spruces

Assessing directional bias in interspecific gene flow might be important in determining the evolutionary trajectory of closely related species pairs. Using a set of 300 single nucleotide polymorphisms (SNPs) having variable propensity to cross species boundary, we evaluated the genomic extent and direction of interspecific gene flow in a progenitor‐derivative spruce species pair (black spruce and red spruce). A higher rate of gene flow was found from black spruce toward red spruce purebreds than vice versa. This asymmetry could reflect the historical gene flow between the two taxa at the time of species inception and during postglacial colonization. A clear asymmetry in introgression was depicted by a greater gene flow between red spruce and hybrids than between black spruce and hybrids. While backcrossing toward red spruce was invariably high across the genome, the actual species boundary is between hybrids and black spruce where gene flow is impeded at those genomic regions impermeable to introgression. Associations between hybrid index and climatic variables (total annual precipitation and mean annual temperature) were tested, as these might indicate a role for exogenous selection in maintaining the species boundary. While an apparent association was found between the hybrid index and precipitation, it collapsed when considered in light of the directional bias in interspecific gene flow. Hence, considering asymmetrical patterns of introgression allowed us to falsify an apparent role for exogenous selection. Although this was not formerly tested here, we suggest that this pattern could result from asymmetrical endogenous selection, a contention that deserves further investigations.     

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

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

Phylogeography of white spruce (Picea glauca) in eastern North America reveals contrasting ecological trajectories

Aim The main objective of this study was to investigate the origin of white spruce (Picea glauca) range discontinuities in eastern North America. A first analysis aimed at uncovering the glacial origin and post-glacial migration route of white spruce in the eastern part of the range. A second analysis aimed at evaluating whether disjunct white spruce populations in central Québec result from a recent expansion or are remnant stands following a northern range contraction. Location The first analysis covered the North American transcontinental boreal distribution of white spruce, and the second focused on the eastern part of its range. Methods A total of 589 P. glauca individuals from 51 populations, 24 of which were from a previous published source, were analysed for variation in three chloroplast DNA regions (trnT/L, trnL/F and ndhK/C). At the continental scale, haplotypic diversity, zones of genetic discontinuities, clustering analysis and analysis of molecular variance (AMOVA) were employed. In the eastern part of the range, haplotypic diversity was assessed with a permutation procedure, demographic history with mismatch analyses, and genetic structure with AMOVA. Results At the continental scale, various lines of evidence (genetic-diversity hotspots, endemic haplotypes, zones of genetic discontinuities, and pattern of molecular variation) show that different lineages migrated from three glacial refugia (Beringian, Mississippian, east Appalachian). Strong genetic discontinuities, low diversity and lack of evidence for demographic expansion were found among stands from isolated high hills in central Québec. Main conclusions White spruce in the eastern part of its range has an east Appalachian origin, a lineage genetically distinct from populations west of the Great Lakes. Previous studies have shown that the northward migration of white spruce is still ongoing as the northernmost populations expand into the tundra. In contrast, mixed white spruce and balsam fir (Abies balsamea) stands from isolated high hills in central Québec have already reached their Holocene maxima and are currently declining. This range contraction is probably a result of increased wildfire disturbance, which gradually creates a pattern of regional white spruce isolation and replacement by fire-prone black spruce (Picea mariana) forest.     

Inheritance of chloroplast and mitochondrial DNA in Picea and composition of hybrids from introgression zones

Summary The cloning of white spruce (Picea glauca) mitochondrial DNA homologous to the cytochrome oxidase II and ATPase genes of maize is described. These probes were used to define restriction fragment length polymorphisms which distinguish the white, Engelmann (P. engelmannii) and Sitka spruce (P. sitchensis) populations that occur in British Columbia. Analysis of progeny from crosses between the species revealed that mitochondrial DNA was maternally inherited in all cases (32 progeny from five independent crosses). The inheritance of chloroplast DNA was determined using a probe described previously; in this case, all progeny exhibited paternal inheritance (27 progeny from four crosses). Mitochondrial and chloroplast probes were used to test trees from zones of introgression between coastal (Sitka) and interior spruces (white and Engelmann). In most cases mitochondria and chloroplasts within individuals were contributed by different species. The data shows that there is a significant Sitka spruce component in trees east of the coastal watershed in British Columbia.     

Molecular cytogenetics of the genes encoding 18s-5.8s-26s rRNA and 5s rRNA in two species of spruce (Picea)

 Molecular cytogenetics is a convenient tool to investigate the organization and evolution of plant genomes. In coniferous trees of the Pinaceae, cytogenetic data is rudimentary since individual chromosomes are difficult to distinguish and karyotypes of related species are poorly differentiated. We determined the chromosomal locations of ribosomal RNA genes in white spruce (Picea glauca) and Sitka spruce (Picea sitchensis) using fluorescence in situ hybridization. The biotin-labeled DNA probes consisted of the 5s ribosomal DNA (rDNA) amplified from white spruce using the polymerase chain reaction and a heterologous 18s-5.8s-26s rDNA sequence. The 5s rDNA was present only on chromosome 5 at a single locus and near to an 18s-5.8s-26s rDNA locus in both species. Additional 18s-5.8s-26s rDNA loci were found at interstitial sites on six and four chromosomes of white and Sitka spruce, respectively, providing potentially useful interspecific differences. Progress in karyotyping both species is presented. A molecular analysis of 5s rDNA of white spruce revealed the presence of two classes of repeating units, one of 221 bp corresponding to the PCR amplification product, and another of approximately 600 bp. The nucleotide sequence and copy number of the 221-bp class is reported. Received: 17 September 1996/Accepted: 20 December 1996     

EVOLUTION OF POPULUS NIGRA (SECT. AIGEIROS):INTROGRESSIVE HYBRIDIZATION AND THE CHLOROPLAST CONTRIBUTION OF POPULUS ALBA (SECT. POPULUS)

Restriction site variation in chloroplast DNA and nuclear ribosomal DNA was examined in 16 accessions from the Salicaceae comprising ten species of Populus and one outgroup species of Salix. Forty-nine restriction site mutations in the chloroplast DNAs were used to generate one most parsimonious phylogenetic tree. This tree indicates that all varieties of P. nigra (black poplars of sect. Aigeiros) have a chloroplast genome, maternally inherited, derived from the clade including the white poplars (P. alba and segregate species of sect. Populus) and divergent from the American cottonwoods of their own section. Twenty-one restriction site mutations in the nuclear ribosomal DNAs generated a single most parsimonious phylogenetic tree that indicates that the nuclear genome ofP. nigra is distinct from both the white poplars and American cottonwoods. The incongruity of these independent molecular phylogenies provides evidence for an unusual origin of the black poplars. Populus alba or its immediate ancestor acted as the maternal parent in a hybridization event with the paternal lineage of P. nigra. Subsequent backcrosses to the paternal species gave rise to the extant P. nigra with a chloroplast genome of P. alba and the nuclear genome of the paternal species. These hybridization and introgression events must have pre-dated the divergence of the black poplar varieties. The biphyletic nature of the P. nigra genomes suggests that dependency on one class of molecular or morphological markers or the merging of the two kinds of data sets to derive accurate estimates of true phylogenies could be misleading in plants.     

Classifying seedlots of Picea sitchensis and P. glauca in zones of introgression using restriction analysis of chloroplast DNA

Summary Chloroplast DNA (cpDNA) restriction analysis was used to classify five reforestation seedlots as to species. The material included two Sitka spruce (Picea sitchensis (Bong.) Carr.), one white spruce (P. glauca (Moench) Voss) from interior British Columbia, and two putative hybrid seedlots from the coast-interior introgression zone in British Columbia. The cpDNA patterns generated by Bam-HI and Bc1-I from individual trees of Sitka spruce, white spruce, western white spruce (P. glauca var. albertiana (S. Brown)), and Engelmann spruce (P. engelmanni (Parry)) were species-specific. They were used as reference patterns for comparisons. In addition, two controlled crosses between white and Sitka spruce were analyzed to demonstrate the paternal inheritance of cpDNA in spruces. The cpDNA restriction patterns for the five seedlots were obtained from composite samples of seedlings from each lot and compared to the typical cpDNA patterns of each species. Restriction patterns for the two Sitka spruce seedlots agreed with those from the Sitka spruce tree, while patterns for the white spruce seedlots from British Columbia agreed with those from the white spruce tree, lacking evidence of any Engelmann spruce component in the sample. On the other hand, one putative hybrid seedlot showed cpDNA patterns similar to white spruce while the other showed fragments unique to both Sitka and white spruce, indicating that this was a hybrid seedlot. The analysis of cpDNA restriction polymorphism has proven to be an effective tool for classifying seedlots in regions of introgression. To our knowledge, these results provide the first demonstration of the use of cpDNA analysis for solving practical forestry problems.     

Natural hybridization between black spruce and red spruce

Using species-specific random amplified polymorphic DNA (RAPD) markers and morphological characters, natural hybridization between the closely related black spruce Picea mariana (Mill.) B.S.P. and red spruce P. rubens Sarg. was evaluated in natural populations from north-eastern North America. Sampling included populations from both areas of allopatry and also 14 populations from part of the area of sympatry located in the province of Québec and covering several thousands of square kilometres. Classification results from RAPD species-specific markers and from a discriminant function based on morphology were compared. Molecular analysis of the allopatric populations indicated a small amount of interspecific gene leakage with no asymmetric directionality to introgression. A high occurrence of hybrid/introgressant individuals was observed within sympatric populations, suggesting weak reproductive isolation. As expected, the detection of such individuals was more efficient using molecular markers than with morphological traits. The hybrid zone appeared extensive with variable species structure and, in some stands, the main component composed of hybrid/introgressant trees. Implications for the genecology and genetic management of these species are discussed.     

Exploring genomic variation associated with drought stress in Picea mariana populations

Predicted increases in drought and heat stress will likely induce shifts in species bioclimatic envelopes. Genetic variants adapted to water limitation may prove pivotal for species response under scenarios of increasing drought. In this study, we aimed to explore this hypothesis by investigating genetic variation in 16 populations of black spruce (Picea mariana) in relation to climate variables in Alaska. A total of 520 single nucleotide polymorphisms (SNPs) were genotyped for 158 trees sampled from areas of contrasting climate regimes. We used multivariate and univariate genotype‐by‐environment approaches along with available gene annotations to investigate the relationship between climate and genetic variation among sampled populations. Nine SNPs were identified as having a significant association with climate, of which five were related to drought stress response. Outlier SNPs with respect to the overall environment were significantly overrepresented for several biological functions relevant for coping with variable hydric regimes, including osmotic stress response. This genomic imprint is consistent with local adaptation of black spruce to drought stress. These results suggest that natural selection acting on standing variation prompts local adaptation in forest stands facing water limitation. Improved understanding of possible adaptive responses could inform our projections about future forest dynamics and help prioritize populations that harbor valuable genetic diversity for conservation.     

Occurrence of somaclonal variation among somatic embryo-derived white spruces (Picea glauca,Pinaceae)

Four different variegata phenotypes were identified among 2270 white spruce plants [Picea glauca (Moench.) Voss.] produced over a period of 2 yr from the time of induction of embryogenic tissues. The four variegated plants differed from each other in the extent and distribution of chlorophyll-deficient needles. Light microscopy showed that variegated leaves of a selected variant consisted of a chimeral mixture of green and white cells. Electron microscopy showed that cells in completely white needles had large nuclei with predominant euchromatin, lacked large cytoplasmic vacuoles, and harbored vacuolized plastids with aberrant morphologies. Various observations suggest that the recovered variegata phenotypes reflect some kind of genetic instability of either chloroplastic or nuclear genomes. To elucidate the genetic basis of these variegata phenotypes in white spruce, three out of the four variants were subjected to randomly amplified polymorphic DNA (RAPD) analysis. Out of more than 250 RAPD markers screened, only one correlated with white needles of variegated plants. The nucleotide sequence of this DNA fragment showed no homology with any known gene, but the amplified sequence appears most likely of nuclear origin.     

Hybrid status of Tibouchina urvilleana Cogn. revealed by ITS sequences of nuclear ribosomal DNA

Tibouchina urvilleana Cogn. is native to southern Brazil and currently cultivated as an important ornamental shrub in frost free areas around the world. Its rapid vegetative growth and sterility suggests that it might be of hybrid origin. In this study, Internal Transcribed Spacer (ITS) of nuclear ribosomal DNA and chloroplast trnL-trnF spacer from T. urvilleana and three other congeneric species were sequenced to test its hybrid status. Cloning sequencing revealed two distinct types of ITS sequences from T. urvilleana, with one type almost identical with Tibouchina aspera. Genetic distance between the two types was much larger than the average interspecific genetic distances calculated from other Tibouchina species. Sequencing of chloroplast trnL-trnF spacer showed that T. urvilleana has identical sequence with T. aspera, but differed from other congeneric species by one nucleotide substitution and two indels. Molecular data demonstrated clearly that T. urvilleana indeed was a hybrid, with T. aspera or closely related species acting as the maternal parent.     

New evidence from mitochondrial DNA of a progenitor-derivative species relationship between black spruce and red spruce (Pinaceae)

Mitochondrial DNA (mtDNA) markers were used to assess the genetic diversity in allopatric populations of black spruce (Picea mariana [Mill.] BSP) and red spruce (P. rubens Sarg.). Patterns of mitochondrial haplotypes (mitotypes) were strikingly different between the two species. All mtDNA markers surveyed were polymorphic in black spruce, revealing four different mitotypes and high levels of mtDNA diversity (P(p) = 100%, A = 2.0, H = 0.496). In contrast, populations of red spruce had only two mitotypes and harbored low levels of ggenetic diversity (P(p) = 13.2%, A = 1.1, H = 0.120). When the southernmost allopatric populations of red spruce were considered, only one mitotype was detected. As previously reported for nuclear gene loci, the diversity observed for mtDNA in red spruce was a subset of that found in black spruce. Comparison of present and previously published data supports the hypothesis of a recent progenitor-derivative relationship between these species, red spruce presumably being derived by allopatric speciation of an isolated population of black spruce during the Pleistocene.     

Frequency of somaclonal variation in plants of black spruce (Picea mariana,Pinaceae) and white spruce (P. glauca,Pinaceae) derived from somatic embryogenesis and identification of some factors involved in genetic instability

Plants of black spruce (Picea mariana, N = 7047 individuals) and white spruce (P. glauca, N = 3995 individuals) were regenerated from a total of 87 clones over a 5-yr period by somatic embryogenesis to study factors that might be associated with the appearance of variant phenotypes. Morphological evaluation of the plants showed several types of variation. These variations were grouped into nine types: dwarfism (type A), reduced height with various form anomalies (types B, C, and D), needle fasciation (type E), abnormality in tree architecture (type F), variegata phenotype (type G), and plants with an overall regular morphology but smaller than normal plants (type H). Plagiotropic plants were also observed (type I). Each plant from types A to H (except type C where no plants survived more than 6 mo) had retained its phenotype over 4–5 yr of growth. Some of the variant types could be related to chromosomic instability: chromosome counts showed aneuploid cells for type-A and type-D plants. The type I (plagiotropism) was not related to genetic instability but rather to physiological disorders. In total, spruce variants of types A–H were obtained at relatively low frequencies, i.e., 1.0% (39/3995) for white spruce and 1.6% (110/7047) for black spruce. Statistical analyses, conducted with family, clone, and time in maintenance as variables, showed that clone was the most important source of genetic instability followed by time in maintenance.     

Barcoding Poplars (Populus L.) from Western China

Background

Populus is an ecologically and economically important genus of trees, but distinguishing between wild species is relatively difficult due to extensive interspecific hybridization and introgression, and the high level of intraspecific morphological variation. The DNA barcoding approach is a potential solution to this problem.

Methodology/Principal Findings

Here, we tested the discrimination power of five chloroplast barcodes and one nuclear barcode (ITS) among 95 trees that represent 21 Populus species from western China. Among all single barcode candidates, the discrimination power is highest for the nuclear ITS, progressively lower for chloroplast barcodes matK (M), trnG-psbK (G) and psbK-psbI (P), and trnH-psbA (H) and rbcL (R); the discrimination efficiency of the nuclear ITS (I) is also higher than any two-, three-, or even the five-locus combination of chloroplast barcodes. Among the five combinations of a single chloroplast barcode plus the nuclear ITS, H+I and P+I differentiated the highest and lowest portion of species, respectively. The highest discrimination rate for the barcodes or barcode combinations examined here is 55.0% (H+I), and usually discrimination failures occurred among species from sympatric or parapatric areas.

Conclusions/Significance

In this case study, we showed that when discriminating Populus species from western China, the nuclear ITS region represents a more promising barcode than any maternally inherited chloroplast region or combination of chloroplast regions. Meanwhile, combining the ITS region with chloroplast regions may improve the barcoding success rate and assist in detecting recent interspecific hybridizations. Failure to discriminate among several groups of Populus species from sympatric or parapatric areas may have been the result of incomplete lineage sorting, frequent interspecific hybridizations and introgressions. We agree with a previous proposal for constructing a tiered barcoding system in plants, especially for taxonomic groups that have complex evolutionary histories (e.g. Populus).