Saturated genic SNP mapping identified functional candidates and selection tools for the Pinus monticola Cr2 locus controlling resistance to white pine blister rust

Transcriptional Profiling of Early Defense Response to White Pine Blister Rust Infection in Pinus albicaulis (Whitebark Pine)

Pathogen perception generates the activation of signal transduction cascades to host defense. White pine blister rust (WPBR) is caused by Cronartium ribicola J.C. Fisch and affects a number of species of Pinus. One of the most severely affected species is Pinus albicaulis Engelm (whitebark pine). WPBR resistance in the species is a polygenic and complex trait that requires an optimized immune response. We identified early responses in 2-year-old seedlings after four days of fungal inoculation and compared the underlying transcriptomic response with that of healthy non-inoculated individuals. A de novo transcriptome assembly was constructed with 56,796 high quality-annotations derived from the needles of susceptible and resistant individuals in a resistant half-sib family. Differential expression analysis identified 599 differentially expressed transcripts, from which 375 were upregulated and 224 were downregulated in the inoculated seedlings. These included components of the initial phase of active responses to abiotic factors and stress regulators, such as those involved in the first steps of flavonoid biosynthesis. Four days after the inoculation, infected individuals showed an overexpression of chitinases, reactive oxygen species (ROS) regulation signaling, and flavonoid intermediates. Our research sheds light on the first stage of infection and emergence of disease symptoms among whitebark pine seedlings. RNA sequencing (RNA-seq) data encoding hypersensitive response, cell wall modification, oxidative regulation signaling, programmed cell death, and plant innate immunity were differentially expressed during the defense response against C. ribicola.

Comparative Association Mapping Reveals Conservation of Major Gene Resistance to White Pine Blister Rust in Southwestern White Pine (Pinus strobiformis) and Limber Pine (P. flexilis)

All native North American white pines are highly susceptible to white pine blister rust (WPBR) caused by Cronartium ribicola. Understanding genomic diversity and molecular mechanisms underlying genetic resistance to WPBR remains one of the great challenges in improvement of white pines. To compare major gene resistance (MGR) present in two species, southwestern white pine (Pinus strobiformis) Cr3 and limber pine (P. flexilis) Cr4, we performed association analyses of Cr3-controlled resistant traits using single nucleotide polymorphism (SNP) assays designed with Cr4-linked polymorphic genes. We found that ∼70% of P. flexilis SNPs were transferable to P. strobiformis. Furthermore, several Cr4-linked SNPs were significantly associated with the Cr3-controlled traits in P. strobiformis families. The most significantly associated SNP (M326511_1126R) almost colocalized with Cr4 on the Pinus consensus linkage group 8, suggesting that Cr3 and Cr4 might be the same R locus, or have localizations very close to each other in the syntenic region of the P. strobiformis and P. flexilis genomes. M326511_1126R was identified as a nonsynonymous SNP, causing amino acid change (Val₃₇₆Ile) in a putative pectin acetylesterase, with coding sequences identical between the two species. Moreover, top Cr3-associated SNPs were further developed as TaqMan genotyping assays, suggesting their usefulness as marker-assisted selection (MAS) tools to distinguish genotypes between quantitative resistance and MGR. This work demonstrates the successful transferability of SNP markers between two closely related white pine species in the hybrid zone, and the possibility for deployment of MAS tools to facilitate long-term WPBR management in P. strobiformis breeding and conservation.

Achievements and Challenges of Genomics-Assisted Breeding in Forest Trees: From Marker-Assisted Selection to Genome Editing

Forest tree breeding efforts have focused mainly on improving traits of economic importance, selecting trees suited to new environments or generating trees that are more resilient to biotic and abiotic stressors. This review describes various methods of forest tree selection assisted by genomics and the main technological challenges and achievements in research at the genomic level. Due to the long rotation time of a forest plantation and the resulting long generation times necessary to complete a breeding cycle, the use of advanced techniques with traditional breeding have been necessary, allowing the use of more precise methods for determining the genetic architecture of traits of interest, such as genome-wide association studies (GWASs) and genomic selection (GS). In this sense, main factors that determine the accuracy of genomic prediction models are also addressed. In turn, the introduction of genome editing opens the door to new possibilities in forest trees and especially clustered regularly interspaced short palindromic repeats and CRISPR-associated protein 9 (CRISPR/Cas9). It is a highly efficient and effective genome editing technique that has been used to effectively implement targetable changes at specific places in the genome of a forest tree. In this sense, forest trees still lack a transformation method and an inefficient number of genotypes for CRISPR/Cas9. This challenge could be addressed with the use of the newly developing technique GRF-GIF with speed breeding.

Fine dissection of limber pine resistance to Cronartium ribicola using targeted sequencing of the NLR family

BACKGROUND

Proteins with nucleotide binding site (NBS) and leucine-rich repeat (LRR) domains (NLR) make up one of most important resistance (R) families for plants to resist attacks from various pathogens and pests. The available transcriptomes of limber pine (Pinus flexilis) allow us to characterize NLR genes and related resistance gene analogs (RGAs) in host resistance against Cronartium ribicola, the causal fungal pathogen of white pine blister rust (WPBR) on five-needle pines throughout the world. We previously mapped a limber pine major gene locus (Cr4) that confers complete resistance to C. ribicola on the Pinus consensus linkage group 8 (LG-8). However, genetic distribution of NLR genes as well as their divergence between resistant and susceptible alleles are still unknown.

RESULTS

To identify NLR genes at the Cr4 locus, the present study re-sequenced a total of 480 RGAs using targeted sequencing in a Cr4-segregated seed family. Following a call of single nucleotide polymorphisms (SNPs) and genetic mapping, a total of 541 SNPs from 155 genes were mapped across 12 LGs. Three putative NLR genes were newly mapped in the Cr4 region, including one that co-segregated with Cr4. The tight linkage of NLRs with Cr4-controlled phenotypes was further confirmed by bulked segregation analysis (BSA) using extreme-phenotype genome-wide association study (XP-GWAS) for significance test. Local tandem duplication in the Cr4 region was further supported by syntenic analysis using the sugar pine genome sequence. Significant gene divergences have been observed in the NLR family, revealing that diversifying selection pressures are relatively higher in local duplicated genes. Most genes showed similar expression patterns at low levels, but some were affected by genetic background related to disease resistance. Evidence from fine genetic dissection, evolutionary analysis, and expression profiling suggests that two NLR genes are the most promising candidates for Cr4 against WPBR.

CONCLUSION

This study provides fundamental insights into genetic architecture of the Cr4 locus as well as a set of NLR variants for marker-assisted selection in limber pine breeding. Novel NLR genes were identified at the Cr4 locus and the Cr4 candidates will aid deployment of this R gene in combination with other major/minor genes in the limber pine breeding program.

Comparative Transcriptomics and RNA-Seq-Based Bulked Segregant Analysis Reveals Genomic Basis Underlying Cronartium ribicola vcr2 Virulence

Breeding programs of five-needle pines have documented both major gene resistance (MGR) and quantitative disease resistance (QDR) to Cronartium ribicola (Cri), a non-native, invasive fungal pathogen causing white pine blister rust (WPBR). WPBR is one of the most deadly forest diseases in North America. However, Cri virulent pathotypes have evolved and can successfully infect and kill trees carrying resistance (R) genes, including vcr2 that overcomes MGR conferred by the western white pine (WWP, Pinus monticola) R gene (Cr2). In the absence of a reference genome, the present study generated a vcr2 reference transcriptome, consisting of about 20,000 transcripts with 1,014 being predicted to encode secreted proteins (SPs). Comparative profiling of transcriptomes and secretomes revealed vcr2 was significantly enriched for several gene ontology (GO) terms relating to oxidation-reduction processes and detoxification, suggesting that multiple molecular mechanisms contribute to pathogenicity of the vcr2 pathotype for its overcoming Cr2. RNA-seq-based bulked segregant analysis (BSR-Seq) revealed genome-wide DNA variations, including about 65,617 single nucleotide polymorphism (SNP) loci in 7,749 polymorphic genes shared by vcr2 and avirulent (Avcr2) pathotypes. An examination of the distribution of minor allele frequency (MAF) uncovered a high level of genomic divergence between vcr2 and Avcr2 pathotypes. By integration of extreme-phenotypic genome-wide association (XP-GWAS) analysis and allele frequency directional difference (AFDD) mapping, we identified a set of vcr2-associated SNPs within functional genes, involved in fungal virulence and other molecular functions. These included six SPs that were top candidate effectors with putative activities of reticuline oxidase, proteins with common in several fungal extracellular membrane (CFEM) domain or ferritin-like domain, polysaccharide lyase, rds1p-like stress responsive protein, and two Cri-specific proteins without annotation. Candidate effectors and vcr2-associated genes provide valuable resources for further deciphering molecular mechanisms of virulence and pathogenicity by functional analysis and the subsequent development of diagnostic tools for monitoring the virulence landscape in the WPBR pathosystems.

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

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

Association Mapping and Development of Marker-Assisted Selection Tools for the Resistance to White Pine Blister Rust in the Alberta Limber Pine Populations

Since its introduction to North America in the early 1900s, white pine blister rust (WPBR) caused by the fungal pathogen Cronartium ribicola has resulted in substantial economic losses and ecological damage to native North American five-needle pine species. The high susceptibility and mortality of these species, including limber pine (Pinus flexilis), creates an urgent need for the development and deployment of resistant germplasm to support recovery of impacted populations. Extensive screening for genetic resistance to WPBR has been underway for decades in some species but has only started recently in limber pine using seed families collected from wild parental trees in the USA and Canada. This study was conducted to characterize Alberta limber pine seed families for WPBR resistance and to develop reliable molecular tools for marker-assisted selection (MAS). Open-pollinated seed families were evaluated for host reaction following controlled infection using C. ribicola basidiospores. Phenotypic segregation for presence/absence of stem symptoms was observed in four seed families. The segregation ratios of these families were consistent with expression of major gene resistance (MGR) controlled by a dominant R locus. Based on linkage disequilibrium (LD)-based association mapping used to detect single nucleotide polymorphism (SNP) markers associated with MGR against C. ribicola, MGR in these seed families appears to be controlled by Cr4 or other R genes in very close proximity to Cr4. These associated SNPs were located in genes involved in multiple molecular mechanisms potentially underlying limber pine MGR to C. ribicola, including NBS-LRR genes for recognition of C. ribicola effectors, signaling components, and a large set of defense-responsive genes with potential functions in plant effector-triggered immunity (ETI). Interactions of associated loci were identified for MGR selection in trees with complex genetic backgrounds. SNPs with tight Cr4-linkage were further converted to TaqMan assays to confirm their effectiveness as MAS tools. This work demonstrates the successful translation and deployment of molecular genetic knowledge into specific MAS tools that can be easily applied in a selection or breeding program to efficiently screen MGR against WPBR in Alberta limber pine populations.

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.

Construction of genetic linkage map and identification of a novel major locus for resistance to pine wood nematode in Japanese black pine (Pinus thunbergii)

BACKGROUND

Pine wilt disease (PWD), which is caused by the pine wood nematode (PWN) Bursaphelenchus xylophilus, is currently the greatest threat to pine forests in Europe and East Asian countries including Japan. Constructing a detailed linkage map of DNA markers and identifying PWD resistance genes/loci lead to improved resistance in Pinus thunbergii, as well as other Pinus species that are also susceptible to PWD.

RESULTS

A total F₁ mapping population of 188 individuals derived from a cross between the PWD-resistant P. thunbergii varieties 'Tanabe 54' (resistant rank 2 to PWD) and 'Tosashimizu 63' (resistant rank 4 to PWD) was inoculated with PWN, and was evaluated for disease symptoms. To perform linkage analysis for PWN resistance, a set of three maps was constructed; two parental maps generated using the integrated two-way pseudo-testcross method, and a consensus map with population-type cross-pollination. The linkage map of 'Tanabe 54' consisted of 167 loci, and covered 14 linkage groups (LGs), with a total genetic distance of 1214.6 cM. The linkage map of 'Tosashimizu 63' consisted of 252 loci, and covered 14 LGs, with a total genetic distance of 1422.1 cM. The integrated consensus map comprised 12 LGs with the basic chromosome number of P. thunbergii, and a total genetic distance of 1403.6 cM. Results from quantitative trait loci (QTL) analysis using phenotype data and linkage maps indicated that PWN resistance is controlled by a single dominant allele, which was derived from the 'Tanabe 54' female parent. This major QTL was located on linkage group 3 and was designated PWD1 for PINE WILT DISEASE 1.

CONCLUSIONS

The PWD1 locus is a major resistance QTL located on the Pinus consensus LG03 that acts in a dominant manner to confer pine wood nematode resistance. Information from the present study will be useful for P. thunbergii breeding programs to improve resistance to PWD, and also to help identify susceptibility genes in Pinus species.

Plant nonhost resistance: paradigms and new environments

Nonhost resistance (NHR) protects plants from a large and diverse array of potential phytopathogens. Each phytopathogen can parasitise some plant species, but most plant species are nonhosts that are innately immune due to a series of physical, chemical and inducible defenses these nonadapted pathogens cannot overcome. New evidence supports the NHR paradigm that posits the inability of potential pathogens to colonise nonhost plants is frequently due to molecular incompatibility between pathogen virulence factors and plant cellular targets. While NHR is durable, it is not insurmountable. Environmental changes can facilitate pathogen host jumps or alternatively result in new encounters between previously isolated plant species and pathogens. Climate change is predicted to substantially alter the current distribution of plants and their pathogens which could result in parasitism of new plant species.

Multiomics analysis reveals a distinct mechanism of oleaginousness in the emerging model alga Chromochloris zofingiensis

Chromochloris zofingiensis, featured due to its capability to simultaneously synthesize triacylglycerol (TAG) and astaxanthin, is emerging as a leading candidate alga for production uses. To better understand the oleaginous mechanism of this alga, we conducted a multiomics analysis by systematically integrating time-resolved transcriptomes, lipidomes and metabolomes in response to nitrogen deprivation. The data analysis unraveled the distinct mechanism of TAG accumulation, which involved coordinated stimulation of multiple biological processes including supply of energy and reductants, carbon reallocation from protein and starch, and 'pushing' and 'pulling' carbon to TAG synthesis. Unlike the model alga Chlamydomonas, de novo fatty acid synthesis in C. zofingiensis was promoted, together with enhanced turnover of both glycolipids and phospholipids, supporting the drastic need of acyls for TAG assembly. Moreover, genomewide analysis identified many key functional enzymes and transcription factors that had engineering potential for TAG modulation. Two genes encoding glycerol-3-phosphate acyltransferase (GPAT), the first committed enzyme for TAG assembly, were found in the C. zofingiensis genome; in vivo functional characterization revealed that extrachloroplastic GPAT instead of chloroplastic GPAT played a central role in TAG synthesis. These findings illuminate distinct oleaginousness mechanisms in C. zofingiensis and pave the way towards rational manipulation of this alga to becone an emerging model for trait improvements.

Identification and Functional Characterization of an Effector Secreted by Cronartium ribicola

Cri-9402 was identified as a protein effector from Cronartium ribicola, based on the effect of its expression on growth of Pseudomonas syringae Psm ES4326 introduced into transiently transformed tobacco leaves and stably transformed Arabidopsis seedlings. In tobacco leaves transiently expressing its coding sequence, growth of P. syringae Psm ES4326 was inhibited. Expression of pathogenesis-related (PR) protein 2 (PR2), PR4a, endochitinase B, hypersensitive-related 201 (HSR201), HSR203J, and proteinase inhibitor 1 was upregulated but expression of PR1, coronatine insensitive 1, and abscisic acid 1 was significantly suppressed. In transformed Arabidopsis seedlings, the effector stimulated growth of P. syringae Psm ES4326; significantly suppressed expression of PR1, PR2, nonexpresser of pathogenesis-related genes 1 (NPR1), NPR3, NPR4, phytoalexin deficient 4, and salicylic acid induction deficient 2; and enhanced expression of plant defensin 1.2 (PDF1.2). The above results showed that the majority of responses to this effector in tobacco leaves were converse to those in transformed Arabidopsis. We could conclude that Cri-9402 promoted disease resistance in tobacco leaves and disease susceptibility in Arabidopsis seedlings. Its transcript was mainly expressed in aeciospores of C. ribicola and was probably involved in production or germination of aeciospores, and it was an effector potentially functioning in white pine-blister rust interactions.

Characterization of the western white pine TIR-NBS-LRR (PmTNL2) gene by transcript profiling and promoter analysis

Proteins with nucleotide-binding site (NBS) and leucine-rich repeats (LRRs) have been reported to play important roles in plant disease resistance, growth, and development. However, no comprehensive analysis of this protein family has been performed in conifers. Here we report that the Pinus monticola PmTNL2 gene is a member of the NBS-LRR superfamily. Quantitative reverse transcription-PCR (qRT-PCR) analysis revealed that the PmTNL2 transcript was expressed in a tissue-specific pattern with extensive regulation by various environmental stimuli in western white pine seedlings, suggesting its wide involvement in stress defense and diverse developmental processes. In silico analysis of the PmTNL2 promoter region revealed multiple cis-regulatory elements characterized with potential functions for development-, light-, and stress-regulated transcript expression. Expression patterns were largely confirmed by PmTNL2 promoter-directed reporter gene expression using stable transgenic Arabidopsis plants. Notably, the PmTNL2 promoter activity was highly expressed in shoot apical and floral meristems and was induced strongly with vascular specificity by pathogen infection. Our data has provided a fundamental insight into both expression regulation and putative functions of the PmTNL2 gene in the context of plant growth and development, as well as in responses to environmental stressors. Promoter application as a potential tool for tree improvement was further discussed.

Pyramiding of nine transgenes in maize generates high-level resistance against necrotrophic maize pathogens

Key message Nine transgenes from different categories, viz. plant defense response genes and anti-apoptosis genes, played combined roles in maize to inhibit the necrotrophic pathogens Rhizoctonia solani and Bipolaris maydis. Maize sheath blight and southern corn leaf blight are major global threats to maize production. The management of these necrotrophic pathogens has encountered limited success due to the characteristics of their lifestyle. Here, we presented a transgenic pyramiding breeding strategy to achieve nine different resistance genes integrated in one transgenic maize line to combat different aspects of necrotrophic pathogens. These nine genes, selected from two different categories, plant defense response genes (Chi, Glu, Ace-AMP1, Tlp, Rs-AFP2, ZmPROPEP1 and Pti4), and anti-apoptosis genes (Iap and p35), were successfully transferred into maize and further implicated in resistance against the necrotrophic pathogens Rhizoctonia solani and Bipolaris maydis. Furthermore, the transgenic maize line 910, with high expression levels of the nine integrated genes, was selected from 49 lines. Under greenhouse and field trial conditions, line 910 showed significant resistance against maize sheath blight and southern corn leaf blight diseases. Higher-level resistance was obtained after the pyramiding of more resistance transgenes from different categories that function via different mechanisms. The present study provides a successful strategy for the management of necrotrophic pathogens.

Pyramiding of nine transgenes in maize generates high-level resistance against necrotrophic maize pathogens

Key message Nine transgenes from different categories, viz. plant defense response genes and anti-apoptosis genes, played combined roles in maize to inhibit the necrotrophic pathogens Rhizoctonia solani and Bipolaris maydis. Maize sheath blight and southern corn leaf blight are major global threats to maize production. The management of these necrotrophic pathogens has encountered limited success due to the characteristics of their lifestyle. Here, we presented a transgenic pyramiding breeding strategy to achieve nine different resistance genes integrated in one transgenic maize line to combat different aspects of necrotrophic pathogens. These nine genes, selected from two different categories, plant defense response genes (Chi, Glu, Ace-AMP1, Tlp, Rs-AFP2, ZmPROPEP1 and Pti4), and anti-apoptosis genes (Iap and p35), were successfully transferred into maize and further implicated in resistance against the necrotrophic pathogens Rhizoctonia solani and Bipolaris maydis. Furthermore, the transgenic maize line 910, with high expression levels of the nine integrated genes, was selected from 49 lines. Under greenhouse and field trial conditions, line 910 showed significant resistance against maize sheath blight and southern corn leaf blight diseases. Higher-level resistance was obtained after the pyramiding of more resistance transgenes from different categories that function via different mechanisms. The present study provides a successful strategy for the management of necrotrophic pathogens.

Profiling methyl jasmonate-responsive transcriptome for understanding induced systemic resistance in whitebark pine (Pinus albicaulis)

RNA-seq analysis on whitebark pine needles demonstrated that methyl jasmonate (MeJA)-triggered transcriptome re-programming substantially overlapped with defense responses against insects and fungal pathogens in Pinus species, increasing current knowledge regarding induced systemic resistance (ISR) to pathogens and pests in whitebark pine. Many whitebark pine populations are in steep decline due to high susceptibility to mountain pine beetle and the non-native white pine blister rust (WPBR). Resistance, including induced systemic resistance (ISR), is not well characterized in whitebark pine, narrowing the current options for increasing the success of restoration and breeding programs. Exogenous jasmonates are known to trigger ISR by activating the plant's immune system through regulation of gene expression to produce chemical defense compounds. This study reports profiles of whitebark pine needle transcriptomes, following methyl jasmonate (MeJA) treatment using RNA-seq. A MeJA-responsive transcriptome was de novo assembled and transcriptome profiling identified a set of differentially expressed genes (DEGs), revealing 1422 up- and 999 down-regulated transcripts with at least twofold change (FDR corrected p < 0.05) in needle tissues in response to MeJA application. GO analysis revealed that these DEGs have putative functions in plant defense signalling, transcription regulation, biosyntheses of secondary metabolites, and other biological processes. Lineage-specific expression of defense-related genes was characterized through comparison with MeJA signalling in model plants. In particular, MeJA-triggered transcriptome re-programming substantially overlapped with defense responses against WPBR and insects in related Pinus species, suggesting that MeJA may be used to improve whitebark pine resistance to pathogens/pests. Our study provides new insights into molecular mechanisms and metabolic pathways involved in whitebark pine ISR. DEGs identified in this study can be used as candidates to facilitate identification of genomic variation contributing to host resistance and aid in breeding selection of elite genotypes with better adaptive fitness to environmental stressors in this endangered tree species.