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.

Quantitative metabolome, proteome and transcriptome analysis of midgut and fat body tissues in the mountain pine beetle, Dendroctonus ponderosae Hopkins, and insights into pheromone biosynthesis.

Bark beetles (Coleoptera: Scolytinae) are pests of many forests around the world. The mountain pine beetle (MPB), Dendroctonus ponderosae Hopkins, is a significant pest of western North American pine forests. The MPB is able to overcome the defences of pine trees through pheromone-assisted aggregation that results in a mass attack of host trees. These pheromones, both male and female produced, are believed to be biosynthesized in the midgut and/or fat bodies of these insects. We used metabolite analysis, quantitative proteomics (iTRAQ) and transcriptomics (RNA-seq) to identify proteins and transcripts differentially expressed between sexes and between tissues when treated with juvenile hormone III. Juvenile hormone III induced frontalin biosynthesis in males and trans-verbenol biosynthesis in females, as well as affected the expression of many proteins and transcripts in sex- and tissue-specific ways. Based on these analyses, we identified candidate genes involved in the biosynthesis of frontalin, exo-brevicomin, and trans-verbenol pheromones.

Histology and cell wall biochemistry of stone cells in the physical defence of conifers against insects.

Conifers possess an array of physical and chemical defences against stem-boring insects. Stone cells provide a physical defence associated with resistance against bark beetles and weevils. In Sitka spruce (Picea sitchensis), abundance of stone cells in the cortex of apical shoots is positively correlated with resistance to white pine weevil (Pissodes strobi). We identified histological, biochemical and molecular differences in the stone cell phenotype of weevil resistant (R) or susceptible (S) Sitka spruce genotypes. R trees displayed significantly higher quantities of cortical stone cells near the apical shoot node, the primary site for weevil feeding. Lignin, cellulose, xylan and mannan were the most abundant components of stone cell secondary walls, respectively. Lignin composition of stone cells isolated from R trees contained a higher percentage of G-lignin compared with S trees. Transcript profiling revealed higher transcript abundance in the R genotype of coumarate 3-hydroxylase, a key monolignol biosynthetic gene. Developing stone cells in current year apical shoots incorporated fluorescent-tagged monolignol into the secondary cell wall, while mature stone cells of previous year apical shoots did not. Stone cell development is an ephemeral process, and fortification of shoot tips in R trees is an effective strategy against insect feeding.

Improved white spruce (Picea glauca) genome assemblies and annotation of large gene families of conifer terpenoid and phenolic defense metabolism.

White spruce (Picea glauca), a gymnosperm tree, has been established as one of the models for conifer genomics. We describe the draft genome assemblies of two white spruce genotypes, PG29 and WS77111, innovative tools for the assembly of very large genomes, and the conifer genomics resources developed in this process. The two white spruce genotypes originate from distant geographic regions of western (PG29) and eastern (WS77111) North America, and represent elite trees in two Canadian tree-breeding programs. We present an update (V3 and V4) for a previously reported PG29 V2 draft genome assembly and introduce a second white spruce genome assembly for genotype WS77111. Assemblies of the PG29 and WS77111 genomes confirm the reconstructed white spruce genome size in the 20 Gbp range, and show broad synteny. Using the PG29 V3 assembly and additional white spruce genomics and transcriptomics resources, we performed MAKER-P annotation and meticulous expert annotation of very large gene families of conifer defense metabolism, the terpene synthases and cytochrome P450s. We also comprehensively annotated the white spruce mevalonate, methylerythritol phosphate and phenylpropanoid pathways. These analyses highlighted the large extent of gene and pseudogene duplications in a conifer genome, in particular for genes of secondary (i.e. specialized) metabolism, and the potential for gain and loss of function for defense and adaptation.

Assembling the 20 Gb white spruce (Picea glauca) genome from whole-genome shotgun sequencing data.

UNLABELLED: White spruce (Picea glauca) is a dominant conifer of the boreal forests of North America, and providing genomics resources for this commercially valuable tree will help improve forest management and conservation efforts. Sequencing and assembling the large and highly repetitive spruce genome though pushes the boundaries of the current technology. Here, we describe a whole-genome shotgun sequencing strategy using two Illumina sequencing platforms and an assembly approach using the ABySS software. We report a 20.8 giga base pairs draft genome in 4.9 million scaffolds, with a scaffold N50 of 20,356 bp. We demonstrate how recent improvements in the sequencing technology, especially increasing read lengths and paired end reads from longer fragments have a major impact on the assembly contiguity. We also note that scalable bioinformatics tools are instrumental in providing rapid draft assemblies. AVAILABILITY: The Picea glauca genome sequencing and assembly data are available through NCBI (Accession#: ALWZ0100000000 PID: PRJNA83435). http://www.ncbi.nlm.nih.gov/bioproject/83435.