Improvement of cordycepin production by an isolated Paecilomyces hepiali mutant from combinatorial mutation breeding and medium screening

ImageXue Cai1,2 · Jie‑Yi Jin1,2 · Bo Zhang1,2 · Zhi‑Qiang Liu1,2 · Yu‑Guo Zheng1,2

Received: 11 May 2021 / Accepted: 30 June 2021
© The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature 2021

Abstract
Cordycepin is a major bioactive compound found in Cordyceps sinensis that exhibits a broad spectrum of biological activi- ties. Here a Paecilomyces hepiali OR-1 strain was initially isolated from plateau soil for the bioproduction of cordycepin. Subsequently, strain modification including 60Co γ-ray and ultraviolet irradiation were employed to increase the cordycepin titer, resulted in a high-yield mutant strain P. hepiali ZJB18001 with the cordycepin content of 0.61 mg/gDCW, showing a 2.3-fold to that from the wild strain (0.26 mg/gDCW). Furthermore, medium screening based on Box-Behnken design and the response surface methodology facilitated the enhancement of cordycepin yield to the value of 0.96 mg/gDCW at 25 °C for 5 days in submerged cultivation with an optimized medium composition. The high cordycepin yield, rapid growth rate and stable genetic characteristics of P. hepiali ZJB18001 are beneficial in terms of costs and time for the industrialization of cordycepin production.
Keywords Cordycepin · 60Co-γ ray irradiation · Paecilomyces hepiali · Medium screening · Response surface methodology

Introduction
Cordycepin is a derivative of the nucleoside adenosine with the hydroxy group absent at the 3’ position of its ribose moiety (3’-deoxyadenosine, C10H10N5O3, Fig. 1) [1]. As the first antibiotic nucleoside isolated from Cordyceps genus, cordycepin is well-known for a variety of biological and pharmacological activities [2, 3]. Applications in cosme- ceutical as a bioactive ingredient are expanding as well as anti-photoaging and anti-pigmentation [4]. However, the market price of cordycepin keeps high over the years due to the depletion and slow growth of the wild resources of Cordyceps sinensis caused by predatory excavation and

* Zhi-Qiang Liu [email protected]
1 The National and Local Joint Engineering Research Center for Biomanufacturing of Chiral Chemicals, Zhejiang University of Technology, Hangzhou 310014, People’s Republic of China
2 Key Laboratory of Bioorganic Synthesis of Zhejiang Province, College of Biotechnology and Bioengineering, Zhejiang University of Technology, Hangzhou 310014, Zhejiang, People’s Republic of China

environmental destruction. How to realize the cordycepin production in an industrial scale has therefore attracted the interests of researchers especially in Asia.Biological production of cordycepin could be achieved by direct extraction from C. militaris fruiting body, or from its fermentation broth through solid-state or liquid fermentation [5, 6]. It takes 60 days for the C. militaris fruiting body to be ready for extraction and the culture conditions are com- plex and hard to control as well. In comparison, the produc- tion cycle of liquid fermentation is shorter (~ 15 days), and the fermentation process is easier to control [2, 7]. Liquid fermentation technology has therefore become the major means of cordycepin production with a bunch of work in cordycepin-producing fungi such as Cordyceps militaris, Cordyceps sinensis, Cordyceps kyushuensis, Aspergillus nidulan, and Paecilomyces hepiali, which exhibited high similarity in components and efficacy of active compounds [2, 4, 7, 8].The low cordycepin productivity has been improved by strategies including physical and chemical mutagenesis breeding, process optimization, genetic and metabolic engineering [9]. Among them, mutagenesis breeding technologies such as ultraviolet, γ-ray irradiation, plasma, and chemical mutagens have been considered as practical

 Chemical structure of adenosine (a) and cordycepin (b)approaches in the industry for decades [10]. γ-Rays pro- duced by the radioisotope 60Co or 137Cs are high-energy electromagnetic waves which are able to penetrate through the cell walls of mycelia, which would cause DNA sequence mutation. Liu et al. [11] reported a successful example of 60Co γ-irradiation on an edible mushroom Aga- ricus brasiliensis, resulting in a mutant with 50% lower Cd accumulation and 40.3% higher biological efficiency than those of the parent. As far as cordycepin production is concerned, proton beam irradiation has been employed on Cordyceps militaris with a cordycepin titer increased to 14.3 g/L [12–14]. Cordyceps militaris KYL05 was also mutated by UV irradiation with a 2.3-fold improvement in cordycepin production [15]. An extensive literature search showed that there has not yet any study reporting using 60Co γ-ray irradiation for cordycepin production improvement.

Medium optimization along with the response surface methodology (RSM) is another effective approach to enhance the yield and to determine key factors influenc- ing the process with limited experiments as it is a com- bination of mathematical and statistical techniques [16, 17]. A recent research has reported additional growth supplements such as nucleosides-hypoxanthine increased the cordycepin titer by 34-fold in mycelial biomass of Ophiocordyceps sinensis [18]. Tang et al. [19] evaluated different vegetable oils as the second carbon source of C. militaris resulting in 3.17 times higher cordycepin titer than control using 20 g/L peanut oil. Casein hydrolysate was also found effective as a nitrogen source with a 2.3- fold improvement to 445 mg/L [15].

The present work isolated and identified a P. hepiali mutant strain from the mutation library of the wild C. sinensis by 60Co γ-ray and ultraviolet irradiation. Further attempt was made to optimize the basic culture medium for achieving a high level of cordycepin by the mutant strain using Box-Behnken design and response surface method- ology in a four-factor three levels experiment.

Materials and methods
Chemical reagentsStandard cordycepin was purchased from Sigma-Aldrich (Shanghai, China). Methanol of HPLC grade was purchased from Honeywell Burdick & Jackson Chemical Ltd. (New Jersey, USA). All other chemical reagents were of analytical grade and commercially available.

Microorganism isolation and cultivation

Wild Cordyceps sinensis samples were collected from Qing- hai, China, washed with tap water (3 times), distilled water (3 times) and 0.1% mercuric chloride solution for surface disinfection, sucked dry with sterile paper on a clean bench. The body parts were dissected with a sterile scalpel and transplanted onto potato dextrose agar plate (PDA plate) containing 80 mg/L of gentamicin. The plates were incu- bated at 10 °C for everyday observation until white mycelia without contamination were selected for transplantation and cultivated on new PDA plates without gentamicin at 25 °C for screening. A single colony was inoculated into 50 mL of seed medium (glucose 20 g/L, peptone 5 g/L, yeast extract 2 g/L, KH2PO4 1 g/L and MgSO4 0.5 g/L) in a 250 mL flask, and incubated at 25 °C for 3 days. The fermentation medium contained glucose 20 g/L, peptone 5 g/L, KH2PO4 1 g/L and MgSO4 0.5 g/L. Unless otherwise specified, fermentations were carried out in a 250 mL flask at 25 °C and an initial pH of 6.0 for 5 days, with 2% (v/v) inoculation of the seed culture. The morphological observations were carried out by an optical microscope (DM3000, Leica, Germany) and an environmental scanning electron microscope (ESEM, SU-8010, Hitachi, Tokyo, Japan) with the samples treated as Abrusci et al. [20] reported for SEM observation.

The wild strain was cultivated in a 5 L bioreactor (Shang-hai Baoxing Biological Equipment Engineering Co. Ltd, China) with 150 mL inoculation and a working volume of 3 L using optimized fermentation media. Cells were incu- bated at 25 °C for 5 days with agitation at 150 rpm and 6 L/ min filter sterilized airflow with automatic pH adjustment by ammonia to keep the pH at 6.0 and foaming control by ster- ile polypropylene glycol addition. Samples were taken every 12 h to measure the dry cell weight and cordycepin titer.

Microorganism identification by 18S rDNA and Biolog system

The 18S rDNA of the newly isolated microorganism was extracted and enzymatically amplified according to a previ- ously described method [21], sequenced and compared in the

NCBI database by using the Basic Local Alignment Search Tool (BLAST, http://www.ncbi.nlm.nih.gov/BLAST). The neighbor-joining tree was constructed using the Clustal Omega program. The utilizations of carbon and nitrogen sources of the isolated microorganism were studied with Biolog system kit following the manual (Biolog, Hayward, CA, USA).

Combinatorial mutation breeding of 60Co γ‑ray and ultraviolet irradiation

The suspension of strain P. hepiali OR-1 with 1.0 × 106 cells/ mL was prepared with 0.1% of Tween80 for 60Co γ-ray and ultraviolet irradiation sequentially. The radiation process was conducted in the Irradiation Center of Zhejiang Academy of Agricultural Sciences (Hangzhou, China) using the HFY-YC type γ-ray radiation device with radiation doses of 100, 200, 400, 600, 800 and 1000 Gy. Eight samples were irradiated for each dose at the same time. The suspension was placed in ice bath until use after irradiation. The retrieved spore solution was diluted to reach an even spore concentration of1.0 × 103 spore/mL and 50 μL of the diluted spore solution was
added onto the PDA plate at 25 °C for 3 days. Colony counts were then performed to determine the death rates. Single colonies were picked for cultivation and cordycepin content determination.

The mutant strain with the highest cordycepin content obtained by 60Co γ-ray irradiation was stored and prepared in 106 spores/mL as a suspension solution, evenly placed onto a sterile plate to form a 2 mm-thick film for ultraviolet irradia- tion with a 15 W UV lamp (emission wavelength of 253.7 nm). The lengths of ultraviolet irradiation were of 5 min, 10 min, 15 min, 20 min, 25 min and 30 min, respectively, and immedi- ately placed on ice for 1 h in dark. A volume of 100 μL diluted samples (100 spore/mL) were coated on PDA plates in tripli- cate and incubated at 25 °C for 3 days in the dark.
The lethal rate, positive and negative mutation rates of 60Co γ-ray and ultraviolet irradiation were calculated accord- ing to Eqs. (1), (2) and (3), respectively. The positive/negative mutation rate was defined as the cordycepin yield increased/ decreased by over 10% comparing to the parent. Equations (1),
(2) and (3) were as follows:N0 indicates the number of viable bacteria before irradiation,
N1 indicates the number of viable bacteria after irradiation.
Study on the genetic stability of the mutant strain

The genetic stability of the obtained mutant strain P. hepiali ZJB18001 was tested by cultivation as the standard cultiva- tion procedure as described in “Microorganism identifica- tion by 18S rDNA and Biolog system” for 10 generations. Samples from each generation were taken for cordycepin quantification by HPLC.
Medium screening

Single factor optimization

Nine carbon sources of glucose, D-fructose, sucrose, malt- ose, mannitol, corn syrup, molasses, peanut oil, methyl oleate in a concentration of 20 g/L and eight nitrogen sources of peptone, yeast extract, beef extract, silkworm chrysalis, peanut powder, yeast extract, bran, soybean meal of 10 g/L were tested, respectively. Besides, metal ions including Ca2+, Mn2+, Zn2+, Fe2+, Cu2+, Al3+ and Fe3+ with sulfate ions of 4 mmol/L and precursors including adenosine, adenine, guanosine, guanine, uridine, uracil, thymidine, inosine and D-ribose of 2 g/L were added separately to further optimize the fermentation medium.

Response surface methodology

Four factors (peanut oil, yeast extract, metal ion Fe3+ and precursor adenosine) at three different levels (− 1, 0, 1) were following the Box-Behnken design (Table 1). A set of 29 experiments was generated with five replicates at the center point for enhanced cordycepin production. A second-order polynomial model was used to fit the responses to the values of the factors.

Validation of the model

Experiment was performed in triplicate, to verify the statisti-

N0 − N1
Lethal rate(%) = × 100%
N0

NΔCordycepin>+10%

(1)
cal model at the optimal values of four factors. The observed

Table 1 Concentration ranges of the four factors used in the central composite design

Positive mutation rate(%) =
× 100%
N1
(2)

Factor Actual levels of coded factors

– 1 0 + 1

NΔCordycepin<−10%
Negative mutation rate(%) = × 100%
N1

(3)
Peanut oil (g/L) 16 20 24
Yeast extract (g/L) 8 10 12
Fe2(SO4)3 (mmol/L) 1.2 1.6 2.0
Adenosine (g/L) 1.5 2 2.5 values of the cordycepin titer and the biomass concentra- tion after 5 days of fermentation were compared with the predicted values.

Determination of mycelia biomass and cordycepin concentration
the significance of the model and coefficients. The results were analyzed by Design Expert software (Version 8.0.5, Stat-Ease Inc., Minneapolis, MN, USA), and fitted with a quadratic regression model established for variable response, as described by Eq. (4):

i
Y = b0 + biXi + biiX2 + bijXiXj
(4)

The fermentation liquid was centrifuged at 8000 rpm for
where Y is the predicted results, X and X for the independ-

10 min at room temperature. The precipitated mycelia
i j
ent variables of the coded value;were washed twice with distilled water and freeze-dried on a Christ Alpha 2–4 LDD plus freeze dryer (Osterode,

linear coefficient;
b0 is intercepted; I is the
bii is square coefficient; bij (i ≠ j) is interac-

Germany) at − 55 °C and 1.2 × 104 mpa for 48 h to obtain the mycelia dry mass. The determination of cordycepin concentration was carried out by HPLC with an Agi- lent 1260 system equipped with a Unitary C18 column (4.6 mm × 250 mm) and a UV detector at 259 nm at 30 °C, gradiently eluted with ultrapure water (A) and methanol (B) at a flow rate of 0.8 mL/min in a process of 0–4 min, 15% B; 4–8 min, 20% B; 8–12 min, 25% B; 12–16 min, 20% B;
and 16–20 min, 15% B [22]. The retention time of the stand- ard cordycepin was 9.8 min. The fermentation broth was re-dissolved in 50% methanol for LC–MS analysis (Model 2020, Shimadzu, Kyoto, Japan).

Statistical analysis

All of the samples were analyzed in triplicate, and mean values were calculated with standard deviations of < 4%. The analysis of variance (ANOVA) was performed to evaluatetion coefficient [23].

Results and discussion
Isolation and identification of wild strain
Paecilomyces hepiali OR‑1

A single colony with a relatively high cordycepin titer was isolated and named as OR-1. The OR-1 colonies on PDA plate were white, flat, villous and dense with a neat edge in the early growth stage, and the color of the colony surface turned into orange, the hyphae were developed and raised in the middle with extended incubation time to the average size of 1.0 cm in diameter (Fig. S1). The morphology of OR-1 was also investigated by ESEM through which typical characteristics of aerial hyphae and spores were observed (Fig. 2). The hyphae of OR-1 were smooth and branched Mycelia and spores of Paecilomyces sp. OR-1 byESEM a mycelia observation chart magnified 2,000 times; b is enlarge by 8000 times; c, d are enlarged views of the myce- lium by 20,000 timeswith a diameter of ~ 2.0 μm (Fig. 2d). Biolog analysis of met- abolic fingerprints showed that OR-1 could make good use of 41 carbon sources, while 45 sources could not be used, and 9 sources were in the borderline (Table S1). According to the new model bacteria category proposed by the state key laboratory of microbiology of the Chinese Academy of Sciences in 2015, strain OR-1 was preliminarily identified as Paecilomyces genus.

To further identify these isolated fungi, 18S rDNA sequencing and BLAST were performed revealing that the isolated strain shares a 100% identity of 560 bp to Paecilo- myces hepiali (KX237743.1), and a 72.4% identity to Hirsu- tella sinensis (Fig. 3). Based on the above physiological and biochemical characteristics, the isolated strain was identified as Paecilomyces hepiali OR-1. In addition, the active com- ponent from the fermentation broth of P. hepiali OR-1 was confirmed by LC–MS referring to the standard cordycepin (MS(ESI−) m/z = 252.02) (Fig. 4).

Fermentation of wild strain P. hepiali OR‑1 for cordycepin production

The submerged fermentation of P. hepiali OR-1 was per- formed for 144 h at 25 °C with a final cordycepin titer of
0.26 mg/g, hard to meet the industrial production require- ment. In the first 84 h, the mycelium accumulation was in logarithmic growth from 3.3 g/L to 10.0 g/L (Fig. 5). During the subsequent stable growth period, dry mycelium weight maintained of ~ 11.4 g/L. Correspondingly, the cordycepin content started to increase significantly from 36 h until 84 h, accumulated from 0.05 mg/g to 0.25 mg/g. Though the dry mycelium weight increased by 1.8 g/L to the maxi- mum of 11.8 g/L during the extended fermentation period (84–144 h), the cordycepin content increased slightly by0.1 mg/g to the highest value of 0.26 mg/g. Further breed- ing mutagenesis would be performed on this wild strain P. hepiali OR-1 to obtain a high-cordycepin production strain as the current cordycepin titer of P. hepiali OR-1 is hard to satisfy the industrial production requirement.

The phylogenetic tree for Paecilomyces sp. OR-1 and related strains based on the 18S rDNA sequence. Numbers in parentheses are accession numbers of published sequences. Bootstrap values were based on 1000 replicates 60Co γ‑ray and ultraviolet irradiation on P. hepiali OR‑1 To improve the cordycepin production of the isolated strain P. hepiali OR-1, 60Co γ-ray and UV irradiation were employed to generate mutant libraries. The suitable irra- diation dose of 60Co γ-ray on mycelia varies substantially among different species, for example, 1.5 kGy was selected for the mycelia of Agaricus brasiliensis strain with 40.3% higher biological efficiency and 14.7% higher crude protein content [11]. The appropriate dose for P. hepiali OR-1 was examined of 0, 100, 200, 400, 600, 800 and 1000 Gy show-ing a death rate of 0, 18.4, 35.0, 51.7, 78.6, 95.2 and 100% respectively. The positive rates were calculated as well of 0, 5.0, 4.3, 5.0, 6.4 and 3.1% correspondingly. Therefore, 600 Gy was selected as the median lethal dose with a lethal- ity rate of 78.6% and a positive mutation rate of 6.4%. A mutant strain with a cordycepin yield of 0.60 mg/gDCW was obtained as M1 from 17 positive strains, which is 2.3-fold to that of the wild strain. UV irradiation was subsequently applied on M1 with a 20 min irradiation time, accompany- ing a lethality rate of 76.4% and a positive mutation rate of 5.98%. Finally, a mutant strain exhibited an average cordycepin yield of 0.61 mg/gDCW was named as ZJB18001, and stored at the China Center for Type Culture Collection (CCTCC M 2018091, Wuhan, China). The adopt of UV irra- diation was not as effective as 60Co γ-ray irradiation which increased by 0.34 mg/gDCW in this study, but still worth trial as it is a simple and fast approach with satisfied results in many other strains [24, 25].

The best mutant strain ZJB18001 was compared with the wild strain P. hepiali OR-1 in morphology, physiological and biochemical characteristics. The single colonies of both strains showed similar morphology of white and villous, and white globular cell clusters formed in liquid medium as well. ZJB18001 was able to use 47 kinds of carbon sources, 44 kinds of which were invalid, and 4 kinds of weak usability (Table S1), indicating that the mutant strain showed wider substrate spectrum than the wild strain. The HPLC result of the fermentation broth of mutant strain ZJB18001 was shown in Fig. S3 indicating a complicated composition and a cordycepin peak at 22 min.

The improvement of cordycepin titer in the fermentation media might be relevant to the chro- mosomal aberration or genetic mutation of key enzymes in the metabolic pathway of cordycepin synthesis in ZJB18001 after exposure to irradiation, however, genome sequencing of ZJB18001 should be further conducted to explore the underlying mechanism of cordycepin enhancement.

Single‑factor optimization of medium composition

As the limited knowledge of nutrient requirements of mutant strain P. hepiali ZJB18001, medium optimization Mass spectrometry results of cordycepin standards( a) and fermentation broth of Paecilomyces hepiali OR-1 (b)is a necessary step to improve the cordycepin content con- cerning the factors such as carbon source, nitrogen source, heavy metal content and precursor content. As a result, pea- nut oil was favorable for mycelia growth and cordycepin production among nine carbon sources by 7% enhancement in comparison with glucose (Fig. 6a). The yields of mycelia and cordycepin were most increased by yeast extract with 10.4% and 6% respectively, comparing to peptone (Fig. 6b). The optimal carbon and nitrogen sources improved the cordycepin titer from 0.61 mg/g to 0.74 mg/g indicating that P. hepiali ZJB18001 favors the acidic condition as pea- nut oil offers a relatively acidic environment. Metal ions have been reported to be participated in the regulation ofactivities of enzymes, and to change the permeability of cell membrane, thus influencing the metabolites out pump- ing [26]. The cordycepin titer was improved by 9.5% to the value of 0.81 mg/g with the addition of 4 mmol/L metal ion Fe3+ (Fig. 6c). Addition of 2 g/L adenosine enabled the cordycepin increasement to 0.93 mg/g by 13.5% and dry mycelia weight by 16.3% to the value of 18.2 g/L (Fig. 6d), confirming that adenosine is a key secondary metabolite in P. hepiali ZJB18001, the same as the synthetic pathway of cordycepin in C. sinensis (Fig. 7). Lin et al. [27] dem- onstrated that 5′-nucleotidase participated in cordycepin synthesis, hexokinase, and glucose phosphate isomerase in cordycepic acid synthesis in Hirsutella sinensis by real-time.

 Fermentation curve of wild strain Paecilomyces hepiali OR-1 in a 5-L fermenter

PCR, accompanying with significant 15.0, 5.3-, and 3.9- fold up-regulation, respectively. Pang et al. [28] speculated 3′,5′-cyclic-nucleotide phosphodiesterase might be related
with cordycepin biosynthesis in Paecilomyces hepialid by transciptome analysis. The metabolic pathway indicated that precursor substances of cordycepin might be as 3′-AMP synthesized from 2′,3′-cyclic AMP, which was a by-product from mRNA degradation of C. militaris (Fig. 7b) [9, 29]. In sum, peanut oil, yeast extract, Fe2(SO4)3 and adenosine were confirmed to be beneficial for cordycepin production and further optimization was carried out according to the response surface methodology.
Culture medium screening for cordycepin production enhancement

The Box-Behnken design-based experiments with four fac- tors in three levels were performed to select the most effec- tive culture medium specifically for P. hepiali ZJB18001 (Table 2). The steepest ascent experiment results showed the center points of four factors were of peanut oil 20 g/L, yeast extract 10 g/L, Fe2(SO4)3 1.6 g/L, and adenosine 2 g/L. The regression coefficients F value, p value and standard error were summarized and selected for ANOVAs of the Single-factor optimization of medium composition for the mutant strain Paecilomyces hepiali ZJB18001. a effects of car- bon sources on mycelia mass and cordycepin production; b effects of nitrogen sources on mycelia mass and cordycepin production; ceffects of metal ion on mycelia mass and cordycepin production; d effects of precursor substance on mycelia mass and cordycepin pro- duction

Imagemodel to analyze the relationship between the cordycepin yield and each variable (Table S2). According to the 3D response surface curve which is a symmetrical mound shape with a coordinate axis, changes on cordycepin yield could be observed by small variations on each factor concentration (Fig. 8), indicating the variables of square interaction (X12, X22, X32, X42) were all significant on the cordycepin yield.
The values of F, R2, R2 , R2 , C.V.% and Adeq Precision The possible synthetic routes of cordycepin. a, b are the syn- thetic pathway found in C. militaris and C. kyushensis; c is the syn- thetic pathway found in C. cicadaewere 23.78, 0.96, 0.92, 0.78, 2.35% and 17.42, respectively, for cordycepin yield (Y). The p value is less than 0.0001 implying the interaction between each factor was significant, and this model can be used to guide the optimization of the mutant strain ZJB18001.By linear regression quadratic polynomial fitting on the experimental data, the following second-order polyno- mial Eq. (5) was established to explain the production of cordycepin:
Experimental design and results of the Box-Behnken design Runs Coded values and real values (g/L) Cordycepin titer (mg/g)

X1
(peanut oil)
X2
(yeast extract)
X3 (Fe2(SO4)3)
X4
(adenosine)
Experimental Predicted

1 16 (-1) 8 (-1) 1.6 (0) 2 (0) 0.71 ± 0.21 0.77
2 24 (1) 8 (-1) 1.6 (0) 2 (0) 0.85 ± 0.32 0.81
3 16 (-1) 12 (1) 1.6 (0) 2 (0) 0.82 ± 0.20 0.74
4 24 (1) 12 (1) 1.6 (0) 2 (0) 0.82 ± 0.14 0.88
5 20 (0) 10 (0) 1.2 (-1) 2.5 (1) 0.84 ± 0.17 0.89
6 20 (0) 10 (0) 2 (1) 2.5 (1) 0.74 ± 0.23 0.71
7 20 (0) 10 (0) 1.2 (-1) 1.5 (-1) 0.60 ± 0.14 0.66
8 20 (0) 10 (0) 2 (1) 1.5 (-1) 0.86 ± 0.18 0.81
9 16 (-1) 10 (0) 1.6 (0) 2.5 (1) 0.84 ± 0.16 0.79
10 24 (1) 10 (0) 1.6 (0) 2.5 (1) 0.72 ± 0.24 0.67
11 16 (-1) 10 (0) 1.6 (0) 1.5 (-1) 0.63 ± 0.19 0.68
12 24 (1) 10 (0) 1.6 (0) 1.5 (-1) 0.83 ± 0.20 0.87
13 20 (0) 8 (-1) 1.2 (-1) 2 (0) 0.63 ± 0.16 0.67
14 20 (0) 12 (1) 1.2 (-1) 2 (0) 0.82 ± 0.31 0.86
15 20 (0) 8 (-1) 2 (1) 2 (0) 0.85 ± 0.25 0.79
16 20 (0) 12 (1) 2 (1) 2 (0) 0.78 ± 0.27 0.71
17 16 (-1) 10 (0) 1.2 (-1) 2 (0) 0.66 ± 0.16 0.62
18 24 (1) 10 (0) 1.2 (-1) 2 (0) 0.70 ± 0.13 0.77
19 16 (-1) 10 (0) 2 (1) 2 (0) 0.83 ± 0.23 0.87
20 24 (1) 10 (0) 2 (1) 2 (0) 0.76 ± 0.14 0.71
21 20 (0) 8 (-1) 1.6 (0) 2.5 (1) 0.90 ± 0.34 0.84
22 20 (0) 12 (1) 1.6 (0) 2.5 (1) 0.72 ± 0.15 0.75
23 20 (0) 8 (-1) 1.6 (0) 1.5 (-1) 0.62 ± 0.17 0.67
24 20 (0) 12 (1) 1.6 (0) 1.5 (-1) 0.77 ± 0.24 0.72
25 20 (0) 10 (0) 1.6 (0) 2 (0) 0.99 ± 0.15 0.91
26 20 (0) 10 (0) 1.6 (0) 2 (0) 0.98 ± 0.16 0.91
27 20 (0) 10 (0) 1.6 (0) 2 (0) 0.99 ± 0.17 0.91
28 20 (0) 10 (0) 1.6 (0) 2 (0) 0.97 ± 0.20 0.91
29 20 (0) 10 (0) 1.6 (0) 2 (0) 0.97 ± 0.11 0.91

 3-D response surface curves of the cordycepin production of Paecilomyces hepiali ZJB18001 as a function of yeast extract powder and peanut oil (a), yeast extract powder and Fe3+ (b), yeast extractpowder and adenosine (c), peanut oil and Fe3+ (d), peanut oil and adenosine (e), Fe3+ and adenosine (f)

Y = 0.98 + 0.022 ∗ X1 + 0.017 ∗ X2 + 0.044 ∗ X3

+ 0.034 ∗ X4 − 0.036 ∗ X1 ∗ X2 − 0.037 ∗ X1
∗ X3 − 0.079 ∗ X1 ∗ X4 − 0.067 ∗ X2 ∗ X3
1
– 0.094 ∗ X2 ∗ X4 − 0.090 ∗ X3 ∗ X4 − 0.099 ∗ X2
2.14 g/L. A mean cordycepin content of 0.96 mg/gDCW was obtained, representing a good agreement with the predicted value of 0.98 mg/gDCW. These results confirmed the validity and precision of the constructed model, and the accumula- tion amount of cordycepin was greatly increased by 1.57-

– 0.096 ∗ X2 − 0.11 ∗ X2 − 0.11 ∗ X2
(5)
fold comparing to that with the original medium.

2 3 4
Besides the high cordycepin production of P. hepiali

where Y represents the predicted response of cordycepin content (mg/g) of mutant ZJB18001, X1, X2, X3, and X4 are the coded values for peanut oil, yeast extract, Fe2(SO4)3 and adenosine, respectively.
The model was validated by conducting experiments in triplicates under the predicted medium: peanut oil 21 g/L, yeast extract 9.8 g/L, Fe2(SO4)3 4.68 mmol/L and adenosine vZJB18001 in comparison with reported strains (Table 3), such as C. militaris [14] and Hirsutella sinensis [27, 30], the cordycepin production of 0.96 mg/g by ZJB18001 in the opti- mized medium was threefold to that from Hirsutella sinensis [27]. P. hepiali is advantageous in a 5-day cultivation time, whereas other strains need 8–20 days, which is favored by industry in the aspect of cost reduction and time-saving.

Summary of cordycepin production, culture medium and cultivation conditions reported in recent yearsStrain Method Medium composition (g/L) Cordycepin Culture condition Ref

Carbon source
Nitrogen source
Metal ion Precursor
production in the broth (g/L)
Temp(°C) Time
(Day)

Cordyceps militaris
Submerged culture
Glucose 86.2
Yeast
extract 93.8
1/10 Vogel’s medium
– 14.3 25 20 [14]

Hirsutella sinensis
Submerged culture
Glucose 10
Molas- ses 10,
silkworm chrysalis pow-
der 5, soybean cake powder 10, yeast
extract 5
MgSO4 0.1,
KH2PO4 0.2
– 3.85
(0.35 mg/g)
16 12 [27]

Hirsutella sinensis
Submerged culture
Glucose 10
Molas- ses 10,
silkworm chrysalis pow-
der 5, soybean cake powder 10, yeast
extract 5
MgSO4 0.1,
KH2PO4 0.2
Adenosine-5′- monophosphate 4
16.35
(1.09 mg/g)
16 12 [27]

Cordyceps militaris CICC 14,014

Cordyceps militaris KYL05
Static
culture

Submerged culture
Glucose 40
Peanut oil 20

Glucose 20
Tryptone 10, yeast
extract 6

Casein hydro- lysate 20
MgSO4·7H2O 0.5, K2HPO4·3H2O 0.5, KH2PO4
0.5
K2HPO4·3H2O 0.2, KH2PO4
0.1, MgSO4·7H2O 0.2
Adenine 1.6 5.29 25 20 [19]

– 0.45 25 6 [15]

Paecilomy- ces hepiali
Submerged culture
Glucose 20
Peptone 5 KH2PO4 1,
MgSO4 0.5
– 3.11
(0.26 mg/g)
25 5 This study

Paecilomy- ces hepiali ZJB18001

Paecilomy- ces hepiali ZJB18001
Submerged culture

Submerged culture
Peanut oil 20

Peanut oil 21
Yeast
extract 10

Yeast
extract 9.8
KH2PO4 1,
MgSO4 0.5, Fe2(SO4)3
4 mmol/L
KH2PO4 1,
MgSO4 0.5, Fe2(SO4)3
4.68 mmol/L
Adenosine 2 8.10
(0.61 mg/g)
Adenosine 2.14 17.86
(0.96 mg/g)
25 5 This study
25 5 This study

Conclusions
The aim of this study was to screen a strain that are capa- ble of producing high content of cordycepin as an alter- native to C. sinensis in an industrial-scale production. Irradiation mutagenesis implemented on the wild strains was a common but effective approach for the enhance- ment of bioactive metabolites that have been used for
decades in industry when the genetic background of a certain strain is not clear yet. The results provide us a combined approach of 60Co γ-ray and ultraviolet irradia- tion and medium screening to improve the productivity of bioactive compounds by strains. Further efforts such as genome and transcriptome analysis should be made to discover the mechanism behind the cordycepin production enhancement that could be possibly transferred to other strains with high cordycepin yield but a slow growth rate for metabolic engineering.

Supplementary Information The online version contains supplemen- tary material available at https://doi.org/10.1007/s00449-021-02611-w.

Acknowledgements None.

Author contribution Xue Cai: Data curation, Writing-original draft and Writing-Review and Editing; Jie-Yi Jin: Investigation and Validation; Bo Zhang: Methodology and Writing-Review and Editing; Zhi-Qiang Liu: Conceptualization, Writing-Review and Editing; Yu-Guo Zheng: Supervision.

Funding This work was financially supported by the National Natural Science Foundation of China (32070099 and 32000898), the Natural Science Foundation of Zhejiang Province (LQ21B060006), the Fun- damental Research Funds for the Provincial Universities of Zhejiang (RF-C2019005).

Declarations

Conflict of interest Authors have no declarations of conflict of interest.

References
1. Verma AK (2020) Cordycepin: a bioactive metabolite of  Cordyceps militaris and polyadenylation inhibitor with therapeu- tic potential against COVID-19. J Biomol Struct Dyn. 1. https://doi. 1. org/10.1080/07391102.2020.1850352
2. Lin LT, Lai YJ, Wu SC, Hsu WH, Tai CJ (2018) Optimal con- ditions for cordycepin production in surface liquid-cultured Cordyceps militaris treated with porcine liver extracts for sup- pression of oral cancer. J Food Drug Anal 26(1):135–144. 2. https:// 2. doi.org/10.1016/j.jfda.2016.11.021
3. Lee CT, Huang KS, Shaw JF, Chen JR, Kuo WS, Shen G, Grumezescu AM, Holban AM, Wang YT, Wang JS, Hsiang YP, Lin YM, Hsu HH, Yang CH (2020) Trends in the immunomodula- tory effects of Cordyceps militaris: total extracts, polysaccharides and cordycepin. Front Pharmacol 11:575704. 3. https://doi.org/10. 3. 3389/fphar.2020.575704
4. Kunhorm P, Chaicharoenaudomrung N, Noisa P (2019) Enrich- ment of cordycepin for cosmeceutical applications: culture sys- tems and strategies. Appl Microbiol Biotechnol 103:1681–1691. https://doi.org/10.1007/s00253-019-09623-3
5. Cui JD (2015) Biotechnological production and applications of Cordyceps militaris, a valued traditional Chinese medicine. Crit Rev Biotechnol 35(4):475–484. 5. https://doi.org/10.3109/07388551. 5. 2014.900604
6. Sari N, Suparmin A, Kato T, Park EY (2016) Improved cordycepin production in a liquid surface culture of Cordyceps militaris iso- lated from wild strain. Biotechnol Bioprocess Eng 21:595–600. https://doi.org/10.1007/s12257-016-0405-0
7. Tang J, Liu Y, Zhu L (2014) Optimization of fermentation condi- tions and purification of cordycepin from Cordyceps militaris. Prep Biochem Biotechnol 44(1):90–106. 7. https://doi.org/10.1080/ 7. 10826068.2013.833111
8. Chang CY, Lue MY, Pan TM (2005) Determination of adeno- sine, cordycepin and ergosterol contents in cultivated Antrodia camphorata by HPLC method. J Food Drug Anal 13(4):338–342. https://doi.org/10.38212/2224-6614.2569
9. Yang L, Li G, Chai Z, Gong Q, Guo J (2020) Synthesis of cordycepin: current scenario and future perspectives. Fungal Genet Biol 143:103431. 9. https://doi.org/10.1016/j.fgb.2020. 9. 103431
10. Zou SP, Zhong W, Xia CJ, Gu YN, Niu K, Zheng YG, Shen YC (2015) Mutagenesis breeding of high echinocandin B producing strain and further titer improvement with culture medium optimi- zation. Bioprocess Biosyst Eng 38(10):1845–1854. 10. https://doi.org/ 10. 10.1007/s00449-015-1425-4
11. Liu P, Yuan J, Jiang Z, Wang Y, Weng B, Li G (2019) A lower cadmium accumulating strain of Agaricus brasiliensis produced by 60Co-γ-irradiation. Lwt-Food Sci Technol 114:108370. 11. https:// 11. doi.org/10.1016/j.lwt.2019.108370
12. Das SK, Masuda M, Hatashita M, Sakurai A, Sakakibara M (2008) A new approach for improving cordycepin productivity in surface liquid culture of Cordyceps militaris using high-energy ion beam irradiation. Lett Appl Microbiol 47(6):534–538. 12. https:// 12. doi.org/10.1111/j.1472-765X.2008.02456.x
13. Masuda M, Das SK, Fujihara S, Hatashita M, Sakurai A (2011) Production of cordycepin by a repeated batch culture of a Cordyceps militaris mutant obtained by proton beam irradiation. J Biosci Bioeng 111(1):55–60. 13. https://doi.org/10.1016/j.jbiosc. 13. 2010.08.018
14. Masuda M, Das SK, Hatashita M, Fujihara S, Sakurai A (2014) Efficient production of cordycepin by the Cordyceps militaris mutant G81–3 for practical use. Process Biochem 49(2):181–187. https://doi.org/10.1016/j.procbio.2013.10.017
15. Lee SK, Lee JH, Kim HR, Chun Y, Lee JH, Yoo HY, Park C, Kim SW (2019) Improved cordycepin production by Cordyceps mili- taris KYL05 using casein hydrolysate in submerged conditions. Biomolecules 9(9):461. https://doi.org/10.3390/biom9090461
16. Kalil SJ, Maugeri F, Rodrigues MI (2000) Response surface analy- sis and simulation as a tool for bioprocess design and optimiza- tion. Process Biochem 35(6):539–550. 16. https://doi.org/10.1016/ 16. S0032-9592(99)00101-6
17. Chaudhary P, Chhokar V, Choudhary P, Kumar A, Beniwal V (2017) Optimization of chromium and tannic acid bioremediation by Aspergillus niveus using Plackett-Burman design and response surface methodology. AMB Express 7:201. 17. https://doi.org/10. 17. 1186/s13568-017-0504-0
18. Kaushik V, Singh A, Arya A, Sindhu SC, Sindhu A, Singh A (2020) Enhanced production of cordycepin in Ophiocordyceps sinensis using growth supplements under submerged conditions. Biotechnol Reports 28:e00557. 18. https://doi.org/10.1016/j.btre. 18. 2020.e00557
19. Tang J, Qian Z, Wu H (2018) Enhancing cordycepin production in liquid static cultivation of Cordyceps militaris by adding veg- etable oils as the secondary carbon source. Bioresour Technol 268:60–67. https://doi.org/10.1016/j.biortech.2018.07.128
20. Abrusci C, Martín-González A, Del Amo A, Catalina F, Col- lado J, Platas G (2005) Isolation and identification of bacteria and fungi from cinematographic films. Int Biodeterior Biodegrad 56(1):58–68. https://doi.org/10.1016/j.ibiod.2005.05.004
21. Liu ZQ, Li Y, Ping LF, Xu YY, Cui FJ, Xue YP, Zheng YG (2007) Isolation and identification of a novel Rhodococcus sp. ML-0004 producing epoxide hydrolase and optimization of enzyme produc- tion. Process Biochem 42(5):889–894. 21. https://doi.org/10.1016/j. 21. procbio.2007.01.009
22. Lee HH, Kang N, Park I, Park J, Kim I, Kim J, Kim N, Lee JY, Seo YS (2017) Characterization of newly bred Cordyceps mili- taris strains for higher production of cordycepin through HPLC and URP-PCR analysis. J Microbiol Biotechnol 27(7):1223–1232. https://doi.org/10.4014/jmb.1701.01043
23. Ye Z, Wang W, Yuan Q, Ye H, Sun Y, Zhang H, Zeng X (2016) Box-Behnken design for extraction optimization, characteriza- tion and in vitro antioxidant activity of Cicer arietinum L. hull

polysaccharides. Carbohydr Polym 147:354–364. https://doi.org/ 10.1016/j.carbpol.2016.03.092
24. Takaku H, Ebina S, Kasuga K, Sato R, Ara S, Kazama H, Mat- suzawa T, Yaoi K, Araki H, Shida Y, Ogasawara W, Ishiya K, Aburatani S, Yamazaki H (2021) Isolation and characterization of Lipomyces starkeyi mutants with greatly increased lipid produc- tivity following UV irradiation. J Biosci Bioeng 131(6):613–621. https://doi.org/10.1016/j.jbiosc.2021.01.006
25. Kang CK, Yang JE, Park HW, Choi YJ (2021) Enhanced lycopene production by UV-C irradiation in radiation-resistant Deinococcus radiodurans R1. J Microbiol Biotechn 30(12):1937–1943. https:// doi.org/10.4014/jmb.2009.09013
26. Cui JD, Zhang YN (2012) Evaluation of metal ions and sur- factants effect on cell growth and exopolysaccharide production in two-stage submerged culture of Cordyceps militaris. Appl Biochem Biotechnol 168(6):1394–1404. https://doi.org/10.1007/ s12010-012-9865-7
27. Lin S, Liu ZQ, Xue YP, Baker PJ, Wu H, Xu F, Teng Y, Brathwaite ME, Zheng YG (2016) Biosynthetic pathway analysis for improv- ing the cordycepin and cordycepic acid production in Hirsutella
sinensis. Appl Biochem Biotechnol 179:633–649. https://doi.org/ 10.1007/s12010-016-2020-0
28. Pang F, Wang L, Jin Y, Guo L, Song L, Liu G, Feng C (2018) Transcriptome analysis of Paecilomyces hepiali at different growth stages and culture additives to reveal putative genes in cordycepin biosynthesis. Genomics 110(3):162–170. https://doi. org/10.1016/j.ygeno.2017.09.008
29. Xia Y, Luo F, Shang Y, Chen P, Lu Y, Wang C (2017) Fungal cordycepin biosynthesis is coupled with the production of the safeguard molecule pentostatin. Cell Chem Biol 24(12):1479– 1489. https://doi.org/10.1016/j.chembiol.2017.09.001
30. Jin LQ, Xu ZW, Men XH, Zhang B, Liu ZQ, Zheng YG (2020) Enhancement of protoplast preparation and regeneration of Hir- sutella sinensis based on process optimization. Biotechnol Lett 42(11):2357–2366. https://doi.org/10.1007/s10529-020-02958-2
Publisher’s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.