A global phylogenomic analysis of the shiitake genus Lentinula

10 p.-6 fig.
Lentinula is a broadly distributed group of fungi that contains the cultivated shiitake mushroom, L. edodes. We sequenced 24 genomes representing eight described species and several unnamed lineages of Lentinula from 15 countries on four continents. Lentinula comprises four major clades that arose in the Oligocene, three in the Americas and one in Asia–Australasia. To expand sampling of shiitake mushrooms, we assembled 60 genomes of L. edodes from China that were previously published as raw Illumina reads and added them to our dataset. Lentinula edodes sensu lato (s. lat.) contains three lineages that may warrant recognition as species, one including a single isolate from Nepal that is the sister group to the rest of L. edodes s. lat., a second with 20 cultivars and 12 wild isolates from China, Japan, Korea, and the Russian Far East, and a third with 28 wild isolates from China, Thailand, and Vietnam. Two additional lineages in China have arisen by hybridization among the second and third groups. Genes encoding cysteine sulfoxide lyase (lecsl) and γ-glutamyl transpeptidase (leggt), which are implicated in biosynthesis of the organosulfur flavor compound lenthionine, have diversified in Lentinula. Paralogs of both genes that are unique to Lentinula (lecsl 3 and leggt 5b) are coordinately up-regulated in fruiting bodies of L. edodes. The pangenome of L. edodes s. lat. contains 20,308 groups of orthologous genes, but only 6,438 orthogroups (32%) are shared among all strains, whereas 3,444 orthogroups (17%) are found only in wild populations, which should be targeted for conservation.
This research was supported by the NSF awards DEB-1456588 to D.H. and DEB-1442113 to A.R., the Department of Energy award 0000238431 to D.H., the GENOBIOREF (BIO2017-86559-R) project of the Spanish Ministry of Science and Innovation (cofinanced by FEDER funds) to F.J.R.-D., S.C., and A.T.M., grant PID2019-103901GB-I00 of the Spanish Ministry of Science and Innovation to P.F., and the NIH award R56AI146096 to A.R. J.L.S. and A.R. were supported by the HHMI and the Burroughs Wellcome Fund. The work (proposal: 10.46936/10.25585/60001048) conducted by the US Department of Energy JGI (https://ror.org/04xm1d337), a DOE Office of Science User Facility, was supported by the Office of Science of the US Department of Energy operated under contract no. DE-AC02-05CH11231.