Skip to main content

Enhancement of ergothioneine production by discovering and regulating its metabolic pathway in Cordyceps militaris



Cordyceps militaris is a traditional medicinal fungus contains a variety of functional ingredients and has been developed as an important mushroom food recently. Ergothioneine, one of the antioxidative compounds in C. militaris, is benefits on aging-related diseases and therefore became a novel functional food nutritive fortifier. Currently, the main diet source of ergothioneine is mushroom food. However, the mushroom farming faces the problems such as rather low ingredient yield and spontaneous degeneration associated fruiting body that restricts large scale production of ergothioneine.


In this study, we excavated the ergothioneine synthetases in mushroom and modified the genes in C. militaris to construct a new ergothioneine synthesis pathway. By further introducing this pathway into C. militaris genome, we succeeded to increase the ingredients’ production of engineering strain, the highest amount of ergothioneine and cordycepin were up to 2.5 g/kg dry weight and 2 g/L, respectively. Additionally, the expression of ergothioneine synthetase genes in the shape-mutated degenerative C. militaris could recover the ability of degenerative strain to produce high amount of ingredients, suggesting the metabolic regulation of ergothioneine might release the symptom of mushroom degeneration.


This study reveals a new pathway to fulfill the market needs of functional mushroom food and food fortifier ergothioneine. It implied the mycelium of C. militaris could be engineered as a novel medicinal mushroom food which could produce higher amount of valuable ingredients.

Graphical abstract


Cordyceps militaris is a well-known edible and medicinal fungi for centuries in Asian [1], it has been developed as an important mushroom food recently [2]. It contains many valuable bioactive compounds [3, 4] such as cordycepin and ergothioneine. Ergothioneine, one of the ingredients with great antioxidant and anti-inflammatory activity [5] in C. militaris [6], is a natural sulfur-containing thiol molecular and has been regarded as longevity vitamins these years [7, 8]. The specific cationic transport protein (OCTN1) has high affinity for ergothioneine, suggesting the important role of ergothioneine in human physiological activities [9, 10]. But ergothioneine could not be synthesized by plants and animals, including human being. The major diet source of ergothioneine for human is mushroom [11] food, which contains the highest amount (average 150–727 mg/kg dry weight [12]) of ergothioneine than other species in nature. Recently, ergothioneine received self-affirmed Generally Recognized As Safe (GRAS) and approval as a food supplement by the European Food Safety Authority (EFSA) [13], which led it to be a promising new functional food nutritive fortifier.

To obtain higher yield of ergothioneine for the commercial market, engineered strains with different ergothioneine synthetases have been developed in model organism Escherichia coli [14,15,16,17], Saccharomyces cerevisiae [18,19,20], Aspergillus oryzae [21] and et al. [22, 23], but ergothioneine production could only reach 687 mg/L or 231 mg/kg media in their continuous fermentation. The production of ergothioneine in E. coli and S. cerevisiae can be further increased by strategies of metabolic engineering and fed-batch fermentation, the highest content of ergothioneine were reported to reach 1.31 g/L [15] and 2.39 g/L [20] respectively. However, even with such delicate optimization and high cost in fermentation, the ergothioneine content is not higher enough to fulfill the market needs. As the applications of high-throughput sequencing analysis [24, 25] and genome editing technologies [26] in Cordyceps, the genetic breeding strategy can be applicated in C. militaris to obtain strain with higher amount of ingredients, therefore C. militaris could be developed as a great engineered mushroom host to improve the production of ergothioneine and fulfill the medicinal mushroom food market.

The content of ergothioneine in wild type C. militaris was 382 to 799 mg/kg dry weight (DW) fruiting body [27] and 130 ± 11 mg/kg DW mycelia [6], respectively. The rather low content in mycelia indicated the production of ergothioneine in C. militaris could be further increased via genetic modification and metabolic engineering. Besides, as one of confirmed factors that caused spontaneous degradation of C. militaris is the cellular accumulation of reactive oxygen species (ROS) during the light-induced stage [28, 29], the construction of C. militaris strain with overexpression of ergothioneine synthetases is necessary for further studying the degeneration protective function of ergothioneine. However, the biosynthesis pathway of ergothioneine in C. militaris has not been discovered yet.

In this study, we identified the ergothioneine synthetases of several mushrooms. The function of C. militaris ergothioneine synthetases was verified in E. coli with a reconstructed de novo pathway. By introducing this pathway into C. militaris genome, we succeeded to increase the production of ergothioneine in engineering strains, the highest amount was 2485.32 ± 50.72 mg/kg DW of the 10 days fermentative mycelium. Additionally, we found that the overproduction of ergothioneine in fruiting body degenerated mutant (shape of fat and dwarf) C. militaris could recover the high productive ability (approximately 1 g/L) of its main bioactive ingredient cordycepin, while overexpression in wild type could further double the cordycepin content. In conclusion, the discovery and regulation of the ergothioneine metabolic pathway greatly improved the yield of ergothioneine and cordycepin in C. militaris, revealing a new way to meet the market needs of functional mushroom food and ergothioneine, and implied that C. militaris mycelium could be engineered as a suitable platform to produce higher amount of valuable ingredients.

Results and discussion

Conserved domain prediction and sequence alignment of edible mushrooms ergothioneine synthetases

Ergothioneine could be synthesized by several pathways [30, 31] in nature (Fig. 1A). These pathways shared similar catalytic route but with different numbers of synthetases (EgtA, B, C, D, E for bacteria [30] and Egt1,2 for fungi [32]). To excavate the biosynthesis pathway of ergothioneine in mushroom, we first analyzed the sequence characterization and functional domain of putative synthetases, which were found by homology sequence analysis. We used the amino acid sequence of NcEgt1/2 as probes to align the sulfoxide synthetases or PLP(pyridoxal 5-phosphatemonohydrate)-depend C-S lyase of ergothioneine with the known genomic sequence of mushrooms (including C. militaris, A. aegerita, A. bisporus, L. edodes, P. ostreatus, F. velutipes, G. frondosa and G. sinense).

As showed in Table 1, the amino acid sequence of Egt1 from edible mushrooms were showed identification of 33.75–56.43% with the amino acid sequence of NcEgt1, while the amino acid sequence of Egt2 were showed identification of 31.10–48.78% with NcEgt2. They all contained the same functional domains and shared similar sequence length to each other (Fig. 1B, C), which implied the edible mushrooms have a similar ergothioneine biosynthesis pathway with N. crassa.

Table 1 Putative ergothioneine biosynthesis genes blast from edible mushrooms by comparing with NcEgt1 and NcEgt2

We labeled the key residues for binding and catalysis of the putative Egt1s (Fig. 1E) and Egt2s (Fig. 1F) of edible mushrooms (Additional file 1: Table S1) based on previously identified motifs, which were responsible for enzymic functions [31, 33]. The amino acid sequences of Egt1s in different edible mushrooms were sharing the same iron binding sites and catalytic residues tyrosine. But the hercynine (TMH) binding sites of putative Egt1 genes from P. ostreatus, G. frondosa and G. sinense were distinguish with others, which suggested that these three species might use a different substrate rather than TMH to perform synthesis of hercynylcysteine sulfoxide (Cys-TMH).

Fig. 1
figure 1

Two different natural biosynthesis pathways of ergothioneine (A). A phylogenetic tree with schematic diagram of functional domains contains putative Egt1 (B) or Egt2 (C) homologous (boxes represented the predicted domains and labeled with corresponding superfamily names). A schematic diagram of genes structure (boxes represented exon and the black line between boxes represented intron) of NcEgt1/2 and CmEgt1/2 (D). Sequence alignment of putative Egt1 (E) or Egt2 (F) homologous in edible mushrooms (predicted key residues for binding were labeled)

Previous study using biochemical and X-ray crystallographic approaches to reveal the NcEgt2 contains a unique bisulfide binding site, but none of Egt2 extracted from edible mushrooms in this study were found to own the similar residues, which implicated that the transsulfuration reaction within ergothioneine biosynthesis pathway in edible mushrooms was different from the reaction in N. crassa. Regarding the substrate hercynylcysteine sulfoxide binding site, FvEgt2 showed different residues from others, which could be the reason why the ergothioneine biosynthesis pathway of F. velutipes requires two copies of Egt2 [16]. For the binding of cofactor PLP, CmEgt2 contains the same 4 PLP binding sites as NcEgt2, but the other 7 edible mushrooms exist one or two different sites. Furthermore, since GfEgt2 showed a lower identity with NcEgt2 and possessed a different active site with others, it implies that they might share different protein structure and chemical reaction mechanism.

In particular, the putative ergothioneine biosynthesis genes CmEgt1 (accession number CCM_07351) and CmEgt2 (accession number CCM_01645) extracted from C. militaris showed the highest identity score (56.43% and 48.78%) with NcEgt1 and NcEgt2. In addition, these genes have similar characteristics of sequence length and exon location (Fig. 1D). All these results implied the ergothioneine synthetases CmEgt1 and CmEgt2 were shared similar catalytic activity and substrate specific with NcEgt1 and NcEgt2, respectively.

Function determination of ergothioneine synthetases CmEgt1 and CmEgt2

The strategy of heterologous expression in E. coli with ergothioneine synthetases in C. militaris

Fig. 2
figure 2

Vector diagram of pET-CmEgt1-CmEgt2 (A) and pET-CmE1b2 (D). The SDS-PAGE analysis of time course experiment of BL21-pET-CmEgt1_Egt2 (B; CmEgt1, 99.4 kDa; CmEgt2_His-tag, 53.3 kDa) and BL21-pET-CmE1b2 (E; CmE1b_His tag, 59.5 kDa; CmEgt2_His tag, 53.8 kDa) (Pre, the sample of fermentative cells before induction; 3–12 h, the sample of fermentative cells during the fermentation after induction; triangle and arrow, pointed out the putative CmE1b and CmEgt2 respectively). Predicted protein structure (C) of CmEgt1 by I-TASSER (the cut site of truncated protein CmE1b was pointed out)

To verify the ergothioneine synthesis function of CmEgt1 and CmEgt2, we cloned their coding genes from the transcript of C. militaris [25], and inserted the genes into E. coli expression vector to build pET-CmEgt1-CmEgt2 (Fig. 2A). For monitoring the expression of these proteins, the vector was transformed into E. coli BL21(DE3) to perform time course experiment. As result showed (Fig. 2B), the CmEgt2 were expressed after the activation of the promoter by inducer. However, we did not observe the overexpression of CmEgt1. We subsequently performed multiple trials, including changing different expression vectors (including pET-22b from Novagen, Merck Millipore, USA and pCold-I from Takara, Japan) and performing optimization of fermentation, but all of them were failed to optimize the expression of CmEgt1. Even though the SDS-PAGE analysis did not show the soluble expression of CmEgt1, we further performed fermentation of BL21-pET-CmEgt1-CmEgt2 to detect the ergothioneine. Unfortunately, we could not detect any ergothioneine in fermentative medium or sonicated cell pellet via HPLC assay (Additional file 1: Fig. S1).

Previous study indicated the NcEgt1 could be expressed in E. coli [32], but CmEgt1 was failed to express in this study. It implicated they may share different protein structure, so we performed protein structure prediction of CmEgt1 with I-Tasser to figure out the problem. As the protein structure showed (Fig. 2C), the CmEgt1 was assembled with three distinguish parts, which contained a putative unstructured signal peptide, a N-terminal S-adenosylmethionine (SAM)-dependent methyltransferase domain and EgtB (a nonheme iron-depended sulfoxide synthase that plays role in ergothioneine biosynthesis in bacterial) [34]-like domain. The signal peptides in fungal proteins usually affect the protein expression in E. coli. But we have tried to fuse a pelB leader sequence (the inherent peptide of pET-22b, performed function of secretion) to the N-terminal of CmEgt1, the pelB leader could not help to increase the expression of soluble recombined CmEgt1. It implied that the signal peptide of CmEgt1 may not be the only reason of expressing obstacle of CmEgt1.

Previous study [32] also reported that the N-terminal methyltransferase part of NcEgt1 showed a lower catalytic rate than the other part (NcEgt1 could catalyze the substrate histidine or TMH to synthesize Cys-TMH but the reaction rate of using TMH as substrate to perform synthesis was 100-fold fast than using the histidine). As the purpose of this research was constructing a C. militaris strain with higher production of ergothioneine, we would rather not to intensive study the function of the N-terminal of CmEgt1 because of its putative low enzymatic efficiency. Moreover, the N-terminal domain could be replaced by a more efficiency isoenzyme such as EgtD (ergothioneine biosynthetic methyltransferase in bacteria).

So, we directly cut off the signal peptide alongside with N-terminal domain at the middle of dividing peptide between the SAM-dependent methyltransferase domain and EgtB-like domain (Fig. 2C, Additional file 1: Table S1) to design the truncated protein CmE1b. Next, the CmE1b, together with CmEgt2, were inserted into the expression vector to build pET-CmE1b2 (Fig. 2D). To increase the production of proteins, the genes of CmE1b and CmEgt2 were driven by two separate T7 promoter respectively. This vector was therefore transformed into E. coli BL21(DE3) to perform enzyme co-expression. The time course experiment (Fig. 2E) obviously indicated the overexpression of CmE1b and CmEgt2 and proved our speculation that the signal peptide and/or N-terminal methyltransferase domain could impede the expression of CmEgt1 in E. coli.

The in vitro reaction verified the ergothioneine biosynthesis function of CmE1b and CmEgt2

Fig. 3
figure 3

A Vector diagram of pET-CmE1b2; B SDS-PAGE analysis of protein expression of E. coli BL21(DE3), BL21-pBAD-EgtD, BL21-pET-CmE1b2 (E1B2) and BL21-EgtD1b2 (E1B2D). (The supernatant (SP) and cell pellet (CP) of samples after the process of cell sonication and centrifuge; CmE1b, 59.5 kDa; CmEgt2, 53.8 kDa; EgtD, 36.4 kDa; putative proteins were pointed out by arrow); C The HPLC analysis (partly chromatogram of the retention time from 36.5 to 58.5 min, all the sample performed at the same condition and showed in same intensity) of in vivo reaction by supernatant from the fermented E. coli samples; D The putative ergothioneine synthesis pathway of BL21-EgtD1B2

The cut off of the methyltransferase domain of CmE1b suggested the pathway could not perform the conversion of histidine to TMH, so we built a pBAD-EgtD vector (Fig. 3A), which contains the sequence of EgtD with E. coli codon optimization, to endow E. coli with the ability to synthesize hercynine as reported before [14]. The pBAD-EgtD, coupling with pET-CmE1b2, were transformed into E. coli BL21(DE3) to build BL21-EgtD1B2. The fermentation of this strain could simultaneously obtain the proteins of CmE1b, CmEgt2, and EgtD (Fig. 3B). To rapid test the function of CmE1b and CmEgt2, we set up the in vitro reaction and directly processed the fermentative strains with sonication to obtain enzymes for whole cells reaction. As the HPLC chromatogram showed (Fig. 3C), only the samples extracted from BL21-EgtD1b2 could catalyze the substrates (contained histidine, cysteine, and SAM) to synthesize ergothioneine.

In this part, we combined the known function EgtD, which could catalysis histidine to TMH, with putative ergothioneine synthetases CmE1B and CmEgt2 to establish an in vitro ergothioneine synthesis pathway (Fig. 3D). Though the function of N-terminal of CmEgt1 was not verified by overexpression or deletion experiment, its high identical with N-terminal of NcEgt1 implied they both share the similar function of catalyzing histidine to TMH.

Construct the recombined C. militaris strains with high productive of ergothioneine

Previous section indicated the assembling of EgtD, CmE1b and CmEgt2 could perform in vitro ergothioneine synthesis, but the bio-function (ergothioneine biosynthesis and degenerative protection) of this new pathway in C. militaris is unclear. The construction of C. militaris strains, which were carrying the overexpression of ergothioneine synthetases, is therefore necessary for biofunction verification. So, the vectors of p390-CmEgt2-E1B and p390-CmEgt2-E1b_EgtD (Fig. 4A) were constructed and transformed into two kinds of C. militaris (contains the normal wild type CM15 and fruiting body degenerated mutant strain CMdf, Fig. 4B) via ATMT (Agrobacterium tumefaciens-mediated transformation) method. After several rounds of unselective subculture with the transformants, the colonies were inoculated on the resistance screening plates (Additional file 1: Fig. S2A). The remaining colonies were therefore selected based on the performance of grown rate and pigment accumulation (Additional file 1: Fig. S2B). The qualified recombined strain was further verified by PCR (Fig. S2C) and picked to perform fermentation (Additional file 1: Fig. S2DEF). Finally, we obtained four recombined C. militaris strains in total, and they were named 15-E1b2, Df-E1b2, 15-E1bD2 and Df-E1bD2.

These recombined strains were inoculated in PSB medium (Fig. 4C, Additional file 1: Fig. S2G), which were optimized to increase the ergothioneine production of edible mushrooms with sucrose as carbon source, to perform flask fermentation of ergothioneine. After 10 days cultivating, the mycelium of CMdf and CM15 produced 511.08 ± 39.69 and 867.72 ± 31.88 mg/kg DW of ergothioneine, respectively. The production were corresponding or higher than the strains (contains S. cerevisiae [18] and A. oryzae [21]) with heterologous expression of NcEgt1 and Egt2s. It indicated the potential of building C. militaris as an ergothioneine high productive host. In addition, the overexpression of CmE1B, CmEgt2, and EgtD in recombined C. militaris strains (Df-E1b2D or 15-E1b2D) could perform a 1.36 to 1.54 times (than CMdf or CM15) increase of ergothioneine production in the fermentation with PSB media (Fig. 4C). Consider the precursor shortage always be the main problem of metabolite production, we supplemented the substrates (histidine and SAM) in PSB media before the fermentation. As the result showed (Fig. 4D and E), the production of ergothioneine were all increased. The highest amount was up to 2485.32 ± 50.72 mg/kg DW, which were produced by 15-E1b2D with PSB-HSs medium (contained 0.3 mM histidine and SAM in PSB). It suggested that the bacteria source EgtD, CmE1B and CmEgt2 could be over-expressed and constructed a high yield ergothioneine synthesis pathway in C. militaris.

Besides, compared the ergothioneine yield of strains in PSB-HSs and PSB-HSl (contained 3 mM histidine and SAM in PSB), they showed that the excessive substrates may decrease the yield of ergothioneine, which suggested that the production of ergothioneine could be further optimized by starting from the adjustment of the concentration of additives.

Meanwhile, the ergothioneine production of Df-E1b2 and 15-E1b2 showed obviously decreased (compared to CMdf or CM15). It suggested the biosynthesis of ergothioneine in C. militaris was not limited by CmE1b and CmEg2. The overexpression of them resulted in a decreasing of ergothioneine implied they may join in other unknown reactions that compete with ergothioneine precursor. In other words, it means that the synthesis of the precursor substrate TMH was the rate-limiting step to synthesize ergothioneine in C. militaris.

Fig. 4
figure 4

A Vector diagram of p390-CmEgt2-E1B and p390-CmEgt2-E1b_EgtD; B Mature fruiting body of C. militaris CMdf and CM15; The ergothioneine concentration of C. militaris fruiting body degenerated mutant CMdf (spot and yellow) and wild type CM15 (slash and orange) and synthetases overexpressed strains (Df-E1b2/15-E1b2: C. militaris CMdf/CM15 with overexpression of CmE1b and CmEgt2; Df-E1bD2/15-E1b2D: C. militaris CMdf/CM15 with overexpression of CmE1b, CmEgt2 and EgtD.) in PSB (C), PSB-HSl (D) and PSB-HSs (E) media. Statistical analyses were performed using t-tests (**p < 0.001, ***p < 0.0001)

In this section, we constructed two expression boxes which contained the coding genes for ergothioneine synthesis (contained CmE1B, CmEgt2, and EgtD), and inserted them into the genome of C. militaris CMdf/15 for the purpose of breeding the ergothioneine high productive recombined strains. We obtained up to 2533.28 mg/kg DW of ergothioneine in PSB-HSs flask fermentation of C. militaris. So far, the fruiting body of edible mushrooms is the major diet source of ergothioneine for human being, because it usually contains the higher concentration of metabolites than the mycelium. But the fruiting body always faces a serious degeneration problem, which restricts the development of industrial bioactive ingredient extraction and medicinal mushroom food. The engineering overexpression strategy in this study revealed that the molecule breeding of C. militaris could engineer mycelium as a suitable platform to promote the yield of high-value ingredients.

Endogenous synthesized ergothioneine highly promote the cordycepin production of degenerated C. militaris

The continuous cultivation of C. militaris always along with high frequency degeneration, which represented as fruiting bodies distortion and secondary metabolites reduction. Previous research indicated that the ROS accumulation was one of the main factors to induce degeneration in C. militaris [28, 29]. Although the degenerate phenotype of C. militaris is irreversible, we speculated that the relief of ROS accumulation would recover the property of producing the main metabolites. Because the ergothioneine was reported to prevent the ROS over-accumulated in fungi cells [35], we investigated the differences of morphology and metabolites yield among wild type CM15, degenerative mutant CMdf and their recombined strains with ergothioneine synthetases (15-E1b2, Df-E1b2, 15-E1bD2 and Df-E1bD2) in C. militaris.

Fig. 5
figure 5

The cordycepin concentration (A) and the total content (B) of C. militaris wild type and its derived in AM medium after 30 days static submerged fermentation. Statistical analyses were performed using t-tests (*p < 0.01, **p < 0.001, ***p < 0.0001)

As Fig. 5A showed, the wild type C. militaris CM15 could produce 0.89 g/L of cordycepin in AM medium, while CMdf only produced 0.32 g/L of cordycepin. The lower yield of cordycepin of CMdf showed that the degeneration of CM15 will change the shape of the fruiting body, and perform a significant reduction in metabolite production, which consisted of the pattern of the C. militaris degeneration [28].

In addition, the cordycepin content (Fig. 5B) of 15-E1b2 showed that overexpression of CmE1B and CmEgt2 in healthy wild type (compared to CM15) could barely affect the cordycepin yield. But the concentration of cordycepin in Df-E1b2 had a double amount increasing (compared to CMdf). It indicated these two genes were not directly relevant to the production of cordycepin, but they were roles of joining the metabolic regulation of degenerated strains.

The results (Fig. 5) showed that the overexpression of CmE1B, CmEgt2 and EgtD obviously improved the cordycepin production, while the increasing patterns were different between wild type CM15 and shape mutant type CMdf. As Fig. 5A showed, the synthetases overexpression promoted a triple increase of cordycepin content in fermented broth (compared with samples of CMdf and Df-E1b2D), while the concentration in the mycelium barely changed. But the cordycepin concentration of 15-E1b2D was increased 1.29-time (in broth) and 2.34-time (in mycelium) than CM15. As result, the total cordycepin amount of 15-E1b2D (Fig. 5B), which reached 0.2 g per grass jar (contained 100 ml of AM medium), was double than its original strain CM15. However, the fermentation of Df-E1b2D performed just equal amount to that of in CM15. It showed that the cordycepin yield of the shape degenerated strain could be raised, but the in vivo cordycepin production will face increasing bottleneck. Regarding to the antioxidant ability of ergothioneine, it also implied that the change of ROS environment might be one of the reasons for the decrease in cordycepin yield in the shape degenerative strain. In summary, the increasing content of cordycepin indicated the ergothioneine synthetases played roles in the excretion or production of cordycepin.

Previous report [29] showed that the generation of C. militaris fruiting body could be recovered by introducing the antioxidase Gpx in genome of degenerative strain. It suggested that the antioxidant ability of C. militaris might also be important to fruiting body germination. We therefore performed fruiting body farming with wild types and mutants to further study the relationship between ergothioneine and fruiting body germination. The experiment was different of previous study [29], the degenerated strain of this study have the ability of generating fruiting body. We aimed to test the ability of recover the fruiting body shape from dwarf to long. Unfortunately, the morphology of the fruiting body between recombined strains and its original wild type did not show difference. Though the ergothioneine content of CMdf was increased, the ability of forming healthy slender fruiting body was failed to recover. It indicated the fruiting body mutant degeneration in C. militaris might be irreversible, but the amount of metabolites in mutant strains could be increased back to normal level.


We excavated the ergothioneine synthetases Egt1 and Egt2 of 8 micro fungi by bioinformatic homologous alignment and verified the ergothioneine synthesized function of CmE1B, CmEgt2, and EgtD in E. coli. By constructing this new pathway in C. militaris wild type, we increased the yield of ergothioneine and cordycepin by up to 2.5 g/kg DW and 2 g/L, respectively. Though the overexpression of these synthetases in degenerated C. militaris strain failed to recover the morphology of fruiting body back to normal, the improving ingredients concentration implied the potential of developing C. militaris as chassis cells for producing high amount of ergothioneine, and breeding a functional mushroom food to fulfill market needs.


DNA manipulation and strains construction

Vectors and strains used in this study were listed in Table 2. The coding sequence of putative ergothioneine synthetases (Additional file 1: Table S1) in C. militaris were amplified from total RNA of C. militaris CM15 [25] by Phanta Max Super-Fidelity DNA polymerase (Vazyme Biotech, Nanjing, China) with default conditions. All the vectors used in this study were constructed by overlap-PCR and ClonExpress II One Step Cloning Kit (Vazyme Biotech, Nanjing, China). Primers and detail constructed methods were list in Additional file 1: Table S2. The backbone vector used for the construction of synthetases overexpression in C. militaris was built before [26]. It contains two expression boxes which were under controlled by strong promoters PtrpC and Pcmgpd respectively. Briefly, we assembled two vectors, the first vector called p390-CmEgt2-E1B, which contained the CmEgt2, controlled by a medium strength PtrpC promoter, and the CmE1b, promoted by a strong strength Pcmgpd promoter. The other one called p390-CmEgt2-E1b_EgtD, which contained the coding gene of EgtD (Additional file 1: Table S3). In particularly, the EgtD coding sequence was located at the end of previous vector’s CmE1b open reading frame by a 2 A peptide linker [36]. The C. militaris strains were constructed with these two vectors by ATMT method as previous described [26]. The wild type of strain with inserting the expression boxes of CmE1b and CmEgt2 was named 15-E1b2, while the degenerated fruiting body mutant strain with the same insertion was named Df-E1b2. Similarly, the wild type CM15 with inserting the CmE1b_EgtD (coding genes of Egt2 and CmE1b were linked together) and CmEgt2 was named 15-E1bD2, while the degenerated mutant CMdf with the same insertion was named 10-E1bD2.

Table 2 Vectors and strains used in this study

Fermentation conditions

The protein overexpression of E. coli was performed in Lysogeny Broth (LB) and supplemented with corresponding antibiotic and cofactors. The precultures for overexpression cultivation were prepared at 37 ℃ with 200 rpm agitation for 16 h. The fermentation of E. coli was performed as follow. Seed culture were inoculated in LB media and cultivated at 37℃ until OD600 was 0.8, 0.2 mM IPTG or/and 0.5% l-arabinose was/were supplemented and continue cultivated at 20 ℃ for 12 h. The overexpression of proteins was verified by SDS-PAGE electrophoresis. In particularly, the samples of time course were collected by each hour while induced and cultivated for three hours.

C. militaris was inoculated in 100 ml PDB media and cultivated at 25℃ with shaking for 5 to 8 days as seed culture of ingredients fermentation. Ergothioneine fermentation of C. militaris was performed in potato sucrose broth (PSB, contains 200 g/L potato, 20 g/L sucrose, 4 g/L NH4Cl, 3 g/L KH2PO4 and 1.5 g/L MgSO4·7H2O; PSB-HSl, supplementary with 3 mM 0.3mM histidine and SAM in PSB; PSB-HSs, supplementary with 0.3mM histidine and SAM in PSB), which could induce C. militaris to produce higher yield of ergothioneine compared to several common macro-fungi cultivation media. 10% of seed culture were inoculated in PSB and cultivated at 25 ℃ for 4 days with shaking in dark, then cultivated with light for 6 days. Cordycepin fermentation were also begun with 10% of seed culture and performed in AM media at 25 ℃ for 30 days as previously described.

Detection and quantification of cordycepin and ergothioneine

After fermenting in AM media, the mycelia were dried off at 50 ℃ and ground to powder. Two hundred milligram of dry powder was weighted and immersed in 20% of methanol to perform ultrasonication. Then, the supernatant of media and the crushing fluid were collected to perform HPLC detection as previously method [25]. The detection of ergothioneine from C. militaris fermentative samples were performed as follow. The mycelia were collected and immersed in 70% of ethanol. After treating by ultrasonication, the supernatant was used for HPLC detection while the precipitation was dried off to obtain the wight of original mycelia. HPLC assay was performed by Ultimate HILIC Amphion II Column (4.6 × 250 mm, 5-Micron, Welch, Shanghai, China). The analysis conditions were as follow: mobile phase, 10% ultra-pure water and 80% acetonitrile (v/v); flow rate, 0.9 ml/min; detection wavelength, 259 nm; column temperature, 30 ℃. A standard ergothioneine curve was generated using 5–60 mg/L ergothioneine standard (Sigma-Aldrich, United States). The yield of cordycepin and ergothioneine were calculated using the detected peak area according to the standard curve. The cordycepin or ergothioneine concentrations of mycelia and/or fermentative media presented in the study were calculated by normalizing in the equal biomass.

Statistical and bioinformatic analysis

The experimental data were represented as mean ± standard error of mean of three replicates. For graphical representation and analysis Microsoft Excel were used. The yield of cordycepin or ergothioneine between two strains were compared based on biological triplicates and significant difference was identified with p-value < 0.01. The protein structure prediction was performed by uploading the amino acid sequence into I-TASSER On-line Server [37] ( with default setting. The sequence alignments were performed by using the online software BLAST41 in NCBI ( and MEGA (version 7.0,

Enzymatic reaction assay

Since the enzyme kinetics of ergothioneine synthetase EgtD [38, 39], EgtB [34], NcEgt1 [32] and NcEgt2 [33] had been studied. We performed enzymatic reaction with ultrasonicated whole cell to verify the ergothioneine biosynthesized function of CmEgt1 and CmEgt2 rather than using purified enzyme because of the overexpression of CmEgt1 in E. coli faced obstacles. After fermenting in LB for 12 h, the collected cells were diluted in Tris Buffer (100 mM Tris, 150 mM NaCl, pH8.0) to perform sonicating. After using high speed centrifuge to separate the cell pellet, all the supernatants were performed enzymatic reaction immediately. The reaction was set up at 25 ℃ for overnight with stirring and contained the following substrates: 1 mM histidine, 3 mM S-adenosylmethionine (SAM), 1 mM Cysteine, 1 mM dithiothreitol and 100 mM ascorbic acid. Last, the reacted samples were lyophilizate and filtrated before submitting to HPLC detection.

Availability of data and materials

All data generated or analysed during this study are included in this published article [and its supplementary information files]. The materials used and/or constructed during the current study are available from the corresponding author on reasonable request.



Dry weight


Reactive oxygen species


Pyridoxal 5-phosphatemonohydrate






Hercynylcysteine sulfoxide




  1. Cui JD. Biotechnological production and applications of Cordyceps militaris, a valued traditional Chinese medicine. Crit Rev Biotechnol. 2015;35:475–84.

    Article  Google Scholar 

  2. Phull A-R, Ahmed M, Park H-J. Cordyceps militaris as a bio functional food source: pharmacological potential, anti-inflammatory actions and related molecular mechanisms. Microorganisms. 2022;10:405.

    Article  CAS  Google Scholar 

  3. Zhang J, Wen C, Duan Y, Zhang H, Ma H. Advance in Cordyceps militaris (Linn) Link polysaccharides: isolation, structure, and bioactivities: a review. Int J Biol Macromol. 2019;132:906–14.

    Article  CAS  Google Scholar 

  4. Chen Y, Wu Y, Li S, Du S, Hao X, Zhang J, et al. Large-scale isolation and antitumor mechanism evaluation of compounds from the traditional Chinese medicine Cordyceps militaris. Eur J Med Chem. 2021;212:113142.

    Article  CAS  Google Scholar 

  5. Cheah IK, Halliwell B. Ergothioneine; antioxidant potential, physiological function and role in disease. Biochimica et Biophysica Acta. 2012;1822:784–93.

    Article  CAS  Google Scholar 

  6. Chan JSL, Barseghyan GS, Asatiani MD, Wasser SP. Chemical composition and medicinal value of fruiting bodies and submerged cultured mycelia of caterpillar medicinal fungus Cordyceps militaris CBS-132098 (Ascomycetes). Int J Med Mushrooms. 2015;17:649–59.

    Article  CAS  Google Scholar 

  7. Beelman RB, Kalaras MD, Phillips AT, Richie JP. Is ergothioneine a ‘longevity vitamin’ limited in the American diet? J Nutr Sci. 2020;9:e52.

    Article  CAS  Google Scholar 

  8. Ames BN. Prolonging healthy aging: longevity vitamins and proteins. Proc Natl Acad Sci USA. 2018;115:10836–44.

    Article  CAS  Google Scholar 

  9. Lam-Sidun D, Peters KM, Borradaile NM. Mushroom-derived medicine? Preclinical studies suggest potential benefits of ergothioneine for cardiometabolic health. Int J Mol Sci. 2021;22:3246.

    Article  CAS  Google Scholar 

  10. Cheah IK, Halliwell B. Ergothioneine, recent developments. Redox Biol. 2021;42:101868.

    Article  CAS  Google Scholar 

  11. Tsiantas K, Tsiaka T, Koutrotsios G, Siapi E, Zervakis GI, Kalogeropoulos N, et al. On the identification and quantification of ergothioneine and lovastatin in various mushroom species: assets and challenges of different analytical approaches. Molecules. 2021;26:1832.

    Article  CAS  Google Scholar 

  12. Kalaras MD, Richie JP, Calcagnotto A, Beelman RB. Mushrooms: a rich source of the antioxidants ergothioneine and glutathione. Food Chem. 2017;233:429–33.

    Article  CAS  Google Scholar 

  13. EFSA Panel on Dietetic Products, Nutrition and Allergies (NDA), Turck D, Bresson J-L, Burlingame B, Dean T, Fairweather-Tait S, et al. Statement on the safety of synthetic l-ergothioneine as a novel food—supplementary dietary exposure and safety assessment for infants and young children, pregnant and breastfeeding women. EFSA J. 2017;15:e05060.

    Google Scholar 

  14. Kamide T, Takusagawa S, Tanaka N, Ogasawara Y, Kawano Y, Ohtsu I, et al. High production of ergothioneine in Escherichia coli using the sulfoxide synthase from methylobacterium strains. J Agr Food Chem. 2020;68:6390–4.

    Article  Google Scholar 

  15. Tanaka N, Kawano Y, Satoh Y, Dairi T, Ohtsu I. Gram-scale fermentative production of ergothioneine driven by overproduction of cysteine in Escherichia coli. Sci Rep-Uk. 2019;9:1895.

    Article  Google Scholar 

  16. Yang X, Lin S, Lin J, Wang Y, Lin J-F, Guo L. The biosynthetic pathway of ergothioneine in culinary-medicinal winter mushroom, Flammulina velutipes (Agaricomycetes). Int J Med Mushrooms. 2020;22:171–81.

    Article  Google Scholar 

  17. Osawa R, Kamide T, Satoh Y, Kawano Y, Ohtsu I, Dairi T. Heterologous and high production of ergothioneine in Escherichia coli. J Agr Food Chem. 2018;66:1191–6.

    Article  CAS  Google Scholar 

  18. van der Hoek SA, Darbani B, Zugaj KE, Prabhala BK, Biron MB, Randelovic M, et al. Engineering the yeast Saccharomyces cerevisiae for the production of L-(+)-ergothioneine. Front Bioeng Biotech. 2019;7:262.

    Article  Google Scholar 

  19. Yu Y-H, Pan H-Y, Guo L-Q, Lin J-F, Liao H-L, Li H-Y. Successful biosynthesis of natural antioxidant ergothioneine in Saccharomyces cerevisiae required only two genes from Grifola frondosa. Microb Cell Fact. 2020;19:164.

    Article  CAS  Google Scholar 

  20. van der Hoek SA, Rusnák M, Wang G, Stanchev LD, de Fátima Alves L, Jessop-Fabre MM, et al. Engineering precursor supply for the high-level production of ergothioneine in Saccharomyces cerevisiae. Metab Eng. 2022;70:129–42.

    Article  CAS  Google Scholar 

  21. Takusagawa S, Satoh Y, Ohtsu I, Dairi T. Ergothioneine production with Aspergillus oryzae. Biosci Biotechnol Biochem. 2019;83:181–4.

    Article  Google Scholar 

  22. van der Hoek SA, Rusnák M, Jacobsen IH, Martínez JL, Kell DB, Borodina I. Engineering ergothioneine production in Yarrowia lipolytica. Febs Lett. 2021;596:1356–64.

    Article  Google Scholar 

  23. Pluskal T, Ueno M, Yanagida M. Genetic and metabolomic dissection of the ergothioneine and selenoneine biosynthetic pathway in the fission yeast, S. pombe, and construction of an overproduction system. PLoS ONE. 2014;9:e97774.

    Article  Google Scholar 

  24. Zheng P, Xia Y, Xiao G, Xiong C, Hu X, Zhang S, et al. Genome sequence of the insect pathogenic fungus Cordyceps militaris, a valued traditional Chinese medicine. Genome Biol. 2011;12:R116.

    Article  CAS  Google Scholar 

  25. Chen B-X, Wei T, Xue L-N, Zheng Q-W, Ye Z-W, Zou Y, et al. Transcriptome analysis reveals the flexibility of cordycepin network in Cordyceps militaris activated by l-Alanine addition. Front Microbiol. 2020;11:577.

    Article  Google Scholar 

  26. Chen B-X, Wei T, Ye Z-W, Yun F, Kang L-Z, Tang H-B, et al. Efficient CRISPR-Cas9 gene disruption system in edible-medicinal mushroom Cordyceps militaris. Front Microbiol Front. 2018;9:1157.

    Article  Google Scholar 

  27. Cohen N, Cohen J, Asatiani MD, Varshney VK, Yu H-T, Yang Y-C, et al. Chemical composition and nutritional and medicinal value of fruit bodies and submerged cultured mycelia of culinary-medicinal higher basidiomycetes mushrooms. Int J Med Mushrooms. 2014;16:273–91.

    Article  CAS  Google Scholar 

  28. Lou H, Lin J, Guo L, Wang X, Tian S, Liu C, et al. Advances in research on Cordyceps militaris degeneration. Appl Microbiol Biot. 2019;103:7835–41.

    Article  CAS  Google Scholar 

  29. Xiong C, Xia Y, Zheng P, Wang C. Increasing oxidative stress tolerance and subculturing stability of Cordyceps militaris by overexpression of a glutathione peroxidase gene. Appl Microbiol Biot. 2013;97:2009–15.

    Article  CAS  Google Scholar 

  30. Seebeck FP. In vitro reconstitution of Mycobacterial ergothioneine biosynthesis. J Am Chem Soc. 2010;132:6632–3.

    Article  CAS  Google Scholar 

  31. Stampfli AR, Blankenfeldt W, Seebeck FP. Structural basis of ergothioneine biosynthesis. Curr Opin Struc Biol. 2020;65:1–8.

    Article  CAS  Google Scholar 

  32. Hu W, Song H, Sae Her A, Bak DW, Naowarojna N, Elliott SJ, et al. Bioinformatic and biochemical characterizations of C-S bond formation and cleavage enzymes in the fungus Neurospora crassa ergothioneine biosynthetic pathway. Organic Letters. 2014;16:5382–5.

    Article  CAS  Google Scholar 

  33. Irani S, Naowarojna N, Tang Y, Kathuria KR, Wang S, Dhembi A, et al. Snapshots of C-S cleavage in Egt2 reveals substrate specificity and reaction mechanism. Cell Chem Biol. 2018;25:1–11.

    Article  Google Scholar 

  34. Stampfli AR, Goncharenko KV, Meury M, Dubey BN, Schirmer T, Seebeck FP. An alternative active site architecture for O2 activation in the ergothioneine biosynthetic EgtB from Chloracidobacterium thermophilum. J Am Chem Soc. 2019;141:5275–85.

    Article  CAS  Google Scholar 

  35. Saini V, Cumming BM, Guidry L, Lamprecht DA, Adamson JH, Reddy VP, et al. Ergothioneine maintains redox and bioenergetic homeostasis essential for drug susceptibility and virulence of Mycobacterium tuberculosis. Cell Rep. 2016;14:572–85.

    Article  CAS  Google Scholar 

  36. Tang X, Liu X, Tao G, Qin M, Yin G, Suo J, et al. “Self-cleaving” 2A peptide from porcine teschovirus-1 mediates cleavage of dual fluorescent proteins in transgenic Eimeria tenella. Vet Res. 2016;47:68.

    Article  Google Scholar 

  37. Yang J, Zhang Y. I-TASSER server: new development for protein structure and function predictions. Nucleic Acids Res. 2015;43:W174–81.

    Article  CAS  Google Scholar 

  38. Vit A, Misson L, Blankenfeldt W, Seebeck FP. Ergothioneine biosynthetic methyltransferase EgtD reveals the structural basis of aromatic amino acid betaine biosynthesis. Chembiochem. 2015;16:119–25.

    Article  CAS  Google Scholar 

  39. Misson L, Burn R, Vit A, Hildesheim J, Beliaeva MA, Blankenfeldt W, et al. Inhibition and regulation of the ergothioneine biosynthetic methyltransferase EgtD. Acs Chem Biol. 2018;13:1333–42.

    Article  CAS  Google Scholar 

Download references


This work was supported by the China Postdoctoral Science Foundation [Grant Number 2020M682731]; the National Natural Science Foundation of China [Grant Number 32072646, 31772373]; the Guangdong Basic and Applied Basic Research Foundation [Grant Number 2019A1515110609]; the Scientific Research Platform Project of Education Department of Guangdong Province [Grant Number 2020KTSCX019]; the project of Natural Science Foundation of Guangdong Province [Grant Number 2022A1515010057]; the project of Guangdong Province Agriculture Research Project & Agricultural Technique Promotion Project [Grant Number 2021KJ103; 2022KJ103].

Author information

Authors and Affiliations



B-XC designed the study, performed major experiments, analyzed data and was a major contributor in writing the manuscript. L-NX performed vector construction, transformation and metabolic detection. TW and Z-WY helped writing- review and editing the manuscript. X-HL performed strains cultivation and data curation. L-QG and J-FL supervised the study, acquired funding, and revised the paper. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Li-Qiong Guo or Jun-Fang Lin.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Table S1. Amino acid sequence and domains prediction of fungi ergothioneine synthetases. Table S2. Primers and detail constructed methods of vector construction in this study. Table S3. DNA sequence of synthesized EgtD. Fig. S1. The HPLC chromatogram of ergothioneine detection of the fermentative medium and cell pellet of BL21-pET-CmEgt1_Egt2. Fig. S2. The verification of qualified recombined C. militaris strains.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chen, BX., Xue, LN., Wei, T. et al. Enhancement of ergothioneine production by discovering and regulating its metabolic pathway in Cordyceps militaris. Microb Cell Fact 21, 169 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: