- 1Department of Experimental Limnology, Leibniz Institute of Freshwater Ecology and Inland Fisheries, Neuglobsow, Germany
- 2Institute of Biological Sciences, Applied Ecology and Phycology, University of Rostock, Rostock, Germany
- 3Botanic Garden and Botanical Museum Berlin-Dahlem, Freie Universität Berlin, Berlin, Germany
- 4Institute for Inorganic and Analytical Chemistry, Friedrich-Schiller-University Jena, Jena, Germany
- 5Department of Coastal Biology, Argentinean Antarctic Institute, Buenos Aires, Argentina
- 6Department of Basic Sciences, National University of Luján, Luján, Buenos Aires, Argentina
- 7Department of Biology, University of Turku, Turku, Finland
- 8Institute of Biochemistry and Biology, University of Potsdam, Potsdam, Germany
Aquatic ecosystems are frequently overlooked as fungal habitats, although there is increasing evidence that their diversity and ecological importance are greater than previously considered. Aquatic fungi are critical and abundant components of nutrient cycling and food web dynamics, e.g., exerting top-down control on phytoplankton communities and forming symbioses with many marine microorganisms. However, their relevance for microphytobenthic communities is almost unexplored. In the light of global warming, polar regions face extreme changes in abiotic factors with a severe impact on biodiversity and ecosystem functioning. Therefore, this study aimed to describe, for the first time, fungal diversity in Antarctic benthic habitats along the salinity gradient and to determine the co-occurrence of fungal parasites with their algal hosts, which were dominated by benthic diatoms. Our results reveal that Ascomycota and Chytridiomycota are the most abundant fungal taxa in these habitats. We show that also in Antarctic waters, salinity has a major impact on shaping not just fungal but rather the whole eukaryotic community composition, with a diversity of aquatic fungi increasing as salinity decreases. Moreover, we determined correlations between putative fungal parasites and potential benthic diatom hosts, highlighting the need for further systematic analysis of fungal diversity along with studies on taxonomy and ecological roles of Chytridiomycota.
Introduction
Fungi are morphologically, phylogenetically, and functionally diverse. They constitute a well-founded component of terrestrial ecology due to more than 100 years of research that has highlighted their role in biogeochemical cycling and promoting biodiversity (Peay et al., 2016). Aquatic ecosystems, in particular the oceans, are frequently overlooked as fungal habitats with a systematic analysis of fungal diversity and their ecological roles still in their infancy. However, there is increasing evidence that fungal diversity in aquatic ecosystems is greater and more important than previously considered (Shearer et al., 2007; Amend et al., 2019; Grossart et al., 2019). Despite sampling efforts of fungi in aquatic habitats being low compared to terrestrial environments (Rojas-Jimenez et al., 2017), molecular analyses of environmental DNA samples reveal a great diversity of novel fungal sequences, the so-called Dark Matter Fungi (Grossart et al., 2016). Although recent advances in DNA-sequencing technology have revealed that fungi are highly abundant in marine environments (Comeau et al., 2016; Taylor and Cunliffe, 2016; Tisthammer et al., 2016; Hassett et al., 2019a; Banos et al., 2020), their ecological functions and interactions with other microorganisms remain largely unexplored and missing from current general concepts (Amend et al., 2019; Richards et al., 2021).
Based on recent studies, marine fungal communities are dominated by members of the phylum Ascomycota (Tisthammer et al., 2016; Amend et al., 2019; Hassett et al., 2019a) and Chytridiomycota (Comeau et al., 2016; Hassett and Gradinger, 2016; Hassett et al., 2017; Richards et al., 2021). Chytridiomycota, frequently referred to as chytrids, are often recognized as zoosporic virulent parasites on phytoplankton which play significant roles in controlling population sizes of their hosts (Ibelings et al., 2004; Kagami et al., 2004; Rasconi et al., 2011) and altering food web structure by transferring carbon and energy between trophic levels (Rasconi et al., 2012, 2020; Sime-Ngando and Twiss, 2012; Klawonn et al., 2021b). The majority of quantitative studies of zoosporic fungi to date have been laboratory-based, with only a small number of field studies conducted to assess their importance into the broader environment. In addition, zoosporic fungi are difficult to identify using solely morphological characteristics. Although direct examination of the diversity of zoosporic fungi using high-throughput sequencing has been increasingly applied (Monchy et al., 2011; Tedersoo et al., 2015, 2018; Seto et al., 2017; Li et al., 2021; Richards et al., 2021), data are still lacking for polar regions.
Microphytobenthos (MPB) comprise phototrophic benthic microalgae living in intertidal and shallow subtidal sediments and hard substrates, among which benthic diatoms are typically the dominant organisms and are known as preferential hosts for fungal parasites such as chytrids (Scholz et al., 2016a). Although the amount of research on fungal parasites infecting pelagic forms is growing, data about microphytobenthic taxa are sparse (Wulff et al., 2009; Scholz et al., 2016b). Generally, benthic microalgae are among the main contributors to primary production in coastal zones, particularly those poor in nutrients and organic matter. They can contribute up to 42% of the marine benthic primary production (Woelfel et al., 2010). Specifically, according to Cahoon (1999), the global production of benthic microalgae ranges from 8.9 to 14.4 Gt C m–2 year–1, representing approximately 20% of the global ocean production. Hence, they are important suppliers of organic carbon to grazers and sediment feeding macro- and meiofauna (Middelburg et al., 2000; Oakes et al., 2010). Except for sediments, benthic microalgae are also a major and ecologically important components of epilithic biofilms (Cahoon, 1999). They exert further important ecological roles in marine shallow water environments as they occur at the sediment-water interface and thus directly influence, for example, vertical nutrient exchange processes (Sundbäck et al., 2000). Dissolved organic carbon is regularly released by benthic diatom excretion of extrapolymeric substances (EPS), which act as stabilizing compounds against hydrodynamic sediment erosion as well as bacterial substrate (de Brouwer et al., 2005; Aslam et al., 2012).
Marine environments, especially polar regions, which have historically experienced minimum human disturbance, are a vast reservoir of microbial diversity. Their clean air, water and ice are of great importance to science for understanding different aspects of the Earth’s environment. Hence, they have been drawing more attention and as such, they have not escaped the negative impacts of human activity. Their unique marine ecosystems are affected on local and regional scales by overfishing, pollution, introduction of invasive species and exploitation of mineral reserves, oil and gas (Aronson et al., 2011). Recently, in the light of global warming, the specific focus is on the loss of sea ice (Stammerjohn et al., 2008). During summer and winter, Antarctic sea ice does not show a significant overall trend (Parkinson and DiGirolamo, 2021). Nevertheless, south and west of the Antarctic Peninsula, one region has shown a persistent decline (Grosfeld et al., 2013). In February 2020 weather stations recorded the hottest temperature on record for Antarctica, reaching 18.3°C (64.9°F) (Rocha Francelino et al., 2021). This heatwave was the third major melt event of the 2019–2020 summer, following warm spells in November 2019 and January 2020. Consequently, polar aquatic ecosystems face local changes that include higher water temperatures and altered patterns in light penetration, dust deposition, sediment load and changes in salinity due to intensified glacier melt and subsequent terrestrial runoff (Eraso and Dominguez, 2007; Hernández et al., 2019). Similar climate-change-driven changes in abiotic factors were shown to affect microbial community structure, e.g., chytrid community composition in Arctic waters (Kilias et al., 2020), but for Antarctica, such data are scarce. Slow-moving or sedentary benthic communities in shallow nearshore areas are most vulnerable to these processes, with their species richness and diversity strongly impacted. Moon et al. (2015) discuss the impact of glacial retreat on epibenthic megafaunal assemblages on the West Antarctic Peninsula (WAP) and note significant differences in assemblages related to distance from the glacier, substrate grain size, and organic content. Furthermore, Braeckman et al. (2021) showed that benthic communities along the WAP shift from net autotrophy to net heterotrophy when affected by glacial meltwater. Nevertheless, more information is needed to gain a comprehensive understanding of the effects of glacial meltwater on benthic, especially fungal, communities. Although a growing amount of research focuses on fungal parasites in aquatic ecosystems, polar regions remain mainly unexplored. Rojas-Jimenez et al. (2017) reported early diverging lineages within Chytridiomycota and Cryptomycota as dominant among fungal sequences in the Antarctic Dry Valley lakes. Yet, in Antarctic marine systems, further research needs to be done to gain more detailed insights into their diversity and abundance changes and their ecological role under current climate change.
Despite the potentially important effect of fungal parasites on benthic diatoms, to our knowledge, no published experimental or field studies have been conducted on their diversity in Antarctic coastal waters. Therefore, in this study, we aim to qualitatively describe fungal community composition in Antarctic benthic habitats. Our main focus being on the diversity of parasitic fungi and the co-occurrence with their benthic microalgal hosts along the salinity gradient in Potter Cove, Antarctic Peninsula, which is impacted by a summer induced glacial meltwater runoff. We hypothesized that members of the fungal phylum Chytridiomycota predominate fungal taxa in Antarctic shallow coastal waters and exhibit a still unexplored biodiversity with many undescribed species.
Materials and Methods
Sites and Sample Collection
The study was conducted at the German-Argentine Dallmann Laboratory in the Argentine Scientific Station “Dr. Carlini,” which is located at Potter Cove, King George Island/25 de Mayo I, Antarctic Peninsula (62°14′S, 58°31′W). This area combines zones of glacier fronts, rocky shores and extensive seabed overlaid with sand or sediment. Sampling took place in January–February 2020 in several locations around Carlini Station and the Antarctic Specially Protected Area (ASPA) 132 (Figure 1). The sites included marine shallow water locations and open water locations down to 60 m depth, limnic locations that included lakes and ponds, and brackish meltwater runoff (Table 1). The area provided diverse sediment and biofilm habitats for the sampling of benthic diatoms. Sediment surface samples were taken using Plexiglas sediment corer, and biofilm samples were taken by scratching off the surface of at least five stones with a sterilized knife. Sampling of the ocean bottom was done with the help of Argentinian army divers. Each sample was divided into two parts and prepared the following way: (i) the part for DNA analysis was fixed with 70% ethanol, and the samples were kept at –20°C; (ii) the part for microscopy was fixed with Lugol’s solution.
Figure 1. Sampling sites in the Potter Cove area (South Shetland Islands, western Antarctic Peninsula).
DNA Extraction and Sequencing
DNA from the microorganisms in ethanol fixed sediment and biofilm samples was extracted according to a modified protocol described by Nercessian et al. (2005). Briefly, cell lysis was achieved using small (0.1–1 mm) zirconia-silica beads that were suspended in cetyltrimethyl ammoniumbromide (CTAB), to which anion surfactants sodium dodecyl sulfate and N-Lauroylsarcosin, proteinase K and phenol–chloroform–isoamylalcohol were added. DNA purification was facilitated by the addition of chloroform–isoamylalcohol and polyethylene glycol (PEG). Finally, DNA was precipitated at 4°C, washed with ethanol, air-dried and dissolved in ultra-pure water. The detailed protocol is available in Supplementary Material. PCR, library preparation and sequencing was undertaken at mrDNA laboratories (Shallow Waters, Texas, United States). The V8 region of the 18S rRNA gene was amplified using primers 18S-82F (5′- GAAACTGCGAATGGCTC-3′) and Ek-516R (5′- ACCAGACTTGCCCTCC-3′) (Hassett et al., 2019b) for molecular characterization of diatoms, followed by library preparation (2 × 300 bp) and sequencing on a MiSeq (Illumina) platform. This dataset will be referred to as EukSSU further in the text. For molecular characterization of fungi we amplified the LSU D1 region of rRNA using primers ITS4ngsF (5′-GCA TAT CAA TAA GCG SAG GAA-3′) and LF402R (5′-TTC CCT TTY ARC AAT TTC AC-3′) (Tedersoo et al., 2015). This dataset will be referred to as FunLSU further in the text.
Taxonomic Identification and Amplicon Sequence Variant Generation
Sequences were processed in R (ver. 4.0.3) using the DADA2 pipeline (ver. 1.8) (Callahan et al., 2016). Primers were removed from demultiplexed reads using cutadapt ver. 3.5 (Martin, 2011). The primer-free sequences were then filtered and trimmed to remove low-quality sequences. The DADA2 algorithm was used to infer amplicon sequence variants (ASVs). Paired-end reads were merged to obtain full denoised sequences. Chimeric sequences were removed, with the exception of 18S rRNA chimeric sequences, which were identified and removed using DECIPHER online tool ver. 11.4 (Wright et al., 2012). 18S rRNA gene sequences were run against the SILVA SSU ver. 132 database (Quast et al., 2013) for identification of all eukaryotic taxa and LSU sequences were run against LSU database using RDP Classifier and Fungal LSU training set 11 (Wang et al., 2007) for identification of fungal taxa.
Data Analysis
Data processing, visualizations and statistical analysis were performed in R. Prior to any analysis, singletons were removed from both EukSSU and FunLSU datasets. Differences in community structure in respect to water types between samples (beta diversity) were calculated using a Bray-Curtis dissimilarity measure, using phyloseq package (McMurdie and Holmes, 2013), and visualized through non-metric multidimensional scaling (NMDS) ordination. Permutational multivariate analysis of variance (PERMANOVA) was used to test the effect of different water types on community structure using the “Adonis” function in the R package vegan (Oksanen et al., 2019). The statistical significance was calculated with a post hoc pairwise t-test using “pairwise adonis” function (Martinez Arbizu, 2020) with 999 permutations. P-values of the pairwise t-test were adjusted with the Bonferroni method. To calculate alpha diversity values, microbiome package was used (Lahti and Shetty, 2017). The graphical representation of results was realized using the R package ggplot2 (Wickham, 2016). Ternary plots were made using ggtern package (Hamilton and Ferry, 2018). Map of sampling points was generated in software QGIS 3.16.14 (QGIS Development Team, 2021) using a geospatial data package and visualization platform Quantarctica (Matsuoka et al., 2021).
Network Analysis
To infer whether the presence of suitable hosts is correlated with the presence of fungal parasites in collected samples, we performed a co-occurrence network analysis using the SparCC algorithm (Friedman and Alm, 2012) implemented in FastSpar (Watts et al., 2019). Prior to the network analysis, ASVs with < 10% prevalence and occurring in less than three samples were removed from both EukSSU and FunLSU datasets. Data was further processed in R using WGCNA package (Langfelder and Horvath, 2008). We maintained the independence of both datasets and converted the correlation coefficients obtained from FastSpar into a topological overlap matrix (TOM) as the first step to identify modules of co-occurring taxa across our samples. Using parallel minimum (pmin) we calculated a consensus TOM from the EukSSU and FunLSU TOMs which was then used as an input for hierarchical average linkage clustering. To identify consensus modules we used cutreeDynamic with a “deepSplit” of two and “minModuleSize” of five. During this step ASVs are assigned to certain modules depending on their TOM-based topology. Further co-occurrence between identified EukSSU and FunLSU modules was determined by extracting their respective Eigen values and calculating Spearman rank correlation coefficients. Co-occurrence between individual EukSSU and FunLSU ASVs was calculated by multiplying for each ASV pair, the between module Spearman correlation, by the module membership (range 0–1) for each ASV. In practice, this ensures that the two ASVs with the highest module membership from each respective module are most strongly correlated, and corr values decrease with decreasing module membership. We visualized for selected module pairs, weighted correlation coefficients (> 0.5) between ASVs using Cytoscape 3.8.2 software.
Microscopy
CalcoFluor White (CFW) staining approach along with epifluorescence microscopy, as suggested by Rasconi et al. (2009), was used to detect and identify fungal parasites attached to an algal hosts. CFW is a non-specific, fluorescent dye that binds to chitin in fungal cell walls but also cellulose, that is often present in cell wall of some algae and fungi-like organisms (Kagami et al., 2007; Priest et al., 2021). Pre-treatment of sediment samples was done using ultrasound (3 × 10 s) for the mechanical disruption of diatom-sand agglomerates (Scholz et al., 2014). Staining was performed as described in Klawonn et al. (2021a). Five microliter of CFW (1 mg/ml) was added to 1 ml sample and incubated for 15 min at room temperature. Samples were then transferred in Utermöhl counting chambers, after which they were incubated for 10 min, allowing the cells to sink to the bottom of the chamber. For further visualization, Nikon Eclipse Ti2 inverted microscope was used and the samples were screened for fungal parasites on 600x magnification.
Results
With this study it was intended to highlight the prevalence of fungi in Antarctic benthic environments, determine the impact of glacial meltwater on benthic diatoms and their associated fungal communities and to identify correlation in their presence. To this end, freshwater, brackish and marine benthic environments were sampled and two datasets obtained: EukSSU containing the whole eukaryotic community and FunLSU containing fungal community. When examining the EukSSU community composition, a total of 3,790 sequence variants (ASVs) were identified of which 145 were classified as fungi. The relative proportion of these 145 fungal ASVs varied between different water types, from 0.97% in marine, 7.4% in brackish, to 8.9% in limnic habitats. LSU primers were not strictly fungi-specific and out of 1,204 identified ASVs, of which 727 were fungi, the remaining identified as Metazoa (262), Stramenopiles (131), Amoebozoa (35), Chlorophyta (33), Rhodophyta (12), and Alveolata (4). Before any downstream analysis, all non-fungal sequences were removed from FunLSU dataset. Nevertheless, relative proportions of fungal ASVs in FunLSU confirmed differences observed in EukSSU and 51.3% of sequences were classified as fungi in marine, 73.2% in brackish and 69.3% in limnic environments. Calculating alpha diversity measures, we note that eukaryotic and fungal diversity were lower in marine samples compared to either brackish or limnic samples (Figure 2). For eukaryotes, the dominance of a few taxa was the largest impact on diversity in marine environments, with Pielou’s evenness higher for eukaryotic communities in brackish and limnic samples (Figure 2A). In contrast, fungal communities had more even species composition amongst all three water types, with marine samples being less diverse compared to brackish or limnic samples (Figure 2B). Due to artificial variation introduced through subsampling, the omission of valid data through loss of sequence counts, or exclusion of samples with small library sizes, sequences were not rarefied prior to calculating alpha diversity (Cameron et al., 2021). As Willis (2019) discussed, environments can be identical with respect to one alpha diversity metric, but the different abundance structures will induce different biases when rarefied.
Figure 2. Distribution of the alpha diversity estimators according to sampled water types in (A) 18S and (B) LSU data. This includes Pielou’s evenness, observed ASV richness and Shannon diversity index.
Fungal Community Composition
To determine community composition, we carried out a non-NMDS analysis and it confirmed a distinct separation based on sample habitat (marine vs. brackish vs. limnic) (PERMANOVA, p = 0.001) in both EukSSU (Figure 3A) and FunLSU (Figure 3B) datasets. Marine samples were separated from limnic (post hoc pairwise t-test, p = 0.003) and brackish samples (p = 0.006), but no significant difference was observed between limnic and brackish sites. We detected a difference based on sample type (biofilm vs. sediment vs. water) (PERMANOVA, p = 0.001). Post hoc tests confirmed a difference between biofilm and sediment samples (p = 0.003). Summary of all the statistical tests performed can be found in Supplementary Table 1 for EukSSU dataset and Supplementary Table 2 for FunLSU dataset.
Figure 3. NMDS multivariate clustering of communities according to location and sample type. In both (A) EukSSU and (B) FunLSU datasets marine samples separated from limnic (post hoc pairwise t-test, p = 0.003) and brackish samples (post hoc pairwise t-test, p = 0.006) as did biofilm from sediment samples (post hoc pairwise t-test, p = 0.003).
Ascomycota and Chytridiomycota were the most abundant taxa (FunLSU dataset) in marine Antarctic benthic habitats, with no apparent distinction based on habitat or sample type (Figure 4A). Ascomycota represented 49% of the fungal reads and 25.4% of the fungal ASVs, while Chytridiomycota represented 22.4% of the fungal reads and 38.9% of the ASVs. In addition, a high diversity of the phylum Chytridiomycota was observed (Figure 4B), with Spizellomycetales and Monoblepharidales being the most abundant. In the EukSSU dataset, it was noted that fungi were more dominant in freshwater and brackish environments (8.9 and 7.4%), relative to marine environments (0.97%) (Figure 5). This was consistent with an increased diversity of fungi, as observed in both the EukSSU and FunSSU datasets.
Figure 4. Community composition in the Antarctic benthic habitats, considering the abundance distribution in composition in sampled water types (A) of the whole fungal community (B) and the phylum Chytridiomycota relative to the total fungal community.
Figure 5. Ternary plot showing relative abundance of fungal ASVs in LSU data. Each point represents an ASV. The position represents the relative abundance of the ASV with respect to each environment.
Co-occurrence of Fungal Parasites and Their Hosts
Using microscopy, we observed Calcofluor White stained sporangia attached to benthic algae which indicated active fungal infections (Supplementary Figure 1). Furthermore, we performed a co-occurrence network analysis to identify correlations between fungi and benthic algae to identify putative hosts. Due to the compositional nature of both the EukSSU and FunLSU datasets, it was not possible to estimate a direct linear or sparse correlation between the two datasets. Rather, we calculated sparse correlations amongst ASVs within each dataset, clustered ASVs into modules, extracted Eigen values for each module and used these values to identify correlated modules. In total, we identified seven modules in FunLSU and eight modules in EukSSU dataset. Modules largely reflected differences between habitats. FunLSU modules typically always contained at least one Chytrid ASV (Supplementary Figure 2), with Mesochytrium, Maunachytrium, Betamyces, and Spizellomyces being most abundant. For four out of the six FunLSU modules containing chytrids, we observed a correlation between that module and a EukSSU module containing a high abundance of diatoms (Supplementary Figure 3). For two FunLSU modules that contained chytrids, we did not observe a correlation with a EukSSU module, although all EukSSU modules contained diatoms (Supplementary Figure 4).
To better illustrate the correlation between diatoms and zoosporic fungi, we visualized weighted correlation coefficients amongst ASVs obtained from two EukSSU and FunSSU module pairs. In EukSSU3 and FunLSU5 module pair (Figure 6A) we visualized only those correlations that had a correlation value higher than 0.2. Both modules contained ASVs enriched from marine benthic environments. Specifically, we observed Fragilariopsis cylindrus, a pennate, benthic sea-ice diatom found both in Arctic and Antarctic waters (Helmcke and Krieger, 1953), three different Navicula spp. (N. phyllepta, N. rhynchocephala, N. perminuta) and Pleurosigma intermedium, all marine benthic diatoms that could serve as potential hosts for parasitic chytrids. Modules EukSSU4 and FunLSU2 (Figure 6B) contained ASVs enriched in brackish/limnic benthic habitats. In benthic environments, we observed a diverse number of benthic diatoms associated with some Chytridiomycota. In this case, we observed Betamyces (Rhizophydiales), chytrids mostly recognized as parasites on different phytoplankton species (Christaki et al., 2017), and Maunachytrium, for which there is no solid evidence for a specific lifestyle. Most abundant in both networks were ASVs classified as unknown Chytridiomycota, and they were highly correlated with certain diatom ASVs.
Figure 6. Co-occurrence network analysis. (A) EukSSU3 and FunLSU5 modules (B) EukSSU4 and FunLSU2 module. Blue—Baccilariophyceae; dark green—Chytridiomycota; yellow—other Fungi; brown—Opisthokonta, red—Stramenopiles; light green—Chlorophyta.
Discussion
Climate warming in Antarctic environments has been associated with glacier retreat and increased ice melting (Lee et al., 2017) which, in turn, change the vertical structure of the water column, especially in Antarctic shallow coastal environments. During the summer months (December—March), meltwaters on the Western Antarctic Peninsula, where this study was undertaken, occur regularly. It has been shown that meltwater transports high particle loads to the coastal water thereby reducing incident light conditions and salinity that affect microbial community dynamics, leading to changes in species composition (Schloss et al., 2008; Antoni et al., 2020). In this study, we showed different proportions of fungal sequences between different water types. Hence, we assume a shaping impact of salinity on both the eukaryotic and fungal diversity in Potter Cove. Brackish sites, receiving strong inputs of fresh meltwater, exhibited a higher eukaryotic and fungal diversity than that of marine sites. Moreover, many of the same fungal taxa/ASVs found in freshwater environments were also determined in brackish habitats, suggesting fresh meltwater runoff, including terrestrial soil particles, might be a direct source of this increased diversity. Furthermore, it is tempting to speculate that the increased richness of fungi and the decreased dominance of individual algal hosts in brackish/freshwater environments, as shown with alpha diversity measures, might be linked. We argue that this observation can be directly linked with fungal parasites playing an active role in maintaining the diversity of algae (van Donk, 1989; Sime-Ngando and Twiss, 2012). Alternatively, this might be indirect in the sense that fertilization (Schloss et al., 2008) opens up new niches for algae, with fungal diversity, of both saprotrophic and/or parasitic fungi, promoted by an increase in the diversity of autochthonous and allochthonous organic matter (Wang et al., 2017). Finally, we conclude that the increased evenness in brackish/freshwater habitats compared to marine habitats may indicate that glacial meltwaters reduce selective pressure and increase niche opportunities for algae, possibly due to increased nutrients as observed by Death et al. (2014). In contrast, observed differences in fungal diversity may have arisen from a decreased richness in marine habitats relative to brackish/freshwater.
Our results are supported by the study of Kilias et al. (2020), who showed that chytrid fungi are primarily encountered at sites influenced by sea ice melt. The evolutionary history of Chytridiomycota suggests that they originated under brackish-/freshwater-like conditions and thus evolved successively in parallel with the host organisms or existing food (Berbee et al., 2017). It remains an open question whether marine Chytridiomycota proliferate due to lower competition in high salinity water or whether specialization of some taxa has occurred over an evolutionary time scale (Yang et al., 2021). Moreover, a substantial proportion of fungal ASVs in our FunLSU dataset were classified as terrestrial fungi. For example, Spizellomycetales, the most abundant genera among chytrids, are described as soil-inhabiting fungi (Barr, 1980). As a result, we hypothesize that spores and fungal hyphae may arrive in Antarctic coastal waters via glacial/terrestrial runoff from ice melting followed by deposition due to sedimentation alongside other particles. In particular, larger areas of bare soil at higher temperatures lead to an increased deposition of soil particles on the melting glaciers, staining them dark and even decreasing their albedo and increasing their melting (Schmitt et al., 2015). This notion is supported by the presence of Monoblepharidales, which are described as typical freshwater fungi (Sparrow, 1933), and constitute the dominant genera among fungal taxa in the sampled marine environments. Such sediment accumulation driven by increased glacier retreat may decrease the rate of microphytobenthic (MPB) primary production. As earliers discussed by Hoffmann et al. (2019), MPB might survive the consequences but their contribution to the overall primary production as a carbon resource may decline. Consequently, food competition of the benthic heterotrophic community will increase, with unpredictable consequences for biomass, density, structure and diversity of the benthic community and food web structure.
Using microscopy and performing co-occurrence network analysis, we found preliminary evidence suggesting that zoosporic fungi in Antarctic waters are involved in parasitic relationships with benthic diatoms. In summary, in most cases, parasitic chytrids were associated with the presence of Bacillariophyceae, but some diatom taxa occurred in some cases independent of parasitic chytrids. Furthermore, unclassified Chytridiomycota and Spizellomyces, which are known as terrestrial parasites on other chytrid species and soil nematodes (Lozupone and Klein, 2002), were non-specifically associated with multiple potential diatom hosts. It is tempting to speculate that these chytrids exhibit a parasitic lifestyle in aquatic habitats with broad host tolerance, considering that these ASVs represent a large hidden diversity. Whilst network approaches may provide hints for interactions between fungi and diatoms, these are ultimately correlations. Cultivation and targeted single cell approaches that identify host-parasite pairs and better discriminate between parasitic and saprotrophic interactions are still needed (van den Wyngaert et al., 2017). Despite this, our study highlights that diatom-parasite interactions are likely in polar benthic environments and may have significant ecological roles.
Conclusion
Our work suggests that salinity has a major impact in shaping eukaryotic microbial communities in shallow coastal waters of Antarctica. In particular, we determined a significant difference in eukaryotic and fungal diversity between sampled habitats, with both eukaryotic and fungal diversity higher in brackish/freshwater environments. On the one hand, considering the potential parasitic lifestyle of dominant fungal taxa, we may conclude that fungal parasites play an active role in maintaining the diversity of benthic algae. On the other hand, increased ice melting during summer months may cause nutrient loading in brackish waters resulting in increased diversity. It remains an open question whether the positive effect on fungal diversity also has a positive or a negative effect on the ecosystem on a broad-scale, considering parasitic interactions. The factor of inputs of meltwater is also supported by the fact that freshwater and brackish environments shared many of the same fungal taxa/ASVs and that a substantial proportion of fungal ASVs in the FunLSU dataset were classified as terrestrial fungi. In terms of parasitic lifestyle, while many of the chytrid ASVs were found to be non-specifically associated with multiple potential diatom hosts, drawing mechanistic conclusion about host-parasite interactions is difficult because the most abundant ASVs in our networks were classified as unknown Chytridiomycota. Systematic analysis of fungal diversity in Antarctic regions is still in its infancy, and the taxonomy and ecological functions of Chytridiomycota species remain largely unexplored.
Our findings are a significant step toward illuminating the diversity of fungal communities and occurrence of, what it may be, fungal parasites in the Antarctic region. Benthic diatoms play a key role in polar food webs. Hence, studying chytrid diversity and their interactions with the algal hosts will provide us with a comprehensive understanding of polar ecosystems. The effects of salinity changes on fungi-microalgae interactions and microalgal communities, as well as its consequences for the coastal ecosystem of Potter Cove, are therefore subject to future studies.
Data Availability Statement
The original contributions presented in the study are publicly available. This data can be found here: (https://www.ebi.ac.uk/ena/browser/view/PRJEB49266/PRJEB49266).
Author Contributions
AL, H-PG, MLQ, GLC, JZ, UK, and TW organized an expedition to Antarctica and collected samples. AL performed microscopy and extracted DNA from collected samples. DI and JW analyzed the data. DI wrote the manuscript. SVdW assisted with her knowledge during the analysis. H-PG, JW, MLQ, GLC, JZ, UK, TW, and SVdW edited and provided feedback on the original draft. All authors approved the submitted version.
Funding
The study was performed in the frame of the Priority Programme of the German Science Foundation (DFG) entitled “Antarktisforschung mit vergleichenden Untersuchungen in arktischen Eisgebieten” (SPP 1158) in the subproject “Polar parasites” (GR1540/33-1) given to H-PG and UK (KA899/39-1). The DFG-SPP 1158 subproject “Benthic diatom diversity and biogeography” ZI 1628/2-1 is given to JZ and UK (KA899/39-1). GLC was supported by grants from IAA-DNA (PIH18), PADI Foundation (47918/2020), MINCYT-BMBF Programme (AL/17/06-01DN18024), and Universidad Nacional de Luján (Disp CDD-CB 296/18, 343/19, and 69/21). GLC and MLQ were supported by grants from ANPCyT-DNA (PICTO 2010-0116 and PICT 2017-2691). This study presents an outcome of Coastcarb (Funding ID 872690, H2020-MSCA-RISE-2019). Publication is funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) – Project Number 491466077.
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.
The handling editor declared a past co-authorship with one of the authors UK.
Publisher’s Note
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.
Acknowledgments
We would like to thank the team of the Argentinian Antarctic Research Station “Carlini” for their support and logistics as well as AWI Logistics department. We also appreciate the support of the MIBI group at IGB Stechlin for their frequent scientific input and discussions. We acknowledge the support of the Deutsche Forschungsgemeinschaft and Open Access Publishing Fund of University of Potsdam.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.805694/full#supplementary-material
Abbreviations
EukSSU, dataset containing whole eukaryotic community; FunLSU, dataset containing only fungal community.
References
Amend, A., Burgaud, G., Cunliffe, M., Edgcomb, V. P., Ettinger, C. L., Gutiérrez, M. H., et al. (2019). Fungi in the marine environment: open questions and unsolved problems. mBio 10:e01189-18. doi: 10.1128/mBio.01189-18
Antoni, J. S., Almandoz, G. O., Ferrario, M. E., Hernando, M. P., Varela, D. E., Rozema, P. D., et al. (2020). Response of a natural Antarctic phytoplankton assemblage to changes in temperature and salinity. J. Exp. Mar. Biol. Ecol. 532:151444. doi: 10.1016/j.jembe.2020.151444
Aronson, R. B., Thatje, S., McClintock, J. B., and Hughes, K. A. (2011). Anthropogenic impacts on marine ecosystems in Antarctica. Ann. N. Y. Acad. Sci. 1223, 82–107. doi: 10.1111/J.1749-6632.2010.05926.X
Aslam, S. N., Cresswell-Maynard, T., Thomas, D. N., and Underwood, G. J. C. (2012). Production and characterization of the intra- and extracellular carbohydrates and polymeric substances (EPS) of three sea-ice diatom species, and evidence for a cryoprotective role for EPS. J. Phycol. 48, 1494–1509. doi: 10.1111/JPY.12004
Banos, S., Gysi, D. M., Richter-Heitmann, T., Glöckner, F. O., Boersma, M., Wiltshire, K. H., et al. (2020). Seasonal dynamics of pelagic mycoplanktonic communities: interplay of taxon abundance. Temporal occurrence, and biotic interactions. Front. Microbiol. 11:1305. doi: 10.3389/fmicb.2020.01305
Barr, D. J. S. (1980). An outline for the reclassification of the Chytridiales, and for a new order, the Spizellomycetales. Can. J. Bot. 58, 2380–2394. doi: 10.1139/b80-276
Berbee, M. L., James, T. Y., and Strullu-Derrien, C. (2017). Early diverging fungi: diversity and impact at the dawn of terrestrial life. Annu. Rev. Microbiol. 71, 41–60. doi: 10.1146/annurev-micro-030117-020324
Braeckman, U., Pasotti, F., Hoffmann, R., Vázquez, S., Wulff, A., Schloss, I. R., et al. (2021). Glacial melt disturbance shifts community metabolism of an Antarctic seafloor ecosystem from net autotrophy to heterotrophy. Commun. Biol. 4:148. doi: 10.1038/s42003-021-01673-6
Cahoon, L. B. (1999). The role of benthic microalgae in neritic ecosystems. Oceanogr. Mar. Biol. Annu. Rev. 37, 55–94. doi: 10.1201/9781482298550-4
Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., and Holmes, S. P. (2016). DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869
Cameron, E. S., Schmidt, P. J., Tremblay, B. J.-M., Emelko, M. B., and Müller, K. M. (2021). To rarefy or not to rarefy: enhancing diversity analysis of microbial communities through next-generation sequencing and rarefying repeatedly. bioRxiv [Preprint]. doi: 10.1101/2020.09.09.290049
Christaki, U., Genitsaris, S., Monchy, S., Li, L. L., Rachik, S., Breton, E., et al. (2017). Parasitic eukaryotes in a meso-eutrophic coastal system with marked Phaeocystis globosa blooms. Front. Mar. Sci. 4:416. doi: 10.3389/fmars.2017.00416
Comeau, A. M., Vincent, W. F., Bernier, L., and Lovejoy, C. (2016). Novel chytrid lineages dominate fungal sequences in diverse marine and freshwater habitats. Sci. Rep. 6:30120. doi: 10.1038/srep30120
de Brouwer, J. F. C., Wolfstein, K., Ruddy, G. K., Jones, T. E. R., and Stal, L. J. (2005). Biogenic stabilization of intertidal sediments: the importance of extracellular polymeric substances produced by benthic diatoms. Microb. Ecol. 49, 501–512. doi: 10.1007/S00248-004-0020-Z
Death, R., Wadham, J. L., Monteiro, F., le Brocq, A. M., Tranter, M., Ridgwell, A., et al. (2014). Antarctic ice sheet fertilises the Southern Ocean. Biogeosciences 11, 2635–2643. doi: 10.5194/BG-11-2635-2014
Eraso, A., and Dominguez, C. (2007). Physicochemical characteristics of the subglacier discharge in Potter Cove, King George Island, Antarctica. Karst Cryokarst Stud. Faculty Earth Sci. 45, 111–122.
Friedman, J., and Alm, E. J. (2012). Inferring correlation networks from genomic survey data. PLoS Comput. Biol. 8:e1002687. doi: 10.1371/JOURNAL.PCBI.1002687
Grosfeld, K., Treffeisen, R., Asseng, J., Bartsch, A., Bräuer, B., Fritzsch, B., et al. (2013). Online sea-ice knowledge and data platform. Polarforschung 85, 143–155. doi: 10.2312/polfor.2016.011
Grossart, H. P., van den Wyngaert, S., Kagami, M., Wurzbacher, C., Cunliffe, M., and Rojas-Jimenez, K. (2019). Fungi in aquatic ecosystems. Nat. Rev. Microbiol. 17, 339–354. doi: 10.1038/s41579-019-0175-8
Grossart, H. P., Wurzbacher, C., James, T. Y., and Kagami, M. (2016). Discovery of dark matter fungi in aquatic ecosystems demands a reappraisal of the phylogeny and ecology of zoosporic fungi. Fungal Ecol. 19, 28–38. doi: 10.1016/j.funeco.2015.06.004
Hamilton, N. E., and Ferry, M. (2018). Ggtern: ternary diagrams using ggplot2. J. Stat. Softw. 87, 1–17. doi: 10.18637/jss.v087.c03
Hassett, B. T., Borrego, E. J., Vonnahme, T. R., Rämä, T., Kolomiets, M. V., and Gradinger, R. (2019a). Arctic marine fungi: biomass, functional genes, and putative ecological roles. ISME J. 13, 1484–1496. doi: 10.1038/s41396-019-0368-1
Hassett, B. T., Ducluzeau, A. L. L., Collins, R. E., and Gradinger, R. (2017). Spatial distribution of aquatic marine fungi across the western Arctic and sub-arctic. Environ. Microbiol. 19, 475–484. doi: 10.1111/1462-2920.13371
Hassett, B. T., and Gradinger, R. (2016). Chytrids dominate arctic marine fungal communities. Environ. Microbiol. 18, 2001–2009. doi: 10.1111/1462-2920.13216
Hassett, B. T., Thines, M., Buaya, A., Ploch, S., and Gradinger, R. (2019b). A glimpse into the biogeography, seasonality, and ecological functions of arctic marine Oomycota. IMA Fungus 10:6. doi: 10.1186/s43008-019-0006-6
Helmcke, J. G., and Krieger, W. (1953). Diatomeenschalen im elektronenmikroskopischen Bild. | Citations - Diatoms of North America. Available online at: https://diatoms.org/citations/helmcke_jg_and_krieger_w-1953-diatomeenschalen_im_elektronenmikroskopischen (accessed August 5, 2021).
Hernández, E. A., Lopez, J. L., Piquet, A. M. T., mac Cormack, W. P., and Buma, A. G. J. (2019). Changes in salinity and temperature drive marine bacterial communities’ structure at Potter Cove, Antarctica. Polar Biol. 42, 2177–2191. doi: 10.1007/S00300-019-02590-5/FIGURES/6
Hoffmann, R., Al-Handal, A. Y., Wulff, A., Deregibus, D., Zacher, K., Quartino, M. L., et al. (2019). Implications of glacial melt-related processes on the potential primary production of a microphytobenthic community in potter cove (Antarctica). Front. Mar. Sci. 6:655. doi: 10.3389/FMARS.2019.00655/BIBTEX
Ibelings, B. W., de Bruin, A., Kagami, M., Rijkeboer, M., Brehm, M., and van Donk, E. (2004). Host parasite interactions between freshwater phytoplankton and chytrid fungi (Chytridiomycota). J. Phycol. 40, 437–453. doi: 10.1111/j.1529-8817.2004.03117.x
Kagami, M., Arnout De Bruin, A. E., Bas, A. E., Ibelings, W., Ellen, A. E., and Donk, V. (2007). Parasitic chytrids: their effects on phytoplankton communities and food-web dynamics. Hydrobiologia 578, 113–129. doi: 10.1007/s10750-006-0438-z
Kagami, M., van Donk, E., de Bruin, A., Rijkeboer, M., and Ibelings, B. W. (2004). Daphnia can protect diatoms from fungal parasitism. Limnol. Oceanogr. 49, 680–685. doi: 10.4319/lo.2004.49.3.0680
Kilias, E. S., Junges, L., Šupraha, L., Leonard, G., Metfies, K., and Richards, T. A. (2020). Chytrid fungi distribution and co-occurrence with diatoms correlate with sea ice melt in the Arctic Ocean. Commun. Biol. 3:183. doi: 10.1038/s42003-020-0891-7
Klawonn, I., van den Wyngaert, S., Parada, A. E., Arandia-Gorostidi, N., Whitehouse, M. J., Grossart, H.-P., et al. (2021b). Characterizing the “fungal shunt”: parasitic fungi on diatoms affect carbon flow and bacterial communities in aquatic microbial food webs. Proc. Natl. Acad. Sci. U.S.A. 118:e2102225118. doi: 10.1073/PNAS.2102225118
Klawonn, I., Dunker, S., Kagami, M., Grossart, H. P., and van den Wyngaert, S. (2021a). Intercomparison of two fluorescent dyes to visualize parasitic fungi (Chytridiomycota) on Phytoplankton. Microb. Ecol. doi: 10.1007/s00248-021-01893-7
Lahti, L., and Shetty, S. (2017). Tools for Microbiome Analysis in R. Version. Available online at: https://microbiome.github.io/ (accessed January 20, 2022).
Langfelder, P., and Horvath, S. (2008). WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics 9:559. doi: 10.1186/1471-2105-9-559
Lee, J. R., Raymond, B., Bracegirdle, T. J., Chadès, I., Fuller, R. A., Shaw, J. D., et al. (2017). Climate change drives expansion of Antarctic ice-free habitat. Nature 547, 49–54. doi: 10.1038/nature22996
Li, Y., Steenwyk, J. L., Chang, Y., Wang, Y., James, T. Y., Stajich, J. E., et al. (2021). A genome-scale phylogeny of the kingdom Fungi. Curr. Biol. 31, 1653–1665.e5. doi: 10.1016/j.cub.2021.01.074
Lozupone, C. A., and Klein, D. A. (2002). Molecular and cultural assessment of chytrid and Spizellomyces populations in grassland soils. Mycologia 94, 411–420. doi: 10.1080/15572536.2003.11833206
Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.J. 17:10. doi: 10.14806/ej.17.1.200
Martinez Arbizu, P. (2020). pairwiseAdonis: Pairwise Multilevel Comparison Using Adonis. R Package Version 0.4.
Matsuoka, K., Skoglund, A., Roth, G., de Pomereu, J., Griffiths, H., Headland, R., et al. (2021). Quantarctica, an integrated mapping environment for Antarctica, the Southern Ocean, and sub-Antarctic islands. Environ. Modell. Softw. 140:105015. doi: 10.1016/j.envsoft.2021.105015
McMurdie, P. J., and Holmes, S. (2013). Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One 8:e61217. doi: 10.1371/journal.pone.0061217
Middelburg, J. J., Barranguet, C., Boschker, H. T. S., Herman, P. M. J., Moens, T., and Heip, C. H. R. (2000). The fate of intertidal microphytobenthos carbon: an in situ 13C-labeling study. Limnol. Oceanogr. 45, 1224–1234. doi: 10.4319/LO.2000.45.6.1224
Monchy, S., Sanciu, G., Jobard, M., Rasconi, S., Gerphagnon, M., Chabé, M., et al. (2011). Exploring and quantifying fungal diversity in freshwater lake ecosystems using rDNA cloning/sequencing and SSU tag pyrosequencing. Environ. Microbiol. 13, 1433–1453. doi: 10.1111/j.1462-2920.2011.02444.x
Moon, H. W., Wan Hussin, W. M. R., Kim, H. C., and Ahn, I. Y. (2015). The impacts of climate change on Antarctic nearshore mega-epifaunal benthic assemblages in a glacial fjord on King George Island: responses and implications. Ecol. Indic. 57, 280–292. doi: 10.1016/J.ECOLIND.2015.04.031
Nercessian, O., Noyes, E., Kalyuzhnaya, M. G., Lidstrom, M. E., and Chistoserdova, L. (2005). Bacterial populations active in metabolism of C1 compounds in the sediment of Lake Washington, a freshwater lake. Appl. Environ. Microbiol. 71, 6885–6899. doi: 10.1128/AEM.71.11.6885-6899.2005
Oakes, J. M., Connolly, R. M., and Revill, A. T. (2010). Isotope enrichment in mangrove forests separates microphytobenthos and detritus as carbon sources for animals. Limnol. Oceanogr. 55, 393–402. doi: 10.4319/LO.2010.55.1.0393
Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., et al. (2019) Vegan: Community Ecology Package, R package version 2.5-6. Available online at: https://CRAN.R-project.org/package=vegan (accessed January 20, 2022).
Parkinson, C. L., and DiGirolamo, N. E. (2021). Sea ice extents continue to set new records: Arctic, Antarctic, and global results. Remote Sens. Environ. 267:112753. doi: 10.1016/J.RSE.2021.112753
Peay, K. G., Kennedy, P. G., and Talbot, J. M. (2016). Dimensions of biodiversity in the Earth mycobiome. Nat. Rev. Microbiol. 14, 434–447. doi: 10.1038/NRMICRO.2016.59
Priest, T., Fuchs, B., Amann, R., and Reich, M. (2021). Diversity and biomass dynamics of unicellular marine fungi during a spring phytoplankton bloom. Environ. Microbiol. 23, 448–463. doi: 10.1111/1462-2920.15331
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/NAR/GKS1219
Rasconi, S., Jobard, M., Jouve, L., and Sime-Ngando, T. (2009). Use of calcofluor white for detection, identification, and quantification of phytoplanktonic fungal parasites. Appl. Environ. Microbiol. 75, 2545–2553. doi: 10.1128/AEM.02211-08
Rasconi, S., Jobard, M., and Sime-Ngando, T. (2011). Parasitic fungi of phytoplankton: ecological roles and implications for microbial food webs. Aquat. Microb. Ecol. 62, 123–137. doi: 10.3354/AME01448
Rasconi, S., Niquil, N., and Sime-Ngando, T. (2012). Phytoplankton chytridiomycosis: community structure and infectivity of fungal parasites in aquatic ecosystems. Environ. Microbiol. 14, 2151–2170. doi: 10.1111/J.1462-2920.2011.02690.X
Rasconi, S., Ptacnik, R., Danner, S., van den Wyngaert, S., Rohrlack, T., Pilecky, M., et al. (2020). Parasitic chytrids upgrade and convey primary produced carbon during inedible algae proliferation. Protist 171:125768. doi: 10.1016/j.protis.2020.125768
Richards, T. A., Leonard, G., Mahé, F., del Campo, J., Romac, S., Jones, M. D. M., et al. (2021). Molecular diversity and distribution of marine fungi across 130 European environmental samples. Proc. R. Soc. B 282:20152243. doi: 10.1098/rspb.2015.2243
Rocha Francelino, M., Schaefer, C., de Los Milagros Skansi, M., Colwell, S., Bromwich, D. H., Jones, P., et al. (2021). WMO evaluation of two extreme high temperatures occurring in February 2020 for the Antarctic Peninsula region. Bull. Am. Meteorol. Soci. 102, 1–20. doi: 10.1175/BAMS-D-21-0040.1
Rojas-Jimenez, K., Wurzbacher, C., Bourne, E. C., Chiuchiolo, A., Priscu, J. C., and Grossart, H. P. (2017). Early diverging lineages within Cryptomycota and Chytridiomycota dominate the fungal communities in ice-covered lakes of the McMurdo Dry Valleys, Antarctica. Sci. Rep. 7:15348. doi: 10.1038/s41598-017-15598-w
Schloss, I., Ferreyra, G., González, Ó, Atencio, A., Fuentes, V., Tosonotto, G., et al. (2008). Long term hydrographic conditions and climate trends in Potter Cove. The Potter Cove coastal ecosystem, Antarctica. Synopsis of research performed 1999-2006 at the Dallmann Laboratory and Jubany Station, King George Island (Isla 25 de Mayo). Ber. Polarforsch. Meeres.
Schmitt, C. G., All, J. D., Schwarz, J. P., Arnott, W. P., Cole, R. J., Lapham, E., et al. (2015). Measurements of light-absorbing particles on the glaciers in the Cordillera Blanca, Peru. Cryosphere 9, 331–340. doi: 10.5194/TC-9-331-2015
Scholz, B., Guillou, L., Marano, A. V., Neuhauser, S., Sullivan, B. K., Karsten, U., et al. (2016a). Zoosporic parasites infecting marine diatoms — A black box that needs to be opened. Fungal Ecol. 19:59. doi: 10.1016/J.FUNECO.2015.09.002
Scholz, B., Küpper, F. C., Vyverman, W., and Karsten, U. (2016b). Effects of eukaryotic pathogens (Chytridiomycota and Oomycota) on marine benthic diatom communities in the Solthörn tidal flat (southern North Sea, Germany). Eur. J. Phycol. 51, 253–269. doi: 10.1080/09670262.2015.1134814
Scholz, B., Küpper, F. C., Vyverman, W., and Karsten, U. (2014). Eukaryotic pathogens (Chytridiomycota and Oomycota) infecting marine microphytobenthic diatoms-a methodological comparison. J. Phycol. 50, 1009–1019. doi: 10.1111/jpy.12230
Seto, K., Kagami, M., and Degawa, Y. (2017). Phylogenetic position of parasitic chytrids on diatoms: characterization of a novel clade in Chytridiomycota. J. Eukaryot. Microbiol. 64, 383–393. doi: 10.1111/jeu.12373
Shearer, C. A., Descals, E., Kohlmeyer, B., Kohlmeyer, J., Marvanová, L., Padgett, D., et al. (2007). Fungal biodiversity in aquatic habitats. Biodivers. Conserv. 16, 49–67. doi: 10.1007/s10531-006-9120-z
Sime-Ngando, T., and Twiss, M. R. (2012). Phytoplankton chytridiomycosis: fungal parasites of phytoplankton and their imprints on the food web dynamics. Front. Microbiol. 3:361. doi: 10.3389/fmicb.2012.00361
Sparrow, F. K. (1933). The Monoblepharidales. Available online at: https://about.jstor.org/terms (accessed May 26, 2021).
Stammerjohn, S. E., Martinson, D. G., Smith, R. C., and Iannuzzi, R. A. (2008). Sea ice in the western Antarctic Peninsula region: spatio-temporal variability from ecological and climate change perspectives. Deep Sea Res. II Top. Stud. Oceanogr. 55, 2041–2058. doi: 10.1016/j.dsr2.2008.04.026
Sundbäck, K., Miles, A., and Göransson, E. (2000). Nitrogen fluxes, denitrification and the role of microphytobenthos in microtidal shallow-water sediments: an annual study. Mar. Ecol. Prog. Ser. 200, 59–76. doi: 10.3354/MEPS200059
Taylor, J. D., and Cunliffe, M. (2016). Multi-year assessment of coastal planktonic fungi reveals environmental drivers of diversity and abundance. ISME J. 10, 2118–2128. doi: 10.1038/ismej.2016.24
QGIS Development Team (2021). QGIS Geographic Information System. Available online at: https://www.qgis.org (accessed December 7, 2021).
Tedersoo, L., Anslan, S., Bahram, M., Põlme, S., Riit, T., Liiv, I., et al. (2015). Shotgun metagenomes and multiple primer pair-barcode combinations of amplicons reveal biases in metabarcoding analyses of fungi. Mycokeys 10, 1–43. doi: 10.3897/mycokeys.10.4852
Tedersoo, L., Tooming-Klunderud, A., and Anslan, S. (2018). PacBio metabarcoding of Fungi and other eukaryotes: errors, biases and perspectives. New Phytol. 217, 1370–1385. doi: 10.1111/nph.14776
Tisthammer, K. H., Cobian, G. M., and Amend, A. S. (2016). Global biogeography of marine fungi is shaped by the environment. Fungal Ecol. 19, 39–46. doi: 10.1016/j.funeco.2015.09.003
van den Wyngaert, S., Seto, K., Rojas-Jimenez, K., Kagami, M., and Grossart, H.-P. (2017). A new parasitic chytrid, Staurastromyces oculus (Rhizophydiales, Staurastromycetaceae fam. nov.), infecting the freshwater desmid Staurastrum sp. Protist 168, 392–407. doi: 10.1016/J.PROTIS.2017.05.001
van Donk, E. (1989). The role of fungal parasites in phyto-plankton succession.. Plankton Ecol. ed. U. Sommer (Berlin: Springer), 171–194. doi: 10.1007/978-3-642-74890-5_5
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07
Wang, Y., Guo, X., Zheng, P., Zou, S., Li, G., and Gong, J. (2017). Distinct seasonality of chytrid-dominated benthic fungal communities in the neritic oceans (Bohai Sea and North Yellow Sea). Fungal Ecol. 30, 55–66. doi: 10.1016/j.funeco.2017.08.008
Watts, S. C., Ritchie, S. C., Inouye, M., and Holt, K. E. (2019). FastSpar: rapid and scalable correlation estimation for compositional data. Bioinformatics 35, 1064–1066. doi: 10.1093/BIOINFORMATICS/BTY734
Wickham, H. (2016). ggplot2: Elegant Graphics for Data Analysis. New York, NY: Springer-Verlag, doi: 10.1007/978-3-319-24277-4
Willis, A. D. (2019). Rarefaction, alpha diversity, and statistics. Front. Microbiol. 10:2407. doi: 10.3389/FMICB.2019.02407/BIBTEX
Woelfel, J., Schumann, R., Peine, F., Flohr, A., Kruss, A., Tegowski, J., et al. (2010). Microphytobenthos of Arctic Kongsfjorden (Svalbard, Norway): biomass and potential primary production along the shore line. Polar Biol. 33, 1239–1253. doi: 10.1007/s00300-010-0813-0
Wright, E. S., Yilmaz, L. S., and Noguera, D. R. (2012). DECIPHER, a search-based approach to chimera identification for 16S rRNA sequences. Appl. Environ. Microbiol. 78, 717–725. doi: 10.1128/AEM.06516-11
Wulff, A., Iken, K., Quartino, M. L., Al-Handal, A., Wiencke, C., and Clayton, M. N. (2009). Biodiversity, biogeography and zonation of marine benthic micro- and macroalgae in the Arctic and Antarctic. Bot. Mar. 52, 491–507. doi: 10.1515/BOT.2009.072
Keywords: Antarctica, aquatic fungi, Chytridiomycota, phytoplankton host, salinity gradient, Illumina amplicon sequencing, Carlini Station
Citation: Ilicic D, Woodhouse J, Karsten U, Zimmermann J, Wichard T, Quartino ML, Campana GL, Livenets A, Van den Wyngaert S and Grossart H-P (2022) Antarctic Glacial Meltwater Impacts the Diversity of Fungal Parasites Associated With Benthic Diatoms in Shallow Coastal Zones. Front. Microbiol. 13:805694. doi: 10.3389/fmicb.2022.805694
Received: 30 October 2021; Accepted: 12 January 2022;
Published: 04 March 2022.
Edited by:
David Velazquez, Autonomous University of Madrid, SpainCopyright © 2022 Ilicic, Woodhouse, Karsten, Zimmermann, Wichard, Quartino, Campana, Livenets, Van den Wyngaert and Grossart. 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: Hans-Peter Grossart, hgrossart@igb-berlin.de