The Transcriptome and Terpene Profile of Eucalyptus grandis Reveals Mechanisms of Defense Against the Insect Pest, Leptocybe invasa

A chromosome-scale genome of Rhus chinensis Mill. provides new insights into plant-insect interaction and gallotannins biosynthesis

Rhus chinensis Mill., an economically valuable Anacardiaceae species, is parasitized by the galling aphid Schlechtendalia chinensis, resulting in the formation of the Chinese gallnut (CG). Here, we report a chromosomal-level genome assembly of R. chinensis, with a total size of 389.40 Mb and scaffold N50 of 23.02 Mb. Comparative genomic and transcriptome analysis revealed that the enhanced structure of CG and nutritional metabolism contribute to improving the adaptability of R. chinensis to S. chinensis by supporting CG and galling aphid growth. CG was observed to be abundant in hydrolysable tannins (HT), particularly gallotannin and its isomers. Tandem repeat clusters of dehydroquinate dehydratase/shikimate dehydrogenase (DQD/SDH) and serine carboxypeptidase-like (SCPL) and their homologs involved in HT production were determined as specific to HT-rich species. The functional differentiation of DQD/SDH tandem duplicate genes and the significant contraction in the phenylalanine ammonia-lyase (PAL) gene family contributed to the accumulation of gallic acid and HT while minimizing the production of shikimic acid, flavonoids, and condensed tannins in CG. Furthermore, we identified one UDP glucosyltransferase (UGT84A), three carboxylesterase (CXE), and six SCPL genes from conserved tandem repeat clusters that are involved in gallotannin biosynthesis and hydrolysis in CG. We then constructed a regulatory network of these genes based on co-expression and transcription factor motif analysis. Our findings provide a genomic resource for the exploration of the underlying mechanisms of plant-galling insect interaction and highlight the importance of the functional divergence of tandem duplicate genes in the accumulation of secondary metabolites.

Anatomical and Metabolome Features of Haloxylon aphyllum and Haloxylon persicum Elucidate the Resilience against Gall-Forming Insects

Globally, gall-forming insects significantly contribute to the degradation of desert ecosystems. Recent studies have demonstrated that Haloxylon persicum suffers less damage from gall-formers compared to Haloxylon aphyllum. However, the mechanisms driving the long-term metabolic responses of these species to gall-forming biotic stress in their natural environment remain unclear. The current study comparatively analyzes the anatomical features and metabolomic changes in H. aphyllum and H. persicum damaged by gall-forming insects. This research aimed to uncover potential metabolic tolerance mechanisms through GC-MS analysis. The study findings indicate that gall-forming insects cause a reduction in nearly all the anatomical structures of Haloxylon shoots, with the effects being less severe in H. persicum than in H. aphyllum. Thus, the metabolic pathways responsible for the biosynthesis of biologically active substances that enhance resistance to gall inducers were different, specifically in H. aphyllum-the biosynthesis of fatty acids (+their derivatives) and γ-tocopherol (vitamin E) and H. persicum-the biosynthesis of fatty acids (+their derivatives), dialkyl ethers, carbohydrates (+their derivatives), aromatic acid derivatives, phytosterols, γ-tocopherol (vitamin E), phenols, and terpenoids. The results suggest that the modulation of metabolic pathways under biotic stress plays a crucial role in the enhanced survival and growth of H. persicum.

Genome-wide identification and expression profile analysis of metal tolerance protein gene family in Eucalyptus grandis under metal stresses

Metal tolerance proteins (MTPs) as Me²⁺/H⁺(K⁺) antiporters participate in the transport of divalent cations, leading to heavy metal stress resistance and mineral utilization in plants. In the present study, to obtain better knowledge of the biological functions of the MTPs family, 20 potential EgMTPs genes were identified in Eucalyptus grandis and classified into seven groups belonging to three cation diffusion facilitator groups (Mn-CDFs, Zn/Fe-CDFs, and Zn-CDFs) and seven groups. EgMTP-encoded amino acids ranged from 315 to 884, and most of them contained 4-6 recognized transmembrane domains and were clearly prognosticated to localize into the cell vacuole. Almost all EgMTP genes experienced gene duplication events, in which some might be uniformly distributed in the genome. The numbers of cation efflux and the zinc transporter dimerization domain were highest in EgMTP proteins. The promoter regions of EgMTP genes have different cis-regulatory elements, indicating that the transcription rate of EgMTP genes can be a controlled response to different stimuli in multiple pathways. Our findings provide accurate perception on the role of the predicted miRNAs and the presence of SSR marker in the Eucalyptus genome and clarify their functions in metal tolerance regulation and marker-assisted selection, respectively. Gene expression profiling based on previous RNA-seq data indicates a probable function for EgMTP genes during development and responses to biotic stress. Additionally, the upregulation of EgMTP6, EgMTP5, and EgMTP11.1 to excess Cd²⁺ and Cu²⁺ exposure might be responsible for metal translocation from roots to leaves.

Transcriptome-Based Construction of the Gibberellin Metabolism and Signaling Pathways in Eucalyptus grandis × E. urophylla, and Functional Characterization of GA20ox and GA2ox in Regulating Plant Development and Abiotic Stress Adaptations

Gibberellins (GAs) are the key regulators controlling plant growth, wood production and the stress responses in perennial woody plants. The role of GA in regulating the above-mentioned processes in Eucalyptus remain largely unclear. There is still a lack of systematic identification and functional characterization of GA-related genes in Eucalyptus. In this study, a total of 59,948 expressed genes were identified from the major vegetative tissues of the E. grandis × E. urophylla using transcriptome sequencing. Then, the key gene families in each step of GA biosynthesis, degradation and signaling were investigated and compared with those of Arabidopsis, rice, and Populus. The expression profile generated using Real-time quantitative PCR showed that most of these genes exhibited diverse expression patterns in different vegetative organs and in response to abiotic stresses. Furthermore, we selectively overexpressed EguGA20ox1, EguGA20ox2 and EguGA2ox1 in both Arabidopsis and Eucalyptus via Agrobacterium tumefaciens or A. rhizogenes-mediated transformation. Though both Arabidopsis EguGA20ox1- and EguGA20ox2-overexpressing (OE) lines exhibited better vegetative growth performance, they were more sensitive to abiotic stress, unlike EguGA2ox1-OE plants, which exhibited enhanced stress resistance. Moreover, overexpression of EguGA20ox in Eucalyptus roots caused significantly accelerated hairy root initiation and elongation and improved root xylem differentiation. Our study provided a comprehensive and systematic study of the genes of the GA metabolism and signaling and identified the role of GA20ox and GA2ox in regulating plant growth, stress tolerance, and xylem development in Eucalyptus; this could benefit molecular breeding for obtaining high-yield and stress-resistant Eucalyptus cultivars.

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.

Transcriptome profiling of Toona ciliata young stems in response to Hypsipyla robusta Moore

Toona ciliata is a traditional woody plant that can be used as a medicinal material in China. The extracts of its roots, stems, leaves, and flowers all have a wide range of bioactive compounds. However, T. ciliata has been facing an unresolved pest problem caused by Hypsipyla robusta Moore (HRM), which seriously affects its growth and development. In this study, the expression level of TcMYB3 gene reached the maximum (28-fold) at 12 h and transcriptome sequencing of young stems eaten by HRM for 0, 3, 12, and 21 h were performed. A large number of differentially expressed genes (DEGs) were identified including jointly up-regulated genes (263) and down-regulated genes (378). JA synthesis and signaling transduction, terpene biosynthesis, and MAPKs signaling pathway were analyzed in depth and found that TcOPR3, TcJAR1, TcJAZs, and TcTPS9 genes possessed anti-insect potential. Moreover, MYB and ERF transcription factor (TF) families were significantly strengthened to the point that they may participate in induced defense mechanisms in T. ciliata. These data not only provide insights into the molecular mechanisms in resistance of T. ciliata to HRM but also helps to explore the new biocontrol strategies against insects in eco-friendly woody plants.

The survival, development, and reproduction of Gonipterus platensis (Coleoptera: Curculionidae) on the main Eucalyptus (Myrtaceae) genotypes planted in Brazil

Background

Gonipterus platensis Marelli (Coleoptera: Curculionidae) is the main defoliating beetle of Eucalyptus L'Hér. (Myrtaceae) plants worldwide. The suitability of Eucalyptus to this pest varies among host plant genotypes. The objective of this study was to evaluate the development, reproduction, and survival of G. platensis on Eucalyptus species and hybrids to assess their suitability to this insect pest in Brazil.

Methods

The survival, development, and reproduction parameters were evaluated with G. platensis feeding leaves of Eucalyptus camaldulensis Dehnh., Eucalyptus grandis W. Hill., Eucalyptus urophylla S.T. Blake and on the hybrids of E. grandis ×E. urophylla 'H13' and 'VR3748' in the laboratory.

Results

The duration of the larval stage of G. platensis was shorter on E. urophylla. The pupal stage and the period from larva to adult were equally shorter on E. urophylla and E. camaldulensis. The viability of instars of this insect was low on both E. grandis and E. camaldulensis. The complete lifespan, oviposition period and reproduction parameters of G. platensis were greater on E. urophylla, lower on E. camaldulensis and E. grandis, and intermediate on both hybrids tested.

Synthesis

Eucalyptus urophylla is the most suitable host for G. platensis survival, development, and reproduction, while E. grandis and E. camaldulensis are the least suitable.

Recent advances in proteomics and metabolomics in plants

Over the past decade, systems biology and plant-omics have increasingly become the main stream in plant biology research. New developments in mass spectrometry and bioinformatics tools, and methodological schema to integrate multi-omics data have leveraged recent advances in proteomics and metabolomics. These progresses are driving a rapid evolution in the field of plant research, greatly facilitating our understanding of the mechanistic aspects of plant metabolisms and the interactions of plants with their external environment. Here, we review the recent progresses in MS-based proteomics and metabolomics tools and workflows with a special focus on their applications to plant biology research using several case studies related to mechanistic understanding of stress response, gene/protein function characterization, metabolic and signaling pathways exploration, and natural product discovery. We also present a projection concerning future perspectives in MS-based proteomics and metabolomics development including their applications to and challenges for system biology. This review is intended to provide readers with an overview of how advanced MS technology, and integrated application of proteomics and metabolomics can be used to advance plant system biology research.

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.

The Diesel Tree Sindora glabra Genome Provides Insights Into the Evolution of Oleoresin Biosynthesis

Sindora glabra is an economically important tree that produces abundant oleoresin in the trunk. Here, we present a high-quality chromosome-scale assembly of S. glabra genome by combining Illumina HiSeq, Pacific Biosciences sequencing, and Hi-C technologies. The size of S. glabra genome was 1.11 Gb, with a contig N50 of 1.27 Mb and 31,944 predicted genes. This is the first sequenced genome of the subfamily Caesalpinioideae. As a sister taxon to Papilionoideae, S. glabra underwent an ancient genome triplication shared by core eudicots and further whole-genome duplication shared by early-legume in the last 73.3 million years. S. glabra harbors specific genes and expanded genes largely involved in stress responses and biosynthesis of secondary metabolites. Moreover, 59 terpene backbone biosynthesis genes and 64 terpene synthase genes were identified, which together with co-expressed transcription factors could contribute to the diversity and specificity of terpene compounds and high terpene content in S. glabra stem. In addition, 63 disease resistance NBS-LRR genes were found to be unique in S. glabra genome and their expression levels were correlated with the accumulation of terpene profiles, suggesting potential defense function of terpenes in S. glabra. These together provide new resources for understanding genome evolution and oleoresin production.

qtlXplorer: an online systems genetics browser in the Eucalyptus Genome Integrative Explorer (EucGenIE)

BACKGROUND

Affordable high-throughput DNA and RNA sequencing technologies are allowing genomic analysis of plant and animal populations and as a result empowering new systems genetics approaches to study complex traits. The availability of intuitive tools to browse and analyze the resulting large-scale genetic and genomic datasets remain a significant challenge. Furthermore, these integrative genomics approaches require innovative methods to dissect the flow and interconnectedness of biological information underlying complex trait variation. The Plant Genome Integrative Explorer (PlantGenIE.org) is a multi-species database and domain that houses online tools for model and woody plant species including Eucalyptus. Since the Eucalyptus Genome Integrative Explorer (EucGenIE) is integrated within PlantGenIE, it shares genome and expression analysis tools previously implemented within the various subdomains (ConGenIE, PopGenIE and AtGenIE). Despite the success in setting up integrative genomics databases, online tools for systems genetics modelling and high-resolution dissection of complex trait variation in plant populations have been lacking.

RESULTS

We have developed qtlXplorer ( https://eucgenie.org/QTLXplorer ) for visualizing and exploring systems genetics data from genome-wide association studies including quantitative trait loci (QTLs) and expression-based QTL (eQTL) associations. This module allows users to, for example, find co-located QTLs and eQTLs using an interactive version of Circos, or explore underlying genes using JBrowse. It provides users with a means to build systems genetics models and generate hypotheses from large-scale population genomics data. We also substantially upgraded the EucGenIE resource and show how it enables users to combine genomics and systems genetics approaches to discover candidate genes involved in biotic stress responses and wood formation by focusing on two multigene families, laccases and peroxidases.

CONCLUSIONS

qtlXplorer adds a new dimension, population genomics, to the EucGenIE and PlantGenIE environment. The resource will be of interest to researchers and molecular breeders working in Eucalyptus and other woody plant species. It provides an example of how systems genetics data can be integrated with functional genetics data to provide biological insight and formulate hypotheses. Importantly, integration within PlantGenIE enables novel comparative genomics analyses to be performed from population-scale data.

Characterization of terpene biosynthesis in Melaleuca quinquenervia and ecological consequences of terpene accumulation during myrtle rust infection

Plants use a wide array of secondary metabolites including terpenes as defense against herbivore and pathogen attack, which can be constitutively expressed or induced. Here, we investigated aspects of the chemical and molecular basis of resistance against the exotic rust fungus Austropuccinia psidii in Melaleuca quinquenervia, with a focus on terpenes. Foliar terpenes of resistant and susceptible plants were quantified, and we assessed whether chemotypic variation contributed to resistance to infection by A. psidii. We found that chemotypes did not contribute to the resistance and susceptibility of M. quinquenervia. However, in one of the chemotypes (Chemotype 2), susceptible plants showed higher concentrations of several terpenes including α-pinene, limonene, 1,8-cineole, and viridiflorol compared with resistant plants. Transcriptome profiling of these plants showed that several TPS genes were strongly induced in response to infection by A. psidii. Functional characterization of these TPS showed them to be mono- and sesquiterpene synthases producing compounds including 1,8-cineole, β-caryophyllene, viridiflorol and nerolidol. The expression of these TPS genes correlated with metabolite data in a susceptible plant. These results suggest the complexity of resistance mechanism regulated by M. quinquenervia and that modulation of terpenes may be one of the components that contribute to resistance against A. psidii.

Insect egg-induced physiological changes and transcriptional reprogramming leading to gall formation

Gall-inducing insects and their hosts present some of the most intricate plant-herbivore interactions. Oviposition on the host is often the first cue of future herbivory and events at this early time point can affect later life stages. Many gallers are devastating plant pests, yet little information regarding the plant-insect molecular interplay exists, particularly following egg deposition. We studied the physiological and transcriptional responses of Eucalyptus following oviposition by the gall-inducing wasp, Leptocybe invasa, to explore potential mechanisms governing defence responses and gall development. RNA sequencing and microscopy were used to explore a susceptible Eucalyptus-L. invasa interaction. Infested and control material was compared over time (1-3, 7 and 90 days post oviposition) to examine the transcriptional and morphological changes. Oviposition induces accumulation of reactive oxygen species and phenolics which is reflected in the transcriptome analysis. Gene expression supports phytohormones and 10 transcription factor subfamilies as key regulators. The egg and oviposition fluid stimulate cell division resulting in gall development. Eucalyptus responses to oviposition are apparent within 24 hr. Putative defences include the oxidative burst and barrier reinforcement. However, egg and oviposition fluid stimuli may redirect these responses towards gall development.

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.

Tree health in South Africa: Retrospect and prospect

South Africa is a country with very limited natural forest cover. Consequently, the timber and fibre needs of the country cannot be provided for from indigenous forest. It is largely for this reason that South Africa initially developed a highly productive plantation forest industry, which today makes a substantial contribution to the local economy. These plantations are based on non-native species of Eucalyptus, Pinus and Australian Acacia. In the early years of establishment, South African plantations were relatively free of pest and pathogen problems. But, over time, an increasing number of insects, fungi and bacteria have emerged as serious threats to the sustainability of the forestry industry. Numerous native pests and pathogens, especially insects, have adapted to these introduced tree species to cause damage or disease. The problem is compounded by the accidental introduction of non-native pests and pathogens, and this has been at a rapidly increasing rate over the past three decades. Some of these introduced pests and pathogens also threaten the fitness and even the survival of many indigenous South African tree species. Fortunately, South Africa has developed an impressive knowledge base and range of integrated management options to deal with these problems. This development was first driven by government programmes, and in more recent years by public–private partnerships between industry, universities and government. It is clear from the pattern of emergence of pests and pathogens in recent years that South Africa will deal with an increasing number of these problems and a continuously changing tree health environment. This requires robust investment in both quarantine and mitigation mechanisms to protect the country’s biodiversity as well as to ensure the sustainability of its wood and fibre industries.

Comparative Transcriptome Analysis Reveals the Mechanism Related to Fluazinam Stress of Panonychus citri (Acarina: Tetranychidae)

Simple Summary

The citrus red mite, Panonychus citri, is an important pest that causes serious citrus production losses in China. The insecticide fluazinam has a good control effect on the pest mites; however, its mechanism of action on mites remains unclear. In this study, we analyzed the transcriptomic sequencing and differential expression genes in P. citri treated with fluazinam, and identified some of the genes potential involved in detoxification metabolism related with the fluazinam exposure. Evaluating the efficacy of fluazinam, and analyzing the transcriptome data of P. citri under fluazinam stress, potentially provide a new agent for prevention and control of P. citri, and also preliminary research results for exploring the mechanism of action of fluazinam on P. citri. Given the up-regulated expression levels of genes for Mn-superoxide dismutase and catalase, we speculate that they play an important role in fluazinam-stress action on P. citri.

Abstract

The use of a large number of chemical acaricides to control these pest mites has led to an increasing problem of pesticide resistance, which has always been the difficulty in integrated pest management (IPM). Fluazinam has a good control effect on Panonychus citri, the serious pest on citrus; however, we only know the mechanism of action of fluazinam as a fungicide and its mechanism of action on mites remains unclear. Through analysis using Illumina high-throughput transcriptomic sequencing and differential expression genes in P. citri treated with fluazinam, 59 cytochrome P450 genes, 23 glutathione s-transferase genes, five carboxylate esterase genes, 11 superoxide dismutase genes and 15 catalase genes were identified. The Gene Ontology enrichment and the enrichment of KEGG results showed that the treatment were enrichment for redox enzyme pathways. Evaluating the efficacy of fluazinam, and analyzing the transcriptome data of P. citri under fluazinam stress, potentially provide a new agent for prevention and control of P. citri, and also preliminary research results for exploring the mechanism of action of fluazinam on P. citri. Given the up-regulated expression levels of genes for Mn-superoxide dismutase and catalase, we speculate that they play an important role in fluazinam-stress action on P. citri.

Genomic evidence of introgression and adaptation in a model subtropical tree species, Eucalyptus grandis

The genetic consequences of adaptation to changing environments can be deciphered using population genomics, which may help predict species' responses to global climate change. Towards this, we used genome-wide SNP marker analysis to determine population structure and patterns of genetic differentiation in terms of neutral and adaptive genetic variation in the natural range of Eucalyptus grandis, a widely cultivated subtropical and temperate species, serving as genomic reference for the genus. We analysed introgression patterns at subchromosomal resolution using a modified ancestry mapping approach and identified provenances with extensive interspecific introgression in response to increased aridity. Furthermore, we describe potentially adaptive genetic variation as explained by environment-associated SNP markers, which also led to the discovery of what is likely a large structural variant. Finally, we show that genes linked to these markers are enriched for biotic and abiotic stress responses.

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.

The interplay between microRNA and alternative splicing of linear and circular RNAs in eleven plant species

MOTIVATION

MicroRNA (miRNA) and alternative splicing (AS)-mediated post-transcriptional regulation has been extensively studied in most eukaryotes. However, the interplay between AS and miRNAs has not been explored in plants. To our knowledge, the overall profile of miRNA target sites in circular RNAs (circRNA) generated by alternative back splicing has never been reported previously. To address the challenge, we identified miRNA target sites located in alternatively spliced regions of the linear and circular splice isoforms using the up-to-date single-molecule real-time (SMRT) isoform sequencing (Iso-Seq) and Illumina sequencing data in eleven plant species.

RESULTS

In total, we identified 399 401 and 114 574 AS events from linear and circular RNAs, respectively. Among them, there were 64 781 and 41 146 miRNA target sites located in linear and circular AS region, respectively. In addition, we found 38 913 circRNAs to be overlapping with 45 648 AS events of its own parent isoforms, suggesting circRNA regulation of AS of linear RNAs by forming R-loop with the genomic locus. Here, we present a comprehensive database of miRNA targets in alternatively spliced linear and circRNAs (ASmiR) and a web server for deposition and identification of miRNA target sites located in the alternatively spliced region of linear and circular RNAs. This database is accompanied by an easy-to-use web query interface for meaningful downstream analysis. Plant research community can submit user-defined datasets to the web service to search AS regions harboring small RNA target sites. In conclusion, this study provides an unprecedented resource to understand regulatory relationships between miRNAs and AS in both gymnosperms and angiosperms.

AVAILABILITY AND IMPLEMENTATION

The readily accessible database and web-based tools are available at http://forestry.fafu.edu.cn/bioinfor/db/ASmiR.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

eCALIBRATOR: A Comparative Tool to Identify Key Genes and Pathways for Eucalyptus Defense Against Biotic Stressors

Many pests and pathogens threaten Eucalyptus plantations. The study of defense responses in this economically important wood and fiber crop enables the discovery of novel pathways and genes, which may be adopted to improve resistance. Various functional genomics experiments have been conducted in Eucalyptus-biotic stress interactions following the availability of the Eucalyptus grandis genome, however, comparisons between these studies were limited largely due to a lack of comparative tools. To this end, we developed eCALIBRATOR http://ecalibrator.bi.up.ac.za, a tool for the comparison of Eucalyptus biotic stress interaction. The tool, which is not limited to Eucalyptus, allows the comparison of various datasets, provides a visual output in the form of Venn diagrams and clustering and extraction of lists for gene ontology enrichment analyses. We also demonstrate the usefulness of the tool in revealing pathways and key gene targets to further functionally characterize. We identified 708 differentially expressed E. grandis genes in common among responses to the insect pest Leptocybe invasa, oomycete pathogen Phytophthora cinnamomi and fungus Chrysoporthe austroafricana. Within this set of genes, one of the Gene Ontology terms enriched was "response to organonitrogen compound," with NITRATE TRANSPORTER 2.5 (NRT2.5) being a key gene, up-regulated under susceptible interactions and down-regulated under resistant interactions. Although previous functional genetics studies in Arabidopsis thaliana support a role in nitrate acquisition and remobilization under long-term nitrate starvation, the importance of NRT2.5 in plant defense is unclear. The T-DNA mutants of AtNRT2.5 were more resistant to Pseudomonas syringae pv. tomato pv tomato DC3000 inoculation than the wild-type counterpart, supporting a direct role for NRT2.5 in plant defense. Future studies will focus on characterizing the Eucalyptus ortholog of NRT2.5.

Foliage inoculation by Burkholderia vietnamiensis CBMB40 antagonizes methyl jasmonate-mediated stress in Eucalyptus grandis

Methyl jasmonate (MeJA) is widely used as a model chemical to study hypersensitive responses to biotic stress impacts in plants. Elevated levels of methyl jasmonate induce jasmonate-dependent defense responses, associated with a decline in primary metabolism and enhancement of secondary metabolism of plants. However, there is no information of how stress resistance of plants, and accordingly the sensitivity to exogenous MeJA can be decreased by endophytic plant growth promoting rhizobacteria (PGPR) harboring ACC (1-aminocyclopropane-1-carboxylate) deaminase. In this study, we estimated stress alleviating potential of endophytic PGPR against MeJA-induced plant perturbations through assessing photosynthetic traits and stress volatile emissions. We used mild (5 mM) to severe (20 mM) MeJA and endophytic plant growth promoting rhizobacteria Burkholderia vietnamiensis CBMB40 and studied how MeJA and B. vietnamiensis treatments influenced temporal changes in photosynthetic characteristics and stress volatile emissions. Separate application of MeJA markedly decreased photosynthetic characteristics and increased lipoxygenase pathway (LOX) volatiles, volatile isoprenoids, saturated aldehydes, lightweight oxygenated compounds (LOC), geranyl-geranyl diphosphate pathway (GGDP) volatiles, and benzenoids. However, MeJA-treated leaves inoculated by endophytic bacteria B. vietnamiensis had substantially increased photosynthetic characteristics and decreased emissions of LOX, volatile isoprenoids and other stress volatiles compared with non-inoculated MeJA treatments, especially at later stages of recovery. In addition, analysis of leaf terpenoid contents demonstrated that several mono- and sesquiterpenes were de novo synthesized upon MeJA and B. vietnamiensis applications. This study demonstrates that foliar application of endophytic bacteria B. vietnamiensis can potentially enhance resistance to biotic stresses and contribute to the maintenance of the integrity of plant metabolic activity.

The Road to Resistance in Forest Trees

In recent years, forests have been exposed to an unprecedented rise in pests and pathogens. This, coupled with the added challenge of climate change, renders forest plantation stock vulnerable to attack and severely limits productivity. Genotypes resistant to such biotic challenges are desired in plantation forestry to reduce losses. Conventional breeding has been a main avenue to obtain resistant genotypes. More recently, genetic engineering has become a viable approach to develop resistance against pests and pathogens in forest trees. Tree genomic resources have contributed to advancements in both these approaches. Genome-wide association studies and genomic selection in tree populations have accelerated breeding tools while integration of various levels of omics information facilitates the selection of candidate genes for genetic engineering. Furthermore, tree associations with non-pathogenic endophytic and subterranean microbes play a critical role in plant health and may be engineered in forest trees to improve resistance in the future. We look at recent studies in forest trees describing defense mechanisms using such approaches and propose the way forward to developing superior genotypes with enhanced resistance against biotic stress.

Transcriptome Analysis of Chinese Chestnut (Castanea mollissima Blume) in Response to Dryocosmus kuriphilus Yasumatsu Infestation

Chinese chestnut (Castanea mollissima Blume) can be infested by Dryocosmus kuriphilus Yasumatsu, resulting in gall formation and yield losses. Research on the control of gall wasps using genomics approaches is rarely reported. We used RNA-seq to investigate the dynamic changes in the genes of a chestnut species (C. mollissima B.) during four gall-formation stages caused by D. kuriphilus. A total of 21,306 genes were annotated by BLAST in databases. Transcriptome comparison between different gall-formation stages revealed many genes that were differentially expressed compared to the control. Among these, 2410, 7373, 6294, and 9412 genes were differentially expressed in four gall-formation stages: initiation stage (A), early growth stage (B), late growth stage (C), and maturation stage (D), respectively. Annotation analysis indicated that many metabolic processes (e.g., phenylpropanoid biosynthesis, secondary metabolism, plant⁻pathogen interaction) were affected. Interesting genes encoding putative components of signal transduction, stress response, and transcription factors were also differentially regulated. These genes might play important roles in response to D. kuriphilus gall formation. These new data on the mechanism by which D. kuriphilus infests chestnuts could help improve chestnut resistance.

Hardwood Tree Genomics: Unlocking Woody Plant Biology

Woody perennial angiosperms (i.e., hardwood trees) are polyphyletic in origin and occur in most angiosperm orders. Despite their independent origins, hardwoods have shared physiological, anatomical, and life history traits distinct from their herbaceous relatives. New high-throughput DNA sequencing platforms have provided access to numerous woody plant genomes beyond the early reference genomes of Populus and Eucalyptus, references that now include willow and oak, with pecan and chestnut soon to follow. Genomic studies within these diverse and undomesticated species have successfully linked genes to ecological, physiological, and developmental traits directly. Moreover, comparative genomic approaches are providing insights into speciation events while large-scale DNA resequencing of native collections is identifying population-level genetic diversity responsible for variation in key woody plant biology across and within species. Current research is focused on developing genomic prediction models for breeding, defining speciation and local adaptation, detecting and characterizing somatic mutations, revealing the mechanisms of gender determination and flowering, and application of systems biology approaches to model complex regulatory networks underlying quantitative traits. Emerging technologies such as single-molecule, long-read sequencing is being employed as additional woody plant species, and genotypes within species, are sequenced, thus enabling a comparative ("evo-devo") approach to understanding the unique biology of large woody plants. Resource availability, current genomic and genetic applications, new discoveries and predicted future developments are illustrated and discussed for poplar, eucalyptus, willow, oak, chestnut, and pecan.

Transcriptome profile of cup-shaped galls in Litsea acuminata leaves

BACKGROUND

Insect galls are atypical plant tissues induced by the invasion of insects. Compared to the host leaf, gall tissues lose photosynthetic ability, but have higher soluble sugar content. Although the physiological and biochemical regulation of gall tissues have been demonstrated, the mechanism of genetic regulation has only been analyzed in few studies.

RESULTS

In this study, the transcriptome of cup-shaped galls and its host leaf were de novo assembled. Cellular functional enrichment and differentially expressed gene groups in the gall tissues were analyzed. The genes associated with primary metabolism, including photosynthesis, cell wall turnover, and sugar degradation, were expressed differently in galls and leaves. The examination of gene expression demonstrated that the genes involved in brassinosteroid synthesis and responses exhibited a remarkable modulation in cup-shaped galls, suggesting a potential role of steroid hormones in regulating gall development.

CONCLUSIONS

This study revealed the genetic responses, including those involved in source-sink reallocation and phytohormone metabolism, of galls induced by a dipteran insect.

Terpenes associated with resistance against the gall wasp, Leptocybe invasa, in Eucalyptus grandis

Leptocybe invasa is an insect pest causing gall formation on oviposited shoot tips and leaves of Eucalyptus trees leading to leaf deformation, stunting, and death in severe cases. We previously observed different constitutive and induced terpenes, plant specialized metabolites that may act as attractants or repellents to insects, in a resistant and susceptible clone of Eucalyptus challenged with L. invasa. We tested the hypothesis that specific terpenes are associated with pest resistance in a Eucalyptus grandis half-sib population. Insect damage was scored over 2 infestation cycles, and leaves were harvested for near-infrared reflectance (NIR) and terpene measurements. We used Bayesian model averaging for terpene selection and obtained partial least squares NIR models to predict terpene content and L. invasa infestation damage. In our optimal model, 29% of the phenotypic variation could be explained by 7 terpenes, and the monoterpene combination, limonene, α-terpineol, and 1,8-cineole, could be predicted with an NIR prediction ability of  .67. Bayesian model averaging supported α-pinene, γ-terpinene, and iso-pinocarveol as important for predicting L. invasa infestation. Susceptibility was associated with increased γ-terpinene and α-pinene, which may act as a pest attractant, whereas reduced susceptibility was associated with iso-pinocarveol, which may act to recruit parasitoids or have direct toxic effects.

De Novo Transcriptome Study Identifies Candidate Genes Involved in Resistance to Austropuccinia psidii (Myrtle Rust) in Syzygium luehmannii (Riberry)

Austropuccinia psidii, causal agent of myrtle rust, was discovered in Australia in 2010 and has since become established on a wide range of species within the family Myrtaceae. Syzygium luehmannii, endemic to Australia, is an increasingly valuable berry crop. Plants were screened for responses to A. psidii inoculation, and specific resistance, in the form of localized necrosis, was determined in 29% of individuals. To understand the molecular basis underlying this response, mRNA was sequenced from leaf samples taken preinoculation, and at 24 and 48 h postinoculation, from four resistant and four susceptible plants. Analyses, based on de novo transcriptome assemblies for all plants, identified significant expression changes in resistant plants (438 transcripts) 48 h after pathogen exposure compared with susceptible plants (three transcripts). Most significantly up-regulated in resistant plants were gene homologs for transcription factors, receptor-like kinases, and enzymes involved in secondary metabolite pathways. A putative G-type lectin receptor-like kinase was exclusively expressed in resistant individuals and two transcripts incorporating toll/interleukin-1, nucleotide binding site, and leucine-rich repeat domains were up-regulated in resistant plants. The results of this study provide the first early gene expression profiles for a plant of the family Myrtaceae in response to the myrtle rust pathogen.

Transcriptome Profiling of Melaleuca quinquenervia Challenged by Myrtle Rust Reveals Differences in Defense Responses Among Resistant Individuals

Plants have developed complex defense mechanisms to protect themselves against pathogens. A wide-host-range fungus, Austropuccinia psidii, which has caused severe damage to ecosystems and plantations worldwide, is a major threat to Australian ecosystems dominated by members of the family Myrtaceae. In particular, the east coast wetland foundation tree species Melaleuca quinquenervia, appears to be variably susceptible to this pathogen. Understanding the molecular basis of host resistance would enable better management of this rust disease. We identified resistant and susceptible individuals of M. quinquenervia and explored their differential gene expression in order to discover the molecular basis of resistance against A. psidii. Rust screening of germplasm showed a varying degree of response, with fully resistant to highly susceptible individuals. We used transcriptome profiling in samples collected before and at 5 days postinoculation (dpi). Differential gene expression analysis showed that numerous defense-related genes were induced in susceptible plants at 5 dpi. Mapping reads against the A. psidii genome showed that only susceptible plants contained fungal-derived transcripts. Resistant plants exhibited an overexpression of candidate A. psidii resistance-related genes such as receptor-like kinases, nucleotide-binding site leucine-rich repeat proteins, glutathione S-transferases, WRKY transcriptional regulators, and pathogenesis-related proteins. We identified large differences in the expression of defense-related genes among resistant individuals.

The first identification of genomic loci in plants associated with resistance to galling insects: a case study in Eucalyptus L'Hér. (Myrtaceae)

Genomic loci related with resistance to gall-inducing insects have not been identified in any plants. Here, association mapping was used to identify molecular markers for resistance to the gall wasp Leptocybe invasa in two Eucalyptus species. A total of 86 simple sequence repeats (SSR) markers were screened out from 839 SSRs and used for association mapping in E. grandis. By applying the mixed linear model, seven markers were identified to be associated significantly (P ≤ 0.05) with the gall wasp resistance in E. grandis, including two validated with a correction of permutation test (P ≤ 0.008). The proportion of the variance in resistance explained by a significant marker ranged from 3.3% to 37.8%. Four out of the seven significant associations in E. grandis were verified and also validated (P ≤ 0.073 in a permutation test) in E. tereticornis, with the variation explained ranging from 24.3% to 48.5%. Favourable alleles with positive effect were also mined from the significant markers in both species. These results provide insight into the genetic control of gall wasp resistance in plants and have great potential for marker-assisted selection for resistance to L. invasa in the important tree genus Eucalyptus.

The first identification of genomic loci in plants associated with resistance to galling insects: a case study in Eucalyptus L'Hér. (Myrtaceae)

Genomic loci related with resistance to gall-inducing insects have not been identified in any plants. Here, association mapping was used to identify molecular markers for resistance to the gall wasp Leptocybe invasa in two Eucalyptus species. A total of 86 simple sequence repeats (SSR) markers were screened out from 839 SSRs and used for association mapping in E. grandis. By applying the mixed linear model, seven markers were identified to be associated significantly (P ≤ 0.05) with the gall wasp resistance in E. grandis, including two validated with a correction of permutation test (P ≤ 0.008). The proportion of the variance in resistance explained by a significant marker ranged from 3.3% to 37.8%. Four out of the seven significant associations in E. grandis were verified and also validated (P ≤ 0.073 in a permutation test) in E. tereticornis, with the variation explained ranging from 24.3% to 48.5%. Favourable alleles with positive effect were also mined from the significant markers in both species. These results provide insight into the genetic control of gall wasp resistance in plants and have great potential for marker-assisted selection for resistance to L. invasa in the important tree genus Eucalyptus.

The gall wasp Leptocybe invasa (Hymenoptera: Eulophidae) stimulates different chemical and phytohormone responses in two Eucalyptus varieties that vary in susceptibility to galling

Gall-inducing insects produce various types of galls on plants, but little is known about the gall-induction mechanism of these galling insects. The gall wasp Leptocybe invasa Fisher & LaSalle (Hymenoptera: Eulophidae) forms galls of different sizes on several Eucalyptus species. To clarify the physiological responses of Eucalyptus to L. invasa infestation, we measured the dynamics of nitrogen (N), carbon (C), total phenolics, total tannins and four types of phytohormones (zeatin [Z] + zeatin riboside [ZR], gibberellins [GA], indole-3-acetic acid [IAA] and abscisic acid [ABA]) in galled and ungalled leaf tissues of two Eucalyptus horticultural varieties (DH201-2 [Eucalyptus grandis × Eucalyptus camaldulensis] and EA [Eucalyptus exserta]) with different susceptibility to galling throughout the larval developmental stages. Nitrogen, total phenolics, tannins and four kinds of phytohormones strongly accumulated in tissues galled by L. invasa (especially during early larval feeding stages). While N, Z + ZR and GA levels were higher, tannins and ABA levels were lower in the galled tissues on the highly susceptible variety. Nitrogen, total phenolics, GA, Z + ZR and IAA levels in the galled tissues gradually decreased during gall development, but ABA and tannins conversely increased in the galled tissues of the less susceptible variety. Our results suggest that the effects of gall-inducing insects on plants depend not only on the susceptibility of the plant infested but also on the developmental stage of galled tissues. Gall formation process is thus synergistically influenced by both gall-inducing insect and plant genotypes.

Insect Gallers and Their Plant Hosts: From Omics Data to Systems Biology

Gall-inducing insects are capable of exerting a high level of control over their hosts' cellular machinery to the extent that the plant's development, metabolism, chemistry, and physiology are all altered in favour of the insect. Many gallers are devastating pests in global agriculture and the limited understanding of their relationship with their hosts prevents the development of robust management strategies. Omics technologies are proving to be important tools in elucidating the mechanisms involved in the interaction as they facilitate analysis of plant hosts and insect effectors for which little or no prior knowledge exists. In this review, we examine the mechanisms behind insect gall development using evidence from omics-level approaches. The secretion of effector proteins and induced phytohormonal imbalances are highlighted as likely mechanisms involved in gall development. However, understanding how these components function within the system is far from complete and a number of questions need to be answered before this information can be used in the development of strategies to engineer or breed plants with enhanced resistance.

The Eucalyptus grandis NBS-LRR Gene Family: Physical Clustering and Expression Hotspots

Eucalyptus grandis is a commercially important hardwood species and is known to be susceptible to a number of pests and pathogens. Determining mechanisms of defense is therefore a research priority. The published genome for E. grandis has aided the identification of one important class of resistance (R) genes that incorporate nucleotide binding sites and leucine-rich repeat domains (NBS-LRR). Using an iterative search process we identified NBS-LRR gene models within the E. grandis genome. We characterized the gene models and identified their genomic arrangement. The gene expression patterns were examined in E. grandis clones, challenged with a fungal pathogen (Chrysoporthe austroafricana) and insect pest (Leptocybe invasa). One thousand two hundred and fifteen putative NBS-LRR coding sequences were located which aligned into two large classes, Toll or interleukin-1 receptor (TIR) and coiled-coil (CC) based on NB-ARC domains. NBS-LRR gene-rich regions were identified with 76% organized in clusters of three or more genes. A further 272 putative incomplete resistance genes were also identified. We determined that E. grandis has a higher ratio of TIR to CC classed genes compared to other woody plant species as well as a smaller percentage of single NBS-LRR genes. Transcriptome profiles indicated expression hotspots, within physical clusters, including expression of many incomplete genes. The clustering of putative NBS-LRR genes correlates with differential expression responses in resistant and susceptible plants indicating functional relevance for the physical arrangement of this gene family. This analysis of the repertoire and expression of E. grandis putative NBS-LRR genes provides an important resource for the identification of novel and functional R-genes; a key objective for strategies to enhance resilience.