Association genetics of oleoresin flow in loblolly pine: discovering genes and predicting phenotype for improved resistance to bark beetles and bioenergy potential

Physiological, Proteomic, and Resin Yield-Related Genes Expression Analysis Provides Insights into the Mechanisms Regulating Resin Yield in Masson Pine

Masson pine (Pinus massoniana Lamb.) is an important resin-producing conifer species in China. Resin yield is a highly heritable trait and varies greatly among different genotypes. However, the mechanisms regulating the resin yield of masson pine remain largely unknown. In this study, physiological, proteomic, and gene expression analysis was performed on xylem tissues of masson pine with high and low resin yield. Physiological investigation showed that the activity of terpene synthase, as well as the contents of soluble sugar, jasmonic acid (JA), methyl jasmonate (MeJA), gibberellins (GA1, GA4, GA9, GA19, and GA20), indole-3-acetic acid (IAA), and abscisic acid (ABA) were significantly increased in the high yielder, whereas sucrose and salicylic acid (SA) were significantly decreased compared with the low one. A total of 2984 differentially expressed proteins (DEPs) were identified in four groups, which were mainly enriched in the biosynthesis of secondary metabolites, protein processing in the endoplasmic reticulum, carbohydrate metabolism, phytohormone biosynthesis, glutathione metabolism, and plant-pathogen interaction. Integrated physiological and proteomic analysis revealed that carbohydrate metabolism, terpenoid biosynthesis, resistance to stress, as well as JA and GA biosynthesis and signaling, play key roles in regulating resin yield. A series of proteins associated with resin yield, e.g., terpene synthase proteins (TPSs), ATP-binding cassette transporters (ABCs), glutathione S-transferase proteins (GSTs), and heat shock proteins (HSPs), were identified. Resin yield-related gene expression was also associated with resin yield. Our study unveils the implicated molecular mechanisms regulating resin yield and is of pivotal significance to breeding strategies of high resin-yielding masson pine cultivars.

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.

Evaluation of Aphid Resistance and Oleoresin Production in Indigenous Tropical Pine (Pinus merkusii Jungh. and de Vriese)

The native tropical pine (Pinus merkusii Jungh. and de Vriese) has been genetically improved in Indonesia since 1977; nevertheless, minor evaluations of aphid resistance have been conducted since 2004. As a result, a progeny test for aphid resistance was established in 2010 in Lawu, Central Java, Indonesia. Subjects in the trial were attacked significantly at the rate of 30.7% after 4 years, but surprisingly, some individuals were found to be healthy without any aphid attack. The observed a 7-year progeny trial comprised 34 families with 4 trees per unitary plot and replicated in 10 blocks. At 7 years, observations during 9 months (April–December) showed that there were differences in the range of resistance across families. The stem diameter, oleoresin production, and resistance to aphid attack were evaluated, and all traits were distinct among families except for oleoresin exudation from the western side of the stem. Five families performed above average for all three traits, while three other families had high diameter and maintained good oleoresin production. These eight families can be included in a forward selection strategy. Cluster analysis revealed that the eight best families were grouped into two of the eight clusters. Phenotypic correlations revealed that all pairs of traits were significantly related, with the highest correlation registered between stem diameter and resistance to aphid attack (0.99). Forward selection ensures the simultaneous improvement of the three traits.

The Approach in Selecting the Best Genetic Resistance against Invasive Aphid for Indigenous Tropical Pinus merkusii Jungh. et de Vriese in Indonesia

Pinus merkusii, a natural tropical pine species of Indonesia, is cultivated as the second most important artificial forest for the industry in Java, after teak, to produce oleoresin. Its genetic improvement began in 1977 because of its critical role in raising community incomes. Meanwhile, the effort for genetic improvement in aphid (Pineus boerneri) resistance has just recently started since its spread was only found broadly in Java by 2004. The second-generation progeny trial for this purpose was established in 2010, with materials from the best growing 34 families of the first generation. This study aimed to obtain the best pine genotypes through screening the existing natural variations found on important characters. The reported incidence of the trigger was when the experiment was attacked significantly at 30.7% after four years, while some 67 individuals were unexpectedly still performing well after six years. The results show that blocks affect differences for all traits of diameter, aphid resistance, and oleoresin productions, and all families differ except for the west-side yield of oleoresin production. Furthermore, heritability values at individual and family levels were moderate for the diameter (h2i = 0.16; h2f = 0.53) and eastern oleoresin (h2i = 0.14; h2f = 0.42). The gain is 4.3% when 30% of families with the best diameters are retained, while the genetic gain reaches 11% for oleoresin production. As one of the important traits in the breeding program, aphid resistance has a weakly inherited trait (h2i =0.07; h2f =0.29). Interestingly, this trait shows positive moderate genetic correlations with the two essential economic values of diameter (rg = 0.66) and oleoresin production (rg = 0.40). Therefore, the selection of the diameter and oleoresin production will not substantially affect the resistance.

A transcriptome-based association study of growth, wood quality, and oleoresin traits in a slash pine breeding population

Slash pine (Pinus elliottii Engelm.) is an important timber and resin species in the United States, China, Brazil and other countries. Understanding the genetic basis of these traits will accelerate its breeding progress. We carried out a genome-wide association study (GWAS), transcriptome-wide association study (TWAS) and weighted gene co-expression network analysis (WGCNA) for growth, wood quality, and oleoresin traits using 240 unrelated individuals from a Chinese slash pine breeding population. We developed high quality 53,229 single nucleotide polymorphisms (SNPs). Our analysis reveals three main results: (1) the Chinese breeding population can be divided into three genetic groups with a mean inbreeding coefficient of 0.137; (2) 32 SNPs significantly were associated with growth and oleoresin traits, accounting for the phenotypic variance ranging from 12.3% to 21.8% and from 10.6% to 16.7%, respectively; and (3) six genes encoding PeTLP, PeAP2/ERF, PePUP9, PeSLP, PeHSP, and PeOCT1 proteins were identified and validated by quantitative real time polymerase chain reaction for their association with growth and oleoresin traits. These results could be useful for tree breeding and functional studies in advanced slash pine breeding program.

Slash pine is an important source of original timber and resin production on commercial forest plantations. It is necessary to implement precise breeding strategies to improve timber quality and resin yield. However, little is known about the species’ molecular genetic basis. Using a transcriptome dataset with sequencing from 240 individuals in the slash pine population, we combined multiple approaches (based on gene variation, expression variation and co-expression network) to dissect the genetic structure for slash pine major breeding traits. We found that the research population could be divided into three genetic groups with a mean heterozygosity of 0.2246. We also found that six genes with important functions in slash pine resin synthesis and timber formation through association studies. Four new SNPs associatation with the average ring width were also discovered. Our results provide new insights into the molecular genetic basis of important traits in slash pine and provide a comprehensive method for association analyses of conifer tree species with large genome.

Identifying and testing marker‐trait associations for growth and phenology in three pine species: implications for genomic prediction

In tree species, genomic prediction offers the potential to forecast mature trait values in early growth stages, if robust marker–trait associations can be identified. Here we apply a novel multispecies approach using genotypes from a new genotyping array, based on 20,795 single nucleotide polymorphisms (SNPs) from three closely related pine species (Pinus sylvestris, Pinus uncinata and Pinus mugo), to test for associations with growth and phenology data from a common garden study. Predictive models constructed using significantly associated SNPs were then tested and applied to an independent multisite field trial of P. sylvestris and the capability to predict trait values was evaluated. One hundred and eighteen SNPs showed significant associations with the traits in the pine species. Common SNPs (MAF > 0.05) associated with bud set were only found in genes putatively involved in growth and development, whereas those associated with growth and budburst were also located in genes putatively involved in response to environment and, to a lesser extent, reproduction. At one of the two independent sites, the model we developed produced highly significant correlations between predicted values and observed height data (YA, height 2020: r = 0.376, p < 0.001). Predicted values estimated with our budburst model were weakly but positively correlated with duration of budburst at one of the sites (GS, 2015: r = 0.204, p = 0.034; 2018: r = 0.205, p = 0.034–0.037) and negatively associated with budburst timing at the other (YA: r = −0.202, p = 0.046). Genomic prediction resulted in the selection of sets of trees whose mean height was taller than the average for each site. Our results provide tentative support for the capability of prediction models to forecast trait values in trees, while highlighting the need for caution in applying them to trees grown in different environments.

Genetic basis of growth, spring phenology, and susceptibility to biotic stressors in maritime pine

Forest ecosystems are increasingly challenged by extreme events, for example, drought, storms, pest attacks, and fungal pathogen outbreaks, causing severe ecological and economic losses. Understanding the genetic basis of adaptive traits in tree species is of key importance to preserve forest ecosystems, as genetic variation in a trait (i.e., heritability) determines its potential for human-mediated or evolutionary change. Maritime pine (Pinus pinaster Aiton), a conifer widely distributed in southwestern Europe and northwestern Africa, grows under contrasted environmental conditions promoting local adaptation. Genetic variation at adaptive phenotypes, including height, spring phenology, and susceptibility to two fungal pathogens (Diplodia sapinea and Armillaria ostoyae) and an insect pest (Thaumetopoea pityocampa), was assessed in a range-wide clonal common garden of maritime pine. Broad-sense heritability was significant for height (0.219), spring phenology (0.165-0.310), and pathogen susceptibility (necrosis length caused by D. sapinea, 0.152; and by A. ostoyae, 0.021, measured on inoculated, excised branches under controlled conditions), but not for pine processionary moth incidence in the common garden. The correlations of trait variation among populations revealed contrasting trends for pathogen susceptibility to D. sapinea and A. ostoyae with respect to height. Taller trees showed longer necrosis length caused by D. sapinea while shorter trees were more affected by A. ostoyae. Moreover, maritime pine populations from areas with high summer temperatures and frequent droughts were less susceptible to D. sapinea but more susceptible to A. ostoyae. Finally, an association study using 4227 genome-wide SNPs revealed several loci significantly associated with each trait (range of 3-26), including a possibly disease-induced translation initiation factor, eIF-5, associated with needle discoloration caused by D. sapinea. This study provides important insights to develop genetic conservation and breeding strategies integrating species responses to biotic stressors.

Resin ducts as resistance traits in conifers: linking dendrochronology and resin-based defences

Conifers have evolved different chemical and anatomical defences against a wide range of antagonists. Resin ducts produce, store and translocate oleoresin, a complex terpenoid mixture that acts as both a physical and a chemical defence. Although resin duct characteristics (e.g., number, density, area) have been positively related to biotic resistance in several conifer species, the literature reporting this association remains inconclusive. Axial resin ducts recorded in annual growth rings are an archive of annual defensive investment in trees. This whole-life record of defence investment can be analysed using standard dendrochronological procedures, which allows us to assess interannual variability and the effect of understudied drivers of phenotypic variation on resin-based defences. Understanding the sources of phenotypic variation in defences, such as genetic differentiation and environmental plasticity, is essential for assessing the adaptive potential of forest tree populations to resist pests under climate change. Here, we reviewed the evidence supporting the importance of resin ducts in conifer resistance, and summarized current knowledge about the sources of variation in resin duct production. We propose a standardized methodology to measure resin duct production by means of dendrochronological procedures. This approach will illuminate the roles of resin ducts in tree defence across species, while helping to fill pivotal knowledge gaps in plant defence theory, and leading to a robust understanding of the patterns of variation in resin-based defences throughout the tree's lifespan.

A transcriptomic variation map provides insights into the genetic basis of Pinus massoniana Lamb. evolution and the association with oleoresin yield

BACKGROUND

Masson pine (Pinus massoniana Lamb.), the dominant native coniferous species in southern China, is commercially important for supplying timber and oleoresin. However, knowledge of the genetic variability of masson pine germplasm is still limited. In this study, the genetic diversity and population structure of masson pine germplasm were assessed using 204 wild accessions from 10 main distribution regions using 94,194 core single-nucleotide polymorphisms (SNPs) obtained from transcriptome sequencing data.

RESULTS

The average expected heterozygosity was 0.2724, implying abundant genetic diversity within masson pine germplasm. Analysis of molecular variance (AMOVA) revealed that 3.29% of the variation was sourced from genetic differentiation. Structure analysis identified two geographically distinct groups. Discriminant analysis of principal components (DAPC) showed that one of those groups was further divided into two clusters. Sichuan and Chongqing provenance is the geographical origin, which diffused outward along two different lines. Oleoresin yield is reflected in the evolution of the two groups, and exhibits two different trends along the two lines of diffusion. The oleoresin yield may be associated with the genes of chitinase, CYP720B, cytochrome P450, ABC transporter, and AP2/ethylene-responsive transcription factor (ERF) based on SNPs and expression.

CONCLUSIONS

SNP markers from transcriptome sequencing are highly capable of evaluating genetic diversity within different species, as well as the genetic control of objective traits. The functions of these genes will be verified in future studies, and those genes strongly associated with oleoresin yield will be used to improve yields by means of early genotype selection and genetic engineering.

Climate drives intraspecific differentiation in the expression of growth-defence trade-offs in a long-lived pine species

Intraspecific variation in plant defences is expected to be the result of adaptive and plastic responses to environmental conditions, where trade-offs between growth and defences are thought to play a key role shaping phenotypic patterns in defensive investment. Axial resin ducts are costly defensive structures that remain imprinted in the tree rings of conifers, therefore being a valuable proxy of defensive investment along the trees' lifespan. We aimed to disentangle climate-driven adaptive clines and plastic responses to both spatial and temporal environmental variation in resin duct production, and to explore growth-defence trade-offs. To that aim, we applied dendrochronological procedures to quantify annual growth and resin duct production during a 31-year-period in a Mediterranean pine species, including trees from nine populations planted in two common gardens. Both genetic factors and plastic responses modulated annual resin duct production. However, we found no evidence of adaptive clines with climate gradients driving population differentiation. Our results revealed a marked physiological trade-off between growth and defences, where the slope of the trade-off was genetically variable and associated with climatic gradients. Our results help to enlighten the evolutionary patterns and genetic basis of defensive allocation within species, particularly revealing a key role of growth-defence trade-offs.

Transcriptomic and metabolomic analyses reveal several critical metabolic pathways and candidate genes involved in resin biosynthesis in Pinus massoniana

Pine resin, which typically consists of terpenoids, is a natural product used in various industrial applications. Oleoresin can be obtained from the xylem tissue by wounding the stem bark. Pinus massoniana (masson pine) is an important resin-tapping tree species that originated in southern China. Masson pines with different genetic backgrounds typically have different resin-yielding capacities (RYCs). However, the mechanisms underlying high resin yield in masson pines are unclear. The aim of this study was to identify the possible genetic regulation pathways and functional genes that influence the resin yield. In this study, we conducted transcriptomic and metabolomic studies of masson pine secondary xylem with high, medium, and low RYCs. A total of 230,068 unigenes and 3894 metabolites were identified from the tissue of the secondary xylem. Several differentially expressed regulation factors, including WRKY, bHLH, and ERF, and functional genes such as PKc and LRR-RLKs, were identified among these masson pines. The Kyoto Encyclopedia of Genes and Genomes pathways were mainly focused on diterpenoid biosynthesis, plant hormone signal transduction, and ABC transporters. Furthermore, integration of the transcriptomic and metabolomic data indicated that the PKc- and LRR-RLK-related regulatory and metabolic pathways may play critical roles in the biosynthesis of terpenoids. These above results improve our understanding of the biosynthesis mechanism of oleoresin in P. massoniana and facilitate further research work into the functional analysis of these candidate genes.

The Resistance of Scots Pine (Pinus sylvestris L.) Half-sib Families to Heterobasidion annosum

This study observed the genotypic variation among Scots pine (Pinus sylvestris L.) half-sib families’ susceptibility to Heterobasidion annosum. Scots pine susceptibility was tested in 12 half-sib families by inoculating them with four different H. annosum strains. At two, six, and ten months after the inoculations, the susceptibility indicators (incidence rate, pathogen spread, and mortality rate) were compared and the total phenolic compounds (TPC) in the inoculated and control groups determined using the Folin–Ciocalteu method. Among half-sib families, significant differences were found for seedling mortality (range: 1.3%–21.2%); however, the differences in incidence rate (range: 54%–77%) and pathogen spread (range: 24–53 mm) were not significant. The incidence rate among half-sib families correlated positively and significantly (r = 0.72, p < 0.05) with the mortality rate, while the pathogen spread correlated negatively with mortality, although the correlation was not significant (r = −0.29, p > 0.05). The TPC comparison with susceptibility indicators showed that the half-sib families with lower susceptibility were characterized by the ability to increase TPC after inoculation compared to the control group. This tendency was most apparent in stems and roots six and two months after inoculation, respectively. Correlation analyses revealed that higher TPC in stems six months after inoculation determined a lower incidence rate (r = −0.32, p < 0.05), while higher Change in concentration of total phenolic compounds (TPCΔ) indicated a lower pathogen spread (r = −0.60, p < 0.05). The lowest incidence of the pathogen was determined in half-sib families with the highest TPCΔ in the roots two months after inoculation. The lower susceptibility of Scots pine half-sib families is based on a combination of enhanced constitutive and inducible phenolic defense mechanisms. The data may facilitate the selection of Scots pine half-sib families with low susceptibility for breeding programs and forest management strategies.

Identification, expression, and phylogenetic analyses of terpenoid biosynthesis-related genes in secondary xylem of loblolly pine (Pinus taeda L.) based on transcriptome analyses

Loblolly pine (Pinus taeda L.) is one of the most important species for oleoresin (a mixture of terpenoids) in South China. The high oleoresin content of loblolly pine is associated with resistance to bark beetles and other economic benefits. In this study, we conducted transcriptome analyses of loblolly pine secondary xylem to gain insight into the genes involved in terpenoid biosynthesis. A total of 372 unigenes were identified as being critical for oleoresin production, including genes for ATP-binding cassette (ABC) transporters, the cytochrome P450 (CYP) protein family, and terpenoid backbone biosynthesis enzymes. Six key genes involved in terpenoid biosynthetic pathways were selected for multiple sequence alignment, conserved motif prediction, and phylogenetic and expression profile analyses. The protein sequences of all six genes exhibited a higher degree of sequence conservation, and upstream genes were relatively more conserved than downstream genes in terpenoid biosynthetic pathways. The N-terminal regions of these sequences were less conserved than the C-terminal ends, as the N-terminals were quite diverse in both length and composition. The phylogenetic analyses revealed that most genes originated from gene duplication after species divergence, and partial genes exhibited incomplete lineage sorting. In addition, the expression profile analyses showed that all six genes exhibited high expression levels during the high-oleoresin-yielding phase.

Assessment of flux through oleoresin biosynthesis in epithelial cells of loblolly pine resin ducts

The shoot system of pines contains abundant resin ducts, which harbor oleoresins that play important roles in constitutive and inducible defenses. In a pilot study, we assessed the chemical diversity of oleoresins obtained from mature tissues of loblolly pine trees (Pinus taeda L.). Building on these data sets, we designed experiments to assess oleoresin biosynthesis in needles of 2-year-old saplings. Comparative transcriptome analyses of single cell types indicated that genes involved in the biosynthesis of oleoresins are significantly enriched in isolated epithelial cells of resin ducts, compared with those expressed in mesophyll cells. Simulations using newly developed genome-scale models of epithelial and mesophyll cells, which incorporate our data on oleoresin yield and composition as well as gene expression patterns, predicted that heterotrophic metabolism in epithelial cells involves enhanced levels of oxidative phosphorylation and fermentation (providing redox and energy equivalents). Furthermore, flux was predicted to be more evenly distributed across the metabolic network of mesophyll cells, which, in contrast to epithelial cells, do not synthesize high levels of specialized metabolites. Our findings provide novel insights into the remarkable specialization of metabolism in epithelial cells.

Exploring the genetic basis of gene transcript abundance and metabolite levels in loblolly pine (Pinus taeda L.) using association mapping and network construction

BACKGROUND

Identifying genetic variations that shape important complex traits is fundamental to the genetic improvement of important forest tree species, such as loblolly pine (Pinus taeda L.), which is one of the most commonly planted forest tree species in the southern U.S. Gene transcripts and metabolites are important regulatory intermediates that link genetic variations to higher-order complex traits such as wood development and drought response. A few prior studies have associated intermediate phenotypes including mRNA expression and metabolite levels with a limited number of molecular markers, but the identification of genetic variations that regulate intermediate phenotypes needs further investigation.

RESULTS

We identified 1841 single nucleotide polymorphisms (SNPs) associated with 191 gene expression mRNA phenotypes and 524 SNPs associated with 53 metabolite level phenotypes using 2.8 million exome-derived SNPs. The identified SNPs reside in genes with a wide variety of functions. We further integrated the identified SNPs and the associated expressed genes and metabolites into networks. We described the SNP-SNP interactions that significantly impacted the gene transcript abundance and metabolite level in the networks. Key loci and genes in the wood development and drought response networks were identified and analyzed.

CONCLUSIONS

This work provides new candidate genes for research on the genetic basis of gene expression and metabolism linked to wood development and drought response in loblolly pine and highlights the efficiency of using association-mapping-based networks to discover candidate genes with important roles in complex biological processes.

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.

Association genetics of phenolic needle compounds in Norway spruce with variable susceptibility to needle bladder rust

KEY MESSAGE

Accumulation of phenolic needle metabolites in Norway spruce is regulated by many genes with small and additive effects and is correlated with the susceptibility against fungal attack. Norway spruce accumulates high foliar concentrations of secondary phenolic metabolites, with important functions for pathogen defence responses. However, the molecular genetic basis underlying the quantitative variation of phenolic compounds and their role in enhanced resistance of spruce to infection by needle bladder rust are unknown. To address these questions, a set of 1035 genome-wide single nucleotide polymorphisms (SNPs) was associated to the quantitative variation of four simple phenylpropanoids, eight stilbenes, nine flavonoids, six related arithmetic parameters and the susceptibility to infection by Chrysomyxa rhododendri in an unstructured natural population of Norway spruce. Thirty-one significant genetic associations for the flavonoids gallocatechin, kaempferol 3-glucoside and quercetin 3-glucoside and the stilbenes resveratrol, piceatannol, astringin and isorhapontin were discovered, explaining 22-59% of phenotypic variation, and indicating a regulation of phenolic accumulation by many genes with small and additive effects. The phenolics profile differed between trees with high and low susceptibility to the fungus, underlining the importance of phenolic compounds in the defence mechanisms of Norway spruce to C. rhododendri. Results highlight the utility of association studies in non-model tree species and may enable marker-assisted selection of Norway spruce adapted to severe pathogen attack.

Linkage and Association Mapping for Two Major Traits Used in the Maritime Pine Breeding Program: Height Growth and Stem Straightness

Background

Increasing our understanding of the genetic architecture of complex traits, through analyses of genotype-phenotype associations and of the genes/polymorphisms accounting for trait variation, is crucial, to improve the integration of molecular markers into forest tree breeding. In this study, two full-sib families and one breeding population of maritime pine were used to identify quantitative trait loci (QTLs) for height growth and stem straightness, through linkage analysis (LA) and linkage disequilibrium (LD) mapping approaches.

Results

The populations used for LA consisted of two unrelated three-generation full-sib families (n = 197 and n = 477). These populations were assessed for height growth or stem straightness and genotyped for 248 and 217 markers, respectively. The population used for LD mapping consisted of 661 founders of the first and second generations of the breeding program. This population was phenotyped for the same traits and genotyped for 2,498 single-nucleotide polymorphism (SNP) markers corresponding to 1,652 gene loci. The gene-based reference genetic map of maritime pine was used to localize and compare the QTLs detected by the two approaches, for both traits. LA identified three QTLs for stem straightness and two QTLs for height growth. The LD study yielded seven significant associations (P ≤ 0.001): four for stem straightness and three for height growth. No colocalisation was found between QTLs identified by LA and SNPs detected by LD mapping for the same trait.

Conclusions

This study provides the first comparison of LA and LD mapping approaches in maritime pine, highlighting the complementary nature of these two approaches for deciphering the genetic architecture of two mandatory traits of the breeding program.

High-density SNP assay development for genetic analysis in maritime pine (Pinus pinaster)

Maritime pine provides essential ecosystem services in the south-western Mediterranean basin, where it covers around 4 million ha. Its scattered distribution over a range of environmental conditions makes it an ideal forest tree species for studies of local adaptation and evolutionary responses to climatic change. Highly multiplexed single nucleotide polymorphism (SNP) genotyping arrays are increasingly used to study genetic variation in living organisms and for practical applications in plant and animal breeding and genetic resource conservation. We developed a 9k Illumina Infinium SNP array and genotyped maritime pine trees from (i) a three-generation inbred (F2) pedigree, (ii) the French breeding population and (iii) natural populations from Portugal and the French Atlantic coast. A large proportion of the exploitable SNPs (2052/8410, i.e. 24.4%) segregated in the mapping population and could be mapped, providing the densest ever gene-based linkage map for this species. Based on 5016 SNPs, natural and breeding populations from the French gene pool exhibited similar level of genetic diversity. Population genetics and structure analyses based on 3981 SNP markers common to the Portuguese and French gene pools revealed high levels of differentiation, leading to the identification of a set of highly differentiated SNPs that could be used for seed provenance certification. Finally, we discuss how the validated SNPs could facilitate the identification of ecologically and economically relevant genes in this species, improving our understanding of the demography and selective forces shaping its natural genetic diversity, and providing support for new breeding strategies.

Nucleotide polymorphisms in a pine ortholog of the Arabidopsis degrading enzyme cellulase KORRIGAN are associated with early growth performance in Pinus pinaster

We have carried out a candidate-gene-based association genetic study in Pinus pinaster Aiton and evaluated the predictive performance for genetic merit gain of the most significantly associated genes and single nucleotide polymorphisms (SNPs). We used a second generation 384-SNP array enriched with candidate genes for growth and wood properties to genotype mother trees collected in 20 natural populations covering most of the European distribution of the species. Phenotypic data for total height, polycyclism, root-collar diameter and biomass were obtained from a replicated provenance-progeny trial located in two sites with contrasting environments (Atlantic vs Mediterranean climate). General linear models identified strong associations between growth traits (total height and polycyclism) and four SNPs from the korrigan candidate gene, after multiple testing corrections using false discovery rate. The combined genomic breeding value predictions assessed for the four associated korrigan SNPs by ridge regression-best linear unbiased prediction (RR-BLUP) and cross-validation accounted for up to 8 and 15% of the phenotypic variance for height and polycyclic growth, respectively, and did not improve adding SNPs from other growth-related candidate genes. For root-collar diameter and total biomass, they accounted for 1.6 and 1.1% of the phenotypic variance, respectively, but increased to 15 and 4.1% when other SNPs from lp3.1, lp3.3 and cad were included in RR-BLUP models. These results point towards a desirable integration of candidate-gene studies as a means to pre-select relevant markers, and aid genomic selection in maritime pine breeding programs.

Comparison of Genomic Selection Models to Predict Flowering Time and Spike Grain Number in Two Hexaploid Wheat Doubled Haploid Populations

Genomic selection (GS) is becoming an important selection tool in crop breeding. In this study, we compared the ability of different GS models to predict time to young microspore (TYM), a flowering time-related trait, spike grain number under control conditions (SGNC) and spike grain number under osmotic stress conditions (SGNO) in two wheat biparental doubled haploid populations with unrelated parents. Prediction accuracies were compared using BayesB, Bayesian least absolute shrinkage and selection operator (Bayesian LASSO / BL), ridge regression best linear unbiased prediction (RR-BLUP), partial least square regression (PLS), and sparse partial least square regression (SPLS) models. Prediction accuracy was tested with 10-fold cross-validation within a population and with independent validation in which marker effects from one population were used to predict traits in the other population. High prediction accuracies were obtained for TYM (0.51–0.84), whereas moderate to low accuracies were observed for SGNC (0.10–0.42) and SGNO (0.27–0.46) using cross-validation. Prediction accuracies based on independent validation are generally lower than those based on cross-validation. BayesB and SPLS outperformed all other models in predicting TYM with both cross-validation and independent validation. Although the accuracies of all models are similar in predicting SGNC and SGNO with cross-validation, BayesB and SPLS had the highest accuracy in predicting SGNC with independent validation. In independent validation, accuracies of all the models increased by using only the QTL-linked markers. Results from this study indicate that BayesB and SPLS capture the linkage disequilibrium between markers and traits effectively leading to higher accuracies. Excluding markers from QTL studies reduces prediction accuracies.

Genome-Wide Identification of Differentially Expressed Genes Associated with the High Yielding of Oleoresin in Secondary Xylem of Masson Pine (Pinus massoniana Lamb) by Transcriptomic Analysis

Masson pine is an important timber and resource for oleoresin in South China. Increasing yield of oleoresin in stems can raise economic benefits and enhance the resistance to bark beetles. However, the genetic mechanisms for regulating the yield of oleoresin were still unknown. Here, high-throughput sequencing technology was used to investigate the transcriptome and compare the gene expression profiles of high and low oleoresin-yielding genotypes. A total of 40,690,540 reads were obtained and assembled into 137,499 transcripts from the secondary xylem tissues. We identified 84,842 candidate unigenes based on sequence annotation using various databases and 96 unigenes were candidates for terpenoid backbone biosynthesis in pine. By comparing the expression profiles of high and low oleoresin-yielding genotypes, 649 differentially expressed genes (DEGs) were identified. GO enrichment analysis of DEGs revealed that multiple pathways were related to high yield of oleoresin. Nine candidate genes were validated by QPCR analysis. Among them, the candidate genes encoding geranylgeranyl diphosphate synthase (GGPS) and (-)-alpha/beta-pinene synthase were up-regulated in the high oleoresin-yielding genotype, while tricyclene synthase revealed lower expression level, which was in good agreement with the GC/MS result. In addition, DEG encoding ABC transporters, pathogenesis-related proteins (PR5 and PR9), phosphomethylpyrimidine synthase, non-specific lipid-transfer protein-like protein and ethylene responsive transcription factors (ERFs) were also confirmed to be critical for the biosynthesis of oleoresin. The next-generation sequencing strategy used in this study has proven to be a powerful means for analyzing transcriptome variation related to the yield of oleoresin in masson pine. The candidate genes encoding GGPS, (-)-alpha/beta-pinene, tricyclene synthase, ABC transporters, non-specific lipid-transfer protein-like protein, phosphomethylpyrimidine synthase, ERFs and pathogen responses may play important roles in regulating the yield of oleoresin. These DEGs are worthy of special attention in future studies.

A Consensus Genetic Map for Pinus taeda and Pinus elliottii and Extent of Linkage Disequilibrium in Two Genotype-Phenotype Discovery Populations of Pinus taeda

A consensus genetic map for Pinus taeda (loblolly pine) and Pinus elliottii (slash pine) was constructed by merging three previously published P. taeda maps with a map from a pseudo-backcross between P. elliottii and P. taeda. The consensus map positioned 3856 markers via genotyping of 1251 individuals from four pedigrees. It is the densest linkage map for a conifer to date. Average marker spacing was 0.6 cM and total map length was 2305 cM. Functional predictions of mapped genes were improved by aligning expressed sequence tags used for marker discovery to full-length P. taeda transcripts. Alignments to the P. taeda genome mapped 3305 scaffold sequences onto 12 linkage groups. The consensus genetic map was used to compare the genome-wide linkage disequilibrium in a population of distantly related P. taeda individuals (ADEPT2) used for association genetic studies and a multiple-family pedigree used for genomic selection (CCLONES). The prevalence and extent of LD was greater in CCLONES as compared to ADEPT2; however, extended LD with LGs or between LGs was rare in both populations. The average squared correlations, r², between SNP alleles less than 1 cM apart were less than 0.05 in both populations and r² did not decay substantially with genetic distance. The consensus map and analysis of linkage disequilibrium establish a foundation for comparative association mapping and genomic selection in P. taeda and P. elliottii.

Discovering candidate genes that regulate resin canal number in Pinus taeda stems by integrating genetic analysis across environments, ages, and populations

Genetically improving constitutive resin canal development in Pinus stems may enhance the capacity to synthesize terpenes for bark beetle resistance, chemical feedstocks, and biofuels. To discover genes that potentially regulate axial resin canal number (RCN), single nucleotide polymorphisms (SNPs) in 4027 genes were tested for association with RCN in two growth rings and three environments in a complex pedigree of 520 Pinus taeda individuals (CCLONES). The map locations of associated genes were compared with RCN quantitative trait loci (QTLs) in a (P. taeda × Pinus elliottii) × P. elliottii pseudo-backcross of 345 full-sibs (BC1). Resin canal number was heritable (h(2) ˜ 0.12-0.21) and positively genetically correlated with xylem growth (rg ˜ 0.32-0.72) and oleoresin flow (rg ˜ 0.15-0.51). Sixteen well-supported candidate regulators of RCN were discovered in CCLONES, including genes associated across sites and ages, unidirectionally associated with oleoresin flow and xylem growth, and mapped to RCN QTLs in BC1. Breeding is predicted to increase RCN 11% in one generation and could be accelerated with genomic selection at accuracies of 0.45-0.52 across environments. There is significant genetic variation for RCN in loblolly pine, which can be exploited in breeding for elevated terpene content.

Molecular proxies for climate maladaptation in a long-lived tree (Pinus pinaster Aiton, Pinaceae)

Understanding adaptive genetic responses to climate change is a main challenge for preserving biological diversity. Successful predictive models for climate-driven range shifts of species depend on the integration of information on adaptation, including that derived from genomic studies. Long-lived forest trees can experience substantial environmental change across generations, which results in a much more prominent adaptation lag than in annual species. Here, we show that candidate-gene SNPs (single nucleotide polymorphisms) can be used as predictors of maladaptation to climate in maritime pine (Pinus pinaster Aiton), an outcrossing long-lived keystone tree. A set of 18 SNPs potentially associated with climate, 5 of them involving amino acid-changing variants, were retained after performing logistic regression, latent factor mixed models, and Bayesian analyses of SNP-climate correlations. These relationships identified temperature as an important adaptive driver in maritime pine and highlighted that selective forces are operating differentially in geographically discrete gene pools. The frequency of the locally advantageous alleles at these selected loci was strongly correlated with survival in a common garden under extreme (hot and dry) climate conditions, which suggests that candidate-gene SNPs can be used to forecast the likely destiny of natural forest ecosystems under climate change scenarios. Differential levels of forest decline are anticipated for distinct maritime pine gene pools. Geographically defined molecular proxies for climate adaptation will thus critically enhance the predictive power of range-shift models and help establish mitigation measures for long-lived keystone forest trees in the face of impending climate change.

Unique features of the loblolly pine (Pinus taeda L.) megagenome revealed through sequence annotation

The largest genus in the conifer family Pinaceae is Pinus, with over 100 species. The size and complexity of their genomes (∼20-40 Gb, 2n = 24) have delayed the arrival of a well-annotated reference sequence. In this study, we present the annotation of the first whole-genome shotgun assembly of loblolly pine (Pinus taeda L.), which comprises 20.1 Gb of sequence. The MAKER-P annotation pipeline combined evidence-based alignments and ab initio predictions to generate 50,172 gene models, of which 15,653 are classified as high confidence. Clustering these gene models with 13 other plant species resulted in 20,646 gene families, of which 1554 are predicted to be unique to conifers. Among the conifer gene families, 159 are composed exclusively of loblolly pine members. The gene models for loblolly pine have the highest median and mean intron lengths of 24 fully sequenced plant genomes. Conifer genomes are full of repetitive DNA, with the most significant contributions from long-terminal-repeat retrotransposons. In depth analysis of the tandem and interspersed repetitive content yielded a combined estimate of 82%.

In situ genetic association for serotiny, a fire-related trait, in Mediterranean maritime pine (Pinus pinaster)

Wildfire is a major ecological driver of plant evolution. Understanding the genetic basis of plant adaptation to wildfire is crucial, because impending climate change will involve fire regime changes worldwide. We studied the molecular genetic basis of serotiny, a fire-related trait, in Mediterranean maritime pine using association genetics. A single nucleotide polymorphism (SNP) set was used to identify genotype : phenotype associations in situ in an unstructured natural population of maritime pine (eastern Iberian Peninsula) under a mixed-effects model framework. RR-BLUP was used to build predictive models for serotiny in this region. Model prediction power outside the focal region was tested using independent range-wide serotiny data. Seventeen SNPs were potentially associated with serotiny, explaining approximately 29% of the trait phenotypic variation in the eastern Iberian Peninsula. Similar prediction power was found for nearby geographical regions from the same maternal lineage, but not for other genetic lineages. Association genetics for ecologically relevant traits evaluated in situ is an attractive approach for forest trees provided that traits are under strong genetic control and populations are unstructured, with large phenotypic variability. This will help to extend the research focus to ecological keystone non-model species in their natural environments, where polymorphisms acquired their adaptive value.