Genome-wide association study identifies SNP markers and putative candidate genes for terpene traits important for Leptocybe invasa resistance in Eucalyptus grandis

Terpenes are an important group of plant specialized metabolites influencing, amongst other functions, defence mechanisms against pests. We used a genome-wide association study to identify single nucleotide polymorphism (SNP) markers and putative candidate genes for terpene traits. We tested 15,387 informative SNP markers derived from genotyping 416 Eucalyptus grandis individuals for association with 3 terpene traits, 1,8-cineole, γ-terpinene, and p-cymene. A multilocus mixed model analysis identified 21 SNP markers for 1,8-cineole on chromosomes 2, 4, 6, 7, 8, 9, 10, and 11, that individually explained 3.0%–8.4% and jointly 42.7% of the phenotypic variation. Association analysis of γ-terpinene found 32 significant SNP markers on chromosomes 1, 2, 4, 5, 6, 9, and 11, explaining 3.4–15.5% and jointly 54.5% of phenotypic variation. For p-cymene, 28 significant SNP markers were identified on chromosomes 1, 2, 3, 5, 6, 7, 10, and 11, explaining 3.4–16.1% of the phenotypic variation and jointly 46.9%. Our results show that variation underlying the 3 terpene traits is influenced by a few minor loci in combination with a few major effect loci, suggesting an oligogenic nature of the traits.

A Genome-Wide Association Study for Resistance to the Insect Pest Leptocybe invasa in Eucalyptus grandis Reveals Genomic Regions and Positional Candidate Defense Genes

The galling insect, Leptocybe invasa, causes significant losses in plantations of various Eucalyptus species and hybrids, threatening its economic viability. We applied a genome-wide association study (GWAS) to identify single-nucleotide polymorphism (SNP) markers associated with resistance to L. invasa. A total of 563 insect-challenged Eucalyptus grandis trees, from 61 half-sib families, were genotyped using the EUChip60K SNP chip, and we identified 15,445 informative SNP markers in the test population. Multi-locus mixed-model (MLMM) analysis identified 35 SNP markers putatively associated with resistance to L. invasa based on four discreet classes of insect damage scores: (0) not infested, (1) infested showing evidence of oviposition but no gall development, (2) infested with galls on leaves, midribs or petioles and (3) stunting and lethal gall formation. MLMM analysis identified three associated genomic regions on chromosomes 3, 7 and 8 jointly explaining 17.6% of the total phenotypic variation. SNP analysis of a validation population of 494 E. grandis trees confirmed seven SNP markers that were also detected in the initial association analysis. Based on transcriptome profiles of resistant and susceptible genotypes from an independent experiment, we identified several putative candidate genes in associated genomic loci including Nucleotide-binding ARC- domain (NB-ARC) and toll-interleukin-1-receptor-Nucleotide binding signal- Leucine rich repeat (TIR-NBS-LRR) genes. Our results suggest that Leptocybe resistance in E. grandis may be influenced by a few large-effect loci in combination with minor effect loci segregating in our test and validation populations.