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

Genome sequencing progenies of magic mushrooms (Psilocybe subaeruginosa) identifies tetrapolar mating and gene duplications in the psilocybin pathway

Knowledge of breeding systems and genetic diversity is critical to select and combine desired traits that advance new cultivars in agriculture and horticulture. Mushrooms that produce psilocybin, magic mushrooms, may potentially be used in therapeutic and wellness industries, and stand to benefit from genetic improvement. We studied haploid siblings of Psilocybe subaeruginosa to resolve the genetics behind mating compatibility and advance knowledge of breeding. Our results show that mating in P. subaeruginosa is tetrapolar, with compatibility controlled at a homeodomain locus with one copy each of HD1 and HD2, and a pheromone/receptor locus with four homologs of the receptor gene STE3. An additional two pheromone/receptor loci homologous to STE3 do not appear to regulate mating compatibility. Alleles in the psilocybin gene cluster did not vary among the five siblings and were likely homozygous in the parent. Psilocybe subaeruginosa and its relatives have three copies of PsiH genes but their impact on production of psilocybin and its analogues is unknown. Genetic improvement in Psilocybe will require access to genetic diversity from the centre of origin of different species, identification of genes behind traits, and strategies to avoid inbreeding depression.

Cytoplasmic Mixing, Not Nuclear Coexistence, Can Explain Somatic Incompatibility in Basidiomycetes

Nonself recognition leading to somatic incompatibility (SI) is commonly used by mycologists to distinguish fungal individuals. Despite this, the process remains poorly understood in basidiomycetes as all current models of SI are based on genetic and molecular research in ascomycete fungi. Ascomycete fungi are mainly found in a monokaryotic stage, with a single type of haploid nuclei, and only briefly during mating do two genomes coexist in heterokaryotic cells. The sister phylum, Basidiomycota, differs in several relevant aspects. Basidiomycete fungi have an extended heterokaryotic stage, and SI is generally observed between heterokaryons instead of between homokaryons. Additionally, considerable nuclear migration occurs during a basidiomycete mating reaction, introducing a nucleus into a resident homokaryon with cytoplasmic mixing limited to the fused or neighboring cells. To accommodate these differences, we describe a basidiomycete model for nonself recognition using post-translational modification, based on a reader-writer system as found in other organisms. This post-translational modification combined with nuclear migration allows for the coexistence of two genomes in one individual while maintaining nonself recognition during all life stages. Somewhat surprisingly, this model predicts localized cell death during mating, which is consistent with previous observations but differs from the general assumptions of basidiomycete mating. This model will help guide future research into the mechanisms behind basidiomycete nonself recognition.

Monokaryotic Pleurotus sapidus Strains with Intraspecific Variability of an Alkene Cleaving DyP-Type Peroxidase Activity as a Result of Gene Mutation and Differential Gene Expression

The basidiomycete Pleurotus sapidus produced a dye-decolorizing peroxidase (PsaPOX) with alkene cleavage activity, implying potential as a biocatalyst for the fragrance and flavor industry. To increase the activity, a daughter-generation of 101 basidiospore-derived monokaryons (MK) was used. After a pre-selection according to the growth rate, the activity analysis revealed a stable intraspecific variability of the strains regarding peroxidase and alkene cleavage activity of PsaPOX. Ten monokaryons reached activities up to 2.6-fold higher than the dikaryon, with MK16 showing the highest activity. Analysis of the PsaPOX gene identified three different enzyme variants. These were co-responsible for the observed differences in activities between strains as verified by heterologous expression in Komagataella phaffii. The mutation S371H in enzyme variant PsaPOX_high caused an activity increase alongside a higher protein stability, while the eleven mutations in variant PsaPOX_low resulted in an activity decrease, which was partially based on a shift of the pH optimum from 3.5 to 3.0. Transcriptional analysis revealed the increased expression of PsaPOX in MK16 as reason for the higher PsaPOX activity in comparison to other strains producing the same PsaPOX variant. Thus, different expression profiles, as well as enzyme variants, were identified as crucial factors for the intraspecific variability of the PsaPOX activity in the monokaryons.

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.

Inheritance of phenotypic traits in the progeny of a Ceratocystis interspecific cross

Ceratocystis fimbriata is a fungal plant pathogen that causes black rot on Ipomoea batatas. Based on inoculation studies on numerous tree species, the pathogen is known to be host specific. The closely related species, Ceratocystis manginecans, causes severe wilt on a broad range of tree hosts, including Mangifera indica, Acacia mangium and other leguminous tree species. The genetic factors underlying the pathogenicity and host specificity of Ceratocystis species have rarely been investigated. In this study, an F₁ population of 70 recombinant progeny from a cross between C. fimbriata and C. manginecans was generated and the inheritance of various phenotypic traits was investigated. Results showed that colony colour, growth rate, asexual spore production and aggressiveness to I. batatas and A. mangium are all quantitative traits with high levels of heritability. However, conidia production and aggressiveness appeared to be regulated by a small number of genes. No correlation could be found between aggressiveness and other phenotypic traits, suggesting that these are inherited independently. This is the first study to consider genetic inheritance of pathogenicity and host specificity in Ceratocystis species and the results will contribute, in future, to the identification of quantitative trait loci and candidate genes associated with the traits investigated.

Somatic incompatibility and genetic structure of fungal crops in sympatric Atta colombica and Acromyrmex echinatior leaf-cutting ants

Obligate mutualistic symbioses rely on mechanisms that secure host-symbiont commitments to maximize host benefits and prevent symbiont cheating. Previous studies showed that somatic incompatibilities correlate with neutral-marker-based genetic distances between fungal symbionts of Panamanian Acromyrmex leaf-cutting ants, but the extent to which this relationship applies more generally remained unclear. Here we showed that genetic distances accurately predicted somatic incompatibility for Acromyrmex echinatior symbionts irrespective of whether neutral microsatellites or AFLP markers were used, but that such correlations were weaker or absent in sympatric Atta colombica colonies. Further analysis showed that the symbiont clades maintained by A. echinatior and A. colombica were likely to represent separate gene pools, so that neutral markers were unlikely to be similarly correlated with incompatibility loci that have experienced different selection regimes. We suggest that evolutionarily derived claustral colony founding by Atta queens may have removed selection for strong incompatibility in Atta fungi, as this condition makes the likelihood of symbiont swaps much lower than in Acromyrmex, where incipient nests stay open because queens have to forage until the first workers emerge.

Crosses between monokaryons of Pleurotus sapidus or Pleurotus florida show an improved biotransformation of (+)-valencene to (+)-nootkatone

Several hundred monokaryotic and new dikaryotic strains derived thereof were established from (+)-valencene tolerant Pleurotus species. When grouped according to their growth rate on agar plates and compared to the parental of Pleurotus sapidus 69, the slowly growing monokaryons converted (+)-valencene more efficiently to the grapefruit flavour compound (+)-nootkatone. The fast growing monokaryons and the slow×slow and the fast×fast dikaryotic crosses showed similar or inferior yields. Some slow×fast dikaryons, however, exceeded the biotransformation capability of the parental dikaryon significantly. The activity of the responsible enzyme, lipoxygenase, showed a weak correlation with the yields of (+)-nootkatone indicating that the determination of enzyme activity using the primary substrate linoleic acid may be misleading in predicting the biotransformation efficiency. This exploratory study indicated that a classical genetics approach resulted in altered and partly improved terpene transformation capability (plus 60%) and lipoxygenase activity of the strains.

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.

Mutualism and asexual reproduction influence recognition genes in a fungal symbiont

Mutualism between microbes and insects is common and alignment of the reproductive interests of microbial symbionts with this lifestyle typically involves clonal reproduction and vertical transmission by insect partners. Here the Amylostereum fungus-Sirex woodwasp mutualism was used to consider whether their prolonged association and predominance of asexuality have affected the mating system of the fungal partner. Nucleotide information for the pheromone receptor gene rab1, as well as the translation elongation factor 1α gene and ribosomal RNA internal transcribed spacer region were utilized. The identification of rab1 alleles in Amylostereum chailletii and Amylostereum areolatum populations revealed that this gene is more polymorphic than the other two regions, although the diversity of all three regions was lower than what has been observed in free-living Agaricomycetes. Our data suggest that suppressed recombination might be implicated in the diversification of rab1, while no evidence of balancing selection was detected. We also detected positive selection at only two codons, suggesting that purifying selection is important for the evolution of rab1. The symbiotic relationship with their insect partners has therefore influenced the diversity of this gene and influenced the manner in which selection drives and maintains this diversity in A. areolatum and A. chailletii.

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.

Quantitative trait locus mapping of yield-related components and oligogenic control of the cap color of the button mushroom, Agaricus bisporus

As in other crops, yield is an important trait to be selected for in edible mushrooms, but its inheritance is poorly understood. Therefore, we have investigated the complex genetic architecture of yield-related traits in Agaricus bisporus through the mapping of quantitative trait loci (QTL), using second-generation hybrid progeny derived from a cross between a wild strain and a commercial cultivar. Yield, average weight per mushroom, number of fruiting bodies per m(2), earliness, and cap color were evaluated in two independent experiments. A total of 23 QTL were detected for 7 yield-related traits. These QTL together explained between 21% (two-flushes yield) and 59% (earliness) of the phenotypic variation. Fifteen QTL (65%) were consistent between the two experiments. Four regions underlying significant QTL controlling yield, average weight, and number were detected on linkage groups II, III, IV, and X, suggesting a pleiotropic effect or tight linkage. Up to six QTL were identified for earliness. The PPC1 locus, together with two additional genomic regions, explained up to 90% of the phenotypic variation of the cap color. Alleles from the wild parent showed beneficial effects for some yield traits, suggesting that the wild germ plasm is a valuable source of variation for several agronomic traits. Our results constitute a key step toward marker-assisted selection and provide a solid foundation to go further into the biological mechanisms controlling productive traits in the button mushroom.

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.

Gene expression associated with vegetative incompatibility in Amylostereum areolatum

In filamentous fungi, vegetative compatibility among individuals of the same species is determined by the genes encoded at the heterokaryon incompatibility (het) loci. The hyphae of genetically similar individuals that share the same allelic specificities at their het loci are able to fuse and intermingle, while different allelic specificities at the het loci result in cell death of the interacting hyphae. In this study, suppression subtractive hybridization (SSH) followed by pyrosequencing and quantitative reverse transcription PCR were used to identify genes that are selectively expressed when vegetatively incompatible individuals of Amylostereum areolatum interact. The SSH library contained genes associated with various cellular processes, including cell-cell adhesion, stress and defence responses, as well as cell death. Some of the transcripts encoded proteins that were previously implicated in the stress and defence responses associated with vegetative incompatibility. Other transcripts encoded proteins known to be associated with programmed cell death, but have not previously been linked with vegetative incompatibility. Results of this study have considerably increased our knowledge of the processes underlying vegetative incompatibility in Basidiomycetes in general and A. areolatum in particular.

Genetic analysis of growth, morphology and pathogenicity in the F(1) progeny of an interspecific cross between Fusarium circinatum and Fusarium subglutinans

Fusarium circinatum and Fusarium subglutinans are two distinct species in the Gibberella fujikuroi species complex. A genetic linkage map produced from an interspecific cross between these species was used to identify quantitative trait loci (QTLs) associated with variation in mycelial growth and morphology of colony margins (CMs) in the 94 F(1) progeny. Mycelial growth was assessed by measuring culture size at 25°C and 30°C, while CM morphology was characterized in the parents and assessed in their F(1) progeny. In order to test the pathogenicity of the progeny, Pinus patula seedlings were inoculated and lesion lengths were measured after 3weeks. Seven putative QTLs were associated with mycelial growth, three for growth at 25°C and four at 30°C. One highly significant QTL (P<0.001) was present at both growth temperatures. For CM morphology, a QTL was identified at the same position (P<0.001) as the QTL responsible for growth at the two temperatures. The putative QTLs accounted for 45 and 41% of the total mycelial growth variation at 25°C and 30°C, respectively, and for 21% of the variation in CM morphology. Only one of the 94 F(1) progeny was pathogenic on P. patula seedlings. This observation could be explained by the genetic constitution of this F(1) isolate, namely that ∼96% of its genome originated from the F. circinatum parent. This F(1) individual also grew significantly faster at 25°C than the F. circinatum parent (P<0.05), as well as more rapidly than the average growth for the remaining 93 F(1) progeny (P<0.05). However, no association was found between mycelial growth and pathogenicity at 25°C. The highly significant QTL associated with growth at two temperatures, suggests that this is a principal genomic region involved in mycelial growth at both temperatures, and that the same region is also responsible for CM morphology.