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

Molecular Profiling of Rice Straw Degradability Discrepancy in Stropharia rugosoannulata Core Germplasm

The rice-S. rugosoannulata pattern is a rapidly growing agricultural practice for straw disposal and mushroom production in China. However, different S. rugosoannulata strains show a large variation in rice straw degradability. Here, we constructed a core collection of S. rugosoannulata containing 14 strains with rich genetic diversity. The molecular profiling of the lignocellulose degradability discrepancy of S. rugosoannulata strains was then explored using enzyme activity assays and transcriptome analysis. The results indicated that mycelial growth rate, lignocellulolytic enzyme activities, and rice straw degradability differed widely among the S. rugosoannulata core strains. The genes encoding lignin modifying and degrading auxiliary enzymes, oxidases, glycoside hydrolases, and detoxification proteins were differentially expressed between two representative S. rugosoannulata strains, resulting in differences in their lignocellulolytic enzyme activities and further causing differences in lignocellulose degradability. This study is useful to improve the production efficiency of S. rugosoannulata and promote the recycling of rice straw.

Near-gapless genome and transcriptome analyses provide insights into fruiting body development in Lentinula edodes

Fruiting body development in macrofungi is an intensive research subject. In this study, high-quality genomes were assembled for two sexually compatible monokaryons from a heterokaryotic Lentinula edodes strain WX1, and variations in L. edodes genomes were analyzed. Specifically, differential gene expression and allele-specific expression (ASE) were analyzed using the two monokaryotic genomes and transcriptome data from four different stages of fruiting body development in WX1. Results revealed that after aeration, mycelia sensed cell wall stress, pheromones, and a decrease in CO2 concentration, leading to up-regulated expression in genes related to cell adhesion, cell wall remodeling, proteolysis, and lipid metabolism, which may promote primordium differentiation. Aquaporin genes and those related to proteolysis, mitosis, lipid, and carbohydrate metabolism may play important roles in primordium development, while genes related to tissue differentiation and sexual reproduction were active in fruiting body. Several essential genes for fruiting body development were allele-specifically expressed and the two nuclear types could synergistically regulate fruiting body development by dominantly expressing genes with different functions. ASE was probably induced by long terminal repeat-retrotransposons. Findings here contribute to the further understanding of the mechanism of fruiting body development in macrofungi.

Dicarboxylic Amino Acid Permease 7219 Regulates Fruiting Body Type of Auricularia heimuer

Auricularia heimuer is a widely cultivated jelly mushroom. The fruiting bodies are categorized into cluster and chrysanthemum types. With changing consumer demands and the need to reduce bio-waste, the demand for clustered fruiting bodies is increasing. Therefore, gene mining for fruiting body types is a matter of urgency. We determined that the A. heimuer locus for fruiting body type was located at one end of the genetic linkage map. The locus was localized between the markers D23860 and D389 by increasing the density of the genetic linkage map. BlastN alignment showed that the marker SCL-18 was also located between D23860 and D389, and a total of 25 coding genes were annotated within this interval. Through parental transcriptome analysis and qRT-PCR verification, the locus g7219 was identified as the gene controlling the fruiting body type. A single-nucleotide substitution in the TATA box of g7219 was detected between the parents. By PCR amplification of the promoter region of g7219, the TATA-box sequences of the cluster- and chrysanthemum-type strains were found to be CATAAAA and TATAAAA, respectively. This study provides a foundation for the breeding of fruiting body types and strain improvement of A. heimuer.

Determination and Analysis of Hyper-Variable A Mating Types in Wild Strains of Lentinula edodes in Korea

The diversity of A mating type in wild strains of Lentinula edodes was extensively analyzed to characterize and utilize them for developing new cultivars. One hundred twenty-three A mating type alleles, including 67 newly discovered alleles, were identified from 106 wild strains collected for the past four decades in Korea. Based on previous studies and current findings, a total of 130 A mating type alleles have been found, 124 of which were discovered from wild strains, indicating the hyper-variability of A mating type alleles of L. edodes. About half of the A mating type alleles in wild strains were found in more than two strains, whereas the other half of the alleles were found in only one strain. About 90% of A mating type combinations in dikaryotic wild strains showed a single occurrence. Geographically, diverse A mating type alleles were intensively located in the central region of the Korean peninsula, whereas only allele A17 was observed throughout Korea. We also found the conservation of the TCCCAC motif in addition to the previously reported motifs, including ATTGT, ACAAT, and GCGGAG, in the intergenic regions of A mating loci. Sequence comparison among some alleles indicated that accumulated mutation and recombination would contribute to the diversification of A mating type alleles in L. edodes. Our data support the rapid evolution of A mating locus in L. edodes, and would help to understand the characteristics of A mating loci of wild strains in Korea and help to utilize them for developing new cultivars.

Chromosome-Wide Characterization of Intragenic Crossover in Shiitake Mushroom, Lentinula edodes

Meiotic crossover plays a critical role in generating genetic variations and is a central component of breeding. However, our understanding of crossover in mushroom-forming fungi is limited. Here, in Lentinula edodes, we characterized the chromosome-wide intragenic crossovers, by utilizing the single-nucleotide polymorphisms (SNPs) datasets of an F₁ haploid progeny. A total of 884 intragenic crossovers were identified in 110 single-spore isolates, the majority of which were closer to transcript start sites. About 71.5% of the intragenic crossovers were clustered into 65 crossover hotspots. A 10 bp motif (GCTCTCGAAA) was significantly enriched in the hotspot regions. Crossover frequencies around mating-type A (MAT-A) loci were enhanced and formed a hotspot in L. edodes. Genome-wide quantitative trait loci (QTLs) mapping identified sixteen crossover-QTLs, contributing 8.5–29.1% of variations. Most of the detected crossover-QTLs were co-located with crossover hotspots. Both cis- and trans-QTLs contributed to the nonuniformity of crossover along chromosomes. On chr2, we identified a QTL hotspot that regulated local, global crossover variation and crossover hotspot in L. edodes. These findings and observations provide a comprehensive view of the crossover landscape in L. edodes, and advance our understandings of conservation and diversity of meiotic recombination in mushroom-forming fungi.

Macrosynteny analysis between Lentinula edodes and Lentinula novae-zelandiae reveals signals of domestication in Lentinula edodes

The basidiomycete fungus Lentinula novae-zelandiae is endemic to New Zealand and is a sister taxon to Lentinula edodes, the second most cultivated mushroom in the world. To explore the biology of this organism, a high-quality chromosome level reference genome of L. novae-zelandiae was produced. Macrosyntenic comparisons between the genome assembly of L. novae-zelandiae, L. edodes and a set of three genome assemblies of diverse species from the Agaricomycota reveal a high degree of macrosyntenic restructuring within L. edodes consistent with signal of domestication. These results show L. edodes has undergone significant genomic change during the course of its evolutionary history, likely a result of its cultivation and domestication over the last 1000 years.

RNA-Seq-based high-resolution linkage map reveals the genetic architecture of fruiting body development in shiitake mushroom, Lentinula edodes

•

We constructed a reference genetic map of Lentinula edodes.

•

We re-assembled a chromosome-level genome of L. edodes.

•

We disclosed three hotspots regions for fruiting body-related traits in shiitake.

•

We scanned candidate genes for fruiting body-related traits.

Fruiting body development (FBD) of mushroom-forming fungi has attracted tremendous interest. However, the genetic and molecular basis of FBD is poorly known. Here, using Lentinula edodes (shiitake) as a model, we deciphered the genetic architecture underlying fruiting body-related traits (FBRTs) by combined genomic, genetic and phenotypic data. Using RNA-Seq of fruiting bodies from 110 dikaryons in a bi-parental mapping population, we constructed an ultra-high-density genetic map of L. edodes (Lemap2.0) with a total length of 810.14 cM, which covered 81.7% of the shiitake genome. A total of 94 scaffolds of the shiitake genome were aligned to Lemap2.0 and re-anchored into nine pseudo-chromosomes. Then via quantitative trait locus (QTL) analysis, we disclosed an outline of the genetic architecture of FBD in shiitake. Twenty-nine QTLs and three main genomic regions associated with FBD of shiitake were identified. Using meta-QTL analysis, seven pleiotropic QTLs for multiple traits were detected, which contributed to the correlations of FBRTs. In the mapped QTLs, the expression of 246 genes were found to significantly correlate with the phenotypic traits. Thirty-three of them were involved in FBD and could represent candidate genes controlling the shape and size of fruiting bodies. Collectively, our findings have advanced our understanding of the genetic regulation of FBD in shiitake and mushroom-forming fungi at large.

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

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

Morphological, cytological and genetic analyzes of the 'sango' mutant with the defects in basidiocarp development in edible mushroom Pleurotus pulmonarius

A spontaneous, morphological variation 'sango' was observed in the progeny of a Pleurotus pulmonarius (Fr.) Quél. wild-type basidiocarp (also known as fruiting body) collected from the field. This variant developed wart- and coral-like structures instead of normal basidiocarps. Microscopic analysis showed that the sango phenotype had defects in the differentiation of the pileus and hymenium. Basidiocarp phenotypic data analysis in the progenies revealed that the sango trait is a heritable mutation character controlled by a single recessive gene. This mutation locus was mapped on linkage group III of a previously constructed genetic linkage map by amplified fragment length polymorphism (AFLP) technique in P. pulmonarius. Four AFLP markers identified by bulked segregant analysis showed linkage to the sango mutation locus, with the genetic distance ranging from 0 to 2.1 cM. Of these markers, one marker was co-segregated with the sango mutation locus. This knowledge will be a useful foundation for practical breeding as well as for elucidating molecular mechanisms in basidiocarp development of main edible mushrooms.

RNA-seq Profiling Showed Divergent Carbohydrate-Active Enzymes (CAZymes) Expression Patterns in Lentinula edodes at Brown Film Formation Stage Under Blue Light Induction

Lentinula edodes (shiitake mushroom) is one of the most important edible mushrooms worldwide. The L. edodes cultivation cycle includes a unique developing stage called brown film formation that directly affects the development of primordium and the quality of fruiting body. Brown film formation is induced by light, especially blue light. To promote our understanding of the role of blue light in brown film formation mechanisms of L. edodes, we used RNA-seq and compared the transcriptomes of L. edodes grown under blue light and in dark, and validated the expression profiles using qRT-PCR. Blue light stimulated the formation of brown film and increased the content of polysaccharides in L. edodes. Blue light also promoted L. edodes to absorb more polysaccharides by enhancing the activities of enzymes. Among the 730 differentially expressed genes (DEGs), 433 genes were up-regulated and 297 were down-regulated. Most of the DEGs were in the oxidoreductase activity group. Pentose and glucuronic acid conversion and starch and sucrose metabolism were the most important pathways in the formation of brown film. A total of 79 genes of DEGs were identified as genes encoding carbohydrate-active enzymes (CAZymes). Fifty-one of the CAZymes genes were up-regulated, suggesting that CAZymes play important roles in brown film formation to provide sufficient nutrition for L. edodes. The results will facilitate future functional investigations of the genes involved in the developmental control of L. edodes.

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.

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

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

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.

High-density genetic mapping identifies the genetic basis of a natural colony morphology mutant in the root rot pathogen Armillaria ostoyae

Filamentous fungi exhibit a broad spectrum of heritable growth patterns and morphological variations reflecting the adaptation of the different species to distinct ecological niches. But also within species, isolates show considerable variation in growth rates and other morphological characteristics. The genetic basis of this intraspecific variation in mycelial growth and morphology is currently poorly understood. By chance, a growth mutant in the root rot pathogen Armillaria ostoyae was discovered. The mutant phenotype was characterized by extremely compact and slow growth, as well as shorter aerial hyphae and hyphal compartments in comparison to the wildtype phenotype. Genetic analysis revealed that the abnormal phenotype is caused by a recessive mutation, which segregates asa single locus in sexual crosses. In order to identify the genetic basis of the mutant phenotype, we performed a quantitative trait locus (QTL) analysis. A mapping population of 198 haploid progeny was genotyped at 11,700 genome-wide single nucleotide polymorphisms (SNPs) making use of double digest restriction site associated DNA sequencing (ddRADseq). In accordance with the genetic analysis, a single significant QTL was identified for the abnormal growth phenotype. The QTL confidence interval spans a narrow, gene dense region of 87kb in the A. ostoyae genome which contains 37 genes. Overall, our study reports the first high-density genetic map for an Armillaria species and shows its successful application in forward genetics by resolving the genetic basis of a mutant phenotype with a severe defect in hyphal growth.

Association Mapping Reveals Genetic Loci Associated with Important Agronomic Traits in Lentinula edodes, Shiitake Mushroom

Association mapping is a robust approach for the detection of quantitative trait loci (QTLs). Here, by genotyping 297 genome-wide molecular markers of 89 Lentinula edodes cultivars in China, the genetic diversity, population structure and genetic loci associated with 11 agronomic traits were examined. A total of 873 alleles were detected in the tested strains with a mean of 2.939 alleles per locus, and the Shannon's information index was 0.734. Population structure analysis revealed two robustly differentiated groups among the Chinese L. edodes cultivars (F_ST = 0.247). Using the mixed linear model, a total of 43 markers were detected to be significantly associated with four traits. The number of markers associated with traits ranged from 9 to 26, and the phenotypic variations explained by each marker varied from 12.07% to 31.32%. Apart from five previously reported markers, the remaining 38 markers were newly reported here. Twenty-one markers were identified as simultaneously linked to two to four traits, and five markers were associated with the same traits in cultivation tests performed in two consecutive years. The 43 traits-associated markers were related to 97 genes, and 24 of them were related to 10 traits-associated markers detected in both years or identified previously, 13 of which had a >2-fold expression change between the mycelium and primordium stages. Our study has provided candidate markers for marker-assisted selection (MAS) and useful clues for understanding the genetic architecture of agronomic traits in the shiitake mushroom.

Population genomic analysis uncovers environmental stress-driven selection and adaptation of Lentinula edodes population in China

The elucidation of genome-wide variations could help reveal aspects of divergence, domestication, and adaptation of edible mushrooms. Here, we resequenced the whole genomes of 39 wild and 21 cultivated strains of Chinese Lentinula edodes, the shiitake mushroom. We identified three distinct genetic groups in the Chinese L. edodes population with robust differentiation. Results of phylogenetic and population structure analyses suggest that the cultivated strains and most of the wild trains of L. edodes in China possess different gene pools and two outlier strains show signatures of hybridization between groups. Eighty-four candidate genes contributing to population divergence were detected in outlier analysis, 18 of which are involved in response to environmental stresses. Gene enrichment analysis of group-specific single nucleotide polymorphisms showed that the cultivated strains were genetically diversified in biological processes related to stress response. As the formation of fruiting bodies is a stress-response process, we postulate that environment factors, such as temperature, drove the population divergence of L. edodes in China by natural or artificial selection. We also found phenotypic variations between groups and identified some wild strains that have potential to diversify the genetic pool for improving agricultural traits of L. edodes cultivars in China.

Genome Sequence of the Edible Cultivated Mushroom Lentinula edodes (Shiitake) Reveals Insights into Lignocellulose Degradation

Lentinula edodes, one of the most popular, edible mushroom species with a high content of proteins and polysaccharides as well as unique aroma, is widely cultivated in many Asian countries, especially in China, Japan and Korea. As a white rot fungus with lignocellulose degradation ability, L. edodes has the potential for application in the utilization of agriculture straw resources. Here, we report its 41.8-Mb genome, encoding 14,889 predicted genes. Through a phylogenetic analysis with model species of fungi, the evolutionary divergence time of L. edodes and Gymnopus luxurians was estimated to be 39 MYA. The carbohydrate-active enzyme genes in L. edodes were compared with those of the other 25 fungal species, and 101 lignocellulolytic enzymes were identified in L. edodes, similar to other white rot fungi. Transcriptome analysis showed that the expression of genes encoding two cellulases and 16 transcription factor was up-regulated when mycelia were cultivated for 120 minutes in cellulose medium versus glucose medium. Our results will foster a better understanding of the molecular mechanism of lignocellulose degradation and provide the basis for partial replacement of wood sawdust with agricultural wastes in L. edodes cultivation.

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.

Genetic dissection of fruiting body-related traits using quantitative trait loci mapping in Lentinula edodes

To provide a better understanding of the genetic architecture of fruiting body formation of Lentinula edodes, quantitative trait loci (QTLs) mapping was employed to uncover the loci underlying seven fruiting body-related traits (FBRTs). An improved L. edodes genetic linkage map, comprising 572 markers on 12 linkage groups with a total map length of 983.7 cM, was constructed by integrating 82 genomic sequence-based insertion-deletion (InDel) markers into a previously published map. We then detected a total of 62 QTLs for seven target traits across two segregating testcross populations, with individual QTLs contributing 5.5 %-30.2 % of the phenotypic variation. Fifty-three out of the 62 QTLs were clustered in six QTL hotspots, suggesting the existence of main genomic regions regulating the morphological characteristics of fruiting bodies in L. edodes. A stable QTL hotspot on MLG2, containing QTLs for all investigated traits, was identified in both testcross populations. QTLs for related traits were frequently co-located on the linkage groups, demonstrating the genetic basis for phenotypic correlation of traits. Meta-QTL (mQTL) analysis was performed and identified 16 mQTLs with refined positions and narrow confidence intervals (CIs). Nine genes, including those encoding MAP kinase, blue-light photoreceptor, riboflavin-aldehyde-forming enzyme and cyclopropane-fatty-acyl-phospholipid synthase, and cytochrome P450s, were likely to be candidate genes controlling the shape of fruiting bodies. The study has improved our understanding of the genetic architecture of fruiting body formation in L. edodes. To our knowledge, this is the first genome-wide QTL detection of FBRTs in L. edodes. The improved genetic map, InDel markers and QTL hotspot regions revealed here will assist considerably in the conduct of future genetic and breeding studies of L. edodes.

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.

Genetic diversity analysis of Hypsizygus marmoreus with target region amplification polymorphism

Hypsizygus marmoreus is an industrialized edible mushroom. In the present paper, the genetic diversity among 20 strains collected from different places of China was evaluated by target region amplification polymorphism (TRAP) analysis; the common fragment of TRAPs was sequenced and analyzed. Six fixed primers were designed based on the analysis of H. marmoreus sequences from GenBank database. The genomic DNA extracted from H. marmoreus was amplified with 28 TRAP primer combinations, which generated 287 bands. The average of amplified bands per primer was 10.27 (mean polymorphism is 69.73%). The polymorphism information content (PIC) value for TRAPs ranged from 0.32 to 0.50 (mean PIC value per TRAP primer combination is 0.48), which indicated a medium level of polymorphism among the strains. A total of 36 sequences were obtained from TRAP amplification. Half of these sequences could encode the known or unknown proteins. According to the phylogenetic analysis based on TRAP result, the 20 strains of H. marmoreus were classified into two main groups.