Analysis of Genetic Diversity and Population Structure of Wild Strains and Cultivars Using Genomic SSR Markers in Lentinula edodes

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.

Japanese “nameko” mushrooms (<i>Pholiota microspora</i>) produced via sawdust-based cultivation exhibit severe genetic bottleneck associated with a single founder

Pholiota microspora ("nameko" in Japanese) is one of the most common edible mushrooms, especially in Japan, where sawdust-based cultivation is the most dominant method accounting for 99% of the production. The current strains for sawdust cultivation in Japan are considered to have been derived from a single wild strain collected from Fukushima, Japan, implying that commercial nameko mushrooms are derived from a severe genetic bottleneck. We tested this single founder hypothesis by developing 14 microsatellite markers for P. microspora to evaluate the genetic diversity of 50 cultivars and 73 wild strains isolated from across Japan. Microsatellite analysis demonstrated that sawdust-cultivated strains from Japan were significantly less genetically diverse than the wild strains, and the former displayed a significant bottleneck signature. Analyzing the genetic relationships among all genotypes also revealed that the sawdust-cultivated samples clustered into one monophyletic subgroup. Moreover, the sawdust-cultivated samples in Japan were more closely related than full-sibs. These results were consistent with the single founder hypothesis that suggests that all commercial nameko mushrooms produced in Japan are descendants of a single ancestor. Therefore, we conclude that cultivated P. microspora originated from a single domestication event that substantially reduced the diversity of commercial nameko mushrooms in Japan.

Chromosomal genome and population genetic analyses to reveal genetic architecture, breeding history and genes related to cadmium accumulation in Lentinula edodes

BACKGROUND

Lentinula edodes (Berk.) is the second most productive mushroom in the world. It contains compounds effective for antiviral, antitumor, antioxidant and immune regulation. Although genomes have previously been reported for this species, a high-quality chromosome-level reference for L. edodes is unavailable. This hinders detailed investigation of population genetics, breeding history of strains and genes related to environmental stress responses.

RESULTS

A high-quality chromosome-level genome was constructed. We separated a monokaryon from protoplasts of the commercial L. edodes strain L808 and assembled the genome of L. edodes using PacBio long-read and Illumina short-read sequencing, along with the high-throughput chromatin conformation capture (Hi-C) technique. We assembled a 45.87 Mb genome, and 99% of the sequences were anchored onto 10 chromosomes. The contig and scaffold N50 length were 2.17 and 4.94 Mb, respectively. Over 96% of the complete Benchmarking Universal Single-Copy Orthologs (BUSCO) were identified, and 9853 protein-coding genes were predicted. We performed population genome resequencing using 34 wild strains and 65 commercial cultivars of L. edodes originating from China, Japan, the United States and Australia. Based on whole-genome variants, we showed substantial differences in the Chinese wild population, which divided into different branches according to the main areas of their geographical distribution. We also determined the breeding history of L. edodes at the molecular level, and demonstrated that the cultivated strains in China mainly originated from wild strains from China and Northeast Asia. Phenotypic analysis showed that 99 strains exhibited differences on the Cd accumulation. Three significant loci in the of L. edodes genome were identified using the genome-wide association study (GWAS) of Cd accumulation traits. Functional genes associated with Cd accumulation traits were related to DNA ligase and aminoacyl tRNA synthetase, indicating that DNA damage repair and in vivo protein translation may be responses to Cd stress.

CONCLUSIONS

A high-quality chromosome-level genome and population genetic data of L. edodes provide genetic resources for functional genomic, evolutionary and artificial breeding studies for L. edodes.

Evaluation of Genetic Diversity and Population Structure Analysis among Germplasm of Agaricus bisporus by SSR Markers

Agaricus bisporus is a popular edible mushroom that is cultivated worldwide. Due to its secondary homothallic nature, cultivated A. bisporus strains have low genetic diversity, and breeding novel strains is challenging. The aim of this study was to investigate the genetic diversity and population structure of globally collected A. bisporus strains using simple sequence repeat (SSR) markers. Agaricus bisporus strains were divided based on genetic distance-based groups and model-based subpopulations. The major allele frequency (MAF), number of genotypes (NG), number of alleles (NA), observed heterozygosity (HO), expected heterozygosity (HE), and polymorphic information content (PIC) were calculated, and genetic distance, population structure, genetic differentiation, and Hardy-Weinberg equilibrium (HWE) were assessed. Strains were divided into two groups by distance-based analysis and into three subpopulations by model-based analysis. Strains in subpopulations POP A and POP B were included in Group I, and strains in subpopulation POP C were included in Group II. Genetic differentiation between strains was 99%. Marker AB-gSSR-1057 in Group II and subpopulation POP C was confirmed to be in HWE. These results will enhance A. bisporus breeding programs and support the protection of genetic resources.