A white spruce gene catalog for conifer genome analyses

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.

Genome-wide identification of late embryogenesis abundant protein family and their key regulatory network in Pinus tabuliformis cold acclimation

Cold acclimation is a crucial biological process that enables conifers to overwinter safely. The late embryogenesis abundant (LEA) protein family plays a pivotal role in enhancing freezing tolerance during this process. Despite its importance, the identification, molecular functions and regulatory networks of the LEA protein family have not been extensively studied in conifers or gymnosperms. Pinus tabuliformis, a conifer with high ecological and economic values and with high-quality genome sequence, is an ideal candidate for such studies. Here, a total of 104 LEA genes were identified from P. tabuliformis, and we renamed them according to their subfamily group: PtLEA1-PtLEA92 (group LEA1-LEA6), PtSMP1-PtSMP6 (group seed maturation protein) and PtDHN1-PtDHN6 (group Dehydrin). While the sequence structure of P. tabuliformis  LEA genes are conserved, their physicochemical properties exhibit unique characteristics within different subfamily groupings. Notably, the abundance of low-temperature responsive elements in PtLEA genes was observed. Using annual rhythm and temperature gradient transcriptome data, PtLEA22 was identified as a key gene that responds to low-temperature induction while conforming to the annual cycle of cold acclimation. Overexpression of PtLEA22 enhanced Arabidopsis freezing tolerance. Furthermore, several transcription factors potentially co-expressed with PtLEA22 were validated using yeast one-hybrid and dual-luciferase assays, revealing that PtDREB1 could directly bind PtLEA22 promoter to positively regulate its expression. These findings reveal the genome-wide characterization of P. tabuliformis  LEA genes and their importance in the cold acclimation, while providing a theoretical basis for studying the molecular mechanisms of cold acclimation in conifers.

What can hornworts teach us?

The hornworts are a small group of land plants, consisting of only 11 families and approximately 220 species. Despite their small size as a group, their phylogenetic position and unique biology are of great importance. Hornworts, together with mosses and liverworts, form the monophyletic group of bryophytes that is sister to all other land plants (Tracheophytes). It is only recently that hornworts became amenable to experimental investigation with the establishment of Anthoceros agrestis as a model system. In this perspective, we summarize the recent advances in the development of A. agrestis as an experimental system and compare it with other plant model systems. We also discuss how A. agrestis can help to further research in comparative developmental studies across land plants and to solve key questions of plant biology associated with the colonization of the terrestrial environment. Finally, we explore the significance of A. agrestis in crop improvement and synthetic biology applications in general.

Reducing Pre- and Post-Treatments in Cryopreservation Protocol and Testing Storage at -80 °C for Norway Spruce Embryogenic Cultures

Somatic embryogenesis (SE) is considered the most effective method for vegetative propagation of Norway spruce (Picea abies L. Karst). For mass propagation, a storage method that is able to handle large quantities of embryogenic tissues (ETs) reliably and at a low cost is required. The aim of the present study was to compare freezing at -80 °C in a freezer to cryopreservation using liquid nitrogen (LN) as a method for storing Norway spruce ETs. The possibility of simplifying both the pre-treatment and thawing processes in cryopreservation was also studied. The addition of abscisic acid (ABA) to the pre-treatment media and using polyethylene glycol PEG4000 instead of PEG6000 in a cryoprotectant solution were tested. Both the pre-and post-treatments on semi-solid media could be simplified by reducing the number of media, without any loss of genotype or embryo production capacity of ETs. On the contrary, the storage of ETs in a freezer at -80 °C instead of using LN was not possible, and the addition of ABA to the pre-treatment media did not provide benefits but increased costs. The lower regeneration rate after using PEG4000 instead of PEG6000 in a cryoprotectant solution in cryovials was unexpected and unwanted. The simplified pre-and post-treatment protocol will remarkably reduce the workload and costs in the mass-cryopreservation of future forest regeneration materials and in thawing the samples for mass propagations, respectively.

Genetic basis of growth reaction to drought stress differs in contrasting high-latitude treeline ecotones of a widespread conifer

Climate change is increasing the frequency and intensity of drought events in many boreal forests. Trees are sessile organisms with a long generation time, which makes them vulnerable to fast climate change and hinders fast adaptations. Therefore, it is important to know how forests cope with drought stress and to explore the genetic basis of these reactions. We investigated three natural populations of white spruce (Picea glauca) in Alaska, located at one drought-limited and two cold-limited treelines with a paired plot design of one forest and one treeline plot. We obtained individual increment cores from 458 trees and climate data to assess dendrophenotypes, in particular the growth reaction to drought stress. To explore the genetic basis of these dendrophenotypes, we genotyped the individual trees at 3000 single nucleotide polymorphisms in candidate genes and performed genotype-phenotype association analysis using linear mixed models and Bayesian sparse linear mixed models. Growth reaction to drought stress differed in contrasting treeline populations. Therefore, the populations are likely to be unevenly affected by climate change. We identified 40 genes associated with dendrophenotypic traits that differed among the treeline populations. Most genes were identified in the drought-limited site, indicating comparatively strong selection pressure of drought-tolerant phenotypes. Contrasting patterns of drought-associated genes among sampled sites and in comparison to Canadian populations in a previous study suggest that drought adaptation acts on a local scale. Our results highlight genes that are associated with wood traits which in turn are critical for the establishment and persistence of future forests under climate change.

Exogenous Glutathione Promotes the Proliferation of Pinus koraiensis Embryonic Cells and the Synthesis of Glutathione and Ascorbic Acid

Somatic embryogenesis (SE), which leads to the formation of embryonic callus (EC) tissue, is the most promising method for large-scale production and selective breeding of woody plants. However, in many species, SE suffers from low proliferation rates, hindering the production of improved plant materials. One way of improving proliferation rates is achieved by improving the redox status of the culture medium. In this study, we investigated the effects of exogenous glutathione (GSH) and L-buthionine sulfoximine (BSO, the inhibitor of glutathione synthase) on the EC proliferation rate in Korean pine (Pinus koraiensis), using cell lines with both high (F: 001#-001) and low (S: 001#-010) proliferation potential. We found that exogenous GSH promoted cell proliferation in both cell lines, while exogenous BSO inhibited proliferation in both cell lines. At 35 d with exogenous GSH treatment, the fresh weight of F and S cell lines increased by 35.48% and 48.39%, respectively, compared with the control. The exogenous application of GSH increased the intracellular levels of GSH, total GSH (T-GSH), oxidized glutathione (GSSG), ascorbic acid (ASA), total ASA (T-ASA), and the ratios of GSH:T-GSH and ASA:T-ASA in both F and S cell lines. Furthermore, exogenous GSH increased the activity of both glutathione reductase (GR) and dehydroascorbate reductase (DHAR) while decreasing the activity of ascorbate peroxidase (APX) in both cell lines. It appears that the application of exogenous GSH promotes a reducing cultural environment, which is conducive to EC proliferation in Korean pine. By helping to reveal the mechanism whereby GSH regulates redox homeostasis in Korean pine EC cells, we have laid the foundation for a large-scale breeding of Korean pine somatic embryogenesis technology system.

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.

Transcriptome Analysis and Intraspecific Variation in Spanish Fir (Abies pinsapo Boiss.)

Spanish fir (Abies pinsapo Boiss.) is an endemic, endangered tree that has been scarcely investigated at the molecular level. In this work, the transcriptome of Spanish fir was assembled, providing a large catalog of expressed genes (22,769), within which a high proportion were full-length transcripts (12,545). This resource is valuable for functional genomics studies and genome annotation in this relict conifer species. Two intraspecific variations of A. pinsapo can be found within its largest population at the Sierra de las Nieves National Park: one with standard green needles and another with bluish-green needles. To elucidate the causes of both phenotypes, we studied different physiological and molecular markers and transcriptome profiles in the needles. "Green" trees showed higher electron transport efficiency and enhanced levels of chlorophyll, protein, and total nitrogen in the needles. In contrast, needles from "bluish" trees exhibited higher contents of carotenoids and cellulose. These results agreed with the differential transcriptomic profiles, suggesting an imbalance in the nitrogen status of "bluish" trees. Additionally, gene expression analyses suggested that these differences could be associated with different epigenomic profiles. Taken together, the reported data provide new transcriptome resources and a better understanding of the natural variation in this tree species, which can help improve guidelines for its conservation and the implementation of adaptive management strategies under climatic change.

Comparative Genomics Analysis of Repetitive Elements in Ten Gymnosperm Species: "Dark Repeatome" and Its Abundance in Conifer and Gnetum Species

Repetitive elements (RE) and transposons (TE) can comprise up to 80% of some plant genomes and may be essential for regulating their evolution and adaptation. The “repeatome” information is often unavailable in assembled genomes because genomic areas of repeats are challenging to assemble and are often missing from final assembly. However, raw genomic sequencing data contain rich information about RE/TEs. Here, raw genomic NGS reads of 10 gymnosperm species were studied for the content and abundance patterns of their “repeatome”. We utilized a combination of alignment on databases of repetitive elements and de novo assembly of highly repetitive sequences from genomic sequencing reads to characterize and calculate the abundance of known and putative repetitive elements in the genomes of 10 conifer plants: Pinus taeda, Pinus sylvestris, Pinus sibirica, Picea glauca, Picea abies, Abies sibirica, Larix sibirica, Juniperus communis, Taxus baccata, and Gnetum gnemon. We found that genome abundances of known and newly discovered putative repeats are specific to phylogenetically close groups of species and match biological taxa. The grouping of species based on abundances of known repeats closely matches the grouping based on abundances of newly discovered putative repeats (kChains) and matches the known taxonomic relations.

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 (Piceaglauca [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.

see also the Perspective by Lars Opgenoorth and Christian Rellstab

Construction of a reference transcriptome for the analysis of male sterility in sugi (Cryptomeria japonica D. Don) focusing on MALE STERILITY 1 (MS1)

Sugi (Cryptomeria japonica D. Don) is an important conifer used for afforestation in Japan. As the genome of this species is 11 Gbps, it is too large to assemble within a short timeframe. Transcriptomics is one approach that can address this deficiency. Here we designed a workflow consisting of three stages to de novo assemble transcriptome using Oases and Trinity. The three transcriptomic stage used were independent assembly, automatic and semi-manual integration, and refinement by filtering out potential contamination. We identified a set of 49,795 cDNA and an equal number of translated proteins. According to the benchmark set by BUSCO, 87.01% of cDNAs identified were complete genes, and 78.47% were complete and single-copy genes. Compared to other full-length cDNA resources collected by Sanger and PacBio sequencers, the extent of the coverage in our dataset was the highest, indicating that these data can be safely used for further studies. When two tissue-specific libraries were compared, there were significant expression differences between male strobili and leaf and bark sets. Moreover, subtle expression difference between male-fertile and sterile libraries were detected. Orthologous genes from other model plants and conifer species were identified. We demonstrated that our transcriptome assembly output (CJ3006NRE) can serve as a reference transcriptome for future functional genomics and evolutionary biology studies.

Champions of winter survival: cold acclimation and molecular regulation of cold hardiness in evergreen conifers

Evergreen conifers are champions of winter survival, based on their remarkable ability to acclimate to cold and develop cold hardiness. Counterintuitively, autumn cold acclimation is triggered not only by exposure to low temperature, but also by a combination of decreasing temperature, decreasing photoperiod and changes in light quality. These environmental cues control a network of signaling pathways that coordinate cold acclimation and cold hardiness in overwintering conifers, leading to cessation of growth, bud dormancy, freezing tolerance and changes in energy metabolism. Advances in genomic, transcriptomic and metabolomic tools for conifers have improved our understanding of how trees sense and respond to changes in temperature and light during cold acclimation and the development of cold hardiness, but there remain considerable gaps deserving further research in conifers. In the first section of this review, we focus on the physiological mechanisms used by evergreen conifers to adjust metabolism seasonally and to protect overwintering tissues against winter stresses. In the second section, we review how perception of low temperature and photoperiod regulate the induction of cold acclimation. Finally, we explore the evolutionary context of cold acclimation in conifers and evaluate challenges imposed on them by changing climate and discuss emerging areas of research in the field.

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.

Genomic selection for resistance to spruce budworm in white spruce and relationships with growth and wood quality traits

With climate change, the pressure on tree breeding to provide varieties with improved resilience to biotic and abiotic stress is increasing. As such, pest resistance is of high priority but has been neglected in most tree breeding programs, given the complexity of phenotyping for these traits and delays to assess mature trees. In addition, the existing genetic variation of resistance and its relationship with productivity should be better understood for their consideration in multitrait breeding. In this study, we evaluated the prospects for genetic improvement of the levels of acetophenone aglycones (AAs) in white spruce needles, which have been shown to be tightly linked to resistance to spruce budworm. Furthermore, we estimated the accuracy of genomic selection (GS) for these traits, allowing selection at a very early stage to accelerate breeding. A total of 1,516 progeny trees established on five sites and belonging to 136 full-sib families from a mature breeding population in New Brunswick were measured for height growth and genotyped for 4,148 high-quality SNPs belonging to as many genes along the white spruce genome. In addition, 598 trees were assessed for levels of AAs piceol and pungenol in needles, and 578 for wood stiffness. GS models were developed with the phenotyped trees and then applied to predict the trait values of unphenotyped trees. AAs were under moderate-to-high genetic control (h 2: 0.43-0.57) with null or marginally negative genetic correlations with other traits. The prediction accuracy of GS models (GBLUP) for AAs was high (PAAC: 0.63-0.67) and comparable or slightly higher than pedigree-based (ABLUP) or BayesCπ models. We show that AA traits can be improved and that GS speeds up the selection of improved trees for insect resistance and for growth and wood quality traits. Various selection strategies were tested to optimize multitrait gains.

Whole-exome sequencing reveals a long-term decline in effective population size of red spruce (Picea rubens)

Understanding the factors influencing the current distribution of genetic diversity across a species range is one of the main questions of evolutionary biology, especially given the increasing threat to biodiversity posed by climate change. Historical demographic processes such as population expansion or bottlenecks and decline are known to exert a predominant influence on past and current levels of genetic diversity, and revealing this demo-genetic history can have immediate conservation implications. We used a whole-exome capture sequencing approach to analyze polymorphism across the gene space of red spruce (Picea rubens Sarg.), an endemic and emblematic tree species of eastern North America high-elevation forests that are facing the combined threat of global warming and increasing human activities. We sampled a total of 340 individuals, including populations from the current core of the range in northeastern USA and southeastern Canada and from the southern portions of its range along the Appalachian Mountains, where populations occur as highly fragmented mountaintop "sky islands." Exome capture baits were designed from the closely relative white spruce (P. glauca Voss) transcriptome, and sequencing successfully captured most regions on or near our target genes, resulting in the generation of a new and expansive genomic resource for studying standing genetic variation in red spruce applicable to its conservation. Our results, based on over 2 million exome-derived variants, indicate that red spruce is structured into three distinct ancestry groups that occupy different geographic regions of its highly fragmented range. Moreover, these groups show small Ne , with a temporal history of sustained population decline that has been ongoing for thousands (or even hundreds of thousands) of years. These results demonstrate the broad potential of genomic studies for revealing details of the demographic history that can inform management and conservation efforts of nonmodel species with active restoration programs, such as red spruce.

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.

Plastome Phylogenomic and Biogeographical Study on Thuja (Cupressaceae)

Investigating the biogeographical disjunction of East Asian and North American flora is key to understanding the formation and dynamics of biodiversity in the Northern Hemisphere. The small Cupressaceae genus Thuja, comprising five species, exhibits a typical disjunct distribution in East Asia and North America. Owing to obscure relationships, the biogeographical history of the genus remains controversial. Here, complete plastomes were employed to investigate the plastome evolution, phylogenetic relationships, and biogeographic history of Thuja. All plastomes of Thuja share the same gene content arranged in the same order. The loss of an IR was evident in all Thuja plastomes, and the B-arrangement as previously recognized was detected. Phylogenomic analyses resolved two sister pairs, T. standishii-T. koraiensis and T. occidentalis-T. sutchuenensis, with T. plicata sister to T. occidentalis-T. sutchuenensis. Molecular dating and biogeographic results suggest the diversification of Thuja occurred in the Middle Miocene, and the ancestral area of extant species was located in northern East Asia. Incorporating the fossil record, we inferred that Thuja likely originated from the high-latitude areas of North America in the Paleocene with a second diversification center in northern East Asia. The current geographical distribution of Thuja was likely shaped by dispersal events attributed to the Bering Land Bridge in the Miocene and subsequent vicariance events accompanying climate cooling. The potential effect of extinction may have profound influence on the biogeographical history of Thuja.

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.

An Axiom SNP genotyping array for Douglas-fir

BACKGROUND

In forest trees, genetic markers have been used to understand the genetic architecture of natural populations, identify quantitative trait loci, infer gene function, and enhance tree breeding. Recently, new, efficient technologies for genotyping thousands to millions of single nucleotide polymorphisms (SNPs) have finally made large-scale use of genetic markers widely available. These methods will be exceedingly valuable for improving tree breeding and understanding the ecological genetics of Douglas-fir, one of the most economically and ecologically important trees in the world.

RESULTS

We designed SNP assays for 55,766 potential SNPs that were discovered from previous transcriptome sequencing projects. We tested the array on ~ 2300 related and unrelated coastal Douglas-fir trees (Pseudotsuga menziesii var. menziesii) from Oregon and Washington, and 13 trees of interior Douglas-fir (P. menziesii var. glauca). As many as ~ 28 K SNPs were reliably genotyped and polymorphic, depending on the selected SNP call rate. To increase the number of SNPs and improve genome coverage, we developed protocols to 'rescue' SNPs that did not pass the default Affymetrix quality control criteria (e.g., 97% SNP call rate). Lowering the SNP call rate threshold from 97 to 60% increased the number of successful SNPs from 20,669 to 28,094. We used a subset of 395 unrelated trees to calculate SNP population genetic statistics for coastal Douglas-fir. Over a range of call rate thresholds (97 to 60%), the median call rate for SNPs in Hardy-Weinberg equilibrium ranged from 99.2 to 99.7%, and the median minor allele frequency ranged from 0.198 to 0.233. The successful SNPs also worked well on interior Douglas-fir.

CONCLUSIONS

Based on the original transcriptome assemblies and comparisons to version 1.0 of the Douglas-fir reference genome, we conclude that these SNPs can be used to genotype about 10 K to 15 K loci. The Axiom genotyping array will serve as an excellent foundation for studying the population genomics of Douglas-fir and for implementing genomic selection. We are currently using the array to construct a linkage map and test genomic selection in a three-generation breeding program for coastal Douglas-fir.

Functional and morphological evolution in gymnosperms: A portrait of implicated gene families

Gymnosperms diverged from their sister plant clade of flowering plants 300 Mya. Morphological and functional divergence between the two major seed plant clades involved significant changes in their reproductive biology, water-conducting systems, secondary metabolism, stress defense mechanisms, and small RNA-mediated epigenetic silencing. The relatively recent sequencing of several gymnosperm genomes and the development of new genomic resources have enabled whole-genome comparisons within gymnosperms, and between angiosperms and gymnosperms. In this paper, we aim to understand how genes and gene families have contributed to the major functional and morphological differences in gymnosperms, and how this information can be used for applied breeding and biotechnology. In addition, we have analyzed the angiosperm versus gymnosperm evolution of the pleiotropic drug resistance (PDR) gene family with a wide range of functionalities in plants' interaction with their environment including defense mechanisms. Some of the genes reviewed here are newly studied members of gene families that hold potential for biotechnological applications related to commercial and pharmacological value. Some members of conifer gene families can also be exploited for their potential in phytoremediation applications.

The large repertoire of conifer NLR resistance genes includes drought responsive and highly diversified RNLs

The NLRs or NBS-LRRs (nucleotide-binding, leucine-rich-repeat) form the largest resistance gene family in plants, with lineage-specific contingents of TNL, CNL and RNL subfamilies and a central role in resilience to stress. The origin, evolution and distribution of NLR sequences has been unclear owing in part to the variable size and diversity of the RNL subfamily and a lack of data in Gymnosperms. We developed, searched and annotated transcriptomes assemblies of seven conifers and identified a resource of 3816 expressed NLR sequences. Our analyses encompassed sequences data spanning the major groups of land plants and determinations of NLR transcripts levels in response to drought in white spruce. We showed that conifers have among the most diverse and numerous RNLs in tested land plants. We report an evolutionary swap in the formation of RNLs, which emerged from the fusion of an RPW8 domain to a NB-ARC domain of CNL. We uncovered a quantitative relationship between RNLs and TNLs across all land plants investigated, with an average ratio of 1:10. The conifer RNL repertoire harbours four distinct groups, with two that differ from Angiosperms, one of which contained several upregulated sequences in response to drought while the majority of responsive NLRs are downregulated.

Functional Analysis of the PgCesA3 White Spruce Cellulose Synthase Gene Promoter in Secondary Xylem

Cellulose is an essential structural component of the plant cell wall. Its biosynthesis involves genes encoding cellulose synthase enzymes and a complex transcriptional regulatory network. Three cellulose synthases have been identified in conifers as being potentially involved in secondary cell wall biosynthesis because of their preferential expression in xylem tissues; however, no direct functional association has been made to date. In the present work, we characterized the white spruce [Picea glauca (Moench) Voss] cellulose synthase PgCesA3 gene and 5' regulatory elements. Phylogenetic analysis showed that PgCesA1-3 genes grouped with secondary cell wall-associated Arabidopsis cellulose synthase genes, such as AtCesA8, AtCesA4, and AtCesA7. We produced transgenic spruce expressing the GUS reporter gene driven by the PgCesA3 promoter. We observed blue staining in differentiating xylem cells from stem and roots, and in foliar guard cells indicating that PgCesA3 is clearly involved in secondary cell wall biosynthesis. The promoter region sequence of PgCesA3 contained several putative MYB cis-regulatory elements including AC-I like motifs and secondary wall MYB-responsive element (SMRE); however, it lacked SMRE4, 7 and 8 that correspond to the sequences of AC-I, II, and III. Based on these findings and results of previous transient trans-activation assays that identified interactions between the PgCesA3 promoter and different MYB transcription factors, we performed electrophoretic mobility shift assays with MYB recombinant proteins and cis-regulatory elements present in the PgCesA3 promoter. We found that PgMYB12 bound to a canonical AC-I element identified in the Pinus taeda PAL promoter and two AC-I like elements. We hypothesized that the PgMYB12 could regulate PgCesA3 in roots based on previous expression results. This functional study of PgCesA3 sequences and promoter opens the door for future studies on the interaction between PgMYBs and the PgCesA3 regulatory elements.

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.

Association genetics of acetophenone defence against spruce budworm in mature white spruce

BACKGROUND

Outbreaks of spruce budworm (SBW, Choristoneura fumiferana Clem.) cause major recurrent damage in boreal conifers such as white spruce (Picea glauca [Moench] Voss) and large losses of forest biomass in North America. Although defensive phenolic compounds have recently been linked to chemical resistance against SBW, their genetic basis remains poorly understood in forest trees, especially in conifers. Here, we used diverse association genetics approaches to discover genes and their variants that may control the accumulation of acetophenones, and dissect the genetic architecture of these defence compounds against SBW in white spruce mature trees.

RESULTS

Out of 4747 single nucleotide polymorphisms (SNPs) from 2312 genes genotyped in a population of 211 unrelated individuals, genetic association analyses identified 35 SNPs in 33 different genes that were significantly associated with the defence traits by using single-locus, multi-locus and multi-trait approaches. The multi-locus approach was particularly effective at detecting SNP-trait associations that explained a large fraction of the phenotypic variance (from 20 to 43%). Significant genes were regulatory including the NAC transcription factor, or they were involved in carbohydrate metabolism, falling into the binding, catalytic or transporter activity functional classes. Most of them were highly expressed in foliage. Weak positive phenotypic correlations were observed between defence and growth traits, indicating little or no evidence of defence-growth trade-offs.

CONCLUSIONS

This study provides new insights on the genetic architecture of tree defence traits, contributing to our understanding of the physiology of resistance mechanisms to biotic factors and providing a basis for the genetic improvement of the constitutive defence of white spruce against SBW.

Tree rings provide a new class of phenotypes for genetic associations that foster insights into adaptation of conifers to climate change

Local adaptation in tree species has been documented through a long history of common garden experiments where functional traits (height, bud phenology) are used as proxies for fitness. However, the ability to identify genes or genomic regions related to adaptation to climate requires the evaluation of traits that precisely reflect how and when climate exerts selective constraints. We combine dendroecology with association genetics to establish a link between genotypes, phenotypes and interannual climatic fluctuations. We illustrate this approach by examining individual tree responses embedded in the annual rings of 233 Pinus strobus trees growing in a common garden experiment representing 38 populations from the majority of its range. We found that interannual variability in growth was affected by low temperatures during spring and autumn, and by summer heat and drought. Among-population variation in climatic sensitivity was significantly correlated with the mean annual temperature of the provenance, suggesting local adaptation. Genotype-phenotype associations using these new tree-ring phenotypes validated nine candidate genes identified in a previous genetic-environment association study. Combining dendroecology with association genetics allowed us to assess tree vulnerability to past climate at fine temporal scales and provides avenues for future genomic studies on functional adaptation in forest trees.

Expansion of the dehydrin gene family in the Pinaceae is associated with considerable structural diversity and drought-responsive expression

Temperatures are expected to increase over the next century in all terrestrial biomes and particularly in boreal forests, where drought-induced mortality has been predicted to rise. Genomics research is helping to develop hypotheses regarding the molecular basis of drought tolerance and recent work proposed that the osmo-protecting dehydrin proteins have undergone a clade-specific expansion in the Pinaceae, a major group of conifer trees. The objectives of this study were to identify all of the putative members of the gene family, trace their evolutionary origin, examine their structural diversity and test for drought-responsive expression. We identified 41 complete dehydrin coding sequences in Picea glauca, which is four times more than most angiosperms studied to date, and more than in pines. Phylogenetic reconstructions indicated that the family has undergone an expansion in conifers, with parallel evolution implicating the sporadic resurgence of certain amino acid sequence motifs, and a major duplication giving rise to a clade specific to the Pinaceae. A variety of plant dehydrin structures were identified with variable numbers of the A-, E-, S- and K-segments and an N-terminal (N1) amino acid motif including assemblages specific to conifers. The expression of several of the spruce dehydrins was tissue preferential under non-stressful conditions or responded to water stress after 7-18 days without watering, reflecting changes in osmotic potential. We found that dehydrins with N1 K2 and N1 AESK2 sequences were the most responsive to the lack of water. Together, the family expansion, drought-responsive expression and structural diversification involving loss and gain of amino acid motifs suggests that subfunctionalization has driven the diversification seen among dehydrin gene duplicates. Our findings clearly indicate that dehydrins represent a large family of candidate genes for drought tolerance in spruces and in other Pinaceae that may underpin adaptability in spatially and temporally variable environments.

Single-Copy Genes as Molecular Markers for Phylogenomic Studies in Seed Plants

Phylogenetic relationships among seed plant taxa, especially within the gymnosperms, remain contested. In contrast to angiosperms, for which several genomic, transcriptomic and phylogenetic resources are available, there are few, if any, molecular markers that allow broad comparisons among gymnosperm species. With few gymnosperm genomes available, recently obtained transcriptomes in gymnosperms are a great addition to identifying single-copy gene families as molecular markers for phylogenomic analysis in seed plants. Taking advantage of an increasing number of available genomes and transcriptomes, we identified single-copy genes in a broad collection of seed plants and used these to infer phylogenetic relationships between major seed plant taxa. This study aims at extending the current phylogenetic toolkit for seed plants, assessing its ability for resolving seed plant phylogeny, and discussing potential factors affecting phylogenetic reconstruction. In total, we identified 3,072 single-copy genes in 31 gymnosperms and 2,156 single-copy genes in 34 angiosperms. All studied seed plants shared 1,469 single-copy genes, which are generally involved in functions like DNA metabolism, cell cycle, and photosynthesis. A selected set of 106 single-copy genes provided good resolution for the seed plant phylogeny except for gnetophytes. Although some of our analyses support a sister relationship between gnetophytes and other gymnosperms, phylogenetic trees from concatenated alignments without 3rd codon positions and amino acid alignments under the CAT + GTR model, support gnetophytes as a sister group to Pinaceae. Our phylogenomic analyses demonstrate that, in general, single-copy genes can uncover both recent and deep divergences of seed plant phylogeny.

Gene expression analysis of primordial shoot explants collected from mature white spruce (Picea glauca) trees that differ in their responsiveness to somatic embryogenesis induction

Within a plantation of clonal somatic embryo-derived white spruce trees that belonged to four genotypes, one genotype (G6) has consistently responded for the last 16 years, to the induction of somatic embryogenesis within primordial shoot explants. Analysis of fourteen individuals within this genotype subsequently revealed a group of clonal trees that were nonresponsive. This in turn provided a unique opportunity to conduct differential gene expression analysis in the absence of genotype-specific factors. Absolute qPCR was first used to expand the analysis of several genes previously identified via microarray analysis to be differentially expressed during SE induction, along with the inclusion of two nonresponsive genotypes. While this demonstrated a high level of repeatability within, and between, responsive and nonresponsive genotypes, it did not support our previous contention that an adaptive stress response plays a role in SE induction responsiveness, at least with respect to the candidate genes we analyzed. RNAseq analysis was then used to compare responsive and nonresponsive G6 primordial shoots during the somatic embryogenesis induction treatment. Although not analyzed in this study, this included samples of callus and embryonal masses previously generated from G6 explants. In addition to revealing a large number of differentially expressed genes, de novo assembly of unmapped reads was used to generate over 25,000 contigs that potentially represent previously unidentified transcripts. This included a MADS-domain gene that was found to be the most highly differentially expressed gene within responsive shoot explants during the first seven days of the induction treatment.

NorWood: a gene expression resource for evo-devo studies of conifer wood development

The secondary xylem of conifers is composed mainly of tracheids that differ anatomically and chemically from angiosperm xylem cells. There is currently no high-spatial-resolution data available profiling gene expression during wood formation for any coniferous species, which limits insight into tracheid development. RNA-sequencing data from replicated, high-spatial-resolution section series throughout the cambial and woody tissues of Picea abies were used to generate the NorWood.conGenIE.org web resource, which facilitates exploration of the associated gene expression profiles and co-expression networks. Integration within PlantGenIE.org enabled a comparative regulomics analysis, revealing divergent co-expression networks between P. abies and the two angiosperm species Arabidopsis thaliana and Populus tremula for the secondary cell wall (SCW) master regulator NAC Class IIB transcription factors. The SCW cellulose synthase genes (CesAs) were located in the neighbourhoods of the NAC factors in A. thaliana and P. tremula, but not in P. abies. The NorWood co-expression network enabled identification of potential SCW CesA regulators in P. abies. The NorWood web resource represents a powerful community tool for generating evo-devo insights into the divergence of wood formation between angiosperms and gymnosperms and for advancing understanding of the regulation of wood development in P. abies.

Gene copy number variations in adaptive evolution: The genomic distribution of gene copy number variations revealed by genetic mapping and their adaptive role in an undomesticated species, white spruce (Picea glauca)

Gene copy number variation (CNV) has been associated with phenotypic variability in animals and plants, but a genomewide understanding of their impacts on phenotypes is largely restricted to human and agricultural systems. As such, CNVs have rarely been considered in investigations of the genomic architecture of adaptation in wild species. Here, we report on the genetic mapping of gene CNVs in white spruce, which lacks a contiguous assembly of its large genome (~20 Gb), and their relationships with adaptive phenotypic variation. We detected 3,911 gene CNVs including de novo structural variations using comparative genome hybridization on arrays (aCGH) in a large progeny set. We inferred the heterozygosity at CNV loci within parents by comparing haploid and diploid tissues and genetically mapped 82 gene CNVs. Our analysis showed that CNVs were distributed over 10 linkage groups and identified four CNV hotspots that we predict to occur in other species of the Pinaceae. Significant relationships were found between 29 of the gene CNVs and adaptive traits based on regression analyses with timings of bud set and bud flush, and height growth, suggesting a role for CNVs in climate adaptation. The importance of CNVs in adaptive evolution of white spruce was also indicated by functional gene annotations and the clustering of 31% of the mapped adaptive gene CNVs in CNV hotspots. Taken together, these results illustrate the feasibility of studying CNVs in undomesticated species and represent a major step towards a better understanding of the roles of CNVs in adaptive evolution.

The gene expression landscape of pine seedling tissues

Conifers dominate vast regions of the Northern hemisphere. They are the main source of raw materials for timber industry as well as a wide range of biomaterials. Despite their inherent difficulties as experimental models for classical plant biology research, the technological advances in genomics research are enabling fundamental studies on these plants. The use of laser capture microdissection followed by transcriptomic analysis is a powerful tool for unravelling the molecular and functional organization of conifer tissues and specialized cells. In the present work, 14 different tissues from 1-month-old maritime pine (Pinus pinaster) seedlings have been isolated and their transcriptomes analysed. The results increased the sequence information and number of full-length transcripts from a previous reference transcriptome and added 39 841 new transcripts. In total, 2376 transcripts were ubiquitously expressed in all of the examined tissues. These transcripts could be considered the core 'housekeeping genes' in pine. The genes have been clustered in function to their expression profiles. This analysis reduced the number of profiles to 38, most of these defined by their expression in a unique tissue that is much higher than in the other tissues. The expression and localization data are accessible at ConGenIE.org (http://v22.popgenie.org/microdisection/). This study presents an overview of the gene expression distribution in different pine tissues, specifically highlighting the relationships between tissue gene expression and function. This transcriptome atlas is a valuable resource for functional genomics research in conifers.

Development of a Traceability System Based on a SNP Array for Large-Scale Production of High-Value White Spruce (Picea glauca)

Biological material is at the forefront of research programs, as well as application fields such as breeding, aquaculture, and reforestation. While sophisticated techniques are used to produce this material, all too often, there is no strict monitoring during the "production" process to ensure that the specific varieties are the expected ones. Confidence rather than evidence is often applied when the time comes to start a new experiment or to deploy selected varieties in the field. During the last decade, genomics research has led to the development of important resources, which have created opportunities for easily developing tools to assess the conformity of the material along the production chains. In this study, we present a simple methodology that enables the development of a traceability system which, is in fact a by-product of previous genomic projects. The plant production system in white spruce (Picea glauca) is used to illustrate our purpose. In Quebec, one of the favored strategies to produce elite varieties is to use somatic embryogenesis (SE). In order to detect human errors both upstream and downstream of the white spruce production process, this project had two main objectives: (i) to develop methods that make it possible to trace the origin of plants produced, and (ii) to generate a unique genetic fingerprint that could be used to differentiate each embryogenic cell line and ensure its genetic monitoring. Such a system had to rely on a minimum number of low-cost DNA markers and be easy to use by non-specialists. An efficient marker selection process was operationalized by testing different classification methods on simulated datasets. These datasets were generated using in-house bioinformatics tools that simulated crosses involved in the breeding program for which genotypes from hundreds of SNP markers were already available. The rate of misidentification was estimated and various sources of mishandling or contamination were identified. The method can easily be applied to other production systems for which genomic resources are already available.

Factors affecting the accuracy of genomic selection for growth and wood quality traits in an advanced-breeding population of black spruce (Picea mariana)

Background

Genomic selection (GS) uses information from genomic signatures consisting of thousands of genetic markers to predict complex traits. As such, GS represents a promising approach to accelerate tree breeding, which is especially relevant for the genetic improvement of boreal conifers characterized by long breeding cycles. In the present study, we tested GS in an advanced-breeding population of the boreal black spruce (Picea mariana [Mill.] BSP) for growth and wood quality traits, and concurrently examined factors affecting GS model accuracy.

Results

The study relied on 734 25-year-old trees belonging to 34 full-sib families derived from 27 parents and that were established on two contrasting sites. Genomic profiles were obtained from 4993 Single Nucleotide Polymorphisms (SNPs) representative of as many gene loci distributed among the 12 linkage groups common to spruce. GS models were obtained for four growth and wood traits. Validation using independent sets of trees showed that GS model accuracy was high, related to trait heritability and equivalent to that of conventional pedigree-based models. In forward selection, gains per unit of time were three times higher with the GS approach than with conventional selection. In addition, models were also accurate across sites, indicating little genotype-by-environment interaction in the area investigated. Using information from half-sibs instead of full-sibs led to a significant reduction in model accuracy, indicating that the inclusion of relatedness in the model contributed to its higher accuracies. About 500 to 1000 markers were sufficient to obtain GS model accuracy almost equivalent to that obtained with all markers, whether they were well spread across the genome or from a single linkage group, further confirming the implication of relatedness and potential long-range linkage disequilibrium (LD) in the high accuracy estimates obtained. Only slightly higher model accuracy was obtained when using marker subsets that were identified to carry large effects, indicating a minor role for short-range LD in this population.

Conclusions

This study supports the integration of GS models in advanced-generation tree breeding programs, given that high genomic prediction accuracy was obtained with a relatively small number of markers due to high relatedness and family structure in the population. In boreal spruce breeding programs and similar ones with long breeding cycles, much larger gain per unit of time can be obtained from genomic selection at an early age than by the conventional approach. GS thus appears highly profitable, especially in the context of forward selection in species which are amenable to mass vegetative propagation of selected stock, such as spruces.

A high-resolution reference genetic map positioning 8.8 K genes for the conifer white spruce: structural genomics implications and correspondence with physical distance

Over the last decade, extensive genetic and genomic resources have been developed for the conifer white spruce (Picea glauca, Pinaceae), which has one of the largest plant genomes (20 Gbp). Draft genome sequences of white spruce and other conifers have recently been produced, but dense genetic maps are needed to comprehend genome macrostructure, delineate regions involved in quantitative traits, complement functional genomic investigations, and assist the assembly of fragmented genomic sequences. A greatly expanded P. glauca composite linkage map was generated from a set of 1976 full-sib progeny, with the positioning of 8793 expressed genes. Regions with significant low or high gene density were identified. Gene family members tended to be mapped on the same chromosomes, with tandemly arrayed genes significantly biased towards specific functional classes. The map was integrated with transcriptome data surveyed across eight tissues. In total, 69 clusters of co-expressed and co-localising genes were identified. A high level of synteny was found with pine genetic maps, which should facilitate the transfer of structural information in the Pinaceae. Although the current white spruce genome sequence remains highly fragmented, dozens of scaffolds encompassing more than one mapped gene were identified. From these, the relationship between genetic and physical distances was examined and the genome-wide recombination rate was found to be much smaller than most estimates reported for angiosperm genomes. This gene linkage map shall assist the large-scale assembly of the next-generation white spruce genome sequence and provide a reference resource for the conifer genomics community.

Single-Step BLUP with Varying Genotyping Effort in Open-Pollinated Picea glauca

Maximization of genetic gain in forest tree breeding programs is contingent on the accuracy of the predicted breeding values and precision of the estimated genetic parameters. We investigated the effect of the combined use of contemporary pedigree information and genomic relatedness estimates on the accuracy of predicted breeding values and precision of estimated genetic parameters, as well as rankings of selection candidates, using single-step genomic evaluation (HBLUP). In this study, two traits with diverse heritabilities [tree height (HT) and wood density (WD)] were assessed at various levels of family genotyping efforts (0, 25, 50, 75, and 100%) from a population of white spruce (Picea glauca) consisting of 1694 trees from 214 open-pollinated families, representing 43 provenances in Québec, Canada. The results revealed that HBLUP bivariate analysis is effective in reducing the known bias in heritability estimates of open-pollinated populations, as it exposes hidden relatedness, potential pedigree errors, and inbreeding. The addition of genomic information in the analysis considerably improved the accuracy in breeding value estimates by accounting for both Mendelian sampling and historical coancestry that were not captured by the contemporary pedigree alone. Increasing family genotyping efforts were associated with continuous improvement in model fit, precision of genetic parameters, and breeding value accuracy. Yet, improvements were observed even at minimal genotyping effort, indicating that even modest genotyping effort is effective in improving genetic evaluation. The combined utilization of both pedigree and genomic information may be a cost-effective approach to increase the accuracy of breeding values in forest tree breeding programs where shallow pedigrees and large testing populations are the norm.

An Evolutionary Framework for Carpel Developmental Control Genes

Carpels are the female reproductive organs of flowering plants (angiosperms), enclose the ovules, and develop into fruits. The presence of carpels unites angiosperms, and they are suggested to be the most important autapomorphy of the angiosperms, e.g., they prevent inbreeding and allow efficient seed dispersal. Many transcriptional regulators and coregulators essential for carpel development are encoded by diverse gene families and well characterized in Arabidopsis thaliana. Among these regulators are AGAMOUS (AG), ETTIN (ETT), LEUNIG (LUG), SEUSS (SEU), SHORT INTERNODE/STYLISH (SHI/STY), and SEPALLATA1, 2, 3, 4 (SEP1, 2, 3, 4). However, the timing of the origin and their subsequent molecular evolution of these carpel developmental regulators are largely unknown. Here, we have sampled homologs of these carpel developmental regulators from the sequenced genomes of a wide taxonomic sampling of the land plants, such as Physcomitrella patens, Selaginella moellendorfii, Picea abies, and several angiosperms. Careful phylogenetic analyses were carried out that provide a phylogenetic background for the different gene families and provide minimal estimates for the ages of these developmental regulators. Our analyses and published work show that LUG-, SEU-, and SHI/STY-like genes were already present in the Most Recent Common Ancestor (MRCA) of all land plants, AG- and SEP-like genes were present in the MRCA of seed plants and their origin may coincide with the ξ Whole Genome Duplication. Our work shows that the carpel development regulatory network was, in part, recruited from preexisting network components that were present in the MRCA of angiosperms and modified to regulate gynoecium development.

CNVs into the wild: screening the genomes of conifer trees (Picea spp.) reveals fewer gene copy number variations in hybrids and links to adaptation

Background

Copy number variations (CNVs) have been linked to different phenotypes in human, including many diseases. A genome-scale understanding of CNVs is available in a few plants but none are wild species, leaving a knowledge gap regarding their genome biology and evolutionary role. We developed a reliable CNV detection method for species lacking contiguous reference genome. We selected multiple probes within 14,078 gene sequences and developed comparative genome hybridization on arrays. Gene CNVs were assessed in three full-sib families from species with 20 Gb genomes, i.e., white and black spruce, and interior spruce - a natural hybrid.

Results

We discovered hundreds of gene CNVs in each species, 3612 in total, which were enriched in functions related to stress and defense responses and narrow expression profiles, indicating a potential role in adaptation. The number of shared CNVs was in accordance with the degree of relatedness between individuals and species. The genetically mapped subset of these genes showed a wide distribution across the genome, implying numerous structural variations. The hybrid family presented significantly fewer CNVs, suggesting that the admixture of two species within one genome reduces the occurrence of CNVs.

Conclusions

The approach we developed is of particular interest in non-model species lacking a reference genome. Our findings point to a role for CNVs in adaptation. Their reduced abundance in the hybrid may limit genetic variability and evolvability of hybrids. We propose that CNVs make a qualitatively distinct contribution to adaptation which could be important for short term change.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3458-8) contains supplementary material, which is available to authorized users.

Identification of Norway Spruce MYB-bHLH-WDR Transcription Factor Complex Members Linked to Regulation of the Flavonoid Pathway

Transcription factors (TFs) forming MYB-bHLH-WDR complexes are known to regulate the biosynthesis of specialized metabolites in angiosperms through an intricate network. These specialized metabolites participate in a wide range of biological processes including plant growth, development, reproduction as well as in plant immunity. Studying the regulation of their biosynthesis is thus essential. While MYB (TFs) have been previously shown to control specialized metabolism (SM) in gymnosperms, the identity of their partners, in particular bHLH or WDR members, has not yet been revealed. To gain knowledge about MYB-bHLH-WDR transcription factor complexes in gymnosperms and their regulation of SW, we identified two bHLH homologs of AtTT8, six homologs of the MYB transcription factor AtTT2 and one WDR ortholog of AtTTG1 in Norway spruce. We investigated the expression levels of these genes in diverse tissues and upon treatments with various stimuli including methyl-salicylate, methyl-jasmonate, wounding or fungal inoculation. In addition, we also identified protein-protein interactions among different homologs of MYB, bHLH and WDR. Finally, we generated transgenic spruce cell lines overexpressing four of the Norway spruce AtTT2 homologs and observed differential regulation of genes in the flavonoid pathway and flavonoid contents.

A transcriptomic resource for Douglas-fir seed development and analysis of transcription during late megagametophyte development

Douglas-fir transcriptomics. Douglas-fir (Pseudotsuga menziesii (Mirb.) Franco) is economically important with extensive breeding programs and seed trade. However, the molecular genetics of its seed development are largely unknown. We developed a transcriptome resource covering key developmental stages of megagametophytes over time: prefertilization, fertilization, embryogenesis, and early, unfertilized abortion. RNA sequencing reads were assembled de novo into 105,505 predicted high-confidence transcripts derived from 34,521 predicted genes. Expression levels were estimated based on alignment of the original reads to the reference. Megagametophytes express a distinct set of genes compared to those of vegetative tissues. Transcripts related to signaling, protein turnover, and RNA biogenesis have lower expression values in vegetative tissues, whereas cell wall remodeling, solute transport, and seed storage protein transcripts have higher expression values in megagametophytes. Seed storage protein transcripts become very abundant in both pollinated and unpollinated megagametophytes over time, even in aborting ovules. However, the absence of protein storage bodies in unfertilized megagametophytes suggests extensive posttranscriptional mechanisms that either inhibit storage protein translation or their aggregation into protein bodies. This novel transcriptome resource provides a foundation for further important insights into conifer seed development.

The challenge of separating signatures of local adaptation from those of isolation by distance and colonization history: The case of two white pines

Accurately detecting signatures of local adaptation using genetic‐environment associations (GEAs) requires controlling for neutral patterns of population structure to reduce the risk of false positives. However, a high degree of collinearity between climatic gradients and neutral population structure can greatly reduce power, and the performance of GEA methods in such case is rarely evaluated in empirical studies. In this study, we attempted to disentangle the effects of local adaptation and isolation by environment (IBE) from those of isolation by distance (IBD) and isolation by colonization from glacial refugia (IBC) using range‐wide samples in two white pine species. For this, SNPs from 168 genes, including 52 candidate genes for growth and phenology, were genotyped in 133 and 61 populations of Pinus strobus and P. monticola, respectively. For P. strobus and using all 153 SNPs, climate (IBE) did not significantly explained among‐population variation when controlling for IBD and IBC in redundancy analyses (RDAs). However, 26 SNPs were significantly associated with climate in single‐locus GEA analyses (Bayenv2 and LFMM), suggesting that local adaptation took place in the presence of high gene flow. For P. monticola, we found no evidence of IBE using RDAs and weaker signatures of local adaptation using GEA and F_ST outlier tests, consistent with adaptation via phenotypic plasticity. In both species, the majority of the explained among‐population variation (69 to 96%) could not be partitioned between the effects of IBE, IBD, and IBC. GEA methods can account differently for this confounded variation, and this could explain the small overlap of SNPs detected between Bayenv2 and LFMM. Our study illustrates the inherent difficulty of taking into account neutral structure in natural populations and the importance of sampling designs that maximize climatic variation, while minimizing collinearity between climatic gradients and neutral structure.

Transcriptome responses to temperature, water availability and photoperiod are conserved among mature trees of two divergent Douglas-fir provenances from a coastal and an interior habitat

BACKGROUND

Local adaptation and phenotypic plasticity are important components of plant responses to variations in environmental conditions. While local adaptation has been widely studied in trees, little is known about plasticity of gene expression in adult trees in response to ever changing environmental conditions in natural habitats. Here we investigate plasticity of gene expression in needle tissue between two Douglas-fir provenances represented by 25 adult trees using deep RNA sequencing (RNA-Seq).

RESULTS

Using linear mixed models we investigated the effect of temperature, soil water availability and photoperiod on the abundance of 59189 detected transcripts. Expression of more than 80 % of all identified transcripts revealed a response to variations in environmental conditions in the field. GO term overrepresentation analysis revealed gene expression responses to temperature, soil water availability and photoperiod that are highly conserved among many plant taxa. However, expression differences between the two Douglas-fir provenances were rather small compared to the expression differences observed between individual trees. Although the effect of environment on global transcript expression was high, the observed genotype by environment (GxE) interaction of gene expression was surprisingly low, since only 21 of all detected transcripts showed a GxE interaction.

CONCLUSIONS

The majority of the transcriptome responses in plant leaf tissue is driven by variations in environmental conditions. The small variation between individuals and populations suggests strong conservation of this response within Douglas-fir. Therefore we conclude that plastic transcriptome responses to variations in environmental conditions are only weakly affected by local adaptation in Douglas-fir.

Different Alleles of a Gene Encoding Leucoanthocyanidin Reductase (PaLAR3) Influence Resistance against the Fungus Heterobasidion parviporum in Picea abies

Despite the fact that fungal diseases are a growing menace for conifers in modern silviculture, only a very limited number of molecular markers for pathogen resistance have been validated in conifer species. A previous genetic study indicated that the resistance of Norway spruce (Picea abies) to Heterobasidion annosum s.l., a pathogenic basidiomycete species complex, is linked to a quantitative trait loci that associates with differences in fungal growth in sapwood (FGS) that includes a gene, PaLAR3, which encodes a leucoanthocyanidin reductase. In this study, gene sequences showed the presence of two PaLAR3 allelic lineages in P. abies. Higher resistance was associated with the novel allele, which was found in low frequency in the four P. abies populations that we studied. Norway spruce plants carrying at least one copy of the novel allele showed a significant reduction in FGS after inoculation with Heterobasidion parviporum compared to their half-siblings carrying no copies, indicating dominance of this allele. The amount of (+) catechin, the enzymatic product of PaLAR3, was significantly higher in bark of trees homozygous for the novel allele. Although we observed that the in vitro activities of the enzymes encoded by the two alleles were similar, we could show that allele-specific transcript levels were significantly higher for the novel allele, indicating that regulation of gene expression is responsible for the observed effects in resistance, possibly caused by differences in cis-acting elements that we observe in the promoter region of the two alleles.

Dissection of expression-quantitative trait locus and allele specificity using a haploid/diploid plant system - insights into compensatory evolution of transcriptional regulation within populations

Regulation of gene expression plays a central role in translating genotypic variation into phenotypic variation. Dissection of the genetic basis of expression variation is key to understanding how expression regulation evolves. Such analyses remain challenging in contexts where organisms are outbreeding, highly heterozygous and long-lived such as in the case of conifer trees. We developed an RNA sequencing (RNA-seq)-based approach for both expression-quantitative trait locus (eQTL) mapping and the detection of cis-acting (allele-specific) vs trans-acting (non-allele-specific) eQTLs. This method can be potentially applied to many conifers. We used haploid and diploid meiotic seed tissues of a single self-fertilized white spruce (Picea glauca) individual to dissect eQTLs according to linkage and allele specificity. The genetic architecture of local eQTLs linked to the expressed genes was particularly complex, consisting of cis-acting, trans-acting and, surprisingly, compensatory cis-trans effects. These compensatory effects influence expression in opposite directions and are neutral when combined in homozygotes. Nearly half of local eQTLs were under compensation, indicating that close linkage between compensatory cis-trans factors is common in spruce. Compensated genes were overrepresented in developmental and cell organization functions. Our haploid-diploid eQTL analysis in spruce revealed that compensatory cis-trans eQTLs segregate within populations and evolve in close genetic linkage.

Genetic architecture of wood properties based on association analysis and co-expression networks in white spruce

Association studies are widely utilized to analyze complex traits but their ability to disclose genetic architectures is often limited by statistical constraints, and functional insights are usually minimal in nonmodel organisms like forest trees. We developed an approach to integrate association mapping results with co-expression networks. We tested single nucleotide polymorphisms (SNPs) in 2652 candidate genes for statistical associations with wood density, stiffness, microfibril angle and ring width in a population of 1694 white spruce trees (Picea glauca). Associations mapping identified 229-292 genes per wood trait using a statistical significance level of P < 0.05 to maximize discovery. Over-representation of genes associated for nearly all traits was found in a xylem preferential co-expression group developed in independent experiments. A xylem co-expression network was reconstructed with 180 wood associated genes and several known MYB and NAC regulators were identified as network hubs. The network revealed a link between the gene PgNAC8, wood stiffness and microfibril angle, as well as considerable within-season variation for both genetic control of wood traits and gene expression. Trait associations were distributed throughout the network suggesting complex interactions and pleiotropic effects. Our findings indicate that integration of association mapping and co-expression networks enhances our understanding of complex wood traits.

Implementation of the Realized Genomic Relationship Matrix to Open-Pollinated White Spruce Family Testing for Disentangling Additive from Nonadditive Genetic Effects

The open-pollinated (OP) family testing combines the simplest known progeny evaluation and quantitative genetics analyses as candidates’ offspring are assumed to represent independent half-sib families. The accuracy of genetic parameter estimates is often questioned as the assumption of “half-sibling” in OP families may often be violated. We compared the pedigree- vs. marker-based genetic models by analysing 22-yr height and 30-yr wood density for 214 white spruce [Picea glauca (Moench) Voss] OP families represented by 1694 individuals growing on one site in Quebec, Canada. Assuming half-sibling, the pedigree-based model was limited to estimating the additive genetic variances which, in turn, were grossly overestimated as they were confounded by very minor dominance and major additive-by-additive epistatic genetic variances. In contrast, the implemented genomic pairwise realized relationship models allowed the disentanglement of additive from all nonadditive factors through genetic variance decomposition. The marker-based models produced more realistic narrow-sense heritability estimates and, for the first time, allowed estimating the dominance and epistatic genetic variances from OP testing. In addition, the genomic models showed better prediction accuracies compared to pedigree models and were able to predict individual breeding values for new individuals from untested families, which was not possible using the pedigree-based model. Clearly, the use of marker-based relationship approach is effective in estimating the quantitative genetic parameters of complex traits even under simple and shallow pedigree structure.

Somatic Embryogenesis in Broad-Leaf Woody Plants: What We Can Learn from Proteomics

Proteomic approaches have been used to understand several regulatory aspects of plant development. Somatic embryogenesis is one of those developmental pathways that have beneficiated from the integration of proteomics data to the understanding of the molecular mechanisms that control embryogenic competence acquisition, somatic embryo development and conversion into viable plants. Nevertheless, most of the results obtained are based on the traditional model systems, very often not easily compared with the somatic embryogenesis systems of economical relevant woody species. The aim of this work is to summarize some of the applications of proteomics in the understanding of particular aspects of the somatic embryogenesis process in broad-leaf woody plants (model and non-model systems).

Conifer genomics and adaptation: at the crossroads of genetic diversity and genome function

Conifers have been understudied at the genomic level despite their worldwide ecological and economic importance but the situation is rapidly changing with the development of next generation sequencing (NGS) technologies. With NGS, genomics research has simultaneously gained in speed, magnitude and scope. In just a few years, genomes of 20-24 gigabases have been sequenced for several conifers, with several others expected in the near future. Biological insights have resulted from recent sequencing initiatives as well as genetic mapping, gene expression profiling and gene discovery research over nearly two decades. We review the knowledge arising from conifer genomics research emphasizing genome evolution and the genomic basis of adaptation, and outline emerging questions and knowledge gaps. We discuss future directions in three areas with potential inputs from NGS technologies: the evolutionary impacts of adaptation in conifers based on the adaptation-by-speciation model; the contributions of genetic variability of gene expression in adaptation; and the development of a broader understanding of genetic diversity and its impacts on genome function. These research directions promise to sustain research aimed at addressing the emerging challenges of adaptation that face conifer trees.