Skip to main content

Overlapping promoter library designed for rational heterogenous expression in Cordyceps militaris



Cordyceps militaris, a kind of edible and medicinal fungus widely accepted in East Asia, has attracted much attention as a potential cell factory for producing adenosine analogs. Despite the rapid development in gene editing techniques and genome modeling, the diversity of DNA elements in C. militaris was too short to achieve rational heterogeneous expression for metabolic engineering studies.


In this study, PtrpC, a kind of promoter with a relatively appropriate expression level and small size, was selected as a monomer for promoter library construction. Through in vitro BioBricks assembly, 9 overlapping PtrpC promoters with different copy numbers as well as reporter gene gfp were connected and subsequently integrated into the genome of C. militaris. Both the mRNA transcription level and the expression level of gene gfp gradually increased along with the copy number of the overlapping promoter NPtrpC and peaked at 7. In the meantime, no significant difference was found in either the biomass or morphological characteristic of engineered and wild-type strains.


This study firstly expanded the overlapping promoter strategy used in model microorganism in C. militaris. It was a proof-of-concept in fungi synthetic biology and provide a general method to pushed the boundary of promoter engineering in edible mushroom.


Cordyceps militaris is generally regarded as an edible and medicinal fungus in East Asia. It has recently been widely reported due to its diverse bio-active component (e.g., cordycepin [1], polysaccharide [2], carotenoid [3]), which have been proved to have anti-cancer, anti-tumor [4, 5], and anti-photoaging [6, 7] effects. It is the only species that can produce a large amount of cordycepin among more than 350 species of the Cordyceps genus. Artificial cultivation techniques of C. militaris have been developed for more than 40 years and most parameters (e.g., culture medium composition, culture conditions) have been fully optimized [8,9,10]. By 2015, the output of large-scale industrialized planting of C. militaris in China has reached 74,000 tons [11].

Diverse synthetic biology systems have been recently developed in C. militaris. In 2018, gene editing techniques based on homologous recombination [12] and CRISPR [13] were established in C. militaris, which made it possible to reconstruct its cordycepin metabolic network. What’s more, the genome metabolism model iNR1329 allowed genome-wide in silico analysis of gene of interest (GOI) in C. militaris [14]. For another, the study of setting up C. militaris as microbial cell factory boomed in recent years with the rapid development in synthetic biology. Zou et al.used C. militaris as chassis cells to transform them into a cell factory with high-yield anticancer drug pentostatin [15]. Since the industrial cultivation techniques for various mushroom was mature, there is growing interest in such attempts to reconstruct mushroom species to overproduce their valuable bio-active compound, and on a more and more complexity level.

Promoter engineering was the most popular and convenient strategy to regulate the expression level of GOI [16, 17]. And constitutive promoter present stable expression level, which was essential to accurate regulation. However, the poor DNA element diversity has limited the present development of fungi synthetic biology to achieve a fixed ratio of multiple gene expression levels. So far, there are only four kinds of constitutive promoters (e.g., Pgpd [18], PtrpC [19], Pcmlsm3 [13], PCaMV [20]) reported in C. militaris. With limited available promoter on hand, the overlapping constitutive promoter was a practicable strategy widely applied in edible mushroom. And BioBricks method was capable to construct gene cluster repeatedly and sequentially by reusing isocaudarners with similar restriction sites [21]. Constitutive promoter PtrpC [19] regulated the glyceraldehyde 3-phosphate dehydrogenase gene in C. militaris was chosen to be overlapped due to its small size and medium expression strength (Additional file 1: Table S1). 1 to 9 PtrpC were overlapped by BioBricks method, ligated to reporter gene gfp and each module was integrated into the genome of C. militaris strain CM10 by Agrobacterium-mediated transformation (ATMT) respectively. Both mRNA and fluorescence intensity increased with the copy number of promoters at the beginning and presented a normal distribution as a whole. This study proved the overlapping promoter strategy worked well in fungi and provided a new path to achieve rational balance among multiple GOI in mushrooms.


Construction of pMD19T (simple)-1-9PtrpC-gfp-Tnos

The size of promoter PtrpC was only 369 bp, which was much smaller than other constitutive promoters (Additional file 1: Table S1) discovered in C. militaris. Therefore, it was chosen to construct an overlapping promoter library with 1–9 copies to achieve quantitative expression of multiple target genes.

After 4 rounds of BioBricks assembly, vectors pMD19T(simple)-1-9PtrpC-gfp-Tnos were constructed (Fig. 1). When these 9 promoters were digested by EcoRI and BglII, the remaining vector backbones should be 2705 bp and the operon with overlapping promoters should be 1362, 1737, 2112, 2487, 2862, 3237, 3612, 3987, and 4362 bp respectively. The experimental result was consistent with the theoretical result (Fig. 2).

Fig. 1
figure 1

Flow diagram of pMD19T(simple)-1-9PtrpC-gfp-Tnos vector construction

Fig. 2
figure 2

Restriction enzyme identification of pMD19T(simple)-1-9PtrpC-gfp-Tnos. (M: DNA marker, 1–9 were: pMD19T(simple)-1-9PtrpC-gfp-Tnos, respectively)

Construction of pCambia0390-blpR-1-9PtrpC-gfp

The cassettes 1-9PtrpC-gfp-Tnos cut from vectors pMD19T(simple)-1-9PtrpC-gfp-Tnos were inserted into shuttle vector pCambia0390 respectively. They were separated by EcoRI and BglII into an 8109 bp vector backbone and a 1362, 1737, 2112, 2487, 2862, 3237, 3612, 3987, and 4362 bp operon respectively. The experimental result was consistent with the theoretical result (Fig. 3).

Fig. 3
figure 3

Restriction enzyme digestion identification of pCambia0390-blpR-1-9PtrpC-gfp. Vector pCambia0390-blpR-1-9PtrpC-gfp were double digested (M: DNA marker; 1–9: pCambia0390-blpR-1-9PtrpC-gfp, respectively)

Colony PCR results showed that vectors pCambia0390-blpR-1-9PtrpC-gfp were successfully transferred into A. tumefaciens strain AGL-1 by electro-transformation. And recombinant C. militaris hyphae infected by A. tumefaciens AGL_ pCambia0390-blpR-1-9PtrpC-gfp were scraped and identified by PCR (Fig. 4). Since the robustness of fungi recombinant was relatively low, the obtained resistant transformed strains were continuously sub-cultured and further identified by PCR (Fig. 4). The stable strains were those with purposeful bands at the sizes of 1325, 1700, 2075, 2450, 2825, 3200, 3575, 3950, and 4325 bp respectively. The 3rd generation of stable recombinants was named as CmNtrpC and selected for subsequent experiments. Result showed that no significant difference in mycelia was found between wild-type and recombinant strains (Fig. 5a and b).

Fig. 4
figure 4

PCR identification of recombinant transformants of C. militaris. (M: DNA marker; 1–9: Cm1-9trpC, respectively)

Fig. 5
figure 5

Mycelia of C. militaris grown on plates. a, the wild-type strains; b, the recombinant strains

Stable expression of gene gfp in C. militaris

The transcription levels of gene gfp in CmNtrpC were detected by qPCR and the mycelial fluorescence intensities were investigated by laser confocal microscope. The ratio of transcription level of engineered strains to Cm1trpC increased with the copy number of overlapping promoters PtrpC and peaked at 7 (Fig. 6). Compared with strain Cm1trpC, the transcription level of strain Cm7trpC increased 63 folds. The transcription levels of strains Cm4-7trpC were relatively high while the others were significantly lower. This result indicated that the transcription strength can be enhanced by increasing the tandem repeats of the PtrpC promoter. However, the increase was in a normal distribution. And even the copy number increased from 1 to 7, the increase of fluorescence intensity was not linear.

Fig. 6
figure 6

qRT-PCR detection of gene gfp in transformant strains. The gene expression of Cm1trpC in control is normalized to 1. The data is the average of three experiments, and the error line shows the standard error. T-test was used for statistical analysis (p < 0.05).

The protein expression level was basically consistent with the result of transcription level. The fluorescence intensity of a single mycelium gradually increased with the copy number of overlapping promoters (Fig. 7 and Additional file 1: Table S3). It also peaked when the copy number reached 7. The fluorescence intensities of strains Cm4-7trpC were relatively higher while the others were too low to be captured by software ImageJ. The fluorescence intensity values of remaining strain Cm4-8trpC showed similar fold change. For example, the mRNA transcription level of Cm7trpC was 1.7-fold higher than Cm6trpC, while fluorescence intensity was 1.6-fold. Results of both mRNA level and fluorescence intensity indicated that the expression level of the target gene could be enhanced with the increase of overlapping promoter number. The increase was not limitless, and it would usually reach a boundary with the increase of promoter copy number.

Fig. 7
figure 7

Fluorescence intensities of GFP in transformant mycelia. Mycelia of the Cm1-9trpC transformants were shown under bright field image and fluorescence microscope (excitation, 395 to 440 nm, and emission, 470 nm). Scale bars represent 20 µm

The biomass of both wild-type strain and engineered strains CmNtrpC were measured. Compared with the wild-type strain CM10, the biomass of engineered strains slightly increased with the copy number of promoters (Fig. 8), indicating that the integration of gene gfp module with multiple overlapping promoters would not burden the growth of C. militaris. Furthermore, biomass data, growth rate, and biomass productivity, which were key indexes of cell growth, showed no significant difference (Additional file 1: Table S2). And since strain CM10 was degenerated and lost its capability to form fruiting body, it was necessary to compare the morphological characteristic of fruiting body when similar promoter library was integrated in an undegenerated strain.

Fig. 8
figure 8

Biomass of wild-type and recombinant strains of C. militaris. The data is the average of three experiments, and the error line shows the standard error. *It indicates that the biomass of recombinant strain is significantly higher when it was compared with that of wild-type strain, *p < 0.05, ** p < 0.01


In this study, modules with 9 overlapping promoters were constructed by the BioBricks method and integrated into the edible fungus C. militaris via ATMT. The result of transcription level detected by RT-PCR was well correlated with the result of fluorescence intensity detected by confocal microscope. The initiation intensity of reporter gene gfp showed a normal distribution while the initiation intensity reached the highest level with 7 monomer promoters. The pattern was consistent with the previous report in E. coli [16]. It decreased when the copy number of promoters passed a particular singularity. It might be may be due to the transcription interference caused by the connection between RNA polymerases and tandem promoters [22]. Yet, it takes further experimental methods, e.g., ChIP, to dig out the truth.

Furthermore, even though this study explored the pattern of overlapping PtrpC promoter, the standard of monomer promoter (e.g., minimum and maximum size, strength, DNA sequence characteristic) was still unclear. Previous study reported that promoter P37 with the size of only 55 bp was a practicable candidate for overlapping promoter engineering in E. coli [23]. Therefore, the minimum size of the repetitive promoter worked in C. militaris might be much smaller than the size of promoter PtrpC. Kanamasa et al. concatenated the promoter region III with 12 copies and used it to enhance the expression of β-Mannosidase in Aspergillus echinosporus [24]. What’s more, Zhang et al. overlapped 4 copies of cbhI inducible promoter to improve the expression level of the β-glucosidase gene (bgl1) in Trichoderma reesei) [25]. The size of promoter cbh1 was larger than 1291 bp. It implied that the maximum size of the monomer promoter could be much larger than the promoter we used in this study. And not only constitutive promoter but inducible promoter could be applied as a candidate in overlapping promoter strategy. For another, the carrying capacity of the general plasmid was limited and usually no more than 50 kb. The maximum size of the monomer promoter might be limited by the current molecular cloning method. Considering the convenience of molecular cloning and better binding between RNA polymerase and promoter sequence, a hypothesis for overlapping promoter strategy in filamentous fungi was provided in this study: a small promoter with medium strength would be a better candidate. Therefore, frequently used promoters in filamentous fungi (e.g., PoliC with only 212 bp, and PacuD with only 286 bp [26]) might be a practical choice.

Diverse metabolic engineering strategies in C. militaris, including CRISPR and split marker homologous recombination to edit genomic DNA as well as genome model iNR1329 to perform in silico metabolic analysis, were developed in recent years [12,13,14]. To make full use of these strategies and further develop fungus synthetic biology, it was essential and urgent to expand the promoter database and achieve accurate gene expression. And to rebalance the metabolic flux of different pathways, it was common to use promoters with different intensities to regulate and control the expression of multiple genes [27]. However, it was knotty to develop diverse promoters with different intensities in undeveloped mushroom species. Previous research generally mutated homologous promoters [28] or utilized heterogeneous promoters, which were either laborious or violate the principle of gene editing products [29]. This research aided the rational construction of promoter library with different intensities in undeveloped mushroom species and provided a one-step method for heterogeneous expression with diverse strength in C. militaris without overusing promoter sources, which saves more available choices for additional editing. It also proved that it was feasible to transform the synthetic biology strategy developed in model species to non-model fungus.

The engineered strains CmNtrpC were sub-cultured and tested by PCR for 3 generations. No significant degeneration was observed. With the addition of applying natural selection and culture preservation, the natural problem of C. militaris degeneration was soluble. Furthermore, even though the molecular mechanism of C. militaris degeneration remained unclear [30], the overlapping promoter strategy has provided one-step repetitive construction of engineered C. militaris strain with fix expression ratio of multiple target strains. Once the engineered strains were degenerated, they could be quickly reconstructed by re-invasion of engineered Agrobacterium strain AGL-1_pCAMBIA0390-NPtrpC-GOI to any wild-type C. militaris strain. The Agrobacterium strain rarely degenerated and the constructed vector was highly stable.

Compared with model species with clear genetic background and diverse promoter elements, such as E. coli and Saccharomyces cerevisiae, C. militaris was the Robinson Crusoe island with restricted choices. Even though promoter PtrpC had appropriate moderate intensity and relatively small size, it was still way too larger than the promoter chosen in the overlapping promoter research in E. coli [16, 23]. The combination of 9 PtrpC promoters was about 3.2 kb and larger than general gene. It made the library lack flexibility in vector construction and module integration. With more and more in silico tools developed for promoter prediction, it was still necessary to discover promoters with smaller sizes and appropriate starting strength in the future.

Materials and methods

Strains and cultivation conditions

E. coli strain DH5α (stored in our laboratory) was used in the construction and transformation of vectors. A. tumefaciens strain AGL-1(Weidi Bio, Shanghai, China) and the pCAMBIA0390 vector (Cambia, Queensland, Australia) were used for fungal transformation. C. militaris strain CM10 was purchased from Haixin Biological Co., LTD (Shandong, China) as the host for gene disruption. Luria- Bertani medium was used for the cultivation of E. coli strain DH5α and A. tumefaciens strain AGL-1. C. militaris strain CM10 was cultured in potato peptone dextrose agar (PDA: 200 g/L potatoes, 20 g/L glucose, 3 g/L KH2PO4, 1.5 g/L MgSO4, 20 g/L agar)) at 25℃. The materials for the strains and plasmids were shown in Table 1.

Table 1 Materials of strains and plasmids

Plasmid construction

Construction of pMD19T (simple)-1-9PtrpC-gfp-Tnos

The terminator Tnos was cloned from pCambia0390 by forward primer Tnos-F with EcoRI and BamHI restriction sites and reverse primer Tnos-R with BglII restriction site to build up BioBricks adaptor. And it was linked to vector pMD19T(Simple)(Takara, Beijing, China) by TA cloning and transformed into the competent cell of E. coli strain DH5α prepared by 0.1 M CaCl2. The recombinant was identified by colony PCR and DNA sequencing from transformants.

Gene gfp coding green fluorescent protein was amplified by PrimeSTAR Max (Takara, Beijing, China) with vector pPrstA-GFP as template and GFP-F/R as primer. EcoRI-BamHI cleavage site and protective base were added to forward primer GFP-F, and BglII cleavage site and protective base were added to reverse primer GFP-R. Gene gfp fragment and vector pMD19T(simple)-Tnos were digested by isocaudarner BamHI and BglII and further linked by T4 DNA ligase (Takara, Beijing, China) to generate a scar sequence (AGATCC). The linked product was transformed into strain DH5α and the transformant was identified.

Promoter PtrpC was amplified largely by the same method that performed to gene gfp. Vector pAg1-H3 was used as template and trpC-F/R was used as primer (Table. 2). What’s more, more copies of promoter PtrpC were inserted and vector library containing pMD19T(simple)-2–9-PtrpC-gfp-Tnos were constructed by the BioBricks method.

Table 2 Primers for PCR and qPCR

Construction of pCambia0390-blpR-1-9PtrpC-gfp-Tnos

To achieve Agrobacterium-mediated transformation, shuttle vector pCambia0390-blpR was chosen for the integration of the repetitive PtrpC promoter library. Digested by EcoRI and BglII, an 8109-bp linear vector backbone was purified by a DNA Gel Extraction Kit (Magen, Beijing, China). To further separate the operon with overlapping promoters, vectors pMD19T(simple)-1-3PtrpC-gfp-Tnos and pMD19T(simple)-7-9PtrpC-gfp-Tnos were digested by EcoRI and BglII, thereinto, ScaI was additionally added to digest backbone vector and obtain NPtrpC-gfp-Tnos module since the size of restriction fragments of pMD19T(simple)-4-6PtrpC-gfp-Tnos were similar, while pMD19T(simple)-4-6PtrpC-gfp-Tnos were digested by EcoRI, BglII, and ScaI. The operon with overlapping promoters with sizes of 1362, 1737, 2112, 2487, 2862, 3237, 3612, 3987, and 4362 bp was respectively purified. Ligation with overlapping promoters was transformed. Vectors pCambia0390-blpR-1-9PtrpC-gfp-Tnos were obtained by colony PCR and enzyme digestion.

Agrobacterium-mediated transformation

A. tumefaciens-mediated transformation(ATMT) was used as described in a previous study [13]. Shuttle vectors pCambia0390-blpR-1-9PtrpC-gfp-Tnos were transformed to the competent cells of A. tumefaciens strain AGL-1. Recombinant was re-inoculated in 5 mL LB medium (with kanamycin 50 mg/L and carbenicillin 50 mg/L) at 30 ℃ for 2–3 days to obtain seed solution. It was further inoculated into IM medium (1.45 g/L KH2PO4,2.05 g/L K2HPO4, 0.15 g/L NaCl, 0.5 g/L MgSO4·7H2O, 66 mg/L CaCl2·2H2O, 2.48 mg/L FeSO4·7H2O, 0.5 g/L (NH4)2SO4, 1.8 g/L Glucose, 5 mL/L Glycerin, pH 5.5, with 200 μM AS) with initial 0.15 OD600 value and kept in suspension culture at 30℃ till the OD600 value reached 0.8. A. tumefaciens suspension preheated at 25℃ and C. militaris conidia suspensions (1 × 104 conidia mL−1) were mixed in diverse certain proportions (1:100, 1:1, 100:1). Mixture was separated on IMA solid medium (1.45 g/L KH2PO4,2.05 g/L K2HPO4, 0.15 g/L NaCl, 0.5 g/L MgSO4·7H2O, 66 mg/L CaCl2·2H2O, 2.48 mg/L FeSO4·7H2O, 0.5 g/L (NH4)2SO4, 0.9 g/L Glucose, 5 mL/L Glycerin, 15 g/L Agar, pH 5.5, with 200 μM AS, pH = 5.5) pasted with glass paper (diameter = 90 mm) and cultured in the dark at 25℃ for 3 days. Glass paper was transferred to PDA medium with cephalosporin (300 mg/L) and Basta (400 mg/L) for 7 days when spores were germinated and white spots were emitted. Once was visible, a single mycelial mass was picked it by toothpick and transferred to the potato PDA with Basta (400 mg/L). When a single mycelium grows to 2 cm in diameter, cut 0.5 cm in diameter for PCR verification of crude mycelium samples. Verified transformants were re-cultivated till passed to the third generation. Mycelia were collected to extract genomic DNA and perform PCR analysis with primer pairs Tu3-F/GFP-R. Putative transformants were tested by quantitative real-time PCR and confocal laser scanning microscope.

qPCR verification of gene gfp transcription level initiated by overlapping promoters

Transformants with overlapping promoters were inoculated in liquid PDB media at 25 ℃ and 150 rpm for about 7 days. Mycelium was washed with ddH2O 3 times and ground with liquid nitrogen. 100 mg powder was separated to extract total RNA with E.Z.N.A. Fungal RNA Miniprep kit (OMEGA Bio-Tek Inc., GA, USA). 1 µg RNA was taken and reverse transcribed into cDNA by HiScript III-RT SuperMix (Vazyme, Nanjing, China). 50 ng of cDNA, 0.2 μL each of relevant primers (Table. 2) and SYBR qPCR Master Mix (Vazyme, Nanjing, China) were combined to build up the 20-µL qPCR system. ABI 7500 Real-Time PCR System (Thermo Fisher Scientific, MA, USA) and 18S rRNA gene as the internal control gene [31] were used for Quantitative Real-Time PCR. The cDNA of wild-type strain CM10 and engineered strains Cm1-9trpC were used as templates and 18S rRNA-F/R and GFP-F/R were used as primers for qRT-PCR amplification. Taking wild-type strain CM10 as the negative control and the gfp gene expression of C. militaris Cm1trpC as the control, the relative transcription level of gene gfp in other 8 strains were calculated respectively. The relative transcription of mRNA was computed by 2−ΔΔCT (livak) [32].

Observation of fluorescence C. militaris by confocal laser scanning microscope

Dry and sterilized cover glass (20 mm × 20 mm) was inserted into the PDA medium inoculated with corresponding transformant blocks at an angle of 45°. Transformant was subsequently cultured at 25℃ till mycelium grew beyond the insertion position. And the cover glass was taken out for observation in Carl Zeiss LSM800 light and fluorescence microscope (Carl Zeiss AG, Oberkochen, Germany). The mycelial fluorescence of engineered strains was observed and the fluorescence intensity of wild-type strain CM10 was regarded as the negative control to remove the influence of its own fluorescence background.

Availability of data and materials

All data generated or analyzed during this study are included in this published article and its additional files.


  1. Cunningham KG, Manson W, Spring FS, Hutchinson SA. Cordycepin, a metabolic product isolated from cultures of Cordyceps militaris (Linn.) Link. Nature. 1950;166:949.

    Article  CAS  Google Scholar 

  2. Huang ZF, Zhang ML, Zhang S, Wang YH, Jiang XW. Structural characterization of polysaccharides from Cordyceps militaris and their hypolipidemic effects in high fat diet fed mice. RSC Adv. 2018;8:41012–22.

    Article  CAS  Google Scholar 

  3. Zhao Y, Li SL, Chen HY, Zou Y, Zheng QW, Guo LQ, Wu GH, Lu J, Lin JF, Ye ZW. Enhancement of carotenoid production and its regulation in edible mushroom Cordyceps militaris by abiotic stresses. Enzyme Microb Technol. 2021;148: 109808.

    Article  CAS  Google Scholar 

  4. Lee HH, Lee S, Lee K, Shin YS, Kang H, Cho H. Anti-cancer effect of Cordyceps militaris in human colorectal carcinoma RKO cells via cell cycle arrest and mitochondrial apoptosis. DARU J Pharm Sci. 2015;23:35.

    Article  Google Scholar 

  5. Dong J, Li Y, Xiao H, Luo D, Zhang S, Zhu C, Jiang M, Cui M, Lu L, Fan S. Cordycepin sensitizes breast cancer cells toward irradiation through elevating ROS production involving Nrf2. Toxicol Appl Pharmacol. 2019;364:12–21.

    Article  CAS  Google Scholar 

  6. Lee YR, Noh EM, Jeong EY, Yun SK, Jeong YJ, Kim JH, Kwon KB, Kim BS, Lee SH, Park CS, et al. Cordycepin inhibits UVB-induced matrix metalloproteinase expression by suppressing the NF-kappaB pathway in human dermal fibroblasts. Exp Mol Med. 2009;41:548–54.

    Article  CAS  Google Scholar 

  7. Kunhorm P, Chaicharoenaudomrung N, Noisa P. Enrichment of cordycepin for cosmeceutical applications: culture systems and strategies. Appl Microbiol Biotechnol. 2019;103:1681–91.

    Article  CAS  Google Scholar 

  8. Mao XB, Eksriwong T, Chauvatcharin S, Zhong JJ. Optimization of carbon source and carbon/nitrogen ratio for cordycepin production by submerged cultivation of medicinal mushroom Cordyceps militaris. Process Biochem. 2005;40:1667–72.

    Article  CAS  Google Scholar 

  9. Masuda M, Urabe E, Honda H, Sakurai A, Sakakibara M. Enhanced production of cordycepin by surface culture using the medicinal mushroom Cordyceps militaris. Enzyme Microb Technol. 2007;40:1199–205.

    Article  CAS  Google Scholar 

  10. Chen YS, Liu BL, Chang YN. Effects of light and heavy metals on Cordyceps militaris fruit body growth in rice grain-based cultivation. Korean J Chem Eng. 2011;28:875–9.

    Article  Google Scholar 

  11. Bao DP, Liu MX, Wang Y, Yang RH, Mao WJ, He HQ. The compatible crosses assisted by the molecular marker of the mating genes of Cordyceps militaris. J Nuclear Agric Sci. 2017;31:2113–20.

    Google Scholar 

  12. Lou HW, Ye ZW, Yun F, Lin JF, Guo LQ, Chen BX, Mu Z. Targeted gene deletion in Cordyceps militaris using the split-marker approach. Mol Biotechnol. 2018;60:380–5.

    Article  CAS  Google Scholar 

  13. Chen BX, Wei T, Ye ZW, Yun F, Kang LZ, Tang HB, Guo LQ, Lin JF. Efficient CRISPR-Cas9 gene disruption system in edible-medicinal mushroom Cordyceps militaris. Front Microbiol. 2018;9:1157.

    Article  Google Scholar 

  14. Raethong N, Wang H, Nielsen J, Vongsangnak W. Optimizing cultivation of Cordyceps militaris for fast growth and cordycepin overproduction using rational design of synthetic media. Comput Struct Biotechnol J. 2020;18:1–8.

    Article  CAS  Google Scholar 

  15. Zou G, Li B, Wang Y, Yin X, Gong M, Shang JJ, Wei YJ, Li XL, Bao DP. Efficient conversion of spent mushroom substrate into a high value-added anticancer drug pentostatin with engineered Cordyceps militaris. Green Chem. 2021;23:10030–8.

    Article  CAS  Google Scholar 

  16. Li M, Wang J, Geng Y, Li Y, Wang Q, Liang Q, Qi QS. A strategy of gene overexpression based on tandem repetitive promoters in Escherichia coli. Microb Cell Fact. 2012;11:19.

    Article  CAS  Google Scholar 

  17. Lu Q, Bu YF, Liu JZ. Metabolic engineering of Escherichia coli for producing astaxanthin as the predominant carotenoid. Mar Drugs. 2017;15:296.

    Article  Google Scholar 

  18. Gong Z, Ying S, Lei H, Lin J, Zhou X. Cloning and analysis of glyceraldehyde-3-phosphate dehydrogenase gene from Cordyceps militaris. Afr J Agric Res. 2009;4:402–8.

    Google Scholar 

  19. Zheng Z, Huang C, Li C, Xie C, Han R. Agrobacterium tumefaciens-mediated transformation as a tool for insertional mutagenesis in medicinal fungus Cordyceps militaris. Fungal Biol. 2011;115:274.

    Google Scholar 

  20. Zhang H, Wang YX, Tong XX, Yokoyama W, Cao J, Wang F, Peng C, Guo JL. Overexpression of ribonucleotide reductase small subunit, RNRM, increases cordycepin biosynthesis in transformed Cordyceps militaris. Chin J Nat Med. 2020;18:393–400.

    CAS  PubMed  Google Scholar 

  21. Knight T. Idempotent vector design for standard assembly of BioBricks. MIT Synthetic Biology Working Group. 2005.

  22. Sneppen K, Dodd IB, Shearwin KE, Palmer AC, Schubert RA, Callen BP, Egan JB. A mathematical model for transcriptional interference by RNA polymerase traffic in Escherichia coli. J Mol Biol. 2005;346:399–409.

    Article  CAS  Google Scholar 

  23. Wei T, Cheng BY, Liu JZ. Genome engineering Escherichia coli for L-DOPA overproduction from glucose. Sci Rep. 2016;6:30080.

    Article  CAS  Google Scholar 

  24. Kanamasa S, Kawaguchi T, Takada G, Kajiwara S, Sumitani J, Arai M. Development of an efficient production method for beta-mannosidase by the creation of an overexpression system in Aspergillus aculeatus. Lett Appl Microbiol. 2007;45:142–7.

    Article  CAS  Google Scholar 

  25. Liu T, Wang TH, Li X, Liu X. Improved heterologous gene expression in Trichoderma reesei by cellobiohydrolase I gene (cbh1) promoter optimization. Acta Biochim Biophys Sin (Shanghai). 2008;40:158–65.

    Article  CAS  Google Scholar 

  26. De Lucas JR, Gregory S, Turner G. Analysis of the regulation of the Aspergillus nidulans acuD gene, encoding isocitrate lyase, by construction of a hybrid promoter. Mol Gen Genet MGG. 1994;243:654–9.

    Article  Google Scholar 

  27. Cox RS 3rd, Surette MG, Elowitz MB. Programming gene expression with combinatorial promoters. Mol Syst Biol. 2007;3:145.

    Article  Google Scholar 

  28. Blazeck J, Alper H. Systems metabolic engineering: genome-scale models and beyond. Biotechnol J. 2010;5:647–59.

    Article  CAS  Google Scholar 

  29. Zhang Y, Liu H, Liu Y, Huo K, Wang S, Liu R, Yang C. A promoter engineering-based strategy enhances polyhydroxyalkanoate production in Pseudomonas putida KT2440. Int J Biol Macromol. 2021;191:608–17.

    Article  CAS  Google Scholar 

  30. Lou HW, Lin JF, Guo LQ, Wang XW, Tian SQ, Liu CX, Zhao Y, Zhao RY. Advances in research on Cordyceps militaris degeneration. Appl Microbiol Biotechnol. 2019;103:7835–41.

    Article  CAS  Google Scholar 

  31. Fan DD, Wang W, Zhong JJ. Enhancement of cordycepin production in submerged cultures of Cordyceps militaris by addition of ferrous sulfate. Biochem Eng J. 2012;60:30–5.

    Article  CAS  Google Scholar 

  32. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods. 2001;25:402–8.

    Article  CAS  Google Scholar 

Download references


Biorender ( was used for the making of Fig. 1.


We are grateful to the Guangdong Basic and Applied Basic Research Foundation (No. 2019A1515110609), Scientific Research Platform Project of Education Department of Guangdong Province (No. 2020KTSCX019) for their financial support.

Author information

Authors and Affiliations



ML wrote the manuscript text and carried out most experiments. JZ, YZ, TZ, JM, and ZW are responsible for performing some experiments. AW, ACG, EGA carried out part of confocal laser microscope experiment. JL provided the resources. TW conceived of the study, participated in its design, coordination and helped to draft the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Tao Wei.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no financial or commercial conflict of interest.

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.

Basic parameters of the promoter commonly used in Cordyceps militaris. Table S2. Biomass, growth rate, and biomass productivity of wild-type and recombinant strains. Table S3. Fluorescence intensities of GFP in transformant mycelia. Figure S1. Plasmid map of pCambia0390-blpR-9PtrpC-gfp.

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

Lyu, M., Zeng, J., Zhou, Y. et al. Overlapping promoter library designed for rational heterogenous expression in Cordyceps militaris. Microb Cell Fact 21, 107 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: