A genetic linkage map for the ectomycorrhizal fungus Laccaria bicolor and its alignment to the whole-genome sequence assemblies

Identification and characterization of eight metallothionein genes involved in heavy metal tolerance from the ectomycorrhizal fungus Laccaria bicolor

Metallothioneins (MTs) are small, cysteine-rich, heavy metal-binding proteins involved in metal homeostasis and detoxification. The increasing numbers of available genomic sequences of ectomycorrhizal (ECM) fungi enable deeper insights into the characteristics of MT genes in these fungi that form the most important symbiosis with the host trees in forest ecosystems. The aim of this study was to establish a comprehensive, genome-wide inventory of MT genes from the ECM fungus Laccaria bicolor. Eight MT genes in L. bicolor were cloned, and the expression patterns of their transcripts at various developmental stages based on expressed sequence tag (EST) counts were analyzed. The expression levels of four MTs were significantly increased during symbiosis stages. Quantitative real-time PCR (qRT-PCR) analysis revealed that transcripts of LbMT1 were dominant in free-living mycelia and strongly induced by excessive copper (Cu), cadmium (Cd), and hydrogen peroxide (H2O2). To determine whether these eight MTs functioned as metal chelators, we expressed them in the Cu- and Cd-sensitive yeast mutants, cup1∆ and yap1∆, respectively. All LbMT proteins provided similar levels of Cu(II) or Cd(II) tolerance, but did not affect by H2O2. Our findings provide novel data on the evolution and diversification of fungal MT gene duplicates, a valuable resource for understanding the vast array of biological processes in which these proteins are involved.

The Assembled and Annotated Genome of the Fairy-Ring Fungus Marasmius oreades

Marasmius oreades is a basidiomycete fungus that grows in so called "fairy rings," which are circular, underground mycelia common in lawns across temperate areas of the world. Fairy rings can be thought of as natural, long-term evolutionary experiments. As each ring has a common origin and expands radially outwards over many years, different sectors will independently accumulate mutations during growth. The genotype can be followed to the next generation, as mushrooms producing the sexual spores are formed seasonally at the edge of the ring. Here, we present new genomic data from 95 single-spore isolates of the species, which we used to construct a genetic linkage map and an updated version of the genome assembly. The 44-Mb assembly was anchored to 11 linkage groups, producing chromosome-length scaffolds. Gene annotation revealed 13,891 genes, 55% of which contained a pfam domain. The repetitive fraction of the genome was 22%, and dominated by retrotransposons and DNA elements of the KDZ and Plavaka groups. The level of assembly contiguity we present is so far rare in mushroom-forming fungi, and we expect studies of genomics, transposons, phylogenetics, and evolution to be facilitated by the data we present here of the iconic fairy-ring mushroom.

Genetic linkage map construction and quantitative trait loci mapping of agronomic traits in Gloeostereum incarnatum

Gloeostereum incarnatum is an edible medicinal mushroom widely grown in China. Using the whole genome of G. incarnatum, simple sequence repeat (SSR) markers were developed and synthetic primers were designed to construct its first genetic linkage map. The 1,048.6 cm map is composed of 10 linkage groups and contains 183 SSR markers. In total, 112 genome assembly sequences were anchored, representing 16.43 Mb and covering 46.41% of the genome. Selfing populations were used for quantitative trait loci (QTL) targeting, and the composite interval mapping method was used to co-localize the mycelium growth rate (potato dextrose agar and sawdust), growth period, yield and fruiting body length, and width and thickness. The 14 QTLs of agronomic traits had LOD values of 3.20-6.51 and contribution rates of 2.22-13.18%. No linkage relationship was found between the mycelium growth rate and the growth period, but a linkage relationship was observed among the length, width and thickness of the fruiting bodies. Using NCBI's BLAST alignment, the genomic sequences corresponding to the QTL regions were compared, and a TPR-like protein candidate gene was selected. Using whole-genome data, 138 candidate genes were found in four sequence fragments of two SSR markers located in the same scaffold. The genetic map and QTLs established in this study will aid in developing selective markers for agronomic traits and identifying corresponding genes, thereby providing a scientific basis for the further gene mapping of quantitative traits and the marker-assisted selection of functional genes in G. incarnatum breeding programs.

Subchromosome-Scale Nuclear and Complete Mitochondrial Genome Characteristics of Morchella crassipes

Morchella crassipes (Vent.) Pers., a typical yellow morel species with high economic value, is mainly distributed in the low altitude plains of Eurasia. However, rare research has been performed on its genomics and polarity, thus limiting its research and development. Here, we reported a fine physical map of the nuclear genome at the subchromosomal-scale and the complete mitochondrial genome of M. crassipes. The complete size of the nuclear genome was 56.7 Mb, and 23 scaffolds were assembled, with eight of them being complete chromosomes. A total of 11,565 encoding proteins were predicted. The divergence time analysis showed that M. crassipes representing yellow morels differentiated with black morels at ~33.98 Mya (million years), with 150 gene families contracted and expanded in M. crassipes versus the two black morels (M. snyderi and M. importuna). Furthermore, 409 CAZYme genes were annotated in M. crassipes, containing almost all plant cell wall degrading enzymes compared with the mycorrhizal fungi (truffles). Genomic annotation of mating type loci and amplification of the mating genes in the monospore population was conducted, the results indicated that M. crassipes is a heterothallic fungus. Additionally, a complete circular mitochondrial genome of M. crassipes was assembled, the size reached as large as 531,195 bp. It can be observed that the strikingly large size was the biggest up till now, coupled with 14 core conserved mitochondrial protein-coding genes, two rRNAs, 31 tRNAs, 51 introns, and 412 ncORFs. The total length of intron sequences accounted for 53.67% of the mitochondrial genome, with 19 introns having a length over 5 kb. Particularly, 221 of 412 ncORFs were distributed within 51 introns, and the total length of the ncORFs sequence accounted for 40.83% of the mitochondrial genome, and 297 ncORFs had expression activity in the mycelium stage, suggesting their potential functions in M. crassipes. Meanwhile, there was a high degree of repetition (51.31%) in the mitochondria of M. crassipes. Thus, the large number of introns, ncORFs and internal repeat sequences may contribute jointly to the largest fungal mitochondrial genome to date. The fine physical maps of nuclear genome and mitochondrial genome obtained in this study will open a new door for better understanding of the mysterious species of M. crassipes.

Construction of a genetic linkage map of Lentinula edodes based on SSR, SRAP and TRAP markers

Genetic mapping is a basic tool for eukaryotic genomic research. It allows the localization of genes or quantitative trait loci (QTLs) and map-based cloning. In this study, we constructed a linkage map based on DNA samples from a commercial strain L808, including two parental monokaryons and 93 single spore isolates considered with segregating to 1:1:1:1 at four mating types (A₁B₁, A₁B₂, A₂B₁ and A₂B₂). Using Simple Sequence Repeats (SSR), Sequence Related Amplified Polymorphism (SRAP), Target Region Amplified Polymorphism (TRAP) molecular markers, 182 molecular markers and two mating factors were located on 11 linkage groups (LGs). The total length of the map was 948.083 centimorgan (cM), with an average marker interval distance of 4.817 cM. Only two gaps spanning more than 20 cM was observed. The probability of 20 cM, 10 cM, 5 cM genetic distance cover one marker was 99.68%, 94.36%, 76.43% in our genetic linkage map, respectively. This is the first linkage map of Lentinula edodes using SSR markers, which provides essential information for quantitative trait analyses and improvement of genome assembly.

Data on the genome analysis of the wild edible mushroom, Russula griseocarnosa

In the present article, we report data on the whole genome sequence of a wild edible and medicinal ectomycorrhizal fungus Russula griseocarnosa. The R. griseocarnosa genome consists of 64.81 Mb with a GC-pair content of 49.41%. The genome assembly consists of 471 scaffolds and 16128 coding protein genes. The coding protein genes was annotated in different databases (GO, KEGG and CAZys), respectively. The whole genome sequence and functional annotation provide important information for ectomycorrhizal fungus, which can be used as a basis for cultivation and breeding of R. griseocarnosa. The Whole Genome project of Russula griseocarnosa has been deposited at DDBJ/ENA/GenBank under the accession RMVF00000000. The version described is RMVF01000000. To further interpretation of the data provided in this article, please refer to the research article 'Whole genome sequencing and genome annotation of the wild edible mushroom, Russula griseocarnosa' [1].

Whole genome sequencing and genome annotation of the wild edible mushroom, Russula griseocarnosa

Russula griseocarnosa is a species of edible ectomycorrhizal fungi with medicinal properties that grows in southern China. Total DNA was isolated from a fresh fruiting body of R. griseocarnosa and subjected to sequencing using Illumina Hiseq with the PacBio RS sequencing platform. Here, we present the 64.81 Mb draft genome map of R. griseocarnosa based on 471 scaffolds and 16,128 coding protein genes. The gene annotation of protein coding genes was used to obtain corresponding annotations by blastp. Phylogenetic analysis revealed a close evolutionary relationship of R. griseocarnosa to Heterobasidion irregulare and Stereum hirsutum in the core Russulales clade. The R. griseocarnosa genome encodes a repertoire of enzymes engaged in carbohydrate and polysaccharide metabolism, along with cytochrome P450s and secondary metabolite biosynthesis. The genome content of R. griseocarnosa provides insights into the genetic basis of its reported medicinal properties and serves as a reference for comparative genomics of fungi.

Construction of a genetic linkage map and QTL mapping of agronomic traits in Auricularia auricula-judae

Auricularia auricula-judae is a traditional edible fungus that is cultivated widely in China. In this study, a genetic linkage map for A. auricula-judae was constructed using a mapping population consisting of 138 monokaryons derived from a hybrid strain (A119-5). The monokaryotic parent strains A14-5 and A18-119 were derived from two cultivated varieties, A14 (Qihei No. 1) and A18 (Qihei No. 2), respectively. In total, 130 simple sequence repeat markers were mapped. These markers were developed using the whole genome sequence of A. auricula-judae and amplified in A14-5, A18- 119, and the mapping population. The map consisted of 11 linkage groups (LGs) spanning 854 cM, with an average interval length of 6.57 cM. A testcross population was derived from crossing between the monokaryon A184-57 (from the wild strain A184 as a tester strain) and the mapping population. Important agronomic trait-related QTLs, including mycelium growth rate on potato dextrose agar for the mapping population, mycelium growth rate on potato dextrose agar and sawdust for the testcross population, growth period (days from inoculation to fruiting body harvesting), and yield for the testcross population, were identified using the composite interval mapping method. Six mycelium growth raterelated QTLs were identified on LG1 and LG4, two growth period-related QTLs were identified on LG2, and three yieldrelated QTLs were identified on LG2 and LG6. The results showed no linkage relationship between mycelium growth rate and growth period. The present study provides a foundation for locating genes for important agronomic characteristics in A. auricula-judae in the future.

Fungal Diversity in Field Mold-Damaged Soybean Fruits and Pathogenicity Identification Based on High-Throughput rDNA Sequencing

Continuous rain and an abnormally wet climate during harvest can easily lead to soybean plants being damaged by field mold (FM), which can reduce seed yield and quality. However, to date, the underlying pathogen and its resistance mechanism have remained unclear. The objective of the present study was to investigate the fungal diversity of various soybean varieties and to identify and confirm the FM pathogenic fungi. A total of 62,382 fungal ITS1 sequences clustered into 164 operational taxonomic units (OTUs) with 97% sequence similarity; 69 taxa were recovered from the samples by internal transcribed spacer (ITS) region sequencing. The fungal community compositions differed among the tested soybeans, with 42 OTUs being amplified from all varieties. The quadratic relationships between fungal diversity and organ-specific mildew indexes were analyzed, confirming that mildew on soybean pods can mitigate FM damage to the seeds. In addition, four potentially pathogenic fungi were isolated from FM-damaged soybean fruits; morphological and molecular identification confirmed these fungi as Aspergillus flavus, A. niger, Fusarium moniliforme, and Penicillium chrysogenum. Further re-inoculation experiments demonstrated that F. moniliforme is dominant among these FM pathogenic fungi. These results lay the foundation for future studies on mitigating or preventing FM damage to soybean.

Structural Variation (SV) Markers in the Basidiomycete Volvariella volvacea and Their Application in the Construction of a Genetic Map

Molecular markers and genetic maps are useful tools in genetic studies. Novel molecular markers and their applications have been developed in recent years. With the recent advancements in sequencing technology, the genomic sequences of an increasingly great number of fungi have become available. A novel type of molecular marker was developed to construct the first reported linkage map of the edible and economically important basidiomycete Volvariella volvacea by using 104 structural variation (SV) markers that are based on the genomic sequences. Because of the special and simple life cycle in basidiomycete, SV markers can be effectively developed by genomic comparison and tested in single spore isolates (SSIs). This stable, convenient and rapidly developed marker may assist in the construction of genetic maps and facilitate genomic research for other species of fungi.

Genome-wide patterns of segregation and linkage disequilibrium: the construction of a linkage genetic map of the poplar rust fungus Melampsora larici-populina

The poplar rust fungus Melampsora larici-populina causes significant yield reduction and severe economic losses in commercial poplar plantations. After several decades of breeding for qualitative resistance and subsequent breakdown of the released resistance genes, breeders now focus on quantitative resistance, perceived to be more durable. But quantitative resistance also can be challenged by an increase of aggressiveness in the pathogen. Thus, it is of primary importance to better understand the genetic architecture of aggressiveness traits. To this aim, our goal is to build a genetic linkage map for M. larici-populina in order to map quantitative trait loci related to aggressiveness. First, a large progeny of M. larici-populina was generated through selfing of the reference strain 98AG31 (which genome sequence is available) on larch plants, the alternate host of the poplar rust fungus. The progeny's meiotic origin was validated through a segregation analysis of 115 offspring with 14 polymorphic microsatellite markers, of which 12 segregated in the expected 1:2:1 Mendelian ratio. A microsatellite-based linkage disequilibrium analysis allowed us to identify one potential linkage group comprising two scaffolds. The whole genome of a subset of 47 offspring was resequenced using the Illumina HiSeq 2000 technology at a mean sequencing depth of 6X. The reads were mapped onto the reference genome of the parental strain and 144,566 SNPs were identified across the genome. Analysis of distribution and polymorphism of the SNPs along the genome led to the identification of 2580 recombination blocks. A second linkage disequilibrium analysis, using the recombination blocks as markers, allowed us to group 81 scaffolds into 23 potential linkage groups. These preliminary results showed that a high-density linkage map could be constructed by using high-quality SNPs based on low-coverage resequencing of a larger number of M. larici-populina offspring.

The influence of phosphorus availability and Laccaria bicolor symbiosis on phosphate acquisition, antioxidant enzyme activity, and rhizospheric carbon flux in Populus tremuloides

Many forest tree species are dependent on their symbiotic interaction with ectomycorrhizal (ECM) fungi for phosphorus (P) uptake from forest soils where P availability is often limited. The ECM fungal association benefits the host plant under P limitation through enhanced soil exploration and increased P acquisition by mycorrhizas. To study the P starvation response (PSR) and its modification by ECM fungi in Populus tremuloides, a comparison was made between nonmycorrhizal (NM) and mycorrhizal with Laccaria bicolor (Myc) seedlings grown under different concentrations of phosphate (Pi) in sand culture. Although differences in growth between NM and Myc plants were small, Myc plants were more effective at acquiring P from low Pi treatments, with significantly lower k m values for root and leaf P accumulation. Pi limitation significantly increased the activity of catalase, ascorbate peroxidase, and guaiacol-dependent peroxidase in leaves and roots to greater extents in NM than Myc P. tremuloides. Phosphoenolpyruvate carboxylase activity also increased in NM plants under P limitation, but was unchanged in Myc plants. Formate, citrate, malonate, lactate, malate, and oxalate and total organic carbon exudation by roots was stimulated by P limitation to a greater extent in NM than Myc plants. Colonization by L. bicolor reduced the solution Pi concentration thresholds where PSR physiological changes occurred, indicating that enhanced Pi acquisition by P. tremuloides colonized by L. bicolor altered host P homeostasis and plant stress responses to P limitation. Understanding these plant-symbiont interactions facilitates the selection of more P-efficient forest trees and strategies for tree plantation production on marginal soils.

A 2nd generation linkage map of Heterobasidion annosum s.l. based on in silico anchoring of AFLP markers

In this study, we present a 2^nd 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 1^st 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.

Characterization of transposable elements in the ectomycorrhizal fungus Laccaria bicolor

Background

The publicly available Laccaria bicolor genome sequence has provided a considerable genomic resource allowing systematic identification of transposable elements (TEs) in this symbiotic ectomycorrhizal fungus. Using a TE-specific annotation pipeline we have characterized and analyzed TEs in the L. bicolor S238N-H82 genome.

Methodology/Principal Findings

TEs occupy 24% of the 60 Mb L. bicolor genome and represent 25,787 full-length and partial copy elements distributed within 171 families. The most abundant elements were the Copia-like. TEs are not randomly distributed across the genome, but are tightly nested or clustered. The majority of TEs exhibits signs of ancient transposition except some intact copies of terminal inverted repeats (TIRS), long terminal repeats (LTRs) and a large retrotransposon derivative (LARD) element. There were three main periods of TE expansion in L. bicolor: the first from 57 to 10 Mya, the second from 5 to 1 Mya and the most recent from 0.5 Mya ago until now. LTR retrotransposons are closely related to retrotransposons found in another basidiomycete, Coprinopsis cinerea.

Conclusions

This analysis 1) represents an initial characterization of TEs in the L. bicolor genome, 2) contributes to improve genome annotation and a greater understanding of the role TEs played in genome organization and evolution and 3) provides a valuable resource for future research on the genome evolution within the Laccaria genus.

Genetic linkage mapping in fungi: current state, applications, and future trends

Genetic mapping is a basic tool for eukaryotic genomic research. Linkage maps provide insights into genome organization and can be used for genetic studies of traits of interest. A genetic linkage map is a suitable support for the anchoring of whole genome sequences. It allows the localization of genes of interest or quantitative trait loci (QTL) and map-based cloning. While genetic mapping has been extensively used in plant or animal models, this discipline is more recent in fungi. The present article reviews the current status of genetic linkage map research in fungal species. The process of linkage mapping is detailed, from the development of mapping populations to the construction of the final linkage map, and illustrated based on practical examples. The range of specific applications in fungi is browsed, such as the mapping of virulence genes in pathogenic species or the mapping of agronomically relevant QTL in cultivated edible mushrooms. Future prospects are finally discussed in the context of the most recent advances in molecular techniques and the release of numerous fungal genome sequences.

Phylogenetic, genomic organization and expression analysis of hydrophobin genes in the ectomycorrhizal basidiomycete Laccaria bicolor

Hydrophobins are morphogenetic, small secreted hydrophobic fungal proteins produced in response to changing development and environmental conditions. These proteins are important in the interaction between certain fungi and their hosts. In mutualistic ectomycorrhizal fungi several hydrophobins form a subclass of mycorrhizal-induced small secreted proteins that are likely to be critical in the formation of the symbiotic interface with host root cells. In this study, two genomes of the ectomycorrhizal basidiomycete Laccaria bicolor strains S238N-H82 (from North America) and 81306 (from Europe) were surveyed to construct a comprehensive genome-wide inventory of hydrophobins and to explore their characteristics and roles during host colonization. The S238N-H82 L. bicolor hydrophobin gene family is composed of 12 genes while the 81306 strain encodes nine hydrophobins, all corresponding to class I hydrophobins. The three extra hydrophobin genes encoded by the S238N-H82 genome likely arose via gene duplication and are bordered by transposon rich regions. Expression profiles of the hydrophobin genes of L. bicolor varied greatly depending on life stage (e.g. free living mycelium vs. root colonization) and on the host root environment. We conclude from this study that the complex diversity and range of expression profiles of the Laccaria hydrophobin multi-gene family have likely been a selective advantage for this mutualist in colonizing a wide range of host plants.

Evaluation of the use of SCAR markers for screening genetic diversity of Lentinula edodes strains

In this study, three molecular marker systems including sequence related amplified polymorphism (SRAP), random amplified polymorphic DNA (RAPD), and inter-simple sequence repeats (ISSR) were screened to select polymorphisms of 24 main commercial strains of Lentinula edodes cultivated widely in China. Twenty-nine sequence characterized amplified region (SCAR) markers were developed to set up a dendrogram using UPMGA based on nucleotide sequences of some SRAP, RAPD, and ISSR polymorphic fragments. The grouping showed that the 24 strains were apparently clustered into five groups at a level of 0.68 similarity coefficient, and those that have similar breeding background clustered preferentially into the same subgroup. Results also revealed that the 24 strains had a low level of genetic diversity, and the breeding source of L. edodes should be broadened by exploiting wild types and introducing exotic strains. In addition, the tested strains of L. edodes could be clearly distinguished and identified from others by using different combinations of SCAR primers. Thus, results of this work demonstrated that SCAR was an excellent genetic marker system to characterize and investigate genetic diversity of L. edodes. Furthermore, this provided an alternative method to identify the genetic relationship of different strains of other fungi.

Construction of a genetic linkage map based on amplified fragment length polymorphism markers and development of sequence-tagged site markers for marker-assisted selection of the sporeless trait in the oyster mushroom (Pleurotus eryngii)

A large number of spores from fruiting bodies can lead to allergic reactions and other problems during the cultivation of edible mushrooms, including Pleurotus eryngii (DC.) Quél. A cultivar harboring a sporulation-deficient (sporeless) mutation would be useful for preventing these problems, but traditional breeding requires extensive time and labor. In this study, using a sporeless P. eryngii strain, we constructed a genetic linkage map to introduce a molecular breeding program like marker-assisted selection. Based on the segregation of 294 amplified fragment length polymorphism markers, two mating type factors, and the sporeless trait, the linkage map consisted of 11 linkage groups with a total length of 837.2 centimorgans (cM). The gene region responsible for the sporeless trait was located in linkage group IX with 32 amplified fragment length polymorphism markers and the B mating type factor. We also identified eight markers closely linked (within 1.2 cM) to the sporeless locus using bulked-segregant analysis-based amplified fragment length polymorphism. One such amplified fragment length polymorphism marker was converted into two sequence-tagged site markers, SD488-I and SD488-II. Using 14 wild isolates, sequence-tagged site analysis indicated the potential usefulness of the combination of two sequence-tagged site markers in cross-breeding of the sporeless strain. It also suggested that a map constructed for P. eryngii has adequate accuracy for marker-assisted selection.

Extensive gene flow over Europe and possible speciation over Eurasia in the ectomycorrhizal basidiomycete Laccaria amethystina complex

Biogeographical patterns and large-scale genetic structure have been little studied in ectomycorrhizal (EM) fungi, despite the ecological and economic importance of EM symbioses. We coupled population genetics and phylogenetic approaches to understand spatial structure in fungal populations on a continental scale. Using nine microsatellite markers, we characterized gene flow among 16 populations of the widespread EM basidiomycete Laccaria amethystina over Europe (i.e. over 2900 km). We also widened our scope to two additional populations from Japan (10(4) km away) and compared them with European populations through microsatellite markers and multilocus phylogenies, using three nuclear genes (NAR, G6PD and ribosomal DNA) and two mitochondrial ribosomal genes. European L. amethystina populations displayed limited differentiation (average F(ST) = 0.041) and very weak isolation by distance (IBD). This panmictic European pattern may result from effective aerial dispersal of spores, high genetic diversity in populations and mutualistic interactions with multiple hosts that all facilitate migration. The multilocus phylogeny based on nuclear genes confirmed that Japanese and European specimens were closely related but clustered on a geographical basis. By using microsatellite markers, we found that Japanese populations were strongly differentiated from the European populations (F(ST) = 0.416), more than expected by extrapolating the European pattern of IBD. Population structure analyses clearly separated the populations into two clusters, i.e. European and Japanese clusters. We discuss the possibility of IBD in a continuous population (considering some evidence for a ring species over the Northern Hemisphere) vs. an allopatric speciation over Eurasia, making L. amethystina a promising model of intercontinental species for future studies.

Population genetics of ectomycorrhizal fungi: from current knowledge to emerging directions

Ectomycorrhizal (EM) fungi are major microbial components of boreal, temperate and Mediterranean forests, as well as some tropical forest ecosystems. Nearly two decades of studies have clarified many aspects of their population biology, based on several model species from diverse lineages of fungi where the EM symbiosis evolved, i.e. among Hymenomycetes and, to a lesser extent, among Ascomycetes. In this review, we show how tools for individual recognition have changed, shifting from the use of somatic incompatibility reactions to dominant and non-specific markers (such as random amplified polymorphic DNA (RAPD) and amplified fragment length polymorphism (AFLP)) and, more recently, to co-dominant and specific markers (such as microsatellites and single nucleotide polymorphisms (SNPs)). At the same time, the theoretical focus has also changed. In earlier studies, a major aim was the description of genet size and popul/ation strategy. For example, we show how some studies supported or challenged the simple, classical model of colonization of new forest stands by ruderal (R) species, propagating by spores and forming small genets, progressively replaced in older forests by more competitive (C) species, propagating by mycelial growth and forming larger genets. By contrast, more recent studies give insights into some genetic traits, such as partners' assortment (allo- versus autogamy), genetic structure of populations and gene flow that turn out to depend both on distance and on whether spores are animal- or wind-dispersed. We discuss the rising awareness that (i) many morphospecies contain cryptic biological species (often sympatric) and (ii) trans- and inter-continental species may often contain several biological species isolated by distance. Finally, we show the emergence of biogeographic approaches and call for some aspects to be developed, such as fine-scale and long-term population monitoring, analyses of subterranean populations of extra-radical mycelia, or more model species from the tropics, as well as from the Ascomycetes (whose genetic idiosyncrasies are discussed). With the rise of the '-omics' sciences, analysis of population structure for non-neutral genes is expected to develop, and forest management and conservation biology will probably profit from published and expected work.

Molecular and histochemical characterisation of two distinct poplar Melampsora leaf rust pathosystems

In this study, we compared interactions of two Melampsora foliar rust species with poplar, which resulted in either limited or abundant pathogen proliferation. In the pathosystem exhibiting limited pathogen growth, a defence response was observed after invasion of poplar leaf tissues by the biotroph, with late and clear production of reactive oxygen species (ROS) and other products. Characterisation of the histological, biochemical and transcriptional events occurring in both pathosystems showed striking similarity with components of plant defence reactions observed during qualitative resistance. Key components associated with development of an active defence response, such as up-regulation of pathogenesis-related (PR) genes, were observed during infection. Moreover, the time course and strength of gene induction appear to be critical determinants for the outcome of the tree-pathogen interaction. This work provides basic biochemical characterisation and expression data for the study of so-called partial resistance in the poplar-rust pathosystem, which is also applicable to other plant-pathogen interactions resulting in quantitative disease resistance.

An expanded genetic linkage map of an intervarietal Agaricus bisporus var. bisporusxA. bisporus var. burnettii hybrid based on AFLP, SSR and CAPS markers sheds light on the recombination behaviour of the species

A genetic linkage map for the edible basidiomycete Agaricus bisporus was constructed from 118 haploid homokaryons derived from an intervarietal A. bisporus var. bisporus x A. bisporus var. burnettii hybrid. Two hundred and thirty-one AFLP, 21 SSR, 68 CAPS markers together with the MAT, BSN, PPC1 loci and one allozyme locus (ADH) were evenly spread over 13 linkage groups corresponding to the chromosomes of A. bisporus. The map covers 1156cM, with an average marker spacing of 3.9cM and encompasses nearly the whole genome. The average number of crossovers per chromosome per individual is 0.86. Normal recombination over the entire genome occurs in the heterothallic variety, burnettii, contrary to the homothallic variety, bisporus, which showed adaptive genome-wide suppressed recombination. This first comprehensive genetic linkage map for A. bisporus provides foundations for quantitative trait analyses and breeding programme monitoring, as well as genome organisation studies.

Phylogenetic analysis, genomic organization, and expression analysis of multi-copper oxidases in the ectomycorrhizal basidiomycete Laccaria bicolor

In forest soils, ectomycorrhizal and saprotrophic Agaricales differ in their strategies for carbon acquisition, but share common gene families encoding multi-copper oxidases (MCOs). These enzymes are involved in the oxidation of a variety of soil organic compounds. The MCO gene family of the ectomycorrhizal fungus Laccaria bicolor is composed of 11 genes divided into two distinct subfamilies corresponding to laccases (lcc) sensu stricto (lcc1 to lcc9), sharing a high sequence homology with the coprophilic Coprinopsis cinerea laccase genes, and to ferroxidases (lcc10 and lcc11) that are not present in C. cinerea. The fet3-like ferroxidase genes lcc10 and lcc11 in L. bicolor are each arranged in a mirrored tandem orientation with an ftr gene coding for an iron permease. Unlike C. cinerea, L. bicolor has no sid1/sidA gene for siderophore biosynthesis. Transcript profiling using whole-genome expression arrays and quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR) revealed that some transcripts were very abundant in ectomycorrhizas (lcc3 and lcc8), in fruiting bodies (lcc7) or in the free-living mycelium grown on agar medium (lcc9 and lcc10), suggesting a specific function of these MCOs. The amino acid composition of the MCO substrate binding sites suggests that L. bicolor MCOs interact with substrates different from those of saprotrophic fungi.

The Laccaria genome: a symbiont blueprint decoded

The first genomic sequence for a representative of symbiotic fungi, the ectomycorrhizal basidiomycete Laccaria bicolor, has been published. The unravelling of this genome provides tantalizing hints about differences between this symbiotic fungus and its saprotrophic and pathogenic relatives. An expansion of several multigene families occurred in L. bicolor, suggesting that adaptation to symbiosis proceeded by gene duplication. Within lineage-specific genes those coding for symbiosis-regulated secreted proteins showed an up-regulated expression in ectomycorrhizas. L. bicolor is lacking enzymes involved in the degradation of plant cell wall components (cellulose, hemicellulose, pectins and pectates), preventing the symbiont from degrading host cells. By contrast, L. bicolor possesses expanded multigene families associated with hydrolysis of bacterial and microfauna polysaccharides and proteins. The genome analysis revealed the dual saprotrophic and biotrophic lifestyle of the mycorrhizal fungus that enables it to grow within both soil and living plant roots. The next stages will involve finer-scale investigation of gene networks to reveal the details of the general patterns now uncovered at the genomic level. The acceptance of L. bicolor as a model organism for symbiosis genetics will, however, depend strongly on the availability of additional genetic, genomic and molecular biological resources, such as gene inactivation procedures.