Mårten Lind

Insight into trade???off between wood decay and parasitism from the genome of a fungal forest pathogen

Parasitism and saprotrophic wood decay are two fungal strategies fundamental for succession and nutrient cycling in forest ecosystems. An opportunity to assess the trade-off between these strategies is provided by the forest pathogen and wood decayer Heterobasidion annosum sensu lato. We report the annotated genome sequence and transcript profiling, as well as the quantitative trait loci mapping, of one member of the species complex: H. irregulare. Quantitative trait loci critical for pathogenicity, and rich in transposable elements, orphan and secreted genes, were identified. A wide range of cellulose-degrading enzymes are expressed during wood decay. By contrast, pathogenic interaction between H. irregulare and pine engages fewer carbohydrate-active enzymes, but involves an increase in pectinolytic enzymes, transcription modules for oxidative stress and secondary metabolite production. Our results show a trade-off in terms of constrained carbohydrate decomposition and membrane transport capacity during interaction with living hosts. Our findings establish that saprotrophic wood decay and necrotrophic parasitism involve two distinct, yet overlapping, processes.

A Genome-Wide Association Study Identifies Genomic Regions for Virulence in the Non-Model Organism Heterobasidion annosum s.s

The dense single nucleotide polymorphisms (SNP) panels needed for genome wide association (GWA) studies have hitherto been expensive to establish and use on non-model organisms. To overcome this, we used a next generation sequencing approach to both establish SNPs and to determine genotypes. We conducted a GWA study on a fungal species, analysing the virulence of Heterobasidion annosum s.s., a necrotrophic pathogen, on its hosts Picea abies and Pinus sylvestris. From a set of 33,018 single nucleotide polymorphisms (SNP) in 23 haploid isolates, twelve SNP markers distributed on seven contigs were associated with virulence (P<0.0001). Four of the contigs harbour known virulence genes from other fungal pathogens and the remaining three harbour novel candidate genes. Two contigs link closely to virulence regions recognized previously by QTL mapping in the congeneric hybrid H. irregulare × H. occidentale. Our study demonstrates the efficiency of GWA studies for dissecting important complex traits of small populations of non-model haploid organisms with small genomes.

Identification of quantitative trait loci affecting virulence in the basidiomycete Heterobasidion annosum s.l.

Identification of virulence factors of phytopathogens is important for the fundamental understanding of infection and disease progress in plants and for the development of control strategies. We have identified quantitative trait loci (QTL) for virulence on 1-year-old Pinus sylvestris and 2-year-old Picea abies seedlings and positioned them on a genetic linkage map of the necrotrophic phytopathogen Heterobasidion annosum sensu lato (s.l.), a major root rot pathogen on conifers. The virulence of 102 progeny isolates was analysed using two measurements: lesion lengths and fungal growth in sapwood from a cambial infection site. We found negative virulence effects of hybridization although this was contradicted on a winter-hardened spruce. On P. abies, both measurements identified several partially overlapping QTLs on linkage group (LG) 15 of significant logarithm of odds (LOD) values ranging from 2.31 to 3.85. On P. sylvestris, the lesion length measurement also identified a QTL (LOD 3.09) on LG 15. Moreover, QTLs on two separate smaller LGs, with peak LOD values of 2.78 and 4.58 were identified for fungal sapwood growth and lesion lengths, respectively. The QTL probably represent loci important for specific as well as general aspects of virulence on P. sylvestris and P. abies.

A Picea abies linkage map based on SNP markers identifies QTLs for four aspects of resistance to Heterobasidion parviporum infection.

A consensus linkage map of Picea abies, an economically important conifer, was constructed based on the segregation of 686 SNP markers in a F1 progeny population consisting of 247 individuals. The total length of 1889.2 cM covered 96.5% of the estimated genome length and comprised 12 large linkage groups, corresponding to the number of haploid P. abies chromosomes. The sizes of the groups (from 5.9 to 9.9% of the total map length) correlated well with previous estimates of chromosome sizes (from 5.8 to 10.8% of total genome size). Any locus in the genome has a 97% probability to be within 10 cM from a mapped marker, which makes the map suited for QTL mapping. Infecting the progeny trees with the root rot pathogen Heterobasidion parviporum allowed for mapping of four different resistance traits: lesion length at the inoculation site, fungal spread within the sapwood, exclusion of the pathogen from the host after initial infection, and ability to prevent the infection from establishing at all. These four traits were associated with two, four, four and three QTL regions respectively of which none overlapped between the traits. Each QTL explained between 4.6 and 10.1% of the respective traits phenotypic variation. Although the QTL regions contain many more genes than the ones represented by the SNP markers, at least four markers within the confidence intervals originated from genes with known function in conifer defence; a leucoanthocyanidine reductase, which has previously been shown to upregulate during H. parviporum infection, and three intermediates of the lignification process; a hydroxycinnamoyl CoA shikimate/quinate hydroxycinnamoyltransferase, a 4-coumarate CoA ligase, and a R2R3-MYB transcription factor.

The pathogenic white-rot fungus Heterobasidion parviporum triggers non-specific defence responses in the bark of Norway spruce

Norway spruce [Picea abies (L.) Karst.] is one of the economically most important conifer species in Europe. The major pathogen on Norway spruce is Heterobasidion parviporum (Fr.) Niemelä & Korhonen. To achieve a better understanding of Norway spruces defence mechanisms, transcriptional responses in bark to H. parviporum infection were compared with the response to wounding using cDNA-amplified fragment length polymorphism. The majority of the recovered transcript-derived fragments (TDFs) showed a similar expression pattern for infection and wounding treatment, although inoculated samples showed an enhanced reaction. Genes related to systemic acquired resistance, e.g., PR1, accumulated after H. parviporum infection. Simultaneously, several transcripts involved in various aspects of jasmonic acid (JA)- and ethylene (ET)-mediated signalling accumulated. Genes involved in the ubiquitin/proteasome system were also regulated. Expression patterns have been confirmed by quantitative polymerase chain reaction. The expression patterns of the isolated TDFs suggest that infection with H. parviporum in Norway spruce induces a broad defence, with many similarities to non-specific defence responses in angiosperms. The parallel induction of salicylic acid- and JA/ET-mediated pathways implies spatially separated responses in different cell layers, with and without hyphal contact. A set of TDFs were analysed in an independent experiment with unrelated material treated with wounding or with inoculation with H. parviporum or Phlebiopsis gigantea, verifying the original observations and underlining the non-specific defence responses. In addition, our data suggest that rerouting of carbon in secondary metabolism is an integral part of Norway spruce induced defence. We report the sequences of three 3-deoxy-d-arabino-heptulosonate 7-phosphate synthase genes (PaDAHP1, PaDAHP2 and PaDAHP3) and their relative expression in response to wounding and infection with H. parviporum and P. gigantea. The results clearly indicate differential regulation of the three DAHPs in the induced defence responses in Norway spruce. This study gives insights into the central mechanisms in the induced defences in Norway spruce.

Distribution and evolution of het gene homologs in the basidiomycota

In filamentous fungi a system known as somatic incompatibility (SI) governs self/non-self recognition. SI is controlled by a regulatory signaling network involving proteins encoded at the het (heterokaryon incompatible) loci. Despite the wide occurrence of SI, the molecular identity and structure of only a small number of het genes and their products have been characterized in the model fungi Neurospora crassa and Podospora anserina. Our aim was to identify and study the distribution and evolution of putative het gene homologs in the Basidiomycota. For this purpose we used the information available for the model fungi to identify homologs of het genes in other fungi, especially the Basidiomycota. Putative het-c, het-c2 and un-24 homologs, as well as sequences containing the NACHT, HET or WD40 domains present in the het-e, het-r, het-6 and het-d genes were identified in certain members of the Ascomycota and Basidiomycota. The widespread phylogenetic distribution of certain het genes may reflect the fact that the encoded proteins are involved in fundamental cellular processes other than SI. Although homologs of het-S were previously known only from the Sordariomycetes (Ascomycota), we also identified a putative homolog of this gene in Gymnopus luxurians (Basidiomycota, class Agaricomycetes). Furthermore, with the exception of un-24, all of the putative het genes identified occurred mostly in a multi-copy fashion, some with lineage and species-specific expansions. Overall our results indicated that gene duplication followed by gene loss and/or gene family expansion, as well as multiple events of domain fusion and shuffling played an important role in the evolution of het gene homologs of Basidiomycota and other filamentous fungi.

A 2nd Generation Linkage Map of Heterobasidion annosum s.l. Based on In Silico Anchoring of AFLP Markers

In this study, we present a 2nd generation genetic linkage map of a cross between the North American species Heterobasidion irregulare and H. occidentale, based on the alignment of the previously published 1st generation map to the parental genomes. We anchored 216 of the original 308 AFLP markers to their respective restriction sites using an in silico-approach. The map resolution was improved by adding 146 sequence-tagged microsatellite markers and 39 sequenced gene markers. The new markers confirmed the positions of the anchored AFLP markers, fused the original 39 linkage groups together into 17, and fully expanded 12 of these to single groups covering entire chromosomes. Map coverage of the genome increased from 55.3% to 92.8%, with 96.3% of 430 markers collinearly aligned with the genome sequence. The anchored map also improved the H. irregulare assembly considerably. It identified several errors in scaffold arrangements and assisted in reducing the total number of major scaffolds from 18 to 15. This denser, more comprehensive map allowed sequence-based mapping of three intersterility loci and one mating type locus. This demonstrates the possibility to utilize an in silico procedure to convert anonymous markers into sequence-tagged ones, as well as the power of a sequence-anchored linkage map and its usefulness in the assembly of a whole genome sequence.

Gene expression associated with intersterility in Heterobasidion

Intersterility (IS) is thought to prevent mating compatibility between homokaryons that belong to different species. Although IS in Heterobasidion is regulated by the genes located at the IS loci, it is not yet known how the IS genes influence sexual compatibility and heterokaryon formation. To increase our understanding of the molecular events underlying IS, we studied mRNA abundance changes during IS compatible and incompatible interactions over time. The clustering of the transcripts into expression profiles, followed by the application of Gene Ontology (GO) enrichment pathway analysis of each of the clusters, allowed inference of biological processes participating in IS. These analyses identified events involved in mating and sexual development (i.e., linked with IS compatibility), which included processes associated with cell-cell adhesion and recognition, cell cycle control and signal transduction. We also identified events potentially involved in overriding mating between individuals belonging to different species (i.e., linked with IS incompatibility), which included reactive oxygen species (ROS) production, responses to stress (especially to oxidative stress), signal transduction and metabolic biosynthesis. Our findings thus enabled detection and characterization of gene expression changes associated with IS in Heterobasidion, as well as identification of important processes and pathways associated with this phenomenon. Overall, the results of this study increase current knowledge regarding the molecular mechanisms underpinning IS in Heterobasidion and allowed for the establishment of a vital baseline for further studies.

Heterobasidion spp, triggers non-specific defence responses in bark of Norway spruce

Norway spruce [Picea abies (L.) Karst.] is one of the economically most important conifer species in Europe. The major pathogen on Norway spruce is Heterobasidion parviporum (Fr.) Niemela & Korhonen. The completed genome sequence of H. irregulare opens up new possibilities to understand the interactions between Heterobasidion spp. and spruce. However to date there are no completed genome sequencing projects of conifer genomes available for complementing studies. To achieve a better understanding of induced transcriptional defence responses in Norway spruce upon Heterobasidion spp. attack, we compared transcriptional responses in bark to H. parviporum infection to the response to wounding using cDNA-AFLP and transcriptome sequencing. In an initial study bark samples were harvested at 3, 7 and 14 days post inoculation (dpi) and untreated bark was used as negative control. About 2500 transcribed derived fragments (TDFs) generated by cDNA-AFLP were screened. 199 TDFs were investigated further based on band intensity in the inoculated bark in relation to either untreated bark or wounded bark. Out of these, 119 TDFs had a putative homology and a consistent band intensity pattern between replications. A majority of these TDFs showed homology to genes known to associate with defence e.g. 3-deoxy-d-arabino-heptulosonate 7-phosphate synthetase (DAHP), Pathogenesis-related protein 1 (PR1), Lipoxygenase (LOX), ACC-synthase (ACS), ACC-oxidase (ACO) and Jasmonate ZIM-domain 1 (JAZ1). Many of these are found in Salicylic acid- or Jasmonic acid/ethylene-signalling pathways. The majority of the TDFs showed a similar expression pattern for all treatments but samples inoculated with H. parviporum generally showed an enhanced reaction (induction/repression) compared to wounding alone. Expression patterns were confirmed by qPCR in material treated with wounding and inoculation with H. parviporum or Phlebiopsis gigantea. Our data suggest that infection with H. parviporum in Norway spruce induces a broad defence, with many similarities to non-specific defence responses in angiosperms. Additionally signs of reallocation of carbon from primary to secondary metabolism were evident. With this information at hand we analysed four Norway spruce genotypes with either high or low susceptibility to Heterobasidion spp. [1] sampled 0, 5, 15 and 28 dpi with H. annosum. The bark phenol-composition was profiled in each sample. The 500,000 454-reads were assembled into 17,228 contigs that assembled in 14,364 putative transcript units (PTU) using the sequence assembler software Newbler™ (http://www.454.com). The assembled reference file was annotated with the software Blast2Go and the PTUs were submitted to differential expression analysis. Data on associations between gene expression levels and phenol composition in bark upon H. annosum inoculation and level of susceptibility will be presented.