Response of living tissues of Pinus sylvestris to the saprotrophic biocontrol fungus Phlebiopsis gigantea

An Experiment in Denmark to Infect Wounded Sitka Spruce with the Rotstop Isolate of Phlebiopsis gigantea, and Its Implications for the Control of Heterobasidion annosum in Britain

The formulation of a Finnish isolate of the saprotrophic wood-rotting fungus Phlebiopsis gigantea into the biocontrol agent (BCA) Rotstop, which is used to prevent infection of Norway spruce stumps by aerial basidiospores of H. annosum, has led to its application to more than 200,000 ha of forest in Scandinavia and North Europe. The success of this treatment opens the possibility of introducing the Rotstop strain into Britain for use on Sitka spruce stumps, which at present (2022) are lacking adequate prophylactic treatment. However, Rotstop is probably non-native to Britain and to North America (the ancestral home of this spruce), and we do not know if this xylem-decaying BCA can invade standing trees. In this paper, we describe a trial into this issue conducted for the U.K. Forestry Commission in Denmark, in a country where both Rotstop and Sitka spruce have been naturalised. It was preliminary to further stump treatment trials, and is relevant to long-term issues surrounding stump treatment in Britain. Inoculations into 44-year-old standing Sitka spruce with 20 mm wooden Scots pine plugs pre-colonised with Rotstop resulted in decay of the S1, S2, S3 and middle lamellae of sapwood above and below the wounds after 11–18 months. In contrast, infection of sapwood occurred in only one of 39 wounds on the same trees challenged with oidial spore inoculants adpressed to undamaged xylem sapwood during the same period. While the results suggest that release of Rotstop into the productive stands of Sitka spruce in Britain would be unlikely to cause long-term commercial losses to wounded trees, the work highlights issues relating to the assessment of risk associated with the introduction of non-native BCAs within the forest environment.

Fourier-transform infrared (FT-IR) spectroscopy analysis discriminates asymptomatic and symptomatic Norway spruce trees

Conifer trees, including Norway spruce, are threatened by fungi of the Heterobasidion annosum species complex, which severely affect timber quality and cause economic losses to forest owners. The timely detection of infected trees is complicated, as the pathogen resides within the heartwood and sapwood of infected trees. The presence of the disease and the extent of the wood decay often becomes evident only after tree felling. Fourier-transform infrared (FT-IR) spectroscopy is a potential method for non-destructive sample analysis that may be useful for identifying infected trees in this pathosystem. We performed FT-IR analysis of 18 phloem, 18 xylem, and 18 needle samples from asymptomatic and symptomatic Norway spruce trees. FT-IR spectra from 1066 - 912 cm^-1 could be used to distinguish phloem, xylem, and needle tissue extracts. FT-IR spectra collected from xylem and needle extracts could also be used to discriminate between asymptomatic and symptomatic trees using spectral bands from 1657 - 994 cm^-1 and 1104 - 994 cm^-1, respectively. A partial least squares regression model predicted the concentration of condensed tannins, a defense-related compound, in phloem of asymptomatic and symptomatic trees. This work is the first to show that FT-IR spectroscopy can be used for the identification of Norway spruce trees naturally infected with Heterobasidion spp.

Comprehensive assembly and analysis of the transcriptome of maritime pine developing embryos

BACKGROUND

There are clear differences in embryo development between angiosperm and gymnosperm species. Most of the current knowledge on gene expression and regulation during plant embryo development has derived from studies on angiosperms species, in particular from the model plant Arabidopsis thaliana. The few published studies on transcript profiling of conifer embryogenesis show the existence of many putative embryo-specific transcripts without an assigned function. In order to extend the knowledge on the transcriptomic expression during conifer embryogenesis, we sequenced the transcriptome of zygotic embryos for several developmental stages that cover most of Pinus pinaster (maritime pine) embryogenesis.

RESULTS

Total RNA samples collected from five zygotic embryo developmental stages were sequenced with Illumina technology. A de novo transcriptome was assembled as no genome sequence is yet published for Pinus pinaster. The transcriptome of reference for the period of zygotic embryogenesis in maritime pine contains 67,429 transcripts, which likely encode 58,527 proteins. The annotation shows a significant percentage, 31%, of predicted proteins exclusively present in pine embryogenesis. Functional categories and enrichment analysis of the differentially expressed transcripts evidenced carbohydrate transport and metabolism over-representation in early embryo stages, as highlighted by the identification of many putative glycoside hydrolases, possibly associated with cell wall modification, and carbohydrate transport transcripts. Moreover, the predominance of chromatin remodelling events was detected in early to middle embryogenesis, associated with an active synthesis of histones and their post-translational modifiers related to increased transcription, as well as silencing of transposons.

CONCLUSIONS

Our results extend the understanding of gene expression and regulation during zygotic embryogenesis in conifers and are a valuable resource to support further improvements in somatic embryogenesis for vegetative propagation of conifer species. Specific transcripts associated with carbohydrate metabolism, monosaccharide transport and epigenetic regulation seem to play an important role in pine early embryogenesis and may be a source of reliable molecular markers for early embryogenesis.

Genome-Wide Exon-Capture Approach Identifies Genetic Variants of Norway Spruce Genes Associated With Susceptibility to Heterobasidion parviporum Infection

Root and butt rot caused by members of the Heterobasidion annosum species complex is the most economically important disease of conifer trees in boreal forests. Wood decay in the infected trees dramatically decreases their value and causes considerable losses to forest owners. Trees vary in their susceptibility to Heterobasidion infection, but the genetic determinants underlying the variation in the susceptibility are not well-understood. We performed the identification of Norway spruce genes associated with the resistance to Heterobasidion parviporum infection using genome-wide exon-capture approach. Sixty-four clonal Norway spruce lines were phenotyped, and their responses to H. parviporum inoculation were determined by lesion length measurements. Afterwards, the spruce lines were genotyped by targeted resequencing and identification of genetic variants (SNPs). Genome-wide association analysis identified 10 SNPs located within 8 genes as significantly associated with the larger necrotic lesions in response to H. parviporum inoculation. The genetic variants identified in our analysis are potential marker candidates for future screening programs aiming at the differentiation of disease-susceptible and resistant trees.

Secondary Metabolites from the Root Rot Biocontrol Fungus Phlebiopsis gigantea

Three cyclopentanoids (phlebiopsin A⁻C), one glycosylated p-terphenyl (methyl-terfestatin A), and o-orsellinaldehyde were isolated from the biocontrol fungus Phlebiopsis gigantea, and their structures were elucidated by 1D and 2D NMR spectroscopic analysis, as well as by LC-HRMS. The biological activity of the compounds against the root rot fungus Heterobasidion occidentale, as well as against Fusarium oxysporum and Penicillium canescens, was also investigated, but only o-orsellinaldehyde was found to have any antifungal activity in the concentration range tested.

Modular organization of the white spruce (Picea glauca) transcriptome reveals functional organization and evolutionary signatures

Transcript profiling has shown the molecular bases of several biological processes in plants but few studies have developed an understanding of overall transcriptome variation. We investigated transcriptome structure in white spruce (Picea glauca), aiming to delineate its modular organization and associated functional and evolutionary attributes. Microarray analyses were used to: identify and functionally characterize groups of co-expressed genes; investigate expressional and functional diversity of vascular tissue preferential genes which were conserved among Picea species, and identify expression networks underlying wood formation. We classified 22 857 genes as variable (79%; 22 coexpression groups) or invariant (21%) by profiling across several vegetative tissues. Modular organization and complex transcriptome restructuring among vascular tissue preferential genes was revealed by their assignment to coexpression groups with partially overlapping profiles and partially distinct functions. Integrated analyses of tissue-based and temporally variable profiles identified secondary xylem gene networks, showed their remodelling over a growing season and identified PgNAC-7 (no apical meristerm (NAM), Arabidopsis transcription activation factor (ATAF) and cup-shaped cotyledon (CUC) transcription factor 007 in Picea glauca) as a major hub gene specific to earlywood formation. Reference profiling identified comprehensive, statistically robust coexpressed groups, revealing that modular organization underpins the evolutionary conservation of the transcriptome structure.

Strategies utilized by trophically diverse fungal species for Pinus sylvestris root colonization

Physiological changes in host plants in response to the broad spectrum of fungal modes of infection are still not well understood. The current study was conducted to better understand the infection of in vitro cultures of Pinus sylvestris L. seedlings by three trophically diverse fungal species, Fusarium oxysporum E. F. Sm. & Swingle, Trichoderma harzianum Rifai and Hebeloma crustuliniforme (Bull.) Quél. Biochemical methods and microscopy were utilized to determine (i) which factors (apoplastic and cellular pH, reactive oxygen species, glutathione and cell death) play a role in the establishment of pathogenic, saprotrophic and mycorrhizal fungi, and (ii) whether cell death is a common response of conifer seedling tissues when they are exposed to trophically diverse fungi. Establishment of the pathogen, F. oxysporum, was observed more frequently in the meristematic region of root tips than in the elongation zone, which was in contrast to T. harzianum and H. crustuliniforme. Ectomycorrhizal (ECM) hyphae, however, were occasionally observed in the studied root zone and caused small changes in the studied factors. Colonization of the meristematic zone occurred due to host cell death. Independently of the zone, changes in cellular pH resulting in an acidic cytoplasm conditioned the establishment of F. oxysporum. Additionally, cell death was negatively correlated with hydrogen peroxide (H2O2) in roots challenged by a pathogenic fungus. Cell death was the only factor uniquely associated with the colonization of host roots by a saprotrophic fungus. The mechanism may differ, however, between the zones since apoplastic pH was negatively correlated with cell death in the elongation zone, whereas in the meristematic zone, none of the studied factors explained cell death. Colonization by the ECM fungus, H. crustuliniforme, was associated with a decreasing number of cells with acidic apoplast and by production of H2O2 in the elongation zone resulting in cell death. Saprotrophic and ECM fungi had a greater effect on cell acidification in the meristematic zone than the pathogenic fungus.

Expression analysis of the impact of culture filtrates from the biocontrol agent, Phlebiopsis gigantea on the conifer pathogen, Heterobasidion annosum s.s. Transcriptome

Phlebiopsis gigantea has been routinely used as the biological control agent for the conifer pathogen Heterobasidion annosum sensu lato, but the actual mechanism for the biocontrol process is not known. To investigate the effect of secreted molecules from culture filtrate produced by P. gigantea on the gene expression profile of H. annosum s.s., microarray analysis was used. Analysis of the differentially expressed genes led to the identification of genes with diverse functions. A major proportion of the up- and downregulated genes were either uncharacterized or genes whose functions were not known. A number of genes coding for proteins involved in metabolism, transport, and signal transduction were differentially downregulated; comparatively lower number of such genes were upregulated. Some genes involved in transport (polyamine transporters, 2573-fold, P = 0.002) and metabolism (endoglucanase, 622.5-fold, P = 0.002, cytochrome P450, 133.2-fold, P = 0.05) showed high transcript fold changes and were statistically significantly upregulated. Genes encoding defense-related proteins such as hydrophobins were either downregulated or expressed at relatively low levels. Further analysis of the effect of the culture filtrate on glucose metabolism showed downregulation of some key enzymes at the early stage of the glycolytic pathway while some genes were upregulated at the later stage of the pathway. A subset of the genes were selected and used to validate the micro-array result by quantitative real time polymerase chain reaction (qPCR) method. Generally, the high transcript levels of genes encoding several biochemically important genes (protein kinases, major facilitator superfamily polyamine transporters, endoglucanase, cytochrome P450, endoglucanase) suggests their potential functional relevance in signal perception, stress tolerance, cell defenses, and detoxification of toxic molecules during competitive interaction. These results have provided further insights into possible molecular and genetic factors underlying the response of H. annosum to metabolites from P. gigantea during interspecific interaction.

The primary module in Norway spruce defence signalling against H. annosum s.l. seems to be jasmonate-mediated signalling without antagonism of salicylate-mediated signalling

A key tree species for the forest industry in Europe is Norway spruce [Picea abies (L.) Karst.]. One of its major diseases is stem and butt rot caused by Heterobasidion parviporum (Fr.) Niemelä & Korhonen, which causes extensive revenue losses every year. In this study, we investigated the parallel induction of Norway spruce genes presumably associated with salicylic acid- and jasmonic acid/ethylene-mediated signalling pathways previously observed in response to H. parviporum. Relative gene expression levels in bark samples of genes involved in the salicylic acid- and jasmonic acid/ethylene-mediated signalling pathways after wounding and inoculation with either the saprotrophic biocontrol fungus Phlebiopsis gigantea or with H. parviporum were analysed with quantitative PCR at the site of the wound and at two distal locations from the wound/inoculation site to evaluate their roles in the induced defence response to H. parviporum in Norway spruce. Treatment of Norway spruce seedlings with methylsalicylate, methyljasmonate and inhibitors of the jasmonic acid/ethylene signalling pathway, as well as the Phenylalanine ammonia lyase inhibitor 2-aminoindan-2-phosphonic acid were conducted to determine the responsiveness of genes characteristic of the different pathways to different hormonal stimuli. The data suggest that jasmonic acid-mediated signalling plays a central role in the induction of the genes analysed in this study irrespective of their responsiveness to salicylic acid. This may suggest that jasmonic acid-mediated signalling is the prioritized module in the Norway spruce defence signalling network against H. parviporum and that there seems to be no immediate antagonism between the modules in this interaction.

RNA-Seq Assembly - Are We There Yet?

Transcriptomic sequence resources represent invaluable assets for research, in particular for non-model species without a sequenced genome. To date, the Next Generation Sequencing technologies 454/Roche and Illumina have been used to generate transcriptome sequence databases by mRNA-Seq for more than fifty different plant species. While some of the databases were successfully used for downstream applications, such as proteomics, the assembly parameters indicate that the assemblies do not yet accurately reflect the actual plant transcriptomes. Two different assembly strategies have been used, overlap consensus based assemblers for long reads and Eulerian path/de Bruijn graph assembler for short reads. In this review, we discuss the challenges and solutions to the transcriptome assembly problem. A list of quality control parameters and the necessary scripts to produce them are provided.