Juan P. Jaramillo-Correa

Variation in mitochondrial DNA reveals multiple distant glacial refugia in black spruce (Picea mariana), a transcontinental North American conifer

Range‐wide genetic variation of black spruce (Picea mariana) was studied using polymerase chain reaction‐random fragment length polymorphism markers of the mitochondrial genome. Four polymorphic mitochondrial DNA (mtDNA) loci were surveyed and two or three alleles were detected at each locus, resulting in 10 multilocus mtDNA types or mitotypes. A significant subdivision of population genetic diversity was detected (GST = 0.671; NST = 0.726), suggesting low levels of gene flow among populations. The distribution of mitotypes was not random (NST > GST; P < 0.05) and revealed four partially overlapping zones, presumably representative of different glacial populations. Comparison of the genetic structure derived from mtDNA markers and the colonization paths previously deduced from the fossil and pollen records allow us to infer at least three southern and one northeastern glacial populations for black spruce. The patterns revealed in this study suggest that black spruce shares its biogeographical history with other forest‐associated North American species.

A mitochondrial DNA minisatellite reveals the postglacial history of jack pine (Pinus banksiana), a broad‐range North American conifer

Jack pine (Pinus banksiana Lamb.) is a broadly distributed North American conifer and its current range was covered by the Laurentian ice sheet during the last glacial maximum. To infer about the history and postglacial colonization of this boreal species, range‐wide genetic variation was assessed using a new and highly variable minisatellite‐like marker of the mitochondrial genome. Among the 543 trees analysed, 14 distinct haplotypes were detected, which corresponded to different repeat numbers of the 32‐nucleotide minisatellite‐like motif. Several haplotypes were rare with limited distribution, suggesting recent mutation events during the Holocene. At the population level, an average of 2.6 haplotypes and a mean haplotype diversity (H) of 0.328 were estimated. Population subdivision of genetic diversity was quite high with GST and RST values of 0.569 and 0.472, respectively. Spatial analyses identified three relatively homogeneous groups of populations presumably representative of genetically distinct glacial populations, one west and one east of the Appalachian Mountains in the United States and a third one presumably on the unglaciated northeastern coastal area in Canada. These results indicate the significant role of the northern part of the US Appalachian Mountains as a factor of vicariance during the ice age. A fourth distinct group of populations was observed in central Québec where the continental glacier retreated last. It included populations harbouring haplotypes present into the three previous groups, and it had higher level of haplotype diversity per population (H = 0.548) and lower population differentiation (GST = 0.265), which indicates a zone of suture or secondary contact between the migration fronts of the three glacial populations. Introgression from Pinus contorta Dougl. var. latifolia Engelm. was apparent in one western population from Alberta. Altogether, these results indicate that the mitochondrial DNA variation of jack pine is geographically highly structured and it correlates well with large‐scale patterns emerging from recent phylogeographical studies of other tree boreal species in North America.

Inferring the past from the present phylogeographic structure of North American forest trees: seeing the forest for the genes.

The study of past historical events that have led to ecological changes is a recurrent topic in many disciplines. Given that many of these events have left a large and long-lasting evolutionary imprint on the extant population genetic structure of species, phylogeographic studies on modern taxa have been largely used to infer the impacts of these events and to complement previous paleoecological and paleobotanical surveys. In spite of the geographical and geological complexity of North America, converging patterns can be observed when comparing the available genetic data for forest trees. Such patterns include the co-location of genetic discontinuities among species and their coincidence with mountain ranges (e.g., the Appalachians, the Rocky Mountains, the Sierra Nevada, or the Transverse Volcanic Belt) and with previously inferred glacial refugia. Using examples drawn from the available literature, we illustrate such shared features and present the contrasting phylogeographic patterns observed among the...

Contrasting evolutionary forces driving population structure at expressed sequence tag polymorphisms, allozymes and quantitative traits in white spruce

Patterns of variation in quantitative characters and genetic markers were compared among six regional populations of white spruce [Picea glauca (Moench) Voss]. Although some phenotypic characters were correlated with latitude (r = 0.791), longitude (r = −0.796) and precipitation during the growing season (r = 0.789), variability at genetic markers was not correlated with geographical or bioclimatic variables, and followed neutral expectations. Estimates of genetic diversity and population differentiation for 14 allozymes (translated regions of coding genes) were essentially indistinguishable from those observed for 11 expressed sequence tag polymorphisms (ESTPs) from untranslated regions of coding genes. Variation among populations for quantitative traits such as eighth year height (QST = 0.082), thirteenth year height (QST = 0.069), total wood density (QST = 0.102) and date of budset (QST = 0.246), was greater than for allozymes (GST = 0.014) and ESTPs (GST = 0.019). These trends suggest a strong adaptive response in quantitative traits, contrasting to allozymes and ESTPs where no selective response could be detected and where populations appeared to be essentially in a migration–drift equilibrium.

Decoupled mitochondrial and chloroplast DNA population structure reveals Holocene collapse and population isolation in a threatened Mexican-endemic conifer.

Chihuahua spruce (Picea chihuahuana Martínez) is a montane subtropical conifer endemic to the Sierra Madre Occidental in northwestern México. Range‐wide variation was investigated using maternally inherited mitochondrial (mtDNA) and paternally inherited chloroplast (cpDNA) DNA markers. Among the 16 mtDNA regions analysed, only two mitotypes were detected, while the study of six cpDNA microsatellite markers revealed eight different chlorotypes. The average cpDNA diversity (H = 0.415) was low but much higher than that for mtDNA (H = 0). The distribution of mitotypes revealed two clear nonoverlapping areas (GST = NST = 1), one including northern populations and the second one including the southern and central stands, suggesting that these two regions may represent different ancestral populations. The cpDNA markers showed lower population differentiation (GST = 0.362; RST = 0.230), implying that the two ancestral populations continued to exchange pollen after their initial geographic separation. A lack of a phylogeographic structure was revealed by different spatial analyses of cpDNA (GST > RST; and samova), and reduced cpDNA gene flow was noted among populations (Nm = 0.873). Some stands deviated significantly from the mutation–drift equilibrium, suggesting recent bottlenecks. Altogether, these various trends are consistent with the hypothesis of a population collapse during the Holocene warming and suggest that most of the modern P. chihuahuana populations are now effectively isolated with their genetic diversity essentially modelled by genetic drift. The conservation efforts should focus on most southern populations and on the northern and central stands exhibiting high levels of genetic diversity. Additional mtDNA sequence analysis confirmed that P. martinezii (Patterson) is not conspecific with P. chihuahuana, and thus deserves separate conservation efforts.

Ancestry and divergence of subtropical montane forest isolates: molecular biogeography of the genus Abies (Pinaceae) in southern México and Guatemala

The genus Abies has a complex history in southern México and Guatemala. In this region, four closely related species, Abies flinckii, A. guatemalensis, A. hickelii, and A. religiosa, are distributed in fragmented and isolated montane populations. Range‐wide genetic variation was investigated across species using cytoplasmic DNA markers with contrasted inheritance. Variation at two maternally inherited mitochondrial DNA markers was low. All species shared two of the nine mitotypes detected, while the remaining seven mitochondrial DNA types were restricted to a few isolated stands. Mitochondrial genetic differentiation across taxa was high (GST = 0.933), it was not related to the taxonomic identity (amova; P > 0.05) of the populations, and it was not phylogeographically structured (GST ≈ NST). In contrast, variation at three paternally inherited chloroplast DNA microsatellites was high. Chloroplast genetic differentiation was lower (GST = 0.402; RST = 0.547) than for mitochondrial DNA, but it was significantly related to taxonomy (amova; P < 0.001), and exhibited a significant phylogeographical structure (GST < RST). Different analyses of population structure indicated that A. flinckii was the most divergent taxon, while the remaining three species formed a relatively homogeneous group. However, a small number of the populations of these three taxa, all located at the limits of their respective ranges or in the Transverse Volcanic Belt, diverged from this main cluster. These trends suggest that the Mesoamerican Abies share a recent common ancestor and that their divergence and speciation is mainly driven by genetic drift and isolation during the warm interglacial periods.

Diverging patterns of mitochondrial and nuclear DNA diversity in subarctic black spruce: imprint of a founder effect associated with postglacial colonization

High‐latitude ecotonal populations at the species margins may exhibit altered patterns of genetic diversity, resulting from more or less recent founder events and from bottleneck effects in response to climate oscillations. Patterns of genetic diversity were investigated in nine populations of the conifer black spruce (Picea mariana[Mill.] BSP.) in northwestern Québec, Canada, using seed‐dispersed mitochondrial (mt) DNA and nuclear (nc) DNA. mtDNA diversity (mitotypes) was assessed at three loci, and ncDNA diversity was estimated for nine expressed sequence tag polymorphism (ESTP) loci. Sampling included populations from the boreal forest and the southern and northern subzones of the subarctic forest‐tundra, a fire‐born ecotone. For ncDNA, populations from all three vegetation zones were highly diverse with little population differentiation (θN = 0.014); even the northernmost populations showed no loss of rare alleles. Patterns of mitotype diversity were strikingly different: within‐population diversity and population differentiation were high for boreal forest populations [expected heterozygosity per locus (HE) = 0.58 and θM = 0.529], but all subarctic populations were fixed for a single mitotype (HE = 0). This lack of variation suggests a founder event caused by long‐distance seed establishment during postglacial colonization, consistent with palaeoecological data. The estimated movement of seeds alone (effective number of migrants per generation, NmM < 2) was much restricted compared to that estimated from nuclear variants, which including pollen movement (NmN > 17). This could account for the conservation of a founder imprint in the mtDNA of subarctic black spruce. After reduction, presumably in the early Holocene, the diversity in ncDNA would have been replenished rapidly by pollen‐mediated gene flow, and maintained subsequently through vegetative layering during the current cooler period covering the last 3000 years.

Development and implementation of a highly- multiplexed SNP array for genetic mapping in maritime pine and comparative mapping with loblolly pine

BackgroundSingle nucleotide polymorphisms (SNPs) are the most abundant source of genetic variation among individuals of a species. New genotyping technologies allow examining hundreds to thousands of SNPs in a single reaction for a wide range of applications such as genetic diversity analysis, linkage mapping, fine QTL mapping, association studies, marker-assisted or genome-wide selection. In this paper, we evaluated the potential of highly-multiplexed SNP genotyping for genetic mapping in maritime pine (Pinus pinaster Ait.), the main conifer used for commercial plantation in southwestern Europe.ResultsWe designed a custom GoldenGate assay for 1,536 SNPs detected through the resequencing of gene fragments (707 in vitro SNPs/Indels) and from Sanger-derived Expressed Sequenced Tags assembled into a unigene set (829 in silico SNPs/Indels). Offspring from three-generation outbred (G2) and inbred (F2) pedigrees were genotyped. The success rate of the assay was 63.6% and 74.8% for in silico and in vitro SNPs, respectively. A genotyping error rate of 0.4% was further estimated from segregating data of SNPs belonging to the same gene. Overall, 394 SNPs were available for mapping. A total of 287 SNPs were integrated with previously mapped markers in the G2 parental maps, while 179 SNPs were localized on the map generated from the analysis of the F2 progeny. Based on 98 markers segregating in both pedigrees, we were able to generate a consensus map comprising 357 SNPs from 292 different loci. Finally, the analysis of sequence homology between mapped markers and their orthologs in a Pinus taeda linkage map, made it possible to align the 12 linkage groups of both species.ConclusionsOur results show that the GoldenGate assay can be used successfully for high-throughput SNP genotyping in maritime pine, a conifer species that has a genome seven times the size of the human genome. This SNP-array will be extended thanks to recent sequencing effort using new generation sequencing technologies and will include SNPs from comparative orthologous sequences that were identified in the present study, providing a wider collection of anchor points for comparative genomics among the conifers.

Mitochondrial genome recombination in the zone of contact between two hybridizing conifers.

Variation in mitochondrial DNA was surveyed at four gene loci in and around the zone of contact between two naturally hybridizing conifers, black spruce (Picea mariana) and red spruce (P. rubens) in northeastern North America. Most of the mtDNA diversity of these species was found in populations next to or into the zone of contact, where some individuals bore rare mitotypes intermediate between the common mitotypes observed in the allopatric areas of each species. Sequence analysis and tests for mtDNA recombination point to this phenomenon, rather than to recurrent mutation, as the most tenable hypothesis for the origin of these rare mitotypes. From the 10 mitotypes observed, at least 4 would be the product of recombination between 4 of the 5 putative ancestral mitotypes. Tests for cytonuclear disequilibrium and geographical structure of the putative recombinant mitotypes suggest that mtDNA recombination is not frequent and relatively recent on the geological time scale. mtDNA recombination would have been promoted by transient heteroplasmy due to leakage of paternal mtDNA since the Holocene secondary contact between the two species.

Development of high-density SNP genotyping arrays for white spruce (Picea glauca) and transferability to subtropical and nordic congeners.

High‐density SNP genotyping arrays can be designed for any species given sufficient sequence information of high quality. Two high‐density SNP arrays relying on the Infinium iSelect technology (Illumina) were designed for use in the conifer white spruce (Picea glauca). One array contained 7338 segregating SNPs representative of 2814 genes of various molecular functional classes for main uses in genetic association and population genetics studies. The other one contained 9559 segregating SNPs representative of 9543 genes for main uses in population genetics, linkage mapping of the genome and genomic prediction. The SNPs assayed were discovered from various sources of gene resequencing data. SNPs predicted from high‐quality sequences derived from genomic DNA reached a genotyping success rate of 64.7%. Nonsingleton in silico SNPs (i.e. a sequence polymorphism present in at least two reads) predicted from expressed sequenced tags obtained with the Roche 454 technology and Illumina GAII analyser resulted in a similar genotyping success rate of 71.6% when the deepest alignment was used and the most favourable SNP probe per gene was selected. A variable proportion of these SNPs was shared by other nordic and subtropical spruce species from North America and Europe. The number of shared SNPs was inversely proportional to phylogenetic divergence and standing genetic variation in the recipient species, but positively related to allele frequency in P. glauca natural populations. These validated SNP resources should open up new avenues for population genetics and comparative genetic mapping at a genomic scale in spruce species.