Abstract
The temperate seagrass species eelgrass Zostera marina can be infected by the wasting disease pathogen Labyrinthula zosterae, which is believed to have killed about 90% of the seagrass in the Atlantic Ocean in the 1930s. It is not known why this opportunistic pathogen sometimes becomes virulent, but the recurrent outbreaks may be due to a weakening of the Z. marina plants from adverse environmental changes. This study investigated the individual and interactive effects of multiple extrinsic factors (temperature, light, and tissue damage) on the host-pathogen interaction between Z. marina and L. zosterae in a fully crossed infection experiment. The degree of infection was measured as both lesion coverage and L. zosterae cell concentration. We also investigated if the treatment factors affect the chemical defense of the host, measured as the inhibitory capacity of seagrass extracts in bioassays with L. zosterae. Finally, gene expression of a set of targeted genes was quantified in order to investigate how the treatments change Z. marina’s response to infection. Light had a pronounced effect on L. zosterae infection measured as lesion coverage, where reduced light conditions increased lesions by 35%. The response to light on L. zosterae cell concentration was more complex and showed significant interaction with the temperature treatment. Cell concentration was also significantly affected by physical damage, where damage surprisingly resulted in a reduced cell concentration of the pathogen. No treatment factor caused detectable decrease in the inhibitory capacity of the seagrass extracts. There were several interactive effects between L. zosterae infection and the treatment factors on Z. marina growth, and on the expression of genes associated with immune defense, phenol synthesis and primary metabolism, showing that the molecular reaction toward L. zosterae infection depends on prevailing environmental conditions. Our study shows that individual or interactive effects of light, temperature and tissue damage can affect multiple aspects of host-pathogen interactions in seagrasses. These results highlight the complexity of marine host-pathogen systems, showing that more multi-factorial investigations are needed to gain a better understanding of disease in marine plants under different environmental conditions.
Introduction
Seagrass meadows are habitats dominated by one or several species of seagrasses, and they constitute one of the most important types of coastal marine ecosystems globally (). The meadows act as nursery grounds for commercially important fish and shellfish species (), function as coastal protection (), and sequester large amounts of carbon (). During the last decades seagrass populations have declined sharply worldwide (Waycott et al., 2009), due to multiple stressors such as sediment and nutrient runoff, physical disturbance, invasive species, commercial fishing practices, aquaculture, overgrazing, algal blooms, global warming, and disease ().
The dominant seagrass species in the northern hemisphere, Zostera marina L. or eelgrass, is regularly infected by the endophytic pathogen Labyrinthula zosterae Porter and Muehlstein (1991) (; ; ). The pathogen is a colony-forming stramenopile that causes black necrotic lesions on its host plant, and is likely the causative organism responsible the wasting disease outbreaks in the 1930s, when about 90% of the Atlantic populations of Z. marina vanished over a span of a few years (reviewed by Sullivan et al., 2013). Labyrinthula spp. has since been detected in areas with seagrass decline on a smaller scale in e.g., New Zealand in the 1960s (), and on the north-eastern coast of the United States in the 1980s (). Presently, L. zosterae and other pathogenic Labyrinthula spp. are omnipresent in seagrass populations around the world (), but it is unknown why they sometimes become virulent pathogens. The recurrent virulent outbreaks of L. zosterae may be due to the weakening of the Z. marina plants because of adverse environmental changes (e.g., Tutin, 1938). Stressors caused by a changing environment, e.g., heat stress, reduced light and salinity, and lowered pH, are all factors that potentially can increase pathogen susceptibility of an organism (; ).
It is well known that plants resist pathogens through production of chemical defenses (; ), but in comparison to many terrestrial plants the chemical ecology of seagrasses is vastly understudied. Many marine organisms, including seagrasses, are chemically defended against herbivores and pathogens (; ; Zidorn, 2016), and extracts from Z. marina can inhibit L. zosterae growth (). Furthermore, external stressors can have an impact on seagrass metabolite production, e.g., increased temperature and reduced light intensity both affect the concentration of phenolic acids in Z. marina (Vergeer et al., 1995). However, it is unknown if and how environmental stressors also affect the production of the compounds that are inhibitory against L. zosterae (). Plant resistance to pathogens consists of multiple defensive strategies, often including the activation of several defense-related genes (), and laboratory experiments show that L. zosterae infection can alter Z. marina’s expression of genes connected to a metabolic defense cascade ().
The influence of temperature on wasting disease dynamics has been poorly investigated, and most data on temperature effects are correlational (Rasmussen, 1977; ; ). For example, mass mortality of eelgrass caused by wasting disease in Denmark during the 1930s was connected to unusually high water temperatures (Rasmussen, 1977). Furthermore, temperature correlates with L. zosterae abundance in natural Z. marina populations, and the pathogen is more active during warmer months of the year (). Increased summer temperature has, however, been reported to reduce lesion size caused by a pathogenic Labyrinthula sp. in laboratory experiments on the Mediterranean seagrass species Posidonia oceanica Delile (1813) and Cymodocea nodosa Acherson (1870) (; ).
A correlation between L. zosterae infection rate and water depth has been detected in natural Z. marina meadows (; ). Shoots growing in shallower waters, and consequently with higher light availability, have a higher degree of L. zosterae infection compared to deeper populations (; ). In contrast, petri dish experiments show increased lesions coverage on Z. marina leaves subjected to complete darkness compared to pieces kept at a diel light cycle (). However, laboratory experiments with whole Z. marina shoots and artificial lights found no change in lesion size or L. zosterae cell concentrations (). Additional experimental studies at ecologically relevant light levels and with intact Z. marina shoots are needed to determine how light can affect L. zosterae infection in Z. marina.
Tissue wounding caused by physical/mechanical damage is correlated to increased risk of pathogen infection in terrestrial plants (), but this issue has not been addressed for seagrasses with one exception. In a study on Thalassia testudinum Koenig (1805), in Florida Bay, showed that plants subjected to simulated herbivory were more likely to be infected by Labyrinthula sp. The Swedish northwestern coast is a popular recreational boating area, and 74% of the docks and 91% of the marinas in this area are located in Z. marina meadows (). The effect of boat traffic and risk of damage by boat propellers on Z. marina plants is therefore potentially high, and shoots with damaged leaves are frequent in some highly trafficked areas (personal observations by ).
Most previous studies have focused on how seagrass-pathogen interactions are affected by a single external factor, even though coastal areas can be subjected to changes in several environmental factors (e.g., temperature and light) simultaneously (e.g., ). Only two studies have, to our knowledge, specifically aimed at investigating interactive effects between two or several factors (; ). Investigations of salinity, temperature, sulfide, and night-time hypoxia showed that the interaction between ambient salinity and high sulfide concentrations results in larger Labyrinthula sp. lesions on T. testudinum (), whereas temperature and light have no interactive effect on lesion size caused by L. zosterae infection in Z. marina in Washington, United States (). This interaction was not detected in a laboratory experiment with L. zosterae and Z. marina in Germany, instead increased temperature in combination with low salinity reduced L. zosterae infection (). In general, interactive effects of climate change and infectious diseases have not been well studied (Sullivan et al., 2018), and more complex studies are needed to predict future outcomes of climate change on pathogen outbreaks.
The overall aim of this study was to investigate how individual and interactive effects of temperature, light, and tissue damage affect L. zosterae infection in Z. marina in a fully crossed infection experiment. The degree of infection was measured both as lesion coverage and as L. zosterae cell concentration. We also investigated if the treatment factors affect growth and the chemical defense (inhibitory capacity of tissue extracts) of the host, and the gene expression pattern of Z. marina in response to pathogen infection. We hypothesized that low light, increased temperature, and tissue damage would increase L. zosterae infection and reduce host resistance, and that interactions among these factors would have an additive or synergistic effect on L. zosterae infection in Z. marina.
Materials and Methods
Experimental Setup
The individual and interactive effects of temperature, light, tissue damage, and L. zosterae infection on Z. marina were investigated in an outdoor laboratory experiment. A total of 320 Z. marina shoots between 22 and 60 cm in length were haphazardly collected at 1–2 m depth from the Ide Fjord, Sweden, by snorkeling in August 2015 and transported to Tjärnö marine laboratory, Sweden. The sampling was carried out in meadows free from L. zosterae infection (). The shoots were planted into individual 6 L aquariums containing 600 mL sterilized sediment and connected to an outdoor flow-through seawater system (19–31 PSU, X̄ = 23 PSU). Surface seawater was pumped from 1.5 m depth, filtered (10 μm), and distributed to the individual aquariums. Four factors (temperature, light, tissue damage, and infection) with two levels each were manipulated in a fully factorial design (n = 20). Shoots received seawater at either ambient (19–20°C) or increased (25–26°C) temperatures, and light at ambient (approximately 600 μmol m–1 s–1) or reduced (approximately 100 μmol m–1 s–1) levels. The latter temperature and light treatments correspond to the expected maximum surface-seawater conditions in shallow bays on the northwest coast of Sweden during a heat wave (), and the measured irradiance at 5 m depth within Z. marina meadows located in the Tjärnö archipelago (). The incoming seawater for the increased temperature treatment was elevated and maintained by a sensory controlled heating system and thereafter distributed into the separate aquariums. Different light levels were achieved by placing one or several layers of fiber cloth above the aquariums. The natural light/dark cycle at Tjärnö, Sweden at the time of the experiment was approximate 14 h light:10 h dark. Shoots were acclimated to their temperature and light conditions for 1 week before the start of the experiment. Shoots receiving increased temperature were initially placed in seawater with ambient temperature that increased to the intended temperature over the course of 2 days. Treatments and replicates were randomized at the start of the experiment.
At the start of the experiment shoots were either left intact or had their leaves cut by scissors to simulate damages by boat propellers. One clean cut was made over all the leaves at approximately 10 cm from the tip of the longest leaf, resulting in all or most leaves of the plant being cut, but at different individual lengths. The L. zosterae strain for the infection treatment was isolated from Z. marina shoots collected in the Tjärnö archipelago in July 2015. To isolate L. zosterae, 3 cm pieces of Z. marina leaves with apparent lesions were surface sterilized in 0.5% sodium hypochlorite, rinsed in distilled water, soaked in filtered (0.2 μm) seawater, and finally placed upon serum seawater agar (SSA) plates. The SSA plates contained filtered (0.2 μm) seawater (30 PSU), 12 g L–1 agar, 1 g L–1 glucose L–1, 0.1 g L–1 peptone, 0.1 g L–1 yeast extract, and 3 mg L–1 germanium dioxide. The medium was autoclaved after which 10 ml L–1 horse serum and 25 ml L–1 streptomycin/penicillin (10,000 units penicillin and 10 mg streptomycin mL–1) were added. The SSA plates with L. zosterae (Figure 1A) cultures were incubated in the dark at 25°C. Inoculation of Z. marina plants in the experiment was performed by attaching L. zosterae-covered agar plugs, taken from the SSA plates (Ø = 5 mm), to the third youngest leaf with a small piece of silicone tubing. Agar plugs without L. zosterae were attached to non-inoculated shoots. All Z. marina shoots were photographed (Figure 1B) before and after the 7-day experiment. At the end of the experiment, Z. marina tissue samples were taken for measurements of gene expression, and the remaining parts of the shoots were placed in −80°C for further assessment of chemical defense and L. zosterae cell concentration (see paragraphs below).
FIGURE 1
Zostera marina Growth and Lesion Coverage
ImageJ software was used to measure Z. marina growth (% of initial shoot area, n = 20) and L. zosterae disease symptoms, i.e., lesion coverage (% of final shoot area, n = 20) from the photographs taken before and after the experiment.
Labyrinthula zosterae Concentration
The pathogen load was investigated by measuring the L. zosterae cell concentration in seagrass tissue by extracting L. zosterae DNA from the whole shoots (n = 20) at the end of the experiment. Each sample was freeze-dried and homogenized with a mortar and pestle, after which 2–3 mg of each sample was extracted for L. zosterae DNA using the Invisorb Spin Plant Mini Kit (Stratec Molecular, Berlin, Germany) according to the manufacturers protocol, with the exception that 1 μL salmon sperm (Invitrogen, Waltham, MA, United States) was added to each extraction to saturate the silica columns with untargeted DNA. DNA was eluted in 100 μL elution buffer and the samples were analyzed using real-time quantitative PCR (RT-qPCR) performed on a StepOnePlus qPCR machine (Applied Biosystems, Waltham, MA, United States) as previously described by . Each sample ran in technical triplicates, and cycle threshold (CT) was calculated with a fixed threshold of 0.06. A standard curve with known L. zosterae concentrations of 0.75 cells μL–1 (CT: 36.75 ± 1.21 mean ± SD), 7.5 cells μL–1 (CT: 32.87 ± 1.04 mean ± SD), 75 cells μL–1 (CT: 29.23 ± 1.04 mean ± SD), and 750 cells μL–1 (CT: 24.59 ± 0.95 mean ± SD) ran along the samples on each plate. Samples with CT values over 39, or where the standard variation of the triplicates exceeded 0.5, were excluded from further analysis.
Chemical Defense Assay
Eelgrass shoots of each treatment (n = 15) were used to investigate if environmental stressors affect Z. marina’s production of inhibitory compounds against L. zosterae. The shoots were freeze-dried and homogenized, and individually extracted in methanol and dichloromethane 1:1 for 1 h on a shaker table. The samples were diluted to 1/16 volumetric concentration to better detect inhibitory differences between extracts from different shoots, since all Z. marina extracts at natural concentration strongly inhibit L. zosterae growth (). The plant material was filtered (1.6 μm) and solvents were evaporated using a SpeedVac.
The inhibitory effects of the Z. marina extracts on the pathogen was tested in a modified L. zosterae growth assay first described by . Each sample was redissolved in 0.7 mL liquid growth media (SSA without agar, see above) containing 1% dimethyl sulfoxide. The samples were transferred to six well plates (Ø = 35 mm, Thermo Fisher Scientific, Waltham, MA, United States), and L. zosterae covered plugs (Ø = 7 mm) taken from the edge of the L. zosterae cultures growing on SSA plates were placed in the center of each well. L. zosterae used in the growth assays was isolated from Z. marina shoots collected in the Tjärnö archipelago in August 2016 as described in ‘Experimental setup’ above. Wells containing growth media without Z. marina extracts (n = 48) served as controls. The well plates were sealed with Parafilm and incubated at 25°C in the dark. After 11 h the size of the L. zosterae colonies were outlined with a marker and the well plates photographed. The total area of the colonies was measured using ImageJ software.
Gene Expression
A 5 cm piece from the middle part of the fourth leaf of each shoot (n = 6) was placed in RNA-later (Thermo Fisher Scientific, Waltham, MA, United States) at the end of the experiment for gene expression measurements. We chose to analyze the middle section of fourth leaf (adjacent to the originally inoculated leaf) for gene expression, thus avoiding the mostly necrotic tissue (i.e., lesions caused through L. zosterae infection) of the originally inoculated leaf. RNA was extracted with the Invitrap® Spin Plant RNA Mini kit (Stratec Molecular, Berlin, Germany) according to the manufactures protocol. Next, cDNA was synthesized using qScriptTM cDNA synthesis kit (Quanta Biosciences, Beverly, MA, United States) in combination with DNA wipe-out treatment by PerfeCTa® DNase I treatment (Quanta Biosciences, Beverly, MA, United States).
Fourteen target genes that previously have been investigated in connection with reduced light intensity, increased temperature, and L. zosterae infection (; Winters et al., 2011; , ; Salo et al., 2015) were selected for gene expression analysis, along with development of two new primer pairs to target the phenol synthesis. See Table 1 for a summary of the genes and their function, and for further details on the newly developed primers.
TABLE 1
| Gene association | Gene code | Gene name | Function | Source | P Primer sequences |
| General stress, redox and | hsp60a | Heat shock protein 60 isoform 1 | Molecular chaperone | F: 5′-GGTTGATTCTGCAAGCGTCT-3′ R: 5′-ATCCCACCCATTCCTCCA-3′ | |
| detoxification | hsp80 | Heat shock protein 80 | Molecular chaperone | F: 5′-GCAATTCTTCCGATGCTCTT-3′ R: 5′-GCGATCGTACCCAAGTTGTT-3′ | |
| BIP | Luminal binding protein | Luminal binding protein | F: 5′-AATGGTTCTTGCTGGATGTTT-3′ R: 5′-CCCATGATGGTGTGATACGA-3′ | ||
| PPIM1 | 70 kDa peptidyl-prolyl isomerase | Catalytic activity | F: 5′-TATGCACCAAGGTGCTGGA-3′ R: 5′-TTCCACTTTAACATCCCTGTTG-3′ | ||
| CAT | Catalase II | Antioxidant | Winters et al., 2011 | F: 5′-ACAAAATTCCGTCCGTCA-3′ R: 5′-GTCCTCAAGGAGTATTGGTCCTC-3′ | |
| GST | Glutathione S-transferase | Detoxification | Winters et al., 2011 | F: 5′-CATGAATCCATTCGGACAAG-3′ R: 5′-CAGCAAGGTGAGTAAGGTCAG-3′ | |
| Immune defense | RppA | NB-ARC domain-containing disease resistance gene | Immune receptor | F: 5′-GCATCACATCGATATCTGATTCTTT-3′ R: 5′-CTGTGGTAATTTCGACCCATC-3′ | |
| pl 206 | Disease resistance-responsive protein 206 | Pathogenesis-related protein | F: 5′-CTCTTCTAGCACGCAATTTGG-3′ R: 5′-CCGAAAATGTCTCCTTCGAG-3′ | ||
| CLT1 | Chitinase like protein 1 | Pathogenesis-related protein | F: 5′-AAACAGCCATCAGCACATGA-3′ R: 5′-GTCAGCAAATCCCTGTCCAC-3′ | ||
| Phenol synthesis | CYP73A | Trans-cinnamate 4-monooxygenase | Phenol synthesis | F: 5′-ATATCCACCTTGTCCATTCCC-3′ R: 5′-CTGACTTCCGATACTTGCCT-3′ | |
| PAL1a | Phenylalanine ammonia lyase 1a | Phenol synthesis | This study | F: 5′-TTTATCTGCCGTGTTCTGTG-3′ R: 5′-GATGCTCGAATAACTTCCACC-3′ | |
| PAL1b | Phenylalanine ammonia lyase 1b | Phenol synthesis | This study | F: 5′-GATAACCGAGTTCACGAATGACC-3′ R: 5′-CAGACACTTAACCAGCACATCC-3′ | |
| Primary metabolism | Chl_synth | Chlorophyll synthase | Enzyme, chlorophyll synthesis | Salo et al., 2015 | F: 5′-GTCTCAGATGCCCTAACCAC-3′ R: 5′-TCCCTCATAACTACTTCACGAC-3′ |
| RuBisCO | Ribulose-1,5-biphosphate carboxylase/oxygenase | Enzyme, photosynthesis | Salo et al., 2015 | F: 5′-CCATCTCTACCGCTATCCCT-3′ R: 5′-GACGACCTCACAACAAACCT-3′ | |
| PSI_IV | Photosystem I, reaction center subunit IV | Photosynthesis, chloroplast precursor | Winters et al., 2011 | F: 5′-GGGAACCAAGGTGAAGATT-3′ R: 5′-GAATCTCACCACAACTGGGTA-3′ | |
| PSII_22kDa | Photosystem II, 22 kDa protein | Photosynthesis, chloroplast precursor | Winters et al., 2011 | F: 5′-TTCCCAAAAAGGTGGTAGTTA-3′ R: 5′-ATAAAGAAGCGGCAAAACC-3′ | |
| Reference | 10 kDa | 10 kDa chaperonin | Reference gene | F: 5′-CGTCACTCCCAAGTATACTACTGTC-3′ R: 5′-CTCCGGTCGGAACACTGA-3′ | |
| APX | L-Ascorbate peroxidase 2 (cytosolic) | Reference gene | Winters et al., 2011 | F: 5′-GGTGATTTCTACCAGCTTGC-3′ R: 5′-GATCCGCACCTTGGGTA-3′ | |
| SOD | Superoxidase dismutase (mitochondrial) | Reference gene | Winters et al., 2011 | F: 5′-ATGGGTGTGGCTTGCTTA-3′ R: 5′-ATGCATGCTCCCATACATCT-3′ | |
| GRFb | Growth-Regulating Factor | Reference gene | This study | F: 5′-CATCTCCATCCCTAACTCATCAC-3′ R: 5′-CCCGTATGTTAGCAGATTCAGAG-3′ | |
| ubiquitin | Ubiquitin | Reference gene | This study | F: 5′-GGATCATCAAGGAAACGCAG-3′ R: 5′-GCAGCCATAGGATACTCTTCAG-3′ |
Information on target genes.
Gene expression was measured with a Biomark HD System (Fluidigm, South San Francisco, CA, United States) on a 96.96 Dynamic Array IFC chip according to published protocols (Salo et al., 2015). Each primer pair ran in technical triplicates. If the standard deviation of a triplicate was >0.5 the value was excluded from further analysis. To normalize gene expression data five reference genes were chosen based on an algorithm geNorm implemented in the software qbase+ (Biogazelle, Zwijnaarde, Belgium). The results were cross validated by two additional algorithms, Bestkeeper (), and NormFinder (), see Supplementary Table S1. Based on the combined results we used 10kDa, SOD, APX, GRFb, and ubiquitin as reference genes. The cycle threshold (Ct) values of target genes were normalized with the geometric mean of the Ct of reference targets. The relative gene expression levels were calculated as: −ΔCT = Geometric mean of reference genes (CT)–CT (target gene).
Statistical Analyses
Lesion coverage (% of total shoot area) and L. zosterae cell concentration (no. of L. zosterae cells mg–1Z. marina dry weight) were statistically analyzed in SuperANOVA (ABACUS Concepts) statistical software using three-way analysis of variance (ANOVA) with temperature (two levels), light (two levels), and tissue damage (two levels) as fixed, orthogonal factors.
Growth of Z. marina shoots (% growth of initial shoot size) and the inhibitory effect of the extracted Z. marina metabolites on L. zosterae growth (% growth of control) in the bioassays were analyzed using four-way ANOVA with temperature (two levels), light (two levels), tissue damage (two levels), and infection (two levels) as fixed, orthogonal factors.
Relative gene expression values (−ΔCT) were statistically analyzed using permutational multivariate ANOVA (Permanova) for each functional gene group (general stress genes and redox and detoxification genes, genes associated to immune response, genes of the phenol synthesis, and genes of primary metabolism), with infection (two levels), temperature (two levels), light (two levels), and tissue damage (two levels) as fixed orthogonal factors. Permanova was conducted, using Euclidian distances and 999 permutations of the package vegan (Oksanen et al., 2016), within R version 3.1.2. In order to illustrate the results of the gene expression for the different treatments a heat map was created using the mean −ΔCT values. Principal component analyses (PCA) were performed for each gene group with significant treatment effects in the Permanova, using the ade4 package () in R (version 3.3.1). The expression of individual genes in comparison to undamaged, non-inoculated shoots at ambient conditions of light and temperature (hereafter referred to as ambient control) is reported as relative gene expression, −ΔΔCT = −ΔCT−(−ΔCT (ambient control)).
Prior to all statistical analyses, data was tested for homogeneity of variances with a Cochran’s test, and transformed when required. Data on L. zosterae growth in the chemical defense assays were square root transformed prior to statistical analysis. However, the order of the mean values changed after transformation for data on Z. marina growth. Due to the large number of groups and replicates, ANOVA is robust against heterogeneous variances (Underwood, 1997) and therefore, the data was not transformed. Multiple comparisons of means were performed using the Student–Newman–Keuls (SNK) procedure.
Results
Labyrinthula zosterae Lesion Coverage and Cell Concentration
All Z. marina shoots developed lesions when inoculated with L. zosterae cells, whereas no lesions were detected on non-inoculated shoots. After the 7 days experiment lesions covered large parts of the inoculated leaf, but also spread to a smaller extent to other leaves within the shoots. There was no significant interaction between the treatment factors (Table 2). However, reduced light conditions resulted in a significant increase of lesion coverage (Table 2), with shoots in low light treatments showing 35.0% ± 0.7 SE (mean ± SE) higher lesion coverage compared to shoots under ambient light (Figure 2). Increased temperature and tissue damage had no significant effect on lesion coverage (Table 2).
TABLE 2
| Lesion coverage | L. zosterae cell concentration | ||||||
| Source of Variation | df | MS | F | P | MS | F | P |
| Temperature (T) | 1 | 0.67 | 0.01 | 0.95 | 1144.56 | 0.52 | 0.47 |
| Light (L) | 1 | 2389.21 | 16.62 | <0.01 | 6093.89 | 2.76 | 0.10 |
| Damage (D) | 1 | 17.64 | 0.01 | 0.73 | 11629.22 | 5.26 | 0.02 |
| T × L | 1 | 26.33 | 0.18 | 0.67 | 17945.52 | 8.12 | <0.01 |
| T × D | 1 | 183.08 | 1.27 | 0.26 | 0.39 | <0.01 | 0.99 |
| L × D | 1 | 236.53 | 1.65 | 0.20 | 3024.51 | 1.37 | 0.24 |
| T × L × D | 1 | 0.10 | <0.01 | 0.98 | 604.37 | 0.27 | 0.60 |
| Residual | 152 | 143.76 | 2209.73 | ||||
Analysis of variance of lesion coverage (%) and Labyrinthula zosterae cell concentration (cells mg–1Zostera marina dry weight) in Z. marina grown under different levels and combinations of temperature (ambient and increased), light (reduced and ambient), and tissue damage (damaged and undamaged).
FIGURE 2
Reduced light conditions significantly (SNK, P < 0.05) increased the L. zosterae cell concentration in Z. marina, but only under ambient temperatures, as shown by the significant interaction between the main factors light and temperature (Table 2). Shoots experiencing ambient temperature showed a slightly higher concentration of L. zosterae cells following infection under reduced light conditions compared to shoots kept at ambient light (Figure 3A). At increased temperatures, there was no significant (SNK, P > 0.05) difference in L. zosterae cell concentration between shoots exposed to different light levels. Contrary to our hypothesis, tissue damage had a significant negative effect on L. zosterae cell concentration in Z. marina leaves (Table 2 and Figure 3B).
FIGURE 3
Zostera marina Growth
Eelgrass growth was observed under all treatment conditions and varied between 8 and 33% of the initial shoot area. When data on Z. marina growth was analyzed, we found a significant interaction between the factors damage and temperature (Table 3). Damaged Z. marina shoots grew significantly (SNK, P < 0.05) more (58.9 ± 2.1% mean ± SE) compared to undamaged shoots under ambient temperature, whereas only a trend toward a similar result was found at increased temperature (Figure 4A). Furthermore, we found a significant interaction between the factors light, temperature, and infection on Z. marina growth (Table 3). Highest growth was observed for non-inoculated Z. marina shoots exposed to increased temperature under ambient light conditions (Figure 4B). However, there was no significant (SNK, P > 0.05) difference in growth between inoculated and non-inoculated Z. marina shoots exposed to different temperatures under reduced light conditions (Figure 4B).
TABLE 3
| Source of variation | df | MS | F | P |
| Temperature (T) | 1 | 206.07 | 1.01 | 0.32 |
| Light (L) | 1 | 20777.05 | 102.01 | <0.01 |
| Damage (D) | 1 | 3216.83 | 15.79 | <0.01 |
| Infection (I) | 1 | 1256.03 | 6.17 | 0.01 |
| T × L | 1 | 87.83 | 0.43 | 0.51 |
| T × D | 1 | 810.16 | 3.98 | 0.05* |
| T × I | 1 | 237.74 | 1.17 | 0.28 |
| L × D | 1 | 64.39 | 0.32 | 0.58 |
| L × I | 1 | 1.93 | 0.01 | 0.92 |
| D × I | 1 | 4.33 | 0.02 | 0.88 |
| T × L × D | 1 | 31.08 | 0.15 | 0.70 |
| T × L × I | 1 | 853.20 | 4.19 | 0.04 |
| T × D × I | 1 | 109.60 | 0.54 | 0.46 |
| L × I × D | 1 | 224.72 | 1.10 | 0.29 |
| T × L × I × D | 1 | 178.05 | 0.87 | 0.35 |
| Residual | 304 | 203.68 |
Analysis of variance of Zostera marina growth (% growth relative to initial size) under different levels and combinations of temperature (ambient and increased), light (reduced and ambient), tissue damage (damaged and undamaged), and Labyrinthula zosterae infection (non-inoculated and inoculated).
*P = 0.047. Bold indicate significance (P < 0.05). Data on mean values and standard error are presented in Figure 4.
FIGURE 4
Chemical Defense Assay
The statistical analysis of data on the inhibitory effect of Z. marina extracts on L. zosterae growth revealed a significant interaction between the factors tissue damage and temperature (Table 4). The inhibitory effect on L. zosterae growth was significantly (SNK, P < 0.05) higher for extracts from undamaged Z. marina shoots exposed to ambient temperature compared to undamaged shoots at increased temperature or damaged shoots at ambient temperature (Figure 5).
TABLE 4
| Source of variation | df | MS | F | P |
| Temperature (T) | 1 | 13.69 | 5.20 | 0.02 |
| Light (L) | 1 | 1.23 | 4.68 | 0.50 |
| Damage (D) | 1 | 0.24 | 0.09 | 0.76 |
| Infection (I) | 1 | 1.93 | 0.73 | 0.39 |
| T × L | 1 | 8.64 | 3.28 | 0.07 |
| T × D | 1 | 27.43 | 10.42 | <0.01 |
| T × I | 1 | 2.97 | 1.13 | 0.29 |
| L × D | 1 | 8.00 | 3.04 | 0.08 |
| L × I | 1 | 6.40 | 2.43 | 0.12 |
| D × I | 1 | 0.37 | 0.14 | 0.71 |
| T × L × D | 1 | 4.85 | 1.84 | 0.18 |
| T × L × I | 1 | 1.44 | 0.55 | 0.46 |
| T × D × I | 1 | 0.61 | 0.23 | 0.63 |
| L × I × D | 1 | 0.69 | 0.26 | 0.61 |
| T × L × I × D | 1 | 3.30 | 1.25 | 0.27 |
| Residual | 224 | 2.63 |
Analysis of variance of inhibition of Labyrinthula zosterae growth by extracts of Zostera marina (% relative to controls) from shoots grown under different levels and combinations of temperature (ambient and increased), light (reduced and ambient), tissue damage (damaged and undamaged), and L. zosterae infection (non-inoculated and inoculated).
Data is square root transformed to meet assumptions for ANOVA. Bold indicate significance (P < 0.05). Data on mean values and standard error are presented in Figure 5.
FIGURE 5
Gene Expression
Expression of the 16 targeted Z. marina genes following L. zosterae infection under different combinations of light, temperature and tissue damage showed a complex pattern with multiple interactive effects between treatments (Supplementary Figure S1). A summary of the expression of the different functional gene groups in relation to non-inoculated shoots at ambient conditions (ambient control) is shown for the groups with significant differences (PERMANOVA, P < 0.05) among treatments (Table 5).
TABLE 5
| Treatment | Genes associated with | ||||
| Light | Temperature | Infection | Immune defense | Phenol synthesis | Primary metabolism |
| Ambient | Ambient | Non-inoculated | → | → | → |
| Inoculated | ↓ | → | → | ||
| Increased | Non-inoculated | → | ↓ | → | |
| Inoculated | → | ↓ | → | ||
| Reduced | Ambient | Non-inoculated | ↓ | ↓ | ↓ |
| Inoculated | ↓ | ↓ | → | ||
| Increased | Non-inoculated | → | → | → | |
| Inoculated | ↓ | → | ↓ | ||
Expression of genes associated with immune defense, phenol synthesis, and primary metabolism in Zostera marina grown under different levels and combinations of light, temperature and Labyrinthula zosterae infection.
Gene expression is based on results from PERMANOVA analyses, and expressed relative to the response of the genes in non-inoculated Z. marina shoots at ambient light and temperature (top row, i.e., ambient control). → = ambient gene expression, ↓ = downregulated gene expression.
The expression pattern of targeted immune genes was affected by a statistically significant interactive effect between temperature, light and infection (Table 6). The PCA showed that 87.7% of the total variation is explained by PC1 (53.5%) and PC2 (34.2%). PC1 is associated with the genes CLT1 and RppA and showed that inoculated shoots had a decreased expression of these genes compared to non-inoculated shoots, in both the treatment combinations ambient temperature/ambient light and increased temperature/reduced light, but not in ambient temperature/reduced light, or increased temperature/ambient light (Table 6). However, ambient temperature/reduced light led to a reduced expression of the targeted immune genes compared to non-infected shoots at ambient temperature and light (Table 6).
TABLE 6
| Stress, redox and detox responses | Immune defense | Phenol synthesis | Primary metabolism | ||||||||||
| Variable | df | MS | F model | Pr (>F) | MS | F model | Pr (>F) | MS | F model | Pr (>F) | MS | F model | Pr (>F) |
| Infection (I) | 1 | 16.36 | 1.05 | 0.34 | 11.36 | 2.48 | 0.07 | 41.22 | 2.38 | 0.09 | 14.93 | 1.74 | 0.15 |
| Temperature (T) | 1 | 26.33 | 1.68 | 0.34 | 14.71 | 3.21 | 0.03 | 11.27 | 0.65 | 0.51 | 13.09 | 1.52 | 0.21 |
| Light (L) | 1 | 2.78 | 0.18 | 0.94 | 26.56 | 5.80 | 0.03 | 41.64 | 2.40 | 0.09 | 96.51 | 11.22 | <0.01 |
| Damage (D) | 1 | 13.54 | 0.87 | 0.43 | 1.17 | 0.26 | 0.83 | 20.89 | 1.21 | 0.28 | 13.90 | 1.62 | 0.18 |
| T × I | 1 | 10.10 | 0.65 | 0.56 | 2.62 | 0.57 | 0.61 | 19.72 | 1.14 | 0.30 | 5.24 | 0.61 | 0.57 |
| L × I | 1 | 9.26 | 0.59 | 0.60 | 5.51 | 1.20 | 0.31 | 0.62 | 0.04 | 0.99 | 9.88 | 1.15 | 0.31 |
| T × L | 1 | 18.75 | 1.20 | 0.27 | 10.02 | 2.19 | 0.10 | 77.22 | 4.46 | 0.02 | 10.01 | 1.16 | 0.31 |
| I × D | 1 | 8.10 | 0.52 | 0.66 | 8.07 | 1.76 | 0.16 | 10.25 | 1.59 | 0.53 | 0.46 | 0.05 | 0.99 |
| T × D | 1 | 9.81 | 0.63 | 0.57 | 0.75 | 1.16 | 0.91 | 7.67 | 0.44 | 0.64 | 12.08 | 1.41 | 0.23 |
| L × D | 1 | 10.20 | 0.65 | 0.56 | 6.05 | 1.32 | 0.26 | 6.39 | 0.37 | 0.70 | 6.78 | 0.79 | 0.46 |
| T × L × I | 1 | 435.15 | 2.25 | 0.09 | 19.06 | 4.16 | 0.01 | 41.14 | 2.37 | 0.10 | 27.62 | 3.21 | 0.04 |
| T × D × I | 1 | 14.34 | 0.92 | 0.40 | 2.84 | 0.62 | 0.56 | 7.52 | 0.43 | 0.63 | 12.03 | 1.40 | 0.23 |
| L × I × D | 1 | 19.12 | 1.22 | 0.27 | 1.04 | 0.24 | 0.88 | 9.08 | 0.52 | 0.57 | 3.01 | 0.35 | 0.75 |
| T × L × D | 1 | 3.15 | 0.20 | 0.92 | 5.66 | 1.24 | 0.29 | 7.57 | 0.44 | 0.66 | 3.32 | 0.39 | 0.73 |
| T × L × I × D | 1 | 9.66 | 0.62 | 0.58 | 1.63 | 0.36 | 0.75 | 1.73 | 0.10 | 0.94 | 4.63 | 0.54 | 0.61 |
| Residual | 72 | 15.65 | 4.58 | 17.33 | 8.60 | ||||||||
Permutational analysis of variance of expression of genes associated with different functions (stress, redox, and detox responses, immune defense, phenol synthesis, and primary metabolism) in Zostera marina in response to different levels and combinations of temperature (ambient and increased), light (reduced and ambient), tissue damage (damaged and undamaged), and Labyrinthula zosterae infection (non-inoculated and inoculated).
Bold indicate significance (P < 0.05).
The expression of genes associated with phenol synthesis was influenced by a significant interactive effect between light and temperature (Table 6). The first two principal components explained 88.6% of the total variation, and accounted for 72.6 and 16.0% of the variation, respectively. PC1 is associated with CYP73A, PAL1a, and PAL1b and showed a decreased expression of these genes in Z. marina when the shoots were subjected to combinations of ambient light/increased temperature, and reduced light/ambient temperature (Table 6).
As for the targeted immune genes, the expression of genes associated with phenol synthesis was influenced by a significant interactive effect between temperature, light and infection (Table 6). PCA showed that 88.3% of the total variation was explained by PC1 (49.3%) and PC2 (39.0%). There was a down regulation of the genes for inoculated shoots, but only at increased temperature/reduced light. In addition, non-inoculated shoots at ambient temperature/reduced light also had a decreased expression of the genes associated with the primary metabolism, compared to non-inoculated shoots at ambient conditions of temperature and light (Table 6 and Figure 6).
FIGURE 6
No significant effect between the treatments was detected for the expression on the general stress, and redox and detoxification genes (Table 6). Mean fold-change expression relative to ambient control (−ΔΔCT) for all individual genes are presented in Figure 6.
Discussion
This study investigated how changes in different extrinsic factors influence L. zosterae infection in Z. marina, and is one of the first studies to include interactive effects of multiple factors on Labyrinthula sp. infection in seagrass systems (see review by Sullivan et al., 2018). We show that changes in extrinsic factors can have an impact on the host-pathogen interaction on a relatively short time scale. In general, light had the largest effect on L. zosterae infection, as well as on Z. marina growth and gene expression, but there were also a large number of interactions between light and the other tested factors (temperature, tissue damage, and infection), highlighting the complexity of this host-pathogen system (Figure 7 summarizes the impact of L. zosterae infection on Z. marina shoots for the different treatment combinations).
FIGURE 7
Contrary to our hypothesis, reduced light, increased temperature, and tissue damage did not have additive or synergistic negative effects on L. zosterae infection in Z. marina. This corroborate the results from previous studies investigating the effect of multiple factors on Labyrinthula sp. infection in seagrass (Z. marina and T. testudinum), where the effect of single factors did not differ from the effect of cumulative stressors (; ). Instead, light intensity had the strongest impact on disease lesions in the current study, and low light significantly increased lesion coverage to a moderate degree. The effects of light intensity on L. zosterae infection found in this study support previous laboratory findings, showing that reduced light increases lesion coverage (). In contrast, laboratory experiments in the western Baltic Sea found no effect of light on L. zosterae infection (). Although our study and the study by are similar in many aspects (model and pathogen species used, length of experiment, and experimental year), different light sources were used for the different experiments () used artificial light, whereas our study used natural light). This could be a possible explanation to why the results differ. It is also important to consider potential differences in virulence between the pathogens used in respective experiment, since infection by more virulent L. zosterae strains can result in 25% larger lesion coverage on Z. marina plants, compared to infection by less virulent strains (). As of yet, little is know about strain diversity of L. zosterae in northern Europe, and it would be beneficial to our understanding of this host-pathogen interaction to include different stains as a factor in future studies.
Whether light intensity is a primary driver for disease severity in natural seagrass populations is not clear. Increased lesion coverage can be detected on shoots in T. testudinum populations at lower light intensities (Trevathan-Tackett et al., 2013), whereas shallow-growing Z. marina shoots, exposed to relatively high light levels, have larger lesion coverage compared to deeper growing shoots receiving less light (). Although the present study shows that light can regulate L. zosterae infection in Z. marina, the results further suggest that other factors, besides light level, significantly influences the different degree of L. zosterae infection between water depths in natural Z. marina populations.
Reduced light also significantly increased pathogen cell concentrations, but only at ambient temperature. Thus, at higher temperatures, lesion coverage did not correlate with L. zosterae cell concentration. Furthermore, at ambient temperatures, the effect of reduced light was approximately four times greater on lesion coverage than on pathogen load, showing that a relatively small increase in cell concentration can have a large impact on the degree of disease symptoms in this system, or that Z. marina is more sensitive to L. zosterae cells in a light limiting environment.
Increased temperature did not increase L. zosterae infection in Z. marina in the present study, and the reported effects of temperature on Labyrinthula spp. infection in other seagrass systems are variable. A temperature increase from 11 to 18°C did not alter lesion size on Z. marina leaves following L. zosterae infection in laboratory experiments (), whereas increased temperature from 22 to 27°C decreases infection, but only in combination with low salinity (12 PSU) (). Increased temperatures also reduced the size of lesions caused by Labyrinthula spp. in laboratory experiments with the two Mediterranean seagrass species P. oceanica and C. nodosa (; ). While a negative correlation between lesion size and temperature was found between 24 and 32°C for C. nodosa, P. oceanica had an upper threshold for infection with a constantly low lesion size at temperatures above 28°C (). In the subtropical seagrass T. testudinum laboratory investigations showed that a temperature increase from 27 to 30°C had no effect on lesion size caused by Labyrinthula sp., whereas exposure to extreme heat (45°C) reduced the lesion size (). Optimum growth temperature varies among different Labyrinthula species (), which could explain the different effects of temperature on the Labyrinthula spp. infection reported in the literature. However, no study has reported an enhancing effect of increased temperature on lesion size. Thus, in consistence with previous investigations (; ; ; ; ), our results suggest that increased temperature is not a sole driver of increased Labyrinthula spp. virulence and wasting disease outbreaks in seagrass.
Lesion coverage was reduced in Z. marina shoots with tissue damage in this study. In contrast, field experiments with T. testudinum showed that infection by Labyrinthula sp. was more prevalent in shoots subjected to tissue damage compared to undamaged shoots (). The seagrass leaves in the study by where cut very close to the rhizome of the seagrass shoots, whereas in our experiments only moderate amounts of leaf tissue (similar to damages by boat propellers) were removed from the Z. marina shoots. The level of tissue removal in the present study potentially enabled Z. marina to manage the L. zosterae infection better than what was possible for T. testudinum in the study by . Our 7-day experiment suggests that moderate tissue loss, mimicking damage by small boat propellers on the Swedish west coast, will not impact Z. marina negatively. Whether repeated tissue damage over a longer time period would affect the interaction between Z. marina and L. zosterae differently warrants further investigation.
Infection by L. zosterae did not increase the production of inhibitory compounds in Z. marina, in accordance with findings from previous field and experimental investigations (; ), where no correlation between pathogen pressure and inhibition strength of extracted metabolites was found. The chemical identity of the inhibitory compounds is presently not known, but phenolic acids, which increase in seagrasses in response to L. zosterae infection in laboratory inoculation experiments (; Trevathan-Tackett et al., 2013), have been suggested to function as chemical defenses against Labyrinthula spp. (e.g., ; Vergeer et al., 1995). However, in contrast to previous studies (, ), infection did not upregulate genes associated with the phenol synthesis in the present study, suggesting that phenolic acids are produced at similar levels in shoots independent of infection status. The expression of trans-cinnamate 4-monooxygenase (CYP73A) as well as phenolic acid concentrations have previously been investigated in the primary leaf of inoculation (Vergeer and Develi, 1997; Steele et al., 2005; , ; ), indicating that upregulation of phenolic acids may be localized to areas of infection. In this study, gene expression was measured in a secondary leaf, which could explain why no upregulation of CYP73A was detected in the Z. marina shoots.
In contrast to results from previous studies on phenolic acids (Vergeer et al., 1995; ; ), increased temperature and reduced light intensity did not decrease Z. marina’s production of inhibitory compounds against L. zosterae. Since we did not find an upregulation of the genes associated with the phenol synthesis, it is likely that compound(s) other than phenolic acids are responsible for the inhibitory effect against L. zosterae found in this study. Little is known about how a changing environment is affecting other secondary metabolites in seagrasses (Zidorn, 2016). A review on how factors connected to global change can affect chemical defenses against insect herbivores in terrestrial plants shows that the outcome is species specific, and will depend on the compound group investigated (). Therefore, identification of the inhibitory compounds found in this study is needed to further investigate the effects of environmental change on chemically mediated interactions between L. zosterae and Z. marina.
Not surprisingly, Z. marina shoots kept under reduced light conditions had significantly decreased growth compared to shoots subjected to ambient light. Furthermore, increased temperature significantly enhanced growth in uninfected plants under ambient light conditions, but this effect disappeared when shoots were inoculated with L. zosterae. These results indicate that under ambient light conditions, increased temperature has a positive, and L. zosterae infection a negative, effect on Z. marina in terms of growth. Furthermore, Z. marina growth increased significantly in response to tissue damage predominantly at ambient temperatures, indicating that Z. marina shows compensatory growth responses when exposed to tissue damage. Previous studies show, similar to the results in this study, that seagrasses can compensate loss of leaf area by increasing their growth (; Valentine et al., 1997; ; Vergés et al., 2008). However, the response to tissue damage depends on the extent of removed tissue, and compensatory growth is not possible for seagrass with extensive physical damage (Vergés et al., 2008).
Gene expression analysis of Z. marina shoots at increased temperature, where no correlation between L. zosterae cell concentration and lesion coverage was detected, shows that shoots with a low degree of disease symptoms at increased temperature/ambient light had a higher expression of the targeted genes associated with immune defense. Whether a higher expression of these genes conveys a higher (or lower) resistance against L. zosterae is not known (). However, these genes encode for pathogenesis related proteins and an immune receptor, and their function in the plant immune response are well described in terrestrial plant pathosystems (van Loon et al., 2006; Spoel and Dong, 2012). Therefore, it is likely that an upregulation of these genes increases the resistance of the plants against L. zosterae, or against secondary infections of other pathogens.
Eelgrass shoots with a higher degree of lesion coverage at increased temperature/reduced light had a downregulation of genes associated with primary metabolism. Heat stress has in previous mesocosm experiments decreased gene expression of PSII_22kDa in Z. marina, although only for one of three populations investigated (Winters et al., 2011). These primary metabolism genes encode essential proteins active in photosynthesis and assimilation of glucose. Most of the proteins have a regular turnover to keep photosynthesis functional, and e. g. the targeted genes PSI_IV (Zosma35g00710) and PSII_22kDa (Zosma60g00630) have been shown to have one of the highest expression within the transcriptome in Z. marina leaves (position 32 and 9, respectively, from 22732 genes in total; ). Thus, if the expression of these genes is decreased, we can assume that this corresponds to a decreased photosynthetic capacity.
As previously detected, expression of genes associated with immune defense is generally lower in Z. marina following L. zosterae infection (, ), although this study shows that interactive effects between L. zosterae inoculation and environmental factors can alter this relationship. This was also observed for genes associated with primary metabolism, where L. zosterae inoculation changes the expression, but only at reduced light conditions. Similar responses have been detected in Arabidopsis thaliana, where heat stress and draught can alter gene expression in virus-treated plants (Prasch and Sonnewald, 2013). It is, however, important to consider that gene expression responses toward microorganisms can be highly localized or systemic, affecting wide parts of the plant (). This study analyzed Z. marina leaves that were not initially inoculated, since inoculated tissue turned largely necrotic after 7 days as a response to the infection and were considered a confounding factor. Instead of sampling directly inoculated tissue we decided to look at a systemic responses in adjacent leaves to the originally inoculated leaves. Significant changes in gene expression as part of a systemic response were detected showing that the gene response toward L. zosterae infection is depending on prevailing environmental conditions. To understand gene expression responses directly in L. zosterae infected tissue and its interaction with stressors will need further studies that carefully control L. zosterae presence in non-necrotic tissue.
Although complex studies with multiple factors and interactive effects are challenging to perform and interpret, our understanding of host-pathogen interactions in the light of a changing environment is becoming increasingly important. This study shows that changing extrinsic factors, in particular light availability, can have a significant impact on the interaction between the seagrass Z. marina and the wasting disease pathogen L. zosterae on a relatively short time scale. In contrast to our experimental study, natural Z. marina meadows exhibit larger and more prevalent lesions where the light intensity is higher (), suggesting the presence of alternate uninvestigated drivers of wasting disease in the field. Along with the large number of interactive effects between the tested factors in this study (light, temperature, tissue damage, and infection) this highlights the complexity of this system, and show that more multi-factorial investigations of abiotic and biotic (e.g., different L. zosterae strains) factors are needed to better understand the impact of extrinsic factors on seagrass-Labyrinthula interactions.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: Figshare (10.6084/m9.figshare.12931247).
Author contributions
SJ-T, HP, and GT designed the study. SJ-T carried out the fieldwork. SJ-T and GT performed experiments. SJ-T, JB, and GT handled sample and data analysis. SJ-T led the writing of the manuscript, with critical input from JB, GT, and HP. All authors contributed to the article and approved the submitted version.
Funding
This work was funded by the Swedish Research Council FORMAS (no. 2011-1193) and additional support was provided by the Linnaeus Centre for Marine Evolutionary Biology (CeMEB) and the Centre for Marine Chemical Ecology (CeMaCE) at the University of Gothenburg, and by Rådman och Fru Ernst Collianders Stiftelse. JB was funded by the Deutsche Bundesstiftung Umwelt (DBU).
Acknowledgments
Thanks to Ylva Durland, Kristie Rigby, and Arielle Zoellin Jilliot for experimental assistance. Parts of the content of this manuscript have been printed in a thesis ().
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: https://www.frontiersin.org/articles/10.3389/fmars.2020.575183/full#supplementary-material
References
1
AndersenC. L.JensenJ. L.ØrntoftT. F. (2004). Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets.Cancer Res.645245–5250. 10.1158/0008-5472.can-04-0496
2
AndersonP. K.CunninghamA. A.PatelN. G.MoralesF. J.EpsteinP. R.DaszakP. (2004). Emerging infectious diseases of plants: pathogen pollution, climate change and agrotechnology drivers.Trends Ecol. Evol.19535–544. 10.1016/j.tree.2004.07.021
3
ArmigerL. C. (1964). An occurrence of Labyrinthula in New Zealand Zostera.N. Z. J. Bot.23–9.
4
BarbierE. B.HackerS. D.KennedyC.KochE. W.StierA. C.SillimanB. R. (2011). The value of estuarine and coastal ecosystem services.Ecol. Monogr.81169–193.
5
BennettR. N.WallsgroveR. M. (1994). Secondary metabolites in plant defence mechanisms.New Phytol.127617–633. 10.1111/j.1469-8137.1994.tb02968.x
6
BergmannN.WintersG.RauchG.EizaguirreC. (2010). Population-specificity of heat stress gene induction in northern and southern eelgrass Zostera marina populations under simulated global warming.Mol. Ecol.192870–2883. 10.1111/j.1365-294x.2010.04731.x
7
Bidart-BouzatM. G.Imeh-NathanielA. (2008). Global change effects on plant chemical defenses against insect herbivores.J. Integr. Plant Biol.501339–1354. 10.1111/j.1744-7909.2008.00751.x
8
BishopN.MartinD. L.RossC. (2017). Effects of multi-stress exposure on the infection dynamics of a Labyrinthula sp.-turtle grass pathosystem.Mar. Ecol. Prog. Ser.581119–133. 10.3354/meps12318
9
BishopN. D. (2013). The effects of Multiple Abiotic Stressors on the Susceptibility of the Seagrass Thalassia Testudinum to Labyrinthula sp., the Causative Agent of Wasting Disease.Master thesis, University of North Florida, Jacksonville, FL.
10
BockelmannA.-C.TamsV.PloogJ.SchubertP. R.ReuschT. B. H. (2013). Quantitative PCR reveals strong spatial and temporal variation of the wasting disease pathogen, Labyrinthula zosterae in northern European eelgrass (Zostera marina) beds.PLoS One8:e62169. 10.1371/journal.pone.0062169
11
BowlesJ. W.BellS. S. (2004). Simulated herbivory and the dynamics of disease in Thalassia testudinum.Mar. Ecol. Prog. Ser.283127–132. 10.3354/meps283127
12
BrakelJ.Jakobsson-ThorS.BockelmannA.-C.ReuschT. B. H. (2019). Modulation of the eelgrass – labyrinthula zosterae interaction under predicted ocean warming, salinity change and light limitation.Front. Mar. Sci.6:268. 10.3389/fmars.2019.00268
13
BrakelJ.ReuschT. B. H.BockelmannA.-C. (2017). Moderate virulence caused by the protist Labyrinthula zosterae in ecosystem foundation species Zostera marina under nutrient limitation.Mar. Ecol. Prog. Ser.57197–108. 10.3354/meps12104
14
BrakelJ.WernerF. J.TamsV.ReuschT. B. H.BockelmannA.-C. (2014). Current European Labyrinthula zosterae are not virulent and modulate seagrass (Zostera marina) defense gene expression.PLoS One9:e92448. 10.1371/journal.pone.0092448
15
BuchsbaumR. N.ShortF. T.CheneyD. P. (1990). Phenolic-nitrogen interactions in eelgrass, Zostera marina L.: possible implications for disease resistance.Aquat. Bot.37291–297. 10.1016/0304-3770(90)90075-v
16
BullJ. C.KenyonE. J.CookK. J. (2012). Wasting disease regulates long-term population dynamics in a threatened seagrass.Oecologia169135–142. 10.1007/s00442-011-2187-6
17
ChakrabortyS.TiedemannA. V.TengP. S. (2000). Climate change: potential impact on plant diseases.Environ. Pollut.108317–326. 10.1016/s0269-7491(99)00210-9
18
CostanzaR.de GrootR.SuttonP.van der PloegS. (2014). Changes in the global value of ecosystem services.Glob. Environ. Change26152–158.
19
DanglJ. L.HorvathD. M.StaskawiczB. J. (2013). Pivoting the plant immune system from dissection to deployment.Science341746–751. 10.1126/science.1236011
20
DawkinsP.EisenlordM.YoshiokaR.FiorenzaE. (2018). Environment, dosage, and pathogen isolate moderate virulence in eelgrass wasting disease.Dis. Aquat. Organ.13051–63. 10.3354/dao03263
21
DrayS.DufourA.-B. (2007). The ade4 Package: implementing the duality diagram for ecologists.J. Stat. Softw.221–20.
22
DuarteC. M.MiddelburgJ. J.CaracoN. (2005). Major role of marine vegetation on the oceanic carbon cycle.Biogeosciences21–8. 10.5194/bg-2-1-2005
23
ErianderL. (2016). Restoration and Management of Eelgrass (Zostera marina) on the West Coast of Sweden.Doctorial thesis, University of Gothenburg, Gothenburg.
24
GiesenW. B. J. T.Van KatwijkM. M.Den HartogC. (1990). Temperature, salinity, insolation and wasting disease of eelgrass (Zostera marina L.) in the Dutch Wadden Sea in the 1930’s.Netherlands J. Sea Res.25395–404. 10.1016/0077-7579(90)90047-k
25
GronerM. L.BurgeC. A.CouchC. S.KimC. J. S. (2014). Host demography influences the prevalence and severity of eelgrass wasting disease.Dis. Aquat. Organ.108165–175. 10.3354/dao02709
26
GronerM. L.BurgeC. A.KimC. J. S.ReesE. (2016). Plant characteristics associated with widespread variation in eelgrass wasting disease.Dis. Aquat. Organ.118159–168. 10.3354/dao02962
27
HayM. E.FenicalW. (1988). Marine plant-herbivore interactions: the ecology of chemical defense.Annu. Rev. Ecol. Syst.19111–145. 10.1146/annurev.es.19.110188.000551
28
HernánG.OrtegaM. J.GándaraA. M.CastejónI.TerradosJ.TomasF. (2017). Future warmer seas: increased stress and susceptibility to grazing in seedlings of a marine habitat-forming species.Glob. Change Biol.234530–4543. 10.1111/gcb.13768
29
InfantesE.ErianderL.MoksnesP.-O. (2016). Eelgrass (Zostera marina) restoration on the west coast of Sweden using seeds.Mar. Ecol. Prog. Ser.54631–45. 10.3354/meps11615
30
Jakobsson-ThorS. (2017). Seagrass Wasting Disease Impact of Abiotic Factors and Chemical Defense.PhD dissertation, University of Gothenburg, Gothenburg.
31
Jakobsson-ThorS.TothG. B.BrakelJ.BockelmannA.-C.PaviaH. (2018). Seagrass wasting disease varies with salinity and depth in natural Zostera marina populations.Mar. Ecol. Prog. Ser.587105–115. 10.3354/meps12406
32
LevinD. A. (1976). The chemical defenses of plants to pathogens and herbivores.Annu. Rev. Ecol. Syst.7121–159. 10.1146/annurev.es.07.110176.001005
33
MaleckK.DietrichR. A. (1999). Defense on multiple fronts: how do plants cope with diverse enemies?Trends Plant Sci.4215–219. 10.1016/s1360-1385(99)01415-6
34
MartinD. L.BooneE.CaldwellM. M.MajorK. M.BoettcherA. A. (2009). Liquid culture and growth quantification of the seagrass pathogen. Labyrinthula spp.Mycologia101632–635. 10.3852/08-171
35
MartinD. L.ChiariY.BooneE.ShermanT. D. (2016). Functional, phylogenetic and host-geographic Signatures of Labyrinthula spp. provide for putative species delimitation and a global-scale view of seagrass wasting disease.Estuar. Coast.391403–1421. 10.1007/s12237-016-0087-z
36
McKoneK. L.TannerC. E. (2009). Role of salinity in the susceptibility of eelgrass Zostera marina to the wasting disease pathogen Labyrinthula zosterae.Mar. Ecol. Prog. Ser.377123–130. 10.3354/meps07860
37
MitchellC. A. (1987). Growth of Halodule wrightii in culture and the effects of cropping, light, salinity and atrazine.Aquat. Bot.2825–37. 10.1016/0304-3770(87)90054-4
38
MoranK. L.BjorndalK. A. (2005). Simulated green turtle grazing affects structure and productivity of seagrass pastures.Mar. Ecol. Prog. Ser.305235–247. 10.3354/meps305235
39
MuehlsteinL. K.PorterD.ShortF. T. (1988). Labyrinthula sp., a marine slime mold producing the symptoms of wasting disease in eelgrass. Zostera marina.Mar. Biol.99465–472. 10.1007/bf00392553
40
OksanenJ.BlanchetF. G.KindtR.LegendreP.MinchinP. R.O’HaraR. B.et al (2016). Vegan: Community Ecology Package. R Package. 2.3–3. Available online at: https://CRAN.R-project.org/package=vegan
41
OlsenJ. L.RouzéP.VerhelstB.LinY.-C. (2016). The genome of the seagrass Zostera marina reveals angiosperm adaptation to the sea.Nature530331–335.
42
OlsenY. S.DuarteC. M. (2015). Combined effect of warming and infection by Labyrinthula sp. on the Mediterranean seagrass Cymodocea nodosa.Mar. Ecol. Prog. Ser.532101–109. 10.3354/meps11343
43
OlsenY. S.PotouroglouM.Garcias-BonetN.DuarteC. M. (2015). Warming reduces pathogen pressure on a climate-vulnerable seagrass species.Estuar. Coas.38659–667. 10.1007/s12237-014-9847-9
44
OrthR. J.CarruthersT. J. B.DennisonW. C.DuarteC. M. (2006). A global crisis for seagrass ecosystems.BioScience56987–996.
45
OrthR. J.HeckK. L.VanmontfransJ. (1984). Faunal communities in seagrass beds - A review of the influence of plant structure and prey characteristics on predator - prey relationship.Estuaries7339–350. 10.2307/1351618
46
PaviaH.BaumgartnerF.CervinG.EngeS. (2012). “Chemical defences against herbivores,” in Chemical Ecology in Aquatic Systems, edsBrönmarkC.HannsonL. A. (New York, NY: Oxford University Press), 210–235. 10.1093/acprof:osobl/9780199583096.003.0016
47
PérezM.RomeroJ. (1992). Photosynthetic response to light and temperature of the seagrass Cymodocea nodosa and the prediction of its seasonality.Aquat. Bot.4351–62. 10.1016/0304-3770(92)90013-9
48
PfafflM. W.TichopadA.PrgometC.NeuviansT. P. (2004). Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: bestkeeper – Excel-based tool using pair-wise correlations.Biotechnol. Lett.26509–515. 10.1023/b:bile.0000019559.84305.47
49
PozoM. J.CordierC.Dumas-GaudotE.GianinazziS.BareaJ. M.Azcón-AguilarC. (2002). Localized versus systemic effect of arbuscular mycorrhizal fungi on defence responses to Phytophthora infection in tomato plants.J. Exp. Bot.53525–534. 10.1093/jexbot/53.368.525
50
PraschC. M.SonnewaldU. (2013). Simultaneous application of heat, drought, and virus to Arabidopsis plants reveals significant shifts in signaling networks.Plant Physiol.1621849–1866. 10.1104/pp.113.221044
51
RasmussenE. (1977). “The wasting disease of eelgrass (Zostera marina) and its effects on environmental factors and fauna,” in Seagrass Ecosystems: A Scientific Perspective, Vol. 4edsMcRoyC. P.HelfferichC. (New York, NY: M. Dekker), 1–51.
52
SaloT.ReuschT. B. H.BoströmC. (2015). Genotype-specific responses to light stress in eelgrass Zostera marina, a marine foundation plant.Mar. Ecol. Prog. Ser.519129–140. 10.3354/meps11083
53
SpoelS. H.DongX. (2012). How do plants achieve immunity? Defence without specialized immune cells.Nat. Rev. Immunol.1289–100. 10.1038/nri3141
54
SteeleL. T.CaldwellM.BoettcherA.ArnoldT. (2005). Seagrass-pathogen interactions:‘pseudo-induction’ of turtlegrass phenolics near wasting disease lesions.Mar. Ecol. Prog. Ser.303123–131. 10.3354/meps303123
55
SullivanB. K.ShermanT. D.DamareV. S.LiljeO.GleasonF. H. (2013). Potential roles of Labyrinthula spp. in global seagrass population declines.Fungal Ecol.6328–338. 10.1016/j.funeco.2013.06.004
56
SullivanB. K.Trevathan-TackettS. M.NeuhauserS.GoversL. L. (2018). Host-pathogen dynamics of seagrass diseases under future global change.Mar. Pollut. Bull.13475–88. 10.1016/j.marpolbul.2017.09.030
57
Trevathan-TackettS. M.LauerN.LoucksK.RossiA. M.RossC. (2013). Assessing the relationship between seagrass health and habitat quality with wasting disease prevalence in the Florida Keys.J. Exp. Mar. Biol. Ecol.449221–229. 10.1016/j.jembe.2013.10.004
58
TutinT. G. (1938). The autecology of Zostera marina in relation to its wasting disease.New Phytol.3750–71. 10.1111/j.1469-8137.1938.tb06926.x
59
UnderwoodA. J. (1997). Experiments in Ecology: Their Logical Design and Interpretation Using Analysis of Variance.Cambridge: Cambridge University Press.
60
ValentineJ. F.HeckK. L.BusbyJ.Jr.WebbD. (1997). Experimental evidence that herbivory increases shoot density and productivity in a subtropical turtlegrass (Thalassia testudinum) meadow.Oecologia112193–200. 10.1007/s004420050300
61
van LoonL. C.RepM.PieterseC. M. J. (2006). Significance of inducible defense-related proteins in infected plants.Annu. Rev. Phytopathol.44135–162. 10.1146/annurev.phyto.44.070505.143425
62
VergeerL. H. T.AartsT. L.de GrootJ. D. (1995). The ‘wasting disease’ and the effect of abiotic factors (light intensity, temperature, salinity) and infection with Labyrinthula zosterae on the phenolic content of Zostera marina shoots.Aquat. Bot.5235–44. 10.1016/0304-3770(95)00480-n
63
VergeerL. H. T.DeveliA. (1997). Phenolic acids in healthy and infected leaves of Zostera marina and their growth-limiting properties towards Labyrinthula zosterae.Aquat. Bot.5865–72. 10.1016/s0304-3770(96)01115-1
64
VergésA.PérezM.AlcoverroT.RomeroJ. (2008). Compensation and resistance to herbivory in seagrasses: induced responses to simulated consumption by fish.Oecologia155751–760. 10.1007/s00442-007-0943-4
65
WaycottM.DuarteC. M.CarruthersT. J. B.OrthR. J. (2009). Accelerating loss of seagrasses across the globe threatens coastal ecosystems.Proc. Natl. Acad. Sci. U.S.A.10612377–12381. 10.1073/pnas.0905620106
66
WintersG.NelleP.FrickeB.RauchG.ReuschT. B. H. (2011). Effects of a simulated heat wave on photophysiology and gene expression of high-and low-latitude populations of Zostera marina.Mar. Ecol. Prog. Ser.43583–95. 10.3354/meps09213
67
ZidornC. (2016). Secondary metabolites of seagrasses (Alismatales and Potamogetonales; Alismatidae): chemical diversity, bioactivity, and ecological function.Phytochemistry1245–28. 10.1016/j.phytochem.2016.02.004
Summary
Keywords
chemical defense, gene expression, eelgrass, infection, Labyrinthula, multifactorial experiment, marine pathogen
Citation
Jakobsson-Thor S, Brakel J, Toth GB and Pavia H (2020) Complex Interactions of Temperature, Light and Tissue Damage on Seagrass Wasting Disease in Zostera marina. Front. Mar. Sci. 7:575183. doi: 10.3389/fmars.2020.575183
Received
22 June 2020
Accepted
23 September 2020
Published
21 October 2020
Volume
7 - 2020
Edited by
Rui Rosa, University of Lisbon, Portugal
Reviewed by
C. Drew Harvell, Cornell University, United States; Hung Manh Nguyen, University of Naples Federico II, Italy
Updates
Copyright
© 2020 Jakobsson-Thor, Brakel, Toth and Pavia.
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) and the copyright owner(s) 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: Stina Jakobsson-Thor, sjakobsson87@gmail.com
This article was submitted to Global Change and the Future Ocean, a section of the journal Frontiers in Marine Science
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.