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

Comparison of ddRADseq and EUChip60K SNP genotyping systems for population genetics and genomic selection in Eucalyptus dunnii (Maiden)

Eucalyptus dunnii is one of the most important Eucalyptus species for short-fiber pulp production in regions where other species of the genus are affected by poor soil and climatic conditions. In this context, E. dunnii holds promise as a resource to address and adapt to the challenges of climate change. Despite its rapid growth and favorable wood properties for solid wood products, the advancement of its improvement remains in its early stages. In this work, we evaluated the performance of two single nucleotide polymorphism, (SNP), genotyping methods for population genetics analysis and Genomic Selection in E. dunnii. Double digest restriction-site associated DNA sequencing (ddRADseq) was compared with the EUChip60K array in 308 individuals from a provenance-progeny trial. The compared SNP set included 8,011 and 19,008 informative SNPs distributed along the 11 chromosomes, respectively. Although the two datasets differed in the percentage of missing data, genome coverage, minor allele frequency and estimated genetic diversity parameters, they revealed a similar genetic structure, showing two subpopulations with little differentiation between them, and low linkage disequilibrium. GS analyses were performed for eleven traits using Genomic Best Linear Unbiased Prediction (GBLUP) and a conventional pedigree-based model (ABLUP). Regardless of the SNP dataset, the predictive ability (PA) of GBLUP was better than that of ABLUP for six traits (Cellulose content, Total and Ethanolic extractives, Total and Klason lignin content and Syringyl and Guaiacyl lignin monomer ratio). When contrasting the SNP datasets used to estimate PAs, the GBLUP-EUChip60K model gave higher and significant PA values for six traits, meanwhile, the values estimated using ddRADseq gave higher values for three other traits. The PAs correlated positively with narrow sense heritabilities, with the highest correlations shown by the ABLUP and GBLUP-EUChip60K. The two genotyping methods, ddRADseq and EUChip60K, are generally comparable for population genetics and genomic prediction, demonstrating the utility of the former when subjected to rigorous SNP filtering. The results of this study provide a basis for future whole-genome studies using ddRADseq in non-model forest species for which SNP arrays have not yet been developed.

Genomic and metabolic profiling of two tomato contrasting cultivars for tolerance to Tuta absoluta

Dissimilar patterns of variants affecting genes involved in response to herbivory, including those leading to difference in VOC production, were identified in tomato lines with contrasting response to Tuta absoluta. Tuta absoluta is one of the most destructive insect pest affecting tomato production, causing important yield losses both in open field and greenhouse. The selection of tolerant varieties to T. absoluta is one of the sustainable approaches to control this invasive leafminer. In this study, the genomic diversity of two tomato varieties, one tolerant and the other susceptible to T. absoluta infestation was explored, allowing us to identify chromosome regions with highly dissimilar pattern. Genes affected by potential functional variants were involved in several processes, including response to herbivory and secondary metabolism. A metabolic analysis for volatile organic compounds (VOCs) was also performed, highlighting a difference in several classes of chemicals in the two genotypes. Taken together, these findings can aid tomato breeding programs aiming to develop tolerant plants to T. absoluta.

Haplotype mining panel for genetic dissection and breeding in Eucalyptus

To improve our understanding of genetic mechanisms underlying complex traits in plants, a comprehensive analysis of gene variants is required. Eucalyptus is an important forest plantation genus that is highly outbred. Trait dissection and molecular breeding in eucalypts currently relies on biallelic single-nucleotide polymorphism (SNP) markers. These markers fail to capture the large amount of haplotype diversity in these species, and thus multi-allelic markers are required. We aimed to develop a gene-based haplotype mining panel for Eucalyptus species. We generated 17 999 oligonucleotide probe sets for targeted sequencing of selected regions of 6293 genes implicated in growth and wood properties, pest and disease resistance, and abiotic stress responses. We identified and phased 195 834 SNPs using a read-based phasing approach to reveal SNP-based haplotypes. A total of 8915 target regions (at 4637 gene loci) passed tests for Mendelian inheritance. We evaluated the haplotype panel in four Eucalyptus species (E. grandis, E. urophylla, E. dunnii and E. nitens) to determine its ability to capture diversity across eucalypt species. This revealed an average of 3.13-4.52 haplotypes per target region in each species, and 33.36% of the identified haplotypes were shared by at least two species. This haplotype mining panel will enable the analysis of haplotype diversity within and between species, and provide multi-allelic markers that can be used for genome-wide association studies and gene-based breeding approaches.

Genomic consequences of artificial selection during early domestication of a wood fibre crop

From its origins in Australia, Eucalyptus grandis has spread to every continent, except Antarctica, as a wood crop. It has been cultivated and bred for over 100 yr in places such as South Africa. Unlike most annual crops and fruit trees, domestication of E. grandis is still in its infancy, representing a unique opportunity to interrogate the genomic consequences of artificial selection early in the domestication process. To determine how a century of artificial selection has changed the genome of E. grandis, we generated single nucleotide polymorphism genotypes for 1080 individuals from three advanced South African breeding programmes using the EUChip60K chip, and investigated population structure and genome-wide differentiation patterns relative to wild progenitors. Breeding and wild populations appeared genetically distinct. We found genomic evidence of evolutionary processes known to have occurred in other plant domesticates, including interspecific introgression and intraspecific infusion from wild material. Furthermore, we found genomic regions with increased linkage disequilibrium and genetic differentiation, putatively representing early soft sweeps of selection. This is, to our knowledge, the first study of genomic signatures of domestication in a timber species looking beyond the first few generations of cultivation. Our findings highlight the importance of intra- and interspecific hybridization during early domestication.

Genetic control of the operculum and capsule morphology of Eucalyptus globulus

BACKGROUND AND AIMS

The petaline operculum that covers the inner whorls until anthesis and the woody capsule that develops after fertilization are reproductive structures of eucalypts that protect the flower and seeds. Although they are distinct organs, they both develop from flower buds and this common ontogeny suggests shared genetic control. In Eucalyptus globulus their morphology is variable and we aimed to identify the quantitative trait loci (QTL) underlying this variation and determine whether there is common genetic control of these ecologically and taxonomically important reproductive structures.

METHODS

Samples of opercula and capsules were collected from 206 trees that belong to a large outcrossed F2E. globulus mapping population. The morphological variation in these structures was characterized by measuring six operculum and five capsule traits. QTL analysis was performed using these data and a linkage map consisting of 480 markers.

KEY RESULTS

A total of 27 QTL were detected for operculum traits and 28 for capsule traits, with the logarithm of odds ranging from 2.8 to 11.8. There were many co-located QTL associated with operculum or capsule traits, generally reflecting allometric relationships. A key finding was five genomic regions where co-located QTL affected both operculum and capsule morphology, and the overall trend for these QTL was to affect elongation of both organs. Some of these QTL appear to have a significant effect on the phenotype, with the strongest QTL explaining 26.4 % of the variation in operculum shape and 16.4 % in capsule shape. Flower bud measurements suggest the expression of these QTL starts during bud development. Several candidate genes were found associated with the QTL and their putative function is discussed.

CONCLUSIONS

Variation in both operculum and capsule traits in E. globulus is under strong genetic control. Our results suggest that these reproductive structures share a common genetic pathway during flower bud development.

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.

Cell death and changes in primary metabolism: the onset of defence in Eucalyptus in the war against Leptocybe invasa

BACKGROUND

Here, we investigated changes in primary metabolism and cell death around oviposition sites in two hybrid clones of Eucalyptus with different degrees of resistance to Leptocybe invasa Fisher & La Salle (Hymenoptera: Eulophidae), as well as tolerance to water deficiency.

RESULTS

We showed that apices of the resistant clone with oviposition had a higher content of amino acids, organic acids and the compound putrescine compared with those of the susceptible clone with oviposition. By contrast, apices of the resistant clone with oviposition had lower sugar and pyruvate organic acid content than those of the susceptible clone with oviposition. Small areas of necrosis were induced around the oviposition sites in the stem apices of Eucalyptus 24 h after infestation. The resistant clone developed larger necrotic areas that showed progressive increases 24-72 h after infestation compared with the susceptible clone, in which cell death was significantly lower and no changes were observed in necrotic area over time. Thus, the programmed death of cells around the egg, modulated by several amino acids, is likely the first defence response of Eucalyptus against L. invasa.

CONCLUSION

Our results serve as the basis for the early identification of key metabolites produced in plants in defence against galling insects. © 2022 Society of Chemical Industry.

Fungal Communities of Eucalyptus grandis Leaves Are Influenced by the Insect Pest Leptocybe invasa

Fungal communities in above-ground tree tissues are hyperdiverse and are influenced by biotic interactions with other organisms living in or on these tissues. These biotic interactions are, however, still poorly understood. In this study, we aimed to understand how insect-associated gall formation on Eucalyptus foliage correlates with the diversity of foliar fungal communities in surrounding healthy leaf tissue, as well as the co-occurrence patterns among the members of the fungal community. We used ITS metabarcoding to characterise the foliar fungal communities of 179 individual E. grandis trees. These trees were assigned to infestation levels of the wasp Leptocybe invasa (Eulophidae: Hymenoptera), which causes gall formation on shoot tips and leaves of its host. Fungal community networks were calculated using a Pearson correlation coefficient. The composition and diversity of fungal communities were influenced by the severity of L. invasa infestations. We identified potential Eucalyptus pathogens with high sequence abundance at all disease severity levels, but network analysis indicated that the co-occurrence of potential pathogens between no to mild and medium to heavy infestation differed significantly. A better understanding of microbial interactions, especially the role of pathogens, can be useful for controlling disease- and beneficial host-associated microbial communities.

Molecular ecology of plant volatiles in interactions with insect herbivores

Plant-derived volatile organic compounds (VOCs) play pivotal roles in interactions with insect herbivores. Individual VOCs can be directly toxic or deterrent, serve as signal molecules to attract natural enemies, and/or be perceived by distal plant tissues as a priming signal to prepare for expected herbivory. Environmental conditions, as well as the specific plant-insect interaction being investigated, strongly influence the observed functions of VOC blends. The complexity of plant-insect chemical communication via VOCs is further enriched by the sophisticated molecular perception mechanisms of insects, which can respond to one or more VOCs and thereby influence insect behavior in a manner that has yet to be fully elucidated. Despite numerous gaps in the current understanding of VOC-mediated plant-insect interactions, successful pest management strategies such as push-pull systems, synthetic odorant traps, and crop cultivars with modified VOC profiles have been developed to supplement chemical pesticide applications and enable more sustainable agricultural practices. Future studies in this field would benefit from examining the responses of both plants and insects in the same experiment to gain a more complete view of these interactive systems. Furthermore, a molecular evolutionary study of key genetic elements of the ecological interaction phenotypes could provide new insights into VOC-mediated plant communication with insect herbivores.

Genetic basis of growth, spring phenology, and susceptibility to biotic stressors in maritime pine

Forest ecosystems are increasingly challenged by extreme events, for example, drought, storms, pest attacks, and fungal pathogen outbreaks, causing severe ecological and economic losses. Understanding the genetic basis of adaptive traits in tree species is of key importance to preserve forest ecosystems, as genetic variation in a trait (i.e., heritability) determines its potential for human-mediated or evolutionary change. Maritime pine (Pinus pinaster Aiton), a conifer widely distributed in southwestern Europe and northwestern Africa, grows under contrasted environmental conditions promoting local adaptation. Genetic variation at adaptive phenotypes, including height, spring phenology, and susceptibility to two fungal pathogens (Diplodia sapinea and Armillaria ostoyae) and an insect pest (Thaumetopoea pityocampa), was assessed in a range-wide clonal common garden of maritime pine. Broad-sense heritability was significant for height (0.219), spring phenology (0.165-0.310), and pathogen susceptibility (necrosis length caused by D. sapinea, 0.152; and by A. ostoyae, 0.021, measured on inoculated, excised branches under controlled conditions), but not for pine processionary moth incidence in the common garden. The correlations of trait variation among populations revealed contrasting trends for pathogen susceptibility to D. sapinea and A. ostoyae with respect to height. Taller trees showed longer necrosis length caused by D. sapinea while shorter trees were more affected by A. ostoyae. Moreover, maritime pine populations from areas with high summer temperatures and frequent droughts were less susceptible to D. sapinea but more susceptible to A. ostoyae. Finally, an association study using 4227 genome-wide SNPs revealed several loci significantly associated with each trait (range of 3-26), including a possibly disease-induced translation initiation factor, eIF-5, associated with needle discoloration caused by D. sapinea. This study provides important insights to develop genetic conservation and breeding strategies integrating species responses to biotic stressors.

Differential development times of galls induced by Leptocybe invasa (Hymenoptera: Eulophidae) reveal differences in susceptibility between two Eucalyptus clones

BACKGROUND

The significance of morphological responses of hosts on susceptibility against gall-inducing insects is relatively unknown, especially in planted forests. Here, we investigate the temporal morphological responses (gall development) induced by the invasive gall wasp Leptocybe invasa and the subsequent insect development in two Eucalyptus clones.

RESULTS

Our results identified a novel stage of gall development, not previously reported, termed here 'brownish ring'. In both hosts similar gall development stages were observed. Although L. invasa oviposited in both clones, comparison of external morphological traits of galls over time revealed a differential response in the number of galls between clones. Comparison of the developmental time of each gall and insect stage between clones suggests that plant defense mechanisms against L. invasa are activated shortly after oviposition by the wasp, yet before gall formation.

CONCLUSION

Gall number is an important parameter that should be used to measure host susceptibility among Eucalyptus clones. To the best of our knowledge, this is the first study showing differential morphological responses induced by a galling insect, even before gall formation, revealing differences in susceptibility between different plant hosts. These findings provide insight into the use of early stages of gall formation by L. invasa to prevent invasion and establishment of this pest.

Genomic Breeding for Diameter Growth and Tolerance to Leptocybe Gall Wasp and Botryosphaeria/Teratosphaeria Fungal Disease Complex in Eucalyptus grandis

Eucalyptus grandis is one of the most important species for hardwood plantation forestry around the world. At present, its commercial deployment is in decline because of pests and pathogens such as Leptocybe invasa gall wasp (Lepto), and often co-occurring fungal stem diseases such as Botryosphaeria dothidea and Teratosphaeria zuluensis (BotryoTera). This study analyzed Lepto, BotryoTera, and stem diameter growth in an E. grandis multi-environmental, genetic trial. The study was established in three subtropical environments. Diameter growth and BotryoTera incidence scores were assessed on 3,334 trees, and Lepto incidence was assessed on 4,463 trees from 95 half-sib families. Using the Eucalyptus EUChip60K SNP chip, a subset of 964 trees from 93 half-sib families were genotyped with 14,347 informative SNP markers. We employed single-step genomic BLUP (ssGBLUP) to estimate genetic parameters in the genetic trial. Diameter and Lepto tolerance showed a positive genetic correlation (0.78), while BotryoTera tolerance had a negative genetic correlation with diameter growth (-0.38). The expected genetic gains for diameter growth and Lepto and BotryoTera tolerance were 12.4, 10, and -3.4%, respectively. We propose a genomic selection breeding strategy for E. grandis that addresses some of the present population structure problems.