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

Conifers Concentrate Large Numbers of NLR Immune Receptor Genes on One Chromosome

Nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes form a major line of defense in plants, acting in both pathogen recognition and resistance machinery activation. NLRs are reported to form large gene clusters in limber pine (Pinus flexilis), but it is unknown how widespread this genomic architecture may be among the extant species of conifers (Pinophyta). We used comparative genomic analyses to assess patterns in the abundance, diversity, and genomic distribution of NLR genes. Chromosome-level whole genome assemblies and high-density linkage maps in the Pinaceae, Cupressaceae, Taxaceae, and other gymnosperms were scanned for NLR genes using existing and customized pipelines. The discovered genes were mapped across chromosomes and linkage groups and analyzed phylogenetically for evolutionary history. Conifer genomes are characterized by dense clusters of NLR genes, highly localized on one chromosome. These clusters are rich in TNL-encoding genes, which seem to have formed through multiple tandem duplication events. In contrast to angiosperms and nonconiferous gymnosperms, genomic clustering of NLR genes is ubiquitous in conifers. NLR-dense genomic regions are likely to influence a large part of the plant's resistance, informing our understanding of adaptation to biotic stress and the development of genetic resources through breeding.

Epigenetic stress memory in gymnosperms

Gymnosperms are long-lived, cone-bearing seed plants that include some of the most ancient extant plant species. These relict land plants have evolved to survive in habitats marked by chronic or episodic stress. Their ability to thrive in these environments is partly due to their phenotypic flexibility, and epigenetic regulation likely plays a crucial part in this plasticity. We review the current knowledge on abiotic and biotic stress memory in gymnosperms and the possible epigenetic mechanisms underlying long-term phenotypic adaptations. We also discuss recent technological improvements and new experimental possibilities that likely will advance our understanding of epigenetic regulation in these ancient and hard-to-study plants.

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.

Long-insert sequence capture detects high copy numbers in a defence-related beta-glucosidase gene βglu-1 with large variations in white spruce but not Norway spruce

Conifers are long-lived and slow-evolving, thus requiring effective defences against their fast-evolving insect natural enemies. The copy number variation (CNV) of two key acetophenone biosynthesis genes Ugt5/Ugt5b and βglu-1 may provide a plausible mechanism underlying the constitutively variable defence in white spruce (Picea glauca) against its primary defoliator, spruce budworm. This study develops a long-insert sequence capture probe set (Picea_hung_p1.0) for quantifying copy number of βglu-1-like, Ugt5-like genes and single-copy genes on 38 Norway spruce (Picea abies) and 40 P. glauca individuals from eight and nine provenances across Europe and North America respectively. We developed local assemblies (Piabi_c1.0 and Pigla_c.1.0), full-length transcriptomes (PIAB_v1 and PIGL_v1), and gene models to characterise the diversity of βglu-1 and Ugt5 genes. We observed very large copy numbers of βglu-1, with up to 381 copies in a single P. glauca individual. We observed among-provenance CNV of βglu-1 in P. glauca but not P. abies. Ugt5b was predominantly single-copy in both species. This study generates critical hypotheses for testing the emergence and mechanism of extreme CNV, the dosage effect on phenotype, and the varying copy number of genes with the same pathway. We demonstrate new approaches to overcome experimental challenges in genomic research in conifer defences.

Assembly and annotation of the black spruce genome provide insights on spruce phylogeny and evolution of stress response

Black spruce (Picea mariana [Mill.] B.S.P.) is a dominant conifer species in the North American boreal forest that plays important ecological and economic roles. Here, we present the first genome assembly of P. mariana with a reconstructed genome size of 18.3 Gbp and NG50 scaffold length of 36.0 kbp. A total of 66,332 protein-coding sequences were predicted in silico and annotated based on sequence homology. We analyzed the evolutionary relationships between P. mariana and 5 other spruces for which complete nuclear and organelle genome sequences were available. The phylogenetic tree estimated from mitochondrial genome sequences agrees with biogeography; specifically, P. mariana was strongly supported as a sister lineage to P. glauca and 3 other taxa found in western North America, followed by the European Picea abies. We obtained mixed topologies with weaker statistical support in phylogenetic trees estimated from nuclear and chloroplast genome sequences, indicative of ancient reticulate evolution affecting these 2 genomes. Clustering of protein-coding sequences from the 6 Picea taxa and 2 Pinus species resulted in 34,776 orthogroups, 560 of which appeared to be specific to P. mariana. Analysis of these specific orthogroups and dN/dS analysis of positive selection signatures for 497 single-copy orthogroups identified gene functions mostly related to plant development and stress response. The P. mariana genome assembly and annotation provides a valuable resource for forest genetics research and applications in this broadly distributed species, especially in relation to climate adaptation.

The genome of the Wollemi pine, a critically endangered "living fossil" unchanged since the Cretaceous, reveals extensive ancient transposon activity

We present the genome of the living fossil, Wollemia nobilis, a southern hemisphere conifer morphologically unchanged since the Cretaceous. Presumed extinct until rediscovery in 1994, the Wollemi pine is critically endangered with less than 60 wild adults threatened by intensifying bushfires in the Blue Mountains of Australia. The 12 Gb genome is among the most contiguous large plant genomes assembled, with extremely low heterozygosity and unusual abundance of DNA transposons. Reduced representation and genome re-sequencing of individuals confirms a relictual population since the last major glacial/drying period in Australia, 120 ky BP. Small RNA and methylome sequencing reveal conservation of ancient silencing mechanisms despite the presence of thousands of active and abundant transposons, including some transferred horizontally to conifers from arthropods in the Jurassic. A retrotransposon burst 8-6 my BP coincided with population decline, possibly as an adaptation enhancing epigenetic diversity. Wollemia, like other conifers, is susceptible to Phytophthora, and a suite of defense genes, similar to those in loblolly pine, are targeted for silencing by sRNAs in leaves. The genome provides insight into the earliest seed plants, while enabling conservation efforts.

Reference-free assembly of long-read transcriptome sequencing data with RNA-Bloom2

Long-read sequencing technologies have improved significantly since their emergence. Their read lengths, potentially spanning entire transcripts, is advantageous for reconstructing transcriptomes. Existing long-read transcriptome assembly methods are primarily reference-based and to date, there is little focus on reference-free transcriptome assembly. We introduce "RNA-Bloom2 [ https://github.com/bcgsc/RNA-Bloom ]", a reference-free assembly method for long-read transcriptome sequencing data. Using simulated datasets and spike-in control data, we show that the transcriptome assembly quality of RNA-Bloom2 is competitive to those of reference-based methods. Furthermore, we find that RNA-Bloom2 requires 27.0 to 80.6% of the peak memory and 3.6 to 10.8% of the total wall-clock runtime of a competing reference-free method. Finally, we showcase RNA-Bloom2 in assembling a transcriptome sample of Picea sitchensis (Sitka spruce). Since our method does not rely on a reference, it further sets the groundwork for large-scale comparative transcriptomics where high-quality draft genome assemblies are not readily available.

Comparative Genome-Wide Analysis of Two Caryopteris x Clandonensis Cultivars: Insights on the Biosynthesis of Volatile Terpenoids

Caryopteris x Clandonensis, also known as bluebeard, is an ornamental plant containing a large variety of terpenes and terpene-like compounds. Four different cultivars were subjected to a principal component analysis to elucidate variations in terpenoid-biosynthesis and consequently, two representative cultivars were sequenced on a genomic level. Functional annotation of genes as well as comparative genome analysis on long read datasets enabled the identification of cultivar-specific terpene synthase and cytochrome p450 enzyme sequences. This enables new insights, especially since terpenoids in research and industry are gaining increasing interest due to their importance in areas such as food preservation, fragrances, or as active ingredients in pharmaceutical formulations. According to BUSCO assessments, the presented genomes have an average size of 355 Mb and about 96.8% completeness. An average of 52,090 genes could be annotated as putative proteins, whereas about 42 were associated with terpene synthases and about 1340 with cytochrome p450 enzymes.

Sirex noctilio infestation led to inevitable pine death despite activating pathways involved in tolerance

Defense-related metabolome traits in pine species after infestation by Sirex noctilio are largely unknown, despite, in most cases, trees being overwhelmed. Using LC-MS-based untargeted metabolomics, we revealed the systemic metabolic changes induced by this insect in 14-year-old Pinus radiata trees, the most affected species worldwide. An immediate metabolome alteration was expressed in needles after infestation, including the up-regulation of flavonols, flavan-3-ols, oxyneolignans, auxins, proline, and tryptophan, among others. The flavan-3-ols (catechin and procyanidin B1) suggested a rapidly induced photoprotection mechanism aided by diverting proline as an alternative substrate for respiration to compensate for the progressive chlorosis that degrades photosystems. Meanwhile, glutathione, glutamate, and ascorbate levels significantly dropped in needles, which may indicate the critical oxidative stress that trees had to face since the onset of the infestation. They were not fully replenished after long-term infestation, and redox homeostasis was probably not achieved, compromising tree survival. Nevertheless, a huge auxins overexpression detected in needles throughout the infestation may reflect tolerance against the premature senescence caused by the woodwasp venom. In contrast, the metabolome of wood tissues remained initially unchanged, although it seems to collapse after three months. Overall, the metabolomics strategy adopted in this work evidenced its usefulness in uncovering the fundamental roles of plants' chemical defense that govern interactions with specific stressors.

Comparative transcriptomic responses of European and Japanese larches to infection by Phytophthora ramorum

BACKGROUND AND OBJECTIVES

Phytophthora ramorum severely affects both European larch (EL) and Japanese larch (JL) trees as indicated by high levels of mortality particularly in the UK. Field observations suggested that EL is less severely affected and so may be less susceptible to P. ramorum than JL; however, controlled inoculations have produced inconsistent or non-statistically significant differences. The present study aimed to compare RNA transcript accumulation profiles in EL and JL in response to inoculation with P. ramorum to improve our understanding of their defence responses.

METHODOLOGY

RNA-sequencing was carried out on bark tissues following the inoculation with P. ramorum of potted saplings in both EL and JL carried out under controlled environment conditions, with sampling at 1, 3, 10, and 25 days post inoculation in infected and control plants.

RESULTS

All of the inoculated trees rapidly developed lesions but no statistically significant differences were found in lesion lengths between EL and JL. RNA-Sequencing comparing control and inoculate saplings identified key differences in differentially expressed genes (DEGs) between the two larch species. European larch had rapid induction of defence genes within 24 hours of infection followed by sustained expression until 25 days after inoculation. Results in JL were more varied; upregulation was stronger but more transient and represented fewer defence pathways. Gene enrichment analyses highlighted differences in jasmonate signalling and regulation including NPR1 upregulation in EL only, and specific aspects of secondary metabolism. Some DEGs were represented by multiple responsive copies including lipoxygenase, chalcone synthase and nucleotide-binding, leucine-rich-repeat genes.

CONCLUSION

The variations between EL and JL in responsive DEGs of interest as potentially related to differences seen in the field and should be considered in the selection of trees for planting and future breeding.

Identification of Genes and Metabolic Pathways Involved in Resin Yield in Masson Pine by Integrative Analysis of Transcriptome, Proteome and Biochemical Characteristics

Masson pine (Pinus massoniana L.) is one of the most important resin-producing tree species in southern China. However, the molecular regulatory mechanisms of resin yield are still unclear in masson pine. In this study, an integrated analysis of transcriptome, proteome, and biochemical characteristics from needles of masson pine with the high and common resin yield was investigated. The results showed that chlorophyll a (Chl a), chlorophyll b (Chl b), total chlorophyll (Chl C), carotenoids (Car), glucose (Glu), gibberellin A9 (GA9), gibberellin A15 (GA15), and gibberellin A53 (GA53) were significantly increased, whereas fructose (Fru), jasmonic acid (JA), jasmonoyl-L-isoleucine (JA-ILE), gibberellin A1 (GA1), gibberellin A3 (GA3), gibberellin A19 (GA19), and gibberellin A24 (GA24) were significantly decreased in the high resin yield in comparison with those in the common one. The integrated analysis of transcriptome and proteome showed that chlorophyll synthase (chlG), hexokinase (HXK), sucrose synthase (SUS), phosphoglycerate kinase (PGK), dihydrolipoamide dehydrogenase (PDH), dihydrolipoamide succinyltransferase (DLST), 12-oxophytodienoic acid reductase (OPR), and jasmonate O-methyltransferases (JMT) were consistent at the transcriptomic, proteomic, and biochemical levels. The pathways of carbohydrate metabolism, terpenoid biosynthesis, photosynthesis, and hormone biosynthesis may play crucial roles in the regulation of resin yield, and some key genes involved in these pathways may be candidates that influence the resin yield. These results provide insights into the molecular regulatory mechanisms of resin yield and also provide candidate genes that can be applied for the molecular-assisted selection and breeding of high resin-yielding masson pine.

A genome for Cissus illustrates features underlying its evolutionary success in dry savannas

Cissus is the largest genus in Vitaceae and is mainly distributed in the tropics and subtropics. Crassulacean acid metabolism (CAM), a photosynthetic adaptation to the occurrence of succulent leaves or stems, indicates that convergent evolution occurred in response to drought stress during species radiation. Here we provide the chromosomal level assembly of Cissus rotundifolia (an endemic species in Eastern Africa) and a genome-wide comparison with grape to understand genome divergence within an ancient eudicot family. Extensive transcriptome data were produced to illustrate the genetics underpinning C. rotundifolia's ecological adaption to seasonal aridity. The modern karyotype and smaller genome of C. rotundifolia (n = 12, 350.69 Mb/1C), which lack further whole-genome duplication, were mainly derived from gross chromosomal rearrangements such as fusions and segmental duplications, and were sculpted by a very recent burst of retrotransposon activity. Bias in local gene amplification contributed to its remarkable functional divergence from grape, and the specific proliferated genes associated with abiotic and biotic responses (e.g. HSP-20, NBS-LRR) enabled C. rotundifolia to survive in a hostile environment. Reorganization of existing enzymes of CAM characterized as diurnal expression patterns of relevant genes further confer the ability to thrive in dry savannas.

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.

Evolution of complex genome architecture in gymnosperms

Gymnosperms represent an ancient lineage that diverged from early spermatophytes during the Devonian. The long fossil records and low diversity in living species prove their complex evolutionary history, which included ancient radiations and massive extinctions. Due to their ultra-large genome size, the whole-genome assembly of gymnosperms has only generated in the past 10 years and is now being further expanded into more taxonomic representations. Here, we provide an overview of the publicly available gymnosperm genome resources and discuss their assembly quality and recent findings in large genome architectures. In particular, we describe the genomic features most related to changes affecting the whole genome. We also highlight new realizations relative to repetitive sequence dynamics, paleopolyploidy, and long introns. Based on the results of relevant genomic studies of gymnosperms, we suggest additional efforts should be made toward exploring the genomes of medium-sized (5–15 gigabases) species. Lastly, more comparative analyses among high-quality assemblies are needed to understand the genomic shifts and the early species diversification of seed plants.

Genetic architecture behind developmental and seasonal control of tree growth and wood properties in Norway spruce

Genetic control of tree growth and wood formation varies depending on the age of the tree and the time of the year. Single-locus, multi-locus, and multi-trait genome-wide association studies (GWAS) were conducted on 34 growth and wood property traits in 1,303 Norway spruce individuals using exome capture to cover ~130K single-nucleotide polymorphisms (SNPs). GWAS identified associations to the different wood traits in a total of 85 gene models, and several of these were validated in a progenitor population. A multi-locus GWAS model identified more SNPs associated with the studied traits than single-locus or multivariate models. Changes in tree age and annual season influenced the genetic architecture of growth and wood properties in unique ways, manifested by non-overlapping SNP loci. In addition to completely novel candidate genes, SNPs were located in genes previously associated with wood formation, such as cellulose synthases and a NAC transcription factor, but that have not been earlier linked to seasonal or age-dependent regulation of wood properties. Interestingly, SNPs associated with the width of the year rings were identified in homologs of Arabidopsis thaliana BARELY ANY MERISTEM 1 and rice BIG GRAIN 1, which have been previously shown to control cell division and biomass production. The results provide tools for future Norway spruce breeding and functional studies.

Annotation of Siberian Larch (Larix sibirica Ledeb.) Nuclear Genome—One of the Most Cold-Resistant Tree Species in the Only Deciduous GENUS in Pinaceae

The recent release of the nuclear, chloroplast and mitochondrial genome assemblies of Siberian larch (Larix sibirica Ledeb.), one of the most cold-resistant tree species in the only deciduous genus of Pinaceae, with seasonal senescence and a rot-resistant valuable timber widely used in construction, greatly contributed to the development of genomic resources for the larch genus. Here, we present an extensive repeatome analysis and the first annotation of the draft nuclear Siberian larch genome assembly. About 66% of the larch genome consists of highly repetitive elements (REs), with the likely wave of retrotransposons insertions into the larch genome estimated to occur 4–5 MYA. In total, 39,370 gene models were predicted, with 87% of them having homology to the Arabidopsis-annotated proteins and 78% having at least one GO term assignment. The current state of the genome annotations allows for the exploration of the gymnosperm and angiosperm species for relative gene abundance in different functional categories. Comparative analysis of functional gene categories across different angiosperm and gymnosperm species finds that the Siberian larch genome has an overabundance of genes associated with programmed cell death (PCD), autophagy, stress hormone biosynthesis and regulatory pathways; genes that may play important roles in seasonal senescence and stress response to extreme cold in larch. Despite being incomplete, the draft assemblies and annotations of the conifer genomes are at a point of development where they now represent a valuable source for further genomic, genetic and population studies.

Water-related innovations in land plants evolved by different patterns of gene cooption and novelty

The origin of land plants and their descendants was marked by the evolution of key adaptations to life in terrestrial environments such as roots, vascular tissue and stomata. Though these innovations are well characterized, the evolution of the genetic toolkit underlying their development and function is poorly understood. We analysed molecular data from 532 species to investigate the evolutionary origin and diversification of genes involved in the development and regulation of these adaptations. We show that novel genes in the first land plants led to the single origin of stomata, but the stomatal closure of seed plants resulted from later gene expansions. By contrast, the major mechanism leading to the origin of vascular tissue was cooption of genes that emerged in the first land plants, enabling continuous water transport throughout the ancestral vascular plant. In turn, new key genes in the ancestors of plants with true leaves and seed plants led to the emergence of roots and lateral roots. The analysis highlights the different modes of evolution that enabled plants to conquer land, suggesting that gene expansion and cooption are the most common mechanisms of biological innovation in plant evolutionary history.

RResolver: efficient short-read repeat resolution within ABySS

BACKGROUND

De novo genome assembly is essential to modern genomics studies. As it is not biased by a reference, it is also a useful method for studying genomes with high variation, such as cancer genomes. De novo short-read assemblers commonly use de Bruijn graphs, where nodes are sequences of equal length k, also known as k-mers. Edges in this graph are established between nodes that overlap by [Formula: see text] bases, and nodes along unambiguous walks in the graph are subsequently merged. The selection of k is influenced by multiple factors, and optimizing this value results in a trade-off between graph connectivity and sequence contiguity. Ideally, multiple k sizes should be used, so lower values can provide good connectivity in lesser covered regions and higher values can increase contiguity in well-covered regions. However, current approaches that use multiple k values do not address the scalability issues inherent to the assembly of large genomes.

RESULTS

Here we present RResolver, a scalable algorithm that takes a short-read de Bruijn graph assembly with a starting k as input and uses a k value closer to that of the read length to resolve repeats. RResolver builds a Bloom filter of sequencing reads which is used to evaluate the assembly graph path support at branching points and removes paths with insufficient support. RResolver runs efficiently, taking only 26 min on average for an ABySS human assembly with 48 threads and 60 GiB memory. Across all experiments, compared to a baseline assembly, RResolver improves scaffold contiguity (NGA50) by up to 15% and reduces misassemblies by up to 12%.

CONCLUSIONS

RResolver adds a missing component to scalable de Bruijn graph genome assembly. By improving the initial and fundamental graph traversal outcome, all downstream ABySS algorithms greatly benefit by working with a more accurate and less complex representation of the genome. The RResolver code is integrated into ABySS and is available at https://github.com/bcgsc/abyss/tree/master/RResolver .

Mountain Pine Beetle Epidemic: An Interplay of Terpenoids in Host Defense and Insect Pheromones

The mountain pine beetle epidemic has highlighted the complex interactions of bark beetles with conifer host defenses. In these interactions, oleoresin terpenoids and volatiles, produced and released by the host tree, can be both harmful and beneficial to the beetle's success in colonizing a tree and completing its life cycle. The insect spends almost its entire life, from egg to adult, within the bark and phloem of a pine host, exposed to large quantities of complex mixtures of oleoresin terpenoids. Conifer oleoresin comprises mostly monoterpenes and diterpene resin acids as well as many different sesquiterpenes. It functions as a major chemical and physical defense system. However, the insect has evolved host colonization behavior and enzymes for terpenoid metabolism and detoxification that allow it to overcome some of the terpenoid defenses and, importantly, to co-opt pine monoterpenes as cues for host search and as a precursor for its own pheromone system. The insect-associated microbiome also plays a role in the metabolism of conifer terpenoids.

Identification of Metabolic Pathways Differentially Regulated in Somatic and Zygotic Embryos of Maritime Pine

Embryogenesis is a complex phase of conifer development involving hundreds of genes, and a proper understanding of this process is critical not only to produce embryos with different applied purposes but also for comparative studies with angiosperms. A global view of transcriptome dynamics during pine somatic and zygotic embryogenesis is currently missing. Here, we present a genome-wide transcriptome analysis of somatic and zygotic embryos at three developmental stages to identify conserved biological processes and gene functions during late embryogenesis. Most of the differences became more significant as the developmental process progressed from early to cotyledonary stages, and a higher number of genes were differentially expressed in somatic than in zygotic embryos. Metabolic pathways substantially affected included those involved in amino acid biosynthesis and utilization, and this difference was already observable at early developmental stages. Overall, this effect was found to be independent of the line (genotype) used to produce the somatic embryos. Additionally, transcription factors differentially expressed in somatic versus zygotic embryos were analyzed. Some potential hub regulatory genes were identified that can provide clues as to what transcription factors are controlling the process and to how the observed differences between somatic and zygotic embryogenesis in conifers could be regulated.

Wood Formation under Changing Environment: Omics Approaches to Elucidate the Mechanisms Driving the Early-to-Latewood Transition in Conifers

The global change scenarios highlight the urgency of clarifying the mechanisms driving the determination of wood traits in forest trees. Coniferous xylem is characterized by the alternation between earlywood (EW) and latewood (LW), on which proportions the wood density depend, one of the most important mechanical xylem qualities. However, the molecular mechanisms triggering the transition between the production of cells with the typical features of EW to the LW are still far from being completely elucidated. The increasing availability of omics resources for conifers, e.g., genomes and transcriptomes, would lay the basis for the comprehension of wood formation dynamics, boosting both breeding and gene-editing approaches. This review is intended to introduce the importance of wood formation dynamics and xylem traits of conifers in a changing environment. Then, an up-to-date overview of the omics resources available for conifers was reported, focusing on both genomes and transcriptomes. Later, an analysis of wood formation studies using omics approaches was conducted, with the aim of elucidating the main metabolic pathways involved in EW and LW determination. Finally, the future perspectives and the urgent needs on this research topic were highlighted.

Integrating genomic information and productivity and climate-adaptability traits into a regional white spruce breeding program

Tree improvement programs often focus on improving productivity-related traits; however, under present climate change scenarios, climate change-related (adaptive) traits should also be incorporated into such programs. Therefore, quantifying the genetic variation and correlations among productivity and adaptability traits, and the importance of genotype by environment interactions, including defense compounds involved in biotic and abiotic resistance, is essential for selecting parents for the production of resilient and sustainable forests. Here, we estimated quantitative genetic parameters for 15 growth, wood quality, drought resilience, and monoterpene traits for Picea glauca (Moench) Voss (white spruce). We sampled 1,540 trees from three open-pollinated progeny trials, genotyped with 467,224 SNP markers using genotyping-by-sequencing (GBS). We used the pedigree and SNP information to calculate, respectively, the average numerator and genomic relationship matrices, and univariate and multivariate individual-tree models to obtain estimates of (co)variance components. With few site-specific exceptions, all traits examined were under genetic control. Overall, higher heritability estimates were derived from the genomic- than their counterpart pedigree-based relationship matrix. Selection for height, generally, improved diameter and water use efficiency, but decreased wood density, microfibril angle, and drought resistance. Genome-based correlations between traits reaffirmed the pedigree-based correlations for most trait pairs. High and positive genetic correlations between sites were observed (average 0.68), except for those pairs involving the highest elevation, warmer, and moister site, specifically for growth and microfibril angle. These results illustrate the advantage of using genomic information jointly with productivity and adaptability traits, and defense compounds to enhance tree breeding selection for changing climate.

Identification and Tissue-Specific Expression Analysis of CYP720B Subfamily Genes in Slash Pine and Loblolly Pine

Diterpene resin acids (DRAs) are major components of pine oleoresin that can effectively resist the invasion of insects and pathogenic microorganisms. The subfamily of cytochrome P450s, CYP720B, catalyzes diterpene products into DRAs. Identifying CYP720B subfamily members and revealing the characteristics of tissue-specific expression would help understand diterpene-rich structures and diverse types. Slash pine and loblolly pine are important pines that provide oleoresin products. In this study, we identified CYP720B candidate genes based on the Pinus taeda V2.0 genome and full-length transcriptome of slash pine by PacBio. A total of 17 genes in slash pine and 19 in loblolly pine were identified and classified into four main clades by phylogenetic analysis. An analysis of cis-acting elements showed that CYP720B genes were closely related to adversity resistance. The gene expression of these candidates in different tissues was quantified by real-time quantitative PCR (RT–qPCR) analysis. Most of the genes showed relatively higher expression levels in roots and stems than in the other tissues, corresponding with the results of DRA component detection by gas chromatography–mass spectrometry (GC–MS), which indicated that stems and roots might be important tissues in oleoresin biosynthesis. These results provide a valuable resource for a better understanding of the biological role of individual CYP720Bs in slash pine and loblolly pine.

Overexpression of geranyl diphosphate synthase (PmGPPS1) boosts monoterpene and diterpene production involved in the response to pine wood nematode invasion

Outbreaks of pine wood nematode (PWN; Bursaphelenchus xylophilus) represent a severe biotic epidemic for the Pinus massoniana in China. When invaded by the PWN, the resistant P. massoniana might secret abundant oleoresin terpenoid to form certain defensive fronts for survival. However, the regulatory mechanisms of this process remain unclear. Here, the geranyl diphosphate synthase (PmGPPS1) gene was identified from resistant P. massoniana. Tissue-specific expression patterns of PmGPPS1 at transcript and protein level in resistant P. massoniana were determined by quantitative real-time polymerase chain reaction (qRT-PCR) and immunohistochemistry. Functional characteristics analysis of PmGPPS1 was performed on transgenic Nicotiana benthamiana by overexpression, as genetic transformation of P. massoniana is, so far, not possible. In summary, we identified and functionally characterized PmGPPS1 from the resistant P. massoniana following PWN inoculation. Tissue-specific expression patterns and localization of PmGPPS1 indicated that it may play a positive role involved in the metabolic and defensive processes of oleoresin terpenes production in response to PWN attack. Furthermore, overexpression of PmGPPS1 may enhance the production of monoterpene, among which limonene reduced the survival of PWN in vitro. In addition, PmGPPS1 upregulated the expression level of key genes involved in mevalonic acid (MVA) pathway, the methylerythritol phosphate (MEP) pathway and gibberellins (GAs) biosynthesis to boost the growth and development of tobacco through a feedback regulation mechanism. Our results offered new insights into the pivotal role of the PmGPPS1 involved in terpene-based defense mechanisms responding to the PWN invasion in resistant P. massoniana and provided a new metabolic engineering scenario to improve monoterpene production in tobacco.

Monoterpene Synthase Genes and Monoterpene Profiles in Pinus nigra subsp. laricio

In the present study, we carried out a quantitative analysis of the monoterpenes composition in different tissues of the non-model conifer Pinus nigra J.F. Arnold subsp. laricio Palib. ex Maire (P. laricio, in short). All the P. laricio tissues examined showed the presence of the same fourteen monoterpenes, among which the most abundant were β-phellandrene, α-pinene, and β-pinene, whose distribution was markedly tissue-specific. In parallel, from the same plant tissues, we isolated seven full-length cDNA transcripts coding for as many monoterpene synthases, each of which was found to be attributable to one of the seven phylogenetic groups in which the d1-clade of the canonical classification of plants' terpene synthases can be subdivided. The amino acid sequences deduced from the above cDNA transcripts allowed to predict their putative involvement in the biosynthesis of five of the monoterpenes identified. Transcripts profiling revealed a differential gene expression across the different tissues examined, and was found to be consistent with the corresponding metabolites profiles. The genomic organization of the seven isolated monoterpene synthase genes was also determined.

NLR diversity and candidate fusiform rust resistance genes in loblolly pine

Resistance to fusiform rust disease in loblolly pine (Pinus taeda) is a classic gene-for-gene system. Early resistance gene mapping in the P. taeda family 10-5 identified RAPD markers for a major fusiform rust resistance gene, Fr1. More recently, single nucleotide polymorphism (SNP) markers associated with resistance were mapped to a full-length gene model in the loblolly pine genome encoding for a nucleotide-binding site leucine-rich repeat (NLR) protein. NLR genes are one of the most abundant gene families in plant genomes and are involved in effector-triggered immunity. Inter- and intraspecies studies of NLR gene diversity and expression have resulted in improved disease resistance. To characterize NLR gene diversity and discover potential resistance genes, we assembled de novo transcriptomes from 92 loblolly genotypes from across the natural range of the species. In these transcriptomes, we identified novel NLR transcripts that are not present in the loblolly pine reference genome and found significant geographic diversity of NLR genes providing evidence of gene family evolution. We designed capture probes for these NLRs to identify and map SNPs that stably cosegregate with resistance to the SC20-21 isolate of Cronartium quercuum f.sp. fusiforme (Cqf) in half-sib progeny of the 10-5 family. We identified 10 SNPs and 2 quantitative trait loci associated with resistance to SC20-21 Cqf. The geographic diversity of NLR genes provides evidence of NLR gene family evolution in loblolly pine. The SNPs associated with rust resistance provide a resource to enhance breeding and deployment of resistant pine seedlings.

Genome-Wide Prediction of Transcription Start Sites in Conifers

The identification of promoters is an essential step in the genome annotation process, providing a framework for gene regulatory networks and their role in transcription regulation. Despite considerable advances in the high-throughput determination of transcription start sites (TSSs) and transcription factor binding sites (TFBSs), experimental methods are still time-consuming and expensive. Instead, several computational approaches have been developed to provide fast and reliable means for predicting the location of TSSs and regulatory motifs on a genome-wide scale. Numerous studies have been carried out on the regulatory elements of mammalian genomes, but plant promoters, especially in gymnosperms, have been left out of the limelight and, therefore, have been poorly investigated. The aim of this study was to enhance and expand the existing genome annotations using computational approaches for genome-wide prediction of TSSs in the four conifer species: loblolly pine, white spruce, Norway spruce, and Siberian larch. Our pipeline will be useful for TSS predictions in other genomes, especially for draft assemblies, where reliable TSS predictions are not usually available. We also explored some of the features of the nucleotide composition of the predicted promoters and compared the GC properties of conifer genes with model monocot and dicot plants. Here, we demonstrate that even incomplete genome assemblies and partial annotations can be a reliable starting point for TSS annotation. The results of the TSS prediction in four conifer species have been deposited in the Persephone genome browser, which allows smooth visualization and is optimized for large data sets. This work provides the initial basis for future experimental validation and the study of the regulatory regions to understand gene regulation in gymnosperms.

From Genome Sequencing to CRISPR-Based Genome Editing for Climate-Resilient Forest Trees

Due to the economic and ecological importance of forest trees, modern breeding and genetic manipulation of forest trees have become increasingly prevalent. The CRISPR-based technology provides a versatile, powerful, and widely accepted tool for analyzing gene function and precise genetic modification in virtually any species but remains largely unexplored in forest species. Rapidly accumulating genetic and genomic resources for forest trees enabled the identification of numerous genes and biological processes that are associated with important traits such as wood quality, drought, or pest resistance, facilitating the selection of suitable gene editing targets. Here, we introduce and discuss the latest progress, opportunities, and challenges of genome sequencing and editing for improving forest sustainability.

Assembled and annotated 26.5 Gbp coast redwood genome: a resource for estimating evolutionary adaptive potential and investigating hexaploid origin

Sequencing, assembly, and annotation of the 26.5 Gbp hexaploid genome of coast redwood (Sequoia sempervirens) was completed leading toward discovery of genes related to climate adaptation and investigation of the origin of the hexaploid genome. Deep-coverage short-read Illumina sequencing data from haploid tissue from a single seed were combined with long-read Oxford Nanopore Technologies sequencing data from diploid needle tissue to create an initial assembly, which was then scaffolded using proximity ligation data to produce a highly contiguous final assembly, SESE 2.1, with a scaffold N50 size of 44.9 Mbp. The assembly included several scaffolds that span entire chromosome arms, confirmed by the presence of telomere and centromere sequences on the ends of the scaffolds. The structural annotation produced 118,906 genes with 113 containing introns that exceed 500 Kbp in length and one reaching 2 Mb. Nearly 19 Gbp of the genome represented repetitive content with the vast majority characterized as long terminal repeats, with a 2.9:1 ratio of Copia to Gypsy elements that may aid in gene expression control. Comparison of coast redwood to other conifers revealed species-specific expansions for a plethora of abiotic and biotic stress response genes, including those involved in fungal disease resistance, detoxification, and physical injury/structural remodeling and others supporting flavonoid biosynthesis. Analysis of multiple genes that exist in triplicate in coast redwood but only once in its diploid relative, giant sequoia, supports a previous hypothesis that the hexaploidy is the result of autopolyploidy rather than any hybridizations with separate but closely related conifer species.

The Chinese pine genome and methylome unveil key features of conifer evolution

Conifers dominate the world's forest ecosystems and are the most widely planted tree species. Their giant and complex genomes present great challenges for assembling a complete reference genome for evolutionary and genomic studies. We present a 25.4-Gb chromosome-level assembly of Chinese pine (Pinus tabuliformis) and revealed that its genome size is mostly attributable to huge intergenic regions and long introns with high transposable element (TE) content. Large genes with long introns exhibited higher expressions levels. Despite a lack of recent whole-genome duplication, 91.2% of genes were duplicated through dispersed duplication, and expanded gene families are mainly related to stress responses, which may underpin conifers' adaptation, particularly in cold and/or arid conditions. The reproductive regulation network is distinct compared with angiosperms. Slow removal of TEs with high-level methylation may have contributed to genomic expansion. This study provides insights into conifer evolution and resources for advancing research on conifer adaptation and development.

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.

Complete Mitochondrial Genome of a Gymnosperm, Sitka Spruce (Picea sitchensis), Indicates a Complex Physical Structure

Plant mitochondrial genomes vary widely in size. Although many plant mitochondrial genomes have been sequenced and assembled, the vast majority are of angiosperms, and few are of gymnosperms. Most plant mitochondrial genomes are smaller than a megabase, with a few notable exceptions. We have sequenced and assembled the complete 5.5-Mb mitochondrial genome of Sitka spruce (Picea sitchensis), to date, one of the largest mitochondrial genomes of a gymnosperm. We sequenced the whole genome using Oxford Nanopore MinION, and then identified contigs of mitochondrial origin assembled from these long reads based on sequence homology to the white spruce mitochondrial genome. The assembly graph shows a multipartite genome structure, composed of one smaller 168-kb circular segment of DNA, and a larger 5.4-Mb single component with a branching structure. The assembly graph gives insight into a putative complex physical genome structure, and its branching points may represent active sites of recombination.

Breeding for Climate Change Resilience: A Case Study of Loblolly Pine (Pinus taeda L.) in North America

Earth's atmosphere is warming and the effects of climate change are becoming evident. A key observation is that both the average levels and the variability of temperature and precipitation are changing. Information and data from new technologies are developing in parallel to provide multidisciplinary opportunities to address and overcome the consequences of these changes in forest ecosystems. Changes in temperature and water availability impose multidimensional environmental constraints that trigger changes from the molecular to the forest stand level. These can represent a threat for the normal development of the tree from early seedling recruitment to adulthood both through direct mortality, and by increasing susceptibility to pathogens, insect attack, and fire damage. This review summarizes the strengths and shortcomings of previous work in the areas of genetic variation related to cold and drought stress in forest species with particular emphasis on loblolly pine (Pinus taeda L.), the most-planted tree species in North America. We describe and discuss the implementation of management and breeding strategies to increase resilience and adaptation, and discuss how new technologies in the areas of engineering and genomics are shaping the future of phenotype-genotype studies. Lessons learned from the study of species important in intensively-managed forest ecosystems may also prove to be of value in helping less-intensively managed forest ecosystems adapt to climate change, thereby increasing the sustainability and resilience of forestlands for the future.

Identification, classification, and characterization of AP2/ERF superfamily genes in Masson pine (Pinus massoniana Lamb.)

Transcription factors (TFs) play crucial regulatory roles in controlling the expression of the target genes in plants. APETALA2/Ethylene-responsive factors (AP2/ERF) are part of a large superfamily of plant-specific TFs whose members are involved in the control of plant metabolism, development and responses to various biotic and abiotic stresses. However, the AP2/ERF superfamily has not been identified systematically in Masson pine (Pinus massoniana), which is one of the most important conifer in southern China. Therefore, we performed systematic identification of the AP2/ERF superfamily using transcriptome sequencing data from Masson pine. In the current study, we obtained 88 members of the AP2/ERF superfamily. All PmAP2/ERF members could be classified into 3 main families, AP2 (7 members), RAV (7 members), ERF (73 members) families, and a soloist protein. Subcellular localization assays suggested that two members of PmAP2/ERF were nuclear proteins. Based on pine wood nematode (PWN) inoculated transcriptome and qPCR analysis, we found that many members of PmAP2/ERF could respond to PWN inoculation and PWN related treatment conditions in vitro. In general, members of the AP2/ERF superfamily play an important role in the response of Masson pine responds to PWN. Furthermore, the roles of the AP2/ERF superfamily in other physiological activities of Masson pine remain to be further studied.

Effect of elevated gibberellic acid application on growth and gene expression patterns in white spruce families from a tree improvement program in Alberta, Canada

Nine open-pollinated families of Picea glauca (Moench) Voss from the Region D1 Controlled Parentage Program (Alberta, Canada) were systematically chosen from fast, medium and slow-growth rankings based on breeding values for height from field progeny tests at age 30 years. Seeds from these families were sown and grown to age 3 years to analyze the performance and correlations of growth, physiological traits and expression of gibberellin-related genes, with and without elevated gibberellic acid 3 (GA3) application, under greenhouse conditions. We observed a significant interaction effect between families and growth groups subjected to 50 μg μl-1 of GA3 treatment, causing a decrease in apical internode length, diameter, volume and absolute transcript level for fast-growing families but an increase for families in the slow-growth group for the same traits. We also observed that in the apical internode, the gene PgGA20ox1 had significantly more relative expression under the elevated GA3 treatment than the control trees. In the stem, PgGA3ox1 showed a significantly higher relative expression under elevated GA3 treatment compared with control trees. Also, the slow-growth group showed more relative expression of PgGA20ox1 (in the apical internode) and PgGA3ox1 (in the stem) than the fast-growth group. The apical internode length and diameter significantly increased by 24% and 16%, respectively, with the hormone treatment in the slow growing group. In general, the PgGID1 and PgDELLA1 genes were upregulated and downregulated respectively, in spruce shoots under the GA3 treatment, meaning a positive feedback regulation by those genes were influencing PgGA20ox1 and PgGA3ox1 expression in that tissue type. Moreover, there was a significant correlation between absolute transcript levels of PgGA20ox1 in the apical internode and apical internode length, and absolute transcript levels of PgGA3ox1 in the stem and the diameter, in the fast-growth group families. This study shows that expression of GA genes is a limiting factor for growth in certain white spruce families with a complex feedback mechanism. Finally, absolute transcript levels of endogenous GA relative to growth parameters in juvenile seedlings could potentially be used to accelerate the early selection of families with inherently rapid apical and radial growth expansion.

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.

Jasmonate induced terpene-based defense in Pinus massoniana depresses Monochamus alternatus adult feeding

BACKGROUND

Symptoms of pine wilt disease (PWD) are caused when pathogenic pine wood nematodes (PWN) invade healthy host trees via wounds created by adult Monochamus alternatus. Exogenous methyl jasmonate (MeJA) can trigger terpene-based induced defense in conifers, which is a core part of the conifer defense system. Herein, we hypothesized that the terpene-based plant defense induced by MeJA could negatively affect the feeding behavior of M. alternatus adults, and may contribute to a new strategy in the control of PWD.

RESULTS

The feeding area for M. alternatus adults feeding on MeJA-treated seedlings was significantly smaller compared with control seedlings. However, no significant difference was detected in the enzyme activities in the midgut of beetles that had fed on these seedlings. Terpenoids were mainly accumulated in traumatic resin duct, whereas the constitutive resin duct accumulated only diterpenoids. Correspondingly, large-scaled responses at the transcriptional level mainly focused on terpenoid and phenolic biosynthesis in the defending trees.

CONCLUSIONS

Breeding tree species with a high resin yield may contribute to control of the spread of PWD by suppressing the feeding of M. alternatus adults. Transcriptome sequencing results provided abundant information for further breeding of highly resistant trees. Based on these findings, a potential push-pull strategy for the control of M. alternatus was discussed.

Impact of Environmental Factors on Stilbene Biosynthesis

Stilbenes are a small family of polyphenolic secondary metabolites that can be found in several distantly related plant species. These compounds act as phytoalexins, playing a crucial role in plant defense against phytopathogens, as well as being involved in the adaptation of plants to abiotic environmental factors. Among stilbenes, trans-resveratrol is certainly the most popular and extensively studied for its health properties. In recent years, an increasing number of stilbene compounds were subjected to investigations concerning their bioactivity. This review presents the most updated knowledge of the stilbene biosynthetic pathway, also focusing on the role of several environmental factors in eliciting stilbenes biosynthesis. The effects of ultraviolet radiation, visible light, ultrasonication, mechanical stress, salt stress, drought, temperature, ozone, and biotic stress are reviewed in the context of enhancing stilbene biosynthesis, both in planta and in plant cell and organ cultures. This knowledge may shed some light on stilbene biological roles and represents a useful tool to increase the accumulation of these valuable compounds.

Challenges and Perspectives in the Epigenetics of Climate Change-Induced Forests Decline

Forest tree species are highly vulnerable to the effects of climate change. As sessile organisms with long generation times, their adaptation to a local changing environment may rely on epigenetic modifications when allele frequencies are not able to shift fast enough. However, the current lack of knowledge on this field is remarkable, due to many challenges that researchers face when studying this issue. Huge genome sizes, absence of reference genomes and annotation, and having to analyze huge amounts of data are among these difficulties, which limit the current ability to understand how climate change drives tree species epigenetic modifications. In spite of this challenging framework, some insights on the relationships among climate change-induced stress and epigenomics are coming. Advances in DNA sequencing technologies and an increasing number of studies dealing with this topic must boost our knowledge on tree adaptive capacity to changing environmental conditions. Here, we discuss challenges and perspectives in the epigenetics of climate change-induced forests decline, aiming to provide a general overview of the state of the art.

Phytostilbenes as agrochemicals: biosynthesis, bioactivity, metabolic engineering and biotechnology

Covering: 1976 to 2020. Although constituting a limited chemical family, phytostilbenes represent an emblematic group of molecules among natural compounds. Ever since their discovery as antifungal compounds in plants and their ascribed role in human health and disease, phytostilbenes have never ceased to arouse interest for researchers, leading to a huge development of the literature in this field. Owing to this, the number of references to this class of compounds has reached the tens of thousands. The objective of this article is thus to offer an overview of the different aspects of these compounds through a large bibliography analysis of more than 500 articles. All the aspects regarding phytostilbenes will be covered including their chemistry and biochemistry, regulation of their biosynthesis, biological activities in plants, molecular engineering of stilbene pathways in plants and microbes as well as their biotechnological production by plant cell systems.

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.

Techniques for Small Non-Coding RNA Analysis in Seeds of Forest Tree Species

In recent years, the scientific community has become aware that epigenetic mechanisms play a more important role in gene regulatory networks (GRNs) than was hitherto thought, as accumulating evidence has shown that changes in epigenetics without genetic variation can affect complex traits over multiple generations. Within the epigenetic machinery, small non-coding RNAs (sRNAs, 18-24 nucleotides in length) are evolutionarily conserved RNA molecules that target mRNAs for deregulation or translational repression. They commonly have high-level regulatory functions in GRNs by mediating DNA and/or histone methylation and gene silencing essential for plant developmental programs and adaptability. Local adaptation enables plants to acquire a high fitness by, for example, properly timing developmental transitions to match plant growth stages with organism's favorable seasons. In particular, the seed represents a key evolutionary adaptation of seed plants that facilitates dispersal and reinitiates the development coupled in time with suitable environmental conditions. With the advent of high-throughput sequencing for sRNAs and computational approaches for sRNA detection and categorization, it is now feasible to unravel how sRNAs contribute to the fitness of tree species that can survive hundreds of years (e.g., conifers). Of particular interest is to disentangle the roles of sRNAs from complex genomic information in tree species with intimidating genomic sizes (commonly 20-30 Gb in conifers) and abundant nongenic components (e.g., >60% transposable elements). In this chapter, we use seeds of the conifer Picea glauca as a study system to describe the methods and protocols we used or have recently updated, from high-quality RNA isolation to sRNA identification, sequence conservation, abundance comparison, and functional analysis.

Development of a highly efficient 50K single nucleotide polymorphism genotyping array for the large and complex genome of Norway spruce (Picea abies L. Karst) by whole genome resequencing and its transferability to other spruce species

Norway spruce (Picea abies L. Karst) is one of the most important forest tree species with significant economic and ecological impact in Europe. For decades, genomic and genetic studies on Norway spruce have been challenging due to the large and repetitive genome (19.6 Gb with more than 70% being repetitive). To accelerate genomic studies, including population genetics, genome‐wide association studies (GWAS) and genomic selection (GS), in Norway spruce and related species, we here report on the design and performance of a 50K single nucleotide polymorphism (SNP) genotyping array for Norway spruce. The array is developed based on whole genome resequencing (WGS), making it the first WGS‐based SNP array in any conifer species so far. After identifying SNPs using genome resequencing data from 29 trees collected in northern Europe, we adopted a two‐step approach to design the array. First, we built a 450K screening array and used this to genotype a population of 480 trees sampled from both natural and breeding populations across the Norway spruce distribution range. These samples were then used to select high‐confidence probes that were put on the final 50K array. The SNPs selected are distributed over 45,552 scaffolds from the P. abies version 1.0 genome assembly and target 19,954 unique gene models with an even coverage of the 12 linkage groups in Norway spruce. We show that the array has a 99.5% probe specificity, >98% Mendelian allelic inheritance concordance, an average sample call rate of 96.30% and an SNP call rate of 98.90% in family trios and haploid tissues. We also observed that 23,797 probes (50%) could be identified with high confidence in three other spruce species (white spruce [Picea glauca], black spruce [P. mariana] and Sitka spruce [P. sitchensis]). The high‐quality genotyping array will be a valuable resource for genetic and genomic studies in Norway spruce as well as in other conifer species of the same genus.

Unbiased Subgenome Evolution in Allotetraploid Species of Ephedra and Its Implications for the Evolution of Large Genomes in Gymnosperms

The evolutionary dynamics of polyploid genomes and consequences of polyploidy have been studied extensively in angiosperms but very rarely in gymnosperms. The gymnospermous genus Ephedra is characterized by a high frequency of polyploidy, and thus provides an ideal system to investigate the evolutionary mode of allopolyploid genomes and test whether subgenome dominance has occurred in gymnosperms. Here, we sequenced transcriptomes of two allotetraploid species of Ephedra and their putative diploid progenitors, identified expressed homeologs, and analyzed alternative splicing and homeolog expression based on PacBio Iso-Seq and Illumina RNA-seq data. We found that the two subgenomes of the allotetraploids had similar numbers of expressed homeologs, similar percentages of homeologs with dominant expression, and approximately equal numbers of isoforms with alternative splicing, showing an unbiased subgenome evolution as in a few polyploid angiosperms, with a divergence of the two subgenomes at ∼8 Ma. In addition, the nuclear DNA content of the allotetraploid species is almost equal to the sum of two putative progenitors, suggesting limited genome restructuring after allotetraploid speciation. The allopolyploid species of Ephedra might have undergone slow diploidization, and the unbiased subgenome evolution implies that the formation of large genomes in gymnosperms could be attributed to even and slow fractionation following polyploidization.

A Reference Genome Sequence for Giant Sequoia

The giant sequoia (Sequoiadendron giganteum) of California are massive, long-lived trees that grow along the U.S. Sierra Nevada mountains. Genomic data are limited in giant sequoia and producing a reference genome sequence has been an important goal to allow marker development for restoration and management. Using deep-coverage Illumina and Oxford Nanopore sequencing, combined with Dovetail chromosome conformation capture libraries, the genome was assembled into eleven chromosome-scale scaffolds containing 8.125 Gbp of sequence. Iso-Seq transcripts, assembled from three distinct tissues, was used as evidence to annotate a total of 41,632 protein-coding genes. The genome was found to contain, distributed unevenly across all 11 chromosomes and in 63 orthogroups, over 900 complete or partial predicted NLR genes, of which 375 are supported by annotation derived from protein evidence and gene modeling. This giant sequoia reference genome sequence represents the first genome sequenced in the Cupressaceae family, and lays a foundation for using genomic tools to aid in giant sequoia conservation and management.

Genomic basis of white pine blister rust quantitative disease resistance and its relationship with qualitative resistance

The genomic architecture and molecular mechanisms controlling variation in quantitative disease resistance loci are not well understood in plant species and have been barely studied in long-generation trees. Quantitative trait loci mapping and genome-wide association studies were combined to test a large single nucleotide polymorphism (SNP) set for association with quantitative and qualitative white pine blister rust resistance in sugar pine. In the absence of a chromosome-scale reference genome, a high-density consensus linkage map was generated to obtain locations for associated SNPs. Newly discovered associations for white pine blister rust quantitative disease resistance included 453 SNPs involved in wide biological functions, including genes associated with disease resistance and others involved in morphological and developmental processes. In addition, NBS-LRR pathogen recognition genes were found to be involved in quantitative disease resistance, suggesting these newly reported genes are qualitative genes with partial resistance, they are the result of defeated qualitative resistance due to avirulent races, or they have epistatic effects on qualitative disease resistance genes. This study is a step forward in our understanding of the complex genomic architecture of quantitative disease resistance in long-generation trees, and constitutes the first step towards marker-assisted disease resistance breeding in white pine species.

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.

Allowing mutations in maximal matches boosts genome compression performance

MOTIVATION

A maximal match between two genomes is a contiguous non-extendable sub-sequence common in the two genomes. DNA bases mutate very often from the genome of one individual to another. When a mutation occurs in a maximal match, it breaks the maximal match into shorter match segments. The coding cost using these broken segments for reference-based genome compression is much higher than that of using the maximal match which is allowed to contain mutations.

RESULTS

We present memRGC, a novel reference-based genome compression algorithm that leverages mutation-containing matches (MCMs) for genome encoding. MemRGC detects maximal matches between two genomes using a coprime double-window k-mer sampling search scheme, the method then extends these matches to cover mismatches (mutations) and their neighbouring maximal matches to form long and MCMs. Experiments reveal that memRGC boosts the compression performance by an average of 27% in reference-based genome compression. MemRGC is also better than the best state-of-the-art methods on all of the benchmark datasets, sometimes better by 50%. Moreover, memRGC uses much less memory and de-compression resources, while providing comparable compression speed. These advantages are of significant benefits to genome data storage and transmission.

AVAILABILITY AND IMPLEMENTATION

https://github.com/yuansliu/memRGC.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Molecular study of drought response in the Mediterranean conifer Pinus pinaster Ait.: Differential transcriptomic profiling reveals constitutive water deficit-independent drought tolerance mechanisms

Adaptation of long-living forest trees to respond to environmental changes is essential to secure their performance under adverse conditions. Water deficit is one of the most significant stress factors determining tree growth and survival. Maritime pine (Pinus pinaster Ait.), the main source of softwood in southwestern Europe, is subjected to recurrent drought periods which, according to climate change predictions for the years to come, will progressively increase in the Mediterranean region. The mechanisms regulating pine adaptive responses to environment are still largely unknown. The aim of this work was to go a step further in understanding the molecular mechanisms underlying maritime pine response to water stress and drought tolerance at the whole plant level. A global transcriptomic profiling of roots, stems, and needles was conducted to analyze the performance of siblings showing contrasted responses to water deficit from an ad hoc designed full-sib family. Although P. pinaster is considered a recalcitrant species for vegetative propagation in adult phase, the analysis was conducted using vegetatively propagated trees exposed to two treatments: well-watered and moderate water stress. The comparative analyses led us to identify organ-specific genes, constitutively expressed as well as differentially expressed when comparing control versus water stress conditions, in drought-sensitive and drought-tolerant genotypes. Different response strategies can point out, with tolerant individuals being pre-adapted for coping with drought by constitutively expressing stress-related genes that are detected only in latter stages on sensitive individuals subjected to drought.