Skip to main content

Inducing secondary metabolite production of Aspergillus sydowii through microbial co-culture with Bacillus subtilis



The co-culture strategy which mimics natural ecology by constructing an artificial microbial community is a useful tool to activate the biosynthetic gene clusters to generate new metabolites. However, the conventional method to study the co-culture is to isolate and purify compounds separated by HPLC, which is inefficient and time-consuming. Furthermore, the overall changes in the metabolite profile cannot be well characterized.


A new approach which integrates computational programs, MS-DIAL, MS-FINDER and web-based tools including GNPS and MetaboAnalyst, was developed to analyze and identify the metabolites of the co-culture of Aspergillus sydowii and Bacillus subtilis. A total of 25 newly biosynthesized metabolites were detected only in co-culture. The structures of the newly synthesized metabolites were elucidated, four of which were identified as novel compounds by the new approach. The accuracy of the new approach was confirmed by purification and NMR data analysis of 7 newly biosynthesized metabolites. The bioassay of newly synthesized metabolites showed that four of the compounds exhibited different degrees of PTP1b inhibitory activity, and compound N2 had the strongest inhibition activity with an IC50 value of 7.967 μM.


Co-culture led to global changes of the metabolite profile and is an effective way to induce the biosynthesis of novel natural products. The new approach in this study is one of the effective and relatively accurate methods to characterize the changes of metabolite profiles and to identify novel compounds in co-culture systems.


Natural products (NPs) are an important historical source of many useful drugs and other chemical agents, of which microbial secondary metabolites represent a significant part [1]. Numerous novel secondary metabolites have been isolated from marine fungi and 70–80% of them have good biological activities such as anti-cancer, anti-bacteria, anti-parasite and free-radical scavenging. Moreover, some compounds have been marketed as commercial drugs through clinical research [2]. However, with the progress of scientific research, researchers have found that repeated discoveries of known metabolites are increasing. The fact that the biosynthetic potential has eluded is mostly explained by the observation that many genes are transcriptionally silent under standard culture conditions, causing their products inaccessible [3]. Moreover, analyses of microbial whole genome sequences indicate that microbes contain many thousands of biosynthetic gene clusters, which encode a plethora of compounds that are not identified when cultured under standard laboratory conditions [4]. To overcome these impasses, several approaches have been developed, such as nontargeted metabolic engineering, epigenetic modification, and chemical synthesis. Among such approaches, the co-culture method draws increasing attention to stimulate the production of novel natural products.

Co-culture of different microorganisms can imitate the natural microbial environment, and the silent biosynthetic gene clusters are transcriptionally activated by environmental stimuli [5]. The chemical cues released by other microbes can cause various defense responses, including the changes of mycelial morphology, synthesis of diverse secondary metabolites, and production of extracellular enzymes. Indeed, these activated defensive metabolites can act as chemical cues that can trigger a series of transcriptional activation [6]. Recently, there have been many successful studies on the use of microbial co-culture to induce new secondary metabolites. For instance, Zuck et al. demonstrated that co-culture of Aspergillus fumigatus and Streptomyces peucetius induced the production of four formyl xanthine analogs that were not generated in pure culture, of which two were compounds with novel structures, and compound 2 showed significant inhibitory activity against several cell lines [7]. Moreover, Wu et al. co-cultured Bacillus amyloliquefaciens and Trichoderma asperellum and found that the production of antibacterial substances was significantly higher than that in pure culture. When the inoculation ratio was 1:1, the production of specific amino acids was improved [8]. Therefore, co-culture is regarded as a useful research method for effectively inducing the production of metabolites.

Aspergillus sydowii can produce various secondary metabolites which are increasingly utilized in pharmaceuticals, food and chemicals, such as endoglucanases with industrial application value, enzymes with inhibitory activity against protein tyrosine phosphatase A of M. tuberculosis, sesquiterpenoids with antimicrobial and antiviral activities, and alkaloids with activity against S. aureus and S. epidermidis [9]. However, the analysis of the whole genome of A. sydowii revealed a number of genes for the biosynthesis of compounds that have not been observed when cultured under standard conditions [10]. The previous study has showed that the addition of 5-azacytidine, an epigenetic modifier, to the broth of A. sydowii induced the production of the (S)-(+)-sydonol which potentiated the insulin-stimulated glucose consumption, suggesting that the metabolites of A. sydowii obtained through nontargeted metabolic engineering might be developed into antidiabetic agents [11]. Recently, the enzymes from the protein tyrosine phosphatases (PTPs) superfamily are emerging as potential new drug targets for type 2 diabetes. For example, protein tyrosine phosphatase 1b (PTP1b) [12] is a negative regulator of insulin action in the insulin receptor signaling pathway, SH2-containing protein tyrosine phos-phatase-1 (SHP1) is a negative regulator in signaling pathways, which regulates glucose homeostasis through the modulation of insulin signaling in liver and muscle [13], and Leukocyte common antigen (CD45) is the receptor for some ligands, which can regulate the recruitment of SHP-1 [14]. Some microbial metabolites, such as varic acid analogues from fungi, showed selective inhibitory activities against PTPs [15]. Thus, the activities of metabolites of co-culture on the PTPs and their potential to be used in diabetes treatment is worth expecting.

The conventional method for the study of metabolites in co-culture systems is to separate and purify the compounds corresponding to each peak newly detected in HPLC, and to analyze their structures by means of MS, UV, IR, and NMR [16]. However, the conventional method is inefficient and time-consuming, and only products with high contents can be identified. The elucidation of trace newly biosynthesized metabolites in co-culture systems is still challenging. Furthermore, the overall changes of metabolite profiles during the co-culture cannot be displayed. In these years, metabolomics, which is mainly aided by the advances in analytical technologies, such as high-resolution mass spectrometry (MS), is primarily associated with comprehensive analysis of small-molecule compounds that can be found in biological samples. Some useful tools, for instance, computation-based MS-DIAL [17] and MS-FINDER [18] programs, and the web-based global natural product social molecular network (GNPS) [19] have also been developed to predict the structure of the metabolites. However, in most cases, only one single tool was used in structure predictions, and the accuracy of the predictions is still a concern. In addition, a web-based tool, MetaboAnalyst, which combined multivariate statistics to identify spectral features that are statistically different between two (or more) different sample populations, is useful in the statistical and functional analysis of metabolomic data [20]. It has been reported that these tools have been used in the analysis of the metabolite profile of microorganisms regulated by epigenetics [21]. However, to our best knowledge, there has been no report on the application of MetaboAnalyst in the co-culture of microbes.

In the present study, the fungus A. sydowii was co-cultured with the bacterium B. subtilis, and an integrated metabolomics approach, composed of MetaboAnalyst, MS-DIAL, MS-FINDER, and GNPS was developed to analyze the MS/MS data of the co-culture. The changes in the metabolite profile were characterized, and the newly biosynthesized compounds were identified. The purification and NMR analysis of part of the newly biosynthesized compounds were performed to verify the accuracy of the new approach. The activities of newly biosynthesized compounds against protein tyrosine phosphatases (PTPs) were also evaluated.


Microbial interaction induced changes of the metabolite profile

The co-culture of twenty microorganisms with A. sydowii on bran medium showed different degrees of induction between the cultures, among which B. subtilis could significantly induce A. sydowii to produce metabolites (Additional file 1: Fig. S1). After 12 days, the color of the hyphae of A. sydowii turned from dark green to light green in co-culture, and red-brown exudate was generated at the junction between B. subtilis and A. sydowii (Fig. 1). Moreover, a significant deadlock model was observed. This phenomenon indicated that during co-culture, A. sydowii and B. subtilis generated compounds due to the stress response at the confrontation zone, inhibiting the growth of the other. In order to further explore this phenomenon, we collected the confrontation zone and analyzed the metabolites.

Fig. 1
figure 1

Microorganisms on bran agar plate. a Pure culture of A. sydowii on Day 12. b Pure culture of B. subtilis on Day 12. c Co-culture of A. sydowii and B. subtilis on Day 12. A. sydowii was in the center and B. subtilis was at both sides. A reddish-brown exudate was observed between the two microorganisms

The extracts from the bran medium of co-culture and pure culture were compared by LC–MS/MS, and 206 strong signal features whose intensity was over 10% of the highest intensity peaks were detected. The partial least squares discriminant analysis (PLS-DA) of these peaks revealed the intrinsic variation in the data set. In the score plot, the samples from the co-culture were clearly separated from the two pure cultures, indicating the changes of metabolite profile (Fig. 2a). The heatmap generated by hierarchical clustering analysis (HCA) of these 206 features based on the MS data showed that co-culture caused global changes in the metabolomes (Fig. 2d). The heatmap also revealed that 25 features were identified only in co-culture, indicating that about 12.1% of the candidate features were newly biosynthesized during co-culture. In addition, 156 features were recorded in both pure culture and co-culture. Among them, 70 features in the co-culture system were significantly decreased, while 4 features in the co-culture system were up-regulated when compared with the pure culture of B. subtilis. On the contrary, there were only 8 features in the co-culture that were significantly decreased, and 37 features in the co-culture that were up-regulated when compared with the pure culture of A. sydowii (Fig. 3). In the loading plot of PLS-DA, the 25 newly biosynthesized features were mainly deviated from the center and clustered into the lower right zone of the plot (Fig. 2b). These features showed good linear correlation. Only the features that had a large contribution to the classification generated by co-culture were distributed on this line. The features that contributed more to the classification were closer to the lower right, while the features that contributed less to the classification were clustered on the upper left of the line and were closer to the origin (Additional file 1: Fig. S2). In the meantime, the variable importance in projection (VIP) score data indicated newly biosynthesized features (N1N4, N7, N13, and N20) with monoisotopic mass of m/z 168.4234, 266.1459, 282.1436, 282.4537, 353.1765, 402.1640, and 480.3248, respectively, were ranked in top features detected by VIP score (Fig. 2c). These data indicated that the newly biosynthesized features in co-culture made important contribution to group classification.

Fig. 2
figure 2

The PLS-DA, VIP and Heatmap of metabolomics data of co-culture and their pure-cultures on Day 12 is shown. a The score plot of the data analyzed by LC-HRMS. b The loading plot of the data analyzed by LC-HRMS. c The top compounds ranked based on the VIP score. The colored boxes on the right indicate the relative concentrations of the corresponding metabolite in each group. d Hierarchical clustering analysis (HCA) of the most significantly variable 206 features among the samples corresponding to the three different groups and represented on a heatmap (ranging from red color for high abundance to blue for low abundance). Data was acquired from three independent biological replicates

Fig. 3
figure 3

Overview of the number of the induced features in the co-culture of A. sydowii and B. subtilis

Metabolomics study of newly biosynthesized metabolites in the co-culture

To understand the structure of the newly biosynthesized metabolites, the 25 features were identified with the integrated approach. Here, we demonstrate our results using four annotation levels (Level 1–4). Level 1: the structures were annotated on MS-DIAL linked MS/MS databases by the characteristic product ions and neutral losses; Level 2: the metabolite ions were converted into structural information and the structures were annotated by the structure elucidation tool (MS-FINDER); Level 3: the structures were annotated putatively by the correlation with the known structures with the assistance of the network analysis tool (GNPS); Level 4: the structures were identified by separation, purification and NMR spectrum analysis. Next, the identification of feature N2 was discussed in detail as an example. The adduct ions of feature N2 were detected as m/z 265.1459 [M−H] and m/z 325.2437 [M+CH3COOH−H], suggesting that the monoisotopic mass was 266.1518. The PLS-DA analysis indicated that this feature was detected only in the co-culture and contributed greatly to the cluster. Firstly, the structure of N2 was identified by MS-DIAL (Level 1), which used deconvolution algorithm to obtain the retention time (RT) and m/z data sets of MS/MS data. The compounds were then annotated by comparing the characteristic products and neutral losses of features with the public MS/MS databases [22]. As no structural candidate was obtained from Level 1, the structure of N2 was annotated through MS-FINDER (Level 2), which was embedded in MS-DIAL software version 3.90. MS-FINDER was a strategy for computational MS/MS fragmentations. In this software, all isomer structures of the predicted formula were retrieved from metabolome databases, and the structure was predicated based on a combined weighting score considering bond dissociation energies, mass accuracies, fragment linkages, and nine rules of hydrogen rearrangement proposed during bond cleavages in low-energy collision-induced-based fragmentation [18]. After comparing in silico spectra and the structures provided by MS-FINDER, N2 was identified as sydonic acid based on the fact that the mass peaks (m/z 265.1459, 253.4632, 249.1126, 180.0446, 137.0338, and 93.0326) matched well with the MS/MS database (Additional file 1: Fig. S3) with the lowest mass error of 1.0942 ppm and the highest structure score of 7.55. Similarly, among the other 20 features, 5 features were identified through Level 1 process and 15 features were identified through Level 2 process. The detailed information of the metabolites was summarized in the Table 2. There were 5 major classes of metabolites induced by co-culture, including sesquiterpenes, macrolides, esters, polyketides, and flavonoids. These metabolites of microorganism, are usually not generated in the normal condition of microorganisms, and can only be synthesized under certain stress. These compounds were reported to participate in the defense and communication between microbial cells, promote metabolism, and have a certain bacteriostatic effect [23].

In order to verify the validity of the identification approach, five compounds (N1, N2, N3, N4, and N13) which had higher VIP scores in PLS-DA analysis indicating larger contributions to the cluster, were isolated and purified through silica gel column chromatography, ODS column chromatography and preparative HPLC from the confrontation zone of the co-culture. According to the results, N1 (20 mg), N2 (60 mg), N3 (19 mg), N4 (13 mg), and N13 (21 mg) with purity over 95% were obtained. After analyzing of the NMR data, the compounds were identified as Orsellinic acid (N1) [24], Sydonic acid (N2) [25], (7S)-(−)-10-Hydroxysydonic acid (N3) [23], (R)-(−)-Hydroxysydonic acid (N4) [26], and Macrolactin A (N13) [27], which were consistent with those identified by the approach above, suggesting the credibility of the approach. The structures of the 7 compounds were shown in Fig. 4 and the detailed NMR data were provided in Additional file 1: compounds information.

Fig. 4
figure 4

Structures of the compounds

Identification of the novel metabolites in the co-culture

There were still 4 features (N6, N7, N9, and N20) that did not match any features in the public MS/MS spectrum library. To elucidate the structures of these potential novel metabolites generated through co-culture, MS/MS data were analyzed with Level 3 process which was assisted by GNPS platform and manual dereplication. The GNPS approach can capture similar structures and analog features into the same cluster regardless of retention time in the LC–MS.

The GNPS data indicated that three induced features, including N6 (m/z 350.1610), N7 (m/z 352.1765), and N9 (m/z 368.1713) were clustered, suggesting that these features have very close structural relationships (Fig. 5). As none of the three features was identified within the LC–MS/MS database, compound N7, with the highest content in LC–MS data, was separated and purified. Compound N7 was obtained as white powder. The UV absorption was at 213 nm, 254 nm, and 298 nm. The molecular formula C18H27O6N was indicated by the ESI-HRMS at m/z 354.1908 [M+H]+ (calculated as 354.1872), indicating 6 degrees of unsaturation. The 1H NMR and 13C NMR of N7, N2 and the reference data of N2 [25] were shown in Table 1. The 13C NMR and HMQC spectra indicated the presence of a total of 18 carbon signals attributable to one carboxyl carbon at δc 172.05 (C-10); one ketone carbon δc 166.14 (C-7); three methyl groups at δc 28.48 (C-8′), δc 22.57 (C-6′) and δc 22.37 (C-7′); four methylene groups at δc 61.35 (C-9), 41.60 (C-2′), 21.33(C-3′) and 38.78(C-4′); five methines including δc 22.37 (C-7′), three aromatic carbons at δc 126.64 (C-6), δc 117.47 (C-5) and δc 115.10 (C-3), one oxygenated carbon at δc 55.53 (C-8); four quaternary carbons including three aromatics at δc 154.74 (C-2), δc 135.86 (C-1) and δc 133.50 (C-4), one oxygenated atom at δc 75.04 (C-1′).

Fig. 5
figure 5

Molecular network analysis of newly discovered features in the co-culture. a A constructed subnetwork to dereplicate co-culture and pure culture; b two newly identified features (m/z 369.171, 351.161) and the confirmed feature (m/z 353.177) were clustered together

Table 1. 1H and 13C NMR spectral data for compound N7 in DMSO-d6, N2 in CD3OH and the reference data of sydonic acid in CD3OH (500 MHz for 1H NMR and 125 MHz for 13C NMR)

The inspection of the 1H NMR spectrum indicated the presence of δH 7.35 (d, J = 8.05 Hz, H-6), 7.28 (dd, J = 8.05, 1.35 Hz, H-5) and 7.24 (d, J = 1.35 Hz, H-3), which was a typical spectrum of 1,3,4-trisubstituted benzene ring. This sub-structure was confirmed by the relationship between H-5 (δH 7.28) to H-6 (δH 7.35) in 1H–1H COSY, and H-3 (δH 7.24) to C-1 (δc 135.86), H-5 (δH 7.28) to C-1, and H-6 (δH 7.35) to C-4 (δc 133.5) in HMBC. The correlations from H-7′ (δH 0.77) to H-2′ (δH 1.94, 1.66) in 1H–1H COSY indicated the existence of a hexane substructure. The HMBC relationship between H-5 (δH 7.28) with C-7 (δc 166.14) and H-6 (δH 7.35) with C-1′ (δc 75.04) implied the linkage of C-4 to the C-1 position. The position of the amine bond was determined by the relationship of H-8 to C-7 in HMBC (Fig. 6; Additional file 1: Figs. S4–S10).

Fig. 6
figure 6

The key COSY and HMBC of compound N7

The absolute configuration of N7 was deduced based on the comparison of experimental data and calculated ECD curves by Gaussian 09. The conformers were optimized using DFT at the B3LYP/6-31G (d) level in methanol. The energies were calculated through the TDDFT methodology at the B3LYP/6-31G (d, p) level in MeOH with PCM model (Additional file 1: Figs. S11–S14; Tables S1–S8). The calculated CD spectrum of N7 (1′R, 8S) agreed well with the experimental CD curve (Fig. 7), indicating that absolute configuration of N7 was 1′R, 8S, and was named as Serine sydonate.

Fig. 7
figure 7

Calculated and experimental ECD spectrums of compound N7 (black, experimental in MeOH; red, blue, pink, green, calculated at the B3LYP/6-31G (d) //B3LYP/6-31G (d, p) level in MeOH)

The structures of compounds N6 and N9 were mainly determined by the LC–MS/MS data from the negative-ion mode and compared with N7. Comparison of the fragment ions of the compounds showed some common fragments of m/z 224.0562, m/z 194.0458 and m/z 150.0563, indicating that these features had the same backbone structure and belonged to a series of structural derivatives (Additional file 1: Fig. S15). The structures of compounds N6 and N9 were determined by comparing the negative-ion mode LC–MS/MS data with N7 because these three compounds shared similar MS/MS pattern and the content of N6 and N9 was low in broth. For the LC–MS/MS data of N7, the abundant fragments were m/z 334.1656 and 322.1660, which arose from molecular ion m/z 352.1765 by the loss of H2O (18 Da) and two methyl groups (30 Da), respectively. The fragments of m/z 304.1557 and 290.1758 were generated by further facile loss of H2O (18 Da) and CH3OH (32 Da) from the fragment of m/z 322.1654, respectively. The fragments of m/z 224.0562, m/z 194.0458 and m/z 150.0563 indicated the substructure of serine substituted hydroxy-benzoic acid. The monoisotopic mass of N6 was 350.1610 [M−H], which was common neutral loss of 2 Da of N7, indicating that the compound N6 was most likely the dihydrogen product of N7. The fragments of m/z 224.0562, m/z 194.0458 and m/z 150.0563 suggested the existence of the substructure of serine substituted hydroxy-benzoic acid in N6. The stable m/z 302.1394 fragment, which represented the conjugated olefin structure formed by the dehydroxylation of hydroxyl methylheptane, suggested that the double bond was located at the 4′ position. Thus, compound N6 was determined as 4′-alkenyl serine sydonate. For compound N9, the fragments of m/z 224.0562, m/z 194.0458 and m/z 150.0563 were also the characteristics of the substructure of serine substituted hydroxy-benzoic acid. The residue mass was 18 Da higher than that of N6, suggesting that N6 was the dehydration product of N9. The stable m/z of 302.1394 in N6 and m/z of 320.1497 in N9 also indicated that dehydration occurred in the substructure of hydroxyl methylheptan. Thus, the structure of N9 was determined as hydroxyl serine sydonate. However, whether this hydroxyl group was located at 4′ or 5′ position could not be determined by LC–MS/MS data alone (Fig. 8). N6 and N9 (4′-hydroxyl or 5′-hydroxyl serine sydonate) were found to be novel compounds after database searching.

Fig. 8
figure 8

Annotated MS/MS spectrum of m/z 369.1713, 351.1610 and the confirmed feature (m/z 353.1765) acquired by LTQ-Orbitrap-XL in negative mode

As compound N20 cannot be connected with the other metabolites in Level 3, it was forwarded to Level 4 for direct isolation and purification. This compound was isolated as pale-yellow creamy solid. The molecular formula C31H44O4 was indicated by the ESI-HRMS at m/z 503.3134 [M+Na]+ (calculated for 480.3240), indicating 10 degrees of unsaturation. The UV absorption was at 212 nm and 273 nm. The 1H NMR (DMSO-d6, 500 MHz) and 13C NMR (DMSO-d6, 125 MHz) data were provided in the Additional file 1: compounds information. The 1D NMR data were consistent with the data of the known compound Macrolactin U identified by Xue et al. [28], and its relative configuration was 4S, 5S. However, the methyl group H3-29 (δ 1.10, d) and H2-6 (δ 2.56, m) were defined as trans due to the NOESY relationship between H3-29 and H-5 (δ 4.25, m) in compound N20 (Additional file 1: Fig. S16-S23). Thus, the relative configuration of C-4 and C-5 were S and R, respectively, (Fig. 9), which was different from S and S in Macrolactin U. Therefore, compound N20 was identified as the isomer of Macrolactin U and named as Macrolactin U′, which is also a novel compound. Unfortunately, the absolute configuration of N20 was not determined yet as no obvious difference was identified between ECD and crystal of N20, which could not be obtained due to the limited amount of the compound.

Fig. 9
figure 9

The key COSY and NOESY of compound N20

Thus, a total of 25 features induced only in the co-culture were identified by the combination of the computational approach (MS-DIAL), the web-based tools (GNPS and MetaboAnalyst) with chemical isolation and purification (Table 2, Additional file 1: Fig. S24). Four of the features were novel metabolites, including two compounds confirmed by NMR. Five known compounds were also purified to verify the validity of the approach.

Table 2 List of induced features only in the co-culture of A. sydowii with B. subtilis analyzed by LC-HRMS

Biological activity assay

The isolated compounds N1N4, N7, N13 and N20 were evaluated for their anti-nematode activity and antidiabetic activity. Among the compounds, compound N3 showed a certain degree of anti-nematode activity with an IC50 of 50 μM. Furthermore, compounds N2N4, N7 and N13 exhibited potent activity against SHP1 and PTP1b, both of which are targets for the development of diabetes (Table 3). In addition, compounds N7 and N13 displayed inhibition activities against CD45 with IC50 values of 16.0 μM and 17.9 μM, respectively.

Table 3 Activities of compounds


Microbial metabolites have always been considered as a very important source of new drugs, due to their various biological activities such as anti-bacteria, anti-oxidation, and anti-tumor. When two microorganisms are co-cultured, new metabolites can be biosynthesized by one or both microorganisms as a result of interspecific crosstalk or induction by biochemical signaling molecules [29]. For example, Akone et al. [30] co-cultured Chaetomium sp. with B. subtilis, obtaining 5 new compounds and 7 known compounds. However, how to characterize the overall changes in the metabolite profile induced by co-culture and to identify the newly biosynthesized metabolites is still a complicated and challenging task.

The fragmentation pattern in the MS/MS spectrum represents a specific feature of a certain compound. The structures and chemical properties of the molecules determine the observable fragmentation patterns in MS/MS data. Therefore, similar fragmentation patterns of related compounds are used as indications of chemical relatedness [31, 32]. In the field of structural prediction of natural products, some useful tools draw great attention. For instance, computational MS-DIAL program can be used to obtain deconvoluted spectra from high-resolution LC–MS data; another computation based MS-FINDER program can be used for structure elucidation of unknown HR-MS spectra through fragment comparison and MS database searching [18, 33]. GNPS provides a visualization approach to detect sets of spectra from related molecules (molecular networks), even when the spectra themselves do not match any known compounds. Using these tools, some metabolites of co-culture were predicted. For example, Ernest et al. [34] analyzed 9 co-cultures of marine-adapted fungi and phytopathogens by GNPS and annotated 18 molecular clusters, 9 of which were exclusively produced in co-cultures. Several clusters contained compounds that could not be annotated to any known compounds, suggesting that they are putatively newly metabolites. However, as only GNPS was mainly involved in these studies, and due to the limited volume of MS/MS library in GNPS, few structures can be predicted. Most structures of newly biosynthesized compounds have not been elucidated yet. Recently, some researchers have also tried to integrate multiple tools to assist structure elucidation. Lai et al. [35] showcased a combined workflow, including GC–MS metabolome database, MS-DIAL and MS-FINDER, to analyze the volatile organic compounds, and three biomarkers and two propofol derivatives were annotated successfully in over 110, 000 biological samples. To our best knowledge, there is still a lack of integrated, effective and accurate strategy to reveal the changes of metabolite profile and characteristics in microorganism following treatment to activate silenced genes, including co-culture.

In this study, MetaboAnalyst, MS-DIAL, MS-FINDER, and GNPS were integrated with the publicly available spectral library to compare the MS/MS data, including common losses of MS and fragmentation similarity, while obtaining the same molecules, analogs, or metabolism families, thereby facilitating structural analysis. Analysis of the co-culture of A. sydowii and B. subtilis by this new approach revealed 206 features induced in the confrontation zone, and 25 features which occupied 12.1% of the detected features were newly induced by co-culture. Especially, 4 features (N7, N20, N9, and N6) were identified as novel compounds. All the 25 newly biosynthesized metabolites were identified by the integrated approach, and the accuracy of the integrated approach was also partially verified by the isolation, purification and spectrum analysis of five newly biosynthesized metabolites with high content. These results suggested that this new approach provided an effective and time-efficient manner to characterize the overall changes of the metabolite profile and to elucidate the structures of metabolites simultaneously. In the integrated approach, N6 and N9, which were derivatives of N7 with a low content, were detected by GNPS molecular network based on the similar fragment pattern. Their structures were further elucidated with the assistance of MS-DIAL and MS-FINDER programs. These results suggested the new approach was effective to discover trace derivatives of metabolites, and would help to understand the global metabolite profile changes in the co-culture system. Interestingly, the newly biosynthesized features were linearly correlated in the loading plot of the PLS-DA analysis in MetaboAnalyst (Additional file 1: Fig. S2). In the previous co-culture study, a similar phenomenon was presented in the PCA analysis using SIMCA-P software when two fungi, Trametes versicolor and Ganoderma applanatum were co-cultured (Fig. 1 in reference Xu et al. [36]), although the authors did not describe this linear correlation. These data suggested that using the linear correlation rule of newly biosynthesized metabolites in PLS-DA or PCA analysis, the new biosynthesized metabolites, especially the newly biosynthesized metabolites with low content, might be easily figured out, although the mechanism of this rule still needs to be clarified.

Analysis of the structural features of the part of newly biosynthesized metabolites also revealed their producing microorganism. Sydonic acid (N2) has been reported as a typical metabolite of Aspergillus sp. [25], suggesting that this compound and its structurally similar compounds N3, N4, N6, N7, and N9 were produced by A. sydowii under the inducing stress of B. subtilis. Similarly, N21 has been reported to be produced by the species of Bacillus [28], suggesting that this compound and its analogues N13 and N20 were supposed to be produced by B. subtilis.

The structures of 25 newly biosynthesized metabolites in the co-culture can be categorized into five classes, including macrolides, sesquiterpenes, esters, polyketides, and flavonoids. N13 and N20 belonged to macrolides. Macrolide antibiotics have multiple conjugated double bonds, hydroxyl side chain groups, and macrolide skeleton structures. This class of antibiotics has no effect on bacteria, but has an inhibitory effect against fungi. Macrolides can interact with sterols on the membrane of fungal cells, causing the leakage of small molecules and ions in the cell content from the transmembrane pores, eventually leading to the death of fungal cells [37]. For instance, Macrolactin A (N13), was reported to display meaningful antifungal activity with MIC values of 0.04–0.3 mM [38]. Compounds N2, N3 and N4 are classified into sesquiterpenoids, which are widely distributed in nature with anti-bacteria, anti-inflammation and immunoregulatory activities [39]. For example, (7S)-(−)-10-Hydroxysydonic acid (N3), was reported to display inhibitory activities against S. aureus with IC50 values ranging from 31.5 to 41.9 μM [40]. (R)-(−)-Hydroxy Sydonic acid (N4), showed broad spectrum activities against S. aureus and B. cereus, with MIC less than 25 μM [26]. These results, together with the analysis of the producing microorganisms of these compounds, indicated that in order to exert the antagonistic effect, A. sydowii and B. subtilis induced the biosynthesis of macrolides and sesquiterpenoids, respectively, to inhibit the growth of the opponent. In addition, inducing the production of sesquiterpenoids by the fungus might play important roles in symbiosis, such as enhancing the immune regulation of bacteria, removing free radicals, enhancing the vitality of bacteria and gaining more nutrients effectively. As the response of antagonistic and symbiosis effects, expression of silent genes was activated and the new metabolites were biosynthesized effectively during the co-culture of A. sydowii and B. subtilis.

The biological assay indicated that purified newly biosynthesized metabolites showed specific inhibitory activities against PTP1b, SHP1 and CD45. The PTP1b assay data of compounds N2N7 in this study indicated that the side chains of these compounds influenced the activities. The addition of hydroxyl groups to the side chain lowered the activity significantly. However, the addition of hydroxyl groups to the side chain had no obvious effect on the inhibitory activity against CD45. Further research is still needed to reveal the structure–activity relationship among these compounds, which will help to design new agents for the treatment of diabetes or immunomodulation.


Co-culture of A. sydowii and B. subtilis increased the diversity of metabolites. The new integrated approach in this study, which includes MetaboAnalyst, MS-DIAL, MS-FINDER, and GNPS, explored the overall changes of microbial metabolites profile of co-culture, and elucidated the structural information of 25 compounds. Four of the compounds identified are novel. The structures of the other 5 compounds were purified and their NMR data were analyzed to verify the accuracy of the new approach. These data suggest that the new approach is effective and reliable for the rapid identification of metabolites. The biological activities of 7 compounds isolated showed relatively strong inhibitory activity of N2 to PTP1b, indicating that the co-culture strategy could induce the production of bioactive secondary metabolites, and provide a valuable platform for the discovery of more novel secondary metabolites. The co-culture strategy will also contribute to the revelation of the metabolic mechanisms that can activate silent genes.


General experimental procedures

HPLC analysis was performed with a Waters HPLC system equipped with a 2998 detector and a 1525 pump. Routine detection wavelengths were at 235, 254, 280, and 340 nm. Twenty (20) μL of the samples was injected to a Shimadzu TC-C18 column (10 × 250 mm, 5 μm), and the following gradient was used (mobile phase A: 0.2% CH3COOH in H2O, mobile phase B: acetonitrile): 0–30 min (20–80% B), 30–35 min (80–100% B), 35–40 min (100% B) at 37℃ with a flow rate of 1 mL/min. Compounds were prepared by silica gel column chromatography and an LC3000 semi-preparative HPLC system (Beijing Chuang Xin Tong Heng Science and Technology Co., Ltd). The analytical, semi-preparative and preparative HPLC was performed using an ODS column from Shimadzu Co. (TC-C18, 10 × 250 mm, 5 μm), an YMC semi-preparative column (YMC-Pack Pro C18 RS, 10 × 250 mm, 5 μm), and an YMC preparative column (YMC-Pack ODS-A, 20 × 250 mm, 10 μm). The NMR data were recorded on a Bruker 500 MHz spectrometer from Bruker Co.

Microorganisms and co-cultivation experiment

Aspergillus sydowii was isolated from a piece of deep-sea mud from Dalian, China. In order to further excavate the secondary metabolites of the strain, A. sydowii was co-cultured with B. subtilis. The fungal and bacterial strains were activated in potato dextrose agar (PDA) medium (200 g Potato/L, 20 g dextrose/L and 15 g agar/L) for 3 days. Then, each strain was suspended in 2 mL of sterile water. To establish individual pure cultures, 80 μL of bacterial suspension was inoculated into a 90 mm petri dish containing 20 mL of bran medium (100 g bran/L, 20 g dextrose/L, 15 g agar/L). For the co-culture, 80 μL of each bacterial suspension of A. sydowii and B. subtilis was inoculated approximately 10 mm apart on the PDA medium. The plates were incubated at 28 °C for 12 days.

Measurement of the metabolome

The extracts were dissolved in 150 μL methanol and centrifuged at 16000×g for 10 min. The supernatants were transferred into HPLC autosampler vials and analyzed on a LTQ Orbitrap XL mass spectrometer (Thermo Fisher Scientific, Hemel Hempstead, UK) at a flow rate of 0.6 mL/min. The ESI conditions were as follows: the spray voltage was fixed at 4200 V; the sheath gas pressure was 35 arb; the auxiliary gas pressure was 10 arb and the heater temperature were 320 °C, and the capillary temperature was 300 °C. The mass scanning range was m/z 50–1200 Da in centroid mode with a scan rate of 1.5 spectra/s. The mass detection was performed by an electrospray source functioning in positive and negative ion mode at 15,000 resolving power. The mass measurement was externally calibrated before the experiment. Each full MS scan was followed by data-dependent MS/MS on the three most intense peaks using stepped collision-induced dissociation (35% normalized collision energy, isolation width 2 Da, activation Q 0.250). Twenty (20) μL of the samples was separated by a Shimadzu TC-C18 column (10 × 250 mm, 5 μm). The mobile phase A was water with 0.1% acetic acid and the mobile phase B was acetonitrile with 0.1% acetic acid. The elution gradient of reversed-phase liquid chromatography was as follows: 0–10 min, 20% B; 10–30 min, 20–80% B; 30–35 min, 80–85% B; 35–40 min, 100% B; 40–45 min, 25% B; All the samples had three independent biological replicates. The solvent (MeOH) and pure culture were injected under the same conditions as controls.

Metabolites profile and structure analysis

In order to fully exploit the differences of the metabolite profile in co-culture and pure cultures, MS-DIAL (Version 3.90), the computational approach which helps to rapidly characterize the structure of the metabolites [22], and MetaboAnalyst [41], the web-based tools for comprehensive metabolomic data analysis and interpretation, were integrated. This approach mainly includes: (1) determination of monoisotopic mass of peaks with MS-DIAL. In MS-DIAL, the adduct ion dictionary were defined as: [M+H]+, [M+Na]+, [M+K]+, [M−H2O+H]+, and [2M+H]+ for data from positive ion mode, and [M−H], [M−H2O−H], [M+HCOOH−H], [M+CH3COOH−H]and [2M−H] for data from negative ion mode. The monoisotopic mass of each peak was determined when at least two adduct ions matched the adduct ion dictionary. (2) Peak list alignment with MS-DIAL. The MS/MS data were converted to abf format by Analysis Base File Converter, and then subjected to MS-DIAL program to find the peak list alignment. The MS tolerance was set as 0.01 Da, the minimum peak height was set as 1 × 107, and the maximum charge was set to 2. (3) Multivariate analysis of the global metabolites profile with MetaboAnalyst. The aligned data were uploaded to MetaboAnalyst, and the data were first normalized by the sum and auto scaled. The data then were analyzed with PLS-DA to reveal the global profile changes, and the heatmap that could show clustering of the features and visualize the differences between groups was also obtained. (4) Structural identification of the metabolites assisted with MS-DIAL and GNPS. This step mainly included four levels. Level 1: structure annotated on MS-DIAL linked MS/MS databases by the characteristic product ions and neutral losses. The MS/MS public databases mainly include ReSpect, BMDMS-NP, MetaboBASE, Fiehn/Vaniya and natural product library in positive and/or negative manner. Level 2: structure annotated on MS-DIAL linked MS-FINDER program. The metabolite ions were converted into structural information with MS-FINDER. The number of carbon atoms and formula can be determined and the structural formula of all substructures were defined. Compared with the public spectral databases including NIST 14, MassBank, Metlin, ReSpect, and MetaboBase, the compounds with monoisotopic mass error within ± 5 ppm and a structure score higher than 5 were screened for mass spectral peak matching. Then the structures were searched on Reaxys and SciFinder database to confirm whether they were derived from natural products. Finally, the ontology for all substructure forms was defined. Level 3: structure annotation assisted by GNPS. LC–MS/MS data was uploaded to GNPS to create the network between the metabolites. Thus, the features whose structure scores was less than 5 or its fingerprints could not match any compounds in MS/MS database, which might be the novel compounds, might be correlated to the other structures. If any structures in the molecular network can be identified in LC–MS/MS database, the structure of other features in the network can be deduced by comparing the difference between the MS/MS spectrum of unknown and available structures. Otherwise, at least one of the compounds in the network would be separated, purified and its NMR spectrum was analyzed to elucidate the structure of the feature, and then the other structure of the network was deduced accordingly. In GNPS analysis, the LC–MS/MS data were first converted to mzXML format by MS Convert and processed by MZmine 2 [42] and then uploaded to GNPS. The parent mass tolerance was set as 2.0 Da. The ion tolerance was set as 0.5 Da. The maximum connected components value was set as 19 and the minimum cluster size was set as 2. All matches between the network spectra and the library spectra were required to have a score above 0.7 and at least six matched peaks. The output of the molecular networks was visualized using Cytoscape (Version 3.6.1) [43]. Level 4: structure identification by separation, purification and NMR spectrum analysis. The features whose structures could not be determined by Level 1–3 were separated and purified by column chromatography, and analyzed with 1D and 2D NMR spectrum. To verify the veracity of the identification approach, some structures with higher VIP scores in PLS-DA analysis were also separated, purified and analyzed with NMR data.

Extraction and isolation of the metabolites

After cultured for 12 days, the confrontation zones of co-culture were collected and soaked in ethyl acetate to extract the compounds of interest. The extract was evaporated under reduced pressure, and 30 g of the crude extract was obtained. Then, the crude extract was subjected to a silica gel column and a gradient elution using N-hexane/dichloromethane (90:10 → 0:100 over 30 min, 0:100 hold for 10 min) at a flow rate of 12 mL/min. Three fractions, F1–F3, were obtained from the separation. The fraction F1 was further purified using DAISO ODS (20% acetonitrile to 100% acetonitrile over 35 min) at a flow rate of 20 mL/min, followed by preparative HPLC with acetonitrile-H2O (30% isocratic) to yield N1 (20 mg), N3 (19 mg), N4 (13 mg) and N7 (50 mg). The fraction F2 was processed in the same manner as F1 with acetonitrile-H2O (60% isocratic) to yield N2 (60 mg) and N13 (21 mg). The fraction F3 from the separation was further purified by a YMC preparative column and a YMC semi-preparative column at 3 mL/min with acetonitrile-H2O (75% isocratic) to yield N20 (8 mg).

Computational details

The theoretical calculations of compound N7 were performed using Gaussian 09. Firstly, the conformations at B3LYP/6-31G (d) level were optimized in MeOH and the theoretical of ECD was determined using Time Dependent Density Functional Theory (TDDFT) at B3LYP/6-31G (d, p) level in MeOH. Secondly, the ECD spectra was simulated using Gaussian function with band width σ = 0.30 eV. Finally, the ECD spectra of compound N7 was obtained by weighing the Boltzmann distribution rate of each geometric conformation.

Protein tyrosine phosphatase 1b inhibitory assay

The PTP1b inhibitory activity of the tested compounds was measured at 37 °C using p-nitrophenyl phosphate (pNPP) as the substrate. The reaction was performed in a 96-well plate (final volume of 150 μL) and incubated for 30 min in the assay buffer (50 mM citrate (pH 6.0), 0.1 M NaCl, 1 mM EDTA, and 1 mM dithiothreitol) at 37 °C. Subsequently, the reaction was terminated by the addition of 10 M NaOH and the amount of p-nitrophenyl was determined by measuring the absorbance at 405 nm.

Availability of data and materials

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



Natural products

A. sydowii :

Aspergillus sydowii

B. subtilis :

Bacillus subtilis


Mass spectrometry


Global Natural Product Social Molecular Network


Protein tyrosine phosphatase


Hierarchical clustering analysis


Variable importance in projection


Leukocyte common antigen


Potato dextrose agar


Partial least squares discriminant analysis


Time Dependent Density Functional Theory


  1. Bertrand S, Schumpp O, Bohni N, Monod M, Wolfender JL. De novo production of metabolites by fungal co-culture of Trichophyton rubrum and Bionectria ochroleuca. J Nat Prod. 2013;76:1157–65.

    Article  CAS  PubMed  Google Scholar 

  2. Raad II, Zakhem AE, Helou GE, Jiang Y, Kontoyiannis DP, Hachem R. Clinical experience of the use of voriconazole, caspofungin or the combination in primary and salvage therapy of invasive aspergillosis in haematological malignancies. Int J Antimicrob Agents. 2015;45:283–8.

    Article  CAS  PubMed  Google Scholar 

  3. Ma YJ, Zheng LP, Wang JW. Inducing perylenequinone production from a bambusicolous fungus Shiraia sp. S9 through co-culture with a fruiting body-associated bacterium Pseudomonas fulva SB1. Microbial Cell Fact. 2019;18:121.

    Article  CAS  Google Scholar 

  4. Chiang YM, Chang SL, Oakley BR, Wang CC. Recent advances in awakening silent biosynthetic gene clusters and linking orphan clusters to natural products in microorganisms. Curr Opin Chem Biol. 2011;15:137–43.

    Article  CAS  PubMed  Google Scholar 

  5. Karuppiah V, Vallikkannu M, Jie T, et al. Simultaneous andsequential based co-fermentations of Trichoderma asperellum GDFS1009 and Bacillus amyloliquefaciens 1841: a strategy to enhance the gene expression and metabolites to improve the bio-control and plant growth promoting activity. Microbial Cell Fact. 2019;18:185–185.

    Article  CAS  Google Scholar 

  6. Shang Z, Salim AA, Capo RJ. Chaunopyran A: co-cultivation of marine mollusk-derived fungi activates a rare class of 2-alkenyl-tetrahydropyran. J Nat Prod. 2017;80:1167–72.

    Article  CAS  PubMed  Google Scholar 

  7. Zuck KM, Shipley S, Newman DJ. Induced production of N-formyl alkaloids from Aspergillus fumigatus by co-culture with Streptomyces peucetius. J Nat Prod. 2011;74:1653–7.

    Article  CAS  PubMed  Google Scholar 

  8. Wu Q, Ni M, Dou K, Tang J, Ren J, Yu C, et al. Co-culture of Bacillus amyloliquefaciens ACCC11060 and Trichoderma asperellum GDFS1009 enhanced pathogen-inhibition and amino acid yield. Microbial Cell Fact. 2018;17:155.

    Article  CAS  Google Scholar 

  9. Liu N, Peng S, Yang J, Cong Z, Lin X, Liao S, et al. Structurally diverse sesquiterpenoids and polyketides from a sponge-associated fungus Aspergillus sydowii SCSIO41301. Fitoterapia. 2019;135:27–32.

    Article  CAS  PubMed  Google Scholar 

  10. Chiang YM, Szewczyk E, Nayak T, Davidson AD, Sanchez JF, Lo HC, et al. Molecular genetic mining of the Aspergillus secondary metabolome: discovery of the emericellamide biosynthetic pathway. Cell Chem Biol. 2008;15:527–32.

    CAS  Google Scholar 

  11. Chung YM, Wei CK, Chuang DW, El-Shazly M, Hsieh CT, Asai T, et al. An epigenetic modifier enhances the production of anti-diabetic and anti-inflammatory sesquiterpenoids from Aspergillus sydowii. Bioorgan Med Chem. 2013;21:3866–72.

    Article  CAS  Google Scholar 

  12. Dodd GT, Xirouchaki CE, Eramo M, Mitchell CA, Andrews ZB, Henry BA, et al. Intranasal targeting of hypothalamic PTP1B and TCPTP reinstates leptin and insulin sensitivity and promotes weight loss in obesity. Cell Rep. 2019;28:2905–22.

    Article  CAS  PubMed  Google Scholar 

  13. Dubois MJ, Bergeron S, Kim HJ, Dombrowski L, Perreault M, Fournès B, et al. The SHP-1 protein tyrosine phosphatase negatively modulates glucose homeostasis. Nat Med. 2006;12:549–56.

    Article  CAS  PubMed  Google Scholar 

  14. Greer SF, Justement LB. CD45 regulates tyrosine phosphorylation of CD22 and its association with the protein tyrosine phosphatase SHP-1. J Immunol. 1999;162:5278–86.

    Article  CAS  PubMed  Google Scholar 

  15. Sun W, Zhang B, Zheng H, Zhuang C, Xiu Z. Trivaric acid, a new inhibitor of PTP1b with potent beneficial effect on diabetes. Life Sci. 2016;169:52–64.

    Article  PubMed  CAS  Google Scholar 

  16. Kamdem RST, Wang H, Wafo P, Ebrahim W, Zkaya FC, Makhloufi G, et al. Induction of new metabolites from the endophytic fungus Bionectria sp. through bacterial co-culture. Fitoterapia. 2017;124:132–6.

    Article  PubMed  CAS  Google Scholar 

  17. Tsugawa H, Ikeda K, Takahashi M, Satoh A, Arita M. A lipidome atlas in MS-DIAL 4. Nat Biotechnol. 2020;38:1159–63.

    Article  CAS  PubMed  Google Scholar 

  18. Tsugawa H, Kind T, Nakabayashi R, Yukihira D, Tanaka W, Cajka T, et al. Hydrogen rearrangement rules: computational MS/MS fragmentation and structure elucidation using MS-FINDER software. Anal Chem. 2016;88:7946–58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Wang M, Carver JJ, Phelan VV, Sanchez LM, Garg N, Peng Y, et al. Sharing and community curation of mass spectrometry data with global natural products social molecular networking. Nat Biotechnol. 2016;34:828–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Xia J, Wishart DS. Using MetaboAnalyst 3.0 for comprehensive metabolomics data analysis. Curr Protoc Bioinform. 2016;55:14.10.11-14.10.91.

    Article  Google Scholar 

  21. Triastuti A, Vansteelandt M, Barakat F, Trinel M, Haddad M. How histone deacetylase inhibitors alter the secondary metabolites of Botryosphaeria mamane, an endophytic fungus isolated from Bixa orellana L. Chem Biodivers. 2019;16:e1800485.

    Article  PubMed  CAS  Google Scholar 

  22. Tsugawa H, Nakabayashi R, Mori T, Yamada Y, Takahashi M, Rai A, et al. A cheminformatics approach to characterize metabolomes in stable-isotope-labeled organisms. Nat Methods. 2019;16:295–8.

    Article  CAS  PubMed  Google Scholar 

  23. Li XD, Li XM, Xu GM, Zhang P, Wang BG. Antimicrobial phenolic bisabolanes and related derivatives from Penicillium aculeatum SD-321, a deep sea sediment-derived fungus. J Nat Prod. 2015;78:844–9.

    Article  CAS  PubMed  Google Scholar 

  24. Fang MJ, Fang H, Li WJ, Huang DM, Wu Z, Zhao YF. A new diphenyl ether from Phoma sp. strain, SHZK-2. Nat Prod Res. 2012;26:1224–8.

    Article  CAS  PubMed  Google Scholar 

  25. Shinji K, Takanori M, Junsuke M, Kazuaki T, Noboru T. Isolation and absolute stereochemistry of optically active sydonic acid from Glonium sp. (Hysteriales, Ascomycota). Biosci Biotechnol Biochem. 2009;73:203–4.

    Article  CAS  Google Scholar 

  26. Wang CY, Liu YF, Cao F, Wang CY. Bisabolane-type sesquiterpenoids from a gorgonian-derived Aspergillus sp. fungus induced by DNA methyltransferase inhibitor. Chem Nat Compd. 2016;52:1129–32.

    Article  CAS  Google Scholar 

  27. He S, Wang H, Yan X, Zhu P, Chen J, Yang R. Preparative isolation and purification of macrolactin antibiotics from marine bacterium Bacillus amyloliquefaciens using high-speed counter-current chromatography in stepwise elution mode. J Chromatogr A. 2013;1272:15–9.

    Article  CAS  PubMed  Google Scholar 

  28. Xue C, Tian L, Xu M, Deng Z, Lin W. A new 24-membered lactone and a new polyene δ-lactone from the marine bacterium Bacillus marinus. J Antibiot. 2008;61:668–74.

    Article  CAS  Google Scholar 

  29. Marmann A, Aly A, Lin W, Wang B, Proksch P. Co-cultivation—a powerful emerging tool for enhancing the chemical diversity of microorganisms. Mar Drugs. 2014;12:1043–65.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  30. Akone SH, Mandi A, Kurtan T, et al. Inducing secondary metabolite production by the endophytic fungus Chaetomium sp. through fungal-bacterial co-culture and epigenetic modification. Tetrahedron. 2016;72:6340–7.

    Article  CAS  Google Scholar 

  31. Winnikoff JR, Glukhov E, Watrous J, Dorrestein PC, Gerwick WH. Quantitative molecular networking to profile marine cyanobacterial metabolomes. J Antibiot. 2014;67:105–12.

    Article  CAS  Google Scholar 

  32. Cabral RSA, Allard PM, Marcourt L, Young MCM, Queiroz EF, Wolfender JL. Targeted isolation of indolopyridoquinazoline alkaloids from Conchocarpus fontanesianus based on molecular networks. J Nat Prod. 2016;79:2270–8.

    Article  CAS  PubMed  Google Scholar 

  33. Böcker S, Dührkop K. Fragmentation trees reloaded. J Cheminform. 2016;8:5.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  34. Oppong-Danquah E, Parrot D, Blümel M, et al. Molecular networking-based metabolome and bioactivity analyses of marine-adapted fungi co-cultivated with phytopathogens. Front Microbiol. 2018;9:2072.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Lai Z, Tsugawa H, Wohlgemuth G, Mehta S, Mueller M, Zheng Y, et al. Identifying metabolites by integrating metabolome databases with mass spectrometry cheminformatics. Nat Methods. 2017;15:53–6.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  36. Xu XY, Shen XT, Yuan XJ, Zhou YM, Fan H, Zhu LP, et al. Metabolomics investigation of an association of induced features and corresponding fungus during the co-culture of Trametes versicolor and Ganoderma applanatum. Front Microbiol. 2018;8:2647.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Yamamoto T, Hanashima S, Haser R, et al. Role of polyol moiety of amphotericin B in ion channel formation and sterol selectivity in bilayer membrane. Bioorgan Med Chem. 2015;23:5782–8.

    Article  CAS  Google Scholar 

  38. Tareq FS, Kim JH, Lee MA, Lee HS, Lee JS, Lee YJ, et al. Antimicrobial gageomacrolactins characterized from the fermentation of the marine-derived bacterium Bacillus subtilis under optimum growth conditions. J Agric Food Chem. 2013;61:3428–34.

    Article  CAS  PubMed  Google Scholar 

  39. Zhang L, Duan J, Lv J. Phytochemistry and bioactivities of sesquiterpenoids from the Artemisia species. J Pharm Sci. 2017;26:317–34.

    CAS  Google Scholar 

  40. Wang P, Yu JH, Zhu K, Wang Y, Cheng ZQ, Jiang CS, et al. Phenolic bisabolane sesquiterpenoids from a Thai mangrove endophytic fungus, Aspergillus sp. xy02. Fitoterapia. 2018;127:322–7.

    Article  CAS  PubMed  Google Scholar 

  41. Jasmine C, Othman S, Carin L, Iurie C, Shuzhao L, Guillaume B, et al. MetaboAnalyst 4.0: towards more transparent and integrative metabolomics analysis. Nucleic Acids Res. 2018;46:486–94.

    Article  Google Scholar 

  42. Holman JD, Tabb DL, Mallick P. Employing ProteoWizard to convert raw mass spectrometry data. Curr Protoc Bioinform. 2014;46:1–9.

    Article  Google Scholar 

  43. Su G, Morris JH, Demchak B, Bader GD. Biological network exploration with Cytoscape 3. Curr Protoc Bioinform. 2014;47:1–24.

    Article  Google Scholar 

Download references


Not applicable.


This work was supported by the grants from the National Natural Science Foundation of China (31670052) and Natural Science Foundation of Liaoning Province of China (2019-ZD-0143).

Author information

Authors and Affiliations



YS and YSD conceived and designed the experiments. YS and HZZ carried out the main work, analyzed the data and drafted the manuscript. WCL, XS, HZZ, XHL, YX, KJ and ML participated in the research. YSD supervised the work and revised the manuscript. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Mei Liu or Yue-Sheng Dong.

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.

Additional figures, tables and compounds information.

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

Sun, Y., Liu, WC., Shi, X. et al. Inducing secondary metabolite production of Aspergillus sydowii through microbial co-culture with Bacillus subtilis. Microb Cell Fact 20, 42 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: