Conifer defence against insects: microarray gene expression profiling of Sitka spruce (Picea sitchensis) induced by mechanical wounding or feeding by spruce budworms (Choristoneura occidentalis) or white pine weevils (Pissodes strobi) reveals large-scale changes of the host transcriptome

Complementary roles of two classes of defense chemicals in white spruce against spruce budworm

MAIN CONCLUSION

Monoterpenes and phenolics play distinct roles in defending white spruce trees from insect defoliators. Monoterpenes contribute to the toxicity of the foliage, deterring herbivory, whereas phenolics impede budworm growth. This study demonstrates the complex interplay between monoterpenes and phenolics and their collective influence on the defense strategy of white spruce trees against a common insect defoliator. Long-lived coniferous trees display considerable variations in their defensive chemistry. The impact of these defense phenotype variations on insect herbivores of the same species remains to be thoroughly studied, mainly due to challenges in replicating the comprehensive defense profiles of trees under controlled conditions. This study methodically examined the defensive properties of foliar monoterpenes and phenolics across 80 distinct white spruce families. These families were subsequently grouped into two chemotypes based on their foliar monoterpene concentrations. To understand the separate and combined effects of these classes on tree defenses to the eastern spruce budworm, we conducted feeding experiments using actual defense profiles from representative families. Specifically, we assessed budworm response when exposed to substrates amended with phenolics alone or monoterpenes. Our findings indicate that the ratios and amounts of monoterpenes and phenolics present in the white spruce foliage influence the survival of spruce budworms. Phenotypes associated with complete larval mortality exhibited elevated ratios (ranging from 0.4 to 0.6) and concentrations (ranging from 1143 to 1796 ng mg-1) of monoterpenes. Conversely, families characterized by higher phenolic ratios (ranging from 0.62 to 0.77) and lower monoterpene concentrations (ranging from 419 to 985 ng mg-1) were less lethal to the spruce budworm. Both classes of defense compounds contribute significantly to the overall defensive capabilities of white spruce trees. Monoterpenes appear critical in determining the general toxicity of foliage, while phenolics play a role in slowing budworm development, thereby underscoring their collective importance in white spruce defenses.

Non-volatile metabolites mediate plant interactions with insect herbivores

Non-volatile metabolites constitute the bulk of plant biomass. From the perspective of plant-insect interactions, these structurally diverse compounds include nutritious core metabolites and defensive specialized metabolites. In this review, we synthesize the current literature on multiple scales of plant-insect interactions mediated by non-volatile metabolites. At the molecular level, functional genetics studies have revealed a large collection of receptors targeting plant non-volatile metabolites in model insect species and agricultural pests. By contrast, examples of plant receptors of insect-derived molecules remain sparse. For insect herbivores, plant non-volatile metabolites function beyond the dichotomy of core metabolites, classed as nutrients, and specialized metabolites, classed as defensive compounds. Insect feeding tends to elicit evolutionarily conserved changes in plant specialized metabolism, whereas its effect on plant core metabolism varies widely based the interacting species. Finally, several recent studies have demonstrated that non-volatile metabolites can mediate tripartite communication on the community scale, facilitated by physical connections established through direct root-to-root communication, parasitic plants, arbuscular mycorrhizae and the rhizosphere microbiome. Recent advances in both plant and insect molecular biology will facilitate further research on the role of non-volatile metabolites in mediating plant-insect interactions.

Performance of two Ips bark beetles and their associated pathogenic fungi on hosts reflects a species-specific association in the beetle-fungus complex

When Ips bark beetles invade and colonize the host plants, their associated pathogenic fungal partners are carried into the phloem of the host trees. Host trees are lethally attacked by the beetle-fungus complex and the collective damage severely limits forestry production worldwide. It is of great importance to verify whether bark beetles and their associated fungi show concordant performance in terms of biology, physiology, and biochemistry on host trees. In this study, the two Ips bark beetles Ips typographus and Ips subelongatus (Coleoptera: Curculionidae, Scolytinae), their respective associated pathogenic fungi Endoconidiophora polonica and Endoconidiophora fujiensis, and their respective host plants Picea jezoensis and Larix olgensis were selected as test material. Cross-inoculation experiments were conducted indoors and outdoors to investigate the differences in reproduction and development of two beetles and infectivity of two fungi on two plants, as well as the differences in physiological responses of two plants to two fungal infections. The results showed that I. typographus and E. polonica had excellent host performance on P. jezoensis; however, neither successfully colonized and infected L. olgensis. In contrast, I. subelongatus and E. fujiensis showed strong host suitability on L. olgensis and some degree of suitability on P. jezoensis, although the host suitability of P. jezoensis for E. polonica was significantly higher than that for E. fujiensis. In addition, we found that the absolute amount of ergosterol accumulated on the lesion was positively correlated with lesion area. The ergosterol amount and lesion area were both strongly correlated with the release of host monoterpenes, but had no obvious correlation with the concentration of fungi-induced phenols on the lesion area and the side-chain oxidation of lignin in the xylem of the infected sites. Based on these results, we confirmed that "I. typographus-E. polonica" and "I. subelongatus-E. fujiensis" complexes both showed the most suitable consistent performances on their own traditional hosts, establishing a stable species-specific association relationship in these two beetle-fungus complexes, with the "I. subelongatus-E. fujiensis" complex showing broader host suitability. From the perspective of physiological responses of plants to fungal infections, monoterpenes are an important indicator of host suitability.

Plagiodera versicolora feeding induces systemic and sexually differential defense responses in poplars

Dioecious plants have evolved effective defense strategies to deal with various biotic and abiotic stresses. However, little is known regarding sexual differences in their defense against herbivores. In this study, we investigated the mechanism of systemic defense responses in male and female Populus cathayana attacked by Plagiodera versicolora Laicharting. The results revealed that P. cathayana exhibits sexually differential responses to a defoliator. The percentage of damaged leaf area was greater in males than in females. Furthermore, the observed saccharide changes imply that males and females exhibit different response times to defoliators. The contents of flavonoids and anthocyanins were significantly increased in both sexes but were higher in females. Specifically, the jasmonic acid (JA) pathway plays an important role. Expression of pest-related genes further revealed that hormones induce changes in downstream genes and metabolites, and upregulation of JA ZIM-domain (JAZ) and CORONATINE INSENSITIVE 1 (COI1) was more significant in females. In the undamaged adjacent leaves, metabolite and gene changes displayed similar patterns to the damaged local leaves, but levels of JA, JAZ1, and COI1 were higher in females. Therefore, our data confirmed that plants initiate the JA pathway to defend against herbivores, that there is systematic signal transduction, and that this ability is stronger in females than in males. This study provides new insights into the resistance of dioecious plants to herbivory and adds a new theoretical basis for the systemic signal transduction of plants in response to biotic stress.

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.

The involvement of extracellular ATP in regulating the stunted growth of Arabidopsis plants by repeated wounding

BACKGROUND

Extracellular ATP (exATP) has been shown to act as a signal molecule for regulating growth, development, and responses of plants to the external environment.

RESULTS

In this study, we investigated the possible involvement of exATP in regulating the stunted growth caused by repeated wounding. The present work showed that the repeated wounding caused the decreases in leaf area, fresh weight, dry weight, and root length of Arabidopsis seedlings, while the exATP level was enhanced by the repeated wounding. Repeated application of exogenous ATP had similar effects on the plant growth, as the repeated wounding. Through the comparison of p2k1-3 mutant (in which T-DNA disrupted the gene coding P2K1, as exATP receptor) and wide type (WT) plants, it was found that the mutation in P2K1 decreased the sensitivity of plant growth to the repeated wounding and exogenous ATP application. Further works showed that the ibuprofen (IBU, an inhibitor of jasmonate biosynthesis) partially rescued the wound-induced growth degradation. In comparison, the P2K1 mutation partly rescued the wound-induced growth degradation, whereas this mutation failed to do so in the wounded seedlings treated with IBU, indicating that the role of exATP in regulating the growth degradation by repeated wounding could be linked to the JA signaling pathway.

CONCLUSIONS

In conclusion, these results indicate that exATP could be a regulator for the stunted growth of plants by repeated wounding.

Full-Length Transcriptome Sequencing-Based Analysis of Pinus sylvestris var. mongolica in Response to Sirex noctilio Venom

Simple Summary

Sirex noctilio, as a devastating international forestry quarantine pest whose venom can cause a series of physiological changes in the host plants, such as needle wilting, yellowing, decreased transpiration rate and increased respiration rate, etc. In this study, a full-length reference transcript of Pinus sylvestris var. mongolica was constructed by combining second- and third-generation transcriptome sequencing technologies. We also identified the specific expression genes and transcription factors of P. sylvestris var. mongolica under S. noctilio venom and wounding stress. S. noctilio venom mainly induced the expression of genes related to ROS, GAPDH and GPX, and mechanical damage mainly induced the photosynthesis−related genes. The results provide a better understanding of the molecular regulation of pine trees in response to S. noctilio venom.

Abstract

Sirex noctilio is a major international quarantine pest that recently emerged in northeast China to specifically invade conifers. During female oviposition, venom is injected into the host together with its symbiotic fungus to alter the normal Pinus physiology and weaken or even kill the tree. In China, the Mongolian pine (Pinus sylvestris var. mongolica), an important wind-proof and sand-fixing species, is the unique host of S. noctilio. To explore the interplay between S. noctilio venom and Mongolian pine, we performed a transcriptome comparative analysis of a 10-year-old Mongolian pine after wounding and inoculation with S. noctilio venom. The analysis was performed at 12 h, 24 h and 72 h. PacBio ISO-seq was used and integrated with RNA-seq to construct an accurate full-length transcriptomic database. We obtained 52,963 high-precision unigenes, consisting of 48,654 (91.86%) unigenes that were BLASTed to known sequences in the public database and 4309 unigenes without any annotation information, which were presumed to be new genes. The number of differentially expressed genes (DEGs) increased with the treatment time, and the DEGs were most abundant at 72 h. A total of 706 inoculation-specific DEGs (475 upregulated and 231 downregulated) and 387 wounding-specific DEGs (183 upregulated and 204 downregulated) were identified compared with the control. Under venom stress, we identified 6 DEGs associated with reactive oxygen species (ROS) and 20 resistance genes in Mongolian pine. Overall, 52 transcription factors (TFs) were found under venom stress, 45 of which belonged to the AP2/ERF TF family and were upregulated. A total of 13 genes related to the photosystem, 3 genes related photo-regulation, and 9 TFs were identified under wounding stress. In conclusion, several novel putative genes were found in Mongolian pine by PacBio ISO seq. Meanwhile, we also identified various genes that were resistant to S. noctilio venom, such as GAPDH, GPX, CAT, FL2, CERK1, and HSP83A, etc.

Analysis of the transcriptome of the needles and bark of Pinus radiata induced by bark stripping and methyl jasmonate

BACKGROUND

Plants are attacked by diverse insect and mammalian herbivores and respond with different physical and chemical defences. Transcriptional changes underlie these phenotypic changes. Simulated herbivory has been used to study the transcriptional and other early regulation events of these plant responses. In this study, constitutive and induced transcriptional responses to artificial bark stripping are compared in the needles and the bark of Pinus radiata to the responses from application of the plant stressor, methyl jasmonate. The time progression of the responses was assessed over a 4-week period.

RESULTS

Of the 6312 unique transcripts studied, 86.6% were differentially expressed between the needles and the bark prior to treatment. The most abundant constitutive transcripts were related to defence and photosynthesis and their expression did not differ between the needles and the bark. While no differential expression of transcripts were detected in the needles following bark stripping, in the bark this treatment caused an up-regulation and down-regulation of genes associated with primary and secondary metabolism. Methyl jasmonate treatment caused differential expression of transcripts in both the bark and the needles, with individual genes related to primary metabolism more responsive than those associated with secondary metabolism. The up-regulation of genes related to sugar break-down and the repression of genes related with photosynthesis, following both treatments was consistent with the strong down-regulation of sugars that has been observed in the same population. Relative to the control, the treatments caused a differential expression of genes involved in signalling, photosynthesis, carbohydrate and lipid metabolism as well as defence and water stress. However, non-overlapping transcripts were detected between the needles and the bark, between treatments and at different times of assessment. Methyl jasmonate induced more transcriptional responses in the bark than bark stripping, although the peak of expression following both treatments was detected 7 days post treatment application. The effects of bark stripping were localised, and no systemic changes were detected in the needles.

CONCLUSION

There are constitutive and induced differences in the needle and bark transcriptome of Pinus radiata. Some expression responses to bark stripping may differ from other biotic and abiotic stresses, which contributes to the understanding of plant molecular responses to diverse stresses. Whether the gene expression changes are heritable and how they differ between resistant and susceptible families identified in earlier studies needs further investigation.

Transcriptomic Profiling Reveals Shared Signalling Networks Between Flower Development and Herbivory-Induced Responses in Tomato

Most flowering plants must defend themselves against herbivores for survival and attract pollinators for reproduction. Although traits involved in plant defence and pollinator attraction are often localised in leaves and flowers, respectively, they will show a diffuse evolution if they share the same molecular machinery and regulatory networks. We performed RNA-sequencing to characterise and compare transcriptomic changes involved in herbivory-induced defences and flower development, in tomato leaves and flowers, respectively. We found that both the herbivory-induced responses and flower development involved alterations in jasmonic acid signalling, suppression of primary metabolism and reprogramming of secondary metabolism. We identified 411 genes that were involved in both processes, a number significantly higher than expected by chance. Genetic manipulation of key regulators of induced defences also led to the expression changes in the same genes in both leaves and flowers. Targeted metabolomic analysis showed that among closely related tomato species, jasmonic acid and α-tomatine are correlated in flower buds and herbivory-induced leaves. These findings suggest that herbivory-induced responses and flower development share a common molecular machinery and likely have coevolved in nature.

Constitutive and insect-induced transcriptomes of weevil-resistant and susceptible Sitka spruce

Spruce weevil (Pissodes strobi) is a significant pest of regenerating spruce (Picea) and pine (Pinus) forests in North America. Weevil larvae feed in the bark, phloem, cambium, and outer xylem of apical shoots, causing stunted growth or mortality of young trees. We identified and characterized constitutive and weevil-induced patterns of Sitka spruce (Picea sitchensis) transcriptomes in weevil-resistant (R) and susceptible (S) trees using RNA sequencing (RNA-seq) and differential expression (DE) analyses. We developed a statistical model for the analysis of RNA-seq data from treatment experiments with a 2 × 3 factorial design to differentiate insect-induced responses from the effects of mechanical damage. Across the different comparisons, we identified two major transcriptome contrasts: A large set of genes that was constitutively DE between R and S trees, and another set of genes that was DE in weevil-induced S-trees. The constitutive transcriptome unique to R trees appeared to be attuned to defense, while the constitutive transcriptome unique to S trees was enriched for growth-related transcripts. Notably, a set of transcripts annotated as "fungal" was detected consistently in the transcriptomes. Fungal transcripts were identified as DE in the comparison of R and S trees and in the weevil-affected DE transcriptome of S trees, suggesting a potential microbiome role in this conifer-insect interaction.

Bruchid beetle ovipositioning mediated defense responses in black gram pods

BACKGROUND

Black gram [Vigna mungo (L)] seeds are a rich source of digestible protein and dietary fibre, both for human and animal consumption. However, the quality and quantity of the Vigna seeds are severely affected by bruchid beetles during storage. Therefore, analyses of the expression of the bruchid induced transcript dynamics in black gram pods would be helpful to understand the underlying defense mechanism against bruchid oviposition.

RESULTS

We used the RNAseq approach to survey the changes in transcript profile in the developing seeds of a moderately resistant cultivar IC-8219 against bruchid oviposition using a susceptible cultivar T-9 as a control. A total of 96,084,600 and 99,532,488 clean reads were generated from eight (4 each) samples of IC-8219 and T-9 cultivar, respectively. Based on the BLASTX search against the NR database, 32,584 CDSs were generated of which 31,817 CDSs were significantly similar to Vigna radiata, a close relative of Vigna mungo. The IC-8219 cultivar had 630 significantly differentially expressed genes (DEGs) of which 304 and 326 genes up and down-regulated, respectively. However, in the T-9 cultivar, only 168 DEGs were identified of which 142 and 26 genes up and down-regulated, respectively. The expression analyses of 10 DEGs by qPCR confirmed the accuracy of the RNA-Seq data. Gene Ontology and KEGG pathway analyses helped us to better understand the role of these DEGs in oviposition mediated defense response of black gram. In both the cultivars, the most significant transcriptomic changes in response to the oviposition were related to the induction of defense response genes, transcription factors, secondary metabolites, enzyme inhibitors, and signal transduction pathways. It appears that the bruchid ovipositioning mediated defense response in black gram is induced by SA signaling pathways and defense genes such as defensin, genes for secondary metabolites, and enzyme inhibitors could be potential candidates for resistance to bruchids.

CONCLUSION

We generated a transcript profile of immature black gram pods upon bruchid ovipositioning by de novo assembly and studied the underlying defense mechanism of a moderately resistant cultivar.

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.

Combining transcriptomics and genetic linkage based information to identify candidate genes associated with Heterobasidion-resistance in Norway spruce

The Heterobasidion annosum s.l species complex comprises the most damaging forest pathogens to Norway spruce. We revisited previously identified Quantitative Trait Loci (QTLs) related to Heterobasidion-resistance in Norway spruce to identify candidate genes associated with these QTLs. We identified 329 candidate genes associated with the resistance QTLs using a gene-based composite map for Pinaceae. To evaluate the transcriptional responses of these candidate genes to H. parviporum, we inoculated Norway spruce plants and sequenced the transcriptome of the interaction at 3 and 7 days post inoculation. Out of 298 expressed candidate genes 124 were differentially expressed between inoculation and wounding control treatment. Interestingly, PaNAC04 and two of its paralogs in the subgroup III-3 of the NAC family transcription factors were found to be associated with one of the QTLs and was also highly induced in response to H. parviporum. These genes are possibly involved in the regulation of biosynthesis of flavonoid compounds. Furthermore, several of the differentially expressed candidate genes were associated with the phenylpropanoid pathway including a phenylalanine ammonia-lyase, a cinnamoyl-CoA reductase, a caffeoyl-CoA O-methyltransferase and a PgMYB11-like transcription factor gene. Combining transcriptome and genetic linkage analyses can help identifying candidate genes for functional studies and molecular breeding in non-model species.

Association genetics identifies a specifically regulated Norway spruce laccase gene, PaLAC5, linked to Heterobasidion parviporum resistance

It is important to improve the understanding of the interactions between the trees and pathogens and integrate this knowledge about disease resistance into tree breeding programs. The conifer Norway spruce (Picea abies) is an important species for the forest industry in Europe. Its major pathogen is Heterobasidion parviporum, causing stem and root rot. In this study, we identified 11 Norway spruce QTLs (Quantitative trait loci) that correlate with variation in resistance to H. parviporum in a population of 466 trees by association genetics. Individual QTLs explained between 2.1 and 5.2% of the phenotypic variance. The expression of candidate genes associated with the QTLs was analysed in silico and in response to H. parviporum hypothesizing that (a) candidate genes linked to control of fungal sapwood growth are more commonly expressed in sapwood, and; (b) candidate genes associated with induced defences are respond to H. parviporum inoculation. The Norway spruce laccase PaLAC5 associated with control of lesion length development is likely to be involved in the induced defences. Expression analyses showed that PaLAC5 responds specifically and strongly in close proximity to the H. parviporum inoculation. Thus, PaLAC5 may be associated with the lignosuberized boundary zone formation in bark adjacent to the inoculation site.

Physical and Physiological Monitoring on Red Palm Weevil-Infested Oil Palms

The red palm weevil (RPW) is a stem boring Coleoptera that decimates host palm trees from within. The challenge of managing this pest is due to a lack of physical symptoms during the early stages of infestation. Investigating the physiological changes that occur within RPW-infested palm trees may be useful in establishing a new approach in RPW detection. In this study, the effects of RPW infestation were investigated in Elaeis guineensis by observing changes in physical and physiological parameters during the progress of infestation by visual inspection and the comparison of growth, gas exchange, stomatal conductance, and chlorophyll content between the non-infested control, physically wounded, and RPW-infested E. guineensis groups. During the study period, four distinct levels of physical infestation were observed and recorded. The RPW-infested group displayed significantly lower maximum photosynthesis activity (A_max) starting from the third week post-infestation. However, growth in terms of change in plant height and stem circumference, leaves’ stomatal conductance, and chlorophyll content were not significantly different between the three groups during the duration of the study. The significant drop in photosynthesis was observed one week before physical changes appeared. This suggests the promising utilisation of photosynthesis activity as a signal for detecting RPW infestation at the early stage of attacks, which could be useful for integration in integrated pest management (IPM).

Characterization and Analysis of the Full-Length Transcriptomes of Multiple Organs in Pseudotaxus chienii (W.C.Cheng) W.C.Cheng

Pseudotaxus chienii, a rare tertiary relict species with economic and ecological value, is a representative of the monotypic genus Pseudotaxus that is endemic to China. P. chienii can adapt well to habitat isolation and ecological heterogeneity under a variety of climate and soil conditions, and is able to survive in harsh environments. However, little is known about the molecular and genetic resources of this long-lived conifer. Herein, we sequenced the transcriptomes of four organs of P. chienii using the PacBio Isoform Sequencing and Illumina RNA Sequencing platforms. Based on the PacBio Iso-Seq data, we obtained 44,896, 58,082, 50,485, and 67,638 full-length unigenes from the root, stem, leaf, and strobilus, respectively, with a mean length of 2692 bp, and a mean N50 length of 3010.75 bp. We then comprehensively annotated these unigenes. The number of organ-specific expressed unigenes ranged from 4393 in leaf to 9124 in strobilus, suggesting their special roles in physiological processes, organ development, and adaptability in the different four organs. A total of 16,562 differentially expressed genes (DEGs) were identified among the four organs and clustered into six subclusters. The gene families related to biotic/abiotic factors, including the TPS, CYP450, and HSP families, were characterized. The expression levels of most DEGs in the phenylpropanoid biosynthesis pathway and plant–pathogen interactions were higher in the root than in the three other organs, suggesting that root constitutes the main organ of defensive compound synthesis and accumulation and has a stronger ability to respond to stress. The sequences were analyzed to predict transcription factors, long non-coding RNAs, and alternative splicing events. The expression levels of most DEGs of C2H2, C3H, bHLH, and bZIP families in the root and stem were higher than those in the leaf and strobilus, indicating that these TFs may play a crucial role in the survival of the root and stem. These results comprise the first comprehensive gene expression profiles obtained for different organs of P. chienii. Our findings will facilitate further studies on the functional genomics, adaptive evolution, and phylogeny of P. chienii, and lay the foundation for the development of conservation strategies for this endangered conifer.

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

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

Hydroxyacetophenone defenses in white spruce against spruce budworm

We review a recently discovered white spruce (Picea glauca) chemical defense against spruce budworm (Choristoneura fumiferana) involving hydroxyacetophenones. These defense metabolites detected in the foliage accumulate variably as the aglycons, piceol and pungenol, or the corresponding glucosides, picein and pungenin. We summarize current knowledge of the genetic, genomic, molecular, and biochemical underpinnings of this defense and its effects on C. fumiferana. We present an update with new results on the ontogenic variation and the phenological window of this defense, including analysis of transcript responses in P. glauca to C. fumiferana herbivory. We also discuss this chemical defense from an evolutionary and a breeding context.

Multi-trait genomic selection for weevil resistance, growth, and wood quality in Norway spruce

Plantation-grown trees have to cope with an increasing pressure of pest and disease in the context of climate change, and breeding approaches using genomics may offer efficient and flexible tools to face this pressure. In the present study, we targeted genetic improvement of resistance of an introduced conifer species in Canada, Norway spruce (Picea abies (L.) Karst.), to the native white pine weevil (Pissodes strobi Peck). We developed single- and multi-trait genomic selection (GS) models and selection indices considering the relationships between weevil resistance, intrinsic wood quality, and growth traits. Weevil resistance, acoustic velocity as a proxy for mechanical wood stiffness, and average wood density showed moderate-to-high heritability and low genotype-by-environment interactions. Weevil resistance was genetically positively correlated with tree height, height-to-diameter at breast height (DBH) ratio, and acoustic velocity. The accuracy of the different GS models tested (GBLUP, threshold GBLUP, Bayesian ridge regression, BayesCπ) was high and did not differ among each other. Multi-trait models performed similarly as single-trait models when all trees were phenotyped. However, when weevil attack data were not available for all trees, weevil resistance was more accurately predicted by integrating genetically correlated growth traits into multi-trait GS models. A GS index that corresponded to the breeders' priorities achieved near maximum gains for weevil resistance, acoustic velocity, and height growth, but a small decrease for DBH. The results of this study indicate that it is possible to breed for high-quality, weevil-resistant Norway spruce reforestation stock with high accuracy achieved from single-trait or multi-trait GS.

A molecular and genomic reference system for conifer defence against insects

Insect pests are part of natural forest ecosystems contributing to forest rejuvenation but can also cause ecological disturbance and economic losses that are expected to increase with climate change. The white pine or spruce weevil (Pissodes strobi) is a pest of conifer forests in North America. Weevil-host interactions with various spruce (Picea) species have been explored as a genomic and molecular reference system for conifer defence against insects. Interactions occur in two major phases of the insect life cycle. In the exophase, adult weevils are free-moving and display behaviour of host selection for oviposition that is affected by host traits. In the endophase, insects live within the host where mobility and development from eggs to young adults are affected by a complex system of host defences. Genetic resistance exists in several spruce species and involves synergism of constitutive and induced chemical and physical defences that comprise the conifer defence syndrome. Here, we review conifer defences that disrupt the weevil life cycle and mechanisms by which trees resist weevil attack. We highlight molecular and genomic aspects and a possible role for the weevil microbiome. Knowledge of this conifer defence system is supporting forest health strategies and tree breeding for insect resistance.

Growth-defense trade-off regulated by hormones in grass plants growing under different grazing intensities

Herbivory creates conflicts between a plant's need to allocate resources for growth and defense. It is not yet clear how plants rebalance resource utilization between growth and defense in response to increasing grazing intensity. We measured characteristics of the primary and secondary metabolism of Leymus chinensis at five levels of grazing intensity (control, light, moderate, heavy and extremely heavy). Furthermore, we evaluated hormone signaling by quantifying the impact of key hormones on plant growth and defense. Under light grazing intensity, indole-3-acetic acid and jasmonates appeared to promote the growth of L. chinensis through a high photosynthetic rate, high water-use efficiency and high soluble protein contents, whereas abscisic acid decreased these properties. Under moderate grazing intensity, L. chinensis had a low photosynthetic capacity but greater production of secondary metabolites (tannins, total flavonoids and total phenols), possibly induced by salicylic acid. When the grazing pressure further intensified, L. chinensis translocated more carbohydrates to its roots in order to survive and regrow. Leymus chinensis therefore exhibited a trade-off between growth and defense in order to survive and reproduce under herbivory. Plants developed different mechanisms to enhance their grazing tolerance by means of hormonal regulation.

Effect of Long-Term Feeding by Borers on the Antibiotic Properties of Corn Stems

Plant long-term response against chewing insects could become stronger than initial reactions and even turn into systemic. The objectives of the present study were 1) to evaluate whether the long-running attack to the stem by corn borers can improve the stem antibiotic properties; 2) to check whether hydroxycinnamic acids could be involved in this antibiotic response; 3) and to check whether elicitation by Sesamia nonagrioides Lef. (Lepidoptera: Noctuidae) regurgitant could activate long-term plant responses. In this sense, we observed that long-term feeding by S. nonagrioides larvae induced genotype-dependent changes in stem antibiosis and phenolic profiles, but the hydroxycinnamate content does not have a significant role in the systemic defense induced by the attack. In addition, response to long-term feeding by larvae could not be fully mimicked by elicitation using S. nonagrioides regurgitant alone. For the first time, it has been demonstrated that 'long-term' attack to the stem by corn borers can increase the stem antibiotic properties, and this has to be considered attending to breeding strategies.

Cytochrome P450 enzymes: A driving force of plant diterpene diversity

In plant terpene biosynthesis, oxidation of the hydrocarbon backbone produced by terpene synthases is typically carried out by cytochrome P450 oxygenases (CYPs). The modifications introduced by CYPs include hydroxylations, sequential oxidations at one position and ring rearrangements and closures. These reactions significantly expand the structural diversity of terpenoids, but also provide anchoring points for further decorations by various transferases. In recent years, there has been a significant increase in reports of CYPs involved in plant terpene pathways. Plant diterpenes represent an important class of metabolites that includes hormones and a number of industrially relevant compounds such as pharmaceutical, aroma or food ingredients. In this review, we provide a comprehensive survey on CYPs reported to be involved in plant diterpene biosynthesis to date. A phylogenetic analysis showed that only few CYP clans are represented in diterpene biosynthesis, namely CYP71, CYP85 and CYP72. Remarkably few CYP families and subfamilies within those clans are involved, indicating specific expansion of these clades in plant diterpene biosynthesis. Nonetheless, the evolutionary trajectory of CYPs of specialized diterpene biosynthesis is diverse. Some are recently derived from gibberellin biosynthesis, while others have a more ancient history with recent expansions in specific plant families. Among diterpenoids, labdane-related diterpenoids represent a dominant class. The availability of CYPs from diverse plant species able to catalyze oxidations in specific regions of the labdane-related backbones provides opportunities for combinatorial biosynthesis to produce novel diterpene compounds that can be screened for biological activities of interest.

Functions of stone cells and oleoresin terpenes in the conifer defense syndrome

Conifers depend on complex defense systems against herbivores. Stone cells (SC) and oleoresin are physical and chemical defenses of Sitka spruce that have been separately studied in previous work. Weevil oviposit at the tip of the previous year's apical shoot (PYAS). We investigated interactions between weevil larvae and trees in controlled oviposition experiments with resistant (R) and susceptible (S) Sitka spruce. R trees have an abundance of SC in the PYAS cortex. SC are mostly absent in S trees. R trees and S trees also differ in the composition of oleoresin terpenes. Transcriptomes of R and S trees revealed differences in long-term weevil-induced responses. Performance of larvae was significantly reduced on R trees compared with S trees under experimental conditions that mimicked natural oviposition behavior at apical shoot tips and may be attributed to the effects of SC. In oviposition experiments designed for larvae to feed below the area of highest SC abundance, larvae showed an unusual feeding behavior and oleoresin appeared to function as the major defense. The results support a role for both SC and oleoresin terpenes and possible synergies between these traits in the defense syndrome of weevil-resistant Sitka spruce.

Inducibility of Plant Secondary Metabolites in the Stem Predicts Genetic Variation in Resistance Against a Key Insect Herbivore in Maritime Pine

Resistance to herbivores and pathogens is considered a key plant trait with strong adaptive value in trees, usually involving high concentrations of a diverse array of plant secondary metabolites (PSM). Intraspecific genetic variation and plasticity of PSM are widely known. However, their ecology and evolution are unclear, and even the implication of PSM as traits that provide direct effective resistance against herbivores is currently questioned. We used control and methyl jasmonate (MJ) induced clonal copies of genotypes within families from ten populations of the main distribution range of maritime pine to exhaustively characterize the constitutive and induced profile and concentration of PSM in the stem phloem, and to measure insect herbivory damage as a proxy of resistance. Then, we explored whether genetic variation in resistance to herbivory may be predicted by the constitutive concentration of PSM, and the role of its inducibility to predict the increase in resistance once the plant is induced. We found large and structured genetic variation among populations but not between families within populations in resistance to herbivory. The MJ-induction treatment strongly increased resistance to the weevil in the species, and the genetic variation in the inducibility of resistance was significantly structured among populations, with greater inducibility in the Atlantic populations. Genetic variation in resistance was largely explained by the multivariate concentration and profile of PSM at the genotypic level, rather than by bivariate correlations with individual PSM, after accounting for genetic relatedness among genotypes. While the constitutive concentration of the PSM blend did not show a clear pattern of resistance to herbivory, specific changes in the chemical profile and the increase in concentration of the PSM blend after MJ induction were related to increased resistance. To date, this is the first example of a comprehensive and rigorous approach in which inducibility of PSM in trees and its implication in resistance was analyzed excluding spurious associations due to genetic relatedness, often overlooked in intraspecific studies. Here we provide evidences that multivariate analyses of PSM, rather than bivariate correlations, provide more realistic information about the potentially causal relationships between PSM and resistance to herbivory in pine trees.

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.

Partial correlation analysis of transcriptomes helps detangle the growth and defense network in spruce

In plants, there can be a trade-off between resource allocations to growth vs defense. Here, we use partial correlation analysis of gene expression to make inferences about the nature of this interaction. We studied segregating progenies of Interior spruce subject to weevil attack. In a controlled experiment, we measured pre-attack plant growth and post-attack damage with several morphological measures, and profiled transcriptomes of 188 progeny. We used partial correlations of individual transcripts (expressed sequence tags, ESTs) with pairs of growth/defense traits to identify important nodes and edges in the inferred underlying gene network, for example, those pairs of growth/defense traits with high mutual correlation with a single EST transcript. We give a method to identify such ESTs. A terpenoid ABC transporter gene showed strongest correlations (P = 0.019); its transcript represented a hub within the compact 166-member gene-gene interaction network (P = 0.004) of the negative genetic correlations between growth and subsequent pest attack. A small 21-member interaction network (P = 0.004) represented the uncovered positive correlations. Our study demonstrates partial correlation analysis identifies important gene networks underlying growth and susceptibility to the weevil in spruce. In particular, we found transcripts that strongly modify the trade-off between growth and defense, and allow identification of networks more central to the trade-off.

Cell-type- and tissue-specific transcriptomes of the white spruce (Picea glauca) bark unmask fine-scale spatial patterns of constitutive and induced conifer defense

Plant defenses often involve specialized cells and tissues. In conifers, specialized cells of the bark are important for defense against insects and pathogens. Using laser microdissection, we characterized the transcriptomes of cortical resin duct cells, phenolic cells and phloem of white spruce (Picea glauca) bark under constitutive and methyl jasmonate (MeJa)-induced conditions, and we compared these transcriptomes with the transcriptome of the bark tissue complex. Overall, ~3700 bark transcripts were differentially expressed in response to MeJa. Approximately 25% of transcripts were expressed in only one cell type, revealing cell specialization at the transcriptome level. MeJa caused cell-type-specific transcriptome responses and changed the overall patterns of cell-type-specific transcript accumulation. Comparison of transcriptomes of the conifer bark tissue complex and specialized cells resolved a masking effect inherent to transcriptome analysis of complex tissues, and showed the actual cell-type-specific transcriptome signatures. Characterization of cell-type-specific transcriptomes is critical to reveal the dynamic patterns of spatial and temporal display of constitutive and induced defense systems in a complex plant tissue or organ. This was demonstrated with the improved resolution of spatially restricted expression of sets of genes of secondary metabolism in the specialized cell types.

Transcriptomic responses of Solanum dulcamara to natural and simulated herbivory

Plants are attacked by diverse herbivores and respond with manifold defence responses. To study transcriptional and other early regulation events of these plant responses, herbivory is often simulated to standardize the temporal and spatial dynamics that vary tremendously for natural herbivory. Yet, to what extent such simulations of herbivory are able to elicit the same plant response as real herbivory remains largely undetermined. We examined the transcriptional response of a wild model plant to herbivory by lepidopteran larvae and to a commonly used herbivory simulation by applying the larvae's oral secretions to standardized wounds. We designed a microarray for Solanum dulcamara and showed that the transcriptional responses to real and to simulated herbivory by Spodoptera exigua overlapped moderately by about 40%. Interestingly, certain responses were mimicked better than others; 60% of the genes upregulated but not even a quarter of the genes downregulated by herbivory were similarly affected by application of oral secretions to wounds. While the regulation of genes involved in signalling, defence and water stress was mimicked well by the simulated herbivory, most of the genes related to photosynthesis, carbohydrate- and lipid metabolism were exclusively regulated by real herbivory. Thus, wounding and application of oral secretions decently mimics herbivory-induced defence responses but likely not the reallocation of primary metabolites induced by real herbivory.

Simulated herbivory in chickpea causes rapid changes in defense pathways and hormonal transcription networks of JA/ethylene/GA/auxin within minutes of wounding

Chickpea (C. arietinum L.) is an important pulse crop in Asian and African countries that suffers significant yield losses due to attacks by insects like H. armigera. To obtain insights into early responses of chickpea to insect attack, a transcriptomic analysis of chickpea leaves just 20 minutes after simulated herbivory was performed, using oral secretions of H. armigera coupled with mechanical wounding. Expression profiles revealed differential regulation of 8.4% of the total leaf transcriptome with 1334 genes up-regulated and 501 down-regulated upon wounding at log₂-fold change (|FC| ≤ −1 and ≥1) and FDR value ≤ 0.05. In silico analysis showed the activation of defenses through up-regulation of genes of the phenylpropanoid pathway, pathogenesis, oxidases and CYTP450 besides differential regulation of kinases, phosphatases and transcription factors of the WRKY, MYB, ERFs, bZIP families. A substantial change in the regulation of hormonal networks was observed with up-regulation of JA and ethylene pathways and suppression of growth associated hormone pathways like GA and auxin within 20 minutes of wounding. Secondary qPCR comparison of selected genes showed that oral secretions often increased differential expression relative to mechanical damage alone. The studies provide new insights into early wound responses in chickpea.

VqDUF642, a gene isolated from the Chinese grape Vitis quinquangularis, is involved in berry development and pathogen resistance

The DUF642 gene VqDUF642 , isolated from the Chinese grape species V. quinquangularis accession Danfeng-2, participates in berry development and defense responses against Erysiphe necator and Botrytis cinerea. The proteins with domains of unknown function 642 (DUF642) comprise a large protein family according to cell wall proteomic analyses in plants. However, the works about functional characterization of DUF642s in plant development and resistance to pathogens are scarce. In this study, a gene encoding a DUF642 protein was isolated from Chinese grape V. quinquangularis accession Danfeng-2, and designated as VqDUF642. Its full-length cDNA contains a 1107-bp open reading frame corresponding to a deduced 368-amino acid protein. Multiple sequence alignments and phylogenetic analysis showed that VqDUF642 is highly homologous to one of the DUF642 proteins (VvDUF642) in V. vinifera. The VqDUF642 was localized to the cell wall of tobacco epidermal cells. Accumulation of VqDUF642 protein and VqDUF642 transcript abundance increased at the later stage of grape berry development in Danfeng-2. Overexpression of VqDUF642 in transgenic tomato plants accelerated plant growth and reduced susceptibility to Botrytis cinerea. Transgenic Thompson Seedless grapevine plants overexpressing VqDUF642 exhibited enhanced resistance to Erysiphe necator and B. cinerea. Moreover, VqDUF642 overexpression affected the expression of a couple of pathogenesis-related (PR) genes in transgenic tomato and grapevine upon pathogen inoculation. Taken together, these results suggest that VqDUF642 is involved in plant development and defense against pathogenic infections.

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

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

Expression Patterns of Genes Involved in Ascorbate-Glutathione Cycle in Aphid-Infested Maize (Zea mays L.) Seedlings

Reduced forms of ascorbate (AsA) and glutathione (GSH) are among the most important non-enzymatic foliar antioxidants in maize (Zea mays L.). The survey was aimed to evaluate impact of bird cherry-oat aphid (Rhopalosiphum padi L.) or grain aphid (Sitobion avenae F.) herbivory on expression of genes related to ascorbate-glutathione (AsA-GSH) cycle in seedlings of six maize varieties (Ambrozja, Nana, Tasty Sweet, Touran, Waza, Złota Karłowa), differing in resistance to the cereal aphids. Relative expression of sixteen maize genes encoding isoenzymes of ascorbate peroxidase (APX1, APX2, APX3, APX4, APX5, APX6, APX7), monodehydroascorbate reductase (MDHAR1, MDHAR2, MDHAR3, MDHAR4), dehydroascorbate reductase (DHAR1, DHAR2, DHAR3) and glutathione reductase (GR1, GR2) was quantified. Furthermore, effect of hemipterans’ attack on activity of APX, MDHAR, DHAR and GR enzymes, and the content of reduced and oxidized ascorbate and glutathione in maize plants were assessed. Seedling leaves of more resistant Z. mays varieties responded higher elevations in abundance of target transcripts. In addition, earlier and stronger aphid-triggered changes in activity of APX, MDHAR, DHAR and GR enzymes, and greater modulations in amount of the analyzed antioxidative metabolites were detected in foliar tissues of highly resistant Ambrozja genotype in relation to susceptible Tasty Sweet plants.

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

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

Transcriptomic Response to Nitric Oxide Treatment in Larix olgensis Henry

Larix olgensis Henry is an important coniferous species found in plantation forests in northeastern China, but it is vulnerable to pathogens. Nitric oxide (NO) is an important molecule involved in plant resistance to pathogens. To study the regulatory role of NO at the transcriptional level, we characterized the transcriptomic response of L. olgensis seedlings to sodium nitroprusside (SNP, NO donor) using Illumina sequencing and de novo transcriptome assembly. A significant number of putative metabolic pathways and functions associated with the unique sequences were identified. Genes related to plant pathogen infection (FLS2, WRKY33, MAPKKK, and PR1) were upregulated with SNP treatment. This report describes the potential contribution of NO to disease resistance in L. olgensis as induced by biotic stress. Our results provide a substantial contribution to the genomic and transcriptomic resources for L. olgensis, as well as expanding our understanding of the involvement of NO in defense responses at the transcriptional level.

The NAC transcription factor family in maritime pine (Pinus Pinaster): molecular regulation of two genes involved in stress responses

BACKGROUND

NAC transcription factors comprise a large plant-specific gene family involved in the regulation of diverse biological processes. Despite the growing number of studies on NAC transcription factors in various species, little information is available about this family in conifers. The goal of this study was to identify the NAC transcription family in maritime pine (Pinus pinaster), to characterize ATAF-like genes in response to various stresses and to study their molecular regulation.

METHODS

We have isolated two maritime pine NAC genes and using a transient expression assay in N. benthamiana leaves estudied the promoter jasmonate response.

RESULTS

In this study, we identified 37 NAC genes from maritime pine and classified them into six main subfamilies. The largest group includes 12 sequences corresponding to stress-related genes. Two of these NAC genes, PpNAC2 and PpNAC3, were isolated and their expression profiles were examined at various developmental stages and in response to various types of stress. The expression of both genes was strongly induced by methyl jasmonate (MeJA), mechanical wounding, and high salinity. The promoter regions of these genes were shown to contain cis-elements involved in the stress response and plant hormonal regulation, including E-boxes, which are commonly found in the promoters of genes that respond to jasmonate, and binding sites for bHLH proteins. Using a transient expression assay in N. benthamiana leaves, we found that the promoter of PpNAC3 was rapidly induced upon MeJA treatment, while this response disappeared in plants in which the transcription factor NbbHLH2 was silenced.

CONCLUSION

Our results suggest that PpNAC2 and PpNAC3 encode stress-responsive NAC transcription factors involved in the jasmonate response in pine. Furthermore, these data also suggest that the jasmonate signaling pathway is conserved between angiosperms and gymnosperms. These findings may be useful for engineering stress tolerance in pine via biotechnological approaches.

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.

Activation of defence pathways in Scots pine bark after feeding by pine weevil (Hylobius abietis)

BACKGROUND

During their lifetime, conifer trees are exposed to numerous herbivorous insects. To protect themselves against pests, trees have developed a broad repertoire of protective mechanisms. Many of the plant's defence reactions are activated upon an insect attack, and the underlying regulatory mechanisms are not entirely understood yet, in particular in conifer trees. Here, we present the results of our studies on the transcriptional response and the volatile compounds production of Scots pine (Pinus sylvestris) upon the large pine weevil (Hylobius abietis) feeding.

RESULTS

Transcriptional response of Scots pine to the weevil attack was investigated using a novel customised 36.4 K Pinus taeda microarray. The weevil feeding caused large-scale changes in the pine transcriptome. In total, 774 genes were significantly up-regulated more than 4-fold (p≤0.05), whereas 64 genes were significantly down-regulated more than 4-fold. Among the up-regulated genes, we could identify genes involved in signal perception, signalling pathways, transcriptional regulation, plant hormone homeostasis, secondary metabolism and defence responses. The weevil feeding on stem bark of pine significantly increased the total emission of volatile organic compounds from the undamaged stem bark area. The emission levels of monoterpenes and sesquiterpenes were also increased. Interestingly, we could not observe any correlation between the increased production of the terpenoid compounds and expression levels of the terpene synthase-encoding genes.

CONCLUSIONS

The obtained data provide an important insight into the transcriptional response of conifer trees to insect herbivory and illustrate the massive changes in the host transcriptome upon insect attacks. Moreover, many of the induced pathways are common between conifers and angiosperms. The presented results are the first ones obtained by the use of a microarray platform with an extended coverage of pine transcriptome (36.4 K cDNA elements). The platform will further facilitate the identification of resistance markers with the direct relevance for conifer tree breeding.

Influence of water deficit on the molecular responses of Pinus contorta × Pinus banksiana mature trees to infection by the mountain pine beetle fungal associate, Grosmannia clavigera

Conifers exhibit a number of constitutive and induced mechanisms to defend against attack by pests and pathogens such as mountain pine beetle (Dendroctonus ponderosae Hopkins) and their fungal associates. Ecological studies have demonstrated that stressed trees are more susceptible to attack by mountain pine beetle than their healthy counterparts. In this study, we tested the hypothesis that water deficit affects constitutive and induced responses of mature lodgepole pine × jack pine hybrids (Pinus contorta Dougl. ex Loud. var. latifolia Engelm. ex S. Wats. × Pinus banksiana Lamb.) to inoculation with the mountain pine beetle fungal associate Grosmannia clavigera (Robinson-Jeffrey and Davidson) Zipfel, de Beer and Wingfield. The degree of stress induced by the imposed water-deficit treatment was sufficient to reduce photosynthesis. Grosmannia clavigera-induced lesions exhibited significantly reduced dimensions in water-deficit trees relative to well-watered trees at 5 weeks after inoculation. Treatment-associated cellular-level changes in secondary phloem were also observed. Quantitative RT-PCR was used to analyze transcript abundance profiles of 18 genes belonging to four families classically associated with biotic and abiotic stress responses: aquaporins (AQPs), dehydration-responsive element binding (DREB), terpene synthases (TPSs) and chitinases (CHIs). Transcript abundance profiles of a TIP2 AQP and a TINY-like DREB decreased significantly in fungus-inoculated trees, but not in response to water deficit. One TPS, Pcb(+)-3-carene synthase, and the Class II CHIs PcbCHI2.1 and PcbCHI2.2 showed increased expression under water-deficit conditions in the absence of fungal inoculation, while another TPS, Pcb(E)-β-farnesene synthase-like, and two CHIs, PcbCHI1.1 and PcbCHI4.1, showed attenuated expression under water-deficit conditions in the presence of fungal inoculation. The effects were observed both locally and systemically. These results demonstrate that both constitutive and induced carbon- and nitrogen-based defenses are affected by water deficit, suggesting potential consequences for mountain pine beetle dynamics, particularly in novel environments.

Transcriptome analysis of soybean lines reveals transcript diversity and genes involved in the response to common cutworm (Spodoptera litura Fabricius) feeding

The interaction between soybeans and the destructive common cutworm insect is complicated. In this paper, the time course of induced responses to common cutworm was characterized in two soybean lines, and the results showed that the induced resistance peaked at different times in the resistant (WX) and susceptible (NN) soybean lines. Two sets of transcriptome profiles from the WX and NN lines at the peak of their induced resistance were compared using microarray analysis. In total, 827 and 349 transcripts were differentially expressed in the WX and NN lines, respectively, with 80 probes common regulated and seven regulated in the opposite direction. All common- and unique-regulated genes were grouped into 10 functional categories based on sequence similarity searches, which showed that most of the genes were related to stress and defence responses. qRT-PCR analysis of 22 genes confirmed the results of the microarray analysis. The spatiotemporal expression patterns of the six genes revealed the consistency of systemic expression levels with the timing of the resistance response observed in the bioassay experiments. In summary, we described the conceptual model of induced resistance in two soybean lines and provided the first large-scale survey of common cutworm-induced defence transcripts in soybean.

Transcriptome analysis of callus from Picea balfouriana

BACKGROUND

Picea likiangensis var. balfouriana (Rehd. et Wils.) Hillier ex Slavin (also known as Picea balfouriana) is an ecologically and economically important conifer that grows rapidly under optimum conditions and produces high-quality wood. It has a wide geographic distribution and is prevalent in southwest and eastern regions of China. Under suboptimal conditions, P. balfouriana grows slowly, which restricts its cultivation. Somatic embryogenesis has been used in the mass propagation of commercial species. However, low initiation rates are a common problem and the mechanisms involved in the induction of somatic embryogenesis are not fully understood. To understand the molecular mechanisms regulating somatic embryogenesis in P. balfouriana, high-throughput RNA-seq technology was used to investigate the transcriptomes of embryogenic and non-embryogenic tissues from three P. balfouriana genotypes. We compared the genes expressed in these tissues to identify molecular markers with embryogenic potential.

RESULTS

A total of 55,078,846 nucleotide sequence reads were obtained for the embryogenic and non-embryogenic tissues of P. balfouriana, and 49.56% of them uniquely matched 22,295 (84.3%) of the 26,437 genes in the Picea abies genome database (Nature 497: 579-584, 2013). Differential gene expression analysis identified 1,418 differentially expressed genes (false discovery rate <0.0001; fold change ≥2) in the embryogenic tissues relative to the non-embryogenic tissues, including 431 significantly upregulated and 987 significantly downregulated genes. KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis revealed that the most significantly altered genes were involved in plant hormone signal transduction, metabolic pathways (starch and sucrose metabolism), and phenylalanine metabolism.

CONCLUSIONS

We found that the initiation of embryogenic tissues affected gene expression in many KEGG pathways, but predominantly in plant hormone signal transduction, plant-pathogen interaction, and starch and sucrose metabolism. The changes in multiple pathways related to induction in the P. balfouriana embryogenic tissues described here, will contribute to a more comprehensive understanding of the mechanisms involved in the initiation of somatic embryogenesis. Additionally, we found that somatic embryogenesis receptor kinase (SERK), arabinogalactan proteins, and members of the WUS-related homeobox protein family may play important roles and could act as molecular markers in the early stage of somatic embryogenesis, as reported previously.

Cloning and characterization of chitinases from interior spruce and lodgepole pine

Chitinases have been implicated in the defence of conifers against insects and pathogens. cDNA for six chitinases were cloned from interior spruce (Picea glauca x engelmannii) and four from lodgepole pine (Pinus contorta). The cloned interior spruce chitinases were annotated class I PgeChia1-1 and PgeChia1-2, class II PgeChia2-1, class IV PgeChia4-1, and class VII PgeChia7-1 and PgeChia7-2; lodgepole pine chitinases were annotated class I PcChia1-1, class IV PcChia4-1, and class VII PcChia7-1 and PcChia7-2. Chitinases were expressed in Escherichia coli with maltose-binding-protein tags and soluble proteins purified. Functional characterization demonstrated chitinolytic activity for the three class I chitinases PgeChia1-1, PgeChia1-2 and PcChia1-1. Transcript analysis established strong induction of most of the tested chitinases, including all three class I chitinases, in interior spruce and lodgepole pine in response to inoculation with bark beetle associated fungi (Leptographium abietinum and Grosmannia clavigera) and in interior spruce in response to weevil (Pissodes strobi) feeding. Evidence of chitinolytic activity and inducibility by fungal and insect attack support the involvement of these chitinases in conifer defense.

Metabolic responses of Quercus ilex seedlings to wounding analysed with nuclear magnetic resonance profiling

Plants defend themselves against herbivory at several levels. One of these is the synthesis of inducible chemical defences. Using NMR metabolomic techniques, we studied the metabolic changes of plant leaves after a wounding treatment simulating herbivore attack in the Mediterranean sclerophyllous tree Quercus ilex. First, an increase in glucose content was observed in wounded plants. There was also an increase in the content of C-rich secondary metabolites such as quinic acid and quercitol, both related to the shikimic acid pathway and linked to defence against biotic stress. There was also a shift in N-storing amino acids, from leucine and isoleucine to asparagine and choline. The observed higher content of asparagine is related to the higher content of choline through serine that was proved to be the precursor of choline. Choline is a general anti-herbivore and pathogen deterrent. The study shows the rapid metabolic response of Q. ilex in defending its leaves, based on a rapid increase in the production of quinic acid, quercitol and choline. The results also confirm the suitability of (1)H NMR-based metabolomic profiling studies to detect global metabolome shifts after wounding stress in tree leaves, and therefore its suitability in ecometabolomic studies.

Overexpression of an isoprenyl diphosphate synthase in spruce leads to unexpected terpene diversion products that function in plant defense

Spruce (Picea spp.) and other conifers employ terpenoid-based oleoresin as part of their defense against herbivores and pathogens. The short-chain isoprenyl diphosphate synthases (IDS) are situated at critical branch points in terpene biosynthesis, producing the precursors of the different terpenoid classes. To determine the role of IDS and to create altered terpene phenotypes for assessing the defensive role of terpenoids, we overexpressed a bifunctional spruce IDS, a geranyl diphosphate and geranylgeranyl diphosphate synthase in white spruce (Picea glauca) saplings. While transcript level (350-fold), enzyme activity level (7-fold), and in planta geranyl diphosphate and geranylgeranyl diphosphate levels (4- to 8-fold) were significantly increased in the needles of transgenic plants, there was no increase in the major monoterpenes and diterpene acids of the resin and no change in primary isoprenoids, such as sterols, chlorophylls, and carotenoids. Instead, large amounts of geranylgeranyl fatty acid esters, known from various gymnosperm and angiosperm plant species, accumulated in needles and were shown to act defensively in reducing the performance of larvae of the nun moth (Lymantria monacha), a conifer pest in Eurasia. These results show the impact of overexpression of an IDS and the defensive role of an unexpected accumulation product of terpenoid biosynthesis with the potential for a broader function in plant protection.

Large-scale analysis of differential gene expression in coffee genotypes resistant and susceptible to leaf miner-toward the identification of candidate genes for marker assisted-selection

BACKGROUND

A successful development of herbivorous insects into plant tissues depends on coordination of metabolic processes. Plants have evolved complex mechanisms to recognize such attacks, and to trigger a defense response. To understand the transcriptional basis of this response, we compare gene expression profiles of two coffee genotypes, susceptible and resistant to leaf miner (Leucoptera coffella). A total of 22000 EST sequences from the Coffee Genome Database were selected for a microarray analysis. Fluorescence probes were synthesized using mRNA from the infested and non-infested coffee plants. Array hybridization, scanning and data normalization were performed using Nimble Scan® e ArrayStar® platforms. Genes with foldchange values +/-2 were considered differentially expressed. A validation of 18 differentially expressed genes was performed in infected plants using qRT-PCR approach.

RESULTS

The microarray analysis indicated that resistant plants differ in gene expression profile. We identified relevant transcriptional changes in defense strategies before insect attack. Expression changes (>2.00-fold) were found in resistant plants for 2137 genes (1266 up-regulated and 873 down-regulated). Up-regulated genes include those responsible for defense mechanisms, hypersensitive response and genes involved with cellular function and maintenance. Also, our analyses indicated that differential expression profiles between resistant and susceptible genotypes are observed in the absence of leaf-miner, indicating that defense is already build up in resistant plants, as a priming mechanism. Validation of selected genes pointed to four selected genes as suitable candidates for markers in assisted-selection of novel cultivars.

CONCLUSIONS

Our results show evidences that coffee defense responses against leaf-miner attack are balanced with other cellular functions. Also analyses suggest a major metabolic reconfiguration that highlights the complexity of this response.

Differential peroxidase activities in three different crops upon insect feeding

Peroxidases are the ubiquitous enzyme and reported to be present in all living genera. They catalyses reduction of peroxide and generate reactive oxygen species. In the present study we demonstrated that insect infestation induces peroxidase activity in sap and total soluble protein (TSP) of plant leaves. Three important crop plants viz. tomato, cowpea and cotton were used for this study. After infestation of chewing insect, Peroxidase activity in the sap and TSP of all the studied plants were enhanced in the range of 1.6 to 3.14 fold. Similar observations were also obtained with feeding of sap sucking insects, in which increment in peroxidase activity of sap and TSP was in the range of 1.8 to 2.53 fold. Enhanced peroxidase activity was reconfirmed by in-gel peroxidase assay. Enzyme kinetic study showed turn over efficiency of peroxidase from cotton (~101.3 min(-1)) was almost similar to tomato (~100.8 min(-1)) but higher than cowpea (~98.21 in(-1)). MS/MS analysis of observed band showed significant similarity with the reported peroxidases in database.

De novo characterization of Larix gmelinii (Rupr.) Rupr. transcriptome and analysis of its gene expression induced by jasmonates

BACKGROUND

Larix gmelinii is a dominant tree species in China's boreal forests and plays an important role in the coniferous ecosystem. It is also one of the most economically important tree species in the Chinese timber industry due to excellent water resistance and anti-corrosion of its wood products. Unfortunately, in Northeast China, L. gmelinii often suffers from serious attacks by diseases and insects. The application of exogenous volatile semiochemicals may induce and enhance its resistance against insect or disease attacks; however, little is known regarding the genes and molecular mechanisms related to induced resistance.

RESULTS

We performed de novo sequencing and assembly of the L. gmelinii transcriptome using a short read sequencing technology (Illumina). Chemical defenses of L. gmelinii seedlings were induced with jasmonic acid (JA) or methyl jasmonate (MeJA) for 6 hours. Transcriptomes were compared between seedlings induced by JA, MeJA and untreated controls using a tag-based digital gene expression profiling system. In a single run, 25,977,782 short reads were produced and 51,157 unigenes were obtained with a mean length of 517 nt. We sequenced 3 digital gene expression libraries and generated between 3.5 and 5.9 million raw tags, and obtained 52,040 reliable reference genes after removing redundancy. The expression of disease/insect-resistance genes (e.g., phenylalanine ammonialyase, coumarate 3-hydroxylase, lipoxygenase, allene oxide synthase and allene oxide cyclase) was up-regulated. The expression profiles of some abundant genes under different elicitor treatment were studied by using real-time qRT-PCR.The results showed that the expression levels of disease/insect-resistance genes in the seedling samples induced by JA and MeJA were higher than those in the control group. The seedlings induced with MeJA elicited the strongest increases in disease/insect-resistance genes.

CONCLUSIONS

Both JA and MeJA induced seedlings of L. gmelinii showed significantly increased expression of disease/insect-resistance genes. MeJA seemed to have a stronger induction effect than JA on expression of disease/insect-resistance related genes. This study provides sequence resources for L. gmelinii research and will help us to better understand the functions of disease/insect-resistance genes and the molecular mechanisms of secondary metabolisms in L. gmelinii.