High-density genetic linkage mapping in Sitka spruce advances the integration of genomic resources in conifers.

In species with large and complex genomes such as conifers, dense linkage maps are a useful resource for supporting genome assembly and laying the genomic groundwork at the structural, populational, and functional levels. However, most of the 600+ extant conifer species still lack extensive genotyping resources, which hampers the development of high-density linkage maps. In this study, we developed a linkage map relying on 21,570 single nucleotide polymorphism (SNP) markers in Sitka spruce (Picea sitchensis [Bong.] Carr.), a long-lived conifer from western North America that is widely planted for productive forestry in the British Isles. We used a single-step mapping approach to efficiently combine RAD-seq and genotyping array SNP data for 528 individuals from 2 full-sib families. As expected for spruce taxa, the saturated map contained 12 linkages groups with a total length of 2,142 cM. The positioning of 5,414 unique gene coding sequences allowed us to compare our map with that of other Pinaceae species, which provided evidence for high levels of synteny and gene order conservation in this family. We then developed an integrated map for P. sitchensis and Picea glauca based on 27,052 markers and 11,609 gene sequences. Altogether, these 2 linkage maps, the accompanying catalog of 286,159 SNPs and the genotyping chip developed, herein, open new perspectives for a variety of fundamental and more applied research objectives, such as for the improvement of spruce genome assemblies, or for marker-assisted sustainable management of genetic resources in Sitka spruce and related species.

Spruce giga-genomes: structurally similar yet distinctive with differentially expanding gene families and rapidly evolving genes.

Spruces (Picea spp.) are coniferous trees widespread in boreal and mountainous forests of the northern hemisphere, with large economic significance and enormous contributions to global carbon sequestration. Spruces harbor very large genomes with high repetitiveness, hampering their comparative analysis. Here, we present and compare the genomes of four different North American spruces: the genome assemblies for Engelmann spruce (Picea engelmannii) and Sitka spruce (Picea sitchensis) together with improved and more contiguous genome assemblies for white spruce (Picea glauca) and for a naturally occurring introgress of these three species known as interior spruce (P. engelmannii × glauca × sitchensis). The genomes were structurally similar, and a large part of scaffolds could be anchored to a genetic map. The composition of the interior spruce genome indicated asymmetric contributions from the three ancestral genomes. Phylogenetic analysis of the nuclear and organelle genomes revealed a topology indicative of ancient reticulation. Different patterns of expansion of gene families among genomes were observed and related with presumed diversifying ecological adaptations. We identified rapidly evolving genes that harbored high rates of non-synonymous polymorphisms relative to synonymous ones, indicative of positive selection and its hitchhiking effects. These gene sets were mostly distinct between the genomes of ecologically contrasted species, and signatures of convergent balancing selection were detected. Stress and stimulus response was identified as the most frequent function assigned to expanding gene families and rapidly evolving genes. These two aspects of genomic evolution were complementary in their contribution to divergent evolution of presumed adaptive nature. These more contiguous spruce giga-genome sequences should strengthen our understanding of conifer genome structure and evolution, as their comparison offers clues into the genetic basis of adaptation and ecology of conifers at the genomic level. They will also provide tools to better monitor natural genetic diversity and improve the management of conifer forests. The genomes of four closely related North American spruces indicate that their high similarity at the morphological level is paralleled by the high conservation of their physical genome structure. Yet, the evidence of divergent evolution is apparent in their rapidly evolving genomes, supported by differential expansion of key gene families and large sets of genes under positive selection, largely in relation to stimulus and environmental stress response.

Combining QTL Mapping and Transcriptomics to Decipher the Genetic Architecture of Phenolic Compounds Metabolism in the Conifer White Spruce.

Conifer forests worldwide are becoming increasingly vulnerable to the effects of climate change. Although the production of phenolic compounds (PCs) has been shown to be modulated by biotic and abiotic stresses, the genetic basis underlying the variation in their constitutive production level remains poorly documented in conifers. We used QTL mapping and RNA-Seq to explore the complex polygenic network underlying the constitutive production of PCs in a white spruce (Picea glauca) full-sib family for 2 years. QTL detection was performed for nine PCs and differentially expressed genes (DEGs) were identified between individuals with high and low PC contents for five PCs exhibiting stable QTLs across time. A total of 17 QTLs were detected for eight metabolites, including one major QTL explaining up to 91.3% of the neolignan-2 variance. The RNA-Seq analysis highlighted 50 DEGs associated with phenylpropanoid biosynthesis, several key transcription factors, and a subset of 137 genes showing opposite expression patterns in individuals with high levels of the flavonoids gallocatechin and taxifolin glucoside. A total of 19 DEGs co-localized with QTLs. Our findings represent a significant step toward resolving the genomic architecture of PC production in spruce and facilitate the functional characterization of genes and transcriptional networks responsible for differences in constitutive production of PCs in conifers.

Connecting tree-ring phenotypes, genetic associations and transcriptomics to decipher the genomic architecture of drought adaptation in a widespread conifer.

As boreal forests face significant threats from climate change, understanding evolutionary trajectories of coniferous species has become fundamental to adapting management and conservation to a drying climate. We examined the genomic architecture underlying adaptive variation related to drought tolerance in 43 populations of a widespread boreal conifer, white spruce (Picea glauca [Moench] Voss), by combining genotype-environment associations, genotype-phenotype associations, and transcriptomics. Adaptive genetic variation was identified by correlating allele frequencies for 6,153 single nucleotide polymorphisms from 2,606 candidate genes with temperature, precipitation and aridity gradients, and testing for significant associations between genotypes and 11 dendrometric and drought-related traits (i.e., anatomical, growth response and climate-sensitivity traits) using a polygenic model. We identified a set of 285 genes significantly associated with a climatic factor or a phenotypic trait, including 110 that were differentially expressed in response to drought under greenhouse-controlled conditions. The interlinked phenotype-genotype-environment network revealed eight high-confidence genes involved in white spruce adaptation to drought, of which four were drought-responsive in the expression analysis. Our findings represent a significant step toward the characterization of the genomic basis of drought tolerance and adaptation to climate in conifers, which is essential to enable the establishment of resilient forests in view of new climate conditions.

Genomic prediction for hastening and improving efficiency of forward selection in conifer polycross mating designs: an example from white spruce.

Genomic selection (GS) has a large potential for improving the prediction accuracy of breeding values and significantly reducing the length of breeding cycles. In this context, the choice of mating designs becomes critical to improve the efficiency of breeding operations and to obtain the largest genetic gains per time unit. Polycross mating designs have been traditionally used in tree and plant breeding to perform backward selection of the female parents. The possibility to use genetic markers for paternity identification and for building genomic prediction models should allow for a broader use of polycross tests in forward selection schemes. We compared the accuracies of genomic predictions of offspring's breeding values from a polycross and a full-sib (partial diallel) mating design with similar genetic background in white spruce (Picea glauca). Trees were phenotyped for growth and wood quality traits, and genotyped for 4092 SNPs representing as many gene loci distributed across the 12 spruce chromosomes. For the polycross progeny test, heritability estimates were smaller, but more precise using the genomic BLUP (GBLUP) model as compared with pedigree-based models accounting for the maternal pedigree or for the reconstructed full pedigree. Cross-validations showed that GBLUP predictions were 22-52% more accurate than predictions based on the maternal pedigree, and 5-7% more accurate than predictions using the reconstructed full pedigree. The accuracies of GBLUP predictions were high and in the same range for most traits between the polycross (0.61-0.70) and full-sib progeny tests (0.61-0.74). However, higher genetic gains per time unit were expected from the polycross mating design given the shorter time needed to conduct crosses. Considering the operational advantages of the polycross design in terms of easier handling of crosses and lower associated costs for test establishment, we believe that this mating scheme offers great opportunities for the development and operational application of forward GS.

A catalog of annotated high-confidence SNPs from exome capture and sequencing reveals highly polymorphic genes in Norway spruce (Picea abies).

BACKGROUND: Norway spruce [Picea abies (L.) Karst.] is ecologically and economically one of the most important conifer worldwide. Our main goal was to develop a large catalog of annotated high confidence gene SNPs that should sustain the development of genomic tools for the conservation of natural and domesticated genetic diversity resources, and hasten tree breeding efforts in this species. RESULTS: Targeted sequencing was achieved by capturing P. abies exome with probes previously designed from the sequenced transcriptome of white spruce (Picea glauca (Moench) Voss). Capture efficiency was high (74.5%) given a high level of exome conservation between the two species. Using stringent criteria, we delimited a set of 61,771 high-confidence SNPs across 13,543 genes. To validate SNPs, a high-throughput genotyping array was developed for a subset of 5571 predicted SNPs representing as many different gene loci, and was used to genotype over 1000 trees. The estimated true positive rate of the resource was 84.2%, which was comparable with the genotyping success rate obtained for P. abies control SNPs recycled from previous genotyping efforts. We also analyzed SNP abundance across various gene functional categories. Several GO terms and gene families involved in stress response were found over-represented in highly polymorphic genes. CONCLUSION: The annotated high-confidence SNP catalog developed herein represents a valuable genomic resource, being representative of over 13 K genes distributed across the P. abies genome. This resource should serve a variety of population genomics and breeding applications in Norway spruce.

Genetic Adaptation to Climate in White Spruce Involves Small to Moderate Allele Frequency Shifts in Functionally Diverse Genes.

Understanding the genetic basis of adaptation to climate is of paramount importance for preserving and managing genetic diversity in plants in a context of climate change. Yet, this objective has been addressed mainly in short-lived model species. Thus, expanding knowledge to nonmodel species with contrasting life histories, such as forest trees, appears necessary. To uncover the genetic basis of adaptation to climate in the widely distributed boreal conifer white spruce (Picea glauca), an environmental association study was conducted using 11,085 single nucleotide polymorphisms representing 7,819 genes, that is, approximately a quarter of the transcriptome.Linear and quadratic regressions controlling for isolation-by-distance, and the Random Forest algorithm, identified several dozen genes putatively under selection, among which 43 showed strongest signals along temperature and precipitation gradients. Most of them were related to temperature. Small to moderate shifts in allele frequencies were observed. Genes involved encompassed a wide variety of functions and processes, some of them being likely important for plant survival under biotic and abiotic environmental stresses according to expression data. Literature mining and sequence comparison also highlighted conserved sequences and functions with angiosperm homologs.Our results are consistent with theoretical predictions that local adaptation involves genes with small frequency shifts when selection is recent and gene flow among populations is high. Accordingly, genetic adaptation to climate in P. glauca appears to be complex, involving many independent and interacting gene functions, biochemical pathways, and processes. From an applied perspective, these results shall lead to specific functional/association studies in conifers and to the development of markers useful for the conservation of genetic resources.

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.

Integrating phylogeography and paleoecology to investigate the origin and dynamics of hybrid zones: insights from two widespread North American firs.

Secondary contact between closely related taxa routinely occurs during postglacial migrations. After initial contact, the location of hybrid zones may shift geographically or remain spatially stable over time in response to various selective pressures or neutral processes. Studying the extent and direction of introgression using markers having contrasted levels of gene flow can help unravel the historical dynamics of hybrid zones. Thanks to their contrasted maternal and paternal inheritance, resulting in different levels of gene flow for mitochondrial and chloroplast DNA (mtDNA and cpDNA), the Pinaceae stand out as a relevant biological model for this purpose. The objective of the study was to assess whether the hybrid zone between Abies balsamea and Abies lasiocarpa (two largely distributed Pinaceae) has moved or remained stable over time by analysing the distribution of cytoplasmic DNA variation as well as published palaeobotanical data. Interspecific gene flow was higher for cpDNA than mtDNA markers; hence, the geographic distribution of mitotypes was more congruent with species distributions than chlorotypes. This genetic signature was contrary to expectations under a moving hybrid zone scenario, as well as empirical observations in other conifers. Genetic evidence for this rare instance of stable hybrid zone was corroborated by the colonization chronology derived from published fossil data, indicating that the two fir species initially came into contact in the area corresponding to the current sympatric zone 11 kyr ago. While an explanatory analysis suggested the putative influence of various environmental factors on the relative abundance of cytoplasmic genome combinations, further research appears necessary to assess the role of both demographic history and selective factors in driving the dynamics of hybrid zones.

Less pollen-mediated gene flow for more signatures of glacial lineages: congruent evidence from balsam fir cpDNA and mtDNA for multiple refugia in eastern and central North America.

The phylogeographic structure and postglacial history of balsam fir (Abies balsamea), a transcontinental North American boreal conifer, was inferred using mitochondrial DNA (mtDNA) and chloroplast DNA (cpDNA) markers. Genetic structure among 107 populations (mtDNA data) and 75 populations (cpDNA data) was analyzed using Bayesian and genetic distance approaches. Population differentiation was high for mtDNA (dispersed by seeds only), but also for cpDNA (dispersed by seeds and pollen), indicating that pollen gene flow is more restricted in balsam fir than in other boreal conifers. Low cpDNA gene flow in balsam fir may relate to low pollen production due to the inherent biology of the species and populations being decimated by recurrent spruce budworm epidemics, and/or to low dispersal of pollen grains due to their peculiar structural properties. Accordingly, a phylogeographic structure was detected using both mtDNA and cpDNA markers and population structure analyses supported the existence of at least five genetically distinct glacial lineages in central and eastern North America. Four of these would originate from glacial refugia located south of the Laurentide ice sheet, while the last one would have persisted in the northern Labrador region. As expected due to reduced pollen-mediated gene flow, congruence between the geographic distribution of mtDNA and cpDNA lineages was higher than in other North American conifers. However, concordance was not complete, reflecting that restricted but nonetheless detectable cpDNA gene flow among glacial lineages occurred during the Holocene. As a result, new cpDNA and mtDNA genome combinations indicative of cytoplasmic genome capture were observed.

Parallel and lineage-specific molecular adaptation to climate in boreal black spruce.

In response to selective pressure, adaptation may follow different genetic pathways throughout the natural range of a species due to historical differentiation in standing genetic variation. Using 41 populations of black spruce (Picea mariana), the objectives of this study were to identify adaptive genetic polymorphisms related to temperature and precipitation variation across the transcontinental range of the species, and to evaluate the potential influence of historical events on their geographic distribution. Population structure was first inferred using 50 control nuclear markers. Then, 47 candidate gene SNPs identified in previous genome scans were tested for relationship with climatic factors using an F(ST) -based outlier method and regressions between allele frequencies and climatic variations. Two main intraspecific lineages related to glacial vicariance were detected at the transcontinental scale. Within-lineage analyses of allele frequencies allowed the identification of 23 candidate SNPs significantly related to precipitation and/or temperature variation, among which seven were common to both lineages, eight were specific to the eastern lineage and eight were specific to the western lineage. The implication of these candidate SNPs in adaptive processes was further supported by gene functional annotations. Multiple evidences indicated that the occurrence of lineage-specific adaptive SNPs was better explained by selection acting on historically differentiated gene pools rather than differential selection due to heterogeneity of interacting environmental factors and pleiotropic effects. Taken together, these findings suggest that standing genetic variation of potentially adaptive nature has been modified by historical events, hence affecting the outcome of recent selection and leading to different adaptive routes between intraspecific lineages.

From glacial refugia to modern populations: new assemblages of organelle genomes generated by differential cytoplasmic gene flow in transcontinental black spruce.

Assessing species' range-wide cytoplasmic diversity provides valuable insights regarding their dispersal and adaptive potential in a changing environment. Transcontinental chloroplast (cpDNA) and mitochondrial DNA (mtDNA) population structures were compared to identify putative ancestral and new cytoplasmic genome assemblages in black spruce (Picea mariana), a North American boreal conifer. Mean within-population diversity and allelic richness for cpSSR markers were 0.80 and 4.21, respectively, and diminished westward. Population differentiation based on G(ST) was lower for cpDNA than for mtDNA (G(ST) =0.104 and 0.645, respectively) but appeared comparable when estimated using Jost differentiation index (D=0.459 and 0.537, respectively). Further analyses resulted in the delineation of at least three genetically distinct cpDNA lineages partially congruent with those inferred from mtDNA data, which roughly corresponded to western, central and eastern Canada. Additionally, the patterns of variation in Alaska for both cpDNA and mtDNA markers suggested that black spruce survived the last glacial maximum in this northern region. The range-wide comparison of the geographic extent of cytoplasmic DNA lineages revealed that extensive pollen gene flow between ancestral lineages occurred preferentially from west to east during the postglacial expansion of the species, while seed-mediated gene flow remained geographically restricted. This differential gene flow promoted intraspecific cytoplasmic capture that generated new assemblages of cpDNA and mtDNA genomes during the Holocene. Hence, black spruce postglacial colonization unexpectedly resulted in an increase in genetic diversity with possible adaptive consequences.