Next Article in Journal
Evaluation of Oxasqualenoids from the Red Alga Laurencia viridis against Acanthamoeba
Next Article in Special Issue
Discovery of Primarolides A and B from Marine Fungus Asteromyces cruciatus Using Osmotic Stress and Treatment with Suberoylanilide Hydroxamic Acid
Previous Article in Journal
Costatone C—A New Halogenated Monoterpene from the New Zealand Red Alga Plocamium angustum
Previous Article in Special Issue
Influence of OSMAC-Based Cultivation in Metabolome and Anticancer Activity of Fungi Associated with the Brown Alga Fucus vesiculosus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rapid Metabolome and Bioactivity Profiling of Fungi Associated with the Leaf and Rhizosphere of the Baltic Seagrass Zostera marina

1
GEOMAR Centre for Marine Biotechnology (GEOMAR-Biotech), Research Unit Marine Natural Products Chemistry, GEOMAR Helmholtz Centre for Ocean Research Kiel, Am Kiel-Kanal 44, 24106 Kiel, Germany
2
Faculty of Mathematics and Natural Sciences, Kiel University, Christian-Albrechts-Platz 4, 24118 Kiel, Germany
*
Author to whom correspondence should be addressed.
Mar. Drugs 2019, 17(7), 419; https://doi.org/10.3390/md17070419
Submission received: 29 May 2019 / Revised: 15 July 2019 / Accepted: 16 July 2019 / Published: 19 July 2019

Abstract

:
Zostera marina (eelgrass) is a marine foundation species with key ecological roles in coastal habitats. Its bacterial microbiota has been well studied, but very little is known about its mycobiome. In this study, we have isolated and identified 13 fungal strains, dominated by Penicillium species (10 strains), from the leaf and the root rhizosphere of Baltic Z. marina. The organic extracts of the fungi that were cultured by an OSMAC (One-Strain–Many-Compounds) regime using five liquid culture media under both static and shaking conditions were investigated for their chemical and bioactivity profiles. All extracts showed strong anti-quorum sensing activity, and the majority of the Penicillium extracts displayed antimicrobial or anti-biofilm activity against Gram-negative environmental marine and human pathogens. HPLC-DAD-MS-based rapid metabolome analyses of the extracts indicated the high influence of culture conditions on the secondary metabolite (SM) profiles. Among 69 compounds detected in all Penicillium sp. extracts, 46 were successfully dereplicated. Analysis of SM relatedness in culture conditions by Hierarchical Cluster Analysis (HCA) revealed generally low similarity and showed a strong effect of medium selection on chemical profiles of Penicillium sp. This is the first study assessing both the metabolite and bioactivity profile of the fungi associated with Baltic eelgrass Z. marina.

Graphical Abstract

1. Introduction

The common eelgrass Zostera marina is a marine flowering plant, which, similar to all seagrasses, evolved from sea to land plants but returned to the sea about 140 million years ago [1]. Today, well-adapted to sea environment, eelgrass meadows dominate the soft sediment shorelines of the temperate northern hemisphere, including Atlantic and Pacific Oceans. Due to its high tolerance to low and variable salinity levels, Z. marina represents the most common seagrass in the brackish Baltic Sea [2]. As a foundation species and a key ecosystem engineer, eelgrass plays important ecological roles in primary production, food web support, production of oxygen, organic matter degradation, and recycling of nutrients in coastal environments. Zostera beds also provide nursery habitat for many animals (juvenile fish, invertebrates, and other larger animals) and stabilize sediments, preventing coastlines from erosion [3]. The recently completed genome-sequencing project [4] showed that, as part of its re-adaptation process into the marine environment, Z. marina had gained genes enabling its tolerance to saltwater. However, it has also lost genes characteristic for terrestrial angiosperms, leading to reduced chemical machinery (e.g., no volatile terpene synthesis), chemical defense, and communication [5]. Nonetheless, a generally low rate of microfouling and biofilm formation in seagrasses [6,7,8] point out defense mechanisms and their ability to control attachment of (micro)foulers.
It is widely accepted that the microbiome is involved in the maintenance of general health and overall performance of plants. The microbes colonizing the above- or below-ground plant tissues are crucial in fitness, disease resistance, growth, and survival of the host plant [9]. Similar to other angiosperms, Z. marina has intimate relationships with both beneficial and harmful microbes in its surrounding environment. Studies on different parts of Z. marina suggest their distinct microbial community [2,10,11]. Particularly, the sediment surrounding the roots of eelgrass has been shown to have a unique and rich microflora [5,11,12]. However, almost all these studies have focused on the bacterial microbiota associated with Z. marina. Being part of the whole microbiome, the fungi, which are talented producers of antibiotics and many other secondary metabolites (SMs) with diverse bioactivities, may also contribute to overall fitness and chemical defense of any plant/marine organisms. Only a recent study by Kirichuk and Pivkin [13] reported the diversity of cultivable filamentous fungi associated with different parts of Z. marina collected in the Sea of Japan (Russia). A few filamentous fungi obtained from Z. marina leaves or rhizomes have been studied for isolation of their chemical constituents [14,15]. To our knowledge, no systematic study has been performed on the chemical profiles of cultivable mycobiota of Z. marina. In this context, marine fungi are considered not a systematical and taxonomical group, but rather a physiological and ecological, including true marine, as well as facultative strains [16].
In this study, we aimed to identify the cultivable mycoflora of the leaf and sediment surrounding the eelgrass roots and investigated their chemical and bioactivity profiles. We used an OSMAC (One-Strain–Many-Compounds)-based culture approach, which led to the isolation and identification of 13 fungal strains from the phyllosphere and the rhizosphere of Zostera marina collected from Kiel Fjord (Baltic Sea, Germany). The organic extracts of all fungi that were grown in five different liquid media under both static and shaking conditions were rapidly investigated for their secondary metabolome (using LC-DAD-MS) and screened for their ecologically as well as medically relevant bioactivities (antibacterial, anti-quorum sensing, and anti-biofilm formation).

2. Results

2.1. Identification and Extraction of Marine Fungi

A total of 13 fungal strains, which were identified based on morphology and ITS sequencing, were isolated. Only one fungus, identified as Phoma macrostoma, was obtained from the leaves of Z. marina, while the sediment surrounding eelgrass roots afforded 12 fungal strains, 10 of which were successfully identified at species level with 99 and 100% sequence similarities to known species (Table 1). This includes one representative of the genus Cladosporium (Cladosporium langeronii), one strain of the genus Trichoderma (Trichoderma harzianum), and 10 Penicillium spp., which were assigned to three different species; as Penicillium antarcticum (four isolates), Penicillium atramentosum (one isolate), Penicillium atrovenetum (five isolates). Although different media were used for the isolation campaign, there was no medium selectively increasing the number of obtained strains.

2.2. Extraction, Metabolite, and Bioactivity Profile of Fungi Associated with Zostera marina

Culture conditions significantly influence SM profiles of microorganisms. To gain a wide chemical diversity, the isolated fungi were grown in five liquid media differing in their carbon and nitrogen sources and salinity (PC (Potato-Carrot), CA (Casamino Acids), WH (Wickerham15), GPY (Glucose-Peptone-Yeast), and AS (Artificial Seawater)) under both static and shaking (“sh”) conditions. The ethyl acetate (EtOAc) extracts of these cultures were submitted to LC-DAD-MS analyses for rapid chemical profiling. All extracts were tested (at 100 µg/mL concentration) for antibacterial activity against two marine Gram-negative environmental pathogens (the fish pathogen, Vibrio anguillarum, and the algal pathogen, Pseudomonas elyakovii), the human pathogen Pseudomonas aeruginosa, as well as Aliivibrio fischeri, which frequently occurs in marine environments [17,18]. The extracts were also evaluated for their biofilm inhibition and anti-quorum sensing efficacy against Ps. aeruginosa and A. fischeri, respectively. Comparison of antibiotic and anti-QS (quorum sensing)/anti-biofilm activity allows discrimination between these activities for later development of these compounds for different applications, such as antibiotic drugs and biofilm preventing agents, e.g., on medical devices or other sensitive surfaces.
The results of the chemical and bioactivity profiling of non-Penicillium strains are presented individually (Figure 1, Figure 2 and Figure 3, Table 2 and Table 3), while the chemical profiles of 10 Penicillii are outlined together to allow comparative analyses (Figure 4, Figure 5 and Figure 6). Only a few selected and striking media profiles are displayed as Figures in the main text, while the full chemical profiles of all culture conditions can be found in the Supplementary file. The abbreviations, e.g., CA-1-sh used in Figures (sometimes also in the text) refer to the medium used for culturing (e.g., CA), strain number of the fungi (1), and the application of shaking (sh) in the culture regime.

2.2.1. Secondary Metabolite Profile of Phoma macrostoma (Strain 1)

Phoma macrostoma was the only fungus isolated from the eelgrass phyllosphere and was the single representative of the genus Phoma. Chemical profile of this strain showed significant variations within different media and culture conditions (Figure 1). Table 2 shows the list of compounds detected in all media and culture regimes (for all non-Penicillium species, strains 1–3), containing information on their UV absorption, molecular weights, and putative identity. The extracts of both Potato-Carrot (PC) and Wickerham15 (WH) static cultures (shown as PC-1 and WH-1, 1 refers to strain number) contained the unknown peaks UP2 (UVmax 220, 275, 325, 360 nm; m/z 273.1 [M + H]+), UP3 (UVmax 218, 296, 364 nm; m/z 289.1 [M + H]+), and compound 1 (UVmax 220, 245, 292, 342 nm; m/z 229.1 [M + H]+). Compound 1 was dereplicated as the furanocoumarin (psoralen) derivative trioxsalen (Table 2). Little to no fungal growth was observed in the statically grown Casamino Acids (CA), Artificial Seawater (AS), and Glucose-Peptone-Yeast (GPY) cultures; hence, their extracts had very poor chemistry (Figure 1, Supplementary Figure S1). In contrast, several compounds were detected in cultures grown under shaking (sh) conditions. The PC-1 medium extract contained three peaks (UP23 and 1), while shaking increased the number of detectable compounds in PC-1-sh to four (UP13 and 1), stimulating the production of another unknown compound UP1 (UVmax 221, 292 nm, m/z 325.2 [M + H]+). In CA-1-sh medium, shaking induced the biosynthesis of UP4 (UVmax 232 nm; m/z 279.3 [M + H]+) and 1. The strongest influence of shaking was observed in the WH medium. The WH-1-sh extract gained two compounds UP4 and 2, which was putatively dereplicated as antibiotic NK-A 17E-233II (UVmax 222, 279 nm; m/z 277.2 [M + H]+) but lacked UP2 (Figure 1). In general, shaking increased the number of detected compounds in strain 1, and PC-1-sh and WH-1-sh cultures exhibited the most diverse chemical composition.
None of the Ph. macrostoma extracts had inhibitory activity on environmental pathogens, V. anguillarum and Ps. elyakovii. Extracts obtained from PC, PC-sh, CA, AS, AS-sh, WH, and WH-sh cultures showed quorum sensing inhibitory (QSI, from 53 to 100% at 100 μg/mL concentration) activity against A. fischeri. The extracts of PC and WH slightly (>20%) inhibited the growth of Ps. aeruginosa, but not the biofilm formation by this Gram-negative bacterium (Table 3).

2.2.2. Secondary Metabolite Profile of Cladosporium langeronii (Strain 2)

Cladosporium langeronii was isolated from the rhizosphere of Z. marina. As shown in Figure 2, only slight differences were visible in its chemical profile upon the use of different media and culture regimes. Two known compounds (3 and 4) were detected in the CA culture extracts of C. langeronii, under both static and shaking conditions (Figure 2). The other media were not favored; hence, no compound was detected in the other four media extracts (Supplementary Figure S2). Compound 3 (UVmax 221, 289, 310 nm; m/z 211.1 [M + H]+) was annotated as the pyranacetal agistatin D, while compound 4 (UVmax 205, 216, 222, 277 nm; m/z 295.1 [M + H]+) was annotated as the orsellinic acid derivative orsellide D (Table 2). Compared to the static CA extracts, compounds 3 and 4 exhibited higher intensities of UV absorption in the CA-2-sh extracts. In general, shaking led minor changes in chemical profiles, and the CA medium appeared to be the best medium for this fungus.
The crude extracts of C. langeronii had no inhibitory activity on the growth of V. anguillarum, Ps. elyakovii, Ps. aeruginosa, or the biofilm formation by the latter. Notably, PC, PC-sh, CA, AS-sh, WH, and WH-sh culture extracts strongly inhibited the QS of A. fischeri (77%–100% at 100 μg/mL concentration), while the growth of this microorganism was not inhibited (Table 3).

2.2.3. Secondary Metabolite Profile of Trichoderma harzianum (Strain 3)

Trichoderma harzianum was the second strain isolated from the rhizosphere of Z. marina. The metabolite profile of this strain possessed distinct diversity within different media and culture regimes (Figure 3, Supplementary Figure S3, Table 2). The static cultures of strain 3 produced altogether nine compounds (UP58 and 59). With seven detected metabolites, statically grown WH culture yielded the chemically richest extracts. It contained four unknown compounds, UP5 (UVmax 215, 303 nm, m/z 327.3 [M + H]+), UP6 (UVmax 218, 297 nm, m/z 312.2 [M + H]+), UP7 (UVmax 220, 266, 276, 287 nm; m/z 393.1 [M + H]+), and UP8 (UVmax 226 nm; m/z 245.2 [M + H]+). The compounds 5 (UVmax 215, 306 nm; m/z 311.3 [M + H]+), 6 (UVmax 223, 275, 288 nm; m/z 599.3 [M + H]+), and 9 (UVmax 223, 261, 271, 281; m/z 282.2 [M + H]+) were putatively identified as orsellide C, the cyclotetrapeptide cyl-2, and 2′-deoxyribofuranosylguanine (2N-Me), respectively. Three peaks were detected in static CA-3 cultures (68, Table 2). Compounds 7 (UVmax 223, 280 nm; m/z 316.1 [M + H]+) and 8 (UVmax 227, 269, 334 nm; m/z 282.1 [M + H]+) were dereplicated as quinolinone B and harzianopyridone alkaloids, respectively. The only detectable compound in the extracts of the static PC-3 and the shaken CA-3-sh cultures was compound 9. No peaks were detected in the statically grown AS and GPY cultures. In general, shaking suppressed the number of compounds produced by strain 3. Strikingly, the overall number of detected peaks in WH cultures declined from nine to two upon shaking (UP9 and 9), while, an unidentified compound UP9 (UVmax 222, 280, 292, 305, 319 nm; m/z 626.2 [M + H]+) and 9 were visible in WH-3-sh extract. Overall, static culture regime using WH medium provided the most favored conditions for yielding the richest secondary metabolite profile (SMP) of T. harzianum.
Several Trichoderma harzianum crude extracts, namely PC, CA, CA-sh, AS-sh, WH, and WH-sh, showed QSI activity against A. fischeri (70% to 97% inhibitory capacity). Also, the growth of A. fischeri was inhibited by extracts derived from the PC (70%) and WH (46%) media cultures of T. harzianum (Table 3).

2.2.4. Secondary Metabolite Profile of Penicillium spp. (Strains 4–13)

Due to their abundance, the chemical profiles of all 10 Penicillium strains associated with the rhizosphere of Z. marina were explored together. Their SM profiles (cultivated in five different liquid media under both shaking and static conditions) were represented in a comparative manner rather than individually. Supplementary Table S1 contains detailed information on UV absorptions, molecular weights, and putative identities of compounds detected in Penicillium extracts in all media and culture regimes. Among 69 compounds detected in Penicillium extracts in total, 46 were assigned to known compounds (Supplementary Figure S14 and Supplementary Table S1). Figure 4 illustrates the total peak number (TPN), as well as the number of putatively dereplicated compounds, for each Penicillium strain fermented under both static and shaking conditions. These numbers represent the sum of detected compounds versus the sum of annotated compounds in all five media.
Figure 4 shows that the most productive Penicillium strain under static conditions (all peaks detected in all static cultures counted together) was strain 7 (P. antarcticum) with 36 detected peaks in all media, followed by strain 12 (32 peaks), 9 (28 peaks), 10 (28 peaks), 5 (27 peaks), 13 (27 peaks), 6 (21 peaks), 4 (20 peaks), 11 (19 peaks), and 8 (15 peaks). In most cases, approximately one-quarter of all detected peaks in extracts of static cultures remained unknown, while three-quarters were dereplicated and assigned to known compounds (see Supplementary Table S1). In general, extracts deriving from statically grown WH and CA cultures were found to be richer in SM than those grown statically in PC, AS, and GPY cultures (Figure 4). The highest TPN in static WH medium cultures was observed in the extract of strain 9, exhibiting 18 detected peaks. Also, strains 7, 10, 11, and 12 produced many compounds in static WH medium. In statically incubated CA cultures, the highest TPNs were detected in strains 7 and 13, both containing 13 peaks. In summary, under static conditions, WH and CA media provide the richest metabolite profiles for Penicillium strains (Figure 4). In contrast, PC, AS, and GPY media yield poor SM profiles in static cultures. Notably, strains 7 and 12 produced high numbers of compounds in all media, including the nutrient-poor AS medium, under static conditions.
Under shaking conditions, the most productive Penicillium sp. was strain 6 (P. antarcticum), which contained altogether 26 peaks detected in all media extracts. This strain was followed by strain 13 (25 peaks), 5 (23 peaks), 9 (21 peaks), 12 (21 peaks), 4 (20 peaks), 8 (20 peaks), 7 (18 peaks), 10 (17 peaks), and 11 (15 peaks). Most compounds found in the extracts of shaking cultures were assigned to known compounds, especially those in the extracts of CA and WH cultures (e.g., strains 6, 7, 8, 10, and 13). However, many compounds remained unknown in the shaking culture extracts of AS medium (e.g., strains 5, 6, 9, 11, and 12). In comparison to static conditions, shaking induced the chemical diversity of the AS and PC cultures, while it was mostly vice versa in the CA and WH cultures. The best performer in shaking AS cultures was strain 6 with a TPN of 11, five of which were annotated to known compounds (Supplementary Table S1). Strains 5, 9, 11, and 12 possessed rich SMPs with many unknown compounds in AS-sh cultures. The highest TPN in shaking PC cultures was found in the extracts of strain 12, possessing seven peaks. In contrast to AS-sh cultures, the peaks in PC-sh cultures were assigned to known compounds (Supplementary Table S1). Except for strains 7, 8, 9, and 13, shaking led to an unfavorable effect in chemical diversity of GPY culture extracts. These findings show that the effect of aeration (shaking) on the SMP of most Penicillium strains is medium-dependent. However, most strains exhibited overall similar sums of TPNs, except for strains 7, 8, and 12.
Although identified as same species, strains 4–7 (P. antarcticum), strain 8 (P. atramentosum), and strains 9–13 (P. atrovenetum, see Table 1) did strongly vary in their SMPs, and their response to stimuli through culture conditions was strain-dependent (Supplementary Figures S4–S13). Hierarchical cluster analysis (HCA) based on Bray-Curtis similarity was applied to compare the similarity of the chemical profiles of all Penicillium strains (based on presence/absence of compounds in different culture extracts) grown in five different media under static and shaking regimes. In this analysis, a similarity value of 1.0 indicates that two (or more) strains have 100% identical SMPs (producing a horizontal line at the top), while a similarity value of 0 shows that their SMPs are completely unrelated (0% similarity, producing a horizontal line at the bottom). Similarly, a value between 0 and one implies partial relatedness among SMPs of the respective strains. Each branch in a dendrogram represents one SMP of a specific strain (or more, if similar) cultivated under two culture regimes. Therefore, strongly branched dendrograms indicate a high variety of SMPs of respective strains in one media and low chemical relatedness between the Penicillium strains. In contrast, weakly branched dendrograms indicate vice versa. Hence, HCA allowed quick assessment of the impact of culture media and culture regimes in a color-coded manner (statically grown cultures in blue, shaken in red) on the chemical diversity of all Penicillium strains. The resulting dendrograms are illustrated in Figure 5. In general, cultivation in CA and WH media (Figure 5B,D) resulted in more complex dendrograms than those in PC, AS, and GPY media. Clusters of the CA and WH cultures were considerably branched, and only small groups of identical SMPs were formed (Figure 5B,D). In contrast, clusters of PC, AS, and GPY cultures exhibited fewer branches but contained more and larger groups of identical similarity. This indicates that SMPs of PC, AS, and GPY cultures are not as diverse as the CA or WH cultures and more sensitive to the culture regime employed. Notably, the strain 8 (P. atramentosum) appeared to be chemically distant to all other Penicillium strains in all five media and both culture regimes, while the remaining Penicillium strains shared similarities of at least 0.05 (Figure 5E).
Specifically, in the PC medium (Figure 5A), the similarity between the static (8) and shaken (8-sh) cultures of strain 8 was approx. 0.65, i.e., showed 65% similarity. However, SMP similarity of these cultures with other Penicillium strains was 0. As shown in Figure 5A, strains 11, 6, 9, 10, 7, and 12 formed a large cluster with a similarity value of 1.0 in static PC cultures. A second cluster with 100% similarity was obvious between strains cultured under both static and shaking conditions (11-sh, 4, 5, 5-sh, 9-sh, 13). Both clusters had a similarity value of approx. 0.65. A third small cluster in Figure 5A with identical SMP was formed by 7-sh and 13-sh cultures. All other strains showed similarities between 0.25 and 0.85, indicating major differences but also overlaps in the SMPs.
The dendrogram of all static and shaking Penicillium CA culture extracts (Figure 5B) was highly branched, indicative of high chemical diversity. As shown in Figure 5B, four small groups with a 100% similarity were formed, each consisting of two representatives (10-sh/12-sh; 9-sh/7-sh; 11/11-sh; 8/8-sh). For strain 8, the static and shaking conditions resulted in an identical SMP, indicating that the culture regime had no impact on its SMP in CA medium. Again, the SMP of strain 8 was chemically unrelated (similarity of 0) to all other Penicillii. The same applies to the SMP of strain 11 in CA medium. The SMPs of all other strains were highly clustered and showed similarities between 0.2 and 0.9. Therefore, some compounds were common in all Penicillium isolates in CA medium, and others were exclusive to specific strains.
In the AS medium (Figure 5C), culture regime had a strong impact on the SMP of strain 8, with 0% similarity. Both 8 and 8-sh cultures were yet unrelated to the remaining Penicillium strains. Statically grown strains 11, 4, 5, 6, 9, and 10 formed a big cluster with 100% similarity. The same applies to strains 10-sh and 13. The other strains exhibited similarity values ranging from 0.35 to 0.95 (Figure 5C), indicating increased variability in their SMPs upon shaking.
The dendrogram of WH cultures (Figure 5D) was highly branched. It contained many clusters and exhibited only one small cluster, consisting of strains 7-sh and 12-sh, with 100% identical SMP. The strain 8 was again chemically unrelated to all other Penicillium spp. but shared approx. 60% SMP similarity between static and shaking cultures. Most of the strains showed similarity values between 0.15 and 0.9 in different culture regimes, indicating high chemical diversity.
In the GPY medium (Figure 5E), the static and shaking culture regimes led to a 65% SMP similarity in strain 8, that was again 100% distant from the remaining Penicillium strains. A few clusters with 100% similarity were visible, i.e., 4/5; 5-sh/6-sh, 9-sh/10-sh/7-sh/12-shs and 11-sh/4-sh/7s in the dendrogram (Figure 5E). The latter cluster, as well as the 8/8-sh cultures, had no chemical relatedness to other Penicillii. The remaining strains formed two big sub-clusters with a similarity value of 0.05. These results highlight the impact of culture regime on the chemical diversity of Penicillium strains in GPY medium.

2.2.5. Bioactivity of Penicillium Species

The crude extracts of all Penicillium isolates were screened against the same panel of microorganisms for their antibacterial activity, as well as their inhibitory effects on quorum sensing (QSI) and biofilm formation. To allow clear and quick assessment of the bioactivity profile of different Penicillium spp., their bioactivities were displayed using matrix plots (Figure 6). Inhibition rates (%, in comparison to standard antibiotics) of fungal extracts are displayed as a function of color intensity ranging from blue (low inhibition) to red (high inhibition). The matrix plots of bioactivities revealed similar patterns of bioactivity against V. anguillarum and Ps. elyakovii for all Penicillii (Figure 6A,B). Almost all PC, PC-sh, AS, and AS-sh media extracts, except for those of strain 8, potently inhibited V. anguillarum and Ps. elyakovii. The crude extracts of WH, CA, and GPY cultures (static and shaking) were inactive against both pathogens. The only exceptions were the strain 9 (in GPY medium) and strains 4, 5, and 13 (in CA medium, Figure 6A,B), all incubated in static conditions. Figure 6C,D displays the activity of all Penicillii extracts for QSI and antibacterial activity against A. fischeri, respectively, as matrix plots. The majority of the crude extracts displayed QSI activity, while only some inhibited the growth of A. fischeri. Similar to the activity observed against V. anguillarum and Ps. elyakovii, almost all PC/PC-sh and AS/AS-sh culture extracts exhibited potent anti-quorum sensing activity. The static cultures of a few additional strains cultured in CA (e.g., strains 4, 5, and 13), CA-sh (strains 8, 13), WH medium (strains 4, 5, 9, 10, and 13), GPY (strains 4, 5, 9–13), and GPY-sh (strains 9 and 13) also showed high QSI potential (Figure 6C). Several extracts obtained from the static AS media (4, 5, 6, 9, 11, and 13) and few extracts obtained under shaking conditions, i.e., strains 9, 11 (PC-sh), and 10 (PC-sh and AS-sh), inhibited the growth of A. fischeri. None of the CA and CA-sh extracts showed reasonable inhibition of A. fischeri.
The bioactivity of Penicillium strains against the growth and biofilm formation by Ps. aeruginosa is depicted in Figure 6E,F, respectively. The highest anti-biofilm activity against Ps. aeruginosa was exerted by four strains grown in static conditions in PC medium (4, 5, 6, 10, and 11) and two strains (7, 13) grown in PC-sh (Figure 6E). The strains 11 and 12 grown in both AS and AS-sh media had high anti-biofilm potential (Figure 6E). Almost all these extracts, as well as strains 4 (in AS), 10 and 11 (in PC), and 12 and 13 (in PC-sh), displayed antibacterial activity against this Gram-negative bacterium (Figure 6F). In comparison to former strains 1–3, Penicillium spp. overall appeared to have higher bioactivity potential, not only showing QSI and anti-biofilm formation but also antibacterial activity against Gram-negative bacteria.

3. Discussions

3.1. Fungi Associated with the Leaf and Rhizosphere of Baltic Z. marina

In continuation of our research interest into Z. marina from the Baltic Sea [6,8], we have investigated fungal communities associated with its phyllosphere and rhizosphere. Dissected eelgrass specimens and seagrass rhizosphere sediment samples led to the isolation of 13 fungi, belonging to four different genera and six different species. Interestingly, only one isolate (Phoma macrostoma) derived from the Zostera leaves. This is lower than that reported in previous studies. Newell (1981), who incubated the green and brown (dead) leaves of Z. marina from Chesapeake Bay (USA) in seawater without an additional carbon source identified fungi belonging to six different genera (by microscopic observation) as well as living mycelia without identifiable propagules [27]. The only detailed mycobiome study performed on Z. marina is from Rifovaya Bay (Sea of Japan) by Kirichuk and Pivkin (2015) that reported 16 fungi from the phylloplane (specifically the leaves, without stems, flowers, or fruits). They belong to diverse genera, i.e., Acremonium, Arthrinium, Aspergillus, Beauveria, Cladosporium, Trichoderma, and Penicillium, the latter being the most abundant one (seven isolates). To our knowledge, no previous study reported the isolation of a Phoma strain from Z. marina leaves. The discrepancy between our work and that of Kirichuk and Pivkin (2015) may stem from potential temporal variations (the different time/season for sampling/isolation of the fungi that was unfortunately not stated), spatial differences (different collection site, i.e., Rifovaya Bay, Sea of Japan versus Baltic Sea), as well as their larger sample size (four stations) [13]. In the present study, the seagrass rhizosphere afforded 12 filamentous fungi belonging to genera Cladosporium, Trichoderma, and Penicillium. With 10 isolates, the genus Penicillium dominated the sediment mycobiome. The previous work by Kirichuk and Pivkin (2015) reported similar numbers (16) from the eelgrass rhizosphere, also predominantly composed of Penicillium spp. [13]. Notably, the same authors have inoculated their samples on Agar wort (prepared with seawater containing unhopped beer wort) or directly on sterile filter papers [13]. Herein, we used five different agar media in Petri dishes, which is regarded as the most suitable inoculation method for fungal strain isolation [28]. The other experimental details (e.g., potential dilution, temperature, total inoculation/incubation time, identification other than morphological traits) are absent in their report to enable comparisons. The common finding in these two studies is the dominance of Penicillium species in rhizosphere samples.
It has been shown that the microbiome (bacterial communities) of the surface (epiphytic), inner tissues (endophytic), and sediments surrounding the roots (rhizopheric) of seagrasses are quite distinct. In particular, the sediment microbiome has been reported to represent a unique and dynamic microsystem [5,11,12,29]. Seagrasses grow in anoxic and sulfidic sediments as a result of deposition of large amounts of organic material in seagrass meadows. The plant-driven gradients of oxygen and dissolved organic carbon determine the distinct microbiome of the rhizosphere [30,31,32]. The availability of various electron acceptors (sulfate) and the exudates released by the seagrass roots promote bacterial growth and microbial processes, hence, allowing a unique niche of microecosystem with complex microbial assemblages around the rhizosphere. Although such high carbon stocks may represent a suitable environment for fungal growth, very little is known on fungal assemblages in Zostera rhizosphere (or in the other parts of Zostera spp.) or how they modulate the plant microbe-interactions, thereby possibly affecting the health of the plant. Furthermore, no attempt has been done for chemical and bioactivity profiling on endophytic and rhizospheric fungal communities of Z. marina. This has been the main aim of this study.

3.2. Dereplication of Secondary Metabolite Profiles of Z. marina Associated Fungi

Despite the advancement of modern analytical instruments and dereplication methods, the identification of known compounds in complex natural extracts remains a challenging endeavor. To decipher the chemical potential of Z. marina associated fungi, we herein followed an OSMAC approach using different media and culture regimes to stimulate the production of putatively novel SMs [33]. For a rapid assessment of the chemical diversity produced in Z. marina associated fungi, we compared information gained from HPLC-DAD-MS measurements (UV/VIS absorption and molecular mass) with entries in large public and commercial databases to conduct dereplication of known compounds.
Strain 1 was genetically identified as Phoma macrostoma. The genus Phoma has never been isolated from Z. marina before, although Phoma spp. are common in Nature, including terrestrial and marine hosts and marine (Arctic) sediments [34,35,36]. Chemically, the crude extracts of this strain displayed moderate chemical diversity. Analyses of the SMP of this strain showed PC and WH media to be the most suitable. The fungal growth, in the static CA, AS, and GPY cultures, was poor, so was the chemical profile. SMPs of the crude PC and PC-sh culture extracts looked alike except for the presence of UP1 with a weak chromophore in PC-sh. UP2 and UP3 occurred in static PC and WH cultures, and the same applied for UP4 in shaken CA and WH cultures. Two known compounds were identified as trioxsalen (1) and “NK-A 17E-233II” (2). Trioxsalen (aka 4,5′,8-trimethylpsoralen) is a phototoxic furanocoumarin originally isolated from domestic celery infected with the phytopathogenic fungus Sclerotinia sclerotium [19]. It is used for the phototherapeutic treatment of vitiligo [37] and is cytotoxic against human melanoma cells [38]. The second compound “NK-A 17E-233II” is obtained from a marine Phoma sp. [20] with reported antitumor activity. In the present study, the extracts of PC, PC-sh, WH, WH-sh, and CA-sh cultures exhibited high QSI activity (Table 3).
The fungal genus Cladosporium has been reported from Z. marina by Kirichuk and Pivkin [13] and Newell [27]. In general, Cladosporium spp. are widespread and have been isolated from numerous marine sources, including algae [39]. The SMP of strain 2 cultivated in five static and shaken media was overall poor. The crude extracts of CA-2 and CA-2-sh cultures contained only two peaks, which were identified as agistatin D (3) and orsellide D (4) from their UVmax and distinct masses (Table 2). This is the first report indicating the presence of agistatin D and orsellide D in the genus Cladosporium. Agistatin D is a pyranacetal originally isolated from a Fusarium sp. that inhibits the cholesterol biosynthesis [21]. Orsellide D is the deoxyhexose glycoside of orsellinic acid (2,4-dihydroxy-6-methylbenzoic acid) that was extracted from an alga-associated Chaetomium sp. It shows moderate antibacterial activity against Escherichia coli, Bacillus subtilis, and Staphylococcus aureus [22]. The chemical profile analyses showed that variations in culture conditions did not much affect the SM production of this strain, although there was a preference for CA medium. The crude extracts of both CA and CA-sh cultures containing agistatin D and orsellide D did not inhibit any of the test microorganisms but inhibited QS with variable potency (95% and 30% inhibition of luminescence in A. fischeri, respectively, see Table 3). It might be of interest to optimize the static culture conditions of this strain to induce the biosynthesis of other compounds and test their QSI potential.
Members of the Trichoderma genus (strain 3) have previously been isolated from Z. marina specimens [13], as well as from the Mediterranean seagrass Posidonia oceanica [40,41]. Trichoderma spp. are common and reported from various marine sources, including algae [42] and terrestrial plant rhizosphere [43]. Strain 3 showed strong medium preferences for static and shaken CA and WH media, while no compounds were detected in AS and GPY media extracts. From 10 peaks detected in different cultures of strain 3, only five were putatively dereplicated. Peak 5 was identified as orsellide C with reported antibacterial activity [40]. This compound was only detectable in the static WH culture. UP5 and UP6 might be potential new derivatives of orsellide C due to shared UV absorption. The compound 6, common to both static WH and CA cultures, was annotated as cyclotetrapeptide (Cyl-2) that originates from the fungal phytopathogen Cylindrocladium scoparium [23] and inhibits the growth of lettuce and rice seedlings [44]. Peak 7 was identified as quinolone B, an alkaloid originally isolated from a soil-derived Penicillium sp. with insecticidal and nematocidal activity [24,45]. Peak 8 was assigned to harzianopyridone (8), a small pyridine derivative with antimicrobial activity [46], including phytopathogens [25]. Compounds 7 and 8 were solely detected in static CA cultures. Peak 9 was identified as 2’-deoxyribofuranosylguanine (2N-Me), a primary metabolite expressed in PC, CA-sh, WH, and WH-sh cultures. All unknown compounds were solely detected in static WH cultures, while UP9 was observed only in the WH-sh extract. The extracts of strain 3 lacked antibacterial activity, except for the static PC and the WH extracts that inhibited A. fischeri. Almost all extracts, particularly those deriving from CA, CA-sh, WH, and WH-sh cultures, strongly inhibited QS (Table 3).
With ten isolated strains, the fungal genus Penicillium showed the highest abundance for Z. marina (Table 1). The strains (4–13) belonged to three species, P. atramentosum (strain 8), P. antarcticum (strains 4–7), and P. atrovenetum (strains 9–13). In general, Penicillium spp. are widely distributed fungi found in all moderate environments [47]. They are also the chemically best-studied genus from which various types of SMs, including polyketides, terpenoids, peptides, and alkaloids, have been reported [48]. Still, they continue providing new SMs [49,50], often due to sampling campaigns conducted at extreme environments, e.g., deep-sea [51]. Afiyatullov and colleagues [14] reported the isolation of a new phthalide derivative from Penicillium claviforme, a surface associated fungus of Z. marina sampled in the Sea of Japan. The same authors [15] isolated four new eudesmane-type sesquiterpenes from Penicillium thomii, from Z. marina sampled from the same location. Three of the four tested compounds induced a significant down-regulation of nitric oxide production in murine macrophages. To our knowledge, none of the previously mentioned studies conducted a systematic approach to investigate the chemical profiles and bioactivities of cultivable mycobiota of Z. marina. In this study, a total number of 69 compounds with different UVmax and molecular masses were detected from all 10 Penicillium spp. Of these, 46 (approx. 70%) were successfully dereplicated, while 23 peaks (approx. 30%) remained unidentified. Known compounds belonged to chemical classes of alkaloids, polyketides, and terpenoids (Supplementary Table S1). Interestingly, this annotation ratio applies largely to every single strain regardless of aeration and independent of taxonomy. As mentioned above, strain 8 appeared to be very unique in chemistry (Figure 5 and Supplementary Table S1). All twelve compounds, identified from this strain, were species-specific, with the only exception being the terpenoid molecule andrastin A that also occurs in strain 4 (P. antarcticum). Several other dereplicated compounds, as well as potential new (unknown) compounds (Supplementary Table S1) also existed in single strain or a single species. For example, atrovenetin (43) was specific to P. atrovenetum and present in all P. atrovenetum strains 9–13, but no other Penicillium species contained this compound. However, the majority of the compounds existed in multiple Penicillium samples, excluding a general strain or species-specificity. For example, all Penicillii, except for strain 8, produced compound 10 in static and shaken PC and AS media. Compound 10 was dereplicated as patulin, a mycotoxin that has been isolated from a variety of Penicillium spp. [52] with a broad range of bioactivities (see Section 3.4). Dictyosphaeric acid B (36), a polyketide, was common to all P. antarcticum and P. atroventumum strains (Supplementary Table S1). Another mycotoxin, dereplicated in almost all extracts of the Penicillium strains, was compound 54. It was assigned to aflatoxin M1, another mycotoxin that shows cytotoxicity and genotoxicity (e.g., DNA damage and gene mutation) in mammalian cells, insects, lower eukaryotes, and bacteria [53]. Besides mycotoxins, the Z. marina associated Penicilli appeared to produce compounds with a broad variety of activities. Examples include the lactone cytosporone E (18) from a Cytospora sp. [54], the terpenoid hongoquercin A (24) from an unidentified fungus [55], and the alkaloid bionectin B (37) that has been reported from the soil or plant-derived fungus Bionectra byssicola [56]. Reportedly, all these compounds show antibacterial activities. Furthermore, podosporin A (35), a quinone from the ascomycete Podospora decipiens [57], and talaroconvolutin B (55), a tetramic acid derivative originating from another fungus Talaromyces convolutus [58], are known to exhibit antifungal activity. Compound 38 (tubingensin A or B from Aspergillus tubingensis) has been shown to display antiviral activity against herpes simplex virus type 1 [59,60]. Compound 43 was dereplicated as atrovenetin, a polyketide isolated from Penicillium paraherquei with antioxidant activity [61]. Adding to this list of bioactivities, the terpenoids territrem C (33) from Aspergillus terreus [62], penitrem E (49) from P. crustosum [63], and compound 53 (penitrem A or pennigritrem from Penicillium sp.) [64] are mycotoxins with tremorgenic activity in mice [63]. Some compounds that are unique to the chemically distinct strain 8 (P. atramentosum) are known for their cytotoxicity and anticancer activity. Examples include globosumone B (12), an ester of orsellinic acid isolated from the fungus Chaetomium globosum with moderate activity against multiple tumor cell lines [65], or compound 51 (penicillenol A1 or A2), which is an anticancer tetramic acid derivative originating from a Penicillium sp. [66]. With their diversity and broad-spectrum bioactivities, Penicillium spp. that are abundantly found in Z. marina rhizosphere may be beneficial for the protection of the eelgrass roots.

3.3. The Influence of Media and Culture Regime on the Secondary Metabolome of Z. marina Associated Fungi

The OSMAC approach has proven to be a very successful strategy for the discovery of novel SMs. In this study, the selection of culture media significantly influenced the chemical diversity of Z. marina associated fungi. Figure 1, Figure 3, Figure 4, and Figure 5D indicate that WH medium provided chemically diverse SMPs in all associated fungi. Dereplication of all fungal extracts annotated metabolites unique to every single medium (independent of culture regime), i.e., WH with 26 unique SMs, AS with 12 unique SMs, CA with 11 unique SMs, PC with 10 unique SMs, and GPY with two unique SMs (Supplementary Table S1). The medium-dependent occurrence of SMs may suggest the activation of silent biosynthetic gene clusters by variation of culture conditions. Figure 1, Figure 3, Figure 4, and Figure 5D indicate that WH medium induced chemically diverse SMPs in all associated fungi. WH medium is rich in nutrients, containing high amounts of carbon sources (glucose and malt) and nitrogen sources (peptone and yeast), as well as sea salt mixture as a salt source. Penicillium spp. generally respond well to changes in growth medium composition, especially to salinity and salt sources [67,68]. The influence of the C source on the SM production has also been reported. Rodríguez-Ortiz and coworkers (2010) observed sucrose to stimulate the production of the antibiotic metabolite bikaverin in Fusarium fujikuroi [69]. Another study by Fan and colleagues (2019) showed that a specific family of linear aminolipids was exclusively detected in the sucrose-containing medium [70]. Besides carbon, the choice of the nitrogen source also plays an important role. An early study on Aspergillus flavus reported that the absence of yeast extract in the respective growth medium led to the silencing of mycotoxin production [71]. However, with 18 produced compounds, strain 9 (P. atrovenetum) was the best producer in static WH medium. The dendrogram of WH medium (Figure 5D) further emphasized the suitability of this medium for high SM production in Penicillium strains. In contrast to PC, AS, and GPY media (Figure 5A,E), the dendrogram of WH medium was highly branched, indicating a high chemical diversity. The AS medium is poor in nutrients, containing low amounts of glucose as a carbon source without additional N source. Therefore, the fungal growth was generally weak, and the AS medium extracts were chemically poor. Nevertheless, cultivation of the Penicillium strains in AS medium resulted in chemically diverse extracts (see Figure 5C). Strain 6 produced 11 SMs with only five known compounds and was the best producer in this medium under shaking conditions. The high number of unknown compounds in shaken AS cultures (see Figure 4) suggests that fungal strains can produce many putatively new natural products in this medium. The CA medium is rich in nitrogen sources (hydrolyzed casein and yeast extract) and supplemented with a sea salt mixture as a salt source and dipotassium phosphate. The PC medium contains a complex mixture of C and N sources without extra salt source. The highest TPN (seven compounds) in PC medium was achieved by strain 12 (P. atrovenetum) under shaking conditions. GPY medium contains glucose as a carbon source and peptone and yeast as a nitrogen source, all in low amounts. The highest TPN (5) in GPY medium was achieved by strain 8 under shaking conditions (Figure 4). In contrast to the three media mentioned above, PC and GPY media did not contain the sea salt mixture. These findings imply the necessity of salt for the cultivation of diverse SMs in Z. marina associated fungi.
The differences in chemical production under different culture conditions might also stem from changes in physical parameters. In general, shaking facilitates aeration and the uptake of oxygen [72]. This is necessary to provide adequate mixing and mass transfer, especially when fungal cells grow in a freely dispersed form, which results in a non-Newtonian broth with high viscosity [73]. In the present study, we applied shaking at 120 rpm to liquid cultures. Bigelis and coworkers (2006) researched fungal antibiotics and showed that fungi formed dispersed mycelia under shaking conditions but mycelial mat under static conditions [74]. Chemical production was significantly affected by different mycelial growth conditions. In the current study, we observed shaking led to the production of more (unknown) compounds in the AS medium (Figure 4). The average number of SMs strongly increased for all Penicillii cultured under shaking regime in AS medium. This positive effect of shaking was not observed for the other four media. For example, cultivation in shaken CA medium resulted in significantly lower numbers of detected SMs in all Penicillium strains compared to static CA medium. The average number of compounds of all Penicillium strains cultured under shaking conditions in CA media, therefore, strongly sank. In any case, the results indicate that the availability of higher levels of oxygen might be a critical factor for the production of SMs in different media. Guo and co-workers (2013) have previously reported the expression of specific fungal SM in a Penicillium sp., where the shift from static to shaking conditions led to the discovery of five new sorbicillinoids [75]. Noteworthy is that strain 8 generally exhibited completely different SMPs than the other Penicillii strains in both static and shaking regimes (Figure 5A–E).

3.4. Bioactivities of the Extracts of Z. marina Associated Fungi

All fungal extracts were tested for bioactivity against V. anguillarum and Ps. Elyakovii, as well as against QS in A. fischeri and biofilm formation of Ps. aeruginosa. V. anguillarum is a Gram-negative bacterium and causative agent of vibriosis, which is pathogenic to a variety of aquatic organisms, such as fish and crustaceans, with importance in aquaculture industry [76,77,78,79]. Ps. elyakovii is another Gram-negative bacterium that is isolated from spot-wounded fronds of the brown alga Laminaria japonica and is believed to cause this disease [80]. Screening for extracts containing active compounds for purification of chemicals with application potential in aquaculture is, therefore, of interest. Multi-drug resistant pathogenic bacteria are on the rise worldwide, and it is of great urgency to develop novel therapeutic measures. A promising approach is to manipulate the cell-to-cell communication (known as QS), which is crucial for microfouling of biological surfaces that increase the vulnerability of marine organisms to water-borne pathogens. QS inhibition is also being regarded as a starting point for the discovery and development of new antibiotics and antipathogenic drugs [81]. In this study, the marine Gram-negative bacterium A. fischeri, emitting bioluminescence through QS [82], was used for anti-quorum sensing assays. Because QS is directly linked to bacterial population-density, tests for growth inhibition in A. fischeri were performed as well. This bacterium can grow without quorum sensing, with the consequence that all QS-coupled properties (e.g., biofilm formation, virulence factors, etc.) will be inactivated causing relevant changes in the ecological functions of the bacterial population. Accordingly, the interference with the QS system by natural products is considered as very relevant for ecological interaction [83]. The Gram-negative bacterium Ps. aeruginosa is an opportunistic human pathogen that exhibits various antibiotic resistance mechanisms and virulence factors, such as biofilm formation. Given the limited antimicrobial arsenal against the resistant strains of this bacterium, the discovery of new active compounds against Ps. aeruginosa is a priority [84]. Pseudomonas species are also widely distributed in marine systems. Hence, activity in this assay can be an indication for activity relevant for bacterial biofilm fouling of marine surfaces.
In general, nearly all extracts of static and shaken PC and AS cultures of Z. marina associated fungi showed potent QSI activity with poor inhibition against A. fischeri. The only exceptions were the AS and AS-sh extracts of strain 8 and the AS extract of strain 2. In contrast, many static and shaken CA, WH, and GPY extracts poorly inhibited QS and the growth of A. fischeri. It is noteworthy that all static and shaken PC and AS crude extracts (except for strain 8) contained the metabolite 10 (patulin). Besides being a mycotoxin, patulin, isolated from a P. coprobium, was found to inhibit quorum sensing in Ps. aeruginosa [85]. Therefore, this compound might have had some impact on the QSI activity. However, as the crude extracts of the non-Penicillium strains lacked compound 10, some other QSI metabolites might also be present. Extracts from static and shaken PC and AS cultures of almost all Penicillium strains strongly inhibited the growth of V. anguillarum and Ps. elyakovii, except those from strain 8. Patulin (10) has been also reported to exhibit antimicrobial activity [86,87], which might explain the activity of the crude extracts containing 10. In contrast to the Penicillium strains, none of the non-Penicillium strains inhibited the growth of V. anguillarum and Ps. elyakovii. Some Penicillium PC and AS culture extracts inhibited the growth of Ps. aeruginosa, although no specific biofilm formation inhibition was recorded. In contrast, the non-Penicillium strains failed to show any significant inhibition of biofilm formation or the growth of Ps. aeruginosa. In general, the bioassays showed that many Z. marina associated fungi (especially those grown in PC and AS media) were active against a variety of Gram-negative bacteria. Furthermore, compounds like patulin (10) and aflatoxin M1 (54) are known for being potent mycotoxins. Due to reduced chemical machinery and defense mechanisms of Z. marina during its re-adaptation from land to sea [4,8], these fungal SMs might have played a role in the chemical protection of this plant through inhibition against growth, QS, and formation of biofilm by marine pathogens. Furthermore, the observed bioactivities might help to generate new drugs against Gram-negative, clinically, and environmentally relevant pathogens in the future.
In conclusion, we have isolated and identified, for the first time by a systematical approach, the cultivable fungal community associated with the Baltic Z. marina. The use of five different media and two culture regimes, as well as the application of a rapid dereplication method using HPLC-DAD-MS, enabled the identification of optimal culture conditions, revealing the chemical diversity and bioactivity of the fungal extracts. Because of their diverse SMPs with putatively new compounds and antimicrobial and/or QS inhibitory activities against Gram-negative bacteria, several strains, e.g., 1, 6, 8, and 11, will be submitted to large-scale cultivation followed by isolation and structure elucidation of their bioactive components.

4. Materials and Methods

4.1. Collection of Zostera marina

Fifteen living specimens of Z. marina were collected by SCUBA diving in a water depth of approx. 3 m in Friedrichsort (54°23′N; 10°11′E) around Kiel Fjord (Baltic Sea, Germany). The eelgrass leaf specimens were stored in sterile plastic bags at 4 °C before being transported and processed in the laboratory. For sampling the sediment surrounding the roots of Z. marina plants, 50 mL sterile Falcon tubes were used for coring the uppermost sediment surrounding the roots (five sediment samples collected). All samples were placed into a cool box, transported to the laboratory, and processed immediately.

4.2. Isolation of Fungi

The Zostera leaf specimens were cut (0.5 cm) and placed into a sterile innuSPEED Lysis Tube (Analytik Jena AG, Jena, Germany) together with 100 µL of sterile saline and homogenized for 6 min at a frequency of 30 s−1 in a RetschMM200 mixer mill (Retsch GmbH, Haan, Germany). Homogenates of Z. marina leaf subsamples and sediment samples were diluted 1:10 and 1:100 with sterile water for a dilution series. The 100 µL of each dilution step (10−1 and 10−2), 15 µL of undiluted homogenate, and a small spoon of sediment surrounding the roots were plated onto five Agar media; four Glucose-Yeast-Malt 4 medium (GYM4, glucose monohydrate 4 g, yeast extract 4 g, malt extract 4 g, CaCO3 2 g, aqua dest. ad 1 L, pH 7.2), Glucose-Peptone-Yeast medium (GPY, glucose monohydrate 1 g, peptone 0.5 g, yeast extract 0.1 g, NaCl 15 g, agar 13 g, aqua dest. ad 1 L, pH 7.3), Marine broth medium (MB, Bacto, agar 15 g, aqua dest. ad 1 L), Tryptic Soy Broth 12 medium (TSB12, 12 g tryptic soy broth, 5 g NaCl, agar 15 g, aqua dest. ad 1 L, pH 7.2), and Hastings medium (HS, Na2HPO4x 12H2O 9.35 g, KH2PO4 1 g, (NH4)2SO4 0.5 g, MgSO4 × 7H2O 0.21 g, NaCl 30 g, trypton 5 g, yeast extract 3 g, glycerol 2 mL, agar 18 g, aqua dest ad 1 L, pH 7.4), for the cultivation of associated microbes. The inoculated Petri dishes were stored at 22 °C in the dark. All plates were checked for different morphotypes. Different morphotypes were isolated by streaking on new agar plates of the corresponding media to obtain axenic cultures, whenever they occurred in a total time frame of 21 days.

4.3. Identification of Fungal Strains

Fungal morphology and ITS sequences were used for taxonomic identification of Z. marina associated fungi. Fungal strains were studied by macroscopic (colony morphology) and microscopic means (hyphae, conidia, conidiophors, septation). Molecular identification was based on the sequence of the ITS1 nrDNA region. Fungal genomic DNA was obtained by simple cell lysis in nuclease-free water with a Retsch mill model MM200 (Retsch GmbH, Haan, Germany). 1 µL of this extract was used as a template for PCR amplifying ITS1-5.8S-ITS2 region with the primers [88] in a total reaction volume of 25 µL (12.5 µL Dream Taq Master Mix (ThermoFisher Scientific, Waltham, MA, USA), 1 µL of each primer (concentration 10 μM, ITS1 [5′-TCC GTA GGT GAA CCT GCG G-3′] and ITS4 [5′-TCC TCC GCT TAT TGA TAT GC-3′] [88]), and 9.5 µL of DNA free water (ThermoFisher Scientific, Waltham, MA, USA). PCR was performed in a T-1 thermocycler (Biometra, Göttingen, Germany) using the following conditions for amplification: denaturation 94 °C for 8 min, 35 cycles of denaturation 94 °C for 30 s, annealing 48 °C for 45 s, elongation 72 °C for 3 min, final elongation 72 °C for 10 min, cooling 4 °C [88]. PCR products were checked for correct length by gel electrophoresis using an SYBR safe (ThermoFisher Scientific, Waltham, MA, USA) pre-stained 1% TBE Agarose gel and Gene ruler 1 Kb DNA ladder (Thermo ScientificGeneRuler™, Schwerte, Germany) as length control. Electrophoresis conditions were 120 V for 20 min. Positive PCR products were analyzed by Sanger sequencing at the Institute of Clinical Molecular Biology (Kiel University) using the primer ITS4. All sequences were processed via the software ChromasPro V1.33 (technelysium Pty Ltd., Brisbane, Australia) and compared to the NCBI database (https://www.ncbi.nlm.nih.gov/genbank/) using the BLAST algorithm.

4.4. Cultivation and Extraction of Marine Fungi

Pre-cultures of fungi were cultivated on the same medium used for isolation of the respective strain for two weeks at 22 °C in the dark. Approximately 1 cm2 piece of overgrown agar was used to inoculate five liquid media (100 mL), including Potato-Carrot (PC, carrot 20 g, potato 20 g, aqua dest. ad 1 L, pH 7.2), Casamino Acids (CA, casein hydrolysate 30 g, yeast extract 4 g, K2HPO4 0.5 g, sea salt mixture (Instant Ocean) 15 g, aqua dest. ad 1 L, pH 7.4), Wickerham15 (WH, glucose monohydrate 10 g, peptone from soymeal 5 g, yeast extract 3 g, malt extract 3 g, sea salt mixture (Instant Ocean) 15 g, aqua dest. ad 1L, pH 7.3), Glucose-Peptone-Yeast (GPY, glucose monohydrate 1 g, peptone 0.5 g, yeast extract 0.1 g, NaCl 15 g, aqua dest. ad 1 L, pH 7.3), Artificial Seawater (AS, glucose monohydrate 1 g, sea salt mixture (Instant Ocean) 15 g, aqua dest. ad 1 L, pH 7.5), for three weeks at 22 °C in the dark either in a shaking incubator (Edmund Bühler GmbH, Bodelshausen, Germany) with 120 rpm (abbreviated as “sh” cultures) or without shaking. The complete culture broth was homogenized at 19,000 rpm for 30 s using a T25 Ultra Turrax (IKA Werke GmbH & Co. KG, Staufen im Breisgau, Germany). The suspension was extracted with ethyl acetate (EtOAc) and centrifuged (4700 rpm) for 10 min at 4 °C. The aqueous phases were discharged, while the respective organic phase was evaporated to dryness at 40 °C by a rotary evaporator. The dried extracts were stored at –20 °C until work-up.

4.5. LC-DAD-MS Analyses and Dereplication

The crude extracts were analyzed by an HPLC (LaChrom Elite, VWR-Hitachi, VWR International, Radnor, PA, USA) system hyphenated to a Diode Array Detector L-2450 (DAD) and an Esquire 4000 Bruker Daltonics mass spectrometer (ESIMS). The extracts were dissolved in HPLC grade MeOH (1 mg/mL), filtered through a 0.2 µm PTFE syringe filter, and injected (20 µL) into the HPLC-DAD-MS system. Chromatographic separations were achieved at 40 °C using a monolithic C18 column (Onyx monolithic C18, 100 × 3.0 mm, Phenomenex, Torrance, CA, USA). The mobile phase used was A: 95.0% MilliQ-water/0.1% formic acid (Fluka) and B: ACN (Actu-All Chemicals BV)/0.1% formic acid, pumped at a rate of 2.0 mL/min with the following gradient: initial, 95% A; 0–4 min 40% A; 6–6.9 min 0% A; 11.5–12.5 min 0% A, and a column reconditioning phase until 8.5 min at a flow of 2.5 mL/min. If not stated otherwise, measurements of crude extracts were conducted at positive ionization mode. As a result, monomer ions ([M + H]+, [M + Na]+) were primarily formed with dimer ions ([2M + H])+, [2M + Na]+) or potassium adducts ([M + K]+) being less commonly found. Both UV-VIS (200–500 nm) and MS data were evaluated using the programs Data Analysis and HyStar Post Processing (Bruker Daltonics, Bruker Corporation, MA, USA). For identification of known compounds, datasets of each crude extract were compared with entries in the database “Dictionary of Natural Products” as well as other publicly available databanks. Data were reconciled with database entries in matters of UV maxima, accurate mass, and biological source.

4.6. Bioassays

All bioassays were performed in a 96-well microtiter format. The crude extracts were dissolved in DMSO (10 mg/mL stock solution) and diluted with the respective media to a final concentration of 0.5% DSMO in test volumes of 200 µL and extract concentrations of 100 µg/mL. Antibiotics (see below) served as positive controls, whereas inoculated medium (containing DMSO but no extracts), blanks (just medium), and DMSO served as negative controls. Assays were prepared in duplicates.

4.6.1. Antibacterial Activity

Crude extracts were screened for their antibacterial activity against Ps. elyakovii and V. anguillarum using resazurin as an indicator. Overnight cultures (28 °C, 160 rpm, Tropic Marine (TM, peptone 5.00 g, yeast extract 1.00 g, sea salt mixture (Instant Ocean) 30.00 g, aqua dest. ad 1 L, pH 7.6) of Ps. elyakovii were diluted to an optical density (OD) of 0.03 with medium, whereas those of V. anguillarum were diluted to an OD of 0.01 and inoculated into the test wells. Chloramphenicol was used as a positive control. Microplates were incubated at 28 °C and 200 rpm for 7 h (Ps. elyakovii) and 5 h (V. anguillarum). Subsequently, 10 µL of resazurin (0.1 mg/mL) solution was added to each well. Microplates were incubated again at RT and 200 rpm in the dark for 20 min. (Ps. elyakovii) and for 15 min. (V. anguillarum), and the fluorescence was measured at 590 nm using an excitation wavelength of 560 nm.

4.6.2. Quorum Sensing Inhibition

Pre-cultures (overnight at 28 °C and 160 rpm in HS medium) of A. fischeri were diluted to an OD of 0.01. Microplates were prepared as above, using 5-(bromomethylene)-2H(5H)-furanone and chloramphenicol as positive controls. Initial cell density was measured at 600 nm. The content of the microplate was transferred into a white microplate suitable for luminescence measurement and incubated at 28 °C and 200 rpm for 5 h. Subsequently, luminescence was measured in an Infinite M200 reader (TECAN Deutschland GmbH, Crailsheim, Germany) for 1000 ms. At a signal intensity of 1 × 106 RLU (relative light unit) in the negative control, cell density was measured at 600 nm in a transparent microplate.

4.6.3. Biofilm Inhibition

TSB12 medium (40 mL) was inoculated with Ps. aeruginosa and incubated for 5 h at 37 °C and 200 rpm. This culture was diluted to an OD of 0.01 and incubated overnight at 37 °C and 200 rpm. For testing biofilm formation, this culture was dispersed into microtiter plates at a volume of 150 µL/well. Polymyxin B served as a positive control. A sterile PCR plate was submerged into the microplates, serving a surface for the biofilm. The constructs consisting of microplates and PCR plates were tightly fixed and incubated for 24 h at 37 °C. After incubation, PCR plates were washed three times with 0.9% sterile saline solution and subsequently dried at 60 °C for 1 h for fixation of biofilms. In the next step, each dried PCR plate was incubated in a 0.1% crystal violet solution at RT for 20 min. All plates were washed again three times with 0.9% sterile saline solution to remove unbound crystal violet. Finally, each PCR plate was placed into a new microplate prefilled with 200 µL 96% EtOH, and both were shaken at 22 °C and 100 rpm for 20 min. In the end, PCR plates were removed, and each microplate was measured for absorbance at 590 nm. As all crystal violet was originally bound to the biofilm, the amount of crystal violet estimated from the absorbance could be related to the amount of formed biofilm.

4.7. Statistics

All statistics were performed, using the soft wares Microsoft Office Excel 2016 and Past 3.10. Multivariate data analyses in the form of UPGMA (Unweighted Pair Group Method with Arithmetic mean) clustering were used for the determination of relatedness between different SMs of all (Penicillium) isolates per medium. UPGMA clustering was based on the calculation of Bray-Curtis dissimilarity. UPGMA clustering was based on the presence-absence matrices of the respective SMs. Only peaks with absorbance and mass information were counted as “present”. If this was not the case, the compound was declared as “absent”. For the construction of matrices used for UPGMA clusters, presence-absence data were normalized using the Hellinger transformation. This transformation was done for giving less weight to rare SMs in the transformed presence-absence table [89].

Supplementary Materials

The following are available online at https://www.mdpi.com/1660-3397/17/7/419/s1, Figures S1–S13: Metabolite profile of strains 1–13 cultured under static and shaking conditions in five different media; Figure S14: Chemical structures of dereplicated compounds in Zostera marina associated fungi; Table S1: Spectral data used for dereplication of all Penicillium spp. extracts.

Author Contributions

Design of the project, D.T. and A.L.; Initial isolation of fungi, M.M.; Growth, chemical analyses, and bioactivity assays, L.E.P.; Write-up of the paper, D.T, L.E.P., and A.L. All authors have read and approved the final version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

We thank Arlette Wenzel-Storjohann for performing the bioactivity assays.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lee, H.; Golicz, A.A.; Bayer, P.E.; Jiao, Y.; Tang, H.; Paterson, A.H.; Sablok, G.; Krishnaraj, R.R.; Chan, C.-K.K.; Batley, J.; et al. The genome of a southern hemisphere seagrass species (Zostera muelleri). Plant Physiol. 2016, 172, 272–283. [Google Scholar] [CrossRef] [PubMed]
  2. Bengtsson, M.M.; Bühler, A.; Brauer, A.; Dahlke, S.; Schubert, H.; Blindow, I. Eelgrass leaf surface microbiomes are locally variable and highly correlated with epibiotic eukaryotes. Front. Microbiol. 2017, 8, 1312. [Google Scholar] [CrossRef] [PubMed]
  3. Boström, C.; Baden, S.; Bockelmann, A.-C.; Dromph, K.; Fredriksen, S.; Gustafsson, C.; Krause-Jensen, D.; Möller, T.; Nielsen, S.L.; Olesen, B.; et al. Distribution, structure and function of Nordic eelgrass (Zostera marina) ecosystems: Implications for coastal management and conservation. Aquat. Conserv. Mar. Freshw. Ecosyst. 2014, 24, 410–434. [Google Scholar] [CrossRef] [PubMed]
  4. Olsen, J.L.; Rouzé, P.; Verhelst, B.; Lin, Y.-C.; Bayer, T.; Collen, J.; Dattolo, E.; De Paoli, E.; Dittami, S.; Maumus, F.; et al. The genome of the seagrass Zostera marina reveals angiosperm adaptation to the sea. Nature 2016, 530, 331–335. [Google Scholar] [CrossRef] [PubMed]
  5. Fahimipour, A.K.; Kardish, M.R.; Lang, J.M.; Green, J.L.; Eisen, J.A.; Stachowicz, J.J. Global-scale structure of the eelgrass microbiome. Appl. Environ. Microbiol. 2017, 83, e03391-16. [Google Scholar] [CrossRef]
  6. Guan, C.; Parrot, D.; Wiese, J.; Sönnichsen, F.D.; Saha, M.; Tasdemir, D.; Weinberger, F. Identification of rosmarinic acid and sulfated flavonoids as inhibitors of microfouling on the surface of eelgrass Zostera marina. Biofouling 2017, 33, 867–880. [Google Scholar] [CrossRef] [PubMed]
  7. Inaba, N.; Trainer, V.L.; Onishi, Y.; Ishii, K.I.; Wyllie-Echeverria, S.; Imai, I. Algicidal and growth-inhibiting bacteria associated with seagrass and macroalgae beds in Puget Sound, WA, USA. Harmful Algae 2017, 62, 136–147. [Google Scholar] [CrossRef] [PubMed]
  8. Papazian, S.; Parrot, D.; Burýšková, B.; Weinberger, F.; Tasdemir, D. Surface chemical defence of the eelgrass Zostera marina against microbial foulers. Sci. Rep. 2019, 9, 3323. [Google Scholar] [CrossRef]
  9. Vandenkoornhuyse, P.; Quaiser, A.; Duhamel, M.; Le Van, A.; Dufresne, A. The importance of the microbiome of the plant holobiont. New Phytol. 2015, 206, 1196–1206. [Google Scholar] [CrossRef]
  10. Ugarelli, K.; Chakrabarti, S.; Laas, P.; Stingl, U. The seagrass holobiont and its microbiome. Microorganisms 2017, 5, 81. [Google Scholar] [CrossRef]
  11. Ettinger, C.L.; Voerman, S.E.; Lang, J.M.; Stachowicz, J.J.; Eisen, J.A. Microbial communities in sediment from Zostera marina patches, but not the Z. marina leaf or root microbiomes, vary in relation to distance from patch edge. PeerJ 2017, 5, e3246. [Google Scholar] [CrossRef] [PubMed]
  12. Ettinger, C.L.; Williams, S.L.; Abbott, J.M.; Stachowicz, J.J.; Eisen, J.A. Microbiome succession during ammonification in eelgrass bed sediments. PeerJ 2017, 5, e3674. [Google Scholar] [CrossRef]
  13. Kirichuk, N.N.; Pivkin, M.V. Filamentous fungi associated with the seagrass Zostera marina Linnaeus, 1753 of Rifovaya Bay (Peter the Great Bay, the Sea of Japan). Russ. J. Mar. Biol. 2015, 41, 351–355. [Google Scholar] [CrossRef]
  14. Afiyatullov, S.S.; Leshchenko, E.V.; Sobolevskaya, M.P.; Gerasimenko, A.V.; Khudyakova, Y.V.; Kirichuk, N.N.; Mikhailov, V.V. New 3-[2′(R)-hydroxybutyl]-7-hydroxyphtalide from marine isolate of the fungus Penicillium claviforme. Chem. Nat. Compd. 2015, 51, 111–115. [Google Scholar] [CrossRef]
  15. Afiyatullov, S.S.; Leshchenko, E.V.; Sobolevskaya, M.P.; Denisenko, V.A.; Kirichuk, N.N.; Khudyakova, Y.V.; Hoai, T.P.T.; Dmitrenok, P.S.; Menchinskaya, E.S.; Pislyagin, E.A.; et al. New eudesmane sesquiterpenes from the marine-derived fungus Penicillium thomii. Phytochem. Lett. 2015, 14, 209–214. [Google Scholar] [CrossRef]
  16. Hyde, K.D.; Sarma, V.V.; Jones, E.B.G. Morphology and taxonomy of higher marine fungi. In Marine Mycology—A practical Approach, 1st ed.; Hyde, K.D., Pointing, S.B., Eds.; Fungal Diversity Press: Hong Kong, China, 2000; pp. 172–204. [Google Scholar]
  17. Ruby, E.G. Lessons from a cooperative, bacterial-animal association: The Vibrio fischeriEuprymna scolopes light organ symbiosis. Annu. Rev. Microbiol. 1996, 50, 591–624. [Google Scholar] [CrossRef]
  18. Ruby, E.G.; Urbanowski, M.; Campbell, J.; Dunn, A.; Faini, M.; Gunsalus, R.; Lostroh, P.; Lupp, C.; McCann, J.; Millikan, D.; et al. Complete genome sequence of Vibrio fischeri: A symbiotic bacterium with pathogenic congeners. Proc. Natl. Acad. Sci. USA 2005, 102, 3004–3009. [Google Scholar] [CrossRef] [PubMed]
  19. Scheel, L.; Perone, V.B.; Larkin, R.L.; Kupel, R.E. The isolation and characterization of two phototoxic furanocoumarins (psoralens) from diseased celery. Biochemistry 1963, 2, 1127–1131. [Google Scholar] [CrossRef] [PubMed]
  20. Nippon Kayaku KK. Lipoxygenase Inhibitor, Antitumour NK-A-17-e-233-II—Produced by Cultivation of Phoma sp. NK-A-17-e-233. Japan Patent Office JP63216894-A, 1988. [Google Scholar]
  21. Zeeck, A. Agistatins, Metabolites from Fusarium sp., Process for Their Preparation and Their Use. Deutsches Patent- und Markenamt EP0492318A2, 13 December 1991. [Google Scholar]
  22. Schlörke, O.; Zeeck, A. Orsellides A–E: An example for 6-deoxyhexose derivatives produced by fungi. Eur. J. Org. Chem. 2006, 2006, 1043–1049. [Google Scholar] [CrossRef]
  23. Hirota, A.; Suzuki, A.; Aizawa, K.; Tamura, S. Structure of Cyl-2, a novel cyclotetrapeptide from Cylindrocladium scoparium. Agric. Biol. Chem. 1973, 37, 955–956. [Google Scholar] [CrossRef]
  24. Hayashi, H.; Nakatani, T.; Inoue, Y.; Nakayama, M.; Nozaki, H. New dihydroquinolinone toxic to Artemia salina produced by Penicillium sp. NTC-47. Biosci. Biotechnol. Biochem. 1997, 61, 914–916. [Google Scholar] [CrossRef]
  25. Dickinson, J.M.; Hanson, J.R.; Hitchcock, P.B.; Claydon, N. Structure and biosynthesis of harzianopyridone, an antifungal metabolite of Trichoderma harzianum. J. Chem. Soc. Perkin Trans. 1 1989, 11, 1885–1887. [Google Scholar] [CrossRef]
  26. Robins, M.J.; Robins, R.K. Purine nucleosides. XXIV. A new method for the synthesis guanine nucleosides. The preparation of 2′-deoxy-α and β-guanosines and the corresponding N2-methyl derivatives. J. Phys. Chem. 1969, 34, 2160–2163. [Google Scholar] [CrossRef]
  27. Newell, S.Y. Fungi and bacteria in or on leaves of eelgrass (Zostera marina L.) from Chesapeake Bay. Appl. Environ. Microbiol. 1981, 41, 1219–1224. [Google Scholar]
  28. Kossuga, M.H.; Romminger, S.; Xavier, C.; Milanetto, M.C.; do Valle, M.C.; Pimenta, E.F.; Morais, R.P.; de Carvalho, E.; Mizuno, C.M.; Coradello, L.F.C.; et al. Evaluating methods for the isolation of marine-derived fungal strains and production of bioactive secondary metabolites. Braz. J. Pharmacogn. 2012, 22, 257–267. [Google Scholar] [CrossRef] [Green Version]
  29. Cúcio, C.; Engelen, A.H.; Costa, R.; Muyzer, G. Rhizosphere microbiomes of European seagrasses are selected by the plant, but are not species specific. Front. Microbiol. 2016, 7, 440. [Google Scholar] [CrossRef]
  30. Pollard, P.C.; Moriarty, D.J.W. Organic carbon decomposition, primary and bacterial productivity, and sulphate reduction, in tropical seagrass beds of the Gulf of Carpentaria, Australia. Mar. Ecol. Prog. Ser. 1991, 69, 149–159. [Google Scholar] [CrossRef]
  31. Brodersen, K.E.; Nielsen, D.A.; Ralph, P.J.; Kühl, M. Oxic microshield and local pH enhancement protects Zostera muelleri from sediment derived hydrogen sulphide. New Phytol. 2015, 205, 1264–1276. [Google Scholar] [CrossRef]
  32. Brodersen, K.E.; Siboni, N.; Nielsen, D.A.; Pernice, M.; Ralph, P.J.; Seymour, J.; Kühl, M. Seagrass rhizosphere microenvironment alters plant-associated microbial community composition. Environ. Microbiol. 2018, 20, 2854–2864. [Google Scholar] [CrossRef]
  33. 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. [Google Scholar] [CrossRef]
  34. Rai, M.; Deshmukh, P.; Gade, A.; Ingle, A.; Kövics, G.J.; Irinyi, L. Phoma saccardo: Distribution, secondary metabolite production and biotechnological applications. Crit. Rev. Microbiol. 2009, 35, 182–196. [Google Scholar] [CrossRef]
  35. Bubnova, E.N. Fungal diversity in bottom sediments of the Kara Sea. Bot. Mar. 2010, 53, 595–600. [Google Scholar] [CrossRef]
  36. Zhang, X.-Y.; Tang, G.-L.; Xu, X.-Y.; Nong, X.-H.; Qi, S.-H. Insights into deep-sea sediment fungal communities from the East Indian Ocean using targeted environmental sequencing combined with traditional cultivation. PLoS ONE 2014, 9, e109118. [Google Scholar] [CrossRef]
  37. Garg, B.J.; Saraswat, A.; Bhatia, A.; Katare, O.P. Topical treatment in vitiligo and the potential uses of new drug delivery systems. Indian J. Dermatol. Venereol. Leprol. 2010, 76, 231–238. [Google Scholar] [CrossRef]
  38. Leite, V.C.; Santos, R.F.; Lee, C.C.; Guillo, L.A. Psoralen derivatives and longwave ultraviolet irradiation are active in vitro against human melanoma cell line. J. Photochem. Photobiol. 2004, 76, 49–53. [Google Scholar] [CrossRef]
  39. Venkatachalam, A.; Govinda Rajulu, M.B.; Thirunavukkarasu, N.; Suryanarayanan, T.S. Endophytic fungi of marine algae and seagrasses: A novel source of chitin modifying enzymes. Mycosphere 2015, 6, 345–355. [Google Scholar] [CrossRef]
  40. Panno, L.; Voyron, S.; Anastasi, A.; Mussat Sartor, R.; Varese, G.C. Biodiversity of marine fungi associated with the seagrass Posidonia oceanica: An ecological and biotechnological perspective. Biol. Mar. Mediterr. 2011, 18, 85–88. [Google Scholar]
  41. Panno, L.; Bruno, M.; Voyron, S.; Anastasi, A.; Gnavi, G.; Miserere, L.; Varese, G.C. Diversity, ecological role and potential biotechnological applications of marine fungi associated to the seagrass Posidonia oceanica. New Biotechnol. 2013, 30, 685–694. [Google Scholar] [CrossRef]
  42. Miao, F.P.; Liang, X.-R.; Yin, X.-L.; Wang, G.; Ji, N.-Y. Absolute configurations of unique harziane diterpenes from Trichoderma species. Org. Lett. 2012, 14, 3815–3817. [Google Scholar] [CrossRef]
  43. Hermosa, R.; Viterbo, A.; Chet, I.; Monte, E. Plant-beneficial effects of Trichoderma and of its genes. Microbiology 2012, 158, 17–25. [Google Scholar] [CrossRef]
  44. Hirota, A.; Suzuki, A.; Suzuki, H.; Tamura, S. Isolation and biological activity of cyl-2, a metabolite of Cylindrocladium scoparium. Agric. Biol. Chem. 1973, 37, 643–647. [Google Scholar] [CrossRef]
  45. Kusano, M.; Koshino, H.; Uzawa, J.; Fujioka, S.; Kawano, T.; Kimura, Y. Nematicidal alkaloids and related compounds produced by the fungus Penicillium cf. simplicissimum. Biosci. Biotechnol. Biochem. 2000, 64, 2559–2568. [Google Scholar] [CrossRef]
  46. Cutler, H.G.; Jacyno, J.M. Biological activity of (–)-harziano-pyridone isolated from Trichoderma harzianum. Agric. Biol. Chem. 1991, 55, 2629–2631. [Google Scholar] [CrossRef]
  47. Roslan, H.A.; Ngo, C.S.; Muid, S. Genetic diversity of Penicillium species isolated from various sources in Sarawak, Malaysia. J. Cell Mol. Biol. 2010, 7, 13–23. [Google Scholar]
  48. Rateb, M.E.; Ebel, R. Secondary metabolites of fungi from marine habitats. Nat. Prod. Rep. 2011, 28, 290–344. [Google Scholar] [CrossRef]
  49. Yurchenko, A.N.; Smetanina, O.F.; Ivanets, E.V.; Kalinovsky, A.I.; Khudyakova, Y.V.; Kirichuk, N.N.; Popov, R.S.; Bokemeyer, C.; von Amsberg, G.; Chingizova, E.A.; et al. Pretrichodermamides D–F from a marine algicolous fungus Penicillium sp. KMM 4672. Mar. Drugs 2016, 14, 122. [Google Scholar] [CrossRef]
  50. Bu, Y.-Y.; Yamazaki, H.; Takahashi, O.; Kirikoshi, R.; Ukai, K.; Namikoshi, M. Penicyrones A and B, an epimeric pair of α-pyrone-type polyketides produced by the marine-derived Penicillium sp. J. Antibiot. 2016, 69, 57–61. [Google Scholar] [CrossRef]
  51. Guo, W.; Zhang, Z.; Zhu, T.; Gu, Q.; Li, D. Penicyclones A–E, antibacterial polyketides from the deep-sea-derived fungus Penicillium sp. F23-2. J. Nat. Prod. 2015, 78, 2699–2703. [Google Scholar] [CrossRef]
  52. Puel, O.; Galtier, P.; Oswald, I.P. Biosynthesis and toxicological effects of patulin. Toxins 2010, 2, 613–631. [Google Scholar] [CrossRef]
  53. Prandini, A.; Tansini, G.; Sigolo, S.; Filippi, L.; Laporta, M.; Piva, G. On the occurrence of aflatoxin M1 in milk and dairy products. Food Chem. Toxicol. 2009, 47, 984–991. [Google Scholar] [CrossRef]
  54. Brady, S.F.; Wagenaar, M.M.; Singh, M.P.; Janso, J.E.; Clardy, J. The cytosporones, new octaketide antibiotics isolated from an endophytic fungus. Org. Lett. 2000, 2, 4043–4046. [Google Scholar] [CrossRef]
  55. Roll, D.M.; Manning, J.K.; Carter, G.T. Hongoquercins A and B, new sesquiterpenoid antibiotics: Isolation, structure elucidation, and antibacterial activity. J. Antibiot. 1998, 51, 635–639. [Google Scholar] [CrossRef]
  56. Zheng, C.-J.; Kim, C.-J.; Bae, K.S.; Kim, Y.-H.; Kim, W.-G. Bionectins A-C, epidithiodioxopiperazines with anti-MRSA activity, from Bionectra byssicola F120. J. Nat. Prod. 2006, 69, 1816–1819. [Google Scholar] [CrossRef]
  57. Weber, H.A.; Baenziger, N.C.; Gloer, J.B. Podosporin A: A novel antifungal metabolite from the coprophilous fungus Podospora decipiens (Wint.) Niessl. J. Org. Chem. 1988, 53, 4567–4569. [Google Scholar] [CrossRef]
  58. Suzuki, S.; Hosoe, T.; Nozawa, K.; Kawai, K.I.; Yaguchi, T.; Udagawa, S.I. Antifungal substances against pathogenic fungi, talaroconvolutins, from Talaromyces convolutes. J. Nat. Prod. 2000, 63, 768–772. [Google Scholar] [CrossRef]
  59. TePaske, M.R.; Gloer, J.B.; Wicklow, D.T.; Dowd, P.F. Tubingensin A: An antiviral carbazole alkaloid from the sclerotia of Aspergillus tubingensis. J. Org. Chem. 1989, 54, 4743–4746. [Google Scholar] [CrossRef]
  60. TePaske, M.R.; Gloer, J.B. The structure of tubingensin B: A cytotoxic carbazole alkaloid from the sclerotia of Aspergillus tubingensis. Tetrahedron Lett. 1989, 30, 5965–5968. [Google Scholar] [CrossRef]
  61. Ishikawa, Y.; Morimoto, K.; Iseki, S. Atrovenetin as a potent antioxidant compound from Penicillium species. J. Am. Oil Chem. Soc. 1991, 68, 666–668. [Google Scholar] [CrossRef]
  62. Ling, K.H.; Liou, H.-H.; Yang, C.-M.; Yang, C.-K. Isolation, chemical structure, acute toxicity, and some physicochemical properties of territrem C from Aspergillus terreus. Appl. Environ. Microbiol. 1984, 47, 98–100. [Google Scholar]
  63. Kyriakidis, N.; Waight, E.S.; Day, J.B.; Mantle, P.G. Novel metabolites from Penicillium crustosum, including penitrem E, a tremorgenic mycotoxin. Appl. Environ. Microbiol. 1981, 42, 61–62. [Google Scholar]
  64. Penn, J.; Biddle, J.R.; Mantle, P.G.; Bilton, J.N.; Sheppard, R.N. Pennigritrem, a naturally-occurring penitrem A analogue with novel cyclisation in the diterpenoid moiety. J. Chem. Soc. Perkin Trans. 1 1992, 1, 23–26. [Google Scholar] [CrossRef]
  65. Bashyal, B.P.; Kithsiri Wijeratne, E.M.; Faeth, S.H.; Leslie Gunatilaka, A.A. Globosumones A-C, cytotoxic orsellinic acid esters from the Sonoran Desert endophytic fungus Chaetomium globosum. J. Nat. Prod. 2005, 68, 724–728. [Google Scholar] [CrossRef]
  66. Lin, Z.-J.; Lu, Z.-Y.; Zhu, T.-J.; Fang, Y.-C.; Gu, Q.-Q.; Zhu, W.-M. Penicillenols from Penicillium sp. GQ-7, an endophytic fungus associated with Aegiceras corniculatum. Chem. Pharm. Bull. 2008, 56, 217–221. [Google Scholar] [CrossRef]
  67. Fill, T.P.; Pallinib, H.F.; da Silva Amaral, L.; da Silva, J.V.; Bidóia, D.L.; Peron, F.; Garcia, F.P.; Nakamura, C.V.; Rodrigues-Filho, E. Copper and manganese cations alter secondary metabolism in the fungus Penicillium brasilianum. J. Braz. Chem. Soc. 2016, 27, 1444–1451. [Google Scholar] [CrossRef]
  68. Gao, S.-S.; Shang, Z.; Li, X.-M.; Li, C.-S.; Cui, C.-M.; Wang, B.-G. Secondary metabolites produced by solid fermentation of the marine-derived fungus Penicillium commune QSD-17. Biosci. Biotechnol. Biochem. 2012, 76, 358–360. [Google Scholar] [CrossRef]
  69. Rodríguez-Ortiz, R.; Mehta, B.J.; Avalos, J.; Limón, M.C. Stimulation of bikaverin production by sucrose and by salt starvation in Fusarium fujikuroi. Appl. Microbiol. Biotechnol. 2010, 85, 1991–2000. [Google Scholar] [CrossRef]
  70. Fan, B.; Parrot, D.; Blümel, M.; Labes, A.; Tasdemir, D. Influence of OSMAC-based cultivation in metabolome and anticancer activity of fungi associated with the brown alga Fucus vesiculosus. Mar. Drugs 2019, 17, 67. [Google Scholar] [CrossRef]
  71. Davis, N.D.; Diener, U.L.; Eldridge, D.W. Production of aflatoxins B1 and G1 by Aspergillus flavus in a semisynthetic medium. Appl. Environ. Microbiol. 1966, 14, 378–380. [Google Scholar]
  72. Chang, S.-C.; Wei, Y.-H.; Wei, D.-L.; Chen, Y.-Y.; Jong, S.-C. Factors affecting the production of eremofortin C and PR toxin in Penicillium roqueforti. Appl. Environ. Microbiol. 1991, 57, 2581–2585. [Google Scholar]
  73. Wang, L.; Ridgway, D.; Gu, T.; Moo-Young, M. Bioprocessing strategies to improve heterologous protein production in filamentous fungal fermentations. Biotechnol. Adv. 2005, 23, 115–129. [Google Scholar] [CrossRef]
  74. Bigelis, R.; He, H.; Yang, H.Y.; Chang, L.-P.; Greenstein, M. Production of fungal antibiotics using polymeric solid supports in solid-state and liquid fermentation. J. Ind. Microbiol. Biotechnol. 2006, 33, 815–826. [Google Scholar] [CrossRef]
  75. Guo, W.; Peng, J.; Zhu, T.; Gu, Q.; Keyzers, R.A.; Li, D. Sorbicillamines A-E, nitrogen-containing sorbicillinoids from the deep-sea-derived fungus Penicillium sp. F23−2. J. Nat. Prod. 2013, 76, 2106–2112. [Google Scholar] [CrossRef]
  76. Austin, B.; Austin, D.; Sutherland, R.; Thompson, F.; Swings, J. Pathogenicity of Vibrios to rainbow trout (Oncorhynchus mykiss, Walbaum) and Artemia nauplii. Environ. Microbiol. 2005, 7, 1488–1495. [Google Scholar] [CrossRef]
  77. Mikkelsen, H.; Lund, V.; Martinsen, L.-C.; Gravningen, K.; Schrøder, M.B. Variability among Vibrio anguillarum O2 isolates from Atlantic cod (Gadus morhua L.): Characterisation and vaccination studies. Aquaculture 2007, 266, 16–25. [Google Scholar] [CrossRef]
  78. Hickey, M.E.; Lee, J.-L. A comprehensive review of Vibrio (Listonella) anguillarum: Ecology, pathology and prevention. Rev. Aquac. 2018, 10, 585–610. [Google Scholar] [CrossRef]
  79. Li, X.; Dierckens, K.; Bossier, P.; Defoirdt, T. The impact of quorum sensing on the virulence of Vibrio anguillarum towards gnotobiotic sea bass (Dicentrarchus labrax) larvae. Aquac. Res. 2018, 49, 3686–3689. [Google Scholar] [CrossRef]
  80. Sawabe, T.; Tanaka, R.; Iqbal, M.M.; Tajima, K.; Ezura, Y.; Ivanova, E.P.; Christen, R. Assignment of Alteromonas elyakovii KMM 162T and five strains isolated from spotwounded fronds of Laminaria japonica to Pseudoalteromonas elyakovii comb. nov. and the extended description of the species. Int. J. Syst. Evol. Microbiol. 2000, 50, 256–271. [Google Scholar] [CrossRef]
  81. Shaaban, M.; Elgam, A.; Habib, E.E. Biotechnological applications of quorum sensing inhibition as novel therapeutic strategies for multidrug resistant pathogens. Microb. Pathog. 2019, 127, 138–143. [Google Scholar] [CrossRef]
  82. Dunlap, P.V. Quorum regulation of luminescence in Vibrio fischeri. J. Molec. Microbiol. Biotechnol. 1999, 1, 5–12. [Google Scholar]
  83. Kalia, V.C. Quorum sensing inhibitors: An overview. Biotechnol. Adv. 2013, 31, 224–245. [Google Scholar] [CrossRef]
  84. Gellatly, S.L.; Hancock, R.E.W. Pseudomonas aeruginosa: New insights into pathogenesis and host defenses. Pathog. Dis. 2013, 67, 159–173. [Google Scholar] [CrossRef]
  85. Rasmussen, T.B.; Skindersoe, M.E.; Bjarnsholt, T.; Phipps, R.K.; Christensen, K.B.; Jensen, P.O.; Andersen, J.B.; Koch, B.; Larsen, T.O.; Hentzer, M.; et al. Identity and effects of quorum-sensing inhibitors produced by Penicillium species. Microbiology 2005, 151, 1325–1340. [Google Scholar] [CrossRef]
  86. Ianiri, G.; Pinedo, C.; Fratianni, A.; Panfili, G.; Castoria, R. Patulin degradation by the biocontrol yeast Sporobolomyces sp. is an inducible process. Toxins 2017, 9, 61. [Google Scholar] [CrossRef]
  87. Kavanagh, F. Activities of twenty-two antibacterial substances against nine species of bacteria. J. Bacteriol. 1947, 54, 761–766. [Google Scholar]
  88. White, T.J.; Bruns, T.; Lee, S.; Taylor, J. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In PCR Protocols: A Guide to Methods and Applications, 1st ed.; Innis, M.A., Gelfand, D.H., Sninsky, J.J., White, T.J., Eds.; Academic Press Inc.: New York, NY, USA, 1990; pp. 315–322. ISBN 978-0123721815. [Google Scholar]
  89. Ramette, A. Multivariate analyses inmicrobial ecology. FEMS Microbiol. Ecol. 2007, 62, 142–160. [Google Scholar] [CrossRef]
Figure 1. Metabolome profiles of strain 1 (Ph. macrostoma) cultured under static and shaking (sh) conditions in selected media. Compounds are displayed as a function of UV absorption (nm) and intensity (mAU) vs. retention time (min). CA: casamino acids medium, WH: Wickerham’s medium. The abbreviations refer to medium (e.g., CA), strain number (1), if used, sh means shaking was used in the culture regime.
Figure 1. Metabolome profiles of strain 1 (Ph. macrostoma) cultured under static and shaking (sh) conditions in selected media. Compounds are displayed as a function of UV absorption (nm) and intensity (mAU) vs. retention time (min). CA: casamino acids medium, WH: Wickerham’s medium. The abbreviations refer to medium (e.g., CA), strain number (1), if used, sh means shaking was used in the culture regime.
Marinedrugs 17 00419 g001
Figure 2. Metabolite profile of strain 2 (C. langeronii) cultured under static and shaking (sh) conditions in CA (casamino acids) medium.
Figure 2. Metabolite profile of strain 2 (C. langeronii) cultured under static and shaking (sh) conditions in CA (casamino acids) medium.
Marinedrugs 17 00419 g002
Figure 3. Metabolite profile of strain 3 (T. harzianum) cultured under static and shaking (sh) conditions in CA (casamino acids) and WH (Wickerham’s) medium.
Figure 3. Metabolite profile of strain 3 (T. harzianum) cultured under static and shaking (sh) conditions in CA (casamino acids) and WH (Wickerham’s) medium.
Marinedrugs 17 00419 g003
Figure 4. Total peak numbers (TPN) of observed and annotated (in brackets) peaks detected in static and shaking cultures of the ten Penicillium strains (e.g., 4 (4) means four peaks were detected, and all four of them were assigned to known compounds). TPNs are sorted by a color code based on the five different liquid media (PC, CA, AS, WH, and GPY). PC: potato-carrot; CA: casamino acids; AS: artificial seawater; WH: Wickerham’s; GPY: glucose-peptone-yeast.
Figure 4. Total peak numbers (TPN) of observed and annotated (in brackets) peaks detected in static and shaking cultures of the ten Penicillium strains (e.g., 4 (4) means four peaks were detected, and all four of them were assigned to known compounds). TPNs are sorted by a color code based on the five different liquid media (PC, CA, AS, WH, and GPY). PC: potato-carrot; CA: casamino acids; AS: artificial seawater; WH: Wickerham’s; GPY: glucose-peptone-yeast.
Marinedrugs 17 00419 g004
Figure 5. Dendrograms of SMPs (secondary metabolite profiles) of all Penicillium strains cultivated under static (blue strain number) and shaking (“sh”, red strain number) conditions in (A): PC (potato-carrot), (B): CA (casamino acids), (C): AS (artificial seawater), (D): WH (Wickerham’s), and (E): GPY (glucose-peptone-yeast) medium. Strains 4–7: P. antarcticum, strain 8: P. atramentosum, strains 9–13: P. atrovenetum (see Table 1). SMPs of all strains are displayed as a function of similarity based on the Bray-Curtis dissimilarity index.
Figure 5. Dendrograms of SMPs (secondary metabolite profiles) of all Penicillium strains cultivated under static (blue strain number) and shaking (“sh”, red strain number) conditions in (A): PC (potato-carrot), (B): CA (casamino acids), (C): AS (artificial seawater), (D): WH (Wickerham’s), and (E): GPY (glucose-peptone-yeast) medium. Strains 4–7: P. antarcticum, strain 8: P. atramentosum, strains 9–13: P. atrovenetum (see Table 1). SMPs of all strains are displayed as a function of similarity based on the Bray-Curtis dissimilarity index.
Marinedrugs 17 00419 g005
Figure 6. Matrix plots of bioactivities of Penicillium spp. (strains 4–13, x-axis) grown in different media (y-axis) against V. anguillarum (A), Ps. elyakovii (B), quorum sensing (C), A. fisheri (D), biofilm formation (E), and growth of Ps. aeruginosa (F). The inhibition levels (%) are displayed in different color intensities (blue: low; red: high). Positive controls: 10 µM chloramphenicol (for A, B, and D); 10 µM furanone (for C), and 10 µM polymyxin B (for E and F). PC: potato-carrot; CA: casamino acids; AS: artificial seawater; WH: Wickerham’s; GPY: glucose-peptone-yeast.
Figure 6. Matrix plots of bioactivities of Penicillium spp. (strains 4–13, x-axis) grown in different media (y-axis) against V. anguillarum (A), Ps. elyakovii (B), quorum sensing (C), A. fisheri (D), biofilm formation (E), and growth of Ps. aeruginosa (F). The inhibition levels (%) are displayed in different color intensities (blue: low; red: high). Positive controls: 10 µM chloramphenicol (for A, B, and D); 10 µM furanone (for C), and 10 µM polymyxin B (for E and F). PC: potato-carrot; CA: casamino acids; AS: artificial seawater; WH: Wickerham’s; GPY: glucose-peptone-yeast.
Marinedrugs 17 00419 g006
Table 1. Identification of fungal strains isolated from the Baltic Z. marina sample based on morphological criteria as well as BLAST searches in GenBank of the Internal Transcribed Spacer (ITS) region. Closest relatives to fungal strains according to BLASTn are presented. L: leaf; S: sediment.
Table 1. Identification of fungal strains isolated from the Baltic Z. marina sample based on morphological criteria as well as BLAST searches in GenBank of the Internal Transcribed Spacer (ITS) region. Closest relatives to fungal strains according to BLASTn are presented. L: leaf; S: sediment.
Strain NumberOriginMorphologyNext Related StrainSimilarity (%)GenBank Accession Number
1LPhoma sp.Phoma macrostoma100MN166393
2SCladosporium sp.Cladosporium langeronii99MN166394
3STrichoderma sp.Trichoderma harzianum100MN166395
4SPenicillium sp.Penicillium antarcticum99MN172369
5SPenicillium sp.Penicillium antarcticum100MN166396
6SPenicillium sp.Penicillium antarcticum100MN166397
7SPenicillium sp.Penicillium antarcticum100MN166398
8SPenicillium sp.Penicillium atramentosum100MN166399
9SPenicillium sp.Penicillium atrovenetum100MN166400
10SPenicillium sp.Penicillium atrovenetum100MN166401
11SPenicillium sp.Penicillium atrovenetum100MN166402
12SPenicillium sp.Penicillium atrovenetum100MN166403
13SPenicillium sp.Penicillium atrovenetum100MN166404
Table 2. Dereplication of extracts from strains 1 (Ph. macrostoma), 2 (C. langeronii), and 3 (T. harzianum). n.d.: not detectable, n.k.: not known, UP: unknown peak.
Table 2. Dereplication of extracts from strains 1 (Ph. macrostoma), 2 (C. langeronii), and 3 (T. harzianum). n.d.: not detectable, n.k.: not known, UP: unknown peak.
Peak/CompoundMediumUVmax (nm)m/z [M + H]+IdentityReference
1PC, PC-sh, CA-sh, WH, WH-sh220, 245, 292, 342229.1Trioxsalen[19]
2WH-sh222, 279277.2NK-A 17E-233II[20]
3CA, CA-sh221, 289, 310211.1Agistatin D[21]
4CA, CA-sh205, 216, 222, 277295.1Orsellide D[22]
5WH215, 306311.3Orsellide C[22]
6CA, WH223, 275, 288599.3Cyl-2[23]
7CA223, 280316.1Quinolinone B[24]
8CA227, 269, 334282.1Harzianopyridone[25]
9PC, CA-sh, WH, WH-sh223, 261, 271, 281282.22′-Deoxyribofuranosylguanine; 2N-Me[26]
UP1PC-sh221, 292325.2n.k.
UP2PC, PC-sh, WH220, 275, 325, 360273.1n.k.
UP3PC, PC-sh, WH, WH-sh218, 296, 364289.1n.k.
UP4CA-sh, WH-sh232279.3n.k.
UP5WH215, 303327.3n.k.
UP6WH218, 297312.2n.k.
UP7WH220, 266, 276, 287393.1n.k.
UP8WH226245.2n.k.
UP9WH-sh222, 280, 292, 305, 319626.2n.k.
Table 3. Observed bioactivities (%) of crude extracts from non-Penicillium strains at 100 µg/mL concentration. Quorum sensing (QS) and biofilm inhibition were compared to growth inhibition of their associated microorganisms, A. fischeri and Ps. aeruginosa, separated by /. Positive control. A: 10 µM 5-(bromomethylene)-2H(5H)-furanone, B: 10 µM chloramphenicol, C: 10 µM polymyxin B.
Table 3. Observed bioactivities (%) of crude extracts from non-Penicillium strains at 100 µg/mL concentration. Quorum sensing (QS) and biofilm inhibition were compared to growth inhibition of their associated microorganisms, A. fischeri and Ps. aeruginosa, separated by /. Positive control. A: 10 µM 5-(bromomethylene)-2H(5H)-furanone, B: 10 µM chloramphenicol, C: 10 µM polymyxin B.
MediumInhibition (%) of A. fischeri
QS/Growth Inhibition
Inhibition (%) of Ps. aeruginosa
Biofilm/Growth Inhibition
Phoma macrostomaCladosporium langeroniiTrichoderma harzianumPhoma macrostomaCladosporium langeroniiTrichoderma harzianum
PC100/0100/3194/705 /240/00/3
PC-sh97/093/3747/00/60/717/3
CA99/095/090/00/160/50/3
CA-sh0/030/081/00/20/10/2
AS73/00/047/50/80/70/9
AS-sh53/077/076/012/10/311/4
WH100/1798/097/460/210/30/1
WH-sh80/095/070/010/50/10/20
GPY34 /01/1320/00/50/00/2
GPY-sh7/2940/00/354/30/10/9
Standard100A/100B100A/100B100A/100B100C/100C100C/100C100C/100C

Share and Cite

MDPI and ACS Style

Petersen, L.-E.; Marner, M.; Labes, A.; Tasdemir, D. Rapid Metabolome and Bioactivity Profiling of Fungi Associated with the Leaf and Rhizosphere of the Baltic Seagrass Zostera marina. Mar. Drugs 2019, 17, 419. https://doi.org/10.3390/md17070419

AMA Style

Petersen L-E, Marner M, Labes A, Tasdemir D. Rapid Metabolome and Bioactivity Profiling of Fungi Associated with the Leaf and Rhizosphere of the Baltic Seagrass Zostera marina. Marine Drugs. 2019; 17(7):419. https://doi.org/10.3390/md17070419

Chicago/Turabian Style

Petersen, Lars-Erik, Michael Marner, Antje Labes, and Deniz Tasdemir. 2019. "Rapid Metabolome and Bioactivity Profiling of Fungi Associated with the Leaf and Rhizosphere of the Baltic Seagrass Zostera marina" Marine Drugs 17, no. 7: 419. https://doi.org/10.3390/md17070419

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop