Abstract
Naturally occurring antagonists toward pathogens play an important role to avoid pathogen outbreaks in ecosystems, and they can be applied as biocontrol agents for crops. Lichens present long-living symbiotic systems continuously exposed to pathogens. To analyze the antagonistic potential in lichens, we studied the bacterial community active against model bacteria and fungi by an integrative approach combining isolate screening, omics techniques, and high resolution mass spectrometry. The highly diverse microbiome of the lung lichen [Lobaria pulmonaria (L.) Hoffm.] included an abundant antagonistic community dominated by Stenotrophomonas, Pseudomonas, and Burkholderia. While antagonists represent 24.5% of the isolates, they were identified with only 7% in the metagenome; which means that they were overrepresented in the culturable fraction. Isolates of the dominant antagonistic genus Stenotrophomonas produced spermidine as main bioactive component. Moreover, spermidine-related genes, especially for the transport, were identified in the metagenome. The majority of hits identified belonged to Alphaproteobacteria, while Stenotrophomonas-specific spermidine synthases were not present in the dataset. Evidence for plant growth promoting effects was found for lichen-associated strains of Stenotrophomonas. Linking of metagenomic and culture data was possible but showed partly contradictory results, which required a comparative assessment. However, we have shown that lichens are important reservoirs for antagonistic bacteria, which open broad possibilities for biotechnological applications.
Introduction
Plant pathogens and the diseases they cause are major threats to humanity. Each year we globally lose over one third of the total harvest to bacterial and fungal pathogens. The past two decades have seen an increasing number of virulent infectious diseases in plants (), and human activity is intensifying pathogen dispersal as well as reducing diversity in agricultural systems (). However, microbial diversity is a key factor in avoiding pathogen outbreaks (; ). Therefore, biocontrol of plant pathogens is a promising solution to control plant pathogens (; ; ) because it was also shown that it enhances general microbial diversity (). Naturally occurring antagonists toward plant pathogens play an important role for biocontrol approaches. In natural ecosystems, which often contain a high proportion of antagonistic microorganisms, such antagonists potentially function in stabilizing the community, but might also protect the community against pathogen outbreaks (; ; ). However, the ecology of naturally occurring antagonistic microorganisms is only partly understood and not yet exploited.
Lichens, which are classic examples of self-sustained symbioses, are interesting models for antagonism studies because within these mini-ecosystems the cooperation between microbial partners facilitates stability and longevity under extreme ecological conditions although they are often attacked by allochthonous bacteria and fungi (; ; ). While the lichen-specific structure is provided by fungal symbionts, which also is the naming component of the symbiosis, green algae, and/or cyanobacteria are incorporated into specific layers or compartments and contribute with photosynthetically fixed carbohydrates to the symbiosis (). Lichen-associated bacteria were only recently shown to be highly diverse and omics approaches have indicated that they are functional contributors to robustness of the lichen holobiome (, ). The intricate association of members of different organismal kingdoms in well-delimited and long-living symbiotic structures – as symbiotic hotspots of terrestrial life – highlights lichens as a veritable treasure chest for interorganismal communication, regulation, and bioactivity in general (; ). Conditioned by the slow growth of many lichens and difficulties in culturing the symbionts, biotechnological exploitation of lichens was lagging behind other natural resources. With the advent of modern technologies, however, the secondary metabolism and antagonistic potentials in lichens receive new impulses, and this will particularly apply to culturable bacterial partners. Although, lichens are equipped with various secondary compounds with antagonistic effects (; ; ), we hypothesize that only a diverse protective microbiome can efficiently maintain stability over longer periods to prevent pathogen attacks.
The objective of this study was to analyze the antagonistic potential of the lichen microbiome against model pathogens by a novel approach bridging metagenomic with culture techniques. Model pathogens associated with human, lichen and plant diseases were accessed to screen for a broad spectrum of antagonistic activity. Furthermore, we utilized the lung lichen Lobaria pulmonaria (L.) Hoffm., which is one of the fastest growing leaf-like lichens () and used as indicator species of undisturbed forests and air pollution (; ). We also characterized the most active as well as the most abundant lichen-associated antagonists Stenotrophomonas, which were already identified as versatile antagonists from plant origin (; ; ). Beneficial Stenotrophomonas strains produced osmoprotectans and spermidine in response to eukaryotic hosts (). In our study we applied multidisciplinary techniques to link metagenomic data with those obtained from bacterial cultures. Moreover, we could show that lichens are important reservoirs for antagonistic bacteria, which can also be used for biological control approaches to protect plants against biotic and abiotic stress.
Materials and Methods
Sampling Strategy and Isolation of Lichen-Associated Bacteria
Lichen thalli of L. pulmonaria were sampled from three different locations in Austria (Tamischbachgraben, N47°32′40′′, E14°37′35′′, Johnsbach, N47°38′07′′, E14°44′45′′, and St. Oswald ob Eibiswald, N46°44′ 50′′, E15° 04′ 26′′) after visual inspection to avoid contamination by lichenicolous fungi and other organisms. Five separate lichen thalli were sampled from each sampling site. The samples were stored on dry ice and were, shortly after, ground with mortar and pestle. A homogenate was prepared using sterile 0.85% NaCl in a 1:10 (w/v) ratio, together with a lab stomacher (BagMixer; Interscience, St Nom, France). Diluted fractions were plated on R2A agar (Carl Roth, Karlsruhe, Germany), R2A agar with 25 μg ml-1 cycloheximide, starch casein agar (SCA; ) and ISP2 agar (). Bacterial colonies were randomly picked within 5 days of incubation at room temperature (RT) and a total of 388 isolates was obtained. The isolates were stored in glycerol stocks at -70°C prior to cultivation-based experiments.
Screening of Isolates for In Vitro Antagonistic Activity Toward Particular Bacteria and Fungi
Dual-culture experiments were carried out as confrontation assays, using different media and target organisms according to and . Lichen-associated isolates were spotted on solid media pre-inoculated with Escherichia coli XL1 and Staphylococcus aureus ATCC 25923 and assessed for inhibition zones after 4 days of incubation at 30°C. Antagonistic activity against the fungus Botrytis cinerea Pers. (in-house culture collection) was tested by dual culture on Waksman agar (WA), according to and assessed after 5–7 days incubation at 20°C. Cultures of the lichen-pathogenic fungus Rhinocladiella sp. (culture collection of Lucia Muggia; Institute of Plant Sciences, University of Graz) were homogenized and re-suspended in sterile 0.85% NaCl. In the following step, 50 μL aliquots from one batch were used to inoculate each well of 24-well plates which contained solid potato dextrose agar (PDA; Carl Roth, Karlsruhe, Germany). Subsequently, 100 μL culture filtrate obtained from each lichen-associated isolate was added to particular wells. After 3 weeks of incubation, the wells were checked for growth reduction. All experiments were conducted with three independent replicates.
Amplicon Library Preparation and Co-Occurrence Analysis
Amplicon libraries obtained by were used to extract distinct taxa for additional studies. The utilized 454-pyrosequencing data was obtained from lichen samples from the same sampling sites that were used for isolation of lichen-associated bacterial cultures. Out of the 454-amplicon dataset 15 thallus samples (five for each sampling site) were used for a co-occurrence analysis. Therefore OTUs (Operational Taxonomic Units) were clustered with UCLUST () at 95% similarity (correlates with the taxonomic genus level). Mitochondrial, chloroplast, and Nostoc sequences were excluded as well as all OTUs with less than three sequences. Co-occurrence patterns were created with calculated Spearman correlations between taxa at family level (>0.6 and <-0.6; R environment version 3.1.21). Only families, which showed a correlation to Pseudomonadaceae, Xanthomonadaceae, and Burkholderiaceae were considered for further analysis and visualized as network with Cytoscape (organic layout; version 3.2.1; ). Node size within the network reflects the sequence abundance of each taxon and nodes were colored according to phylum affiliation.
16S rRNA-Based Identification of Antagonistic Bacteria and Phylogenetic Analysis
Primer pair 27F/1492r was used to amplify specific 16S rRNA gene fragments from antagonistic bacterial cultures. Subsequent sequencing and BLASTn searches within the 16S ribosomal RNA sequence database (NCBI) were conducted for identification of antagonists. These sequences were later trimmed to the hypervariable V4 region to allow alignments with 454-pyrosequencing data of the same rRNA region. QIIME 1.6.0 (Quantitative Insights Into Microbial Ecology, ) and the implemented Greengenes database () was used to search for bacterial OTUs in the corresponding amplicon dataset () that were assigned to the three most dominant proteobacterial genera within prior identified antagonists. OTUs which could not be assigned to any particular phylum or to species level within Proteobacteria were additionally analyzed with Seqmatch (RDP database) for taxonomic assignment. The representative sequences of these additionally identified OTUs were used for further phylogenetic analyses. The phylogenetic tree was constructed with V4-trimmed 16S rRNA sequences from antagonistic bacteria cultures and an amplicon subset. Sequences were aligned with MEGA6 () and processed for bootstrapped neighbor-joining with PHYLIP work package v.3.6952. Confidence levels for the internal branches were assessed by bootstrap analysis with 100 re-samplings. FigTree v.1.4.03 was used for annotation and final graphic visualization of the phylogenetic tree. All utilized 16S rRNA gene fragment sequences from isolate and amplicon sequencing were deposited at GenBank4 (accession numbers: KP739786–KP739797 and KR611621–KR611709).
Metagenomic Mining for Specific Genes of Interest
All metagenome-based analyzes were carried out on the assembled dataset described in a previous study by . CLUSTER CONTROL () was used to search with the blastn algorithm for specific spermidine synthase matches (NCBI accession numbers: NC_010943.1, NC_011071.1, NC_015947.1, and NC_017671.1) within the dataset (368,424 contigs). MEGAN (v4.70.4) was used to retrieve taxonomic classification and relevant SEED functions.
Metabolite Extraction from Bacterial Cultures Grown on Solid Medium
Bacterial cultures were washed from several densely colonized Nutrient Agar (NA; Sifin, Berlin, Germany) plates after 48 h incubation at 30°C and homogenized in 9 mL 0.85% NaCl solution. The homogenate was centrifuged for 20 min, 2,000 g at 4°C. This step was repeated two times to remove residual media from bacterial cells. The pellet was re-suspended in 2 mL ddH2O followed by centrifugation for 15 min, 18,000 g at 4°C. Precooled 90% methanol at -70°C was used for reproducible extraction and to avoid further degradation of metabolites. Subsequently, 1 mL was added to each pellet and the bacterial cells were mechanically disrupted with glass beads for 2 × 45 s at 6 m/s. Followed by a final centrifugation step for 15 min, 18,000 g at 4°C, 100 μl of each supernatant was collected and immediately placed in a deep freezer at -70°C until further analysis. Three independent biological replicates were prepared for each isolate.
Preparation of Culture Supernatants for Spermidine Quantification
Bacterial cultures were used to prepare overnight cultures (ONC) in fluid Nutrient Broth II (NBII; Sifin, Berlin, Germany) medium. These ONCs were used to inoculate 50 mL NBII flasks, which were then incubated at 30°C, 120 rpm for 48 h. In the following step, 2 mL aliquots were taken from the cultures and centrifuged for 20 min, 18,000 g at 4°C. The supernatants were filtered with 0.25 μm filters and immediately placed in a deep freezer at -70°C until further analysis. Three independent biological replicates were prepared for each isolate.
Quantification of Specific Bacterial Metabolites with High Resolution Mass Spectrometry
Samples were analyzed in nine biological/technical replicates with a combined HPLC-hybrid quadrupole-orbitrap mass spectrometer (Q Exactive; Thermo Scientific, Bremen, Germany). A Luna 5u NH2 100A 250 × 4.6 column (Phenomenex, Aschaffenburg, Germany) was used to separate different metabolites from the cell extracts. Formic acid (0.1%, v/v) in acetonitrile was used as solvent A and aqueous formic acid (0.1%, v/v) as solvent B. Starting conditions for the gradient elution were 10% A and 90% B. The conditions were gradually changed to 80% A and 20% B within 15 min. This step was followed by 5 min at 10% A and 90% B for readjustment to initial conditions. The eluent flow was maintained at 0.8 mL/min together with a column temperature of 25°C. Sample analysis was carried out with negative ion ESI detection. ESI conditions were set to 3.2 kV spray voltage and 350°C capillary temperature. Scans were recorded in the range 100.0–300.0 m/z with the AGC target set to 500,000 and maximal accumulation time of 200 ms. The resolution was adjusted to 200,000. Altering full MS-SIM and targeted MS2 cycles were employed and a specific inclusion mass of 146.16517 amu was selected. Standard calibration was obtained with 0, 0.02, 0.03, 0.04, 0.05, 0.1, and 0.2 μM spermidine standard (Duchefa Biochemie, Haarlem, The Netherlands) diluted in 0.2 mM HCL.
Plant Growth Experiments with Stenotrophomonas-Primed Seeds
Overnight cultures with selected Stenotrophomonas isolates were used to inoculate main cultures in fluid NBII. In addition to the lichen-associated isolates, a plant-associated isolate Stenotrophomonas rhizophila P69 (; ) was also utilized for comparisons. After 2 h of growth at 30°C and 120 rpm the fluid cultures were diluted to 5 × 106–1 × 107 cells per mL in sterile 0.85% NaCl solution. Tomato (Solanum lycopersicum L. cv. Kremser Perle; Austrosaat, Graz, Austria) seeds were surface sterilized with 4% NaHClO for 10 min followed by drying at RT. The sterilized seeds were put into the respective bacterial suspensions and incubated for 4 h at 120 rpm and RT. The control samples were put in 0.85% NaCl solution without bacteria. Ground and homogenized seeds from each inoculum were plated on NA to test priming efficiency after the incubation time. The remaining seeds were planted in sterile soil (150 g/tray) with vermiculite (1:3 ratio) and watered with 30 mL sterile H2O. Beside the control with non-primed seeds an additional control with P69-primed seeds was added. Therefore, 60 μL 5-sec-butyl-2,3-dimethylpyrazine (Sigma–Aldrich, Steinheim, Germany) was supplemented into 30 mL sterile H2O used for irrigation after planting. The closed trays were placed without further irrigation for 2 weeks into a greenhouse with 12 h day/night cycles and a constant temperature of 24°C. Root (n = 64) and stem (n = 66) lengths of the plantlets were assessed separately for all samples.
Statistical Analysis
Statistical analyses were performed with SPSS v.20.0.0 (SPSS Inc, Chicago, IL, USA). Data were tested for normal distribution with the Kolmogorov–Smirnov test. Sets with normally distributed data were analyzed with univariate ANOVA and Duncan tests at p < 0.05. The t-test was employed for statistical analysis of data that was not normally distributed (p < 0.05 and p < 0.1).
Results
Antagonistic Bacteria within the Lobaria Microbiome
Cultivable bacteria, which have been isolated from 15 L. pulmonaria samples, were tested in dual-culture assays against the bacterial model pathogens E. coli, S. aureus, the plant pathogen B. cinerea and the lichenicolous fungus Rhinocladoniella sp. to determine the general antagonistic potential. In these experiments, lichen-associated antagonists were shown to primarily target lichen and plant pathogenic fungi (20.1% of all isolates) while a lower proportion was directed against bacterial model pathogens (7.7% of all isolates). A total of 95 isolates (24.5%) showed inhibition of growth of at least one of the model pathogens (Figure 1). Singular antagonists (active against only one target microorganism) were dominated by Stenotrophomonas spp. (31% of singular antagonists) followed by Pseudomonas spp. (19%) and Burkholderia spp. (12%; Supplementary Figure S1). Dual antagonists (active against two microorganisms) comprised nine isolates. Stenotrophomonas and Micrococcus were represented by three isolates and Chryseobacterium, Microbacterium, and Paenibacillus by only one isolate, respectively. Nine bacterial strains inhibited the growth of either three or four model pathogens simultaneously. These cultures were identified at genus level as Bacillus (five isolates), Micrococcus (one isolate), and Paenibacillus (three isolates). A complete taxonomic breakdown for all identified antagonists was visualized in Figure 2.
FIGURE 1
FIGURE 2
Comparison of the hypervariable V4 rRNA region from the most abundant proteobacterial antagonistic isolates: Stenotrophomonas spp., Pseudomonas spp., and Burkholderia spp. with filtered OTUs from an amplicon library constructed with Lobaria samples from the same sampling sites, revealed high homology of sequences within the same genus (Supplementary Figure S2). Specific branches (bootstrap values > 70%) were detected for amplicon- and isolate-based sequences. OTUs that were assigned to Stenotrophomonas sp., Pseudomonas sp., and Burkholderia sp. comprised, respectively, 0.06, 0.56, and 0.09% of all analyzed OTUs in the amplicon library. Interestingly three of the highly active antagonistic genera (Bacillus, Micrococcus, and Paenibacillus) were not substantially represented in the amplicon library.
A co-occurrence pattern between different taxa at family level was created with the calculated Spearman correlations based on 15 lichen thallus samples. In total, 24 correlations between the families Pseudomonadaceae, Burkholderiaceae, and Xanthomonadaceae to other taxa within the microbiome could be detected and were visualized as a co-occurrence network (Figure 3). Most correlations (10 out of 24) were found within the phylum Proteobacteria followed by Actinobacteria (3). The strongest positive correlations (Spearman correlation > ±0.7) showed Xanthomonadaceae with Pseudomonadaceae and Alteromonadaceae, all assigned to Gammaproteobacteria, whereas the strongest negative correlations were found between Burkholderiaceae (Betaproteobacteria) and Phyllobacteriaceae (Alphaproteobacteria) and a family within the class Chloracidobacteria (Acidobacteria) which was not further classified in the utilized database.
FIGURE 3
In addition, the abundance of antagonistic taxa was extracted from the Lobaria metagenome. The proportion of retrievable antagonistic genera was determined for Stenotrophomonas (0.22% of all bacteria within the metagenome), Pseudomonas (1.14%), Burkholderia (2.81%), Xanthomonas (0.43%), Nocardiodes (0.10%), Rhodococcus (0.18%), Bacillus (0.08%), and Staphylococcus (0.02%). The remaining antagonistic taxa could not be retrieved at genus level. However, they comprised at family level together 2.14% of all bacteria. Altogether, antagonistic taxa comprised 7.12% of the total bacterial community. Neither the genus Cellulomonas nor the family Cellulomonadaceae was present in the assembled metagenomic dataset.
Spermidine Production In Vitro and Spermidine-Related Genes within the Lobaria-Associated Metagenome
Genes coding for spermidine synthases were analyzed from the Lobaria metagenome and taxonomically assigned. In addition, spermidine production and secretion was analyzed in vitro.
Seven antagonistic Stenotrophomonas sp. isolates were cultivated on solid agar plates and in liquid media prior to the extraction of spermidine. The detection limit for spermidine on the utilized instruments was determined to be <30 nM. Externalized spermidine levels detected in liquid cultivation media were in the range between 8.2 and 10.5 μmol/g fresh weight. Extracellular spermidine concentration differences between utilized Stenotrophomonas isolates were not statistically significant (Supplementary Figure S3). Conversely, the same isolates were shown to contain different internal spermidine concentrations after cultivation on solid media (Figure 4). The lowest internal spermidine concentration was found to be 168 nmol/g fresh weight, which was around fourfold lower than the highest observed concentration. Isolate 165P3RAB was found to contain significantly higher spermidine concentrations than all other isolates. In contrast, Stenotrophomonas isolate 329P5R contained the lowest spermidine concentrations.
FIGURE 4
For BLASTn searches, reference sequences for spermidine synthases from four different Stenotrophomonas strains were utilized. Additionally, SEED assignments were searched for related functions. Stenotrophomonas-specific contigs that contain known spermidine synthases were not found in the utilized metagenome, while three other bacteria-derived contigs were present. Two spermidine synthase contigs that were retrieved with SEED-based analysis were assigned to Proteobacteria. One contig was assigned to Burkholderiaceae and the other contig to Acetobacteraceae. A third contig could not be assigned to any taxon. Conversely, spermidine putrescine transporter permeases were more abundant in the metagenome. A total of 50 contigs were assigned to this specific transporter protein (Supplementary Figure S4). More than a half of these contigs were assigned to bacteria (58%), while 42% remained unassigned to a specific kingdom. The hits were predominantly associated with Proteobacteria (52%) and more specifically to Rhizobiales (28%).
Stenotrophomonas Treatments Increased Plant Growth of Tomato Under Greenhouse Conditions
Tomato (Solanum lycopersicum L.) seeds were inoculated with three lichen-associated and one plant-associated Stenotrophomonas isolate and stress protecting agent as reference () to analyze the effect of bacterial inoculants on plants. The primed seeds were grown with limited irrigation for 2 weeks in sterile soil. Two control types were implemented to evaluate growth promotion effects by the inoculants. One control (P69_Py) was supplemented with 60 μL 5-sec-butyl-2,3-dimethylpyrazine per tray during the initial irrigation of Stenotrophomonas P69-primed seeds. This heterocyclic compound was found to limit the growth of Stenotrophomonas isolates in previous experiments. Correspondingly, the growth of P69-primed samples that were treated with 5-sec-butyl-2,3-dimethylpyrazine (R_P69_Py and S_P69_Py), was similar to non-primed control samples (Figure 5). Also, the plant growth was not enhanced by isolate 165P3RAB, which was shown to contain the highest internal spermidine concentrations when compared to both implemented controls. Stenotrophomonas isolates 329P5R, 401P2, and P69 enhanced the plant growth significantly when compared to both controls. These isolates were shown to produce low internal spermidine concentrations in previous experiments. Treatments with the lichen-associated isolate 401P2 and the plant-associated isolate P69 resulted in similar plant growth.
FIGURE 5
Discussion
The lichen symbiosis was discovered as reservoir for antagonistic bacteria. Interestingly, it was possible to transfer selected isolates from lichens to cultivated plants while maintaining beneficial effects. In addition, we have shown the usability as well as limits of various applied techniques to efficiently screen for specific characteristics and how to reasonably couple classic microbiology with high-end techniques in a comprehensive approach. Starting from a culture collection and dual-culture experiments to screen for active antagonists, the approach was expanded with detailed specification of continuously filtered isolates.
The microbiome involved in the lichen symbiosis is highly diverse () and was identified as bioresource for antagonistic bacteria. L. pulmonaria is predominately colonized by Alphaproteobacteria (), in particular by various members of Rhizobiales (), which harbor many bacterial genera known for a beneficial host-microbe interaction especially with plants. Interestingly, all antagonistic genera identified for lichens – Stenotrophomonas, Pseudomonas, Burkholderia, Micrococcus, Chryseobacterium, Microbacterium, and Paenibacillus – are well-known from plant studies (; ; Rybakova et al., submitted). This is an interesting finding because it shows that these bacteria have the same redundant function independent of the habitat. This observation is underlined by the greenhouse experiments, which have shown that lichen-associated antagonists are active on plants. This also supports the hypothesis that natural ecosystems are interesting reservoirs for biotechnologically relevant bacteria. The present study depicts that cultivable bacterial taxa with lower occurrence on lichens are mainly responsible for the protection against biotic disturbance. A highly diversified bacterial microbiome enhances the available functional repertoire, which might play a crucial role for the stability and longevity of the lichen symbiosis. Comparison of isolate-derived 16S rRNA gene fragments and amplicon-based sequences of abundant antagonists has indicated that a high proportion of Burkholderia spp., Pseudomonas spp., and Stenotrophomonas spp. can be retrieved from lichen symbioses by cultivation experiments on conventional growth media. It was also demonstrated that several isolated antagonists, such as Bacillus, Paenibacillus, and Micrococcus, were not detectable in the amplicon library but partially in the metagenomic dataset. The most reasonable explanation is that these antagonistic species occur with low abundance within this lichen microbiome and therefore these species might be below the detection limit of the utilized 454 pyrosequencing approach. Other methods with higher coverage might be more suitable to uncover all present bacterial colonizers. Further studies that address this question should preferably subject less multiplexed samples to high-throughput sequencing platforms to obtain a higher read number per sample. This would allow more accurate characterizations of the rare microbial population.
The majority of the antagonistic isolates was assigned to the genus Stenotrophomonas. These bacteria have been reported to protect plants against unaffordable conditions like drought and elevated salinity by exudation of protective compounds like spermidine and different osmolytes (; ). Corresponding to this, polyamines which also include spermidine were shown to be involved in plant response to abiotic stress in prior studies (; ). Environmental strains of S. maltophilia and S. rhizophila were reported to exert a certain degree of tolerance toward salinity of up to 9% (w/v) NaCl which was correlated with the ability to produce the osmolytes trehalose and glucosylglycerol (). Even though Stenotrophomonas-specific spermidine synthases were not present in the analyzed metagenome, it was evident that utilization of spermidine is widely distributed among various detected organisms. Since lichens presents a habitat that is frequently subjected to drought, the association with bacteria having protective properties appears favorable. However, lichens themselves have mechanisms to account for desiccation, which also includes osmolytes (). Moreover, protection mechanisms against oxidative stress-related damage act in a mutual manner among the eukaryotic partners (). Therefore, we consider stress-protective functions of lichen-associated bacteria to act as an enhancer, which might react more flexibly to local fluctuations of the conditions than their hosts.
Although synthesized and excreted in different amounts, for all selected Stenotrophomonas isolates, in vitro spermidine production was detected. According to the result, we assume that the function of lichen-associated bacteria includes assistance in the protection against pathogens as well as against damage caused by desiccation. Stenotrophomonas is a well-known antagonist of plant-associated origin () and connected with a beneficial effect on plant hosts (; ; ). Stenotrophomonas strains might also have the same function independent of the habitat. However, in the last two decades, they have received additional attention for opportunistic infections in humans (). It is difficult to identify specific factors of pathogenicity of the virulent isolates but the ability of persistence, resistance and survival – also essential to colonize lichens – allowed a colonization of immunocompromised patients with predisposition. Similar to the plant rhizosphere (), lichens may act as a reservoir for facultative human pathogenic bacteria, or close relatives thereof.
The host-microbiome balance as well as indigenous diversity is essential for functional stability in ecosystems. This balance also depends on the mutual effects among bacteria within the microbiome. In our co-occurrence analyses we found indications for positive and negative correlations among bacterial groups on samples of the same host lichen symbiosis in the same habitat, which is a strong indication of antagonistic and synergistic effect in the lichen habitat. Further analyses are required for clarifying the mechanisms responsible for these effects, as these might involve direct interactions or diffusible metabolites in the system or both. The present results have already shown that lichen symbioses are valuable bioresources to discover bacteria with antagonistic potential and we suggest that a systematic screening of a broader range of lichens may be useful for finding biocontrol solutions that are specifically tailored for ecologically different plant habitats.
Statements
Acknowledgments
We appreciate the help of Armin Erlacher (Graz) with statistical analyses and his generous contributions to the study. Christoph W. Sensen (Graz) and Jung Soh (Calgary) provided valuable support with metagenomics analyses. We also want to thank Angelika Schäfer (Graz) and Kenneth B. Jensen (Graz) for their support during preparation and realization of LC/MS experiments. We are grateful to Lucia Muggia (Graz) for providing the culture of Rhinocladiella sp. We also would like to thank Rocel A. Ortega (Graz) for proof-reading the manuscript and her helpful remarks. This work was supported by a grant of the Austrian Science Fund (FWF) to GB and MG (FWF Project I882). Parts of the study were conducted in the Central Lab ‘Environmental, Plant and Microbial Metabolomics’ funded by NAWI Graz.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb.2015.00620
Footnotes
2.^http://evolution.genetics.washington.edu/phylip
References
1
AlaviP.StarcherM. R.ThallingerG. G.ZachowC.MüllerH.BergG. (2014). Stenotrophomonas comparative genomics reveals genes and functions that differentiate beneficial and pathogenic bacteria.BMC Genomics15:482. 10.1186/1471-2164-15-482
2
AlaviP.StarcherM. R.ZachowC.MüllerH.BergG. (2013). Root-microbe systems: the effect and mode of interaction of stress protecting agent (SPA) Stenotrophomonas rhizophila DSM14405T.Front. Plant Sci.4:141. 10.3389/fpls.2013.00141
3
AlcazarR.MarcoF.CuevasJ. C.PatronM.FerrandoA.CarrascoP.et al (2006). Involvement of polyamines in plant response to abiotic stress.Biotechnol. Lett.281867–1876. 10.1007/s10529-006-9179-3
4
AschenbrennerI. A.CardinaleM.BergG.GrubeM. (2014). Microbial cargo: do bacteria on symbiotic propagules reinforce the microbiome of lichens?.Environ. Microbiol.163743–3752. 10.1111/1462-2920.12658
5
BatesS. T.CropseyG. W.CaporasoJ. G.KnightR.FiererN. (2011). Bacterial communities associated with the lichen symbiosis.Appl. Environ. Microbiol.771309–1314. 10.1128/AEM.02257-10
6
BerendsenR. L.PieterseC. M.BakkerP. A. (2012). The rhizosphere microbiome and plant health.Trends Plant Sci.17478–486. 10.1016/j.tplants.2012.04.001
7
BergG. (2015). Beyond borders: investigating microbiome interactivity and diversity for advanced biocontrol technologies.Microb. Biotechnol.85–7. 10.1111/1751-7915.12235
8
BergG.EberlL.HartmannA. (2005). The rhizosphere as a reservoir for opportunistic human pathogenic bacteria.Environ. Microbiol.71673–1685. 10.1111/j.1462-2920.2005.00891.x
9
BergG.EgamberdievaD.LugtenbergB.HagemannM. (2010). “Symbiotic plant–microbe interactions: stress protection, plant growth promotion, and biocontrol by Stenotrophomonas,” inSymbioses and StressedsSeckbachJ. M. G.GrubeM. (Netherlands: Springer) 445–460. 10.1007/978-90-481-9449-0_22
10
BergG.MartenP.BallinG. (1996). Stenotrophomonas maltophilia in the rhizosphere of oilseed rape - occurrence, characterization and interaction with phytopathogenic fungi.Microbiol. Res.15119–27. 10.1016/S0944-5013(96)80051-6
11
BergG.MartinezJ. L. (2015). Friends or foes: can we make a distinction between beneficial and harmful strains of the Stenotrophomonas maltophilia complex?Front. Microbiol.6:241. 10.3389/fmicb.2015.00241
12
BergG.RoskotN.SteidleA.EberlL.ZockA.SmallaK. (2002). Plant-dependent genotypic and phenotypic diversity of antagonistic rhizobacteria isolated from different Verticillium host plants.Appl. Environ. Microbiol.683328–3338. 10.1128/AEM.68.7.3328-3338.2002
13
BergG.ZachowC.MüllerH.PhillipsJ.TilcherR. (2013). Next-generation bio-products sowing the seeds of success for sustainable agriculture.Agronomy3648–656. 10.3390/agronomy3040648
14
BoustieJ.GrubeM. (2005). Lichens—a promising source of bioactive secondary metabolites.Plant Genet. Resour.3273–287. 10.1079/PGR200572
15
BoustieJ.TomasiS.GrubeM. (2011). Bioactive lichen metabolites: alpine habitats as an untapped source.Phytochem. Rev.10287–307. 10.1007/s11101-010-9201-1
16
CaporasoJ. G.KuczynskiJ.StombaughJ.BittingerK.BushmanF. D.CostelloE. K.et al (2010). QIIME allows analysis of high-throughput community sequencing data.Nat. Methods7335–336. 10.1038/nmeth.f.303
17
DeSantisT. Z.HugenholtzP.LarsenN.RojasM.BrodieE. L.KellerK.et al (2006). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB.Appl. Environ. Microbiol.725069–5072. 10.1128/AEM.03006-05
18
EdgarR. C. (2010). Search and clustering orders of magnitude faster than BLAST.Bioinformatics262460–2461. 10.1093/bioinformatics/btq461
19
EgamberdievaD.KucharovaZ.DavranovK.BergG.MakarovaN.AzarovaT.et al (2011). Bacteria able to control foot and root rot and to promote growth of cucumber in salinated soils.Biol. Fertil. Soils.47197–205. 10.1007/s00374-010-0523-3
20
ErlacherA.CardinaleM.GrubeM.BergG. (2015a). Biotic stress shifted structure and abundance of enterobacteriaceae in the lettuce microbiome.PLoS ONE10:e0118068. 10.1371/journal.pone.0118068
21
ErlacherA.CernavaT.CardinaleM.SohJ.SensenC. W.GrubeM.et al (2015b). Rhizobiales as functional and endosymbiontic members in the lichen symbiosis of Lobaria pulmonaria L.Front. Microbiol.6:53. 10.3389/fmicb.2015.00053
22
FisherM. C.HenkD. A.BriggsC. J.BrownsteinJ. S.MadoffL. C.McCrawS. L.et al (2012). Emerging fungal threats to animal, plant and ecosystem health.Nature484186–194. 10.1038/nature10947
23
GreenT. G. A.SanchoL. G.PintadoA. (2011). “Ecophysiology of desiccation/rehydration cycles in mosses and lichens,” inPlant Desiccation ToleranceedsLüttgeU.BeckE.BartelsD. (Berlin: Springer) 89–120. 10.1007/978-3-642-19106-0_6
24
GrubeM.CardinaleM.de CastroJ. V.MüllerH.BergG. (2009). Species-specific structural and functional diversity of bacterial communities in lichen symbioses.ISME J.31105–1115. 10.1038/ismej.2009.63
25
GrubeM.CernavaT.SohJ.FuchsS.AschenbrennerI.LassekC.et al (2015). Exploring functional contexts of symbiotic sustain within lichen-associated bacteria by comparative omics.ISME J.9412–424. 10.1038/ismej.2014.138
26
HaasD.DéfagoG. (2005). Biological control of soil-borne pathogens by fluorescent pseudomonads.Nat. Rev. Microbiol.3307–319. 10.1038/nrmicro1129
27
KrannerI.CramW. J.ZornM.WornikS.YoshimuraI.StabentheinerE.et al (2005). Antioxidants and photoprotection in a lichen as compared with its isolated symbiotic partners.Proc. Natl. Acad. Sci. U.S.A.1023141–3146. 10.1073/pnas.0407716102
28
KüsterE.WilliamsS. T. (1964). Selection of media for isolation of Streptomycetes.Nature202928–929. 10.1038/202928a0
29
LawreyJ. D. (2009). “Chemical defense in lichen symbioses,” inDefensive Mutualism in Microbial SymbiosisedsWhiteJ. F.Jr.TorresM. S. (Boca Raton, FL: CRC Press) 167–176. 10.1201/9781420069327.pt3
30
LawreyJ. D.DiederichP. (2003). Lichenicolous fungi: interactions, evolution, and biodiversity.Bryologist10681–120. 10.1639/0007-2745(2003)106[0080:LFIEAB]2.0.CO;2
31
LiuJ. H.KitashibaH.WangJ.BanY.MoriguchiT. (2007). Polyamines and their ability to provide environmental stress tolerance to plants.Plant Biotechnol.24117–126. 10.5511/plantbiotechnology.24.117
32
MacDonaldA.CoxsonD. (2013). A comparison of Lobaria pulmonaria population structure between subalpine fir (Abies lasiocarpa) and mountain alder (Alnus incana) host-tree species in British Columbia’s inland temperate rainforest.Botany91535–544. 10.1139/cjb-2013-0025
33
MendesR.KruijtM.de BruijnI.DekkersE.van der VoortM.SchneiderJ. H.et al (2012). Deciphering the rhizosphere microbiome for disease-suppressive bacteria.Science3321097–1100. 10.1126/science.1203980
34
MinkwitzA.BergG. (2001). Comparison of antifungal activities and 16S ribosomal DNA sequences of clinical and environmental isolates of Stenotrophomonas maltophilia.J. Clin. Microbiol.39139–145. 10.1128/JCM.39.1.139-145.2001
35
MushegianA. A.PetersonC. N.BakerC. C. M.PringleA. (2011). Bacterial diversity across individual lichens.Appl. Environ. Microbiol.774249–4252. 10.1128/AEM.02850-10
36
NashT. H.III. (2008). “Introduction,” inLichen Biologyed.NashT. H. (Cambridge: Cambridge University Press) 1–3. 10.1017/cbo9780511790478.002
37
OksanenI. (2006). Ecological and biotechnological aspects of lichens.Appl. Microbiol. Biotechnol.73723–734. 10.1007/s00253-006-0611-3
38
OpeltK.BergC.BergG. (2007). The bryophyte genus Sphagnum is a reservoir for powerful and extraordinary antagonists and potentially facultative human pathogens.FEMS Microbiol. Ecol.6138–53. 10.1111/j.1574-6941.2007.00323.x
39
Ribbeck-BuschK.RoderA.HasseD.De BoerW.MartínezJ. L.HagemannM.et al (2005). A molecular biological protocol to distinguish potentially human pathogenic Stenotrophomonas maltophilia from plant-associated Stenotrophomonas rhizophila.Environ. Microbiol.71853–1858. 10.1111/j.1462-2920.2005.00928.x
40
RoseF. (1976). “Lichenological indicators of age and environmental continuity in woodlands,” inLichenology: Progress and ProblemsedsBrownD. H.HawksworthD. L.BaileyR. H. (London: Academic Press) 279–307.
41
RyanR. P.MonchyS.CardinaleM.TaghaviS.CrossmanL.AvisonM. B.et al (2009). The versatility and adaptation of bacteria from the genus Stenotrophomonas.Nat. Rev. Microbiol.7514–525. 10.1038/nrmicro2163
42
SaitoR.SmootM. E.OnoK.RuscheinskiJ.WangP. L.LotiaS.et al (2012). A travel guide to Cytoscape plugins.Nat. Methods91069–1076. 10.1038/nmeth.2212
43
ScheideggerC. (1995). Early development of transplanted isidioid soredia of Lobaria pulmonaria in an endangered population.Lichenologist27361–374. 10.1017/s0024282995000478
44
SchmidF.MoserG.MüllerH.BergG. (2011). Functional and structural microbial diversity in organic and conventional viticulture: organic farming benefits natural biocontrol agents.Appl. Environ. Microbiol.772188–2191. 10.1128/AEM.02187-10
45
ShirlingE. T.GottliebD. (1966). Method for characterization of Streptomyces species.Int. J. Syst. Bacteriol.16313–340. 10.1099/00207713-16-3-313
46
StockerG.RiederD.TrajanoskiZ. (2004). Cluster Control: a web interface for distributing and monitoring bioinformatics applications on a Linux cluster.Bioinformatics20805–807. 10.1093/bioinformatics/bth014
47
TamuraK.StecherG.PetersonD.FilipskiA.KumarS. (2013). MEGA6: molecular evolutionary genetics analysis version 6.0.Mol. Biol. Evol.302725–2729. 10.1093/molbev/mst197
48
van ElsasJ. D.ChiurazziM.MallonC. A.ElhottovāD.KrištufekV.SallesJ. F. (2012). Microbial diversity determines the invasion of soil by a bacterial pathogen.Proc. Natl. Acad. Sci. U.S.A.1091159–1164. 10.1073/pnas.1109326109
49
WolfA.FritzeA.HagemannM.BergG. (2002). Stenotrophomonas rhizophila sp. nov., a novel plant-associated bacterium with antifungal properties.Int. J. Sys. Evol. Microbiol.521937–1944. 10.1099/ijs.0.02135-0
50
ZachowC.TilcherR.BergG. (2008). Sugar beet-associated bacterial and fungal communities show a high indigenous antagonistic potential against plant pathogens.Microb. Ecol.55119–129. 10.1007/s00248-007-9257-7
Summary
Keywords
lichen, antagonistic bacteria, plant growth promotion, Stenotrophomonas, spermidine
Citation
Cernava T, Müller H, Aschenbrenner IA, Grube M and Berg G (2015) Analyzing the antagonistic potential of the lichen microbiome against pathogens by bridging metagenomic with culture studies. Front. Microbiol. 6:620. doi: 10.3389/fmicb.2015.00620
Received
03 April 2015
Accepted
08 June 2015
Published
22 June 2015
Volume
6 - 2015
Edited by
Jesús Mercado-Blanco, Consejo Superior de Investigaciones Científicas, Spain
Reviewed by
Yusuke Saijo, Max Planck Institute for Plant Breeding Research, Germany; Diogo Neves Proença, University of Aveiro and CEMUC – University of Coimbra, Portugal
Copyright
© 2015 Cernava, Müller, Aschenbrenner, Grube and Berg.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Gabriele Berg, Institute of Environmental Biotechnology, Graz University of Technology, Petersgasse 12, 8010 Graz, Austria, gabriele.berg@tugraz.at
This article was submitted to Plant-Microbe Interaction, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.