Partial intersterility in Heterobasidion annosum

Introgression between highly divergent fungal sister species

To understand how species evolve and adapt to changing environments, it is important to study gene flow and introgression due to their influence on speciation and radiation events. Here, we apply a novel experimental system for investigating these mechanisms using natural populations. The system is based on two fungal sister species with morphological and ecological similarities occurring in overlapping habitats. We examined introgression between these species by conducting whole genome sequencing of individuals from populations in North America and Europe. We assessed genome-wide nucleotide divergence and performed crossing experiments to study reproductive barriers. We further used ABBA-BABA statistics together with a network analysis to investigate introgression, and conducted demographic modelling to gain insight into divergence times and introgression events. The results revealed that the species are highly divergent and incompatible in vitro. Despite this, small regions of introgression were scattered throughout the genomes and one introgression event likely involves a ghost population (extant or extinct). This study demonstrates that introgression can be found among divergent species and that population histories can be studied without collections of all the populations involved. Moreover, the experimental system is shown to be a useful tool for research on reproductive isolation in natural populations.

Multiple rearrangements and low inter- and intra-species mitogenome sequence variation in the Heterobasidion annosum s.l. species complex

Introduction

Mitochondria are essential organelles in the eukaryotic cells and responsible for the energy production but are also involved in many other functions including virulence of some fungal species. Although the evolution of fungal mitogenomes have been studied at some taxonomic levels there are still many things to be learned from studies of closely related species.

Methods

In this study, we have analyzed 60 mitogenomes in the five species of the Heterobasidion annosum sensu lato complex that all are necrotrophic pathogens on conifers.

Results and Discussion

Compared to other fungal genera the genomic and genetic variation between and within species in the complex was low except for multiple rearrangements. Several translocations of large blocks with core genes have occurred between the five species and rearrangements were frequent in intergenic areas. Mitogenome lengths ranged between 108 878 to 116 176 bp, mostly as a result of intron variation. There was a high degree of homology of introns, homing endonuclease genes, and intergenic ORFs among the five Heterobasidion species. Three intergenic ORFs with unknown function (uORF6, uORF8 and uORF9) were found in all five species and was located in conserved synteny blocks. A 13 bp long GC-containing self-complementary palindrome was discovered in many places in the five species that were optional in presence/absence. The within species variation is very low, among 48 H. parviporum mitogenomes, there was only one single intron exchange, and SNP frequency was 0.28% and indel frequency 0.043%. The overall low variation in the Heterobasidion annosum sensu lato complex suggests a slow evolution of the mitogenome.

Multilocus phylogenies reveal three new truffle-like taxa and the traces of interspecific hybridization in Octaviania (Boletaceae, Boletales)

Among many convergently evolved sequestrate fungal genera in Boletaceae (Boletales, Basidiomycota), the genus Octaviania is the most diverse. We recently collected many specimens of Octaviania subg. Octaviania, including several undescribed taxa, from Japan and the Americas. Here we describe two new species in subgenus Octaviania, O. tenuipes and O. tomentosa, from temperate to subtropical evergreen Fagaceae forests in Japan based on morphological observation and robust multilocus phylogenetic analyses (nrDNA ITS and partial large subunit [LSU], translation elongation factor 1-α gene [TEF1] and the largest subunit of RNA polymerase II gene [RPB1]). Based on specimens from the Americas as well as studies of the holotype, we also taxonomically re-evaluate O. asterosperma var. potteri. Our analysis suggests that O. asterosperma var. potteri is a distinct taxon within the subgenus Octaviania so we recognize this as O. potteri stat. nov. We unexpectedly collected O. potteri specimens from geographically widespread sites in the USA, Japan and Colombia. This is the first verified report of Octaviania from the South American continent. Our molecular analyses also revealed that the RPB1 sequence of one O. tenuipes specimen was identical to that of a closely related species, O. japonimontana, and that one O. potteri specimen from Minnesota had an RPB1 sequence of an unknown species of O. subg. Octaviania. Additionally, one O. japonimontana specimen had an unusually divergent TEF1 sequence. Gene-tree comparison and phylogenetic network analysis of the multilocus dataset suggest that these heterogenous sequences are most likely the result of previous inter- and intra-specific hybridization. We hypothesize that frequent hybridization events in Octaviania may have promoted the high genetic and species diversity found within the genus.

An Updated Global Species Diversity and Phylogeny in the Forest Pathogenic Genus Heterobasidion (Basidiomycota, Russulales)

Heterobasidion species are amongst the most intensively studied polypores because several species are aggressive white rot pathogens of managed coniferous forests mainly in Europe and North America. In the present study, both morphological and multilocus phylogenetic analyses were carried out on Heterobasidion samples from Asia, Oceania, Europe and North America. Three new taxa were found, i.e., H. armandii, H. subinsulare, and H. subparviporum are from Asia and are described as new species. H. ecrustosum is treated as a synonym of H. insulare. So far, six taxa in the H. annosum species complex are recognized. Heterobasidion abietinum, H. annosum, and H. parviporum occur in Europe, H. irregulare, and H. occidentale in North America, and H. subparviporum in East Asia. The North American H. irregulare was introduced to Italy during the Second World War. Species in the H. annosum complex are pathogens of coniferous trees, except H. subparviporum that seems to be a saprotroph. Ten species are found in the H. insulare species complex, all of them are saprotrophs. The pathogenic species are distributed in Europe and North America; the Asian countries should consider the European and North American species as entry plant quarantine fungi. Parallelly, European countries should consider the American H. occidentale and H. irregulare as entry plant quarantine fungi although the latter species is already in Italy, while North America should treat H. abietinum, H. annosum s.s., and H. parviporum as entry plant quarantine fungi. Eight Heterobasidion species found in the Himalayas suggest that the ancestral Heterobasidion species may have occurred in Asia.

Overexpression of PaNAC03, a stress induced NAC gene family transcription factor in Norway spruce leads to reduced flavonol biosynthesis and aberrant embryo development

BACKGROUND

The NAC family of transcription factors is one of the largest gene families of transcription factors in plants and the conifer NAC gene family is at least as large, or possibly larger, as in Arabidopsis. These transcription factors control both developmental and stress induced processes in plants. Yet, conifer NACs controlling stress induced processes has received relatively little attention. This study investigates NAC family transcription factors involved in the responses to the pathogen Heterobasidion annosum (Fr.) Bref. sensu lato.

RESULTS

The phylogeny and domain structure in the NAC proteins can be used to organize functional specificities, several well characterized stress-related NAC proteins are found in III-3 in Arabidopsis (Jensen et al. Biochem J 426:183-196, 2010). The Norway spruce genome contain seven genes with similarity to subgroup III-3 NACs. Based on the expression pattern PaNAC03 was selected for detailed analyses. Norway spruce lines overexpressing PaNAC03 exhibited aberrant embryo development in response to maturation initiation and 482 misregulated genes were identified in proliferating cultures. Three key genes in the flavonoid biosynthesis pathway: a CHS, a F3'H and PaLAR3 were consistently down regulated in the overexpression lines. In accordance, the overexpression lines showed reduced levels of specific flavonoids, suggesting that PaNAC03 act as a repressor of this pathway, possibly by directly interacting with the promoter of the repressed genes. However, transactivation studies of PaNAC03 and PaLAR3 in Nicotiana benthamiana showed that PaNAC03 activated PaLAR3A, suggesting that PaNAC03 does not act as an independent negative regulator of flavan-3-ol production through direct interaction with the target flavonoid biosynthetic genes.

CONCLUSIONS

PaNAC03 and its orthologs form a sister group to well characterized stress-related angiosperm NAC genes and at least PaNAC03 is responsive to biotic stress and appear to act in the control of defence associated secondary metabolite production.

The Role of Hybridization in the Evolution and Emergence of New Fungal Plant Pathogens

Hybridization in fungi has recently been recognized as a major force in the generation of new fungal plant pathogens. These include the grass pathogen Zymoseptoria pseudotritici and the powdery mildew pathogen Blumeria graminis triticale of triticale. Hybridization also plays an important role in the transfer of genetic material between species. This process is termed introgressive hybridization and involves extensive backcrossing between hybrid and the parental species. Introgressive hybridization has contributed substantially to the successful spread of plant pathogens such as Ophiostoma ulmi and O. novo-ulmi, the causal agents of Dutch elm disease, and other tree pathogens such as the rust pathogen Melampsora. Hybridization occurs more readily between species that have previously not coexisted, so-called allopatric species. Reproductive barriers between allopatric species are likely to be more permissive allowing interspecific mating to occur. The bringing together of allopatric species of plant pathogens by global agricultural trade consequently increases the potential for hybridization between pathogen species. In light of global environmental changes, agricultural development, and the facilitated long-distance spread of fungal plant pathogens, hybridization should be considered an important mechanism whereby new pathogens may emerge. Recent studies have gained insight into the genetics and biology of fungal hybrids. Here I summarize current knowledge about hybrid speciation and introgressive hybridization. I propose that future studies will benefit greatly from the availability of large genome data sets and that genome data provide a powerful resource in combination with experimental approaches for analyses of hybrid species.

Selection processes in simple sequence repeats suggest a correlation with their genomic location: insights from a fungal model system

BACKGROUND

Adaptive processes shape the evolution of genomes and the diverse functions of different genomic regions are likely to have an impact on the trajectory and outcome of this evolution. The main underlying hypothesis of this study is that the evolution of Simple Sequence Repeats (SSRs) is correlated with the evolution of the genomic region in which they are located, resulting in differences of motif size, number of repeats, and levels of polymorphisms. These differences should be clearly detectable when analyzing the frequency and type of SSRs within the genome of a species, when studying populations within a species, and when comparing closely related sister taxa. By coupling a genome-wide SSR survey in the genome of the plant pathogenic fungus Heterobasidion irregulare with an analysis of intra- and interspecific variability of 39 SSR markers in five populations of the two sibling species H. irregulare and H. annosum, we investigated mechanisms of evolution of SSRs.

RESULTS

Results showed a clear dominance of trirepeats and a selection against other repeat number, i.e. di- and tetranucleotides, both in regions inside Open Reading Frames (ORFs) and upstream 5' untranslated region (5'UTR). Locus per locus AMOVA showed SSRs both inside ORFs and upstream 5'UTR were more conserved within species compared to SSRs in other genomic regions, suggesting their evolution is constrained by the functions of the regions they are in. Principal coordinates analysis (PCoA) indicated that even if SSRs inside ORFs were less polymorphic than those in intergenic regions, they were more powerful in differentiating species. These findings indicate SSRs evolution undergoes a directional selection pressure comparable to that of the ORFs they interrupt and to that of regions involved in regulatory functions.

CONCLUSIONS

Our work linked the variation and the type of SSRs with regions upstream 5'UTR, putatively harbouring regulatory elements, and shows that the evolution of SSRs might be affected by their location in the genome. Additionally, this study provides a first glimpse on a possible molecular basis for fast adaptation to the environment mediated by SSRs.

Transcriptional Responses Associated with Virulence and Defence in the Interaction between Heterobasidion annosum s.s. and Norway Spruce

Heterobasidion annosum sensu lato is a serious pathogen causing root and stem rot to conifers in the northern hemisphere and rendering the timber defective for sawing and pulping. In this study we applied next-generation sequencing to i) identify transcriptional responses unique to Heterobasidion-inoculated Norway spruce and ii) investigate the H. annosum transcripts to identify putative virulence factors. To address these objectives we wounded or inoculated 30-year-old Norway spruce clones with H. annosum and 454-sequenced the transcriptome of the interaction at 0, 5 and 15 days post inoculation. The 491860 high-quality reads were de novo assembled and the relative expression was analysed. Overall, very few H. annosum transcripts were represented in our dataset. Three delta-12 fatty acid desaturase transcripts and one Clavaminate synthase-like transcript, both associated with virulence in other pathosystems, were found among the significantly induced transcripts. The analysis of the Norway spruce transcriptional responses produced a handful of differentially expressed transcripts. Most of these transcripts originated from genes known to respond to H. annosum. However, three genes that had not previously been reported to respond to H. annosum showed specific induction to inoculation: an oxophytodienoic acid–reductase (OPR), a beta–glucosidase and a germin-like protein (GLP2) gene. Even in a small data set like ours, five novel highly expressed Norway spruce transcripts without significant alignment to any previously annotated protein in Genbank but present in the P. abies (v1.0) gene catalogue were identified. Their expression pattern suggests a role in defence. Therefore a more complete survey of the transcriptional responses in the interactions between Norway spruce and its major pathogen H. annosum would probably provide a better understanding of gymnosperm defence than accumulated until now.

Secondary metabolite comparison of the species within the Heterobasidion annosum s.l. complex

The metabolite production of the five members of the fungal species complex Heterobasidion annosum s.l., i.e. H. annosum s.s., H. abietinum, H. parviporum, H. occidentale and H. irregulare, was analyzed by LC-HRMS. The five members are described to have differences in host preferences: H. annosum s.s. and H. irregulare are pine infecting species, and H. parviporum, H. occidentale and H. abietinum are non-pine infecting. Principal component analysis (PCA) of the LC-HRMS data showed that samples from the five species could be separated into five groups and in accordance with the differences in host preferences. Twenty-three compounds, important to the observed PCA grouping, were isolated and identified. The main contributor to the separation of the pine infecting species from the non-pine infecting species in PC1 was the benzohydrofuran fomannoxin, which was only detected in the pine infecting species H. annosum s.s. and H. irregulare. These two species were further separated in PC3, and one major contributor here was the sesquiterpene deoxyfomannosin A. The three non-pine infecting species were separated in PC2, by epoxydrimenol that was detected in only two of the species and further in PC4, where a few fomannoxin related compounds were important for the grouping. During the work, three unknown compounds were isolated and described: 3-hydroxy-2-(1,3-dihydroxypropan-2-yl)-2,3-dihydrobenzofuran-5-carbaldehyde, 3-hydroxy-2-(1,2,3-trihydroxypropan-2-yl)-2,3-dihydrobenzofuran-5-carbaldehyde and 3-hydroxy-2,3-dihydrobenzofuran-5-carboxylic acid.

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.

Population genetic analyses provide insights on the introduction pathway and spread patterns of the North American forest pathogen Heterobasidion irregulare in Italy

A population genetics approach is used to identify the most likely introduction site and introduction pathway for the North American forest pathogen Heterobasidion irregulare using 101 isolates from six sites in Italy and 34 isolates from five sites in North America. Diversity indices based on sequences from ten loci indicate the highest diversity in Italy is found in Castelfusano/Castelporziano and that diversity progressively decreases with increasing distance from that site. AMOVA, Bayesian clustering and principal coordinates analyses based on 12 SSR loci indicate high levels of gene flow among sites, high frequency of admixing, and fail to identify groups of genotypes exclusive to single locations. Cumulatively, these analyses suggest the current infestation is the result of multiple genotypes expanding their range from a single site. Based on two sequenced loci, a single source site in North America could provide enough variability to explain the variability observed in Italy. These results support the notion that H. irregulare was introduced originally in Castelporziano: because Castelporziano has been sealed off from the rest of the world for centuries except for a camp set up by the US military in 1944, we conclude the fungus may have been transported in infected wood used by the military. Finally, spatial autocorrelation analyses using SSR data indicate a significant under-dispersion of alleles up to 0.5-10 km, while a significant overdispersion of alleles was detected at distances over 80 km: these ranges can be used to make predictions on the likely dispersal potential of the invasive pathogen.

Extensive trans-specific polymorphism at the mating type locus of the root decay fungus Heterobasidion

Incompatibility systems in which individuals bearing identical alleles reject each other favor the maintenance of a diversity of alleles. Mushroom mating type loci (MAT) encode for dozens or hundreds of incompatibility alleles whose loss from the population is greatly restricted through negative frequency selection, leading to a system of alleles with highly divergent sequences. Here, we use DNA sequences of homeodomain (HD) encoding genes at the MAT locus of five closely related species of the root rot basidiomycete Heterobasidion annosum sensu lato to show that the extended coalescence time of MAT alleles greatly predates speciation in the group, contrasting loci outside of MAT that show allele divergences largely consistent with the species phylogeny with those of MAT, which show rampant trans-species polymorphism. We observe a roughly 6-fold greater genealogical depth and polymorphism of MAT compared with non-MAT that argues for the maintenance of balanced polymorphism for a minimum duration of 24 My based on a molecular-clock calibrated species phylogeny. As with other basidiomycete HD genes, balancing selection appears to be concentrated at the specificity-determining region in the N-terminus of the protein based on identification of codons under selection and the absence of recombination within the region. However, the elevated polymorphism extends into the nonspecificity determining regions as well as a neighboring non-MAT gene, the mitochondrial intermediate peptidase (MIP). In doing so, increased divergence should decrease recombination among alleles and as a by-product create incompatibilities in the functional domains not involved in allele recognition but in regulating sexual development.

Biology, epidemiology, and control of Heterobasidion species worldwide

Heterobasidion annosum sensu lato is a species complex comprising five species that are widely distributed in coniferous forests of the Northern Hemisphere and are each characterized by a distinct host preference. More than 1,700 papers have been published on these fungi in the past four decades, making them perhaps the most widely studied forest fungi. Heterobasidion species are at different levels on the saprotroph-necrotroph gradient, and the same individual can switch from one mode to the other. This offers a unique opportunity to study how genomic structure, gene expression, and genetic trade-offs may all interact with environmental factors to determine the life mode of the organism. The abilities of Heterobasidion spp. to infect stumps as saprotrophs and to spread to neighboring trees as pathogens have resulted in significant damages to timber production in managed forests. This review focuses on the current knowledge of the biology, ecology, evolution, and management of these species and is based on classical and modern studies.

Chemical and transcriptional responses of Norway spruce genotypes with different susceptibility to Heterobasidion spp. infection

BACKGROUND

Norway spruce [Picea abies (L.) Karst.] is one of the most important conifer species in Europe. The wood is economically important and infections by wood-rotting fungi cause substantial losses to the industry.The first line of defence in a Norway spruce tree is the bark. It is a very efficient barrier against infection based on its mechanical and chemical properties. Once an injury or an infection is recognized by the tree, induced defences are activated. In this study we examined transcriptional response, using 454-sequencing, and chemical profiles in bark of Norway spruce trees with different susceptibility to Heterobasidion annosum s.l. infection. The aim was to find associations between the transcriptome and chemical profiles to the level of susceptibility to Heterobasidion spp. in Norway spruce genotypes.

RESULTS

Both terpene and phenol compositions were analysed and at 28 days post inoculation (dpi) high levels of 3-carene was produced in response to H. annosum. However, significant patterns relating to inoculation or to genotypes with higher or lower susceptibility could only be found in the phenol fraction. The levels of the flavonoid catechin, which is polymerized into proanthocyanidins (PA), showed a temporal variation; it accumulated between 5 and 15 dpi in response to H. annosum infection in the less susceptible genotypes. The transcriptome data suggested that the accumulation of free catechin was preceded by an induction of genes in the flavonoid and PA biosynthesis pathway such as leucoanthocyanidin reductase. Quantitative PCR analyses verified the induction of genes in the phenylpropanoid and flavonoid pathway. The qPCR data also highlighted genotype-dependent differences in the transcriptional regulation of these pathways.

CONCLUSIONS

The varying dynamics in transcriptional and chemical patterns displayed by the less susceptible genotypes suggest that there is a genotypic variation in successful spruce defence strategies against Heterobasidion. However, both high levels of piceasides and flavonoids in the less susceptible genotypes suggested the importance of the phenolic compounds in the defence. Clearly an extended comparison of the transcriptional responses in the interaction with Heterobasidion between several independent genotypes exhibiting reduced susceptibility is needed to catalogue mechanisms of successful host defence strategies.

A novel putative partitivirus of the saprotrophic fungus Heterobasidion ecrustosum infects pathogenic species of the Heterobasidion annosum complex

We characterized the bisegmented genome of a putative double-stranded (ds) RNA virus from a Chinese isolate of the fungus Heterobasidion ecrustosum, a member of the Heterobasidion insulare species complex. The larger genomic segment of 1885bp encoded a putative RNA dependent RNA polymerase (RdRp, 585aa), and the smaller one for a putative coat protein of 521aa (1826bp). Phylogenetic analyses suggest that this novel virus species, named as 'Heterobasidion RNA virus 3 from H. ecrustosum, strain 1' (HetRV3-ec1), can be assigned to the family Partitiviridae, being most similar to the Helicobasidium mompa dsRNA mycovirus with RdRp amino acid similarity of 54%. The similarity to known viruses of other Heterobasidion species was notably low (25-39%). The virus could be experimentally transmitted to members of the Heterobasidion annosum complex: the European Heterobasidion abietinum and North American Heterobasidion occidentale, and the original host strain could be cured from the virus by thermal treatment. Microscopical observations showed that hyphae of H. ecrustosum anastomosed occasionally with H. abietinum and H. occidentale, and suggested a possible route for horizontal transmission between these sexually incompatible species. The virus infection seemed to cause variable effects on the growth rate of its fungal hosts, but the results were strongly dependent on fungal strain, growth medium and incubation temperature.

Heterobasidion occidentale sp. nov. and Heterobasidion irregulare nom. nov.: a disposition of North American Heterobasidion biological species

The genus Heterobasidion includes some of the most important pathogens of conifers in the world, and as such it is one of the most intensely studied genera of fungi. Because of the remarkable paucity of distinguishing morphological traits, the taxonomy of species within this genus has always been problematic. A partial resolution of the taxonomic issues regarding this genus was achieved by defining the most important and first described species within it, Heterobasidion annosum, as a species complex containing at least two partially intersterile biological species defined as intersterility groups (ISGs). With time, the number of ISGs has increased to include at least two distinct North American and three distinct Eurasian ISGs. Two additional, yet unnamed, taxonomic groups within Heterobasidion have been recently described in Japan. ISGs are distinguishable either by minor morphological differences, by partial intersterility, by ecological traits including host preference, and/or by their geographic range. Several studies employing a variety of molecular tools and analyses have confirmed the distinct genetic divergence among ISGs, identifying each of them as a monophyletic group. Using genetic markers, genotypes can always be unambiguously assigned to one ISG, and very few inter-ISG hybrids have been identified. In this paper, we summarize the available information, both genetic and ecological, that differentiates the two North American ISGs from each other and from other taxonomic units within the genus. We demonstrate that morphometric characteristics such as pore density and pore shape differentiate the two ISGs. Based on the cumulative genetic, ecological, and morphological evidence, we propose a disposition of ISGs of the North American H. annosum by replacing the P ISG with Heterobasidion irregulare, and the S ISG with Heterobasidion occidentale.

Existence of a pattern of reproductive character displacement in Homobasidiomycota but not in Ascomycota

Generally, stronger reproductive isolation is expected between sympatric than between allopatric sibling species. Such reproductive character displacement should predominantly affect premating reproductive isolation and can be due to several mechanisms, including population extinction, fusion of insufficiently isolated incipient species and reinforcement of reproductive isolation in response to low hybrid fitness. Experimental data on several taxa have confirmed these theoretical expectations on reproductive character displacement, but they are restricted to animals and a few plants. Using results reported in the literature on crossing experiments in fungi, we compared the degree and the nature of reproductive isolation between allopatric and sympatric species pairs. In accordance with theoretical expectations, we found a pattern of enhanced premating isolation among sympatric sibling species in Homobasidiomycota. By contrast, we did not find evidence for reproductive character displacement in Ascomycota at similar genetic distances. Both allopatric and sympatric species of Ascomycota had similarly low levels of reproductive isolation, being mostly post-zygotic. This suggests that some phylogeny-dependent life-history trait may strongly influence the evolution of reproductive isolation between closely related species. A significant correlation was found between degree of reproductive isolation and genetic divergence among allopatric species of Homobasidiomycota, but not among sympatric ones or among Ascomycota species.

Inferences on the phylogeography of the fungal pathogen Heterobasidion annosum, including evidence of interspecific horizontal genetic transfer and of human-mediated, long-range dispersal

Fungi in the basidiomycete species complex Heterobasidion annosum are significant root-rot pathogens of conifers throughout the northern hemisphere. We utilize a multilocus phylogenetic approach to examine hypotheses regarding the evolution and divergence of two Heterobasidion taxa associated with pines: the Eurasian H. annosum sensu stricto and the North American H. annosum P intersterility group (ISG). Using DNA sequence information from portions of two nuclear and two mitochondrial loci, we infer phylogenetic relationships via parsimony, Bayesian and median-joining network analysis. Analysis of isolates representative of the entire known geographic range of the two taxa results in monophyletic sister Eurasian and North American lineages, with North America further subdivided into eastern and western clades. Genetically anomalous isolates from the Italian presidential estate of Castelporziano are always part of a North American clade and group with eastern North America, upholding the hypothesis of recent, anthropogenically mediated dispersal. P ISG isolates from Mexico have phylogenetic affinity with both eastern and western North America. Results for an insertion in the mitochondrial rDNA suggest this molecule was obtained from the Heterobasidion S ISG, a taxon sympatric with the P ISG in western North America. These data are compatible with an eastern Eurasian origin of the species, followed by dispersal of two sister taxa into western Eurasia and into eastern North America over a Beringean land bridge, a pattern echoed in the phylogeography of other conifer-associated basidiomycetes.

The origins of plant pathogens in agro-ecosystems

Plant pathogens can emerge in agricultural ecosystems through several mechanisms, including host-tracking, host jumps, hybridization and horizontal gene transfer. High-throughput DNA sequencing coupled with new analytical approaches make it possible to differentiate among these mechanisms and to infer the time and place where pathogens first emerged. We present several examples to illustrate the different mechanisms and timescales associated with the origins of important plant pathogens. In some cases pathogens were domesticated along with their hosts during the invention of agriculture approximately 10,000 years ago. In other cases pathogens appear to have emerged very recently and almost instantaneously following horizontal gene transfer or hybridization. The predominant unifying feature in these examples is the environmental and genetic uniformity of the agricultural ecosystem in which the pathogens emerged. We conclude that agro-ecosystems will continue to select for new pathogens unless they are re-engineered to make them less conducive to pathogen emergence.

Invasion of European pine stands by a North American forest pathogen and its hybridization with a native interfertile taxon

It was recently reported that North American (NA) individuals of the forest pathogen Heterobasidion annosum were found in a single pine stand near Rome, in association with the movement of US troops during World War II. Here, we report on some aspects of the invasion biology of this pathogen in Italian coastal pinewoods, and on its interaction with native (EU) Heterobasidion populations. Spores of Heterobasidion were sampled using woody traps in pine stands along 280 km of coast around Rome. DNA of single-spore colonies was characterized by two sets of nuclear and one set of mitochondrial taxon-specific polymerase chain reaction primers. NA spores were found not only in a single site, but in many locations over a wide geographic area. Invasion occurred at an estimated rate of 1.3 km/year through invasion corridors provided by single trees, and not necessarily by sizable patches of forests. Within the 100-km long range of expansion, the NA taxon was dominant in all pure pine stands. Because abundance of the EU taxon is low and identical among stands within and outside the area invaded by NA individuals, we infer that the exotic population has invaded habitats mostly unoccupied by the native species. Discrepancy between a mitochondrial and a nuclear marker occurred in 3.8% of spores from one site, a mixed oak-pine forest where both taxa were equally represented. Combined phylogenetic analyses on nuclear and mitochondrial loci confirmed these isolates were recombinant. The finding of hybrids indicates that genetic interaction between NA and EU Heterobasidion taxa is occurring as a result of their current sympatry.

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.

An AFLP-markers based genetic linkage map of Heterobasidion annosum locating intersterility genes

A genetic linkage map of the basidiomycete Heterobasidion annosum, casual agent of root rot in conifers, was constructed from a compatible mating between isolates from the North American S and P intersterility groups. In a population consisting of 102 progeny isolates, 358 AFLP markers were scored. The linkage analysis generated 19 large linkage groups, containing 6 or more markers, which covered 1468 cM. The physical size to genetic distance was approximately 11.1 kbp/cM. Segregation of three intersterility gene loci were analysed through mating of the progeny isolates with three tester strains carrying known intersterility genotypes. The loci for the two intersterility genes (S and P) were successfully located in the map. Segregation of the mating type locus was analysed by backcrossing the progeny isolates with their parental strains. The mating type locus could not be located in the map.

Mechanisms of fungal speciation

The objective of this review is to provide a synthesis of speciation theory, of what is known about mechanisms of speciation in fungi and from this, what is expected, and of ideas on how speciation can be elucidated in more fungal systems. The emphasis is on process rather than pattern. Phylogeographic studies in some groups, such as the agarics, demonstrate predominantly allopatric speciation, often through vicariance, as seen in many plants and animals. The variety of life history factors in fungi suggests, however, a diversity in speciation mechanisms that is borne out in comparison of some key examples. Life history features in fungi with a bearing on speciation include genetic mechanisms for intra- and interspecies interactions, haploidy as monokaryons, dikaryons, or coenocytes, distinctive types of propagules with distinctive modes of dispersal, as well as characteristic relationships to the substrate or host as specialized or generalist saprotrophs, parasites or mutualists with associated opportunities and selective pressures for hybridization. Approaches are proposed for both retrospective, phylogeographic determination of speciation mechanisms, and experimental studies with the potential for genomic applications, particularly in examining the relationship between adaptation and reproductive isolation.

Reproductive isolation and phylogenetic divergence in Neurospora: comparing methods of species recognition in a model eukaryote

We critically examined methods for recognizing species in the model filamentous fungal genus Neurospora by comparing traditional biological species recognition (BSR) with more comprehensive applications of both BSR and phylogenetic species recognition (PSR). Comprehensive BSR was applied to a set of 73 individuals by performing extensive crossing experiments and delineating biological species based on patterns of reproductive success. Within what were originally considered two species, N. crassa and N. intermedia, we recognized four reproductively isolated biological species. In a concurrent study (Dettman et al. 2003), we used genealogical concordance of four independent nuclear loci to recognize phylogenetic species in Neurospora. Overall, the groups of individuals identified as species were similar whether recognized by reproductive success or by phylogenetic criteria, and increased genetic distance between parents was associated with decreased reproductive success of crosses, suggesting that PSR using genealogical concordance can be used to reliably recognize species in organisms that are not candidates for BSR. In one case, two phylogenetic species were recognized as a single biological species, indicating that significant phylogenetic divergence preceded the development of reproductive isolation. However, multiple biological species were never recognized as a single phylogenetic species. Each of the putative N. crassa x N. intermedia hybrids included in this study was confidently assigned to a single species, using both PSR and BSR. As such, no evidence for a history of hybridization in nature among Neurospora species was observed. By performing reciprocal mating tests, we found that mating type, parental role, and species identity of parental individuals could all influence the reproductive success of matings. We also observed sympatry-associated sexual dysfunction in interspecific crosses, which was consistent with the existence of reinforcement mechanisms.

Molecular markers reveal genetic isolation and phylogeography of the S and F intersterility groups of the wood-decay fungus Heterobasidion annosum

The root-rot fungus Heterobasidion annosum (Fr.) Bref. species complex consists of three intersterility groups (S, F, and P), separated by their host affinity. The phylogenetic relationship of the species complex was studied, with the focus on the S and F groups, by comparing DNA sequences of four nuclear gene fragments: calmodulin, glyceraldehyde 3-phosphate dehydrogenase, heat stress protein 80-1, and elongation factor 1-alpha, and one anonymous locus, from 29 fungal isolates originating from Europe, Asia, and North America. The phylogeny of each separate gene locus as well as the combined dataset consisted of three main clades: European F group isolates, Euroasian S group isolates, and North American S group isolates, suggesting them to be separated into phylogenetic species. The results also support the hypothesis of an early separation between the S and F groups, indicating that their distribution have followed their host tree species for a considerable time period.

Interspecific hybridization in plant‐associated fungi and oomycetes: a review

Fungi (kingdom Mycota) and oomycetes (kingdom Stramenopila, phylum Oomycota) are crucially important in the nutrient cycles of the world. Their interactions with plants sometimes benefit and sometimes act to the detriment of humans. Many fungi establish ecologically vital mutualisms, such as in mycorrhizal fungi that enhance nutrient acquisition, and endophytes that combat insects and other herbivores. Other fungi and many oomycetes are plant pathogens that devastate natural and agricultural populations of plant species. Studies of fungal and oomycete evolution were extraordinarily difficult until the advent of molecular phylogenetics. Over the past decade, researchers applying these new tools to fungi and oomycetes have made astounding new discoveries, among which is the potential for interspecific hybridization. Consequences of hybridization among pathogens include adaptation to new niches such as new host species, and increased or decreased virulence. Hybrid mutualists may also be better adapted to new hosts and can provide greater or more diverse benefits to host plants.

Relative abundance and potential dispersal range of intersterility groups of Heterobasidion annosum in pure and mixed forests

In Europe the forest pathogen Heterobasidion annosum (Fr.) Bref. includes the S, P, and F intersterility groups (ISGs), each displaying a preferential specialization on Norway spruce (Picea abies (L.) Karst.), pine, and silver fir (Abies alba Mill.), respectively. In this paper, we present data about (i) H. annosum ISGs frequency in different forest types, (ii) the degree of host specificity of each ISG, (iii) the significance of the potential movement of airborne spores among forests, and (iv) the occurrence of S–P chimeras in the northwestern Alps. Using woody spore traps, we sampled natural pure spruce and fir forests and a mixed spruce-fir forest. The ISG of 582 spores was determined by ISG-diagnostic taxon-specific competitive priming (TSCP) polymerase chain reaction (PCR) combined with PCR-mediated detection of ISG-specific introns in the ML5–ML6 DNA region of the mitochondrial large ribosomal RNA (mt LrRNA). All three ISGs were found, and a strong correlation was observed between the F ISG and fir and the S ISG and spruce. In the mixed forest, no clear relationship between tree host species and host-specialized ISGs was found. In spite of a relative dominance of fir in the overstory of the mixed stand, the fir-associated F ISG represented only 11% of the total number of spores collected. This discrepancy was explained by the recent establishment of firs at this site. No S–P nuclear-mitochondrial chimeras were found. This suggests limited gene flow between these ISGs.Key words: Heterobasidion annosum, host specificity, ISGs, gene flow, PCR, Alps.

Genetic Structure of Heterobasidion annosum in White Fir Mortality Centers in California

ABSTRACT The structure of Heterobasidion annosum populations was studied in 15 mixed-conifer sites in central and northern California. Study sites displayed mortality of white fir trees in enlarging discrete patches (mortality centers). At each site, fungal genotypes were defined by somatic compatibility tests. In two sites, further genetic and molecular analyses were performed on field genotypes and on homokaryons obtained by dedikaryotization of field heterokaryons. Isolates were found to be colonizing mostly the roots and the bole sapwood of white fir trees, and no significant infections of other tree species were observed. Each mortality center was characterized by the presence of several fungal genotypes, all belonging to the S intersterility group. Both homokaryotic and heterokaryotic strains were present in all sites. Multiple genotypes were retrieved in individual trees or stumps. Out of 228 fungal genotypes, 86% were found only within a single tree or stump, while 14% had spread to adjacent trees. The two largest genotypes had diameters of 9 and 10 m, and had colonized five and nine trees, stumps, or both, respectively. The maximum distance between two adjacent trees colonized by the same genotype was 6 m, and a highly significant correlation was found between tree diameter and distance of fungal "vegetative" spread. The largest clones were found in areas characterized by high tree and stump densities, and secondary spread of the fungus was more significant in denser stands. In most cases, original infection courts of existing genotypes could be traced to standing trees and not to stumps. The genetic analysis performed in two mortality centers revealed that most local genotypes had different mating alleles, and thus originated from unrelated basidiospores. In a few cases, the same mating allele was shared by two heterokaryons (n+n genome) or by a homokaryon (n genome) and a heterokaryon. Molecular analysis showed that nuclei bearing the same mating allele were identical, providing evidence that the two nuclei forming heterokaryons can act independently in the field and can be shared among isolates, presumably via di-mon mating or by separate matings of different portions of widespread homokaryons.

Distribution patterns of number of nuclei in conidia from heterokaryons of Heterobasidion annosum (Fr.) Bref. and their interpretation in terms of genomic conflict

Ratios of nuclear genotypes observed in conidia from heterokaryotic strains of Heterobasidion annosum (Fr.) Bref., obtained from pairings between sympatrically derived, sib-related and non-sib-related homokaryons, commonly deviated from 1:1. Ratios were temporally stable, and the genotypes examined could be ranked in a strict dominance hierarchy, linked both to the relatedness of the association partners and to the growth rates of the parental homokaryons. Parental homokaryons and sibrelated heterokaryons produced conidia with a mean number of nuclei of about two, whereas non-sib-related-d heterokaryons produced conidia that were predominantly uninucleate. Moreover, whereas conidia containing more than one nucleus germinated most rapidly when derived from homokaryons or sib-related heterokaryons, uninucleate conidia germinated more readily if derived from non-sib-related heterokaryons. In a study of naturally occurring heterokaryons, distribution patterns of the number of nuclei in conidia were found to be similar to those of the homokaryons. The possible interpretation of these findings in terms of inter-nuclear genomic conflict is discussed.

Heterobasidion annosum 5.8s ribosomal DNA and internal transcribed spacer sequence: rapid identification of European intersterility groups by ribosomal DNA restriction polymorphism

Using conserved fungal ribosomal gene sequences the internal transcribed spacer (ITS) regions one and two (ITS1, ITS2) and the 5.8s ribosomal RNA gene (rRNA) of Heterobasidion annosum were amplified by the polymerase chain reaction (PCR). The nucleotide sequence was determined in three European intersterility groups (ISG-S, -F and -P). Three sequence variants of the ITS were found in ISG-S isolates. The sequence of the ITS of ISG-F differed by two residues from the major ISG-S sequence variant. The ISG-P sequence differed from ISG-S and ISG-F at 15-16 and 16 residues, respectively. Amplified intergenic spacer elements were informative for ISG fingerprinting following digestion with various 4-cutter restriction endonucleases. All differences in the restriction fragments between the ISGs were because of sequence differences in the ITS regions. The fingerprint patterns of isolates from the same intersterility group but from different European localities were identical. These results show that ribosomal DNA fingerprinting is a rapid technique to identify ISGs in Heterobasidion annosum.