Comparative linkage mapping in the white button mushroom Agaricus bisporus provides foundation for breeding management

Mating-Type Locus Organization and Mating-Type Chromosome Differentiation in the Bipolar Edible Button Mushroom Agaricus bisporus

In heterothallic basidiomycete fungi, sexual compatibility is restricted by mating types, typically controlled by two loci: PR, encoding pheromone precursors and pheromone receptors, and HD, encoding two types of homeodomain transcription factors. We analysed the single mating-type locus of the commercial button mushroom variety, Agaricus bisporus var. bisporus, and of the related variety burnettii. We identified the location of the mating-type locus using genetic map and genome information, corresponding to the HD locus, the PR locus having lost its mating-type role. We found the mip1 and β-fg genes flanking the HD genes as in several Agaricomycetes, two copies of the β-fg gene, an additional HD2 copy in the reference genome of A. bisporus var. bisporus and an additional HD1 copy in the reference genome of A. bisporus var. burnettii. We detected a 140 kb-long inversion between mating types in an A. bisporus var. burnettii heterokaryon, trapping the HD genes, the mip1 gene and fragments of additional genes. The two varieties had islands of transposable elements at the mating-type locus, spanning 35 kb in the A. bisporus var. burnettii reference genome. Linkage analyses showed a region with low recombination in the mating-type locus region in the A. bisporus var. burnettii variety. We found high differentiation between β-fg alleles in both varieties, indicating an ancient event of recombination suppression, followed more recently by a suppression of recombination at the mip1 gene through the inversion in A. bisporus var. burnettii and a suppression of recombination across whole chromosomes in A. bisporus var. bisporus, constituting stepwise recombination suppression as in many other mating-type chromosomes and sex chromosomes.

Induced mutation in Agaricus bisporus by gamma ray to improve genetic variability, degradation enzyme activity, and yield

PURPOSE

Genetic variability in white button mushroom cultivars is very low due to the life cycle. Induction mutations using gamma irradiation is a useful way to generate diversity in white button mushrooms to obtain genotype(s) with desirable traits.

METHODS

Gamma irradiation Cobalt-60 was used for inducting genetic diversity in white button mushroom to obtain genotype(s) with desirable traits. Gamma irradiation with doses of 0-500 Gy was conducted on spores on Potato Dextrose Agar medium.

RESULTS

The results showed significant differences in days to pin production and harvest, fruit body number, fresh and dry weight, yield, laccase, and manganese peroxidase enzyme activity. After isolating variants, 15 variants were selected on the base of their high yield and enzyme degradation activity. Their genetic variation was confirmed by Sequence Related Amplified Polymorphism (SRAP) markers, and then incubated on three types of substrates (50:50, 75:25, and 100:0 % compost: raw straw). The results showed that all variants, except GR18, colonized in 75:25, and GR3, GR4, GR9, GR61, GR72, and GR74 variants colonized in 50:50. In 100:0 substrate, GR55 and GR63 were the earliest variants, and GR9 produced the highest fruit body number. In 75:25 substrate, GR9, GR3, GR61, GR4, GR74, GR4, GR61, and GR72 showed higher yields. The highest laccase and manganese peroxidase activity were recorded in GR3, GR4, GR9, GR72, and GR61. The isolated 15 variants were clustered into two main groups by cluster analysis and genetic variation was confirmed by SRAP markers.

CONCLUSION

The results showed that the diversity in the white button mushroom could be improved using gamma rays, and the variation would be useful for the development of future breeding programs.

Recombination suppression and evolutionary strata around mating-type loci in fungi: documenting patterns and understanding evolutionary and mechanistic causes

Genomic regions determining sexual compatibility often display recombination suppression, as occurs in sex chromosomes, plant self-incompatibility loci and fungal mating-type loci. Regions lacking recombination can extend beyond the genes determining sexes or mating types, by several successive steps of recombination suppression. Here we review the evidence for recombination suppression around mating-type loci in fungi, sometimes encompassing vast regions of the mating-type chromosomes. The suppression of recombination at mating-type loci in fungi has long been recognized and maintains the multiallelic combinations required for correct compatibility determination. We review more recent evidence for expansions of recombination suppression beyond mating-type genes in fungi ('evolutionary strata'), which have been little studied and may be more pervasive than commonly thought. We discuss testable hypotheses for the ultimate (evolutionary) and proximate (mechanistic) causes for such expansions of recombination suppression, including (1) antagonistic selection, (2) association of additional functions to mating-type, such as uniparental mitochondria inheritance, (3) accumulation in the margin of nonrecombining regions of various factors, including deleterious mutations or transposable elements resulting from relaxed selection, or neutral rearrangements resulting from genetic drift. The study of recombination suppression in fungi could thus contribute to our understanding of recombination suppression expansion across a broader range of organisms.

Variation of the meiotic recombination landscape and properties over a broad evolutionary distance in yeasts

Meiotic recombination is a major factor of genome evolution, deeply characterized in only a few model species, notably the yeast Saccharomyces cerevisiae. Consequently, little is known about variations of its properties across species. In this respect, we explored the recombination landscape of Lachancea kluyveri, a protoploid yeast species that diverged from the Saccharomyces genus more than 100 million years ago and we found striking differences with S. cerevisiae. These variations include a lower recombination rate, a higher frequency of chromosomes segregating without any crossover and the absence of recombination on the chromosome arm containing the sex locus. In addition, although well conserved within the Saccharomyces clade, the S. cerevisiae recombination hotspots are not conserved over a broader evolutionary distance. Finally and strikingly, we found evidence of frequent reversal of commitment to meiosis, resulting in return to mitotic growth after allele shuffling. Identification of this major but underestimated evolutionary phenomenon illustrates the relevance of exploring non-model species.

Meiotic recombination promotes accurate chromosome segregation and genetic diversity. To date, the mechanisms and rules lying behind recombination were dissected using model organisms such as the budding yeast Saccharomyces cerevisiae. To assess the conservation and variation of this process over a broad evolutionary distance, we explored the meiotic recombination landscape in Lachancea kluyveri, a budding yeast species that diverged from S. cerevisiae more than 100 million years ago. The meiotic recombination map we generated revealed that the meiotic recombination landscape and properties significantly vary across distantly related yeast species, raising the yet to confirm possibility that recombination hotspots conservation across yeast species depends on synteny conservation. Finally, the frequent meiotic reversions we observed led us to re-evaluate their evolutionary importance.

Developments in breeding of Agaricus bisporus var. bisporus: progress made and technical and legal hurdles to take

True breeding of button mushrooms has hardly been done in the last decades, despite this species being one of the most cultivated mushrooms worldwide. Research done in the last 20 years has identified and characterised new germplasm and improved our understanding of the genetic base for some traits. A substantial collection of wild-collected strains is now available and partly characterised for a number of important traits such as disease resistance and yield. Most of the variations found in a number of important agronomic traits have a considerable heritability and are thus useful for breeding. Genetic marker technology has also developed considerably for this mushrooms in the last decade and used to identify quantitative trait loci (QTL) for important agronomic traits. This progress has, except for one example, not resulted so far into new commercially varieties. One of the reasons lies in the typical life cycle of the button mushroom Agaricus bisporus var. bisporus which hampers breeding. Joint investment is needed to solve technical problems in breeding. Special attention is needed for the protection of new varieties. Due to its typical life cycle, it is very easy to generate so called "look-a-likes" from protected cultivars by screening fertile single spore cultures. A consensus has been reached within the mushroom (breeding) industry to consider this method as the generation of essentially derived varieties as defined in plant breeding.

A detailed analysis of the recombination landscape of the button mushroom Agaricus bisporus var. bisporus

The button mushroom (Agaricus bisporus) is one of the world's most cultivated mushroom species, but in spite of its economic importance generation of new cultivars by outbreeding is exceptional. Previous genetic analyses of the white bisporus variety, including all cultivars and most wild isolates revealed that crossing over frequencies are low, which might explain the lack of introducing novel traits into existing cultivars. By generating two high quality whole genome sequence assemblies (one de novo and the other by improving the existing reference genome) of the first commercial white hybrid Horst U1, a detailed study of the crossover (CO) landscape was initiated. Using a set of 626 SNPs in a haploid offspring of 139 single spore isolates and whole genome sequencing on a limited number of homo- and heterokaryotic single spore isolates, we precisely mapped all COs showing that they are almost exclusively restricted to regions of about 100kb at the chromosome ends. Most basidia of A. bisporus var. bisporus produce two spores and pair preferentially via non-sister nuclei. Combined with the COs restricted to the chromosome ends, these spores retain most of the heterozygosity of the parent thus explaining how present-day white cultivars are genetically so close to the first hybrid marketed in 1980. To our knowledge this is the first example of an organism which displays such specific CO landscape.

The Genetic Linkage Map of the Medicinal Mushroom Agaricus subrufescens Reveals Highly Conserved Macrosynteny with the Congeneric Species Agaricus bisporus

Comparative linkage mapping can rapidly facilitate the transfer of genetic information from model species to orphan species. This macrosynteny analysis approach has been extensively used in plant species, but few example are available in fungi, and even fewer in mushroom crop species. Among the latter, the Agaricus genus comprises the most cultivable or potentially cultivable species. Agaricus bisporus, the button mushroom, is the model for edible and cultivable mushrooms. We have developed the first genetic linkage map for the basidiomycete A. subrufescens, an emerging mushroom crop known for its therapeutic properties and potential medicinal applications. The map includes 202 markers distributed over 16 linkage groups (LG), and covers a total length of 1701 cM, with an average marker spacing of 8.2 cM. Using 96 homologous loci, we also demonstrated the high level of macrosynteny with the genome of A. bisporus The 13 main LG of A. subrufescens were syntenic to the 13 A. bisporus chromosomes. A disrupted synteny was observed for the three remaining A. subrufescens LG. Electronic mapping of a collection of A. subrufescens expressed sequence tags on A. bisporus genome showed that the homologous loci were evenly spread, with the exception of a few local hot or cold spots of homology. Our results were discussed in the light of Agaricus species evolution process. The map provides a framework for future genetic or genomic studies of the medicinal mushroom A. subrufescens.

Constructing a new integrated genetic linkage map and mapping quantitative trait loci for vegetative mycelium growth rate in Lentinula edodes

The most saturated linkage map for Lentinula edodes to date was constructed based on a monokaryotic population of 146 single spore isolates (SSIs) using sequence-related amplified polymorphism (SRAP), target region amplification polymorphism (TRAP), insertion-deletion (InDel) markers, and the mating-type loci. Five hundred and twenty-four markers were located on 13 linkage groups (LGs). The map spanned a total length of 1006.1 cM, with an average marker spacing of 2.0 cM. Quantitative trait loci (QTLs) mapping was utilized to uncover the loci regulating and controlling the vegetative mycelium growth rate on various synthetic media, and complex medium for commercial cultivation of L. edodes. Two and 13 putative QTLs, identified respectively in the monokaryotic population and two testcross dikaryotic populations, were mapped on seven different LGs. Several vegetative mycelium growth rate-related QTLs uncovered here were clustered on LG4 (Qmgr1, Qdgr1, Qdgr2 and Qdgr9) and LG6 (Qdgr3, Qdgr4 and Qdgr5), implying the presence of main genomic areas responsible for growth rate regulation and control. The QTL hotspot region on LG4 was found to be in close proximity to the region containing the mating-type A (MAT-A) locus. Moreover, Qdgr2 on LG4 was detected on different media, contributing 8.07 %-23.71 % of the phenotypic variation. The present study provides essential information for QTL mapping and marker-assisted selection (MAS) in L. edodes.

Transcriptomic analysis of the interactions between Agaricus bisporus and Lecanicillium fungicola

Agaricus bisporus is susceptible to a number of diseases, particularly those caused by fungi, with Lecanicillium fungicola being the most serious. Control of this disease is important for the security of crop production, however given the lack of knowledge about fungal-fungal interactions, such disease control is rather limited. Exploiting the recently released genome sequence of A. bisporus, here we report studies simultaneously investigating both the host and the pathogen, focussing on transcriptional changes associated with the cap spotting lesions typically seen in this interaction. Forward-suppressive subtractive hybridisation (SSH) analysis identified 68 A. bisporus unigenes induced during infection. Chitin deacetylase showed the strongest response, with almost 1000-fold up-regulation during infection, so was targeted for down-regulation by silencing to see if it was involved in defence against L. fungicola. Transgenic lines were made expressing hairpin RNAi constructs, however no changes in susceptibility to L. fungicola were observed. Amongst the other up-regulated genes there were none with readily apparent roles in resisting infection in this susceptible interaction. Reverse-SSH identified 72 unigenes from A. bisporus showing reduced expression, including two tyrosinases, several genes involved in nitrogen metabolism and a hydrophobin. The forward-SSH analysis of infected mushrooms also yielded 64 unigenes which were not of A. bisporus origin and thus derived from L. fungicola. An EST analysis of infection-mimicking conditions generated an additional 623 unigenes from L. fungicola including several oxidoreductases, cell wall degrading enzymes, ABC and MFS transporter proteins and various other genes believed to play roles in other pathosystems. Together, this analysis shows how both the pathogen and the host modify their gene expression during an infection-interaction, shedding some light on the disease process, although we note that some 40% of unigenes from both organisms encode hypothetical proteins with no ascribed function which highlights how much there is still to discover about this interaction.

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.

Genome sequence of the button mushroom Agaricus bisporus reveals mechanisms governing adaptation to a humic-rich ecological niche

Agaricus bisporus is the model fungus for the adaptation, persistence, and growth in the humic-rich leaf-litter environment. Aside from its ecological role, A. bisporus has been an important component of the human diet for over 200 y and worldwide cultivation of the "button mushroom" forms a multibillion dollar industry. We present two A. bisporus genomes, their gene repertoires and transcript profiles on compost and during mushroom formation. The genomes encode a full repertoire of polysaccharide-degrading enzymes similar to that of wood-decayers. Comparative transcriptomics of mycelium grown on defined medium, casing-soil, and compost revealed genes encoding enzymes involved in xylan, cellulose, pectin, and protein degradation are more highly expressed in compost. The striking expansion of heme-thiolate peroxidases and β-etherases is distinctive from Agaricomycotina wood-decayers and suggests a broad attack on decaying lignin and related metabolites found in humic acid-rich environment. Similarly, up-regulation of these genes together with a lignolytic manganese peroxidase, multiple copper radical oxidases, and cytochrome P450s is consistent with challenges posed by complex humic-rich substrates. The gene repertoire and expression of hydrolytic enzymes in A. bisporus is substantially different from the taxonomically related ectomycorrhizal symbiont Laccaria bicolor. A common promoter motif was also identified in genes very highly expressed in humic-rich substrates. These observations reveal genetic and enzymatic mechanisms governing adaptation to the humic-rich ecological niche formed during plant degradation, further defining the critical role such fungi contribute to soil structure and carbon sequestration in terrestrial ecosystems. Genome sequence will expedite mushroom breeding for improved agronomic characteristics.

Polyploidy in fungi: evolution after whole-genome duplication

Polyploidy is a major evolutionary process in eukaryotes-particularly in plants and, to a less extent, in animals, wherein several past and recent whole-genome duplication events have been described. Surprisingly, the incidence of polyploidy in other eukaryote kingdoms, particularly within fungi, remained largely disregarded by the scientific community working on the evolutionary consequences of polyploidy. Recent studies have significantly increased our knowledge of the occurrence and evolutionary significance of fungal polyploidy. The ecological, structural and functional consequences of polyploidy in fungi are reviewed here and compared with the knowledge acquired with conventional plant and animal models. In particular, the genus Saccharomyces emerges as a relevant model for polyploid studies, in addition to plant and animal models.

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.

Relationship between yield components and partial resistance to Lecanicillium fungicola in the button mushroom, Agaricus bisporus, assessed by quantitative trait locus mapping

Dry bubble, caused by Lecanicillium fungicola, is one of the most detrimental diseases affecting button mushroom cultivation. In a previous study, we demonstrated that breeding for resistance to this pathogen is quite challenging due to its quantitative inheritance. A second-generation hybrid progeny derived from an intervarietal cross between a wild strain and a commercial cultivar was characterized for L. fungicola resistance under artificial inoculation in three independent experiments. Analysis of quantitative trait loci (QTL) was used to determine the locations, numbers, and effects of genomic regions associated with dry-bubble resistance. Four traits related to resistance were analyzed. Two to four QTL were detected per trait, depending on the experiment. Two genomic regions, on linkage group X (LGX) and LGVIII, were consistently detected in the three experiments. The genomic region on LGX was detected for three of the four variables studied. The total phenotypic variance accounted for by all QTL ranged from 19.3% to 42.1% over all traits in all experiments. For most of the QTL, the favorable allele for resistance came from the wild parent, but for some QTL, the allele that contributed to a higher level of resistance was carried by the cultivar. Comparative mapping with QTL for yield-related traits revealed five colocations between resistance and yield component loci, suggesting that the resistance results from both genetic factors and fitness expression. The consequences for mushroom breeding programs are discussed.