Supplementary Figure 1. Boxplots of the height and diameter at the base of the
stem of inoculated plants (half-sib family 440) and t-test results for the comparison
of these parameters’ means between resistant (res) and susceptible (sus) plants.
(A) Boxplot of height (cm) measurements. (B) Boxplot of diameter at the base of
the stem (mm) measurements. Both measurements were made before inoculations. (C) t-Test results for heights comparison. (D) t-Test results for
diameter comparison. N, number of samples. SD, standard deviation.
Supplementary Figure 2. Heatmaps representing the expression patterns of
genes involved in secondary metabolism. (A) Flavonoid biosynthesis pathway. (B)
Terpenoid biosynthesis pathways, including terpenoid backbone biosynthesis
(Terp. Backbone), monoterpenoid biosynthesis, sesquiterpenoid biosynthesis, and
diterpenoid biosynthesis pathways. The color gradient represents mean
expression levels (logTPM) of each gene for control (C), susceptible (S), and
resistant (R) samples.
Supplementary Figure 3. Heatmaps representing the expression patterns of
genes involved in the synthesis of hydrogen peroxide (A) and response to
oxidative stress (B). The color gradient represents mean expression levels
(logTPM) of each gene for control (C), susceptible (S), and resistant (R) samples.
Supplementary Figure 4. Heatmaps representing the expression patterns of
hormone responsive transcription factors (TFs). Jasmonate responsive TFs
JAZ/Tify (A) and ERF (B), salicylic acid responsive TFs WRKY (C), and abscisic
acid responsive TFs NAC (D). The color gradient represents mean expression
levels (logTPM) of each gene for control (C), susceptible (S), and resistant (R)
samples.
Supplementary Table 1. De novo assembly and P. pinaster reference
transcriptome statistics.
Supplementary Table 2. Genes selected for quantitative RT-qPCR, respective
primer sequences, amplicon size, and annealing temperatures used.
Supplementary Table 3. Genes expressed by Bursaphelenchus xylophilus in
inoculated samples.
Supplementary Table 4. Differential expressed genes in susceptible plants
when compared to controls (log2 fold change ≥ |2|, FDR corrected p-value (padj)
≤ 0.05). InterPro, KEGG, and blastx annotations are presented.
Supplementary Table 5. Differential expressed genes in resistant plants when
compared to controls (log2 fold change ≥ |2|, FDR corrected p-value (padj) ≤
0.05). InterPro, KEGG, and blastx annotations are presented.
Supplementary Table 6. GO terms enriched in the upregulated genes in
susceptible samples when compared with controls, after trimming for redundancy.
Supplementary Table 7. GO terms enriched in the upregulated genes in
resistant samples when compared with controls, after trimming for redundancy.
Supplementary Table 8. Complete list of genes annotated with DRAGO 2 tool.
Supplementary Table 9. Genes used for calculating average log2(fold change)
expression levels for Figures 5, 6.
Data Sheet 1. Fasta file with de novo assembled Pinus pinaster transcripts.