# Does Osmotic Stress Affect Natural Product Expression in Fungi?

David Overy, Hebelin Correa, Catherine C Roullier, Wei-Chiung Chi, Ka-Lai Pang, Mostafa Rateb, Rainer Ebel, Zhuo Shang, Rob Capon, Gerald Bills

► To cite this version:

David Overy, Hebelin Correa, Catherine C Roullier, Wei-Chiung Chi, Ka-Lai Pang, et al.. Does Osmotic Stress Affect Natural Product Expression in Fungi?. Marine drugs, 2017, 15 (8), pp.254. 10.3390/md15080254. [https://hal.science/hal-03015848](https://hal.science/hal-03015848)

HAL Id: hal-03015848

https://hal.science/hal-03015848

Submitted on 27 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_1.png" /></figure>

Article

# Does Osmotic Stress Affect Natural Product Expression in Fungi?

David Overy <fnref n="1" />,<fnref n="2" />,<fnref n="3" />, Hebelin Correa <fnref n="1" />,<fnref n="3" /> ID, Catherine Roullier <fnref n="4" /> ID, Wei-Chiung Chi <fnref n="5" />,<fnref n="6" />, Ka-Lai Pang <fnref n="6" />, Mostafa Rateb <fnref n="7" />,<fnref n="8" />, Rainer Ebel <fnref n="7" /> ID, Zhuo Shang <fnref n="9" />, Rob Capon <fnref n="9" />, Gerald Bills <fnref n="10" /> and Russell Kerr <fnref n="1" />,<fnref n="3" />,<fnref n="11" />,*

Received: 12 June 2017; Accepted: 8 August 2017; Published: 13 August 2017

**Abstract:** The discovery of new natural products from fungi isolated from the marine environment has increased dramatically over the last few decades, leading to the identification of over 1000 new metabolites. However, most of the reported marine-derived species appear to be terrestrial in origin yet at the same time, facultatively halo- or osmotolerant. An unanswered question regarding the apparent chemical productivity of marine-derived fungi is whether the common practice of fermenting strains in seawater contributes to enhanced secondary metabolism? To answer this question, a terrestrial isolate of *Aspergillus aculeatus* was fermented in osmotic and saline stress conditions in parallel across multiple sites. The ex-type strain of *A. aculeatus* was obtained from three different culture collections. Site-to-site variations in metabolite expression were observed, suggesting that subculturing of the same strain and subtle variations in experimental protocols can have pronounced effects upon metabolite expression. Replicated experiments at individual sites indicated that secondary metabolite production was divergent between osmotic and saline treatments. Titers of some metabolites increased or decreased in response to increasing osmolite (salt or glycerol) concentrations. Furthermore, in some cases, the expression of some secondary metabolites in relation to osmotic and saline stress was attributed to specific sources of the ex-type strains.

**Keywords:** fungi; metabolite expression; LC-MS; metabolome; osmotic stress

<footnote n="1">¹ Department of Chemistry, University of Prince Edward Island, 550 University Avenue, Charlottetown, PE C1A 4P3, Canada; david.overy@canada.ca (D.O.); hcorreav@gmail.com (H.C.)</footnote>

<footnote n="2">² Department of Pathology and Microbiology, Atlantic Veterinary College, University of PEI, 550 University Avenue, Charlottetown, PE C1A 4P3, Canada</footnote>

<footnote n="3">³ Nautilus Biosciences Canada Inc., 550 University Avenue, Charlottetown, PE C1A 4P3, Canada</footnote>

<footnote n="4">⁴ Mer Molécules Santé—EA 2160, UFR des Sciences Pharmaceutiques et Biologiques, Université de Nantes, 9 Rue Bias, 44035 Nantes, France; catherine.roullier@univ-nantes.fr</footnote>

<footnote n="5">5 Department of Food Science, National Quemoy University, Kinmen County 89250, Taiwan; joan@nqu.edu.tw</footnote>

<footnote n="6">⁶ Institute of Marine Biology and Centre of Excellence for the Oceans, National Taiwan Ocean University, 2 Pei-Ning Road, Keelung 20224, Taiwan; klpang@mail.ntou.edu.tw</footnote>

<footnote n="7">7 Marine Biodiscovery Centre, Department of Chemistry, University of Aberdeen, Meston Walk, Aberdeen AB24 3UE, Scotland, UK; m rateb@abdn.ac.uk (M.R.); r.ebel@abdn.ac.uk (R.E.)</footnote>

<footnote n="8">⁸ School of Science & Sport, University of the West of Scotland, Paisley PA1 2BE, UK</footnote>

<footnote n="9">⁹ Division of Chemistry and Structural Biology, Institute for Molecular Bioscience, The University of Queensland, 306 Carmody Road, St. Lucia QLD 4072, Australia; z3shang@ucsd.edu (Z.S.); r.capon@uq.edu.au (R.C.)</footnote>

<footnote n="10">¹⁰ Texas Therapeutics Institute, The Brown Foundation Institute of Molecular Medicine, University of Texas Health Science Center, 1881 East Rd., Houston, TX 77054, USA; Gerald.F.Bills@uth.tmc.edu</footnote>

<footnote n="11">¹¹ Department of Biomedical Sciences, Atlantic Veterinary College, University of Prince Edward Island, 550 University Avenue, Charlottetown, PE C1A 4P3, Canada</footnote>

<footnote>* Correspondence: rkerr@upei.ca; Tel.: +1-902-566-0565</footnote>
<figure><img src="image_2.png" /></figure>

Article

# Does Osmotic Stress Affect Natural Product Expression in Fungi?

David Overy <fnref n="1" />,<fnref n="2" />,<fnref n="3" />, Hebelin Correa <fnref n="1" />,<fnref n="3" /> ID, Catherine Roullier <fnref n="4" /> ID, Wei-Chiung Chi <fnref n="5" />,<fnref n="6" />, Ka-Lai Pang <fnref n="6" />, Mostafa Rateb <fnref n="7" />,<fnref n="8" />, Rainer Ebel <fnref n="7" /> ID, Zhuo Shang <fnref n="9" />, Rob Capon <fnref n="9" />, Gerald Bills <fnref n="10" /> and Russell Kerr <fnref n="1" />,<fnref n="3" />,<fnref n="11" />,*

Received: 12 June 2017; Accepted: 8 August 2017; Published: 13 August 2017

**Abstract:** The discovery of new natural products from fungi isolated from the marine environment has increased dramatically over the last few decades, leading to the identification of over 1000 new metabolites. However, most of the reported marine-derived species appear to be terrestrial in origin yet at the same time, facultatively halo- or osmotolerant. An unanswered question regarding the apparent chemical productivity of marine-derived fungi is whether the common practice of fermenting strains in seawater contributes to enhanced secondary metabolism? To answer this question, a terrestrial isolate of *Aspergillus aculeatus* was fermented in osmotic and saline stress conditions in parallel across multiple sites. The ex-type strain of *A. aculeatus* was obtained from three different culture collections. Site-to-site variations in metabolite expression were observed, suggesting that subculturing of the same strain and subtle variations in experimental protocols can have pronounced effects upon metabolite expression. Replicated experiments at individual sites indicated that secondary metabolite production was divergent between osmotic and saline treatments. Titers of some metabolites increased or decreased in response to increasing osmolite (salt or glycerol) concentrations. Furthermore, in some cases, the expression of some secondary metabolites in relation to osmotic and saline stress was attributed to specific sources of the ex-type strains.

**Keywords:** fungi; metabolite expression; LC-MS; metabolome; osmotic stress

<footnote n="1">¹ Department of Chemistry, University of Prince Edward Island, 550 University Avenue, Charlottetown, PE C1A 4P3, Canada; david.overy@canada.ca (D.O.); hcorreav@gmail.com (H.C.)</footnote>

<footnote n="2">² Department of Pathology and Microbiology, Atlantic Veterinary College, University of PEI, 550 University Avenue, Charlottetown, PE C1A 4P3, Canada</footnote>

<footnote n="3">³ Nautilus Biosciences Canada Inc., 550 University Avenue, Charlottetown, PE C1A 4P3, Canada</footnote>

<footnote n="4">⁴ Mer Molécules Santé—EA 2160, UFR des Sciences Pharmaceutiques et Biologiques, Université de Nantes, 9 Rue Bias, 44035 Nantes, France; catherine.roullier@univ-nantes.fr</footnote>

<footnote n="5">5 Department of Food Science, National Quemoy University, Kinmen County 89250, Taiwan; joan@nqu.edu.tw</footnote>

<footnote n="6">⁶ Institute of Marine Biology and Centre of Excellence for the Oceans, National Taiwan Ocean University, 2 Pei-Ning Road, Keelung 20224, Taiwan; klpang@mail.ntou.edu.tw</footnote>

<footnote n="7">7 Marine Biodiscovery Centre, Department of Chemistry, University of Aberdeen, Meston Walk, Aberdeen AB24 3UE, Scotland, UK; m rateb@abdn.ac.uk (M.R.); r.ebel@abdn.ac.uk (R.E.)</footnote>

<footnote n="8">⁸ School of Science & Sport, University of the West of Scotland, Paisley PA1 2BE, UK</footnote>

<footnote n="9">⁹ Division of Chemistry and Structural Biology, Institute for Molecular Bioscience, The University of Queensland, 306 Carmody Road, St. Lucia QLD 4072, Australia; z3shang@ucsd.edu (Z.S.); r.capon@uq.edu.au (R.C.)</footnote>

<footnote n="10">¹⁰ Texas Therapeutics Institute, The Brown Foundation Institute of Molecular Medicine, University of Texas Health Science Center, 1881 East Rd., Houston, TX 77054, USA; Gerald.F.Bills@uth.tmc.edu</footnote>

<footnote n="11">¹¹ Department of Biomedical Sciences, Atlantic Veterinary College, University of Prince Edward Island, 550 University Avenue, Charlottetown, PE C1A 4P3, Canada</footnote>

<footnote>* Correspondence: rkerr@upei.ca; Tel.: +1-902-566-0565</footnote>
# 1. Introduction

Genome sequencing confirms that fungi possess many more secondary metabolite gene clusters than are expressed in fermentations employing conventional methods and feedstocks. The impact of fermentation conditions and nutritional composition on secondary metabolite expression has been well documented and even minor variations in medium composition, dissolved oxygen, or temperature may modulate secondary metabolite biosynthesis [1,2]. Secondary metabolite-annotated gene clusters that remain transcriptionally silent even when strains are grown on a wide variety of media have been termed as "orphan", "silent", or "cryptic"; many of these have the potential to encode for new natural products [3]. Transcriptomic studies have revealed that exposure to stress can lead to transcription of "orphan" gene clusters [4]. Gene regulation studies involving halophilic fungi indicate that, although individual fungi exhibit different strategies for osmoregulation [5], differentially expressed genes encoding ion and metabolite transporters were found to be upregulated as a shared response due to salt stress (halotolerance) [6–8]. A link between secondary metabolite production as a response to salt stress is evident from comparative studies where differential expression in secondary metabolite production was observed from a strain of *Aspergillus terreus* (isolated from a saltern) and a marine-derived strain of *Mariannaea elegans* (as *Spicaria elegans*) due to increasing salinity of the cultivation medium compared to water controls [9,10].

The last two decades have seen a rapid increase in the discovery of new natural products from fungi isolated from substrata originating in marine habitats. In such reports, the natural products are often referred to as being "marine-derived". The majority of most frequently cited marine-derived taxa are associated with ruderal substratum relationships, belonging to well-known osmotolerant terrestrial genera, such as *Aspergillus* and *Penicillium* [11–13]. Inversely correlated with the increased attention given to chemistry of marine-derived fungi has been a remarkable decline in natural products exploration of terrestrial soil fungi, as they are unjustly perceived as an exhausted resource of new natural products. However, genomics studies clearly indicate this bias is unfounded [14]. These points raise the question as to whether the use of seawater and artificial seawater supplements in media formulations affects the perceived distinctiveness of the chemistry of marine-derived isolates relative to their terrestrial counterparts? Few authors describing new natural products from marine-derived fungi have questioned whether genotypically related strains of the same species from a terrestrial habitat might produce the same compounds under identical fermentation conditions [15,16].

In August 2014, a workshop held at the University of Prince Edward Island, Charlottetown, Canada brought together an international group of mycologists and marine fungal natural product chemists to discuss topics relevant to marine fungal natural product research. The MaFNaP (Marine Fungal Natural Product) consortium was consequently established and some of the initial members are co-authors of this report. The question of whether culturing a fungus in the presence of seawater has an effect upon the secondary metabolic expression of the organism was debated. We hypothesize that the addition of salt in the culture medium exerts a selective physiological bias on the metabolic output of the fungus, as stress regulatory systems and fungal development are interconnected and directly affect the regulation of secondary metabolism. We also used this as an opportunity to rigorously examine the reproducibility of the "same experiment" conducted by investigators following a defined protocol in different parts of the world. To these ends, experiments were carried out at six different sites, where ex-type strains of *Aspergillus aculeatus* (a tropical soil isolate) were purchased directly from different culture collections or shared by participating labs. The strains were cultured under a standardized set of osmotic stress conditions. The effect upon secondary metabolite output was compared using five different osmotic conditions: natural seawater equivalent to 50% and 100% seawater salinity, glycerol at a water potential equivalent to 50% and 100% seawater salinity, and a 0% osmolite control.
# 1. Introduction

Genome sequencing confirms that fungi possess many more secondary metabolite gene clusters than are expressed in fermentations employing conventional methods and feedstocks. The impact of fermentation conditions and nutritional composition on secondary metabolite expression has been well documented and even minor variations in medium composition, dissolved oxygen, or temperature may modulate secondary metabolite biosynthesis [1,2]. Secondary metabolite-annotated gene clusters that remain transcriptionally silent even when strains are grown on a wide variety of media have been termed as "orphan", "silent", or "cryptic"; many of these have the potential to encode for new natural products [3]. Transcriptomic studies have revealed that exposure to stress can lead to transcription of "orphan" gene clusters [4]. Gene regulation studies involving halophilic fungi indicate that, although individual fungi exhibit different strategies for osmoregulation [5], differentially expressed genes encoding ion and metabolite transporters were found to be upregulated as a shared response due to salt stress (halotolerance) [6–8]. A link between secondary metabolite production as a response to salt stress is evident from comparative studies where differential expression in secondary metabolite production was observed from a strain of *Aspergillus terreus* (isolated from a saltern) and a marine-derived strain of *Mariannaea elegans* (as *Spicaria elegans*) due to increasing salinity of the cultivation medium compared to water controls [9,10].

The last two decades have seen a rapid increase in the discovery of new natural products from fungi isolated from substrata originating in marine habitats. In such reports, the natural products are often referred to as being "marine-derived". The majority of most frequently cited marine-derived taxa are associated with ruderal substratum relationships, belonging to well-known osmotolerant terrestrial genera, such as *Aspergillus* and *Penicillium* [11–13]. Inversely correlated with the increased attention given to chemistry of marine-derived fungi has been a remarkable decline in natural products exploration of terrestrial soil fungi, as they are unjustly perceived as an exhausted resource of new natural products. However, genomics studies clearly indicate this bias is unfounded [14]. These points raise the question as to whether the use of seawater and artificial seawater supplements in media formulations affects the perceived distinctiveness of the chemistry of marine-derived isolates relative to their terrestrial counterparts? Few authors describing new natural products from marine-derived fungi have questioned whether genotypically related strains of the same species from a terrestrial habitat might produce the same compounds under identical fermentation conditions [15,16].

In August 2014, a workshop held at the University of Prince Edward Island, Charlottetown, Canada brought together an international group of mycologists and marine fungal natural product chemists to discuss topics relevant to marine fungal natural product research. The MaFNaP (Marine Fungal Natural Product) consortium was consequently established and some of the initial members are co-authors of this report. The question of whether culturing a fungus in the presence of seawater has an effect upon the secondary metabolic expression of the organism was debated. We hypothesize that the addition of salt in the culture medium exerts a selective physiological bias on the metabolic output of the fungus, as stress regulatory systems and fungal development are interconnected and directly affect the regulation of secondary metabolism. We also used this as an opportunity to rigorously examine the reproducibility of the "same experiment" conducted by investigators following a defined protocol in different parts of the world. To these ends, experiments were carried out at six different sites, where ex-type strains of *Aspergillus aculeatus* (a tropical soil isolate) were purchased directly from different culture collections or shared by participating labs. The strains were cultured under a standardized set of osmotic stress conditions. The effect upon secondary metabolite output was compared using five different osmotic conditions: natural seawater equivalent to 50% and 100% seawater salinity, glycerol at a water potential equivalent to 50% and 100% seawater salinity, and a 0% osmolite control.
## 2. Materials and Methods

### 2.1. *Strain and Medium Selection*

The ex-type strain of *A. aculeatus* was either obtained from three different culture collections (NRRL 20623, ATCC 16872, BCRC 32190) or distributed among the seven participating laboratories. This strain was chosen because it was easily accessible in public collections, and the draft genome sequence of this strain was available at the Joint Genome Institute [17], and a bioinformatic analysis of its genome revealed an extraordinarily rich secondary metabolism [18]. A basal medium was selected for fermentation studies that was comprised of: mannitol (40 g), maltose (40 g), yeast autolysate (Sigma Aldrich 73145, 10 g), K₂HPO₄ (2 g), MgSO₄·7H₂O (0.5 g), FeSO₄·7H₂O (0.01 g) in 1 L of distilled water. The basal medium was then used to prepare five different formulations of osmotic and seawater-enriched media designated as: 100% natural seawater (filtered seawater obtained locally), 50% natural seawater (a 1:1 dilution using deionized water), 100% glycerol (55.2 g/L of glycerol; an equivalent osmotic potential to the average salinity of seawater), 50% glycerol (27.6 g/L of glycerol), and 0% salt/glycerol (no salt/no glycerol control; basal medium formulation only). All flasks containing media were autoclaved at 121 °C for 20 min and left to cool to room temperature prior to inoculation. At each site, natural seawater was collected at the nearest access to the ocean and was stored at 4 °C until used. Seed inoculum was generated for fermentation studies for the three ex-type strains using the basal medium, incubated at 22 °C, shaking at 200 rpm for 5 days.

### 2.2. *Multi-Site Study*

For the multi-site study, a single ex-type strain subculture was used per site for fermentation studies. The ex-type strain NRRL 20623 was distributed between the sites located in Nantes (France), Aberdeen (Scotland), Charlottetown (PE, Canada), and Houston (TX, USA); while the ex-type strain ATCC 16872 was obtained by the group in Brisbane (Australia), and the ex-type strain BCRC 32190 was obtained from the Bioresource Collection and Research Center by the group in Keelung (Taiwan). At each site, 1 mL of seed inoculum was aseptically transferred into three Erlenmeyer flasks for each medium formulation (n = 3, 250 mL flasks each containing 30 mL of medium). Three uninoculated flasks were also prepared for each medium formulation as a medium "blank" (control). Inoculated flasks were incubated at 22 °C, standing, in the dark for a period of 14 days.

### 2.3. *Single-Site Study*

For the single-site study, seed inoculum was prepared for each of the three different *A. aculeatus* ex-type strains (NRRL 20623, ATCC 16872, and BCRC 32190) and 1 mL of inoculum was aseptically transferred into nine Erlenmeyer flasks for each medium formulation (n = 9, 250 mL flasks each containing 30 mL of medium). A total of five uninoculated flasks (of each fermentation medium) were also carried through the fermentation conditions as medium "blanks" to be used as controls in the data analysis. Inoculated flasks were incubated at 22 °C, standing, in the dark for a period of 14 days. Two weeks later, a second round of seed inoculum of *A. aculeatus* strain BCRC 32190 was prepared and used to inoculate nine replicate 250 mL flasks (n = 9, for each of the five fermentation media) that were incubated for 14 days (along with medium "blanks") for a data comparison to assess for experimental reproducibility between fermentation conditions.

### 2.4. *Extraction and UPLC-HRMS Analysis*

Each fermentation flask was extracted with 30 mL of ethyl acetate for 1 h, shaking gently at 150 rpm. A 12 mL aliquot of the organic phase was then removed and liquid partitioned with an equivalent volume of deionized H₂O. From the organic phase, a 5 mL aliquot was then transferred into pre-weighed glass vials and then dried under a stream of air or nitrogen.

Based on determined extract weights, samples were resuspended in methanol to yield a final concentration of 500 µg/mL. Prior to UPLC-HRMS analysis, an aliquot of internal standard
## 2. Materials and Methods

### 2.1. *Strain and Medium Selection*

The ex-type strain of *A. aculeatus* was either obtained from three different culture collections (NRRL 20623, ATCC 16872, BCRC 32190) or distributed among the seven participating laboratories. This strain was chosen because it was easily accessible in public collections, and the draft genome sequence of this strain was available at the Joint Genome Institute [17], and a bioinformatic analysis of its genome revealed an extraordinarily rich secondary metabolism [18]. A basal medium was selected for fermentation studies that was comprised of: mannitol (40 g), maltose (40 g), yeast autolysate (Sigma Aldrich 73145, 10 g), K₂HPO₄ (2 g), MgSO₄.7H₂O (0.5 g), FeSO₄.7H₂O (0.01 g) in 1 L of distilled water. The basal medium was then used to prepare five different formulations of osmotic and seawater-enriched media designated as: 100% natural seawater (filtered seawater obtained locally), 50% natural seawater (a 1:1 dilution using deionized water), 100% glycerol (55.2 g/L of glycerol; an equivalent osmotic potential to the average salinity of seawater), 50% glycerol (27.6 g/L of glycerol), and 0% salt/glycerol (no salt/no glycerol control; basal medium formulation only). All flasks containing media were autoclaved at 121 °C for 20 min and left to cool to room temperature prior to inoculation. At each site, natural seawater was collected at the nearest access to the ocean and was stored at 4 °C until used. Seed inoculum was generated for fermentation studies for the three ex-type strains using the basal medium, incubated at 22 °C, shaking at 200 rpm for 5 days.

### 2.2. *Multi-Site Study*

For the multi-site study, a single ex-type strain subculture was used per site for fermentation studies. The ex-type strain NRRL 20623 was distributed between the sites located in Nantes (France), Aberdeen (Scotland), Charlottetown (PE, Canada), and Houston (TX, USA); while the ex-type strain ATCC 16872 was obtained by the group in Brisbane (Australia), and the ex-type strain BCRC 32190 was obtained from the Bioresource Collection and Research Center by the group in Keelung (Taiwan). At each site, 1 mL of seed inoculum was aseptically transferred into three Erlenmeyer flasks for each medium formulation (n = 3, 250 mL flasks each containing 30 mL of medium). Three uninoculated flasks were also prepared for each medium formulation as a medium "blank" (control). Inoculated flasks were incubated at 22 °C, standing, in the dark for a period of 14 days.

### 2.3. *Single-Site Study*

For the single-site study, seed inoculum was prepared for each of the three different *A. aculeatus* ex-type strains (NRRL 20623, ATCC 16872, and BCRC 32190) and 1 mL of inoculum was aseptically transferred into nine Erlenmeyer flasks for each medium formulation (n = 9, 250 mL flasks each containing 30 mL of medium). A total of five uninoculated flasks (of each fermentation medium) were also carried through the fermentation conditions as medium "blanks" to be used as controls in the data analysis. Inoculated flasks were incubated at 22 °C, standing, in the dark for a period of 14 days. Two weeks later, a second round of seed inoculum of *A. aculeatus* strain BCRC 32190 was prepared and used to inoculate nine replicate 250 mL flasks (n = 9, for each of the five fermentation media) that were incubated for 14 days (along with medium "blanks") for a data comparison to assess for experimental reproducibility between fermentation conditions.

### 2.4. *Extraction and UPLC-HRMS Analysis*

Each fermentation flask was extracted with 30 mL of ethyl acetate for 1 h, shaking gently at 150 rpm. A 12 mL aliquot of the organic phase was then removed and liquid partitioned with an equivalent volume of deionized H₂O. From the organic phase, a 5 mL aliquot was then transferred into pre-weighed glass vials and then dried under a stream of air or nitrogen.

Based on determined extract weights, samples were resuspended in methanol to yield a final concentration of 500 µg/mL. Prior to UPLC-HRMS analysis, an aliquot of internal standard
(dioctyl phthalate) was added to each sample to be used as a quality control for data preprocessing. UPLC-HRMS analysis was carried out using a Kinetex 1.7 µm C₁₈ column (Phenomenex, 50 × 2.1 mm) and an Accela Thermo UPLC with a Thermo Exactive mass spectrometer (ThermoFisher Scientific, Waltham, MA, USA) in ESI+ mode (with a *m/z* 190–2000 mass range).

## 2.5. *Data Processing*

Data preprocessing was carried out using mzMine2 (Cell Unit, Okinawa Institute of Science and Technology (OIST), Onna, Okinawa, Japan). UPLC-HRMS data files were examined (including media controls and MeOH blanks along with treatment data) to determine a minimum noise level threshold. The Full Width at Half Maximum (FWHM) algorithm was used for mass detection with a noise level cut off of 2.0 E⁴. UPLC-HRMS data were deconvoluted, aligned, isotopes removed and then converted into a data matrix of discriminate variables (based on RT and *m/z*) based on peak area measurements. Peak area was normalized by dividing by the total ion current for each sample. Peak alignment was made using the Join Alignment function (with a *m/z* tolerance of 5.0 ppm and a RT tolerance of 0.1 min with a 20:10 weight for *m/z* vs. RT). A gap-filling algorithm was then applied to fill in gaps in the data matrix with detected variable values that fell below the noise limit threshold (2.0 E⁴) during mass detection (using a *m/z* tolerance of 5 ppm and RT tolerance of 0.05 min).

The data matrix was then visualized by plotting a heatmap in Excel. Redundant variables associated with medium "blank" controls and MeOH blanks were removed from the data matrix. To further visualize the data structure, the data matrix was uploaded into the R environment and the data was centered (variable fluctuations around the mean were centered to fluctuations around zero) and scaled using pareto scaling (using the square root of the standard deviation as the scaling factor). Hierarchical cluster analysis was then performed on the scaled data matrix using the Euclidean distance to compute the distance matrix and single, complete and UPGMA (Unweighted Pair Group Method with Arithmetic Mean) as the linkage criteria. Both multivariate and univariate statistical analyses were completed with the *muma* package (metabolomics univariate and multivariate analysis; [19]). Univariate statistical analysis was carried out using a decisional tree algorithm where Shapiro Wilk's test for normality was performed to determine if each variable was normally distributed or not and to decide whether a parametric (Welch's *t*-test) or a non-parametric test (Wilcoxon-Mann Whitney test) was warranted. The univariate function in *muma* then merged the *p*-values deriving from Welch's *t*-test and from Wilcoxon-Mann Whitney Test, according to the distribution of each variable (Shapiro Wilk's *p*-value > 0.05 for normal distribution, Welch's *t*-test; *p*-value < 0.05, Wilcoxon-Mann Whitney Test); in this way, for each pairwise comparison a file was created with merged *p*-values from the two tests. The univariate function also calculated log fold change increases between pairwise comparisons (between treatments, e.g., 0% control vs. 100% seawater) and carried out an ANOVA analysis for each variable.

# 3. Results

## 3.1. *Characterization of A. aculeatus* Ex-Type Strains

Robust growth was observed for each of the three *A. aculeatus* ex-type strains when exposed to increasing concentrations of glycerol and seawater salinity. All standing cultures were visually inspected and no apparent differences in colony morphology were observed between medium controls and treatment conditions. Mycelial mats extended to the walls of the flasks, covering the medium surface, with abundant sporulation. Ethyl acetate extracts of the cultures were chemically profiled using UPLC-HRMS to assess for secondary metabolite differences between culture conditions. A list of known secondary metabolites associated with *A. aculeatus* was compiled from the literature in order to associate pseudomolecular ion signals (variables) with known molecules. Signal annotation was possible for fourteen metabolites (see Figure 1) based on reported metabolite elution order in RP chromatography and associated pseudomolecular ion *m/z* data: aculene A–D, aspergillusol A, asperaculane A, aculin A, acucalbistrin A and B, acu-dioxomorpholine, CJ-15183 (and a novel
(dioctyl phthalate) was added to each sample to be used as a quality control for data preprocessing. UPLC-HRMS analysis was carried out using a Kinetex 1.7 µm C₁₈ column (Phenomenex, 50 × 2.1 mm) and an Accela Thermo UPLC with a Thermo Exactive mass spectrometer (ThermoFisher Scientific, Waltham, MA, USA) in ESI+ mode (with a *m/z* 190–2000 mass range).

## 2.5. *Data Processing*

Data preprocessing was carried out using mzMine2 (Cell Unit, Okinawa Institute of Science and Technology (OIST), Onna, Okinawa, Japan). UPLC-HRMS data files were examined (including media controls and MeOH blanks along with treatment data) to determine a minimum noise level threshold. The Full Width at Half Maximum (FWHM) algorithm was used for mass detection with a noise level cut off of 2.0 E⁴. UPLC-HRMS data were deconvoluted, aligned, isotopes removed and then converted into a data matrix of discriminate variables (based on RT and *m/z*) based on peak area measurements. Peak area was normalized by dividing by the total ion current for each sample. Peak alignment was made using the Join Alignment function (with a *m/z* tolerance of 5.0 ppm and a RT tolerance of 0.1 min with a 20:10 weight for *m/z* vs. RT). A gap-filling algorithm was then applied to fill in gaps in the data matrix with detected variable values that fell below the noise limit threshold (2.0 E⁴) during mass detection (using a *m/z* tolerance of 5 ppm and RT tolerance of 0.05 min).

The data matrix was then visualized by plotting a heatmap in Excel. Redundant variables associated with medium "blank" controls and MeOH blanks were removed from the data matrix. To further visualize the data structure, the data matrix was uploaded into the R environment and the data was centered (variable fluctuations around the mean were centered to fluctuations around zero) and scaled using pareto scaling (using the square root of the standard deviation as the scaling factor). Hierarchical cluster analysis was then performed on the scaled data matrix using the Euclidean distance to compute the distance matrix and single, complete and UPGMA (Unweighted Pair Group Method with Arithmetic Mean) as the linkage criteria. Both multivariate and univariate statistical analyses were completed with the *muma* package (metabolomics univariate and multivariate analysis; [19]). Univariate statistical analysis was carried out using a decisional tree algorithm where Shapiro Wilk's test for normality was performed to determine if each variable was normally distributed or not and to decide whether a parametric (Welch's *t*-test) or a non-parametric test (Wilcoxon-Mann Whitney test) was warranted. The univariate function in *muma* then merged the *p*-values deriving from Welch's *t*-test and from Wilcoxon-Mann Whitney Test, according to the distribution of each variable (Shapiro Wilk's *p*-value > 0.05 for normal distribution, Welch's *t*-test; *p*-value < 0.05, Wilcoxon-Mann Whitney Test); in this way, for each pairwise comparison a file was created with merged *p*-values from the two tests. The univariate function also calculated log fold change increases between pairwise comparisons (between treatments, e.g., 0% control vs. 100% seawater) and carried out an ANOVA analysis for each variable.

# 3. Results

## 3.1. *Characterization of A. aculeatus* Ex-Type Strains

Robust growth was observed for each of the three *A. aculeatus* ex-type strains when exposed to increasing concentrations of glycerol and seawater salinity. All standing cultures were visually inspected and no apparent differences in colony morphology were observed between medium controls and treatment conditions. Mycelial mats extended to the walls of the flasks, covering the medium surface, with abundant sporulation. Ethyl acetate extracts of the cultures were chemically profiled using UPLC-HRMS to assess for secondary metabolite differences between culture conditions. A list of known secondary metabolites associated with *A. aculeatus* was compiled from the literature in order to associate pseudomolecular ion signals (variables) with known molecules. Signal annotation was possible for fourteen metabolites (see Figure 1) based on reported metabolite elution order in RP chromatography and associated pseudomolecular ion *m/z* data: aculene A–D, aspergillusol A, asperaculane A, aculin A, acucalbistrin A and B, acu-dioxomorpholine, CJ-15183 (and a novel
dihydro-analog), and secalonic acids B, D and F [20–26]. Assignment of the [M + H]$^+$ pseudomolecular ion was supported by observed adduct ions and known fragment ions sharing exact chromatographic peak shapes and retention times using single ion monitoring. Pseudomolecular ions associated with four additional metabolites (asperaculane B, okaramines A, H and J) were also observed; however, abundance levels for these ions fell below the threshold cut-off for peak assignment during data preprocessing and therefore were not included in the metabolomic analysis of the data set.

<figure><img src="image_2.png" /><figcaption>**Figure 1.** Annotated total ion current chromatogram of *A. aculeatus* ex-type strain ATCC 16872 (0% alt/glycerol control).</figcaption></figure>

## 3.2. Multi-Site Study: Summary

A “first pass” cluster analysis was performed on the multi-site data set (using both single linkage and complete linkage criteria) to inspect the “structure” of the data. Looking at the overall data structure, obvious structure relating to a consistent secondary metabolism response to seawater or glycerol treatments was absent (Figure 2). Interestingly, the data structure modeled differences between sites as evident from the formation of site-specific clades (whereas if the data structure were to indicate differences between treatments, clades would consist of treatments representing multiple sites). To account for contaminant ions (MS contaminants) commonly detected by mass spectrometry (pseudomolecular ions associated with plasticizers and detergents), the variable peak list was cross-referenced against a published list of common MS contaminants [27]. Site-specific
dihydro-analog), and secalonic acids B, D and F [20–26]. Assignment of the [M + H]$^+$ pseudomolecular ion was supported by observed adduct ions and known fragment ions sharing exact chromatographic peak shapes and retention times using single ion monitoring. Pseudomolecular ions associated with four additional metabolites (asperaculane B, okaramines A, H and J) were also observed; however, abundance levels for these ions fell below the threshold cut-off for peak assignment during data preprocessing and therefore were not included in the metabolomic analysis of the data set.

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Annotated total ion current chromatogram of *A. aculeatus* ex-type strain ATCC 16872 (0% alt/glycerol control).</figcaption></figure>

## 3.2. Multi-Site Study: Summary

A “first pass” cluster analysis was performed on the multi-site data set (using both single linkage and complete linkage criteria) to inspect the “structure” of the data. Looking at the overall data structure, obvious structure relating to a consistent secondary metabolism response to seawater or glycerol treatments was absent (Figure 2). Interestingly, the data structure modeled differences between sites as evident from the formation of site-specific clades (whereas if the data structure were to indicate differences between treatments, clades would consist of treatments representing multiple sites). To account for contaminant ions (MS contaminants) commonly detected by mass spectrometry (pseudomolecular ions associated with plasticizers and detergents), the variable peak list was cross-referenced against a published list of common MS contaminants [27]. Site-specific
differences in the presence of MS contaminants were observed (Figure S1) and these variables were subsequently removed from the data set. However, following removal of the MS contaminant variables, the site-specific clustering observed in Figure 2 remained. To reduce the complexity of the data, cluster analysis was performed using data sets consisting of only one treatment. Cluster analysis results were consistent between each of the individual treatment data sets. Triplicates grouped together into individual clades demonstrating reproducibility between replicates and that the observed differences in metabolism were site based (e.g., Figure 3: dendrogram of the 0% osmolite control samples; see Figure S2 for representative UPLC-HRMS chromatographic profiles). Distinct clades were also observed representing the culture collection from which the *A. aculeatus* ex-type strain was obtained indicating that differences in metabolite expression could not be contributed to site specific differences alone, but also relate to the source of the ex-type strain (ATCC 16872 was obtained by the Australian group, BCRC 32190 was obtained by the Taiwanese group, and NRRL 20623 was shared between the Scottish, Canadian, French and American groups). Follow-up correspondence between the six collaborating research groups led to the discovery that slight deviations in the experimental protocol were implemented at the different research sites (e.g., flasks were covered by cotton bungs, milk filters, and screw caps that would affect air exchange in the culture flasks; and different sources of natural seawater were used).

<figure><img src="image_3.png" /><figcaption>**Figure 2.** Single linkage cluster analysis of the complete *A. aculeatus* data set (all sites, all treatments).</figcaption></figure>
differences in the presence of MS contaminants were observed (Figure S1) and these variables were subsequently removed from the data set. However, following removal of the MS contaminant variables, the site-specific clustering observed in Figure 2 remained. To reduce the complexity of the data, cluster analysis was performed using data sets consisting of only one treatment. Cluster analysis results were consistent between each of the individual treatment data sets. Triplicates grouped together into individual clades demonstrating reproducibility between replicates and that the observed differences in metabolism were site based (e.g., Figure 3: dendrogram of the 0% osmolite control samples; see Figure S2 for representative UPLC-HRMS chromatographic profiles). Distinct clades were also observed representing the culture collection from which the *A. aculeatus* ex-type strain was obtained indicating that differences in metabolite expression could not be contributed to site specific differences alone, but also relate to the source of the ex-type strain (ATCC 16872 was obtained by the Australian group, BCRC 32190 was obtained by the Taiwanese group, and NRRL 20623 was shared between the Scottish, Canadian, French and American groups). Follow-up correspondence between the six collaborating research groups led to the discovery that slight deviations in the experimental protocol were implemented at the different research sites (e.g., flasks were covered by cotton bungs, milk filters, and screw caps that would affect air exchange in the culture flasks; and different sources of natural seawater were used).

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Single linkage cluster analysis of the complete *A. aculeatus* data set (all sites, all treatments).</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Single linkage cluster analysis of the *A. aculeatus* 0% osmolite control treatment (all sites, single treatment).</figcaption></figure>

3.3. *Single-Site Study: Summary*

To compensate for site-to-site variability, the experiment was repeated at a single site and this time
nine rather than three replicates were used with each of the three different *A. aculeatus* ex-type strains
(ATCC 16872, BCRC 32190, and NRRL 20623). To visualize associated relationships between ex-type
strains and culture conditions, hierarchical cluster analysis was performed on the single-site study
data set and the resulting dendrogram using the pareto scaling is presented in Figure 4A. In Figure 4,
samples were highlighted by color to assess the effect that strain differences had upon the data model:
green representing BCRC 32190, red representing NRRL 20623, and black representing ATCC 16872.
Strain associated sample coloration within the dendrogram highlighted distinct clades containing
samples of the same strain associated with different treatments (clade nodes are highlighted by arrows
colored with respect to sample strain, Figure 4A). However several clades of the dendrogram consisted
of samples originating from different stains sharing the same fermentation treatments, indicating
that data structure alone cannot be explained by strain differences. Principal component analysis
was carried out on the pareto scaled data matrix (Figure 4B). The first two principal components
explained a combined 57% of the variance in the data model (with approximately 30% of the variance
explained by PC1 and 26% of the variance explained by PC2). NRRL 20623 samples associated
with different treatments (0% seawater/glycerol control, 50% and 100% seawater) separated from
the remaining samples along PC1 indicating differences in secondary metabolite expression of these
samples compared to the ATCC 16872 and BCRC 32190 ex-type strains. The PCA data supported the
observed formation of a distinct clade within the cluster analysis (red arrow in Figure 4A), confirming
that some of the variance in the data model can be explained by differences in secondary metabolite
production between the ex-type strains. The grouping of the remaining ATCC 16872 samples (50%
glycerol and 100% glycerol) with similar treatment conditions for the other two ex-type strains in
the PCA (similar to the cluster analysis results), reinforced the conclusion that the variance in the
data model could not be explicitly attributed to differences in secondary metabolite production
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Single linkage cluster analysis of the *A. aculeatus* 0% osmolite control treatment (all sites, single treatment).</figcaption></figure>

3.3. *Single-Site Study: Summary*

To compensate for site-to-site variability, the experiment was repeated at a single site and this time
nine rather than three replicates were used with each of the three different *A. aculeatus* ex-type strains
(ATCC 16872, BCRC 32190, and NRRL 20623). To visualize associated relationships between ex-type
strains and culture conditions, hierarchical cluster analysis was performed on the single-site study
data set and the resulting dendrogram using the pareto scaling is presented in Figure 4A. In Figure 4,
samples were highlighted by color to assess the effect that strain differences had upon the data model:
green representing BCRC 32190, red representing NRRL 20623, and black representing ATCC 16872.
Strain associated sample coloration within the dendrogram highlighted distinct clades containing
samples of the same strain associated with different treatments (clade nodes are highlighted by arrows
colored with respect to sample strain, Figure 4A). However several clades of the dendrogram consisted
of samples originating from different stains sharing the same fermentation treatments, indicating
that data structure alone cannot be explained by strain differences. Principal component analysis
was carried out on the pareto scaled data matrix (Figure 4B). The first two principal components
explained a combined 57% of the variance in the data model (with approximately 30% of the variance
explained by PC1 and 26% of the variance explained by PC2). NRRL 20623 samples associated
with different treatments (0% seawater/glycerol control, 50% and 100% seawater) separated from
the remaining samples along PC1 indicating differences in secondary metabolite expression of these
samples compared to the ATCC 16872 and BCRC 32190 ex-type strains. The PCA data supported the
observed formation of a distinct clade within the cluster analysis (red arrow in Figure 4A), confirming
that some of the variance in the data model can be explained by differences in secondary metabolite
production between the ex-type strains. The grouping of the remaining ATCC 16872 samples (50%
glycerol and 100% glycerol) with similar treatment conditions for the other two ex-type strains in
the PCA (similar to the cluster analysis results), reinforced the conclusion that the variance in the
data model could not be explicitly attributed to differences in secondary metabolite production
between the three ex-type strains; rather media composition (treatment) also influenced secondary
metabolite production.

<figure><img src="image_5.png" /><figcaption>**Figure 4.** (A) Cluster analysis dendrogram of complete data set (pareto scaling, Eucledian distance function, UPGMA linkage criterion). Arrows indicate relevant nodes and coloration denote ex-type strain (black for ATCC 16872, green for BCRC 32190, red for NRRL 20623); (B) Principal component analysis score plot for PC1 vs. PC2 (PCA of complete data set).</figcaption></figure>

To visualize the effect that medium composition played upon the single-site hierarchical cluster
analysis result in Figure 4, samples were colored based on treatment: black/grey for 0% control (0%
salt/glycerol), red for 50% glycerol, green for 100% glycerol, dark blue for 50% seawater, and light
blue for 100% seawater (as depicted in Figure 5). Two distinct clades based on seawater vs. glycerol
treatment were observed; however, 0% control samples also formed individual subclades based on
ex-type strain (black arrows in Figure 5). Within the large salt and glycerol treatment clades, samples
intermingled between gradations of seawater or glycerol (i.e., samples within the salt clade do not
between the three ex-type strains; rather media composition (treatment) also influenced secondary
metabolite production.

<figure><img src="image_6.png" /><figcaption>**Figure 4.** (A) Cluster analysis dendrogram of complete data set (pareto scaling, Eucledian distance function, UPGMA linkage criterion). Arrows indicate relevant nodes and coloration denote ex-type strain (black for ATCC 16872, green for BCRC 32190, red for NRRL 20623); (B) Principal component analysis score plot for PC1 vs. PC2 (PCA of complete data set).</figcaption></figure>

To visualize the effect that medium composition played upon the single-site hierarchical cluster
analysis result in Figure 4, samples were colored based on treatment: black/grey for 0% control (0%
salt/glycerol), red for 50% glycerol, green for 100% glycerol, dark blue for 50% seawater, and light
blue for 100% seawater (as depicted in Figure 5). Two distinct clades based on seawater vs. glycerol
treatment were observed; however, 0% control samples also formed individual subclades based on
ex-type strain (black arrows in Figure 5). Within the large salt and glycerol treatment clades, samples
intermingled between gradations of seawater or glycerol (i.e., samples within the salt clade do not
group into uniform subclades based on 50% salt and 100% salt) due to the influence of secondary
metabolite production differences between ex-type strains. Therefore, to further explore differences
in secondary metabolite expression due to different culture conditions, data sets of the three ex-type
strains were independently analyzed.

<figure><img src="image_6.png" /><figcaption>**Figure 5.** (A) Cluster analysis dendrogram of complete data set (pareto scaling, Euclidian distance function, UPGMA linkage criterion). Black arrows indicate clustering of ex-type 0% control treatments. Coloration highlights denote treatment (black: 0% control; red: 50% glycerol; green: 100% glycerol; dark blue: 50% salt; light blue: 100% salt); (B) Principal component analysis score plot for PC1 vs. PC2 (PCA of complete data set).</figcaption></figure>

Principal component analysis of the individual ex-type strains (Figure 6A–C) yielded models with a greater cumulative explained variance for PC1 and 2 (63% for ATCC 16872, 73.4% for NRRL 20623, and 58.5% for BCRC 32190) compared to the previous PCA generated from the combined ex-type strain dataset (Figure 4B). Similar trends were observed in PCA score plots when comparing ex-type strains ATCC 16872 and NRRL 20623 (Figure 6A,C), where salt treatment samples had a positive directional movement along PC1 and glycerol treatments a negative directional movement. In both PCA models of ex-type strains ATCC 16872 and NRRL 20623, 50% salt and 100% salt treatments separate from each
group into uniform subclades based on 50% salt and 100% salt) due to the influence of secondary
metabolite production differences between ex-type strains. Therefore, to further explore differences
in secondary metabolite expression due to different culture conditions, data sets of the three ex-type
strains were independently analyzed.

<figure><img src="image_7.png" /><figcaption>**Figure 5.** (A) Cluster analysis dendrogram of complete data set (pareto scaling, Euclidian distance function, UPGMA linkage criterion). Black arrows indicate clustering of ex-type 0% control treatments. Coloration highlights denote treatment (black: 0% control; red: 50% glycerol; green: 100% glycerol; dark blue: 50% salt; light blue: 100% salt); (B) Principal component analysis score plot for PC1 vs. PC2 (PCA of complete data set).</figcaption></figure>

Principal component analysis of the individual ex-type strains (Figure 6A–C) yielded models with a greater cumulative explained variance for PC1 and 2 (63% for ATCC 16872, 73.4% for NRRL 20623, and 58.5% for BCRC 32190) compared to the previous PCA generated from the combined ex-type strain dataset (Figure 4B). Similar trends were observed in PCA score plots when comparing ex-type strains ATCC 16872 and NRRL 20623 (Figure 6A,C), where salt treatment samples had a positive directional movement along PC1 and glycerol treatments a negative directional movement. In both PCA models of ex-type strains ATCC 16872 and NRRL 20623, 50% salt and 100% salt treatments separate from each
other along PC2. In comparison, in the PCA of the BCRC 32190 samples (Figure 6B) both glycerol treatments group together with those of the 0% control samples, with no movement along PC1 and only slight negative movement along PC2. However, saline treatment samples separate from the other treatments (in this case with positive movement along PC2) and from each other (due to negative movement along PC1 for 100% salt and positive movement along PC1 for 50% salt).

<figure><img src="image_7.png" /><figcaption>**Figure 6.** Principal component analysis score plots of PC1 and 2 from ex-type strain data sets: (A) ATCC 16872; (B) BCRC 32190; (C) NRRL 20623.</figcaption></figure>

### 3.4. Single-Site Study: Effects of Osmolite Concentration upon Secondary Metabolism

Univariate data analyses were carried out to compare the differences in metabolite expression both between treatments for a given ex-type strain as well as between ex-type strains. Significant variables (*p*-value < 0.05) in pairwise comparisons between treatment vs. 0% control for individual ex-type strains were determined and associated log fold changes were calculated (see Table S1). A short list of variables associated with significant fold increases in treatment samples vs. control samples were
other along PC2. In comparison, in the PCA of the BCRC 32190 samples (Figure 6B) both glycerol treatments group together with those of the 0% control samples, with no movement along PC1 and only slight negative movement along PC2. However, saline treatment samples separate from the other treatments (in this case with positive movement along PC2) and from each other (due to negative movement along PC1 for 100% salt and positive movement along PC1 for 50% salt).

<figure><img src="image_8.png" /><figcaption>**Figure 6.** Principal component analysis score plots of PC1 and 2 from ex-type strain data sets: (A) ATCC 16872; (B) BCRC 32190; (C) NRRL 20623.</figcaption></figure>

### 3.4. Single-Site Study: Effects of Osmolite Concentration upon Secondary Metabolism

Univariate data analyses were carried out to compare the differences in metabolite expression both between treatments for a given ex-type strain as well as between ex-type strains. Significant variables (*p*-value < 0.05) in pairwise comparisons between treatment vs. 0% control for individual ex-type strains were determined and associated log fold changes were calculated (see Table S1). A short list of variables associated with significant fold increases in treatment samples vs. control samples were
compiled to reduce the dimensionality of the data set and highlight those variables that contributed the most to the variance in the data model. Significant variables associated with relevant PCA loadings values contributing the most to the variance in the data model, driving sample separation along PC1 and PC2 (based on treatment), were tabulated for each of the three ex-type strains (Figure 7; purple highlights indicate positive loadings values and yellow indicate negative loadings). Trends of significant increased metabolite expression were observed; with particular metabolites have an association with an increase in osmotic potential (red and green highlights; 50% glycerol and 100% glycerol respectively) while others are associated with an increase in salinity (dark and light blue highlights; 50% salt and 100% salt respectively).

<table><caption>**Figure 7.** Summary of relevant variables having significant increases in treatment vs. control pairwise comparisons for all three *A. aculeatus* ex-type strains when cultured at a single site (single site comparison). Highlighted cells represent statistically significant (*p*-value < 0.05) log fold increases in peak area (red: 50% glycerol; green: 100% glycerol; dark blue: 50% salt; light blue: 100% salt). Corresponding PCA loading values for PC1 and PC2 associated with modeling of individual ex-type strain data sets (purple represents positive and yellow represents negative loadings, color intensity based on greatest loading value per PC).</caption><thead><tr><th rowspan="3">Variable (RT_m/z)</th><th rowspan="3">Pseudomolecular Ion Annotation</th><th colspan="4">ATCC 16872</th><th colspan="4">BCRC 32190</th><th colspan="4">NRRL 20623</th><th colspan="2">ATCC 16872</th><th colspan="2">BCRC 32190</th><th colspan="2">NRRL 20623</th></tr><tr><th colspan="3">Glycerol Treatment</th><th rowspan="2">Glycerol Treatment</th><th colspan="3">Glycerol Treatment</th><th rowspan="2">Glycerol Treatment</th><th colspan="3">Glycerol Treatment</th><th rowspan="2">PCA PC1</th><th rowspan="2">PCA PC2</th><th rowspan="2">PCA PC1</th><th rowspan="2">PCA PC2</th><th rowspan="2">PCA PC1</th><th rowspan="2">PCA PC2</th></tr><tr><th>50%G</th><th>100%G</th><th>Salt Treatment</th><th>50%G</th><th>100%G</th><th>Salt Treatment</th><th>50%G</th><th>100%G</th><th>Salt Treatment</th></tr></thead><tbody><tr><td>RT2.45_201.127732</td><td>Aculene A [fragment]+</td><td>0.36</td><td>0.14</td><td>-0.32</td><td>-0.17</td><td>0.08</td><td>0.18</td><td>-0.18</td><td>0.03</td><td>0.16</td><td>0.32</td><td>-0.24</td><td>-0.11</td><td>-0.25</td><td>0.01</td><td>-0.16</td><td>-0.27</td><td>-0.22</td><td>0.13</td></tr><tr><td>RT2.45_316.190751</td><td>Aculene A [M+H]+</td><td>0.16</td><td>0.14</td><td>-0.38</td><td>-0.19</td><td>0.06</td><td>0.16</td><td>-0.21</td><td>0.02</td><td>0.16</td><td>0.30</td><td>-0.33</td><td>-0.17</td><td>-0.31</td><td>0.02</td><td>-0.20</td><td>-0.33</td><td>-0.28</td><td>0.18</td></tr><tr><td>RT2.55_318.206391</td><td>Aculene B [M+H]+</td><td>0.10</td><td>0.10</td><td>-0.57</td><td>-0.28</td><td>-0.09</td><td>0.01</td><td>-0.42</td><td>-0.09</td><td>0.07</td><td>0.20</td><td>-0.38</td><td>-0.26</td><td>-0.20</td><td>0.04</td><td>-0.14</td><td>-0.20</td><td>-0.14</td><td>0.13</td></tr><tr><td>RT2.55_477.150017</td><td>Aspergillus [M+H]+</td><td>0.09</td><td>0.15</td><td>0.64</td><td>0.60</td><td>0.11</td><td>0.13</td><td>0.25</td><td>0.25</td><td>0.11</td><td>0.25</td><td>0.25</td><td>0.34</td><td>0.31</td><td>0.10</td><td>0.25</td><td>0.25</td><td>0.03</td><td>0.27</td></tr><tr><td>RT2.55_499.132347</td><td>Aspergillus [M+Na]+</td><td>0.11</td><td>0.13</td><td>0.60</td><td>0.75</td><td>0.11</td><td>0.13</td><td>0.25</td><td>0.24</td><td>0.11</td><td>0.25</td><td>0.25</td><td>0.33</td><td>0.31</td><td>0.11</td><td>0.25</td><td>0.25</td><td>0.04</td><td>0.28</td></tr><tr><td>RT2.68_291.123173</td><td>Unknown A</td><td>-0.48</td><td>0.30</td><td>-3.18</td><td>-1.48</td><td>-0.77</td><td>0.09</td><td>0.52</td><td>1.70</td><td>-0.25</td><td>0.18</td><td>-1.40</td><td>-0.76</td><td>-0.08</td><td>0.06</td><td>-0.14</td><td>0.06</td><td>-0.02</td><td>0.04</td></tr><tr><td>RT2.71_201.127768</td><td>Aculene C [M+H-H2O]+</td><td>-0.24</td><td>-0.27</td><td>0.05</td><td>0.25</td><td>-0.11</td><td>-0.08</td><td>-0.10</td><td>0.03</td><td>-0.18</td><td>-0.67</td><td>0.03</td><td>0.16</td><td>0.21</td><td>0.01</td><td>-0.08</td><td>0.10</td><td>0.10</td><td>0.10</td></tr><tr><td>RT2.75_241.119986</td><td>Aculene C [M+Na]+</td><td>-0.27</td><td>-0.32</td><td>0.10</td><td>0.28</td><td>-0.10</td><td>-0.08</td><td>-0.09</td><td>0.03</td><td>-0.19</td><td>-0.75</td><td>0.04</td><td>0.16</td><td>0.21</td><td>0.25</td><td>-0.04</td><td>0.10</td><td>0.36</td><td>0.07</td></tr><tr><td>RT2.75_345.131426</td><td>Unknown C</td><td>-1.23</td><td>0.11</td><td>-2.28</td><td>-1.56</td><td>-0.93</td><td>-0.01</td><td>0.34</td><td>1.59</td><td>0.18</td><td>-1.06</td><td>-1.05</td><td>-0.86</td><td>-0.05</td><td>0.05</td><td>-0.14</td><td>0.04</td><td>-0.01</td><td>0.04</td></tr><tr><td>RT2.75_459.2515</td><td>Aculene C [2M+Na]+</td><td>-0.40</td><td>-0.42</td><td>0.03</td><td>0.34</td><td>-0.18</td><td>-0.14</td><td>-0.23</td><td>-0.06</td><td>-0.23</td><td>-1.06</td><td>-0.12</td><td>0.13</td><td>0.08</td><td>0.15</td><td>-0.01</td><td>0.02</td><td>0.17</td><td>0.10</td></tr><tr><td>RT2.77_219.138014</td><td>Aculene C [M+H]+</td><td>-0.33</td><td>-0.34</td><td>0.02</td><td>0.28</td><td>-0.18</td><td>-0.14</td><td>-0.21</td><td>-0.05</td><td>-0.21</td><td>-0.88</td><td>-0.05</td><td>0.13</td><td>0.08</td><td>0.14</td><td>-0.02</td><td>0.02</td><td>0.17</td><td>0.08</td></tr><tr><td>RT3.01_203.142363</td><td>Aculene D [M+H-H2O]+</td><td>-0.34</td><td>-0.39</td><td>-0.16</td><td>0.14</td><td>-0.31</td><td>-0.25</td><td>-0.39</td><td>-0.19</td><td>-0.31</td><td>-0.88</td><td>-0.03</td><td>0.06</td><td>0.07</td><td>0.16</td><td>-0.04</td><td>-0.03</td><td>0.19</td><td>0.05</td></tr><tr><td>RT3.03_201.127843</td><td>Aculene analog [frag]+</td><td>-0.21</td><td>-0.04</td><td>-0.08</td><td>0.05</td><td>-0.12</td><td>-0.09</td><td>-0.13</td><td>0.12</td><td>-0.12</td><td>-0.46</td><td>0.02</td><td>0.10</td><td>0.03</td><td>0.17</td><td>-0.12</td><td>0.07</td><td>0.19</td><td>0.08</td></tr><tr><td>RT3.17_470.217443</td><td>Unknown B [M+H]+</td><td>-0.41</td><td>0.16</td><td>-1.22</td><td>-1.19</td><td>-0.55</td><td>-0.44</td><td>-0.47</td><td>0.42</td><td>0.22</td><td>0.03</td><td>-0.89</td><td>-0.78</td><td>-0.19</td><td>0.02</td><td>-0.23</td><td>0.15</td><td>-0.12</td><td>0.25</td></tr><tr><td>RT3.18_492.200173</td><td>Unknown B [M+Na]+</td><td>-0.44</td><td>0.13</td><td>-1.11</td><td>-1.06</td><td>-0.43</td><td>-0.37</td><td>-0.37</td><td>0.39</td><td>0.25</td><td>0.02</td><td>-0.57</td><td>-0.55</td><td>-0.12</td><td>0.02</td><td>-0.14</td><td>0.09</td><td>-0.07</td><td>0.16</td></tr><tr><td>RT3.79_661.153413</td><td>Secalonic acid [D7] [M+Na]+</td><td>-0.03</td><td>-0.16</td><td>0.09</td><td>0.13</td><td>0.00</td><td>-0.03</td><td>0.04</td><td>0.04</td><td>0.06</td><td>0.03</td><td>0.03</td><td>0.00</td><td>0.15</td><td>0.07</td><td>0.01</td><td>0.14</td><td>0.01</td><td>0.25</td></tr><tr><td>RT3.8_639.171562</td><td>Secalonic acid [D7] [M+H]+</td><td>-0.05</td><td>-0.15</td><td>0.00</td><td>0.05</td><td>-0.02</td><td>-0.06</td><td>-0.02</td><td>-0.01</td><td>0.05</td><td>0.05</td><td>0.03</td><td>0.04</td><td>0.16</td><td>0.18</td><td>-0.03</td><td>0.10</td><td>-0.05</td><td>0.08</td></tr><tr><td>RT3.92_600.264904</td><td>CJ-15183 [M+NH4]+</td><td>0.47</td><td>0.03</td><td>-0.33</td><td>-0.68</td><td>0.52</td><td>0.42</td><td>0.03</td><td>-0.33</td><td>1.00</td><td>1.45</td><td>0.22</td><td>-0.27</td><td>-0.07</td><td>-0.09</td><td>0.03</td><td>-0.10</td><td>-0.14</td><td>-0.08</td></tr><tr><td>RT3.92_605.219954</td><td>CJ-15183 [M+Na]+</td><td>0.44</td><td>0.03</td><td>-0.32</td><td>-0.66</td><td>0.52</td><td>0.42</td><td>0.04</td><td>-0.30</td><td>0.94</td><td>1.35</td><td>0.17</td><td>-0.26</td><td>-0.07</td><td>-0.10</td><td>0.03</td><td>-0.11</td><td>-0.15</td><td>-0.09</td></tr><tr><td>RT4.15_602.281542</td><td>Dihydro-CJ-15183 [M+NH4]+</td><td>0.26</td><td>-0.21</td><td>-0.62</td><td>-0.96</td><td>0.26</td><td>0.16</td><td>-0.25</td><td>-0.67</td><td>0.84</td><td>1.05</td><td>0.07</td><td>-0.54</td><td>-0.12</td><td>-0.13</td><td>0.04</td><td>-0.19</td><td>-0.10</td><td>-0.10</td></tr><tr><td>RT4.15_607.236763</td><td>Dihydro-CJ-15183 [M+Na]+</td><td>0.26</td><td>-0.20</td><td>-0.60</td><td>-0.98</td><td>0.26</td><td>0.16</td><td>-0.24</td><td>-0.65</td><td>0.84</td><td>1.05</td><td>0.09</td><td>-0.51</td><td>-0.12</td><td>-0.14</td><td>0.05</td><td>-0.20</td><td>-0.20</td><td>-0.11</td></tr><tr><td>RT4.28_767.393453</td><td>Unknown D [M+H]+</td><td>-0.22</td><td>0.13</td><td>-0.84</td><td>-0.68</td><td>-0.41</td><td>-0.27</td><td>-0.43</td><td>0.24</td><td>0.04</td><td>-0.34</td><td>-0.81</td><td>-0.45</td><td>-0.08</td><td>-0.01</td><td>-0.09</td><td>0.04</td><td>-0.03</td><td>0.16</td></tr><tr><td>RT4.48_349.154545</td><td>Acu-dioxomorpholine [frag]+</td><td>-0.01</td><td>-0.08</td><td>-0.23</td><td>-0.23</td><td>0.02</td><td>-0.03</td><td>0.03</td><td>-0.71</td><td>-0.11</td><td>-0.19</td><td>-0.04</td><td>-0.40</td><td>0.06</td><td>-0.20</td><td>0.17</td><td>-0.07</td><td>0.01</td><td>-0.06</td></tr><tr><td>RT4.49_361.154289</td><td>Acu-dioxomorpholine [frag]+</td><td>-0.01</td><td>-0.08</td><td>-0.22</td><td>-0.22</td><td>0.02</td><td>-0.03</td><td>0.03</td><td>-0.73</td><td>-0.10</td><td>-0.19</td><td>-0.03</td><td>-0.42</td><td>0.08</td><td>-0.25</td><td>0.22</td><td>-0.08</td><td>0.01</td><td>-0.08</td></tr><tr><td>RT4.5_417.216965</td><td>Acu-dioxomorpholine [M+H]+</td><td>0.01</td><td>-0.08</td><td>-0.23</td><td>-0.03</td><td>0.02</td><td>-0.03</td><td>0.03</td><td>-0.72</td><td>-0.10</td><td>-0.17</td><td>-0.03</td><td>-0.41</td><td>0.12</td><td>-0.43</td><td>-0.15</td><td>-0.15</td><td>0.01</td><td>-0.13</td></tr></tbody></table>

Aculene A (1) and B (2) (Scheme 1) both increased in abundance in response to increased glycerol concentrations and decreased when exposed to an increase in salinity (Figures S3, S5 and S6). Pseudomolecular ions associated with these metabolites had high loading values associated with the separation of osmotic stress treatment samples away from saline treatment samples in PCA models for all three ex-type strains (Figures 6 A–C and 7). The opposite however was observed for aculene C (3) (a presumed biosynthetic precursor to 1 and 2) and a second metabolite (sharing a characteristic fragment ion *m/z* 201.1278 for the aculenes; Figures S3, S5 and S6), which exhibited significant increases in 100% seawater salt treatments for ex-type strains ATCC 16872 and NRRL 20623 and had high loadings values associated with the separation of 100% seawater treatment samples in the PCA models (Figures 6A,B and 7). Decrease in aculene A (1) and B (2) production in saline conditions may be related to the increase in aculene C (3) (which was more pronounced in the ATCC 16872 and NRRL 20623 ex-type strains).
compiled to reduce the dimensionality of the data set and highlight those variables that contributed the most to the variance in the data model. Significant variables associated with relevant PCA loadings values contributing the most to the variance in the data model, driving sample separation along PC1 and PC2 (based on treatment), were tabulated for each of the three ex-type strains (Figure 7; purple highlights indicate positive loadings values and yellow indicate negative loadings). Trends of significant increased metabolite expression were observed; with particular metabolites have an association with an increase in osmotic potential (red and green highlights; 50% glycerol and 100% glycerol respectively) while others are associated with an increase in salinity (dark and light blue highlights; 50% salt and 100% salt respectively).

<table><caption>**Figure 7.** Summary of relevant variables having significant increases in treatment vs. control pairwise comparisons for all three *A. aculeatus* ex-type strains when cultured at a single site (single site comparison). Highlighted cells represent statistically significant (*p*-value < 0.05) log fold increases in peak area (red: 50% glycerol; green: 100% glycerol; dark blue: 50% salt; light blue: 100% salt). Corresponding PCA loading values for PC1 and PC2 associated with modeling of individual ex-type strain data sets (purple represents positive and yellow represents negative loadings, color intensity based on greatest loading value per PC).</caption><thead><tr><th rowspan="3">Variable (RT_m/z)</th><th rowspan="3">Pseudomolecular Ion Annotation</th><th colspan="4">ATCC 16872</th><th colspan="4">BCRC 32190</th><th colspan="4">NRRL 20623</th><th colspan="2">ATCC 16872</th><th colspan="2">BCRC 32190</th><th colspan="2">NRRL 20623</th></tr><tr><th colspan="3">Glycerol Treatment</th><th rowspan="2">Glycerol Treatment</th><th colspan="3">Glycerol Treatment</th><th rowspan="2">Glycerol Treatment</th><th colspan="3">Glycerol Treatment</th><th rowspan="2">PCA Loadings</th><th rowspan="2">PCA Loadings</th><th rowspan="2">PCA Loadings</th><th rowspan="2">PCA Loadings</th></tr><tr><th>50%G</th><th>100%G</th><th>Salt Treatment</th><th>50%G</th><th>100%G</th><th>Salt Treatment</th><th>50%G</th><th>100%G</th><th>Salt Treatment</th><th>50%G</th><th>100%G</th><th>Salt Treatment</th></tr></thead><tbody><tr><td>RT2.45_201.127732</td><td>Aculene A [fragment]+</td><td>0.36</td><td>0.14</td><td>-0.32</td><td>-0.17</td><td>0.08</td><td>0.18</td><td>-0.18</td><td>0.03</td><td>0.16</td><td>0.32</td><td>-0.24</td><td>-0.11</td><td>-0.25</td><td>0.01</td><td>-0.16</td><td>-0.27</td><td>-0.22</td><td>0.13</td></tr><tr><td>RT2.45_316.190751</td><td>Aculene A [M+H]+</td><td>0.36</td><td>0.14</td><td>-0.38</td><td>-0.19</td><td>0.06</td><td>0.36</td><td>-0.21</td><td>0.02</td><td>0.31</td><td>0.30</td><td>-0.33</td><td>-0.17</td><td>-0.31</td><td>0.02</td><td>-0.20</td><td>-0.33</td><td>-0.28</td><td>0.18</td></tr><tr><td>RT2.55_318.206391</td><td>Aculene B [M+H]+</td><td>0.10</td><td>0.10</td><td>-0.57</td><td>-0.28</td><td>-0.09</td><td>0.01</td><td>-0.42</td><td>-0.09</td><td>0.07</td><td>0.20</td><td>-0.38</td><td>-0.26</td><td>-0.20</td><td>0.04</td><td>-0.14</td><td>-0.20</td><td>-0.14</td><td>0.13</td></tr><tr><td>RT2.55_477.150017</td><td>Aspergillus [M+H]+</td><td>0.09</td><td>0.15</td><td>0.64</td><td>0.60</td><td>0.11</td><td>0.13</td><td>0.25</td><td>0.25</td><td>0.11</td><td>0.25</td><td>0.25</td><td>0.34</td><td>0.31</td><td>0.10</td><td>0.25</td><td>0.25</td><td>0.03</td><td>0.27</td></tr><tr><td>RT2.55_499.132347</td><td>Aspergillus [M+Na]+</td><td>0.11</td><td>0.13</td><td>0.60</td><td>0.75</td><td>0.25</td><td>0.25</td><td>0.24</td><td>0.25</td><td>0.11</td><td>0.25</td><td>0.25</td><td>0.33</td><td>0.31</td><td>0.11</td><td>0.25</td><td>0.25</td><td>0.04</td><td>0.28</td></tr><tr><td>RT2.68_291.123173</td><td>Unknown A</td><td>-0.48</td><td>0.30</td><td>-3.18</td><td>-1.48</td><td>-0.77</td><td>0.09</td><td>0.52</td><td>1.70</td><td>-0.25</td><td>0.18</td><td>-1.40</td><td>-0.76</td><td>-0.08</td><td>0.06</td><td>-0.14</td><td>0.06</td><td>-0.02</td><td>0.04</td></tr><tr><td>RT2.71_201.127768</td><td>Aculene C [M+H-H2O]+</td><td>-0.24</td><td>-0.27</td><td>0.05</td><td>0.25</td><td>-0.11</td><td>-0.08</td><td>-0.10</td><td>0.03</td><td>-0.18</td><td>-0.67</td><td>0.03</td><td>0.16</td><td>0.21</td><td>0.01</td><td>-0.08</td><td>0.10</td><td>0.10</td><td>0.10</td></tr><tr><td>RT2.75_241.119986</td><td>Aculene C [M+Na]+</td><td>-0.27</td><td>-0.32</td><td>0.10</td><td>0.28</td><td>-0.10</td><td>-0.08</td><td>-0.09</td><td>0.03</td><td>-0.19</td><td>-0.75</td><td>0.04</td><td>0.16</td><td>0.21</td><td>0.25</td><td>-0.04</td><td>0.10</td><td>0.36</td><td>0.07</td></tr><tr><td>RT2.75_345.131426</td><td>Unknown C</td><td>-1.23</td><td>0.11</td><td>-2.28</td><td>-1.56</td><td>-0.93</td><td>-0.01</td><td>0.34</td><td>1.59</td><td>0.18</td><td>-1.06</td><td>-1.05</td><td>-0.86</td><td>-0.05</td><td>0.05</td><td>-0.14</td><td>0.04</td><td>-0.01</td><td>0.04</td></tr><tr><td>RT2.75_459.2515</td><td>Aculene C [2M+Na]+</td><td>-0.40</td><td>-0.42</td><td>0.03</td><td>0.34</td><td>-0.18</td><td>-0.14</td><td>-0.23</td><td>-0.06</td><td>-0.23</td><td>-1.06</td><td>-0.12</td><td>0.13</td><td>0.08</td><td>0.15</td><td>-0.01</td><td>0.02</td><td>0.17</td><td>0.10</td></tr><tr><td>RT2.77_219.138014</td><td>Aculene C [M+H]+</td><td>-0.33</td><td>-0.34</td><td>0.02</td><td>0.28</td><td>-0.18</td><td>-0.14</td><td>-0.21</td><td>-0.05</td><td>-0.21</td><td>-0.88</td><td>-0.05</td><td>0.13</td><td>0.08</td><td>0.14</td><td>-0.02</td><td>0.02</td><td>0.17</td><td>0.08</td></tr><tr><td>RT3.01_203.142363</td><td>Aculene D [M+H-H2O]+</td><td>-0.34</td><td>-0.39</td><td>-0.16</td><td>0.14</td><td>-0.31</td><td>-0.25</td><td>-0.39</td><td>-0.19</td><td>-0.31</td><td>-0.88</td><td>-0.03</td><td>0.06</td><td>0.07</td><td>0.16</td><td>-0.04</td><td>-0.03</td><td>0.19</td><td>0.05</td></tr><tr><td>RT3.03_201.127843</td><td>Aculene analog [frag]+</td><td>-0.21</td><td>-0.04</td><td>-0.08</td><td>0.05</td><td>-0.12</td><td>-0.09</td><td>-0.13</td><td>0.12</td><td>-0.12</td><td>-0.46</td><td>0.02</td><td>0.10</td><td>0.03</td><td>0.17</td><td>-0.12</td><td>0.07</td><td>0.19</td><td>0.08</td></tr><tr><td>RT3.17_470.217443</td><td>Unknown B [M+H]+</td><td>-0.41</td><td>0.16</td><td>-1.22</td><td>-1.19</td><td>-0.55</td><td>-0.44</td><td>-0.47</td><td>0.42</td><td>0.22</td><td>0.03</td><td>-0.89</td><td>-0.78</td><td>-0.19</td><td>0.02</td><td>-0.23</td><td>0.15</td><td>-0.12</td><td>0.25</td></tr><tr><td>RT3.18_492.200173</td><td>Unknown B [M+Na]+</td><td>-0.44</td><td>0.13</td><td>-1.11</td><td>-1.06</td><td>-0.43</td><td>-0.37</td><td>-0.37</td><td>0.39</td><td>0.25</td><td>0.02</td><td>-0.57</td><td>-0.55</td><td>-0.12</td><td>0.02</td><td>-0.14</td><td>0.09</td><td>-0.07</td><td>0.16</td></tr><tr><td>RT3.79_661.153413</td><td>Secalonic acid [D7] [M+Na]+</td><td>-0.03</td><td>-0.16</td><td>0.09</td><td>0.13</td><td>0.00</td><td>-0.03</td><td>0.04</td><td>0.04</td><td>0.06</td><td>0.03</td><td>0.03</td><td>0.00</td><td>0.15</td><td>0.07</td><td>0.01</td><td>0.14</td><td>0.01</td><td>0.25</td></tr><tr><td>RT3.8_639.171562</td><td>Secalonic acid [D7] [M+H]+</td><td>-0.05</td><td>-0.15</td><td>0.00</td><td>0.05</td><td>-0.02</td><td>-0.06</td><td>-0.02</td><td>-0.01</td><td>0.05</td><td>0.05</td><td>0.03</td><td>0.04</td><td>0.16</td><td>0.18</td><td>-0.03</td><td>0.10</td><td>-0.05</td><td>0.08</td></tr><tr><td>RT3.92_600.264904</td><td>CJ-15183 [M+NH4]+</td><td>0.47</td><td>0.03</td><td>-0.33</td><td>-0.68</td><td>0.52</td><td>0.42</td><td>0.03</td><td>-0.33</td><td>1.00</td><td>1.45</td><td>0.22</td><td>-0.27</td><td>-0.07</td><td>-0.09</td><td>0.03</td><td>-0.10</td><td>-0.14</td><td>-0.08</td></tr><tr><td>RT3.92_605.219954</td><td>CJ-15183 [M+Na]+</td><td>0.44</td><td>0.03</td><td>-0.32</td><td>-0.66</td><td>0.52</td><td>0.42</td><td>0.04</td><td>-0.30</td><td>0.94</td><td>1.35</td><td>0.17</td><td>-0.26</td><td>-0.07</td><td>-0.10</td><td>0.03</td><td>-0.11</td><td>-0.15</td><td>-0.09</td></tr><tr><td>RT4.15_602.281542</td><td>Dihydro-CJ-15183 [M+NH4]+</td><td>0.26</td><td>-0.21</td><td>-0.62</td><td>-0.96</td><td>0.26</td><td>0.16</td><td>-0.25</td><td>-0.67</td><td>0.84</td><td>1.05</td><td>0.07</td><td>-0.54</td><td>-0.12</td><td>-0.13</td><td>0.04</td><td>-0.19</td><td>-0.10</td><td>-0.10</td></tr><tr><td>RT4.15_607.236763</td><td>Dihydro-CJ-15183 [M+Na]+</td><td>0.26</td><td>-0.20</td><td>-0.60</td><td>-0.98</td><td>0.26</td><td>0.16</td><td>-0.24</td><td>-0.65</td><td>0.84</td><td>1.05</td><td>0.09</td><td>-0.51</td><td>-0.12</td><td>-0.14</td><td>0.05</td><td>-0.20</td><td>-0.20</td><td>-0.11</td></tr><tr><td>RT4.28_767.393453</td><td>Unknown D [M+H]+</td><td>-0.22</td><td>0.13</td><td>-0.84</td><td>-0.68</td><td>-0.41</td><td>-0.27</td><td>-0.43</td><td>0.24</td><td>0.04</td><td>-0.34</td><td>-0.81</td><td>-0.45</td><td>-0.08</td><td>-0.01</td><td>-0.09</td><td>0.04</td><td>-0.03</td><td>0.16</td></tr><tr><td>RT4.48_349.154545</td><td>Acu-dioxomorpholine [frag]+</td><td>-0.01</td><td>-0.08</td><td>-0.23</td><td>-0.23</td><td>0.02</td><td>-0.03</td><td>0.03</td><td>-0.71</td><td>-0.11</td><td>-0.19</td><td>-0.04</td><td>-0.40</td><td>0.06</td><td>-0.20</td><td>0.17</td><td>-0.07</td><td>0.01</td><td>-0.06</td></tr><tr><td>RT4.49_361.154289</td><td>Acu-dioxomorpholine [frag]+</td><td>-0.01</td><td>-0.08</td><td>-0.22</td><td>-0.22</td><td>0.02</td><td>-0.03</td><td>0.03</td><td>-0.73</td><td>-0.10</td><td>-0.19</td><td>-0.03</td><td>-0.42</td><td>0.08</td><td>-0.25</td><td>0.22</td><td>-0.08</td><td>0.01</td><td>-0.08</td></tr><tr><td>RT4.5_417.216965</td><td>Acu-dioxomorpholine [M+H]+</td><td>0.01</td><td>-0.08</td><td>-0.23</td><td>-0.02</td><td>0.02</td><td>-0.03</td><td>0.03</td><td>-0.72</td><td>-0.10</td><td>-0.17</td><td>-0.03</td><td>-0.41</td><td>0.12</td><td>-0.43</td><td>-0.15</td><td>-0.15</td><td>0.01</td><td>-0.13</td></tr></tbody></table>

Aculene A (1) and B (2) (Scheme 1) both increased in abundance in response to increased glycerol concentrations and decreased when exposed to an increase in salinity (Figures S3, S5 and S6). Pseudomolecular ions associated with these metabolites had high loading values associated with the separation of osmotic stress treatment samples away from saline treatment samples in PCA models for all three ex-type strains (Figures 6 A–C and 7). The opposite however was observed for aculene C (3) (a presumed biosynthetic precursor to 1 and 2) and a second metabolite (sharing a characteristic fragment ion *m/z* 201.1278 for the aculenes; Figures S3, S5 and S6), which exhibited significant increases in 100% seawater salt treatments for ex-type strains ATCC 16872 and NRRL 20623 and had high loadings values associated with the separation of 100% seawater treatment samples in the PCA models (Figures 6A,B and 7). Decrease in aculene A (1) and B (2) production in saline conditions may be related to the increase in aculene C (3) (which was more pronounced in the ATCC 16872 and NRRL 20623 ex-type strains).
<figure><img src="image_8.png" /><figcaption>**Scheme 1.** Structures of aculenes A–D.</figcaption></figure>

CJ-15,183 (5) (Scheme 2) and an unknown dihydro-CJ-15,183 analog were also observed to increase in conditions of increased osmotic potential. Production differed slightly between ex-type strains in response to osmotic stress. For ex-type strain ATCC 16872, a significant increase in production was observed in 50% glycerol in comparison with observed levels of production at 100% glycerol and in the 0% control (Figure 7 and Figures S3–S12). Both ex-type strains BCRC 32190 and NRRL 20623 demonstrated a significant increase in production under osmotic stress compared to controls for both 50% and 100% glycerol treatments. However for BCRC 32190, production was found to decrease with 50% glycerol and 100% glycerol treatments (similar to ATCC 16872), while production by NRRL 20623 was found to increase significantly with 50% glycerol and increase further with 100% glycerol treatment (Figure 7 and Figures S3–S12).

<figure><img src="image_9.png" /><figcaption>**Scheme 2.** Structure of CJ-15,183.</figcaption></figure>

The metabolite aspergillusol (6) (Scheme 3) was significantly increased under saline conditions from all three ex-type strains. Significant increase in production of 6 was more pronounced in 50% and 100% seawater treatments compared to glycerol treatments for ex-type strain ATCC 16872 (Figure 7, Figures S3 and S10). For this reason, the [M + H]$^+$ and [M + Na]$^+$ pseudomolecular ions for aspergillusol had the highest loading values associated with the separation of saline treatments away from the osmotic stress and control treatments in the ATCC 16872 PCA data model (Figures 6A and 7). Significant increases in 6 production compared to controls for ex-type strain BCRC 32190 were observed for both osmotic stress treatments (50% glycerol and 100% glycerol), albeit greater production was observed in saline treatments (especially in 50% seawater; Figure 7 and Figure S3). This observation was reflected in the PCA analysis of the BCRC 32190 ex-type strain dataset as both pseudomolecular ions observed for aspergillusol had substantial loadings in the PCA model associated with the separation of saline treatments (specifically 50% seawater; Figures 6B and 7). Aspergillusol production by ex-type strain NRRL 20623 was similar to that of BCRC 32190, with significant increases in production associated with 100% glycerol, 50% and 100% seawater treatments (Figure 7, Figures S3 and S15) and PCA loading values associated with the differentiation of 50% and 100% seawater treatments along PC 2 in the NRRL 20623 PCA data model (Figures 6C and 7).
<figure><img src="image_9.png" /><figcaption>**Scheme 1.** Structures of aculenes A–D.</figcaption></figure>

CJ-15,183 (5) (Scheme 2) and an unknown dihydro-CJ-15,183 analog were also observed to increase in conditions of increased osmotic potential. Production differed slightly between ex-type strains in response to osmotic stress. For ex-type strain ATCC 16872, a significant increase in production was observed in 50% glycerol in comparison with observed levels of production at 100% glycerol and in the 0% control (Figure 7 and Figures S3–S12). Both ex-type strains BCRC 32190 and NRRL 20623 demonstrated a significant increase in production under osmotic stress compared to controls for both 50% and 100% glycerol treatments. However for BCRC 32190, production was found to decrease with 50% glycerol and 100% glycerol treatments (similar to ATCC 16872), while production by NRRL 20623 was found to increase significantly with 50% glycerol and increase further with 100% glycerol treatment (Figure 7 and Figures S3–S12).

<figure><img src="image_10.png" /><figcaption>**Scheme 2.** Structure of CJ-15,183.</figcaption></figure>

The metabolite aspergillusol (6) (Scheme 3) was significantly increased under saline conditions from all three ex-type strains. Significant increase in production of 6 was more pronounced in 50% and 100% seawater treatments compared to glycerol treatments for ex-type strain ATCC 16872 (Figure 7, Figures S3 and S10). For this reason, the [M + H]$^+$ and [M + Na]$^+$ pseudomolecular ions for aspergillusol had the highest loading values associated with the separation of saline treatments away from the osmotic stress and control treatments in the ATCC 16872 PCA data model (Figures 6A and 7). Significant increases in 6 production compared to controls for ex-type strain BCRC 32190 were observed for both osmotic stress treatments (50% glycerol and 100% glycerol), albeit greater production was observed in saline treatments (especially in 50% seawater; Figure 7 and Figure S3). This observation was reflected in the PCA analysis of the BCRC 32190 ex-type strain dataset as both pseudomolecular ions observed for aspergillusol had substantial loadings in the PCA model associated with the separation of saline treatments (specifically 50% seawater; Figures 6B and 7). Aspergillusol production by ex-type strain NRRL 20623 was similar to that of BCRC 32190, with significant increases in production associated with 100% glycerol, 50% and 100% seawater treatments (Figure 7, Figures S3 and S15) and PCA loading values associated with the differentiation of 50% and 100% seawater treatments along PC 2 in the NRRL 20623 PCA data model (Figures 6C and 7).
<figure><img src="image_10.png" /><figcaption>**Scheme 3.** Structure of aspergillusol.</figcaption></figure>

Acu-dioxomorpholine (7) (Scheme 4) was found to increase significantly only in the 50% seawater treatment with ex-type strain ATCC 16872 (Figure 7 and Figure S3). However significant decreases in production were observed in the 100% seawater treatment for all three ex-type strains. Pseudomolecular ions associated with acu-dioxomorpholine had an important role in the separation of the 50% seawater treatment samples away from the 100% seawater treatment samples along PC2 in the ATCC 16872 ex-type strain PCA data model (due to high negative loading values; Figures 6A and 7). No significant increase or decrease was observed for acu-dioxomorpholine in either of the glycerol treatments compared to the control for each of the three ex-type strains.

<figure><img src="image_11.png" /><figcaption>**Scheme 4.** Structure of acu-dioxomorpholine.</figcaption></figure>

A significant increase in secalonic acid (Scheme 5) production (likely secalonic acid D as this has been previously reported as being the more abundant secalonic acid produced by *A. aculeatus* [23]) was observed for ex-type strain ATCC 16872 in the 100% seawater treatment, while the 100% glycerol treatment led to a significant decrease in production in comparison to 0% control samples (Figure 7, Figures S3 and S14). The increase in production observed in the 100% seawater treatment had a moderate effect upon the ATCC 16872 ex-type strain PCA data model causing the separation of saline conditions from glycerol conditions along PC1 (positive loading values) and the separation of 100% seawater from 50% seawater treatment samples along PC2 (Figures 6A and 7). A more pronounced increase in both [M + H]$^+$ and [M + Na]$^+$ pseudomolecular ions were observed in 50% seawater treatments for ex-type strain NRRL 20623 (Figures 4 and 7). These variables in particular strongly influenced the separation of the 50% seawater treatment samples from the 100% seawater in the NRRL 20623 ex-type PCA data model due to high negative loading values along PC2 (Figures 6C and 7). A slight increase in production of secalonic acid was observed under saline treatment conditions for ex-type strain BCRC 32190; however these increases were not statistically significant.
<figure><img src="image_11.png" /><figcaption>**Scheme 3.** Structure of aspergillusol.</figcaption></figure>

Acu-dioxomorpholine (7) (Scheme 4) was found to increase significantly only in the 50% seawater treatment with ex-type strain ATCC 16872 (Figure 7 and Figure S3). However significant decreases in production were observed in the 100% seawater treatment for all three ex-type strains. Pseudomolecular ions associated with acu-dioxomorpholine had an important role in the separation of the 50% seawater treatment samples away from the 100% seawater treatment samples along PC2 in the ATCC 16872 ex-type strain PCA data model (due to high negative loading values; Figures 6A and 7). No significant increase or decrease was observed for acu-dioxomorpholine in either of the glycerol treatments compared to the control for each of the three ex-type strains.

<figure><img src="image_12.png" /><figcaption>**Scheme 4.** Structure of acu-dioxomorpholine.</figcaption></figure>

A significant increase in secalonic acid (Scheme 5) production (likely secalonic acid D as this has been previously reported as being the more abundant secalonic acid produced by *A. aculeatus* [23]) was observed for ex-type strain ATCC 16872 in the 100% seawater treatment, while the 100% glycerol treatment led to a significant decrease in production in comparison to 0% control samples (Figure 7, Figures S3 and S14). The increase in production observed in the 100% seawater treatment had a moderate effect upon the ATCC 16872 ex-type strain PCA data model causing the separation of saline conditions from glycerol conditions along PC1 (positive loading values) and the separation of 100% seawater from 50% seawater treatment samples along PC2 (Figures 6A and 7). A more pronounced increase in both [M + H]$^+$ and [M + Na]$^+$ pseudomolecular ions were observed in 50% seawater treatments for ex-type strain NRRL 20623 (Figures 4 and 7). These variables in particular strongly influenced the separation of the 50% seawater treatment samples from the 100% seawater in the NRRL 20623 ex-type PCA data model due to high negative loading values along PC2 (Figures 6C and 7). A slight increase in production of secalonic acid was observed under saline treatment conditions for ex-type strain BCRC 32190; however these increases were not statistically significant.
<figure><img src="image_12.png" /><figcaption>**Scheme 5.** Structures of secalonic acids B, D and F.</figcaption></figure>

**4. Discussion**

When a fungal species is formally described, often a living strain is deposited in a recognized public culture collection in order to fix the name to a living strain; this strain is designated as the type strain and ideally, multiple copies are stored in a metabolically inactive state, either in a lyophilized or ultra-cold frozen state. Descendent of the original living tissues obtained from the type strain are designated as ex-type strains and are often distributed to multiple public service collections. Traditionally, genetic resource centers (such as culture collections) support scientific research by providing cultures, including ex-type cultures, to scientists [28]. Proper culture collection curation and handling of strains is essential to preserve ex-type strains that remain representative of the original type specimen. Phenotypic deviations are frequently reported by researchers when investigating ex-type and conspecific strains obtained from various different culture collections [29]. In our case, *A. aculeatus* is a black-spored species (*Aspergillus* section *Nigri*) that is morphologically similar to at least 18 different black-spored species [30]. PCR-fingerprinting is a reliable means of confirming the correct designation of ex-type strains obtained from different culture collections [29]. In our experiment, comparison of DNA sequences (ITS1-5.8s-ITS2 and the D1/D2 variable region of the nLSU rDNA) confirmed that the *A. aculeatus* ex-type strains ATCC 16872, BCRC 32190 and NRRL 20623 were the same species (100% similarity, data not shown). Therefore, differences in the observed secondary metabolite profiles were assumed to be ex-type strain dependent, likely due to somatic changes and degeneration or mutations occurring during successive transfers between periods of long-term storage of the ex-type stock strains.

**4.1. Site-to-Site Variability**

Although a standardized protocol was circulated among the participating laboratories to minimize operational variations, site-specific variations in secondary metabolism were observed (Figure 2; comparison of UPLC-HRMS profiles in Figure S2). During initial data analysis, various pseudomolecular ions representing plasticizers and detergents (“MS contaminants”), were found in all extracts, and the constituents varied between sites. Such contaminants originate from detergents and plasticware used in the laboratory (e.g., sample bottles, vials, pipette tips, pipette bulbs and filter membranes) that differ from lab-to-lab [27]. After removal of known contaminant ions from the data set, site-to-site variation was still observed during data analysis.

Variations in secondary metabolite production likely resulted from both biotic and abiotic factors. Variations in fermentation conditions and nutritional composition affect secondary metabolite expression [31,32]. Even for labs that exchanged the same ex-type strain (NRRL 20623), secondary
<figure><img src="image_13.png" /><figcaption>**Scheme 5.** Structures of secalonic acids B, D and F.</figcaption></figure>

**4. Discussion**

When a fungal species is formally described, often a living strain is deposited in a recognized public culture collection in order to fix the name to a living strain; this strain is designated as the type strain and ideally, multiple copies are stored in a metabolically inactive state, either in a lyophilized or ultra-cold frozen state. Descendent of the original living tissues obtained from the type strain are designated as ex-type strains and are often distributed to multiple public service collections. Traditionally, genetic resource centers (such as culture collections) support scientific research by providing cultures, including ex-type cultures, to scientists [28]. Proper culture collection curation and handling of strains is essential to preserve ex-type strains that remain representative of the original type specimen. Phenotypic deviations are frequently reported by researchers when investigating ex-type and conspecific strains obtained from various different culture collections [29]. In our case, *A. aculeatus* is a black-spored species (*Aspergillus* section *Nigri*) that is morphologically similar to at least 18 different black-spored species [30]. PCR-fingerprinting is a reliable means of confirming the correct designation of ex-type strains obtained from different culture collections [29]. In our experiment, comparison of DNA sequences (ITS1-5.8s-ITS2 and the D1/D2 variable region of the nLSU rDNA) confirmed that the *A. aculeatus* ex-type strains ATCC 16872, BCRC 32190 and NRRL 20623 were the same species (100% similarity, data not shown). Therefore, differences in the observed secondary metabolite profiles were assumed to be ex-type strain dependent, likely due to somatic changes and degeneration or mutations occurring during successive transfers between periods of long-term storage of the ex-type stock strains.

**4.1. Site-to-Site Variability**

Although a standardized protocol was circulated among the participating laboratories to minimize operational variations, site-specific variations in secondary metabolism were observed (Figure 2; comparison of UPLC-HRMS profiles in Figure S2). During initial data analysis, various pseudomolecular ions representing plasticizers and detergents (“MS contaminants”), were found in all extracts, and the constituents varied between sites. Such contaminants originate from detergents and plasticware used in the laboratory (e.g., sample bottles, vials, pipette tips, pipette bulbs and filter membranes) that differ from lab-to-lab [27]. After removal of known contaminant ions from the data set, site-to-site variation was still observed during data analysis.

Variations in secondary metabolite production likely resulted from both biotic and abiotic factors. Variations in fermentation conditions and nutritional composition affect secondary metabolite expression [31,32]. Even for labs that exchanged the same ex-type strain (NRRL 20623), secondary
metabolite production differed sufficiently between sites to form distinct site-specific clades as recognized by hierarchical cluster analysis (Figure 2). Media composition may have played a minimal role because most of the ingredients were defined reagent grade chemicals, and all used the same brand of yeast autolysate (Sigma Aldrich 73145). However, the source of seawater varied locally (and likely varied in dissolved organic and inorganic matter content and possibly pH). The recognition of site-specific clades during hierarchical clustering of the control conditions (where seawater and glycerol were absent; Figure 3) for sites that fermented the same ex-type strain indicated that additional factors influenced on secondary metabolism. Sources of undefined medium ingredients (e.g., yeast and malt extracts), may differ in composition even between lot numbers from the same manufacturer, and these subtle difference have a direct impact on secondary metabolite expression from the same strain [33]. Abiotic factors such as aeration may have also been a factor as sites differed in how the Erlenmeyer flasks were closed (some pressure-tight caps, cellulose stoppers, or cotton). While it is impossible to identify specific factors that caused the differences in secondary metabolite production, these data highlight the subtleties in culture conditions that may influence natural product biosynthesis.

## 4.2. *Differences in Ex-Type Strain Secondary Metabolism*

When each of the ex-type strains originating from different culture collections were fermented at a single site, only the fermentation condition (saline and osmotic stress) and ex-type strain were found to contribute to the variation in secondary metabolism. Ex-type strains did not segregate uniformly based on treatment following PCA and cluster analysis (Figure 4). Theoretically, genetically equivalent ex-type strains should express a consistent pattern of secondary metabolites under a given fermentation condition as all originated from the same parent. However, this was not the case in this experiment. For example, metabolites 5'-hydroxyasperentin 8-O-methyl ether and alantyleunone increased under osmotic stress for ex-type strains ATCC 16872 and NRRL 20623, while production increased under saline stress for ex-type strain BCRC 32190 (Figure S3).

During prolonged culturing and transfers, filamentous fungi undergo spontaneous somatic changes, referred to as sectoring. Sectors are regions that commonly form in agar culture that can differ in morphology (e.g., lack sporulation), secondary metabolism, and enzyme production from the parent inoculum. Characteristics of a culture sector may remain stable after subculturing. However, sector formation is often associated with physiological instability and frequently results in changes in secondary metabolite production where gains as well as losses have been reported [34–37]. Although the production of some secondary metabolites can be correlated with morphological changes, e.g., loss of sporulation [38], degeneration in sporulation is not necessarily associated with a loss in secondary metabolite production [33,36]. Similarly, heavily mutated factory production strains may have aberrant morphology, including loss of hydrophobicity, pigmentation, aerial hyphae, and sporulation [39,40].

The *A. aculeatus* ex-type strains used in this study originated from the type strain WB 5094 [41]. As the type strain material has been in the collection at least sixty years, repeated handling of ex-type progeny over time at different locations and with different transfer histories has resulted in three ex-type strains that differ in secondary metabolite expression. The stability of secondary metabolite profiles deteriorate over longer storage periods, which can be minimized (but not entirely prevented) by storage in a metabolically inactive state such as lyophilization or cryostorage [28,35].

Annotation of the genome of *A. aculeatus* ex-type strain ATCC 16872 indicates >60 gene clusters encoding non-ribosomal peptides, polyketides, terpenoids, and other secondary metabolites, not all of which have been linked with known molecular structure [18]. However, some of the secondary metabolites that have previously been reported by this strain and for which the gene cluster was present, were not observed in this experiment, e.g., aculeacin [42]. For experimental simplicity (and to minimize variability between labs), the extraction protocol used was not exhaustive and would have favored the extraction of metabolites that were present in the culture broth. Secondary metabolites that accumulate within the mycelium or spores of the fungus would therefore have been less prevalent in the extract and might have been missed in the analysis due to their low abundance in the overall extract
metabolite production differed sufficiently between sites to form distinct site-specific clades as recognized by hierarchical cluster analysis (Figure 2). Media composition may have played a minimal role because most of the ingredients were defined reagent grade chemicals, and all used the same brand of yeast autolysate (Sigma Aldrich 73145). However, the source of seawater varied locally (and likely varied in dissolved organic and inorganic matter content and possibly pH). The recognition of site-specific clades during hierarchical clustering of the control conditions (where seawater and glycerol were absent; Figure 3) for sites that fermented the same ex-type strain indicated that additional factors influenced on secondary metabolism. Sources of undefined medium ingredients (e.g., yeast and malt extracts), may differ in composition even between lot numbers from the same manufacturer, and these subtle difference have a direct impact on secondary metabolite expression from the same strain [33]. Abiotic factors such as aeration may have also been a factor as sites differed in how the Erlenmeyer flasks were closed (some pressure-tight caps, cellulose stoppers, or cotton). While it is impossible to identify specific factors that caused the differences in secondary metabolite production, these data highlight the subtleties in culture conditions that may influence natural product biosynthesis.

## 4.2. *Differences in Ex-Type Strain Secondary Metabolism*

When each of the ex-type strains originating from different culture collections were fermented at a single site, only the fermentation condition (saline and osmotic stress) and ex-type strain were found to contribute to the variation in secondary metabolism. Ex-type strains did not segregate uniformly based on treatment following PCA and cluster analysis (Figure 4). Theoretically, genetically equivalent ex-type strains should express a consistent pattern of secondary metabolites under a given fermentation condition as all originated from the same parent. However, this was not the case in this experiment. For example, metabolites 5'-hydroxyasperentin 8-O-methyl ether and alantyleunone increased under osmotic stress for ex-type strains ATCC 16872 and NRRL 20623, while production increased under saline stress for ex-type strain BCRC 32190 (Figure S3).

During prolonged culturing and transfers, filamentous fungi undergo spontaneous somatic changes, referred to as sectoring. Sectors are regions that commonly form in agar culture that can differ in morphology (e.g., lack sporulation), secondary metabolism, and enzyme production from the parent inoculum. Characteristics of a culture sector may remain stable after subculturing. However, sector formation is often associated with physiological instability and frequently results in changes in secondary metabolite production where gains as well as losses have been reported [34–37]. Although the production of some secondary metabolites can be correlated with morphological changes, e.g., loss of sporulation [38], degeneration in sporulation is not necessarily associated with a loss in secondary metabolite production [33,36]. Similarly, heavily mutated factory production strains may have aberrant morphology, including loss of hydrophobicity, pigmentation, aerial hyphae, and sporulation [39,40].

The *A. aculeatus* ex-type strains used in this study originated from the type strain WB 5094 [41]. As the type strain material has been in the collection at least sixty years, repeated handling of ex-type progeny over time at different locations and with different transfer histories has resulted in three ex-type strains that differ in secondary metabolite expression. The stability of secondary metabolite profiles deteriorate over longer storage periods, which can be minimized (but not entirely prevented) by storage in a metabolically inactive state such as lyophilization or cryostorage [28,35].

Annotation of the genome of *A. aculeatus* ex-type strain ATCC 16872 indicates >60 gene clusters encoding non-ribosomal peptides, polyketides, terpenoids, and other secondary metabolites, not all of which have been linked with known molecular structure [18]. However, some of the secondary metabolites that have previously been reported by this strain and for which the gene cluster was present, were not observed in this experiment, e.g., aculeacin [42]. For experimental simplicity (and to minimize variability between labs), the extraction protocol used was not exhaustive and would have favored the extraction of metabolites that were present in the culture broth. Secondary metabolites that accumulate within the mycelium or spores of the fungus would therefore have been less prevalent in the extract and might have been missed in the analysis due to their low abundance in the overall extract
composition. Pseudomolecular ions relating to several metabolites such as aculeacin, calibistrins, auerospurones were not observed in culture extracts when resulting chromatograms were examined using single ion monitoring. Pseudomolecular ions for other metabolites such as the okaramines and asperaculane B were observed using single ion monitoring but occurred at ion intensities below the noise cutoff threshold used for data analysis. For other metabolites such as the acucalibistrins and asperaculane A where pseudomolecular ions were detected at an intensity above the noise threshold, metabolite changes between fermentation conditions were insignificant. Only half of the variables in the data model could be linked with secondary metabolites previously reported in the *A. aculeatus* literature; the remaining variables represent a considerable number of secondary metabolites present in the metabolome of *A. aculeatus* that have yet to be structurally confirmed.

### 4.3. Effect of Increased Salinity/Osmotic Stress on Secondary Metabolism

Media compositions containing either increasing concentrations of glycerol or sea salts create an environment of low water potential and thus provide a stress condition for fungi. Some fungi have evolved cellular mechanisms to cope with environments with low water potential, and for this reason, are considered as osmotolerant. In the marine environment, increased salinity exerts an additional stress due to the presence of relatively high concentrations of sodium and chloride ions, which potentially exert toxic effects on cellular processes [43]. Genome characteristics of halophilic fungi include enrichment in genes related to stress response proteins, ion/metabolite transporters, hydrophobins, and polyol biosynthesis (including the HOG pathway) [6–8]. Although several mechanisms such as increased polyol biosynthesis are a common response shared between osmotolerant and halotolerant fungi to low water potential, adaptations involving maintenance and efflux of metal ions are only expressed under saline stress induced by exposure to sea salts. Transcriptome studies have demonstrated differences in expression of both primary and secondary metabolites due to salt stress; however, differential induction of several secondary metabolites were found to occur in seawater fermentations compared to fermentations using NaCl alone, where authors attributed induction to occur as a result of upregulation of metal ion efflux [44]. In our project, upregulation of secondary metabolites was found to differ between osmotic stress (induced by glycerol) and saline stress (induced by sea salts). Specifically, production of the metabolites CJ-15,183, a dihydro-CJ-15,183 derivative and aculenes A and B were found to increase in response to increasing glycerol concentration, while production under saline stress conditions was consistent with that of non-stress conditions (0% seawater/0% glycerol control). Other secondary metabolites such as aspergillusol, secalonic acid D, aculene C and another aculene analog were found to be upregulated in saline conditions. Interestingly, aspergillusol was first described from a marine-derived strain (CRI323-04) of *A. aculeatus* fermented in potato dextrose broth using seawater [22]. At the time, Ingavat and collaborators attributed the discovery of aspergillusol to the marine-derived nature of isolate, as it was believed to be adapted to the marine environment and therefore a unique source of new natural products. Based on our results with a terrestrial strain of *A. aculeatus*, expression of aspergillusol was confirmed to be associated with a halotolerant response to saline conditions involving seawater/sea salts.

Many fungal genera such as *Aspergillus*, *Penicillium*, *Trichoderma* and *Cladosporium* have the genetic capacity for prolific production of secondary metabolites, yet generally do not express their full potential for secondary metabolite biosynthesis when cultured in the laboratory. Historically, secondary metabolites associated with saline stress have been overlooked by natural product discovery programs using terrestrial isolates, as the majority of investigators rarely reported fermentation of terrestrial strains using a salt water medium. However, media with high solute concentrations, especially sugars are frequently used. On the other hand, researchers working on marine-derived strains often deliberately include seawater in the fermentation media. This is evident from the numerous reports of new secondary metabolites discovered from commonly occurring, halotolerant fungal genera that have been designated as "marine-derived" [9]. Our findings support the hypothesis that media amended
composition. Pseudomolecular ions relating to several metabolites such as aculeacin, calibistrins, auerospurones were not observed in culture extracts when resulting chromatograms were examined using single ion monitoring. Pseudomolecular ions for other metabolites such as the okaramines and asperaculane B were observed using single ion monitoring but occurred at ion intensities below the noise cutoff threshold used for data analysis. For other metabolites such as the acucalibistrins and asperaculane A where pseudomolecular ions were detected at an intensity above the noise threshold, metabolite changes between fermentation conditions were insignificant. Only half of the variables in the data model could be linked with secondary metabolites previously reported in the *A. aculeatus* literature; the remaining variables represent a considerable number of secondary metabolites present in the metabolome of *A. aculeatus* that have yet to be structurally confirmed.

### 4.3. Effect of Increased Salinity/Osmotic Stress on Secondary Metabolism

Media compositions containing either increasing concentrations of glycerol or sea salts create an environment of low water potential and thus provide a stress condition for fungi. Some fungi have evolved cellular mechanisms to cope with environments with low water potential, and for this reason, are considered as osmotolerant. In the marine environment, increased salinity exerts an additional stress due to the presence of relatively high concentrations of sodium and chloride ions, which potentially exert toxic effects on cellular processes [43]. Genome characteristics of halophilic fungi include enrichment in genes related to stress response proteins, ion/metabolite transporters, hydrophobins, and polyol biosynthesis (including the HOG pathway) [6–8]. Although several mechanisms such as increased polyol biosynthesis are a common response shared between osmotolerant and halotolerant fungi to low water potential, adaptations involving maintenance and efflux of metal ions are only expressed under saline stress induced by exposure to sea salts. Transcriptome studies have demonstrated differences in expression of both primary and secondary metabolites due to salt stress; however, differential induction of several secondary metabolites were found to occur in seawater fermentations compared to fermentations using NaCl alone, where authors attributed induction to occur as a result of upregulation of metal ion efflux [44]. In our project, upregulation of secondary metabolites was found to differ between osmotic stress (induced by glycerol) and saline stress (induced by sea salts). Specifically, production of the metabolites CJ-15,183, a dihydro-CJ-15,183 derivative and aculenes A and B were found to increase in response to increasing glycerol concentration, while production under saline stress conditions was consistent with that of non-stress conditions (0% seawater/0% glycerol control). Other secondary metabolites such as aspergillusol, secalonic acid D, aculene C and another aculene analog were found to be upregulated in saline conditions. Interestingly, aspergillusol was first described from a marine-derived strain (CRI323-04) of *A. aculeatus* fermented in potato dextrose broth using seawater [22]. At the time, Ingavat and collaborators attributed the discovery of aspergillusol to the marine-derived nature of isolate, as it was believed to be adapted to the marine environment and therefore a unique source of new natural products. Based on our results with a terrestrial strain of *A. aculeatus*, expression of aspergillusol was confirmed to be associated with a halotolerant response to saline conditions involving seawater/sea salts.

Many fungal genera such as *Aspergillus*, *Penicillium*, *Trichoderma* and *Cladosporium* have the genetic capacity for prolific production of secondary metabolites, yet generally do not express their full potential for secondary metabolite biosynthesis when cultured in the laboratory. Historically, secondary metabolites associated with saline stress have been overlooked by natural product discovery programs using terrestrial isolates, as the majority of investigators rarely reported fermentation of terrestrial strains using a salt water medium. However, media with high solute concentrations, especially sugars are frequently used. On the other hand, researchers working on marine-derived strains often deliberately include seawater in the fermentation media. This is evident from the numerous reports of new secondary metabolites discovered from commonly occurring, halotolerant fungal genera that have been designated as "marine-derived" [9]. Our findings support the hypothesis that media amended
with sea salts (or glycerol) applies a selective pressure on the secondary metabolic output of a fungus, regardless of whether the isolate was of marine or terrestrial origin. Moreover, upregulated production of particular secondary metabolites differs depending on the selective pressure applied (increasing glycerol concentrations vs. increasing seawater salinity). Implications of this discovery suggest that deliberate manipulation of halotolerant fungal strains by additions of solutes, e.g., seawater, salts, or other osmolites, to fermentation media will amplify the kinds of secondary metabolites discovered from these species. Importantly, our findings also highlight the value of multi-investigator, multi-site studies to rigorously assess the effects of the manipulation of culture conditions on natural product production in fungi. It also documents the effect of unforeseen variables on the metabolome extracted from a fungal fermentation, e.g., plasticizers present in common lab plastic-ware, sources of seawater, and degree of culture aeration.

**5. Conclusions**

This exercise yielded notable conclusions. Firstly, the *A. aculeatus* ex-type strain obtained from different culture collections produced different secondary metabolite expression profiles when fermented at different sites despite using the same fermentation protocols. Secondly, the secondary metabolite profile of *A. aculeatus* fermentations was demonstrated to be responsive to osmotic stress, whether induced by glycerol or seawater. Thirdly, *A. aculeatus* ex-type strains obtained from different culture collections produced different secondary metabolite expression profiles when fermented at the same site using identical fermentation protocols.

**Supplementary Materials:** The following are available online at www.mdpi.com/1660-3397/15/8/254/s1.

**Acknowledgments:** Russell Kerr acknowledges the assistance of Nadia Prigoda-Lee, Marius Grote, Kate McQuillan and Stephanie Duffy, and generous financial support from NSERC, the Canada Research Chair program, the Jeanne and Jean-Louis Lévesque Foundation and the Atlantic Canada Opportunities Agency. Ka-Lai Pang thanks the president of National Taiwan Ocean University, Ching-Fong Chang, for a special fund to attend the workshop held in Charlottetown, Canada in 2014 where this work was discussed. Rob Capon and Zhuo Shang acknowledge support from the University of Queensland, and the UQ Institute for Molecular Bioscience. Zhuo Shang acknowledges the provision of an International Postgraduate Research Scholarship (IPRS) and a Centennial Scholarship by the University of Queensland. Catherine Roullier acknowledges the assistance of Marie-Claude Boumard and Thibaut Robiou du Pont, and support from Region Pays de la Loire, France.

**Author Contributions:** Russell Kerr, Catherine Roullier, Ka-Lai Pang, Rainer Ebel, Rob Capon, Gerald Bills and David Overy conceived and designed the experiments; David Overy, Hebelin Correa, Wei-Chiung Chi, Mostafa Rateb and Zhuo Shang performed the experiments; David Overy analyzed the data; David Overy wrote the paper and Catherine Roullier, Ka-Lai Pang, Rainer Ebel, Rob Capon, Gerald Bills provided editorial assistance.

**Conflicts of Interest:** The authors declare no conflict of interest.

**References**

1. *Secondary Metabolism and Differentiation in Fungi*; Bennett, J.W.; Ciegler, A., Eds.; Marcel Dekker: New York, NY, USA, 1983.

2. Demain, A.L. Regulation of secondary metabolism in fungi. *Pure Appl. Chem.* **1986**, *58*, 219–226. [CrossRef]

3. Lim, F.; Sanchez, J.; Wang, C.C.; Keller, N. Toward awakening cryptic secondary metabolite gene clusters in filamentous fungi. *Methods Enzymol.* **2012**, *517*, 303–324. [PubMed]

4. Keller, N.P. Translating biosynthetic gene clusters into fungal armor and weaponry. *Nat. Chem. Biol.* **2015**, *11*, 671–677. [CrossRef] [PubMed]

5. Pang, K.-L.; Overy, D.P.; Jones, E.B.G.; Calado, M.; Burgaud, G.; Walker, A.K.; Johnson, J.A.; Kerr, R.G.; Cha, H.-J.; Bills, G.F. 'Marine fungi' and 'marine-derived fungi' in natural product chemistry research: Toward a new consensual definition. *Fungal Biol. Rev.* **2016**, *30*, 163–175. [CrossRef]

6. Lenassi, M.; Gostinčar, C.; Jackman, S.; Turk, M.; Sadowski, I.; Nislow, C.; Jones, S.; Birol, I.; Gunde-Cimerman, N. Whole genome duplication and enrichment of metal cation transporters revealed by de novo genome sequencing of extremely halotolerant black yeast *Hortaea werneckii*. *PLoS ONE* **2013**, *8*, e71328. [CrossRef] [PubMed]
with sea salts (or glycerol) applies a selective pressure on the secondary metabolic output of a fungus, regardless of whether the isolate was of marine or terrestrial origin. Moreover, upregulated production of particular secondary metabolites differs depending on the selective pressure applied (increasing glycerol concentrations vs. increasing seawater salinity). Implications of this discovery suggest that deliberate manipulation of halotolerant fungal strains by additions of solutes, e.g., seawater, salts, or other osmolites, to fermentation media will amplify the kinds of secondary metabolites discovered from these species. Importantly, our findings also highlight the value of multi-investigator, multi-site studies to rigorously assess the effects of the manipulation of culture conditions on natural product production in fungi. It also documents the effect of unforeseen variables on the metabolome extracted from a fungal fermentation, e.g., plasticizers present in common lab plastic-ware, sources of seawater, and degree of culture aeration.

**5. Conclusions**

This exercise yielded notable conclusions. Firstly, the *A. aculeatus* ex-type strain obtained from different culture collections produced different secondary metabolite expression profiles when fermented at different sites despite using the same fermentation protocols. Secondly, the secondary metabolite profile of *A. aculeatus* fermentations was demonstrated to be responsive to osmotic stress, whether induced by glycerol or seawater. Thirdly, *A. aculeatus* ex-type strains obtained from different culture collections produced different secondary metabolite expression profiles when fermented at the same site using identical fermentation protocols.

**Supplementary Materials:** The following are available online at www.mdpi.com/1660-3397/15/8/254/s1.

**Acknowledgments:** Russell Kerr acknowledges the assistance of Nadia Prigoda-Lee, Marius Grote, Kate McQuillan and Stephanie Duffy, and generous financial support from NSERC, the Canada Research Chair program, the Jeanne and Jean-Louis Lévesque Foundation and the Atlantic Canada Opportunities Agency. Ka-Lai Pang thanks the president of National Taiwan Ocean University, Ching-Fong Chang, for a special fund to attend the workshop held in Charlottetown, Canada in 2014 where this work was discussed. Rob Capon and Zhuo Shang acknowledge support from the University of Queensland, and the UQ Institute for Molecular Bioscience. Zhuo Shang acknowledges the provision of an International Postgraduate Research Scholarship (IPRS) and a Centennial Scholarship by the University of Queensland. Catherine Roullier acknowledges the assistance of Marie-Claude Boumard and Thibaut Robiou du Pont, and support from Region Pays de la Loire, France.

**Author Contributions:** Russell Kerr, Catherine Roullier, Ka-Lai Pang, Rainer Ebel, Rob Capon, Gerald Bills and David Overy conceived and designed the experiments; David Overy, Hebelin Correa, Wei-Chiung Chi, Mostafa Rateb and Zhuo Shang performed the experiments; David Overy analyzed the data; David Overy wrote the paper and Catherine Roullier, Ka-Lai Pang, Rainer Ebel, Rob Capon, Gerald Bills provided editorial assistance.

**Conflicts of Interest:** The authors declare no conflict of interest.

**References**

1. *Secondary Metabolism and Differentiation in Fungi*; Bennett, J.W.; Ciegler, A., Eds.; Marcel Dekker: New York, NY, USA, 1983.

2. Demain, A.L. Regulation of secondary metabolism in fungi. *Pure Appl. Chem.* **1986**, *58*, 219–226. [CrossRef]

3. Lim, F.; Sanchez, J.; Wang, C.C.; Keller, N. Toward awakening cryptic secondary metabolite gene clusters in filamentous fungi. *Methods Enzymol.* **2012**, *517*, 303–324. [PubMed]

4. Keller, N.P. Translating biosynthetic gene clusters into fungal armor and weaponry. *Nat. Chem. Biol.* **2015**, *11*, 671–677. [CrossRef] [PubMed]

5. Pang, K.-L.; Overy, D.P.; Jones, E.B.G.; Calado, M.; Burgaud, G.; Walker, A.K.; Johnson, J.A.; Kerr, R.G.; Cha, H.-J.; Bills, G.F. 'Marine fungi' and 'marine-derived fungi' in natural product chemistry research: Toward a new consensual definition. *Fungal Biol. Rev.* **2016**, *30*, 163–175. [CrossRef]

6. Lenassi, M.; Gostinčar, C.; Jackman, S.; Turk, M.; Sadowski, I.; Nislow, C.; Jones, S.; Birol, I.; Gunde-Cimerman, N. Whole genome duplication and enrichment of metal cation transporters revealed by de novo genome sequencing of extremely halotolerant black yeast *Hortaea werneckii*. *PLoS ONE* **2013**, *8*, e71328. [CrossRef] [PubMed]
7. Zajc, J.; Liu, Y.; Dai, W.; Yang, Z.; Hu, J.; Gostincar, C.; Gunde-Cimerman, N. Genome and transcriptome sequencing of the halophilic fungus *Wallolemia ichthyophaga*: Loadadaptations present and absent. *BMC Genom.* **2013**, *14*, 617. [CrossRef] [PubMed]

8. Kis-Papo, T.; Weig, A.R.; Riley, R.; Peršoh, D.; Salamov, A.; Sun, H.; Lipzen, A.; Wasser, S.P.; Rambold, G.; Grigoriev, I.V.; et al. Genomic adaptations of the halophilic Ded Sea filamentous fungus *Eurotium rubrum*. *Nat. Commun.* **2014**, *5*, 3745. [CrossRef] [PubMed]

9. Wang, Y.; Lu, Z.; Sun, K.; Zhu, W. Effects of high salt stress on secondary metabolite production in the marine-derived fungus *Spicaria elegans*. *Marine Drugs* **2011**, *9*, 535–542. [CrossRef] [PubMed]

10. Wang, Y.; Zheng, J.; Liu, P.; Wang, W.; Zhu, W. Three new compounds from *Aspergillus terreus* PT06-2 grown in a high salt medium. *Marine Drugs* **2011**, *9*, 1368–1378. [CrossRef] [PubMed]

11. Bugni, T.; Ireland, C. Marine-derived fungi: A chemically and biologically diverse group of microorganisms. *Nat. Prod. Rep.* **2004**, *21*, 143–163. [CrossRef] [PubMed]

12. Rateb, M.E.; Ebel, R. Secondary metabolites of fungi from marine habitats. *Nat. Prod. Rep.* **2012**, *28*, 290–344. [CrossRef] [PubMed]

13. Overy, D.P.; Bayman, P.; Kerr, R.G.; Bills, G.F. An assessment of natural product discovery from marine (*sensu strictu*) and marine-derived fungi. *Mycology* **2014**, *5*, 145–167. [CrossRef] [PubMed]

14. Bérdy, J. Thoughts and facts about antibiotics: Where we are now and where we are heading. *J. Antibiot.* **2012**, *65*, 385–395. [CrossRef] [PubMed]

15. Capon, R.J.; Ratnayake, R.; Stewart, M.; Lacey, E.; Tennant, S.; Gill, J.H. Aspergillaxines A-E: Novel heterocyclic dipeptides from an Australian strain of *Aspergillus unilateralis*. *Org. Biomol. Chem.* **2005**, *3*, 123–129. [CrossRef] [PubMed]

16. Wang, G.Y.S.; Arbell, L.M.; Avelar, A.; Borgeson, B.M.; Crews, P. New hisutane based sesquiterpenes from salt water cultures of a marine sponge-derived fungus and the terrestrial fungus *Coriolus consors*. *Tetrahedron* **1998**, *54*, 7335–7342. [CrossRef]

17. The draft genome sequence of *A. aculeatus*. Available online: http://genome.jgi.doe.gov/Aspac1/Aspac1.home.html (accessed on 15 August 2016).

18. De Vries, R.P.; Riley, R.; Wiebenga, A.; Aguilar-Osorio, G.; Amillis, S.; Uchima, C.A.; Anderluh, G.; Asadollahi, M.; Askin, M.; Barry, K.; et al. Comparative genomics reveals high biological diversity and specific adaptations in the industrially and medically important fungal genus *Aspergillus*. *Genome Biol.* **2017**, *18*, 28. [CrossRef] [PubMed]

19. Gaude, E.; Chignola, F.; Spiliotopoulus, D.; Spitaleri, A.; Ghitti, M.; Garcia-Manteiga, J.M.; Mari, S.; Musco, G. *muma*, an R package for metabolomics univariate and multivariate statistical analysis. *Curr. Metab.* **2013**, *1*, 180–189. [CrossRef]

20. Anderson, R.; Buchi, G.; Kobbe, B.; Demain, A.L. Secalonic acids D and F are toxic metabolites of *Aspergillus aculeatus*. *J. Org. Chem.* **1977**, *42*, 352–353. [CrossRef]

21. Gao, Y.-Q.; Guo, C.-J.; Zhang, Q.; Zhou, W.-M.; Wang, C.C.C.; Gao, J.-M. Asperaculanes A and, B.; two sesquiterpenoids from the fungus *Aspergillus aculeatus*. *Molecules* **2015**, *20*, 325–334. [CrossRef] [PubMed]

22. Ingavat, N.; Dobereiner, J.; Wiyakrutta, S.; Mahidol, C.; Ruchirawat, S.; Kittakoop, P. Aspergillusol A., a α-glucosidase inhibitor from the marine-derived fungus *Aspergillus aculeatus*. *J. Nat. Prod.* **2009**, *72*, 2049–2052. [CrossRef] [PubMed]

23. Kurobane, I.; Vining, L.C. Biosynthetic realtionships among the secalonic acids. Isolation of emodin, endocrocin and secalonic acids from *Pyrenochaeta terrestris* and *Aspergillus aculeatus*. *J. Antibiot.* **1979**, *32*, 1256–1266. [CrossRef] [PubMed]

24. Petersen, L.M.; Hoeck, C.; Frisvad, J.C.; Gotfredsen, C.H.; Larsen, T.O. Dereplication guided discovery of secondary metabolites of mixed biosynthetic origin from *Aspergillus aculeatus*. *Molecules* **2014**, *19*, 10898–10921. [CrossRef] [PubMed]

25. Petersen, L.M.; Holm, D.K.; Gotfredsen, C.H.; Mortensen, U.H.; Larsen, T.O. Investigation of a 6-MSA synthase gene cluster in *Aspergillus aculeatus* reveals 6-MSA-derived aculinic acid, aculins A-B and epi-aculin, *A. ChemBioChem* **2015**, *16*, 2200–2204. [CrossRef] [PubMed]

26. Watanabe, S.; Hirai, H.; Ishiguro, M.; Kambara, T.; Kojima, Y.; Matsunaga, T.; Nishida, H.; Suzuki, Y.; Sugiura, A.; Harwood, H.J.; et al. CJ-15,183, a new inhibitor of squalene synthase produced by a fungus, *Aspergillus aculeatus*. *J. Antibiot.* **2001**, *54*, 904–910. [CrossRef] [PubMed]
7. Zajc, J.; Liu, Y.; Dai, W.; Yang, Z.; Hu, J.; Gostincar, C.; Gunde-Cimerman, N. Genome and transcriptome sequencing of the halophilic fungus *Wallolemia ichthyophaga*: Loadadaptations present and absent. *BMC Genom.* **2013**, *14*, 617. [CrossRef] [PubMed]

8. Kis-Papo, T.; Weig, A.R.; Riley, R.; Peršoh, D.; Salamov, A.; Sun, H.; Lipzen, A.; Wasser, S.P.; Rambold, G.; Grigoriev, I.V.; et al. Genomic adaptations of the halophilic Ded Sea filamentous fungus *Eurotium rubrum*. *Nat. Commun.* **2014**, *5*, 3745. [CrossRef] [PubMed]

9. Wang, Y.; Lu, Z.; Sun, K.; Zhu, W. Effects of high salt stress on secondary metabolite production in the marine-derived fungus *Spicaria elegans*. *Marine Drugs* **2011**, *9*, 535–542. [CrossRef] [PubMed]

10. Wang, Y.; Zheng, J.; Liu, P.; Wang, W.; Zhu, W. Three new compounds from *Aspergillus terreus* PT06-2 grown in a high salt medium. *Marine Drugs* **2011**, *9*, 1368–1378. [CrossRef] [PubMed]

11. Bugni, T.; Ireland, C. Marine-derived fungi: A chemically and biologically diverse group of microorganisms. *Nat. Prod. Rep.* **2004**, *21*, 143–163. [CrossRef] [PubMed]

12. Rateb, M.E.; Ebel, R. Secondary metabolites of fungi from marine habitats. *Nat. Prod. Rep.* **2012**, *28*, 290–344. [CrossRef] [PubMed]

13. Overy, D.P.; Bayman, P.; Kerr, R.G.; Bills, G.F. An assessment of natural product discovery from marine (*sensu strictu*) and marine-derived fungi. *Mycology* **2014**, *5*, 145–167. [CrossRef] [PubMed]

14. Bérdy, J. Thoughts and facts about antibiotics: Where we are now and where we are heading. *J. Antibiot.* **2012**, *65*, 385–395. [CrossRef] [PubMed]

15. Capon, R.J.; Ratnayake, R.; Stewart, M.; Lacey, E.; Tennant, S.; Gill, J.H. Aspergillaxines A-E: Novel heterocyclic dipeptides from an Australian strain of *Aspergillus unilateralis*. *Org. Biomol. Chem.* **2005**, *3*, 123–129. [CrossRef] [PubMed]

16. Wang, G.Y.S.; Arbell, L.M.; Avelar, A.; Borgeson, B.M.; Crews, P. New hisutane based sesquiterpenes from salt water cultures of a marine sponge-derived fungus and the terrestrial fungus *Coriolus consors*. *Tetrahedron* **1998**, *54*, 7335–7342. [CrossRef]

17. The draft genome sequence of *A. aculeatus*. Available online: http://genome.jgi.doe.gov/Aspac1/Aspac1.home.html (accessed on 15 August 2016).

18. De Vries, R.P.; Riley, R.; Wiebenga, A.; Aguilar-Osorio, G.; Amillis, S.; Uchima, C.A.; Anderluh, G.; Asadollahi, M.; Askin, M.; Barry, K.; et al. Comparative genomics reveals high biological diversity and specific adaptations in the industrially and medically important fungal genus *Aspergillus*. *Genome Biol.* **2017**, *18*, 28. [CrossRef] [PubMed]

19. Gaude, E.; Chignola, F.; Spiliotopoulus, D.; Spitaleri, A.; Ghitti, M.; Garcia-Manteiga, J.M.; Mari, S.; Musco, G. *muma*, an R package for metabolomics univariate and multivariate statistical analysis. *Curr. Metab.* **2013**, *1*, 180–189. [CrossRef]

20. Anderson, R.; Buchi, G.; Kobbe, B.; Demain, A.L. Secalonic acids D and F are toxic metabolites of *Aspergillus aculeatus*. *J. Org. Chem.* **1977**, *42*, 352–353. [CrossRef]

21. Gao, Y.-Q.; Guo, C.-J.; Zhang, Q.; Zhou, W.-M.; Wang, C.C.C.; Gao, J.-M. Asperaculanes A and, B.; two sesquiterpenoids from the fungus *Aspergillus aculeatus*. *Molecules* **2015**, *20*, 325–334. [CrossRef] [PubMed]

22. Ingavat, N.; Dobereiner, J.; Wiyakrutta, S.; Mahidol, C.; Ruchirawat, S.; Kittakoop, P. Aspergillusol A., a α-glucosidase inhibitor from the marine-derived fungus *Aspergillus aculeatus*. *J. Nat. Prod.* **2009**, *72*, 2049–2052. [CrossRef] [PubMed]

23. Kurobane, I.; Vining, L.C. Biosynthetic realtionships among the secalonic acids. Isolation of emodin, endocrocin and secalonic acids from *Pyrenochaeta terrestris* and *Aspergillus aculeatus*. *J. Antibiot.* **1979**, *32*, 1256–1266. [CrossRef] [PubMed]

24. Petersen, L.M.; Hoeck, C.; Frisvad, J.C.; Gotfredsen, C.H.; Larsen, T.O. Dereplication guided discovery of secondary metabolites of mixed biosynthetic origin from *Aspergillus aculeatus*. *Molecules* **2014**, *19*, 10898–10921. [CrossRef] [PubMed]

25. Petersen, L.M.; Holm, D.K.; Gotfredsen, C.H.; Mortensen, U.H.; Larsen, T.O. Investigation of a 6-MSA synthase gene cluster in *Aspergillus aculeatus* reveals 6-MSA-derived aculinic acid, aculins A-B and epi-aculin, *A. ChemBioChem* **2015**, *16*, 2200–2204. [CrossRef] [PubMed]

26. Watanabe, S.; Hirai, H.; Ishiguro, M.; Kambara, T.; Kojima, Y.; Matsunaga, T.; Nishida, H.; Suzuki, Y.; Sugiura, A.; Harwood, H.J.; et al. CJ-15,183, a new inhibitor of squalene synthase produced by a fungus, *Aspergillus aculeatus*. *J. Antibiot.* **2001**, *54*, 904–910. [CrossRef] [PubMed]
27. Keller, B.O.; Sui, J.; Young, A.B.; Whittal, R.M. Interferences and contaminants encountered in modern mass spectrometry. *Anal. Chim. Acta* **2008**, *627*, 71–81. [CrossRef] [PubMed]

28. Ryan, M.J.; Smith, D. Fungal genetic resource centres and the genomic challenge. *Mycol. Res.* **2004**, *1088*, 1351–1362. [CrossRef]

29. Kuhls, K.; Lieckfeldt, E.; Börner, T. PCR-fingerprinting used for comparison of ex type strains of *Trichoderma* species deposited in different culture collections. *Microbiol. Res.* **1995**, *150*, 363–371. [CrossRef]

30. Samson, R.A.; Noonim, P.; Meijer, M.; Houbraken, J.; Frisvad, J.C.; Varga, J. Diagnostic tools to identify black aspergilli. *Stud. Mycol.* **2007**, *59*, 129–145. [CrossRef] [PubMed]

31. Bode, H.B.; Bethe, B.; Höfs, R.; Zeeck, A. Big effects from small changes: Possible ways to explore nature's chemical diversity. *ChemBioChem* **2002**, *3*, 619–627. [CrossRef]

32. Bills, G.F.; Platas, G.; Fillola, A.; Jiménez, M.R.; Collado, J.; Vicente, F.; Martín, J.; González, A.; Bur-Zimmermann, J.; Tormo, J.R.; et al. Enhancement of antibiotic and secondary metabolite detection from filamentous fungi by growth on nutritional arrays. *J. Appl. Microbiol.* **2008**, *104*, 1644–1658. [CrossRef] [PubMed]

33. Filtenborg, O.; Frisvad, J.C.; Thrane, U. The significance of yeast extract composition on metabolite production in *Penicillium*. In *Modern Concepts in Penicillium and Aspergillus Classification*; Samson, R.A., Pitt, J.I., Eds.; Plenum Press: New York, NY, USA, 1990; pp. 433–441.

34. Ryan, M.J.; Bridge, P.D.; Smith, D.; Jeffries, P. Phenotypic degeneration occurs during sector formation in *Metarhizium anisopliae*. *J. Appl. Microbiol.* **2002**, *93*, 163–168. [CrossRef] [PubMed]

35. Ryan, M.J.; Smith, D.; Bridge, P.D.; Jeffries, P. The relationship between fungal preservation method and secondary metabolite production in *Metarhizium anisopliae* and *Fusarium oxysporum*. *World J. Microbiol. Biotechnol.* **2003**, *19*, 839–844. [CrossRef]

36. Santos, I.M.; Abrunhosa, L.; Venancio, A.; Lima, N. The effect of culture preservation techniques on patulin and citrinin production by *Penicillium expansum* Link. *Lett. Appl. Microbiol.* **2002**, *35*, 272–275. [CrossRef] [PubMed]

37. Wing, N.; Burgess, L.W.; Bryden, W.L. Cultural degeneration in two Fusarium species and its effects on toxigenicity and cultural morphology. *Mycol. Res.* **1995**, *99*, 615–620. [CrossRef]

38. Wicklow, D.T. Metabolites in the coevolution of fungal chemical defence systems. In *Coevolution of Fungi with Plants and Animals*; Pirozynski, K.A., Hawksworth, D.L., Eds.; Academic Press Limited: London, UK, 1988; pp. 173–201.

39. Anderson, R.F.; Whitmore, L.M.; Brown, W.E.; Peterson, W.H.; Churchill, B.W.; Roegner, F.R.; Campbell, T.H.; Backus, M.P.; Stauffer, J.F. Penicillin production by pigment-free molds. *Ind. Eng. Chem.* **1953**, *45*, 768–773. [CrossRef]

40. Backus, M.P.; Stauffer, J.F. The production and selection of a family of strains in *Penicillium chrysogenum*. *Mycologia* **1955**, *47*, 429–463. [CrossRef]

41. Raper, K.B.; Fennel, D.I. *The Genus Aspergillus*; Williams & Wilkins Co.; Baltimore, MD, USA, 1965; pp. 1–686.

42. Yue, Q.; Chen, L.; Zhang, X.; Li, K.; Sun, J.; Liu, X.; An, Z.; Bills, G.F. Evolution of chemical diversity in the echinocandin lipopeptide antifungal metabolites. *Eukaryot. Cell* **2015**, *14*, 698–718. [CrossRef] [PubMed]

43. Jennings, D.H. Fungal growth in the sea. In *The Biology of Marine Fungi*; Moss, S.T., Ed.; Cambridge University Press: Cambridge, UK, 1986; pp. 1–10.

44. Liu, S.; Li, J.; Wu, Y.; Ren, Y.; Liu, Q.; Wang, Q.; Zhou, X.; Cai, M.; Zhang, Y. De novo transcriptome sequencing of marine-derived *Aspergillus glaucus* and comparative analysis of metabolic and developmental variations in response to salt stress. *Genes Genom.* **2017**, *39*, 317–329. [CrossRef]

<figure><img src="image_13.png" /></figure>
27. Keller, B.O.; Sui, J.; Young, A.B.; Whittal, R.M. Interferences and contaminants encountered in modern mass spectrometry. *Anal. Chim. Acta* **2008**, *627*, 71–81. [CrossRef] [PubMed]

28. Ryan, M.J.; Smith, D. Fungal genetic resource centres and the genomic challenge. *Mycol. Res.* **2004**, *1088*, 1351–1362. [CrossRef]

29. Kuhls, K.; Lieckfeldt, E.; Börner, T. PCR-fingerprinting used for comparison of ex type strains of *Trichoderma* species deposited in different culture collections. *Microbiol. Res.* **1995**, *150*, 363–371. [CrossRef]

30. Samson, R.A.; Noonim, P.; Meijer, M.; Houbraken, J.; Frisvad, J.C.; Varga, J. Diagnostic tools to identify black aspergilli. *Stud. Mycol.* **2007**, *59*, 129–145. [CrossRef] [PubMed]

31. Bode, H.B.; Bethe, B.; Höfs, R.; Zeeck, A. Big effects from small changes: Possible ways to explore nature's chemical diversity. *ChemBioChem* **2002**, *3*, 619–627. [CrossRef]

32. Bills, G.F.; Platas, G.; Fillola, A.; Jiménez, M.R.; Collado, J.; Vicente, F.; Martín, J.; González, A.; Bur-Zimmermann, J.; Tormo, J.R.; et al. Enhancement of antibiotic and secondary metabolite detection from filamentous fungi by growth on nutritional arrays. *J. Appl. Microbiol.* **2008**, *104*, 1644–1658. [CrossRef] [PubMed]

33. Filtenborg, O.; Frisvad, J.C.; Thrane, U. The significance of yeast extract composition on metabolite production in *Penicillium*. In *Modern Concepts in Penicillium and Aspergillus Classification*; Samson, R.A., Pitt, J.I., Eds.; Plenum Press: New York, NY, USA, 1990; pp. 433–441.

34. Ryan, M.J.; Bridge, P.D.; Smith, D.; Jeffries, P. Phenotypic degeneration occurs during sector formation in *Metarhizium anisopliae*. *J. Appl. Microbiol.* **2002**, *93*, 163–168. [CrossRef] [PubMed]

35. Ryan, M.J.; Smith, D.; Bridge, P.D.; Jeffries, P. The relationship between fungal preservation method and secondary metabolite production in *Metarhizium anisopliae* and *Fusarium oxysporum*. *World J. Microbiol. Biotechnol.* **2003**, *19*, 839–844. [CrossRef]

36. Santos, I.M.; Abrunhosa, L.; Venancio, A.; Lima, N. The effect of culture preservation techniques on patulin and citrinin production by *Penicillium expansum* Link. *Lett. Appl. Microbiol.* **2002**, *35*, 272–275. [CrossRef] [PubMed]

37. Wing, N.; Burgess, L.W.; Bryden, W.L. Cultural degeneration in two Fusarium species and its effects on toxigenicity and cultural morphology. *Mycol. Res.* **1995**, *99*, 615–620. [CrossRef]

38. Wicklow, D.T. Metabolites in the coevolution of fungal chemical defence systems. In *Coevolution of Fungi with Plants and Animals*; Pirozynski, K.A., Hawksworth, D.L., Eds.; Academic Press Limited: London, UK, 1988; pp. 173–201.

39. Anderson, R.F.; Whitmore, L.M.; Brown, W.E.; Peterson, W.H.; Churchill, B.W.; Roegner, F.R.; Campbell, T.H.; Backus, M.P.; Stauffer, J.F. Penicillin production by pigment-free molds. *Ind. Eng. Chem.* **1953**, *45*, 768–773. [CrossRef]

40. Backus, M.P.; Stauffer, J.F. The production and selection of a family of strains in *Penicillium chrysogenum*. *Mycologia* **1955**, *47*, 429–463. [CrossRef]

41. Raper, K.B.; Fennel, D.I. *The Genus Aspergillus*; Williams & Wilkins Co.; Baltimore, MD, USA, 1965; pp. 1–686.

42. Yue, Q.; Chen, L.; Zhang, X.; Li, K.; Sun, J.; Liu, X.; An, Z.; Bills, G.F. Evolution of chemical diversity in the echinocandin lipopeptide antifungal metabolites. *Eukaryot. Cell* **2015**, *14*, 698–718. [CrossRef] [PubMed]

43. Jennings, D.H. Fungal growth in the sea. In *The Biology of Marine Fungi*; Moss, S.T., Ed.; Cambridge University Press: Cambridge, UK, 1986; pp. 1–10.

44. Liu, S.; Li, J.; Wu, Y.; Ren, Y.; Liu, Q.; Wang, Q.; Zhou, X.; Cai, M.; Zhang, Y. De novo transcriptome sequencing of marine-derived *Aspergillus glaucus* and comparative analysis of metabolic and developmental variations in response to salt stress. *Genes Genom.* **2017**, *39*, 317–329. [CrossRef]

<figure><img src="image_14.png" /></figure>

© 2017 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).