Genetic linkage map of the edible basidiomycete Pleurotus ostreatus

Quantitative trait loci analysis for molecular markers linked to agricultural traits of Pleurotus ostreatus

Pleurotus ostreatus is a global mushroom crop with nutritional and medicinal benefits. However, the genetic basis of several commercial traits remains unknown. To address this, we analyzed the quantitative trait loci (QTLs) for two representative cultivars, "Heuktari" and "Miso," with apparently distinct alleles. A genetic map with 11 linkage groups was constructed, in which 27 QTLs were assigned to 14 traits. The explained phenotypic variations in QTLs ranged from 7.8% to 22.0%. Relatively high LOD values of 6.190 and 5.485 were estimated for the pinheading period and the number of valid stipes, respectively. Some QTL-derived molecular markers showed potential enhancement rates of selection precision in inbred lines, especially for cap shape (50%) and cap thickness (30%). Candidate genes were inferred from the QTL regions and validated using qRT-PCR, particularly for the cysteine and glutathione pathway, in relation to cap yellowness. The molecular markers in this study are expected to facilitate the breeding of the Heuktari and Miso lines and provide probes to identify related genes in P. ostreatus.

Pleurotus ostreatus as a model mushroom in genetics, cell biology, and material sciences

Pleurotus ostreatus, also known as the oyster mushroom, is a popular edible mushroom cultivated worldwide. This review aims to survey recent progress in the molecular genetics of this fungus and demonstrate its potential as a model mushroom for future research. The development of modern molecular genetic techniques and genome sequencing technologies has resulted in breakthroughs in mushroom science. With efficient transformation protocols and multiple selection markers, a powerful toolbox, including techniques such as gene knockout and genome editing, has been developed, and numerous new findings are accumulating in P. ostreatus. These include molecular mechanisms of wood component degradation, sexual development, protein secretion systems, and cell wall structure. Furthermore, these techniques enable the identification of new horizons in enzymology, biochemistry, cell biology, and material science through protein engineering, fluorescence microscopy, and molecular breeding. KEY POINTS: • Various genetic techniques are available in Pleurotus ostreatus. • P. ostreatus can be used as an alternative model mushroom in genetic analyses. • New frontiers in mushroom science are being developed using the fungus.

A genetic linkage map and improved genome assembly of the termite symbiont Termitomyces cryptogamus

BACKGROUND

The termite-fungus symbiosis is an ancient stable mutualism of two partners that reproduce and disperse independently. With the founding of each termite colony the symbiotic association must be re-established with a new fungus partner. Complementarity in the ability to break down plant substrate may help to stabilize this symbiosis despite horizontal symbiont transmission. An alternative, non-exclusive, hypothesis is that a reduced rate of evolution may contribute to stabilize the symbiosis, the so-called Red King Effect.

METHODS

To explore this concept, we produced the first linkage map of a species of Termitomyces, using genotyping by sequencing (GBS) of 88 homokaryotic offspring. We constructed a highly contiguous genome assembly using PacBio data and a de-novo evidence-based annotation. This improved genome assembly and linkage map allowed for examination of the recombination landscape and its potential effect on the mutualistic lifestyle.

RESULTS

Our linkage map resulted in a genome-wide recombination rate of 22 cM/Mb, lower than that of other related fungi. However, the total map length of 1370 cM was similar to that of other related fungi.

CONCLUSIONS

The apparently decreased rate of recombination is primarily due to genome expansion of islands of gene-poor repetitive sequences. This study highlights the importance of inclusion of genomic context in cross-species comparisons of recombination rate.

Analysis of Biochemical and Genetic Variability of Pleurotus ostreatus Based on the β-Glucans and CDDP Markers

Oyster mushroom (Pleurotus ostreatus) is still one of the most cultivated edible and medicinal mushrooms. Despite its frequent cultivation around the world, there is currently just a little information available on the variability of strains in terms of the content of β-glucans in them. This work presents an extensive study of 60 strains in terms of the content of α-glucans and β-glucans in their caps and stipes. The authenticity of the production strains based on an analysis of the variability of their genome by CDDP (Conserved DNA-derived polymorphism) markers was confirmed, whereas identical CDDP profiles were identified between samples 45, 89, 95, and 96. Genetic variability of the analyzed production strains showed a high polymorphism and effective discriminative power of the used marking technique. Medium positive correlations were found among the CDDP profiles and β-glucan content in the group of strains that generated the same CDDP profiles, and low negative correlation was found among these profiles in the group of low β-glucan content strains. For the determination of glucans content, Mushroom and Yeast analytical enzymatic kit (Megazyme, Bray, Co. Wicklow, Ireland) were used. The results clearly showed that the stipe contains on average 33% more β-glucans than the cap. The minimum detected β-glucan content in the stipe was in strain no. 72, specifically 22%, and the maximum in strain no. 43, specifically 56%, which after the conversion represents a difference of 155%. From the point of view of β-glucan content, the stated strain no. 43 appears to be very suitable for the commercial production of β-glucans under certain conditions.

The spore killers, fungal meiotic driver elements

During meiosis, both alleles of any given gene should have equal chances of being inherited by the progeny. There are a number of reasons why, however, this is not the case, with one of the most intriguing instances presenting itself as the phenomenon of meiotic drive. Genes that are capable of driving can manipulate the ratio of alleles among viable meiotic products so that they are inherited in more than half of them. In many cases, this effect is achieved by direct antagonistic interactions, where the driving allele inhibits or otherwise eliminates the alternative allele. In ascomycete fungi, meiotic products are packaged directly into ascospores; thus, the effect of meiotic drive has been given the nefarious moniker, "spore killing." In recent years, many of the known spore killers have been elevated from mysterious phenotypes to well-described systems at genetic, genomic, and molecular levels. In this review, we describe the known diversity of spore killers and synthesize the varied pieces of data from each system into broader trends regarding genome architecture, mechanisms of resistance, the role of transposable elements, their effect on population dynamics, speciation and gene flow, and finally how they may be developed as synthetic drivers. We propose that spore killing is common, but that it is under-observed because of a lack of studies on natural populations. We encourage researchers to seek new spore killers to build on the knowledge that these remarkable genetic elements can teach us about meiotic drive, genomic conflict, and evolution more broadly.

SNP-Based Genetic Linkage Map and Quantitative Trait Locus Mapping Associated with the Agronomically Important Traits of Hypsizygus marmoreus

White strains of Hypsizygus marmoreus are more difficult to cultivate than are brown strains; therefore, new white strain breeding strategies are required. Accordingly, we constructed the genetic map of H. marmoreus with 1996 SNP markers on 11 linkage groups (LGs) spanning 1380.49 cM. Prior to analysis, 82 backcrossed strains (HM8 lines) were generated by mating between KMCC03106-31 and the progenies of the F1 hybrid (Hami-18 × KMCC03106-93). Using HM8, the first 23 quantitative trait loci (QTLs) of yield-related traits were detected with high limit of detection (LOD) scores (1.98-9.86). The length, thickness, and hardness of the stipe were colocated on LG 1. Especially, length of stipe and thickness of stipe were highly correlated given that the correlation coefficients were negative (-0.39, p value ≤ .01). And a typical biomodal distribution was observed for lightness of the pileus and the lightness of the pileus trait belonged to the LG 8, as did traits of earliness and mycelial growth in potato dextrose agar (PDA) medium. Therefore, results for color traits can be suggested that color is controlled by a multi-gene of one locus. The yield trait was highly negatively correlated with the traits for thickness of the stipe (-0.45, p value ≤ .01). Based on additive effects, the white strain was confirmed as recessive; however, traits of mycelial growth, lightness, and quality were inherited by backcrossed HM8 lines. This new genetic map, finely mapped QTLs, and the strong selection markers could be used in molecular breeding of H. marmoreus.

Strain Degeneration in Pleurotus ostreatus: A Genotype Dependent Oxidative Stress Process Which Triggers Oxidative Stress, Cellular Detoxifying and Cell Wall Reshaping Genes

Strain degeneration has been defined as a decrease or loss in the yield of important commercial traits resulting from subsequent culture, which ultimately leads to Reactive Oxygen Species (ROS) production. Pleurotus ostreatus is a lignin-producing nematophagous edible mushroom. Mycelia for mushroom production are usually maintained in subsequent culture in solid media and frequently show symptoms of strain degeneration. The dikaryotic strain P. ostreatus (DkN001) has been used in our lab as a model organism for different purposes. Hence, different tools have been developed to uncover genetic and molecular aspects of this fungus. In this work, strain degeneration was studied in a full-sib monokaryotic progeny of the DkN001 strain with fast (F) and slow (S) growth rates by using different experimental approaches (light microscopy, malondialdehyde levels, whole-genome transcriptome analysis, and chitosan effect on monokaryotic mycelia). The results obtained showed that: (i) strain degeneration in P. ostreatus is linked to oxidative stress, (ii) the oxidative stress response in monokaryons is genotype dependent, (iii) stress and detoxifying genes are highly expressed in S monokaryons with symptoms of strain degeneration, (iv) chitosan addition to F and S monokaryons uncovered the constitutive expression of both oxidative stress and cellular detoxifying genes in S monokaryon strains which suggest their adaptation to oxidative stress, and (v) the overexpression of the cell wall genes, Uap1 and Cda1, in S monokaryons with strain degeneration phenotype indicates cell wall reshaping and the activation of High Osmolarity Glycerol (HOG) and Cell Wall Integrity (CWI) pathways. These results could constitute a hallmark for mushroom producers to distinguish strain degeneration in commercial mushrooms.

Genetic Linkage and Physical Mapping for an Oyster Mushroom (Pleurotus cornucopiae) and Quantitative Trait Locus Analysis for Cap Color

Oyster mushrooms are grown commercially worldwide, especially in many developing countries, for their easy cultivation and high biological efficiency. Pleurotus cornucopiae is one of the main oyster mushroom species because of its gastronomic value and nutraceutical properties. Cap color is an important trait, since consumers prefer dark mushrooms, which are now represented by only a small portion of the commercial varieties. Breeding efforts are required to improve quality-related traits to satisfy various demands of consumers. Here, we present a saturated genetic linkage map of P. cornucopiae constructed by using a segregating population of 122 monokaryons and 3,449 single nucleotide polymorphism (SNP) markers generated by the 2b-RAD approach. The map contains 11 linkage groups covering 961.6 centimorgans (cM), with an average marker spacing of 0.27 cM. The genome of P. cornucopiae was de novo sequenced, resulting in 425 scaffolds (>1,000 bp) with a total genome size of 35.1 Mb. The scaffolds were assembled to the pseudochromosome level with the assistance of the genetic linkage map. A total of 97% SNP markers (3,357) were physically localized on 140 scaffolds that were assigned to 11 pseudochromosomes, with a total of 32.5 Mb, representing 92.5% of the whole genome. Six quantitative trait loci (QTL) controlling cap color of P. cornucopiae were detected, accounting for a total phenotypic variation of 65.6%, with the highest value for the QTL on pseudochromosome 5 (18%). The results of our study provide a solid base for marker-assisted breeding for agronomic traits and especially for studies on biological mechanisms controlling cap color in oyster mushrooms. IMPORTANCE Oyster mushrooms are produced and consumed all over the world. Pleurotus cornucopiae is one of the main oyster mushroom species. Dark-cap oyster mushrooms are becoming more and more popular with consumers, but dark varieties are rare on the market. Prerequisites for efficient breeding programs are the availability of high-quality whole genomes and genetic linkage maps. Genetic studies to fulfill some of these prerequisites have hardly been done for P. cornucopiae. In this study, we de novo sequenced the genome and constructed a saturated genetic linkage map for P. cornucopiae. The genetic linkage map was effectively used to assist the genome assembly and identify QTL that genetically control the trait cap color. As well, the genome characteristics of P. cornucopiae were compared to the closely related species Pleurotus ostreatus. The results provided a basis for understanding the genetic background and marker-assisted breeding of this economically important mushroom species.

Genomic analyses of two Italian oyster mushroom Pleurotus pulmonarius strains

Pleurotus mushrooms are among the most cultivated fungi in the world and are highly valuable for food, medicine, and biotechnology industries. Furthermore, Pleurotus species are carnivorous fungi; they can rapidly paralyze and kill nematodes when nutrient-deprived. The predator-prey interactions between Pleurotus and nematodes are still widely unexplored. Moreover, the molecular mechanisms and the genes involved in the carnivorous behavior of Pleurotus mushrooms remain a mystery. We are attempting to understand the interactions between Pleurotus mushrooms and their nematode prey through genetic and genomic analyses. Two single spores (ss2 and ss5) isolated from a fruiting body of Pleurotus pulmonarius exhibited significant differences in growth and toxicity against nematodes. Thus, using PacBio long reads, we assembled and annotated two high-quality genomes for these two isolates of P. pulmonarius. Each of these assemblies contains 23 scaffolds, including 6 (ss2) and 8 (ss5) telomere-to-telomere scaffolds, and they are among the most complete assembled genomes of the Pleurotus species. Comparative analyses identified the genomic differences between the two P. pulmonarius strains. In sum, this work provides a genomic resource that will be invaluable for better understanding the Italian oyster mushroom P. pulmonarius.

Construction of a genetic linkage map and detection of quantitative trait locus for the ergothioneine content in tamogitake mushroom (Pleurotus cornucopiae var. citrinopileatus)

Developing high-content strains of L-ergothioneine (EGT), an antioxidant amino acid, is an important breeding target for tamogitake mushroom, Pleurotus cornucopiae var. citrinopileatus. We constructed a genetic linkage map based on segregation analysis of markers in 105 F1 progenies. The loci of 245 markers, including 10 AFLP markers, 195 Rad markers, 2 mating type factors, and 38 gene markers, were mapped. The map contained 12 linkage groups with a total genetic distance of 906.8 cM, and an average marker interval of 4.0 cM. The population from crossing between tester monokaryon and F1 progenies was used to characterize quantitative trait loci (QTL) for EGT content. With composite interval mapping (CIM) method, QTL of EGT content were found to be located in linkage group 10, having a Logarithm of the odds (LOD) score of 2.53 with a 10.1% contribution rate. Moreover, a single nucleotide polymorphism (SNP), A/T, was identified in a gene region of the genome in the neighborhood where the QTL peak existed. This SNP genotype was in good agreement with the EGT phenotypes of each strain in the both QTL population and wild population. Thus, this SNP would have great potential value to use the marker-assisted selection (MAS) for this mushroom with high EGT content.

A Resequencing-Based Ultradense Genetic Map of Hericium erinaceus for Anchoring Genome Sequences and Identifying Genetic Loci Associated With Monokaryon Growth

Hericium erinaceus has attracted tremendous interest owing to its compelling health-promoting properties. However, breeding of elite cultivars of H. erinaceus is hindered by the lack of a genetic and molecular toolbox. Here, we performed resequencing analysis of 127 F1 single-spore isolates and constructed the first high-resolution genetic map of H. erinaceus. With the use of recombination bins as markers, an ultradense genetic map consisting of 1,174 bins (including 37,082 single-nucleotide polymorphisms) was generated. This newly developed genetic map covered 1,096.5 cM, with an average bin spacing of 0.95 cM. High collinearity between genetic map and H. erinaceus genome assembly was revealed by aligning scaffolds to this genetic map using bin markers as anchors. The application of this newly developed genetic map in quantitative trait locus (QTL) mapping was also elucidated, and four QTLs for monokaryon growth were recovered. One QTL, mgr1, which contributes 12.1% of growth variations, was located near the mating type A (MAT-A) loci. Overall, this newly constructed high-resolution genetic map (or bin map) could be used as reference in future genetic, genomic, and breeding studies on H. erinaceus.

The genome sequence of the commercially cultivated mushroom Agrocybe aegerita reveals a conserved repertoire of fruiting-related genes and a versatile suite of biopolymer-degrading enzymes

BACKGROUND

Agrocybe aegerita is an agaricomycete fungus with typical mushroom features, which is commercially cultivated for its culinary use. In nature, it is a saprotrophic or facultative pathogenic fungus causing a white-rot of hardwood in forests of warm and mild climate. The ease of cultivation and fructification on solidified media as well as its archetypal mushroom fruit body morphology render A. aegerita a well-suited model for investigating mushroom developmental biology.

RESULTS

Here, the genome of the species is reported and analysed with respect to carbohydrate active genes and genes known to play a role during fruit body formation. In terms of fruit body development, our analyses revealed a conserved repertoire of fruiting-related genes, which corresponds well to the archetypal fruit body morphology of this mushroom. For some genes involved in fruit body formation, paralogisation was observed, but not all fruit body maturation-associated genes known from other agaricomycetes seem to be conserved in the genome sequence of A. aegerita. In terms of lytic enzymes, our analyses suggest a versatile arsenal of biopolymer-degrading enzymes that likely account for the flexible life style of this species. Regarding the amount of genes encoding CAZymes relevant for lignin degradation, A. aegerita shows more similarity to white-rot fungi than to litter decomposers, including 18 genes coding for unspecific peroxygenases and three dye-decolourising peroxidase genes expanding its lignocellulolytic machinery.

CONCLUSIONS

The genome resource will be useful for developing strategies towards genetic manipulation of A. aegerita, which will subsequently allow functional genetics approaches to elucidate fundamentals of fruiting and vegetative growth including lignocellulolysis.

A genetic linkage map of Pleurotus tuoliensis integrated with physical mapping of the de novo sequenced genome and the mating type loci

Background

Pleurotus tuoliensis (Bailinggu) is a commercially cultivated mushroom species with an increasing popularity in China and other Asian countries. Commercial profits are now low, mainly due to a low yield, long cultivation period and sensitivity to diseases. Breeding efforts are thus required to improve agronomical important traits. Developing saturated genetic linkage and physical maps is a start for applying genetic and molecular approaches to accelerate the precise breeding programs.

Results

Here we present a genetic linkage map for P. tuoliensis constructed by using 115 haploid monokaryons derived from a hybrid strain H6. One thousand one hundred and eighty-two SNP markers developed by 2b–RAD (type IIB restriction-site associated DNA) approach were mapped to 12 linkage groups. The map covers 1073 cM with an average marker spacing of 1.0 cM. The genome of P. tuoliensis was de novo sequenced as 40.8 Mb and consisted of 500 scaffolds (>500 bp), which showed a high level of colinearity to the genome of P. eryngii var. eryngii. A total of 97.4% SNP markers (1151) were physically localized on 78 scaffolds, and the physical length of these anchored scaffolds were 33.9 Mb representing 83.1% of the whole genome. Mating type loci A and B were mapped on separate linkage groups and identified physically on the assembled genomes. Five putative pheromone receptors and two putative pheromone precursors were identified for the mating type B locus.

Conclusions

This study reported a first genetic linkage map integrated with physical mapping of the de novo sequenced genome and the mating type loci of an important cultivated mushroom in China, P. tuoliensis. The de novo sequenced and annotated genome, assembled using a 2b–RAD generated linkage map, provides a basis for marker-assisted breeding of this economic important mushroom species.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4421-z) contains supplementary material, which is available to authorized users.

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.

Somatic transposition and meiotically driven elimination of an active helitron family in Pleurotus ostreatus

Helitrons constitute a superfamily of DNA transposons that were discovered in silico and are widespread in most eukaryotic genomes. They are postulated to mobilize through a "rolling-circle" mechanism, but the experimental evidence of their transposition has been described only recently. Here, we present the inheritance patterns of HELPO1 and HELPO2 helitron families in meiotically derived progeny of the basidiomycete Pleurotus ostreatus. We found distorted segregation patterns of HELPO2 helitrons that led to a strong under-representation of these elements in the progeny. Further investigation of HELPO2 flanking sites showed that gene conversion may contribute to the elimination of such repetitive elements in meiosis, favouring the presence of HELPO2 vacant loci. In addition, the analysis of HELPO2 content in a reconstructed pedigree of subclones maintained under different culture conditions revealed an event of helitron somatic transposition. Additional analyses of genome and transcriptome data indicated that P. ostreatus carries active RNAi machinery that could be involved in the control of transposable element proliferation. Our results provide the first evidence of helitron mobilization in the fungal kingdom and highlight the interaction between genome defence mechanisms and invasive DNA.

Construction of a genetic linkage map and analysis of quantitative trait loci associated with the agronomically important traits of Pleurotus eryngii

Breeding new strains with improved traits is a long-standing goal of mushroom breeders that can be expedited by marker-assisted selection (MAS). We constructed a genetic linkage map of Pleurotus eryngii based on segregation analysis of markers in postmeiotic monokaryons from KNR2312. In total, 256 loci comprising 226 simple sequence-repeat (SSR) markers, 2 mating-type factors, and 28 insertion/deletion (InDel) markers were mapped. The map consisted of 12 linkage groups (LGs) spanning 1047.8cM, with an average interval length of 4.09cM. Four independent populations (Pd3, Pd8, Pd14, and Pd15) derived from crossing between four monokaryons from KNR2532 as a tester strain and 98 monokaryons from KNR2312 were used to characterize quantitative trait loci (QTL) for nine traits such as yield, quality, cap color, and earliness. Using composite interval mapping (CIM), 71 QTLs explaining between 5.82% and 33.17% of the phenotypic variations were identified. Clusters of more than five QTLs for various traits were identified in three genomic regions, on LGs 1, 7 and 9. Regardless of the population, 6 of the 9 traits studied and 18 of the 71 QTLs found in this study were identified in the largest cluster, LG1, in the range from 65.4 to 110.4cM. The candidate genes for yield encoding transcription factor, signal transduction, mycelial growth and hydrolase are suggested by using manual and computational analysis of genome sequence corresponding to QTL region with the highest likelihood odds (LOD) for yield. The genetic map and the QTLs established in this study will help breeders and geneticists to develop selection markers for agronomically important characteristics of mushrooms and to identify the corresponding genes.

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.

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.

Highly expressed captured genes and cross-kingdom domains present in Helitrons create novel diversity in Pleurotus ostreatus and other fungi

BACKGROUND

Helitrons are class-II eukaryotic transposons that transpose via a rolling circle mechanism. Due to their ability to capture and mobilize gene fragments, they play an important role in the evolution of their host genomes. We have used a bioinformatics approach for the identification of helitrons in two Pleurotus ostreatus genomes using de novo detection and homology-based searching. We have analyzed the presence of helitron-captured genes as well as the expansion of helitron-specific helicases in fungi and performed a phylogenetic analysis of their conserved domains with other representative eukaryotic species.

RESULTS

Our results show the presence of two helitron families in P. ostreatus that disrupt gene colinearity and cause a lack of synteny between their genomes. Both putative autonomous and non-autonomous helitrons were transcriptionally active, and some of them carried highly expressed captured genes of unknown origin and function. In addition, both families contained eukaryotic, bacterial and viral domains within the helitron's boundaries. A phylogenetic reconstruction of RepHel helicases using the Helitron-like and PIF1-like helicase conserved domains revealed a polyphyletic origin for eukaryotic helitrons.

CONCLUSION

P. ostreatus helitrons display features similar to other eukaryotic helitrons and do not tend to capture host genes or gene fragments. The occurrence of genes probably captured from other hosts inside the helitrons boundaries pose the hypothesis that an ancient horizontal transfer mechanism could have taken place. The viral domains found in some of these genes and the polyphyletic origin of RepHel helicases in the eukaryotic kingdom suggests that virus could have played a role in a putative lateral transfer of helitrons within the eukaryotic kingdom. The high similarity of some helitrons, along with the transcriptional activity of its RepHel helicases indicates that these elements are still active in the genome of P. ostreatus.

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.

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.

Non-additive transcriptional profiles underlie dikaryotic superiority in Pleurotus ostreatus laccase activity

Background

The basidiomycete Pleurotus ostreatus is an efficient producer of laccases, a group of enzymes appreciated for their use in multiple industrial processes. The aim of this study was to reveal the molecular basis of the superiority of laccase production by dikaryotic strains compared to their parental monokaryons.

Methodology/Principal Findings

We bred and studied a set of dikaryotic strains starting from a meiotic population of monokaryons. We then completely characterised the laccase allelic composition, the laccase gene expression and activity profiles in the dikaryotic strain N001, in two of its meiotic full-sib monokaryons and in the dikaryon formed from their mating.

Conclusions/Significance

Our results suggested that the dikaryotic superiority observed in laccase activity was due to non-additive transcriptional increases in lacc6 and lacc10 genes. Furthermore, the expression of these genes was divergent in glucose- vs. lignocellulose-supplemented media and was highly correlated to the detected extracellular laccase activity. Moreover, the expression profile of lacc2 in the dikaryotic strains was affected by its allelic composition, indicating a putative single locus heterozygous advantage.

Induction of laccase activity in the white rot fungus Pleurotus ostreatus using water polluted with wheat straw extracts

The purpose of this work was to explore the use of polluted water effluents from wheat straw using industries as inducers of lignocellulolytic enzymatic activities in cultures of white rot basidiomycetes. For this purpose, we studied the effect of a wheat straw water extract on the evolution of the laccase activity recovered from submerged cultures of Pleurotus ostreatus made in different media and under various culture conditions. Our results demonstrated an accumulative induction effect in all the cultures and conditions tested. This induction is parallel to changes in the laccase electrophoretic profiles recovered from the culture supernatants. The isoenzyme that appeared to be mainly responsible for the laccase activity under these conditions was laccase 10, as confirmed by sequencing the induced protein. These results support the idea of using wheat straw effluents as inducers in liquid cultures of P. ostreatus mycelia for the production of ligninolytic enzymatic cocktails.

Molecular characterization of sexual diversity in a population of Serpula lacrymans, a tetrapolar basidiomycete

Different mating systems have evolved in the fungal kingdom, including a tetrapolar multiallelic mating system in many basidiomycetes. In tetrapolar species, the presence of different alleles at two mating loci (MAT A and MAT B) is necessary for mating to occur. The tetrapolar fungus Serpula lacrymans causes wood-decay in buildings in temperate regions worldwide and is present in Europe with a genetically homogeneous founder population. Using genome sequence data, we annotated the two mating type loci for S. lacrymans and found the expected synteny with other basidiomycetes, except for a retrotransposon being present in one locus (MAT A). We developed markers linked to the MAT A and MAT B regions and used these to investigate the mating type diversity in the European population. Moreover, we found a good match between the genetic markers and functional mating types as revealed by segregation and mating studies. A low diversity of mating types is present in the European S. lacrymans population caused by the founder event where a limited number of genotypes were introduced. This finding contrasts the situation in natural fungal populations where a high diversity of mating types is normally present. Although S. lacrymans has a large and viable population in Europe, we argue that the low mating type diversity restrains the dispersal and establishment of the fungus.

Fungal laccase, manganese peroxidase and lignin peroxidase: gene expression and regulation

Extensive research efforts have been dedicated to characterizing expression of laccases and peroxidases and their regulation in numerous fungal species. Much attention has been brought to these enzymes broad substrate specificity resulting in oxidation of a variety of organic compounds which brings about possibilities of their utilization in biotechnological and environmental applications. Research attempts have resulted in increased production of both laccases and peroxidases by the aid of heterologous and homologous expression. Through analysis of promoter regions, protein expression patterns and culture conditions manipulations it was possible to compare and identify common pathways of these enzymes' production and secretion. Although laccase and peroxidase proteins have been crystallized and thoroughly analyzed, there are still a lot of questions remaining about their evolutionary origin and the physiological functions. This review describes the present understanding of promoter sequences and correlation between the observed regulatory effects on laccase, manganese peroxidase and lignin peroxidase genes transcript levels and the presence of specific response elements.

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.

Transcriptional and enzymatic profiling of Pleurotus ostreatus laccase genes in submerged and solid-state fermentation cultures

The genome of the white rot basidiomycete Pleurotus ostreatus includes 12 phenol oxidase (laccase) genes. In this study, we examined their expression profiles in different fungal strains under different culture conditions (submerged and solid cultures) and in the presence of a wheat straw extract, which was used as an inducer of the laccase gene family. We used a reverse transcription-quantitative PCR (RT-qPCR)-based approach and focused on determining the reaction parameters (in particular, the reference gene set for the normalization and reaction efficiency determinations) used to achieve an accurate estimation of the relative gene expression values. The results suggested that (i) laccase gene transcription is upregulated in the induced submerged fermentation (iSmF) cultures but downregulated in the solid fermentation (SSF) cultures, (ii) the Lacc2 and Lacc10 genes are the main sources of laccase activity in the iSmF cultures upon induction with water-soluble wheat straw extracts, and (iii) an additional, as-yet-uncharacterized activity (Unk1) is specifically induced in SSF cultures that complements the activity of Lacc2 and Lacc10. Moreover, both the enzymatic laccase activities and the Lacc gene family transcription profiles greatly differ between closely related strains. These differences can be targeted for biotechnological breeding programs for enzyme production in submerged fermentation reactors.

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.

A genomic map enriched for markers linked to Avr1 in Cronartium quercuum f.sp. fusiforme

A novel approach is presented to map avirulence gene Avr1 in the basidiomycete Cronartium quercuum f.sp. fusiforme, the causal agent of fusiform rust disease in pines. DNA markers tightly linked to resistance gene Fr1 in loblolly pine tree 10-5 were used to classify 10-5 seedling progeny as either resistant or susceptible. A single dikaryotic isolate (P2) heterozygous at the corresponding Avr1 gene was developed by crossing Fr1 avirulent isolate SC20-21 with Fr1 virulent isolate NC2-40. Bulk basidiospore inoculum derived from isolate P2 was used to challenge the pine progeny. The ability to unambiguously marker classify 10-5 progeny as resistant (selecting for virulence) or susceptible (non-selecting) permitted the genetic mapping of the corresponding Avr1 gene by bulked segregant analysis. Using this approach, 14 genetic markers significantly linked to Avr1 were identified and placed within the context of a genome-wide linkage map produced for isolate P2 using samples from susceptible seedlings.

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.

A genetic linkage map of Venturia inaequalis, the causal agent of apple scab

BACKGROUND

Venturia inaequalis is an economically-important disease of apple causing annual epidemics of scab worldwide. The pathogen is a heterothallic ascomycete with an annual cycle of sexual reproduction on infected apple leaf litter, followed by several cycles of asexual reproduction during the apple growing season. Current disease control is achieved mainly through scheduled applications of fungicides. Genetic linkage maps are essential for studying genome structure and organisation, and are a valuable tool for identifying the location of genes controlling important traits of interest such as avirulence, host specificity and mating type in V. inaequalis. In this study, we performed a wide cross under in vitro conditions between an isolate of V. inaequalis from China and one from the UK to obtain a genetically diverse mapping population of ascospore progeny isolates and produced a map using AFLP and microsatellite (SSR) markers.

FINDINGS

Eighty-three progeny were obtained from the cross between isolates C0154 (China) x 01/213 (UK). The progeny was screened with 18 AFLP primer combinations and 31 SSRs, and scored for the mating type locus MAT. A linkage map was constructed consisting of 294 markers (283 AFLPs, ten SSRs and the MAT locus), spanning eleven linkage groups and with a total map length of 1106 cM. The length of individual linkage groups ranged from 30.4 cM (Vi-11) to 166 cM (Vi-1). The number of molecular markers per linkage group ranged from 7 on Vi-11 to 48 on Vi-3; the average distance between two loci within each group varied from 2.4 cM (Vi-4) to 7.5 cM (Vi-9). The maximum map length between two markers within a linkage group was 15.8 cM. The MAT locus was mapped to a small linkage group and was tightly linked to two AFLP markers. The map presented is over four times longer than the previously published map of V. inaequalis which had a total genetic distance of just 270 cM.

CONCLUSION

A genetic linkage map is an important tool for investigating the genetics of important traits in V. inaequalis such as virulence factors, aggressiveness and mating type. The linkage map presented here represents a significant improvement over currently published maps for studying genome structure and organisation, and for mapping genes of economic importance on the V. inaequalis genome.

A genetic linkage map of Pleurotus pulmonarius based on AFLP markers, and localization of the gene region for the sporeless mutation

In the cultivation of edible mushrooms, including Pleurotus pulmonarius (Fr.) Quel., the enormous number of spores produced by fruiting bodies can adversely affect mushroom growers' health, mushroom cultivation facilities, and the genetic diversity of natural populations. In this study, we constructed a primary genetic linkage map and identified the locus associated with the sporulation-deficient (sporeless) mutation of P. pulmonarius using 150 progeny isolates derived from a cross between sporeless and wild-type isolates. Based on the segregation of 300 AFLP markers, two mating-type factors, and the sporeless trait, a linkage map was generated consisting of 12 linkage groups. The map covered a total genetic distance of 971 cM, with an average marker interval of 5.2 cM. The gene region responsible for the sporeless mutation was located in linkage group II including 40 AFLP markers and the A mating-type factor locus. Of these markers, the nearest marker to the sporeless locus was located 1.4 cM away. Construction of this P. pulmonarius genetic linkage map and identification of markers that are closely linked to the sporeless locus will facilitate marker-assisted selective breeding of a sporeless strain with economically important traits.

Genetic linkage map for Amylostereum areolatum reveals an association between vegetative growth and sexual and self-recognition

Amylostereum areolatum is a filamentous fungus that grows through tip extension, branching and hyphal fusion. In the homokaryotic phase, the hyphae of different individuals are capable of fusing followed by heterokaryon formation, only if they have dissimilar allelic specificities at their mating-type (mat) loci. In turn, hyphal fusion between heterokaryons persists only when they share the same alleles at all of their heterokaryon incompatibility (het) loci. In this study we present the first genetic linkage map for A. areolatum, onto which the mat and het loci, as well as quantitative trait loci (QTLs) for mycelial growth rate are mapped. The recognition loci (mat-A and het-A) are positioned near QTLs associated with mycelial growth, suggesting that the genetic determinants influencing recognition and growth rate in A. areolatum are closely associated. This was confirmed when isolates associated with specific mat and het loci displayed significantly different mycelial growth rates. Although the link between growth and sexual recognition has previously been observed in other fungi, this is the first time that an association between growth and self-recognition has been shown.

Telomere organization in the ligninolytic basidiomycete Pleurotus ostreatus

Telomeres are structural and functional chromosome regions that are essential for the cell cycle to proceed normally. They are, however, difficult to map genetically and to identify in genome-wide sequence programs because of their structure and repetitive nature. We studied the telomeric and subtelomeric organization in the basidiomycete Pleurotus ostreatus using a combination of molecular and bioinformatics tools that permitted us to determine 19 out of the 22 telomeres expected in this fungus. The telomeric repeating unit in P. ostreatus is TTAGGG, and the numbers of repetitions of this unit range between 25 and 150. The mapping of the telomere restriction fragments to linkage groups 6 and 7 revealed polymorphisms compatible with those observed by pulsed field gel electrophoresis separation of the corresponding chromosomes. The subtelomeric regions in Pleurotus contain genes similar to those described in other eukaryotic systems. The presence of a cluster of laccase genes in chromosome 6 and a bipartite structure containing a Het-related protein and an alcohol dehydrogenase are especially relevant; this bipartite structure is characteristic of the Pezizomycotina fungi Neurospora crassa and Aspergillus terreus. As far as we know, this is the first report describing the presence of such structures in basidiomycetes and the location of a laccase gene cluster in the subtelomeric region, where, among others, species-specific genes allowing the organism to adapt rapidly to the environment usually map.

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

A genetic linkage map for the ectomycorrhizal basidiomycete Laccaria bicolor was constructed from 45 sib-homokaryotic haploid mycelial lines derived from the parental S238N strain progeny. For map construction, 294 simple sequence repeats (SSRs), single-nucleotide polymorphisms (SNPs), amplified fragment length polymorphisms (AFLPs) and random amplified polymorphic DNA (RAPD) markers were employed to identify and assay loci that segregated in backcross configuration. Using SNP, RAPD and SSR sequences, the L. bicolor whole-genome sequence (WGS) assemblies were aligned onto the linkage groups. A total of 37.36 Mbp of the assembled sequences was aligned to 13 linkage groups. Most mapped genetic markers used in alignment were colinear with the sequence assemblies, indicating that both the genetic map and sequence assemblies achieved high fidelity. The resulting matrix of recombination rates between all pairs of loci was used to construct an integrated linkage map using JoinMap. The final map consisted of 13 linkage groups spanning 812 centiMorgans (cM) at an average distance of 2.76 cM between markers (range 1.9-17 cM). The WGS and the present linkage map represent an initial step towards the identification and cloning of quantitative trait loci associated with development and functioning of the ectomycorrhizal symbiosis.

Construction of a genetic linkage map of the fungal pathogen of banana Mycosphaerella fijiensis, causal agent of black leaf streak disease

A genetic linkage map of the fungal plant pathogen Mycosphaerella fijiensis, causal agent of black leaf streak disease of banana was developed. A cross between the isolates CIRAD86 (from Cameroon) and CIRAD139A (from Colombia) was analyzed using molecular markers and the MAT locus. The genetic linkage map consists of 298 AFLP and 16 SSR markers with 23 linkage groups, containing five or more markers, covering 1,879 cM. Markers are separated on average by around 5.9 cM. The MAT locus was shown to segregate in a 1:1 ratio but could not be successfully mapped. An estimate of the relation between physical size and genetic distance was approximately 39.0 kb/cM. The estimated total haploid genome size was calculated using the genetic mapping data at 4,298.2 cM. This is the first genetic linkage map reported for this important foliar pathogen of banana. The great utility of the map will be for anchoring contigs in the genome sequence, evolutionary studies in comparison with other fungi, to identify quantitative trait loci (QTLs) associated with aggressiveness or oxidative stress resistance and with the recently available genome sequence, for positional cloning.

The first linkage map of the plant-pathogenic basidiomyceteTyphula ishikariensis

Speckled snow mold, caused by the basidiomycete Typhula ishikariensis Imai, is one of the most prominent winter diseases on perennial grasses and cereal crops in the northern hemisphere. The first linkage map of T. ishikariensis was constructed using a population of 93 sibling monokaryons derived from a single dikaryotic hybrid isolate that was created by a hyphal fusion of two monokaryotic parental isolates. The parental isolates were produced from a pathogenic dikaryotic isolate collected from a golf course in Wisconsin. The two parents exhibit significant differences in the production of aerial mycelium and sclerotia, and in their aggressiveness on creeping bentgrass ( Agrostis stolonifera L.). A total of 251 loci were mapped, comprising 89 inter-simple sequence repeat (ISSR) and 160 random amplified polymorphic DNA (RAPD) markers along with 2 phenotype-based mating-type (MAT) loci. The MAT loci were mapped on linkage groups (LGs) 1 and 7. The markers were evenly distributed over 7 LGs, covering 436 cM with an average marker interval of 2.2 cM. Seven chromosomes were cytologically observed using germ tube bursting methods with acetocarmine staining. This reference linkage map of T. ishikariensis should provide a framework for the mapping of quantitatively controlled traits such as fungal growth, survival, and virulence/avirulence under low temperatures. The map should also be utilized for studying the genome organization of the cold-loving plant-pathogenic Typhula spp. and for comparative genome analysis among fungal taxa.