Exploring the Genetic Basis of Calonectria spp. Resistance in Eucalypts

Selecting high-quality varieties with disease resistance by artificial crossbreeding is the most fundamental way to address the damage caused by Calonectria spp. in eucalypt plantations. However, understanding the mechanism of disease-resistant heterosis occurrence in eucalypts is crucial for successful crossbreeding. Two eucalypt hybrids, the susceptible EC333 (H1522 × unknown) and the resistant EC338 (W1767 × P9060), were screened through infection with Calonectria isolates, a pathogen that causes eucalypt leaf blight. RNA-Seq was performed on the susceptible hybrid, the disease-resistant hybrid, and their parents. The gene differential expression analysis showed that there were 3912 differentially expressed genes between EC333 and EC338, with 1631 up-regulated and 2281 down-regulated genes. The expression trends of the differential gene sets in P9060 and EC338 were similar. However, the expression trend of W1767 was opposite that of EC338. The similarity of the expression and the advantage of stress resistance in E. pellita suggested that genes with significant differences in expression likely relate to disease resistance. A GSEA based on GO annotations revealed that the carbohydrate binding pathway genes were differentially expressed between EC338 and EC333. The gene pathways that were differentially expressed between EC338 and EC333 revealed by the GSEA based on KEGG annotations were the sesquiterpenoid and triterpenoid biosynthesis pathways. The alternative splicing analysis demonstrated that an AS event between EC338 and EC333 occurred in LOC104426602. According to our SNP analysis, EC338 had 626 more high-impact mutation loci than the male parent P9060 and 396 more than the female parent W1767; W1767 had 259 more mutation loci in the downstream region than EC338, while P9060 had 3107 fewer mutation loci in the downstream region than EC338. Additionally, EC338 had 9631 more mutation loci in the exon region than EC333. Modules were found via WGCNA that were strongly and oppositely correlated with EC338 and EC333, such as module MEsaddlebrown, likely associated with leaf blight resistance. The present study provides a detailed explanation of the genetic basis of eucalypt leaf blight resistance, providing the foundation for exploring genes related to this phenomenon.