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

Medicinal Value, Genetic Diversity, and Genetic Relationship Analysis of Auricularia cornea (Agaricomycetes) Based on ITS, ISSR, and SRAP Markers

Wild resources of Auricularia cornea (A. polytricha) are abundant in China, and genetic diversity and genetic relationships analysis of A. cornea can provide basis for germplasm resource utilization and innovation and molecular marker-assisted breeding. In this study, 22 Auricularia strains collected were identified as A. cornea based on ITS sequence analysis, and its genetic diversity was examined by ISSR and SRAP markers. The results showed that a total of 415 bands were amplified by 11 selected ISSR primers, with an average amplification of 37.73 bands per primer, and the mean values of Ne, I, and H were 1.302, 0.368, and 0.219, respectively. A total of 450 bands were amplified by 10 SRAP primers, with an average of 45 bands per primer, and the average of Ne, I, and H were 1.263, 0.302, and 0.183, respectively. The unweighted pair-group method with arithmetic means analysis based on ISSR-SRAP marker data revealed that the genetic similarity coefficient between the tested strains was 0.73-0.97, and the strains could be divided into five groups at 0.742, which had a certain correlation with regional distribution. The results of PCOA and population structure analysis based on ISSR-SRAP data also produced similar results. These results demonstrate the genetic diversity and distinctness among wild A. cornea and provide a theoretical reference for the classification, breeding, germplasm innovation, utilization, and variety protection of A. cornea resources.

The polysaccharides from Auricularia auricula alleviate non-alcoholic fatty liver disease via modulating gut microbiota and bile acids metabolism

The polysaccharides from Auricularia auricula (AAPs), containing a large number of O-acetyl groups that are related to the physiological and biological properties, seem to be potential prebiotics like other edible fungus polysaccharides. In the present study, therefore, the alleviating effects of AAPs and deacetylated AAPs (DAAPs, prepared from AAPs by alkaline treatment) on nonalcoholic fatty liver disease (NAFLD) induced by high-fat and high-cholesterol diet combined with carbon tetrachloride were investigated. The results revealed that both AAPs and DAAPs could effectively relieve liver injury, inflammation and fibrosis, and maintain intestinal barrier function. Both AAPs and DAAPs could modulate the disorder of gut microbiota and altered the composition of gut microbiota with enrichment of Odoribacter, Lactobacillus, Dorea and Bifidobacterium. Further, the alteration of gut microbiota, especially enhancement of Lactobacillus and Bifidobacterium, was contributed to the changes of bile acids (BAs) profile with increased deoxycholic acid (DCA). Farnesoid X receptor could be activated by DCA and other unconjugated BAs, which participated the BAs metabolism and alleviated the cholestasis, then protected against hepatitis in NAFLD mice. Interestingly, it was found that the deacetylation of AAPs negatively affected the anti-inflammation, thereby reducing the health benefits of A. auricula-derived polysaccharides.

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.

Analysis of the Genome Sequence of Strain GiC-126 of Gloeostereum incarnatum with Genetic Linkage Map

Gloeostereum incarnatum has edible and medicinal value and was first cultivated and domesticated in China. We sequenced the G. incarnatum monokaryotic strain GiC-126 on an Illumina HiSeq X Ten system and obtained a 34.52-Mb genome assembly sequence that encoded 16,895 predicted genes. We combined the GiC-126 genome with the published genome of G. incarnatum strain CCMJ2665 to construct a genetic linkage map (GiC-126 genome) that had 10 linkage groups (LGs), and the 15 assembly sequences of CCMJ2665 were integrated into 8 LGs. We identified 1912 simple sequence repeat (SSR) loci and detected 700 genes containing 768 SSRs in the genome; 65 and 100 of them were annotated with gene ontology (GO) terms and KEGG pathways, respectively. Carbohydrate-active enzymes (CAZymes) were identified in 20 fungal genomes and annotated; among them, 144 CAZymes were annotated in the GiC-126 genome. The A mating-type locus (MAT-A) of G. incarnatum was located on scaffold885 at 38.9 cM of LG1 and was flanked by two homeodomain (HD1) genes, mip and beta-fg. Fourteen segregation distortion markers were detected in the genetic linkage map, all of which were skewed toward the parent GiC-126. They formed three segregation distortion regions (SDR1-SDR3), and 22 predictive genes were found in scaffold1920 where three segregation distortion markers were located in SDR1. In this study, we corrected and updated the genomic information of G. incarnatum. Our results will provide a theoretical basis for fine gene mapping, functional gene cloning, and genetic breeding the follow-up of G. incarnatum.

Prediction of hyaluronic acid target on sucrase-isomaltase (SI) with reverse docking and molecular dynamics simulations for inhibitors binding to SI

Auricularia cornea (E.) polysaccharide is an important component of A. cornea Ehrenb, a white mutant strain of Auricularia with biological activities, such as enhancement of human immune function and cancer prevention. The hyaluronic acids (HAs) are important components of the A. cornea polysaccharide and have extremely high medicinal value. In this study, we used HA to search the target protein sucrase-isomaltase (SI). In addition, we also performed molecular dynamics (MD) simulations to explore the binding of three inhibitors (HA, acarbose and kotalanol) to SI. The MD simulations indicated that the binding of the three inhibitors may induce the partial disappearance of α helix in residues 530–580. Hence, the hydrogen bond for Gly570-Asn572, which was near the catalytic base Asp471 in SI, was broken during the binding of the three inhibitors. We reveal a new inhibitor for SI and provide reasonable theoretical clues for inhibitor binding to SI.

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.

Genome Sequence Analysis of Auricularia heimuer Combined with Genetic Linkage Map

Auricularia heimuer is one of the most popular edible fungi in China. In this study, the whole genome of A. heimuer was sequenced on the Illumina HiSeq X system and compared with other mushrooms genomes. As a wood-rotting fungus, a total of 509 carbohydrate-active enzymes (CAZymes) were annotated in order to explore its potential capabilities on wood degradation. The glycoside hydrolases (GH) family genes in the A. heimuer genome were more abundant than the genes in the other 11 mushrooms genomes. The A. heimuer genome contained 102 genes encoding class III, IV, and V ethanol dehydrogenases. Evolutionary analysis based on 562 orthologous single-copy genes from 15 mushrooms showed that Auricularia formed an early independent branch of Agaricomycetes. The mating-type locus of A. heimuer was located on linkage group 8 by genetic linkage analysis. By combining the genome sequence analysis with the genetic linkage map, the mating-type locus of A. heimuer was located on scaffold45 and consisted of two subunits, α and β. Each subunit consisted of a pair of homeodomain mating-type protein genes HD1 and HD2. The mapping revealed conserved synteny at the whole mating-type loci and mirror symmetry relations near the β subunit between A. heimuer and Exidia glandulosa. This study proposed the potential for the bioethanol production by consolidated bioprocessing of A. heimuer. It will promote understanding of the lignocellulose degradation system and facilitate more efficient conversion of the agricultural wastes used for mushroom cultivation. It also will advance the research on the fruiting body development and evolution of A. heimuer.

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.

Comparative transcriptome analysis reveals relationship of three major domesticated varieties of Auricularia auricula-judae

Auricularia auricula-judae is an edible mushroom and a traditional medicine in China as well as the fourth largest cultivated mushroom species in the world. Here for the first time, we present comparative transcriptome analyses of the fruiting bodies of three morphologically distinguishable A. auricula-judae cultivated varieties (Wujin, smooth; Banjin, partially wrinkled; and Quanjin, fully wrinkled) collected from Jilin Province, China. Biological triplicates were performed to determine the expression levels of 13,937 unigenes. Among them, only 13 unigenes were annotated to A. auricula-judae, highlighting the lack of publicly available reference sequences for this economically important species. Principal component analysis (PCA) determined that the gene expression profile of Quanjin was unique when compared to those of Banjin and Wujin. Such relationships were further supported by analyses of annotated and unannotated unigenes, differentially expressed unigenes, gene ontology functions, and the family of peroxidase genes. Using the KEGG database, significant alternations in biological pathways were detected among the three cultivars. This work contributes a large set of A. auricula-judae sequences to public database, establishes the relationships among major cultivars, and provides molecular guidance for breeding and cultivation.