Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Microbiol., 23 July 2025

Sec. Terrestrial Microbiology

Volume 16 - 2025 | https://doi.org/10.3389/fmicb.2025.1633308

This article is part of the Research TopicPlant Mineral Microbe Interactions, Vol IIView all 4 articles

Mass development of a filamentous and likely nitrophilous aerophytic green alga on tree bark: Apatococcus ammoniophilus sp. nov. (Chlorophyta, Trebouxiophyceae)

  • 1Section Ecology and Evolution, Department of Biology, University of Copenhagen, Copenhagen, Denmark
  • 2Department Experimental Phycology and Culture Collection of Algae (EPSAG), Albrecht-von-Haller-Institute for Plant Sciences, Georg August University, Göttingen, Germany
  • 3Section Marine Biology, Department of Biology, University of Copenhagen, Copenhagen, Denmark
  • 4Grainger Bioinformatics Center, Field Museum, Chicago, IL, United States
  • 5Institute of Pharmacy/Pharmacognosy, University of Innsbruck, Center for Chemistry and Biomedicine, Innsbruck, Austria
  • 6Faculty for Chemistry, Physics and Biosciences, Biology/Ecology, Freiberg, Germany
  • 7Department of Applied Ecology and Phycology, Institute of Biological Sciences, University of Rostock, Rostock, Germany
  • 8Department of Maritime Systems, Interdisciplinary Faculty, University of Rostock, Rostock, Germany

Introduction: A filamentous green alga forming significant biomass on twigs and needles was observed to have increased invasively in Denmark in recent decades. It was particularly abundant in coniferous plantations in western parts of Denmark that experience the highest modelled concentration of atmospheric nitrogen deposition. However, its species identity and taxonomy remained unknown.

Material and methods: Selected algal samples from various substrates were analyzed for their ribosomal DNA sequences, metagenomic, and biochemical compounds [polyols and mycosporine-like amino acids (MAAs)].

Results: Phylogenetic analyses revealed the alga’s position within the Trebouxiophyceae (Chlorophyta), forming an independent lineage within Apatococcus. Though it was associated with various other Trebouxiophyceae species, the metagenome showed exceptionally high coverage of the Apatococcus contigs, proving its predominance, consistent with the amplicon-based approach. The low molecular weight carbohydrates, arabitol, erythritol, and trehalose – with erythritol displaying the highest concentrations—were recovered. The presence of erythritol provided chemotaxonomic support for the classification in Apatococcus. Additionally, a unique UV-absorbing mycosporine amino acid (MAA), likely new for the Trebouxiophyceae, was found. The species is described here as A. ammoniophilus, and the observed morphological features leave no doubt that it has been recorded from Denmark more than a 100 years ago. Morphological features are shared with its closer relatives, such as the presence of a ring of particles surrounding the nucleus and the formation of two-celled units.

Discussion: The presence of low molecular weight carbohydrates and the unique MAA in A. ammoniophilus well explain the biochemical basis for its aeroterrestrial lifestyle, as these organic compounds protect against desiccation and UV-radiation, respectively. Even though the genotype of A. ammoniophilus has also been found in inconspicuous biofilms devoid of filamentous stages on various substrates with presumably low ammonia deposition, the very invasive colonization in recent decades in western Denmark is assumed to be due to ammonia deposition. Consequently, A. ammoniophilus is suggested to be a potential biological indicator of air borne nitrogen deposition. A possible connection between filamentous growth and nitrogen accumulation needs further investigation, including culture experiments.

Introduction

Eukaryotic microalgae and cyanobacteria inhabit numerous terrestrial habitats across all biogeographic regions, forming phototrophic biofilms at the interface between the atmosphere and natural or anthropogenic solid substrates. Consequently, these ecologically unique taxa are classified as aerophytic. In temperate regions of the Northern Hemisphere, eukaryotic green microalgae from the Chlorophyta and Streptophyta are the most abundant terrestrial microorganisms (Häubner et al., 2006; Rindi, 2007; Holzinger and Karsten, 2013). These aerophytic, phototrophic biofilms cover natural surfaces such as rocks, tree bark, soil, and plant surfaces, as well as anthropogenic substrates like concrete, glass, and metal (Karsten et al., 2007), often forming conspicuous communities.

Water availability is a key ecological driver for these aerophytic algae, which are poikilohydric organisms and thus cannot actively regulate their water content (Kranner et al., 2008). As a result, they can easily undergo desiccation stress under water-limited conditions, although they often display adaptive traits to survive partial or even complete water loss (Holzinger and Karsten, 2013).

Desiccation typically results in a significant decline in photosynthesis and growth, with both physiological processes potentially being completely blocked under severe stress conditions (Holzinger and Karsten, 2013). Bertsch (1966) investigated the impact of desiccation on carbon assimilation by reducing air humidity in the aeroterrestrial green alga Apatococcus lobatus. This taxon is among the most abundant green algae in temperate Europe, forming conspicuous biofilms on tree bark (Barkmann, 1958), roof tiles, and building surfaces. Bertsch (1966) demonstrated that these cell packets or biofilms maintain a hydration equilibrium with the vapor pressure of the air, showing optimum carbon assimilation at 97–98% relative air humidity (rah); carbon assimilation decreased to 50% at 90% relative humidity and ceased at 68% relative humidity.

Aerophytic microalgae have developed a variety of morphological, physiological, and biochemical adaptive traits to the harsh terrestrial environment (Holzinger and Karsten, 2013). Photosynthesis is a key physiological process and hence efficient control of light absorption and energy distribution in the photosynthetic apparatus during dehydration is crucial to minimize or prevent photoinhibition. The biochemical capability to synthesize and accumulate organic osmolytes is regarded as a fundamental biochemical mechanism to provide desiccation tolerance (Holzinger and Karsten, 2013). By accumulating such compounds, aerophytic microalgae maintain turgor, membrane integrity, and macromolecule structure by compensating for changes in water potential without incurring metabolic damage (Yancey, 2005; and references therein). Many aerophytic microalgae, such as Apatococcus, Diplosphaera, Stichococcus, and other members of the green algal class Trebouxiophyceae (Chlorophyta), synthesize and accumulate polyols like glycerol, arabitol, ribitol, mannitol, and sorbitol (Gustavs et al., 2011, 2016; Hotter et al., 2018; Medwed et al., 2021). Polyols can also act as antioxidants or heat protectants, leading to the stabilization of proteins (Karsten et al., 2007); due to their multiple functions, polyols serve as one of the biochemical explanations for compensating terrestrial stress conditions.

Besides desiccation, aerophytic microalgae are regularly confronted with solar ultraviolet radiation (UVR), which exerts many harmful effects on cells and can even cause DNA damage (Karsten, 2008). Hence, many aerophytic microalgae are capable of synthesizing and accumulating UV-absorbing mycosporine-like amino acids (MAAs), which are low-molecular weight compounds with numerous chemical variations in their side groups and substituents, resulting in different absorption maxima between 310 and 360 nm (Karsten, 2008; Hartmann et al., 2016, 2020; Hotter et al., 2018). These colorless sunscreens shield the algal cells from UVR by absorbing the harmful radiation energy and converting it into harmless heat (Bandaranayake, 1998).

In temperate regions, the bark of tree trunks is most often covered by lichens or green algal biofilms which are dominated by members of Trebouxiophyceae genera, e.g., Apatococcus, Desmococcus, Diplosphaera, and Trebouxia (Medwed et al., 2021; Karsten et al., 2022), as well as Trentepohlia (Ulvophyceae, Chlorophyta) (Holzinger et al., 2023).

In Denmark, the formation of pulverulent algal biofilms on the bark of trees is extremely common. However, in the 1980s, a filamentous alga suddenly became much more frequent, particularly in Jutland (Søchting, 1997) (Figure 1). This alga occupied the trunks, branches, dead twigs, and needles of conifers including Christmas trees, causing aesthetically unacceptable discoloration, which led to economic problems in selling these products. The occurrence of a filamentous green alga on the twigs of spruce in Denmark was already reported by Petersen (1915) in his studies on Danish ‘aerophilic’ algae. He identified it as Hormidium crenulatum and provided a detailed description and illustration of the species. He also observed that H. crenulatum was growing only where ammonia influence from farming was present, and its substantial expansion during the 1980s was suggested to result from increased ammonia deposition from husbandry farming (Søchting, 1997). In recent years, this alga has been recorded from all parts of Denmark, and in agricultural regions of Jutland it forms almost monospecific filamentous mats on coniferous and deciduous bark, as well as on spruce and fir needles and on dwarf shrubs in heathlands (Figure 2). H. crenulatum was transferred to Klebsormidium and lectotypified by Lokhorst and Star (1985) and later epitypified by Mikhailyuk et al. (2015). However, the morphology of the epitype of K. crenulatum differs from the Petersen (1915) alga (see Discussion below), raising serious doubts about whether the latter could be a member of Klebsormidiophyceae (Streptophyta). It follows that this abundant aerophytic filamentous green alga represents a species for which no name is yet available.

Figure 1
Map of Denmark indicating ammonia concentration levels over land areas, with a color gradient from blue (low concentration) to purple (high concentration). Symbols indicate sample locations: circles show type, diamonds show sequenced samples, and black diamonds denote metagenomics samples. North Sea and Baltic Sea label positions.

Figure 1. Locations from which the samples of the filamentous Apatococcus ammoniophilus have been collected (2018–2024; Table 1) superimposed on the map of Denmark with calculated annual mean values of ammonia-N concentration in the air in 2020 of Ellermann et al. (2021).

Figure 2
Images of tree branches and bark covered with filamentous green algae. Panels a, b, and c show close-ups of branches with green growth. Panel d highlights a vertical section with texture detail. Panel e displays a microscopic view of the surface structure. Panel f features branches with conifer needles.

Figure 2. Apatococcus ammoniophilus biofilms on various tree substrates in Denmark. (a,b) On the trunk of Pinus mugo, Kikkebjerg, Fanø (voucher US13007). (c,d) On the twig of Pinus sp., Råbjerg Hede, Vendsyssel, Jutland (voucher US13044). (e) On the bark of Pinus sp., Mosevrå Kirke, Jutland, SEM (voucher US12952). (f) On the twigs and needles of Pinus sp., Løbners Plantage, Jutland.

We studied a selection of environmental samples of the filamentous alga obtained from various substrates to analyze their ribosomal DNA sequences and perform a metagenomic community profiling, aiming to uncover the alga’s phylogenetic position. Additionally, we examined the samples for specific biochemical compounds, i.e., low molecular weight carbohydrates (polyols) and mycosporine amino acids (MAAs), which may form important adaptive traits for the alga’s subaerial lifestyle.

Materials and methods

Sampling

This study is based primarily on intensive collecting of bark, twigs, branches and needles of mostly pine trees in Denmark, with addition of material from Germany, Norway, Sweden, England and the Netherlands. Studied specimens are listed in Table 1. Vouchers of all collections are deposited in the Natural History Museum of University of Copenhagen (C). Samples for molecular and biochemical analyses were air dried and stored in the dark at ambient room temperature.

Table 1
www.frontiersin.org

Table 1. List of the studied samples of Apatococcus ammoniophilus.

Microscopy and attempts for culture

Microscopic observations were made using an Olympus BX60 microscope with Nomarski DIC optics; micrographs were taken with a Gryphax camera (Jenoptik) and processed with the corresponding Gryphax software (Jenoptik). Scanning Electron Microscopy was performed on freshly dried material from nature (US12952) that was attached to double-sided adhesive tape on a SEM stub and coated with gold–palladium. It was examined at 7 kV in a Jeol (Tokyo, Japan) JEM JSM 6335F field emission scanning electron microscope at the Zoological Museum, University of Copenhagen.

In an attempt to culture the filamentous alga, ten environmental samples were selected for culturing. Only one sample, US13013 from a Pinus twig, succeeded in growing, but the culture stopped growing after about 6 weeks. All further attempts to establish a stable culture failed. The US13013 material was kept in a dark and dry condition for about 1 week after collection. Material was scraped off the substrate with a sterile scalpel and carefully distributed on the surface of agarized (1.5%) Bold Basal culture medium with triple nitrate and vitamins (3NBBM + V, Available at: www.epsag.uni-goettingen.de; accessed 04 December 2024), a modification of the Bold 3 N medium of Starr and Zeikus (1993). The initial culturing was performed on agar plates maintained at 20°C under a light:dark (L: D) regime of 14:10 h and a photosynthetic photon flux density (PPFD) of 23 μmol m−2 s−1.

DNA extraction, PCR, and sequencing

Algal material was scraped off the substrate (wooden twigs) into a sterile petri dish using sterilized scalpels and pooled into 100 μL reaction tubes containing the lysis buffer of the Spin Plant Mini Kit (Invisorb). The cells were carefully mechanically broken at 5,000 rpm twice for 30 s in the presence of equivalent amounts of 100–200 μm and 425–600 μm diameter glass beads (Sigma-Aldrich) using a PowerLyzer® 24 instrument (MO BIO Laboratories, Inc., Carlsbad, California, USA). Genomic DNA was extracted with the same kit. DNA concentration was measured with an Invitrogen™ Qubit™ 3 fluorometer (Fisher Scientific). About 9–14 ng of the extracted DNA was used as a template for PCR amplifications in reaction mixtures of 50 μL volume. The 18S rRNA gene region was amplified with primer pair NS1 (White et al., 1990) and 18 L (Hamby et al., 1988) or NS1 and 1800R the reverse complement of primer 1800F (Friedl, 1996). For amplification of the ITS1–5.8S– ITS2 region with adjacent segments of the 18S and 26S rRNA gene regions served the primer pairs AL1500af (Helms et al., 2001) and LR3 (Vilgalys and Hester, 1990), as modified by Friedl and Rokitta (1997). PCR was performed with MyTaq™ DNA Polymerase (Meridian Bioscience) in a thermocycler Biometra TProfessional basic. The initial denaturation was at 95°C for 5 min, followed by 35 cycles of denaturation at 94°C for 45 s, annealing at 50°C for 45 s, extension at 72°C for 90 s, and a final extension at 72°C for 10 min. Amplicons were checked by agarose (1.5%) electrophoresis and purified with the MSB® Spin PCRapace Kit (Invitek). The PCR products were sequenced at Microsynth Seqlab GmbH (Göttingen, Germany) using the PCR primers. The sequences were edited, and contigs assembled using the DNA sequence contig assembler software DNA Dragon (SequentiX - Digital DNA Processing, Klein Raden, Germany; Available at: www.dna-dragon.com accessed 04 December 2024). The 18S, ITS1, 5.8S, ITS2, and LSU regions were extracted from the long sequences with ITSx version 1.1b (Bengtsson-Palme et al., 2013; Rivers et al., 2018). The new sequences were compared with a cloned sequence, DZK5, which was obtained from the surface of a roof tile of a private house in a residential area within the city of Göttingen, Germany (51.558588 N, 9.939069E). The long sequences are available from the DDBJ/EMBL/GenBank databases under accession numbers PQ763401 - PQ763403.

Paired-end ITS2 metabarcoding and sequence processing

For amplicon sequencing, the ITS2 rDNA region was targeted using the primer pairs ITS3-KYO2 (Toju et al., 2012) and ITS4 (White et al., 1990), which allow for the amplification of green algae and fungi. About 9–14 ng of the extracted DNA was used as a template for PCR amplifications in reaction mixtures of 50 μL volume. The PCR with barcoded primers [P7/P5 indices from Kircher et al. (2012)] and the PhusionTM High–Fidelity DNA Polymerase (Thermo Scientific) was done using the following program in a Biometra TProfessional basic thermal cycler: a first denaturation step at 98°C for 1 min followed by 25 cycles at 98°C for 45 s, annealing at 48°C for 45 s, and extension at 72°C for 30 s, and final extension at 72°C for 5 min. To capture and purify DNA fragments in the size range of 300–600 base pairs, MagSi-NGSPREP Plus beads (Magtivio; Nuth, Netherlands) on an Ambion Invitrogen™ Magnetic Stand (Fisher Scientific), with a ratio of reagent/sample of 0.8, were used. Next, the samples were pooled, and the pools DNA concentration was adjusted to the requirements of the paired-end sequencing (≥ 50 ng/μL) using ethanol precipitation in the presence of ammonium acetate (Zeugin and Hartley, 1985). Illumina paired-end sequencing was performed at Novogene GmbH (Martinsried, Germany) on a Novaseq6000 platform on SP flow cell with v1.5 (500 cycles) chemistry (2 × 250 pair-end reads).

The high-throughput sequencing yielded 5,072,212 paired-end reads, which were demultiplexed using cutadapt [Version 4.9, (Martin, 2011)] with combinatorial dual indexing. We allowed a maximum error rate of 0.15 and a minimum overlap of 4 bases for index matching, resulting in 4,033,349 assigned reads. A more detailed description of the demultiplexing procedure can be found in the Supplementary file S1. Amplicon sequence variants (ASVs) were identified using nf-core/ampliseq [Version 2.11.0; (Straub et al., 2020)], a standardized bioinformatics pipeline implemented in Nextflow (Di Tommaso et al., 2017) within the nf-core (Ewels et al., 2020) framework. Primers were trimmed using cutadapt, and untrimmed sequences were discarded. Primer-free sequences were processed as one pool with DADA2 (Callahan et al., 2016) to remove PhiX contamination, truncate sequences (forward reads and reverse reads each at 220 bp, sequences shorter than this were discarded), discard sequences with more than two expected errors, merge read pairs and remove PCR chimeras. A total of 4,307 ASVs was obtained across all samples. Between 67.55 and 87.56% of reads per sample (average 82.4%) were retained. The ASV count table contained 3,241,892 counts, at least 75,953 and 746,194 per sample (average 294,717). Taxonomic assignment of ASVs was conducted using BLASTN [Version 2.14.0; (Altschul et al., 1990)] against the GenBank Nucleotide database (NCBI-GenBank Flat File Release 259.0 of December 15, 2023) using default parameters (Altschul et al., 1997). Taxonomic labels were assigned to ASVs using the consensus-based protocol of Rybalka et al. (2023). To ensure accuracy, we validated the taxonomic assignments through manual examination of consensus taxonomy, using normalized BLAST bitscore (NB) values as a quantitative criterion for acceptance (Rybalka et al., 2023). ASVs abundance data were filtered using R (Version 4.3.1, R Core Team, 2024) with phyloseseq (Version 1.46.0; (McMurdie and Holmes, 2014)) and tidyverse (Version 2.0.0; (Wickham et al., 2019)). After converting read counts to relative abundances, ASVs were retained if they exceeded a 0.1% relative abundance threshold in at least one sample (236 ASVs). A subset of 59 ASVs were clustered at a 98% sequence similarity threshold using VSEARCH (v2.29.1) (Rognes et al., 2016) after sequence subset extraction with seqkit [Version 2.9.0; (Shen et al., 2016)], with the most abundant ASV selected as representative for multi-ASV clusters and abundances aggregated within each cluster. A heatmap, computed with tidyverse and phyloseq in R (McMurdie and Holmes, 2013), to show the relative abundances of clustered ASVs. Colors represent relative abundance values transformed using the shrinkage function.

f ( x ) = x + ε 1 + 2 ε ( ε = 10 5 ) followed by a quantile logistic transformation to enhance the visibility of low-abundance values.

Of each of the 34 multi-ASV clusters, one representative, as listed in the caption of the heatmap (Figure 3), is available from the DDBJ/EMBL/GenBank databases under the accession numbers PQ784906 - PQ784939.

Figure 3
Heatmap showing the abundance of various ASVs (Amplicon Sequence Variants) across different samples labeled from US12921 to agar sample. Darker blue indicates higher abundance, with a scale ranging from 0.00 to 0.75. Labels on the left include Apatococcus, Chloridium, and Trebouxia among others. Percentage values at the top range from 61.0% to 98.8%.

Figure 3. Heat map of relative abundances of the 34 clusters of green algal ASVs recovered by the ITS2 amplicon-based metabarcoding. It shows their distribution across the ten environmental samples and the subculture of sample US13013 on an agar plate (agar sample). The numbers on top are the relative abundances of ASV_00001 out of all algal reads per sample. An ASV ID is with its genus identification and the normalized score of pairing significance to its closest reference sequences (NB); (Rybalka et al., 2023). Scale, color brightness within the matrix indicates the relative abundance of sequence reads. T, Trebouxiophyceae (Chlorophyta); K, Klebsormidiophyceae (Streptophyta).

ITS2 phylogenetic analyses

The ITS2 region of all selected ASVs and reference sequences was annotated using the ITS2 database (Keller et al., 2009). ITS2 reference sequences were selected from (Zahradníková et al., 2017) and from the top BLAST hit for each ASV. Secondary structures of ITS2 sequences were predicted with RNAfold (Available at http://rna.tbi.univie.ac.at accessed Dec 4, 2024; Gruber et al., 2015) setting the temperature at 25°C and using minimum free energy (MFE) criteria. ITS2 sequences, including their secondary structures, were aligned with 4SALE 1.7.1 (Wolf et al., 2014). The maximum likelihood (ML) tree was estimated by 1,000 bootstrap replicates with the R script available at http://4sale.bioapps.biozentrum.uni-wuerzburg.de/mlseqstr.html in phangorn (Schliep, 2011) using one-letter encoded sequence-structure data exported from 4SALE (Wolf et al., 2014). Nine alternative topologies were constructed using the Interactive Tree of Life (iTOL) Version 6; Available at https://itol.embl.de accessed Dec 4, 2024; (Letunic and Bork, 2024) by repositioning the ASV_00001 lineage in different arrangements relative to the ML tree referred to as the “best tree.” A Shimodaira-Hasegawa (SH) test was conducted using the SpeciesTopoTestR (Adams and DeGiorgio, 2023) in R 4.4.1 (R Core Team, 2024), comparing each alternative topology against the best tree. Almost full 18S sequences were first aligned with a large set of reference sequences (143) that represented major lineages of the Trebouxiophyceae (Chlorophyta) using MAFFT (Katoh and Standley, 2013) and then curated, i.e., selecting regions in the multiple sequence alignment that are suited for phylogenetic inference, with BMGE (Criscuolo and Gribaldo, 2010) on the NGPhylogeny.fr web portal (Lemoine et al., 2019). The set of references included 26 different sequences longer than 1,000 base pairs representing the various lineages of the genus Apatococcus. The final dataset was 1747 nucleotides (columns) long with 753/545 variable/parsimony informative sites. A maximum likelihood phylogenetic tree was generated using IQ-TREE 2 (Minh et al., 2020). Branch support was assessed with the ultrafast bootstrap approximation [1,000 replicates, (Hoang et al., 2017)] implemented in IQ-TREE. IQ-TREE selected the TNe + R4 model to describe the best substitution pattern. The final tree was visualized and rooted using iTOL [Version 6; (Letunic and Bork, 2024)].

Metagenomic approach

For the metagenomic analyses, algal material from a Larix twig (US13019) was scraped off to obtain a powdery tissue from which total DNA was extracted using the Qiagen DNeasy Plant Mini Kit (Qiagen, Germantown). At the Research Resource Center at the University of Illinois, Chicago, the DNA was size-selected to target fragments with an insert size of 200 bp using a Pippin Prep Size Selection system (Sage Science, Beverly). A bead cleanup subsequently followed this process. The selected fragments were then used to create a sequencing library using the Nextera XT kit (Illumina, San Diego) before sequencing approximately 10 Gbp of 150 bp paired-end reads on an Illumina HiSeq4000.

To gather initial genomic information, we assembled all raw reads with SPAdes 3.11.1 (Prjibelski et al., 2020) using the default settings. All reads were inspected by using the program FastQC v0.11.3.1 Further, to improve the metagenomic assembly, we trimmed all reads with Trimmomatic v0.33 (Bolger et al., 2014), setting a quality threshold of 10 (LEADING:10 TRAILING:10) and a minimum read length of 25 bp (MINLEN:25). The surviving paired-end reads were used for another assembly with the metaSPAdes (Nurk et al., 2017) assembler implemented in SPAdes v3.13.0 and designed to perform better on metagenome assemblies. We used Quast v5.2.0 (Gurevich et al., 2013) to compare the quality of the two assemblies resulting from SPAdes and metaSPAdes. To receive coverage information for each scaffold in addition to the k-mer coverage provided by SPAdes, we used minimap2 v2.28-r1209 (Li, 2018) optimized for short reads (−ax sr -L --MD), and samtools v1.13 (Li et al., 2009) for sequence mapping, sorting, and counting of sequence coverage. All metagenomic raw reads were uploaded to NCBI with the accession number SRR32987805 under BioProject PRJNA1246228 and Biosample SAMN47770417.

We identified all contigs in the assembly that contain a 28S, 18S, 5.8S, and 5S ribosomal gene using the program barrnap v0.8 with the built-in “Eukaryota” database.2 In addition, we used the program SSU-ALIGN (Nawrocki, 2009), which identifies all contigs that contain the 18S ribosomal gene. We used all identified gene regions for a reciprocal megaBLAST (Altschul et al., 1990) search against NCBI’s non-redundant nucleotide (nt) database with Geneious Prime v2024.0. We then identified and retained all contigs classified as “Trebouxiophyceae” or “Klebsormidiophyceae” along with the top three results from each search for further review of each scaffold. To visualize the coverage and integrity of specific genomic regions, we mapped all metagenomic NGS reads back to the reference sequences using Geneious Prime v2025.0.3. Mapping was performed with Custom Sensitivity settings to ensure the highest possible stringency. The following parameters were applied: perfect matches only were allowed by setting the maximum mismatches per read to 0%, the minimum overlap identity to 100%, and disallowing gaps in the alignment. Ambiguous bases were limited to a maximum of 1. Reads were trimmed for paired-end overhangs, and the minimum mapping quality was set to 100 to include only high-confidence alignments. The mapping process did not include iterative reference extension or fine-tuning, focusing on fast and precise alignment.

Polyol analysis

For polyol extraction, 10–15 mg of the air-dried biofilm material was carefully scratched off various bark samples from Pinus species in Jutland, Denmark (US12923, 12,924, 12,925, 12,927). The resulting green powder was extracted in 70% aqueous ethanol (v/v) at 70°C for 4 h. After centrifugation at 8,300 × g for 5 min, 800 μL of the supernatant was evaporated to dryness and dissolved in 800 μL of autoclaved ultrapure water. After re-dissolution, vortexing and further centrifugation steps, the supernatant was stored in a 2 mL screw cap glass vial with silicone/PTFE septum (Wicom, Heppenheim, Germany) until high performance liquid chromatography (HPLC) analysis. HPLC was performed on an Agilent 1,260 Infinity HPLC system with RI-detector (Agilent Technologies, Santa Clara, CA, USA). Fast Carbohydrate Column (Bio-Rad, Feldkirchen, Germany) with a Carbo-Pb2+ guard column (Phenomenex, Aschaffenburg, Germany) was used for separation at 70°C and a flow rate of 1.0 mL min−1 at a pressure of 45 bar (Karsten et al., 1991). All polyols for calibration were purchased from Carl Roth (Karlsruhe, Germany). Polyol concentrations were calculated as μg per mg dry weight.

MAA analysis and structural elucidation

MS-grade acetonitrile was acquired from Merck (Darmstadt, Germany) and ammonium formate from Serva (SERVA Electrophoresis GmbH, Heidelberg, Germany). MS-grade formic acid was purchased from VWR International (Vienna, Austria). Ultra-pure water was prepared in-house with a Sartorius Arium 611 UV purification system (Sartorius AG, Göttingen, Germany). Standards of MAAs with confirmed chemical structure and purity were available from previous projects (Orfanoudaki et al., 2020; Zwerger and Ganzera, 2022). The dried algal biofilms were scratched off numerous bark and twig samples from different locations to reach sufficient biomass; they were pooled and milled using a rotating blade coffee grinder from Bosch (Munich, Germany) and sieved through a 710 μm steel mesh from Retsch (S/N 21241665, Body: 200 mm × 50 mm; Haan, Germany). An aliquot of the ground sample (m = 544.98 mg) was then extracted with HPLC-grade water (V = 10 mL) using ultrasonication (t = 15 min). Thereafter, the sample was centrifuged (RCF = 2,700 x g, t = 5 min) and the supernatant filtered through a paper filter with a particle retention of 10–20 μm from VWR International (Vienna, Austria). This extraction procedure was repeated four more times, and the filtered supernatants were combined. Using rotary evaporation under reduced pressure (T water bath = 40°C), the combined aqueous extract was dried. Any solvent residues were removed by freeze-drying to finally obtain 29.5 mg (5.4% dry weight) of dried extract. HPLC-DAD analyses were performed on an Agilent Technologies 1,260 Infinity II system equipped with a quaternary pump, vial sampler, column thermostat, and diode-array detector (Agilent Technologies, Santa Clara, USA). Separation conditions were according to the method described by Orfanoudaki et al. (2020), but the column oven was set to 12°C. The UHPLC analysis of the aqueous extract (c = 5 mg mL−1, solvent for dissolution = HPLC-grade water), a MAA standard mix and several MAA-containing fractions available from previous projects was performed on a Vanquish system (Thermo Scientific, Waltham, MA, USA) consisting of a quaternary pump, an auto-sampler, a column oven, and a variable wavelength detector connected to a Thermo Scientific Exploris 120 Orbitrap HRMS unit. Separation was carried out on a Phenomenex Luna Omega C18 100 Å column (100 mm × 2.1 mm; particle size = 1.6 μm) protected by a Phenomenex SecurityGuard ULTRA guard cartridge system (i.e., a UHPLC C18 pre-column). The mobile phase comprised water with 0.25% formic acid and 20 mM ammonium formate (A) and acetonitrile (B). The applied gradient was as follows: 0 min, 0% B; 10 min 0% B; 11 min, 90% B; 13 min, 90% B; 13.1 min, 0% B. Finally, the column was re-equilibrated with the original solvent composition (i.e., 0% B) for 19.9 min., which corresponds to a total run time of 33 min. The flow rate, column oven temperature, auto-sampler temperature, and injection volume were adjusted to 0.3 mL min−1, 17°C, 20°C, and 1 μL, respectively. The detection wavelengths were set to 330 and 350 nm, the data collection rate to 2.0 Hz, the response time to 2.0 s, and the peak with to 0.2 min. The system was controlled by Thermo Scientific Xcalibur 4.4 software. Calibration of the mass analyzer was done via the Thermo Scientific proprietary calibration mix and the respective automatic calibration function. The mass spectrometric parameters were as follows: heated-ESI ionization source, static spray voltage (positive and negative: 3500 V), sheath gas (N2): 30 arbitrary units, auxiliary gas (N2): 17 arbitrary units, sweep gas (N2): 0 arbitrary units. Temperature of the ion transfer tube and vaporizer was adjusted to 370 and 420°C, respectively. MS data (range 70–1,000 m/z) were recorded from 0 to 13 min with a resolution of 60,000 FWHM for MS1. The RF lens parameter was set to 70%. Data-dependent experiments were conducted with stepped collision energy mode and normalized collision energy type using HCD collision energies of 15, 30, and 45% at a resolution of 15,000 FWHM. The number of dependent scans was set to 3. The following selection of filters was employed: intensity threshold filter (1.0E5), dynamic exclusion (parent ions were placed in the exclusion list for 2 s after detection), isotope exclusion, charge state (perform dependent scans on singly charged precursors only), and apex filter (desired apex window: 75%). In addition, a specific exclusion list was created for the measurement using ultra-pure water as a background extract with an IODA Mass Spec notebook (Zuo et al., 2021).

Results

Distribution and environmental data

The filamentous alga was collected from living or dead bark of coniferous and deciduous trees and shrubs throughout Denmark (Figure 1), as well as from some locations in England, Sweden and the Netherlands (Table 1). Furthermore, it was collected from spruce and fir needles (Søchting et al., 1992), and from dwarf shrubs in highly ammonia-exposed habitats. The alga developed considerable filamentous biofilms (Figure 2), particularly in ammonia-producing farmland, which in Denmark is concentrated in the western part of the country (Figure 1). On twigs, it formed a continuous, light green biofilm on the upper side, while dominating the NNW-exposed sides of standing trunks. Generally, the alga formed a monoculture without being intermixed with lichens. Microscopic observation revealed a few cells of additional coccoid green algae, which were never dominant in the algal biofilms and comprised approximately one to 5 % of the algal community.

Morphology

In the uncultured environmental material, the filaments appeared to be formed by units containing two cells each (Figure 4). New walls encircled the two daughter cells formed within the same mother cell wall (Figures 4dh). We interpret these observations as units that are expanding through cell growth; the old (mother) cell wall stretches and eventually ruptures, leaving cap-like remnants at both ends of a unit (e.g., Figures 4c,d,f,g). Small remnants of the ruptured old wall also persisted on the surface of the new wall (Figures 4e,g). The two-celled units remained attached at both ends, forming uniseriate filaments (Figures 4eh). The filaments broke off between the units (Figures 4c,h), facilitating the reproduction of the alga. The cells contained a flat chloroplast with several lobes (Figures 4dh). Droplets of storage products were found between the chloroplast and the cell wall (Figures 4eg). Remarkably, many cells in the environmental material were not viable (Figures 4a,b); they exhibited damaged chloroplasts, protoplasts detached from the cell wall (Figure 4b), pale remnants of protoplasts within cells, and, ultimately, empty cell walls (Figures 4a,b). In parts of the filaments with empty cells, the walls appeared particularly thick, sometimes even swollen (Figure 4b); likely, they degraded very slowly. In culture, filaments consisting of units with two cells (autospores) each were observed (Figures 5ag). Each two-celled unit may represent an autosporangium, and through their adhesion to each other at their ends, uniseriate filaments may be formed. The filaments easily disintegrated into shorter fragments comprising a few two-celled units (Figures 5d,e). Cap-like remnants of previous walls remained at the ends of the units and were kept distant by the growing filament. In the growing cultures, a tendency to form cell packages was also observed (Figures 5d,f,g). At some points in the filaments, the division plane shifted to the upright of the previous direction, thus forming initial cell packages (Figures 5c,f,g). Cultured cells were almost as long as wide (average 12.5 × 13.0 μm). The cultured cells exhibited a remarkable feature: a collar of globules, possibly representing tiny lipid droplets or large Golgi vesicles, around the nucleus (Figures 5c,f,g). The chloroplast was divided into several lobes, which were appressed to the wall (Figures 5f,g). No pyrenoids were seen. The combination of the latter two features is described as characteristic of the green algal genus Apatococcus (Gärtner and Ingolić, 1989). Unfortunately, after about 2 weeks of growth, the filaments ceased growing and disintegrated into short fragments of a few two-celled units (Figures 5dg). The cells bleached, and empty walls became abundant (Figures 5d,f).

Figure 4
Microscopic images show various green algae filaments labeled from a to h. Each panel features close-up views of intertwined, segmented filaments with distinct green pigmentation. The scales range from 20 to 50 micrometers, highlighting the varying levels of magnification and details within the cellular structures.

Figure 4. Microscopy of environmental samples of the filamentous Apatococcus ammoniophilus. (a) Low magnification overview of filaments, cells with thickened walls, and numerous empty cells. (b) Filaments consisting of two-celled units (autosporangia), thickened cell walls with remnants of ruptured old walls. Note the many pale non-vital cells. (c) Filaments of vital cells and uniseriate arrangement of two-celled units. Note the cap-like structure of the thickened wall (left filament). (d) Short fragments of filaments, cells with lobed chloroplast. (e-h) Closer view of cell structures and thickened cell walls. Note that storage droplets are close to the cell walls. (e) Units of two cells with remnants of old walls on the walls. (f) Two two-celled units in detail. Note cap-like wall structures between two units and lobed chloroplast. (g) Detail of a filament with two-celled units with thickenings of remnants of old cell walls. Filamentous fungus attached to the algal filament. (h) Vital cells (with lobed chloroplast) and damaged cells, filamentous fungus attached to the algal filament. (a,b): sample US13027; (c,d,f): sample US13029; (e,g,h): sample US13028.

Figure 5
Microscopic images of filamentous green algae, segmented into labeled panels (a-g). The algae display various chain formations, with each panel showing different angles and close-ups, highlighting cell structure and formation. Scale bars indicate 20 micrometers.

Figure 5. Microscopy of cultured filaments of Apatococcus ammoniophilus (sample US13013). (a) grown filaments without considerable wall thickenings. Note the filaments constitute a characteristic arrangement of pairs of two cells (autospores). Note cell package formation (lower right). Cells with a chloroplast of several flat lobes. (b); Fragmented filament of a few two-celled units. (c,f,g); Views at cellular details. Note collar of globules around nucleus with nucleolus (arrow heads). (C) Dividing cells in a filament and formation of cell packages. Flat chloroplast with several lobes. (d,e) Disintegrating filaments with cells that undergo division after about 3 weeks of maintenance on agar plates. Arrangement of pairs of cells (autospores). (f,g) Chloroplast with several flat lobes, collar arrangement of globules around the nucleus (arrowheads). (f) Note cell package formation (upper right). (g) Filament with cell package formation at its lower end.

Identification and composition of the algal community

The paired-end ITS2 metabarcoding analysis of ten environmental samples revealed that one ASV was predominant in a community of various green algae, accounting for 61.6 to 97.2% of all read counts per sample. The total number of algal reads per sample varied from 72,804 (sample US13012) to 590,258 (sample US12969). A total of 60 green algal ASVs, clustered into 34 ASV clusters, was found (Figure 3). Interestingly, 37 ASVs exhibited high similarities (NB ≥ 1.76) with 19 closest reference sequences assigned “uncultured alga” (accessions in the range ON119209 - ON119407). These were from a previous metabarcoding analysis of algal communities on tree bark (Dreyling et al., 2022). The consensus-based protocol that we applied could assign taxonomic labels to these ASVs. They represented 11 green algal genera from the class Trebouxiophyceae, two from the Klebsormidiophyceae, and two more remained unidentified (Figure 3). The predominant ASV_00001, along with seven other ASVs, represented Apatococcus. This is a significant finding because it suggests that the dominant mass-forming filamentous alga is a species of that genus. An agar sample grown from US13013 was also included, where the same ASV_00001 was predominant (67.2% of read counts). This proves that the same alga that dominated the filamentous algal masses in the environmental sample US13013 also developed on the agar plate. In addition, most other green algae ASVs appeared to be amplified through culturing on the agar plate (Figure 3), apart from the other Apatococcus ASVs, which may indicate that those algae are hard to culture. Though the two Klebsormidiophyceae ASVs, representing the filamentous genera Interfilum and Klebsormidium, were recovered by the amplicon-based metabarcoding approach from only a single environmental sample, US12921, they also developed on the agar sample, but at a lower relative abundance than ASV_00001. Thus, Klebsormidiophyceae had a rather low abundance in the studied samples.

Phylogenetic position of the filamentous alga

Taking advantage of ASV_00001’s predominant presence in the ten samples selected for the metabarcoding approach, the nearly complete 18S rRNA genes from two of these samples were sequenced to enable the phylogenetic analysis of the filamentous alga. Long contigs of 2,976 and 3,020 base pairs were obtained from environmental samples US12924 and US12925 using the Sanger sequencing method. Both samples were also sequenced through the metabarcoding approach (Figure 3). The long sequences included the ITS1-5.8-ITS2 rRNA regions and the 5′-end of the 26S rRNA gene as well. The highly variable ITS2 of the two long contigs and that of ASV_00001 were identical; it was 211 base pairs long. It confirms that the 18S rRNA gene of the filamentous alga, which was dominant in the environmental samples, was sequenced. Direct comparison showed 100% sequence identity with the 18S rRNA gene and the ITS2 region of the cloned sequence DZK5, which we previously obtained from the surface of a roof tile in an urban region. Additionally, the ITS2 region of all three sequences matched the ITS2 region found in a sequence obtained from beech tree bark, accession number ON119340 (Dreyling et al., 2022). These findings indicate that the geographic distribution of the same genotype of the filamentous alga is broader than the samples we collected (Table 1) may suggest, and that it may colonize also other substrates besides tree bark. Phylogenetic analyses of the 18S rRNA gene placed the filamentous alga of the two environmental samples within a clade of the Trebouxiophyceae representing the genus Apatococcus (Figure 6). Thus, the analyses confirmed the sequence comparisons using the ASV_00001 sequence and proved the filamentous alga to be a species of Apatococcus.

Figure 6
Phylogenetic tree illustrating the evolutionary relationships among various algal species and clades, including Apatococcus, Viridiella, and others. Branches are labeled with species names and clone identifiers, with some highlighted sections such as “A. ammoniophilus” and letters indicating different clades or groups (A, B, C, D, E2). Bootstraps are indicated with purple circles of varying sizes, demonstrating confidence levels.

Figure 6. Maximum likelihood phylogenetic analysis (IQtree) of 18S rRNA gene sequences of Apatococcus ammoniophilus compared to 26 reference sequences representing Apatococcus and 117 other reference sequences representing major clades of the Trebouxiophyceae. The outgroup (two sequences of Nephroselmidales) has been pruned away from the graphic. Filled circles mark internal branches supported by >80% of 1,000 bootstrap replicates. The left figure shows the whole phylogeny with the Apatococcus clade highlighted. The right figure shows the internal structure of the Apatococcus clade with A. ammoniophilus highlighted. Marked lineages within Apatococcus are also resolved in the ITS2 rDNA phylogeny of Figure 7. Names in bold are those references containing 18S and ITS2 rDNA sequences.

The ITS2 of ASV_00001 and clone DZK5 were subjected to secondary structure-based phylogenetic analyses, including various reference sequences and the ITS2 sequences of the 34 ASV clusters as shown in Figure 3. The analyses confirmed the phylogenetic position of the filamentous alga in Apatococcus. Within the genus, the two new sequences and the sequence ON119340 formed an independent lineage (Figure 7) which suggests that the three sequences may represent a species of Apatococcus that has not been sequenced before. The Apatococcus clade was divided into subclades that roughly corresponded to the Apatococcus subclades discussed in Zahradníková et al. (2017). The subclade E of Zahradníková et al. (2017) was split into two lineages, E1 and E2, and the subclades A, B, and C (Zahradníková et al., 2017) were not separated from each other (Figure 7). The phylogenetic analyses also confirmed the identification of the other ASV clusters, i.e., their assignment to clades of Trebouxiophyceae genera was the same as obtained from the consensus-based protocol in the metabarcoding approach (Figures 3, 7). The two unidentified ASVs, ASV_00309 and ASV_00167, were together with their closest references, which have been assigned “Unidentified alga” (Dreyling et al., 2022) in a well-supported lineage, distinct from the other genera of Trebouxiophyceae and, therefore, remained unidentified also in the phylogenetic analyses. The lack of sufficiently close named reference sequence for those ASVs may indicate that the lineage may represent an yet unrecognized genus of Trebouxiophyceae. The phylogenetic position of the filamentous alga studied here within a clade that includes Apatococcus culture strains that exhibit a sarcinoid growth form, i.e., the culture strains SAG 2037, SAG 2145, and SAG 2359, was unexpected. Therefore, we tested the statistical significance of alternative tree topologies in which the filamentous alga, ASV_00001, occupied six other positions compared to the “best” tree from the maximum likelihood analysis presented in Figure 7 (Supplementary Figure S1). The two alternative tree topologies, where the ASV_00001 lineage was outside Apatococcus but adjacent to it (Supplementary figure S1), had high Delta log-likelihood (D-ML) values (p > 0.05), indicating they were significantly worse than the best tree and, consequently, unlikely to represent the true topology. The tests also rejected the other four alternative topologies with the ASV_00001 lineage positioned deeper within the Trebouxiophyceae, implying that any position of the alga other than inside Apatococcus was highly unlikely. Only the different arrangements within Apatococcus displayed low D-ML values (p < 0.05), suggesting they were plausible alternatives not significantly different from the best tree.

Figure 7
Phylogenetic tree of green algae shows relationships among various genera, including Trebouxia, Apatococcus, Coccomyxa, and Symbiochloris. Bootstraps are indicated with circles of varying sizes, demonstrating confidence levels.

Figure 7. Maximum likelihood phylogenetic analysis based on ITS2 secondary structure models of the 34 multi-ASV clusters from the paired-end approach green algal ASVs (Trebouxiophyceae) recovered by the amplicon-based metabarcoding analysis, represented in Figure 3 and closest reference sequence. They are distributed on ten well-supported clades representing different genera of Trebouxiophyceae and one unidentified clade. Highlighted is ASV_00001 and its closest related sequences representing Apatococcus ammoniophilus. Filled circles mark internal branches supported by >80% of 1,000 bootstrap replicates.

Metagenomic identification

Independent of the amplicon-based metabarcoding approach, a metagenomic identification of the filamentous alga was achieved. Illumina sequencing of the algal material from environmental sample US13019 resulted in 36,938,013 paired-end sequences of a size of 150 base pairs, totaling over 11 billion base pairs. Although FastQC results confirmed good sequencing results, we used trimmomatic to improve the read data, which removed 5,286,889 paired-end reads after strict filtering, resulting in 31.651,124 paired-end reads for a metagenomic assembly. We first assembled all untrimmed reads using SPAdes but later attempted another assembly with metaSPAdes due to the metagenomic nature of our data. A comparison between SPAdes and metaSPAdes revealed that SPAdes generated more complete assemblies (Supplementary Table S1); therefore, we continued all further investigations with the contigs of the initial SPAdes assembly.

The program barrnap identified 196 contigs containing 26S, 18S, 5.8S, and/or 5S genes and the program SSU-ALIGN identified 158 contigs among the 482,599 contigs in the SPAdes assembly. When combining and deduplicating these results, we receive 324 contigs, indicating that both programs identify mostly different contigs with ribosomal genes. These contigs showed considerable variation in length and coverage, with an average length of 5,114.1 bp (SD = 16,773.47 bp), average k-mer coverage of 29.31x (SD = 94.06 x), and average per-site coverage of 446.54 x (SD = 2,917.59 x) (Supplementary Table S1). On the barrnap contigs, a total of 694 rDNA genes were identified and subjected to reciprocal BLASTn searches against NCBI’s non-redundant nucleotide (nt) database. This analysis identified 34 contigs of Trebouxiophyceae origin in the metagenomic assembly, with an average length of 14,967.35 bp (SD = 22,342.96 bp), average k-mer coverage of 98.25 x (SD = 205.28 x), and average per-site coverage of 283.39 x (SD = 459.89 x) (Supplementary Table S1). Barrnap did not identify any contigs with ribosomal genes of Klebsormidiophyceae origin. SSU-ALIGN identified 18 contigs with 18S ribosomal genes of Trebouxiophyceae origin, with an average length of 4,256.17 bp (SD = 5,781.62 bp), average k-mer coverage of 147.40 x (SD = 294.06 x), and average per-site coverage of 396.89 x (SD = 599.38 x) (Supplementary Table S1). SSU-ALIGN also identified one contig with parts of the 18S ribosomal gene of Klebsormidiophyceae origin of 306 bp length, a k-mer coverage of 2.11 x, and an average per-site coverage of 5.4 x (Supplementary Table S1). All identified contigs were ranked by k-mer and per-site coverage to assess genomic abundance and identify the algal species.

Seven contigs categorized as Trebouxiophyceae displayed high coverage, exceeding 100 x. Both SSU-ALIGN and barrnap identified 18S ribosomal genes on four of these contigs, which showed high alignment identity (97.9–100%) in a megaBLAST search with Apatococcus (NODE_230161, NODE_1182952, NODE_205589) or Symbiochloris/Dictyochloropsis (NODE_109413) (Supplementary Table S1). Additionally, barrnap detected the 5.8S ribosomal gene on NODE_205589 and NODE_109413. The 5.8S gene on NODE_205589 aligned with 98.1% identity to Apatococcus, while the gene on NODE_109413 matched 100% identity to Symbiochloris/Dictyochloropsis. Furthermore, barrnap identified the 28S ribosomal genes on three high-coverage contigs, NODE_233586, NODE_95039, and NODE_16933, which blasted with lower alignment identity (80.9 to 94.5%) to Coccomyxa/Chlorella, Coccomyxa, and Lobospheara/Chlorella, respectively (Supplementary Table S1). The highest coverage for all Apatococcos-identified contigs was recorded at 902.76 x, 751.16 x, and 677.58 x k-mer coverage, as well as for one of the Coccomyxa contigs containing 26S, which had 716.92 x k-mer coverage. All other identified contigs had k-mer coverages of 354.6 x or less. Using SSU-ALIGN, we identified one contig with parts of a Klebsormidium 18S ribosomal gene on contig NODE_1112260, which had a blast hit with a low identity of 90.8% and a low k-mer coverage of 2.11 x, while barrnap did not find any Klebsormidiophyceae ribosomal genes in the metagenome assembly.

The Apatococcus-identified contigs NODE_1182952, NODE_230161, and NODE_205589 were compared to the partial ribosomal repeat of samples US12924 (acc. no. PQ763402) and US12925 (acc. Number PQ763401) from the amplicon-based approach, which contained parts of the 18S, full ITS1, 5.8S, ITS2, and parts of the 26S rRNAgenes. The NGS sequences were consistent with those from the amplicon-based approach, including the highly variable ITS1 and ITS2 regions (Supplementary Figure S2). Mapping all identical metagenomic NGS reads to the sequence PQ763402 demonstrated continuous and uniform coverage across the reference sequence, with an average depth of 1,043.3 x (SD = 235.0 x), highlighting the robustness and depth of the Apatococcus ribosomal repeat region in the metagenome (Supplementary Figure S2).

Polyols and mycosporine-like amino acids of the filamentous alga

Biochemical analysis of low molecular weight carbohydrates of the biofilm samples dominated by the filamentous alga resulted in the identification of the polyols arabitol and erythritol, as well as the disaccharide trehalose. Erythritol is a C4 polyol and occurred as the quantitatively dominant carbohydrate with 26.6 μg mg−1 dry weight, while arabitol as C5 polyol was measured at 2.8 μg mg−1 dry weight (Figure 8). Trehalose concentration was determined at 10.4 μg mg−1 dry weight (Figure 9).

Figure 8
Analytical data for Apatococcus ammoniophilus extract is shown in a series of graphs and chemical structures. Panel (a) features chromatograms with various peaks labeled with specific compounds. Panel (b) illustrates a mass spectrum with labeled peaks and fragmentation patterns for precursor ions. Panel (c) depicts chemical structures for mycosporine-serine, showing potential fragmentation pathways with annotations for loss of CO2, H2O, and a methyl group.

Figure 8. (a) HPLC-chromatograms of the Apatococcus ammoniophilus water extract (c = 5 mg mL−1), a standard mix with five mycosporine-like amino acids (MAAs) (shinorine (1 | RT = 1.24 min), palythine (2 | RT = 1.54 min), asterina-330 (3 | RT = 1.85 min), porphyra-334 (4 | RT = 2.31 min), aplysiapalythine A (5 | RT = 3.13 min)), and enriched fractions containing various MAAs aplysiapalythine B (6 | RT = 4.99 min), mycosporine-glycine–alanine (7 | RT = 3.50 min), mycosporine-methylamine-threonine (8 | RT = 5.31 min), palythene (9 | RT = 11.60 min), usujirene (10 | RT = 11.58 min) recorded at a wavelength of 330 nm. An on-line UV–Vis spectrum (λ = 200–600 nm) of the main peak in the A. ammoniophilus extract (RT = 4.03 min) is given in the chromatogram in the upper right section. (b) Mass spectra (MS2) of the main peak in the A. ammoniophilus extract (RT = 4.03 min) in positive (black) and negative (red) ESI mode. Characteristic fragments are highlighted with an arrow and the corresponding mass losses are indicated. (c) Hypothetical fragmentation pattern of mycosporine-serine, explaining several ions seen in the mass spectrum (B, negative ESI mode).

Figure 9
Bar graph showing concentrations of low molecular weight carbohydrates in micrograms per milligram dry weight for trehalose, erythritol, and arabitol. Erythritol has the highest concentration around 30, trehalose is about 10, and arabitol is the lowest under 5. Error bars indicate variability.

Figure 9. Low molecular weight carbohydrate pattern in Apatococcus ammoniophilus scratched off various bark samples from Pinus species in Jutland, Denmark (US12923, 12,924, 12,925, 12,927). Concentrations are given as μg mg−1 dry weight and represent mean values ± standard deviation (n = 3).

Chromatographic analysis of the biofilm water extract revealed the presence of a highly polar UV-absorbing compound that appeared as a prominent peak at a retention time of 4.03 min. and exhibited absorption characteristics consistent with those of mycosporine-like amino acids (MAAs) (Figure 8a). The online UV–Vis spectrum showed two maxima, at 224 and 324 nm, which fits well into their frequently stated absorption range of 310–360 nm (Karsten, 2008). Since the maximum at 324 nm is closer to that of MAAs with only one amino-acid substituent such as mycosporine-glycine (λmax = 310 nm) and palythine (λmax = 320 nm) and clearly below that of the aminocyclohexenimines porphyra 334 (λmax = 334 nm), usujirene, and palythene (both λmax = 359 nm), an aminocyclohexenone scaffold of the unknown compound seemed likely. Comparative liquid chromatographic studies showed that the retention time of the compound in the studied biofilm extract (RT = 4.03 min) differed significantly from the retention times of known MAAs listed in Figure 8a. The unknown compound eluted between mycosporine-glycine–alanine (7) and aplysiapalythine B (6) with retention times of 3.50 and 4.99 min. In the high-resolution MS experiments, the main peak of the biofilm extract was characterized by signals at m/z values of 274.0925 [M-H]- (ESI-) and 276.1075 [M + H] + (ESI+) (Figure 8b), suggesting a molecular formula of C11H17NO7 and consequently a neutral mass of 275 Da. The absorption properties, as well as the observed mass of the still unknown compound, suggest that it could be mycosporine-serine, which also complies with the hypothesis of an aminocyclohexenone core. Fragmentation reactions observed for the precursor ion in the positive mode (Figure 8B) include the elimination of water (Δm/z = 18, 276.11→258.10), followed by the loss of two methyl groups (Δm/z = 2 × 15, 258.10→228.09) and a C2H3O2∙moiety (Δm/z = 59, 228.09→169.07). This fragmentation pattern is consistent with that of mycosporine-serine reported by Volkmann and Gorbushina (2006). Moreover, the fragmentation pattern of the negatively charged precursor ion (m/z 274.09 [M-H]-) can also be explained by reactions starting from mycosporine-serine, as illustrated in Figure 8c.

Description of a new species

Both the amplicon-based and the metagenomic approaches identified Apatococcus as the most abundant sequence in the environmental samples, with the identified contigs from Apatococcus showing the highest coverage across the metagenomic assembly. This result demonstrates consistency across different methodologies and validates the dominance of Apatococcus in the sampled material. The exceptionally high coverage of these contigs further underscores the abundance of Apatococcus within the environmental sample. Additionally, the congruence of ribosomal gene regions, including the ITS1 and ITS2 sequences, supports the robust classification of Apatococcus as the predominant species in the microbial community analyzed. Therefore, we describe a new species of Apatococcus for the filamentous alga that prevails in the observed field material here:

Apatococcus ammoniophilus Søchting, Friedl & Moestrup sp. nov.

Diagnosis

Thallus forming up to 2 mm thick, light-green mats consisting of filamentous, unbranched, ca. 10 μm thick, uniseriate aplanospores with flat, lobed chloroplast and lack of pyrenoid. Contains erythritol. Differs from all other known species of Apatococcus by formation of filamentous thallus. On living or dead bark of mostly coniferous trees (Figures 25).

Type

Denmark, West Jutland, 10 km SW of Varde, 55.5459°N, 8.5780°E, on Scandosorbus intermedia bark, collected on 24 April 2024, Søchting 13,005 (C); dried material deposited at Statens Naturvidenskabelige Museum, Copenhagen as no C-A-100671 (C-holotype; BM, GOET, LD, L-isotypes).

Description

Thallus forming up to 2 mm thick, light green mats consisting of filamentous, unbranched, uniseriate, 9–12 μm thick aplanospores formed by the adhesion of two-celled units. Cell lumina 6–7 μm thick, with 0.3 to 0.5 μm thick walls that are stratified and irregularly grooved. Cell lumina typically 3–9 μm long with a single flat chloroplast with irregularly lobed edges appressed to the cell wall; without pyrenoid. Cells with a large, central nucleus with a clearly visible nucleolus and a ring of small particles surrounding it. The cytoplasm close to the cell wall is granular and often filled with oil drops. The reserve product may be oil, not starch. Zoospores not observed, but filaments have intercalary growth and vegetative reproduction is achieved through disruption of filaments and formation of non-motile autospores. Cells containing erythritol.

In culture the filaments disintegrate shortly when the cell walls of previous generations dissociate, leaving behind units where the cells are grouped in pairs (two-celled units). A wall thickening may be present at one apical end of the two-celled units. Packages of cells (sarcinoid stages) can be present alongside the filaments and form through a cell division plane that is perpendicular to the one that results in the filaments; the filament cells are slightly elongated. Oil droplets vanish in culture. The cultured cells also feature a large central nucleus with a clearly visible nucleolus, surrounded by a ring (collar) of small particles.

Etymology

Named after its observed preference for strong atmospheric ammonia deposition.

Ecology

Terrestrial and subaerial, forming algal mats on the living and dead bark of coniferous and deciduous trees and shrubs, but may also grow on artificial hard substrates.

Distribution

Recorded so far from Denmark, Sweden, the Netherlands, Germany, and Great Britain. Very frequent in England, and also occurring in Scotland and Wales (A. Pentecost, in lit.). Most likely found in most of Northwestern Europe. The list of studied specimens is presented in Table 1.

DNA sequences

The ITS2 rDNA sequences are deposited in the NCBI GenBank database under accession number PQ784906 (ASV_00001), the 18S 3′-partial, ITS1, 5.8S, ITS2, and 26S 5′-partial sequences under the accession numbers PQ763401 (from voucher Søchting 12,925) and PQ763402 (from voucher Søchting 12,924). Metagenomic raw reads are also deposited at NCBI GenBank with the accession number SRR32987805 under BioProject PRJNA1246228 and Biosample SAMN47770417.

Discussion

Taxonomy

The filamentous alga studied here has earlier been assigned to Hormidium crenulatum Kütz., when a comprehensive morphological description with illustrations of cellular details was provided by Petersen (1915). Our observations of the studied filamentous alga, mostly from tree bark samples of Danish origin, fully match his descriptions. Therefore, our morphological observations on the fresh environmental (Figures 2-4) and the cultured material (Figure 5) leave no doubt that the alga recovered in our study here is the same as that of Petersen (1915). The first taxonomic description of H. crenulatum (Kützing, 1845) is incomplete because it refers only to the rough feature of the cell wall and the cell dimensions. In addition, Kützing’s description does not include drawings from which further cellular details could be inferred. Hormidium crenulatum Kütz. was lectotypified by Lokhorst in Lokhorst and Star (1985) because it was impossible to find any holotype material relating to Hormidium crenulatum Kütz. The lectotype in Leiden (L): “L 939.67–834, annotated as “1/280” Hormidium crenulatum, Patavii, on wet and warm wall of a bath house, Padua, leg. G. Meneghini” has not been available for microscopic examination. Based on many cultured isolates, mostly from soil samples, Lokhorst (1996) published a new combination, Klebsormidium crenulatum (Kütz.) Lokhorst. Superficially, K. crenulatum resembles A. ammoniophilus as it can form filaments in culture with cell walls that appear “thickened and pronouncedly rough, (slightly) dented or crenulated” (Lokhorst, 1996). However, unlike the alga studied here and described by Petersen (1915), K. crenulatum regularly forms characteristic H-wall pieces. The K. crenulatum chloroplast is “plate- to girdle-shaped” and closely appressed to the cell wall; it possesses a pyrenoid with a prominent shell of starch grains (Lokhorst, 1996). In K. crenulatum the nucleus is within a central bridge of the cytoplasm, connecting opposite sides of the cell wall, a feature characteristic of many Klebsormidiophyceae (e.g., Interfilum, Mikhailyuk et al., 2008). Later, a culture strain, SAG 37.86, isolated from Alpine soil, was designated epitype of K. crenulatum to unambiguously link the name to a sequenced specimen (Mikhailyuk et al., 2015). The phylogenetic position of K. crenulatum as a distinct lineage within the genus Klebsormidium in the Klebsormidiophyceae (Streptophyta) was established (Mikhailyuk et al., 2015). Considering the morphological resemblance between A. ammoniophilus and K. crenulatum, the images of the filamentous K. crenulatum displayed in AlgaeBase (Guiry and Guiry, 2025) depict the same alga we studied here, but not K. crenulatum Kütz., underscoring the potential for confusion between the two. It follows that the alga we studied here, corresponding to the one reported by Petersen (1915), lacks a valid taxonomic name. Our study also uncovers that the alga is phylogenetically distant from K. crenulata because it is a member of the Trebouxiophyceae (Chlorophyta) and within the monophyletic clade representing the genus Apatococcus. Therefore, the alga needed to be described as a new species of Apatococcus. Gustavs et al. (2011) used the occurrence of specific polyols as an additional chemosystematic marker for different Trebouxiophyceae and showed that erythritol was exclusively detected in various Apatococcus strains. To the best of our knowledge, we are not aware of any other Trebouxiophyceae taxa containing this specific polyol, and hence the taxonomic position of A. ammoniophilus is well supported by this chemical marker. The only other algae known to also synthesize and accumulate erythritol are members of the aeroterrestrial genus Trentepohlia, which is in the Ulvophyceae (Chlorophyta) (Holzinger et al., 2023).

It may seem unexpected that a species of Apatococcus exhibits filamentous growth. Cell packages, i.e., a sarcinoid growth form, are typically associated with the species A. lobatus. The “parenchyma-like cell complex” or „multicellular mass “(Kinutani et al., 2015) or” aggregate building “(Gustavs et al., 2010; Gustavs et al., 2011) has often been described for A. lobatus. This species is a common member of algal communities on tree bark (e.g., Geitler, 1942; Vischer, 1960; Gärtner and Ingolić, 1989; Dreyling et al., 2022; Handa and Nakano, 1988). It is also very abundant in aeroterrestrial biofilms, particularly in urban areas, where it can cause extensive discolorations on walls and roofs (Rindi, 2007; Hallmann et al., 2013, 2016; Karsten et al., 2022). A. ammoniophilus filaments form through the adhesion of two-celled units (autosporangia). This feature corresponds to A. lobatus, where two autospores (immotile daughter cells) are formed within an autosporangium (Kinutani et al., 2015). The formation of only two cells by the division of the initially unicellular vegetative cells of A. lobatus has often been observed in earlier studies (Brand and Stockmayer, 1925; Geitler, 1942; Vischer, 1960; Gärtner and Ingolić, 1989). Characteristic of A. lobatus are cell packages that consist of three cells. The study of Kinutani et al. (2015) demonstrates that the A. lobatus three-celled packages arise from a delay in the formation of autospores, i.e., while two autospores are formed in one autosporangium, no division occurs in the other. In A. ammoniophilus, too, sarcinoid growth with packages of autosporangia can be observed in culture (Figure 5f), where the two daughter cells are formed at different times. However, in A. ammoniophilus, the cell division plane persists, resulting in chains of two-celled units that adhere to each other and result in intercalary elongation of the filaments. One could speculate that certain environmental conditions may suppress the change of the division plane, favoring filament formation.

One more remarkable feature connects A. ammoniophilus with A. lobatus: in the cultured cells, there is a centrally located ring of small particles (“Zentralkranz”) surrounding the nucleus (Brand and Stockmayer, 1925; Gärtner and Ingolić, 1989), which is clearly visible with Nomarski optics (Figure 5). According to the original diagnosis of the genus (Brand and Stockmayer, 1925, p. 351, figs 48, 50, 51), this feature distinguishes Apatococcus in culture from other terrestrial green algae found on tree bark.

Adaptive traits of A. ammoniophilus for the aerophytic lifestyle

Aerophytic microalgae exhibit a variety of morphological, physiological, and biochemical traits for acclimation and adaptation to seasonally fluctuating environmental conditions (Holzinger and Karsten, 2013). Water availability and desiccation, along with high solar radiation, are key ecological factors determining the abundance and diversity of green algae in terrestrial biofilms. A. ammoniophilus contains three low molecular weight carbohydrates, i.e., arabitol, erythritol, and trehalose, with erythritol showing the highest concentrations (Figure 9). Both polyols are important protective compounds against desiccation-induced stress (Gustavs et al., 2010). However, erythritol and arabitol exhibit different chemical properties and energy requirements for biosynthesis, as erythritol is the smallest polyol with four carbon atoms, while arabitol consists of five carbon atoms. Trehalose is a disaccharide known to be formed particularly by desiccation-tolerant organisms, as it stabilizes biomolecules in a dehydrated state. The underlying mechanisms are still hypothesized rather than experimentally verified. These mechanisms include, for example, vitrification, as trehalose forms a glassy matrix that acts as a barrier, presumably physically shielding proteins or membranes from water loss (Kageyama and Waditee-Sirisattha, 2023).

Besides this interesting low molecular weight carbohydrate pattern, A. ammoniophilus also contained a probably new UV-absorbing MAA for Trebouxiophycean algae, as most members of this group contain the chemically elucidated MAA prasiolin (Hotter et al., 2018). Volkmann and Gorbushina (2006) published the only traceable report on mycosporine-serine, although the authors also refer to previous studies. As a consequence, the true nature of this interesting compound in A. ammoniophilus seems questionable, especially as its identification in the fruiting bodies of the fungus Stereum hirsutum was purely based on MS and UV–Vis data. Although the results presented herein (absorption maximum, retention time, mass, MS fragmentation pattern) strongly indicate the tentative identity of the unknown compound as mycosporine-serine, this still needs to be verified in additional experiments aimed at isolation and structural elucidation. This requires large amounts of clonal biomaterial, which makes a follow-up study necessary. In any case, final proof can only be provided by submitting unambiguous NMR and high-resolution MS data of the pure compound isolated from a clearly identified organism. The presence and properties of three low molecular weight carbohydrates, along with a rather unique MAA in A. ammoniophilus, well explain the biochemical basis for an aeroterrestrial lifestyle, i.e., the capability to synthesize anti-stress metabolites for compensation of desiccation and UV.

The filamentous morphology with intercalary growth may be advantageous as it allows a continuous and unlimited vertical increase of the filamentous mat. The formation of a thick mat of filaments will allow the species to expand over the substrate and overgrow more crustose competitors, resulting in almost monospecific ecological communities, but it may also result in shading of the lower part of the mat, eventually resulting in flaking off. Coccoid aerophytic algae have a low ability to take up free water (Barkmann, 1958), but are often successful on protected surfaces because they are able to absorb moisture from the air humidity. On the contrary, A. ammoniophilus grows on exposed surfaces on the twigs. It is able to instantly absorb moisture from dew or precipitation in the spongy mat, and additionally have a very large surface exposed to ammonia from the air. Interestingly, their light exposed habitat may also demand special protection by, e.g., MAA’s compared with their shade-preferring relatives. Field studies in combination with laboratory studies on clonal material are needed to investigate growth rates and other physiological traits under different environmental conditions, with a focus on ammonia requirements.

A. ammoniophilus expansion in relation to increased ammonia deposition

Apatococcus ammoniophilus passed fairly unnoticed until it became abundant, most probably in response to decreasing acidification from acid rain in the seventies, and subsequent increasing ammonia deposition due to increased industrial farming. Estimates by Alveteg et al. (1998) and Ellermann et al. (2007) show a continuous average increase in nitrogen deposition from about 1 kg N ha−1 in 1800 to about 16 kg N ha−1 in 1990 (Figure 10). Over the subsequent 30 years nitrogen deposition has declined to about 12 kg N ha−1 in 2023, which is still 12-fold above pristine conditions (Ellermann et al., 2024). The species was unknown to the late Danish phycologist Dr. Tyge Christensen (Christensen, in litt), an extremely experienced phycologist and field biologist, indicating that it had expanded its occurrence only in recent years. Already Petersen (1915) noticed the preference of A. ammoniophilus for nitrogen enrichment. The locations of the studied samples of A. ammoniophilus were superimposed on the map of modelled total nitrogen deposition in Denmark in 2020 (Ellermann et al., 2021; EMEP, 2024) (Figure 1). However, the locations of the >100 A. ammoniophilus collections do not properly illustrate the substantial algal biomass abundance observed in the high-deposition regions in southern and western Jutland and the very sparse biomass in the low-deposition eastern part of the country, where A. ammoniophilus is only found if thoroughly searched for. The predominant occurrence in western Jutland is accentuated by sandy soil resulting in dominance of coniferous plantations, which appear to provide the most favorable substrate for A. ammoniophilus. Intensive farming areas in Denmark are those that have in recent years also received increased precipitation. However, there is no reason to expect this change to explain lichen disappearance; increased precipitation may even have increased wet deposition of ammonium.

Figure 10
Line graph showing annual nitrogen deposition in kilograms of nitrogen per hectare from 1800 to 2000. Values remain stable at around 2 kg N ha-1 until about 1950, after which they sharply increase, reaching approximately 16 kg N ha-1 by 2000.

Figure 10. Average nitrogen deposition in Denmark 1800–1990 (kg N/ha) based on modelled estimates from Alveteg et al. (1998) and Ellermann et al. (2007).

Despite the lack of evidence from correlation data, careful field observations over many years and the abundance pattern concur that the occurrence and abundance of A. ammoniophilus might be related to atmospheric nitrogen deposition. Accordingly, we assume that it may even serve as a biological indicator of elevated ammonia deposition. The poor growth of A. lobatus in culture could be overcome by the addition of organic nutrient sources indicating mixotrophic behavior (Gustavs et al., 2016), which can also be expected for A. ammoniophilus.

A. ammoniophilus occurrence and distribution

Apatococcus ammoniophilus is associated with other terrestrial green algae of the Trebouxiophyceae. The amplicon-based DNA metabarcoding revealed a distinct pattern with only slight variation across ten replicates, specifically nine tree bark samples and one coniferous needle sample (Figure 3). In all the samples, in addition to A. ammoniophilus, additional genotypes of Apatococcus were present. Those and most other genotypes recovered were the same as those revealed by a recent DNA metabarcoding-based study of tree bark microorganisms (Dreyling et al., 2024). It seems that there is a widespread terrestrial green algal community commonly found on tree barks in Central Europe, consisting of various genotypes and species from the Trebouxiophyceae, mainly from the genera Apatococcus, Chloroidium, Coccomyxa, Symbiochloris, and Trebouxia. A. ammoniophilus also belongs to that widespread community.

Apatococcus ammoniophilus has frequently been reported to one of us (US) by colleagues, who identified it as K. crenulatum based on field material observations from various locations in Great Britain. There, it is often viewed as a nuisance alga because of the significant biomass it forms on twigs of diverse trees lacking lichen cover. These observations support a broader distribution of the filamentous alga in Northern European regions with a more humid and cooler climate. The alga was also observed in England, where its appearance correlated with the disappearance of epiphytic lichens (A. Pentecost, in litt.). It is not studied whether the disappearance of lichens is caused by competition from A. ammoniophilus or is a direct effect of ammonia on the lichen thallus or its substrate, but when present A. ammoniophilus leaves very little space for lichens to colonize. From the United Kingdom, there are numerous observations on social media about the thick cover of a filamentous alga called Klebsormidium crenulatum on rocks and other solid surfaces in regions with heavy nitrogen deposition.3

Our DNA-based analyses, at the very fine level of sequence variants, revealed evidence of the same alga occurring in even more southern regions of Europe with a warmer and drier climate. The analyses recovered the same A. ammoniophilus genotype from locations presumably without elevated atmospheric nitrogen levels, i.e., a roof tile from a building in an urban environment and the bark of beech trees in forests. There, only inconspicuous green biofilms devoid of filamentous green algal growth were found.

Conclusion

A microscopic green alga is forming significant biomass, becoming so obvious that it is regarded as a nuisance. This demonstrates that changes in terrestrial ecosystems, such as increasing deposition of eutrophicating nitrogen compounds from the atmosphere or climate warming, may strongly influence the algal communities exposed to open air. As an essential prerequisite for further studies on this phenomenon, our study provided insights into the phylogenetic position, taxonomy, and biochemical adaptive traits of a filamentous green alga that has expanded its occurrence only in recent years. We show that the alga, which was known under the name Klebsormidium crenulatum for more than 100 years, is, however, phylogenetically and structurally distinct from the latter species. We clarified the green alga’s phylogenetic position within the Trebouxiophyceae (Chlorophyta) and its close relationship within the genus Apatococcus. These efforts resulted in the description of a new species, A. ammoniophilus. There is consistent support for its classification across different DNA sequence-based methods and through an independent chemotaxonomic marker, i.e., the presence of the low molecular weight carbohydrate, erythritol. Future studies should focus on establishing culture strains of A. ammoniophilus, which, like its close relative A. lobatus, is difficult to maintain in culture. Furthermore, it is needed to experimentally test its apparent requirement for high levels of available ammonia.

Data availability statement

All sequence data generated from environmental samples as listed in Figure 3 and Table 1 using the paired-end approach (amplicon-based metabarcoding) have been deposited in NCBI in the Sequence Read Archive (SRA) under the BioProject accession number PRJNA1277513, with the BioSample accessions SAMN49111103 - SAMN49111113. One representative of each of the 34 multi-ASV clusters from the paired-end approach can be accessed from DDBJ/EMBL/GenBank databases with the accession numbers PQ784906 - PQ784939, the long 18S-ITS1, 2-26S rDNA sequences obtained from vouchers US12924, US12925, and clone DZK5 from the DDBJ/EMBL/GenBank databases with the accession numbers PQ763401 - PQ763403. All metagenomic raw reads can be accessed from DDBJ/EMBL/GenBank database with the accession number SRR32987805 under BioProject PRJNA1246228 and Biosample SAMN47770417.

Author contributions

US: Data curation, Project administration, Writing – original draft, Writing – review & editing, Funding acquisition, Resources. TF: Funding acquisition, Investigation, Writing – original draft, Writing – review & editing, Formal analysis, Visualization. ØM: Formal analysis, Investigation, Validation, Writing – review & editing. FG: Formal analysis, Funding acquisition, Investigation, Software, Writing – original draft, Writing – review & editing. YS: Formal analysis, Writing – review & editing, Conceptualization. YÇ: Writing – review & editing, Investigation, Visualization. MG: Investigation, Writing – review & editing, Funding acquisition. KG: Formal analysis, Investigation, Methodology, Writing – review & editing. SH: Formal analysis, Investigation, Writing – review & editing. FH: Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft. DN: Investigation, Writing – review & editing. BO: Formal analysis, Investigation, Writing – review & editing. UK: Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. The study was supported by the German Science Foundation (DFG), grant number 512414255, within the program “Biodiversity Exploratories (SPP 1374),” extended to TF, and by the FWF-DFG DACH project “UVision” under grant numbers KA899/45–1 (extended to UK) and FWF I 6122-B (extended to MG) (Grant-DOI 10.55776/I6122, https://www.fwf.ac.at/forschungsradar/10.55776/I6122). FG and YS received funding from the Grainger Foundation. Dept. of Biology, University of Copenhagen funded the printing.

Acknowledgments

Anne Grete Lenz, Bjarke Huus Jensen, Thomas Læssøe, Marco Brodde and Allan Pentecost are thanked for providing collections. Niklas Plag performed polyol analyses and Bjørn Hermansen produced the distribution map; Christian Lange of the herbarium in Copenhagen (C) was instrumental with loans, Nicolien Sol of Naturalis (L) gave information on the lectotype of Hormidium crenulatum, and António José Calado gave linguistic advice. They are all warmly thanked. We thank Kerstin Holst and Elke Zufall-Roth for determining the clone DZK5 sequence. Thorsten Lumbsch is thanked for fruitful discussions.

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.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

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.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2025.1633308/full#supplementary-material

Footnotes

References

Adams, R., and Degiorgio, M. (2023). Likelihood-based tests of species tree hypotheses. Mol. Biol. Evol. 40:159. doi: 10.1093/molbev/msad159

PubMed Abstract | Crossref Full Text | Google Scholar

Altschul, S. F., Gish, W., Miller, W., Myers, E. W., and Lipman, D. J. (1990). Basic local alignment search tool. J. Mol. Biol. 215, 403–410. doi: 10.1016/S0022-2836(05)80360-2

PubMed Abstract | Crossref Full Text | Google Scholar

Altschul, S. F., Madden, T. L., Schäffer, A. A., Zhang, J., Zhang, Z., Miller, W., et al. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389–3402. doi: 10.1093/nar/25.17.3389

PubMed Abstract | Crossref Full Text | Google Scholar

Alveteg, M., Walse, C., and Warfsvinge, P. (1998). Reconstructing historic atmospheric deposition and nutrient uptake from present day values using Makedep. Water Air Soil Pollut. 104, 269–283.

Google Scholar

Bandaranayake, W. M. (1998). Mycosporines: are they nature’s sun-screens? Nat. Prod. Rep. 15, 159–172. doi: 10.1039/A815159Y

PubMed Abstract | Crossref Full Text | Google Scholar

Barkmann, J. J. (1958). On the ecology of cryptogamic epiphytes with special reference to the Netherlands. Assen: Van Gorcum & Co. N.V.

Google Scholar

Bengtsson-Palme, J., Ryberg, M., Hartmann, M., Branco, S., Wang, Z., Godhe, A., et al. (2013). Improved software detection and extraction of ITS1 and ITS2 from ribosomal ITS sequences of fungi and other eukaryotes for analysis of environmental sequencing data. Methods Ecol. Evol. 4, 914–919. doi: 10.1111/2041-210X.12073

Crossref Full Text | Google Scholar

Bertsch, A. (1966). CO2 exchange and water relations in the aerophytic green alga Apatococcus lobatus. Planta 70, 46–72. doi: 10.1007/BF00539910

PubMed Abstract | Crossref Full Text | Google Scholar

Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170

PubMed Abstract | Crossref Full Text | Google Scholar

Brand, F., and Stockmayer, S. (1925). Analyse der aerophilen Grünalgenanflüge, insbesondere der proto- pleurococcoiden Formen. Arch. Protistenkunde 52, 265–355.

Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

Criscuolo, A., and Gribaldo, S. (2010). BMGE (block mapping and gathering with entropy): a new software for selection of phylogenetic informative regions from multiple sequence alignments. BMC Evol. Biol. 10:210. doi: 10.1186/1471-2148-10-210

PubMed Abstract | Crossref Full Text | Google Scholar

Di Tommaso, P., Chatzou, M., Floden, E. W., Barja, P. P., Palumbo, E., and Notredame, C. (2017). Nextflow enables reproducible computational workflows. Nat. Biotechnol. 35, 316–319. doi: 10.1038/nbt.3820

PubMed Abstract | Crossref Full Text | Google Scholar

Dreyling, L., Penone, C., Schenk, N. V., Schmitt, I., and Dal Grande, F. (2024). Biotic interactions outweigh abiotic factors as drivers of bark microbial communities in central European forests. ISME Commun. 4:ycae012. doi: 10.1093/ismeco/ycae012

PubMed Abstract | Crossref Full Text | Google Scholar

Dreyling, L., Schmitt, I., and Dal Grande, F. (2022). Tree size drives diversity and community structure of microbial communities on the bark of beech (Fagus sylvatica). Front. For. Glob. Change 5:858382. doi: 10.3389/ffgc.2022.858382

Crossref Full Text | Google Scholar

Ellermann, T., Andersen, H., Bossi, R., Christensen, J., Frohn, L. M., and Geels, C. (2007) Atmosfærisk deposition 2006. Danmarks Miljøundersøgelser. 62 s. – Faglig rapport fra DMU, nr. 645. Available online at: http://www.dmu.dk/Pub/FR645.pdf (Accessed June 01, 2025).

Google Scholar

Ellermann, T., Bossi, R., Sørensen, M. O. B., Christensen, J., Lansø, A. S., Monies, C., et al. (2024). Atmosfærisk deposition 2023 : NOVANA. Aarhus Universitet, DCE – Nationalt Center for Miljø og Energi. 68 – Videnskabelig rapport fra DCE – Nationalt Center for Miljø og Energi nr. 626.

Google Scholar

Ellermann, T., Bossi, R., Sørensen, M. O. B., Christensen, J., Løfstrøm, P., and Lansø, A. S. (2021) Atmosfærisk deposition 2020. NOVANA. Aarhus Universitet, DCE – Nationalt Center for Miljø og Energi. 95s. – Videnskabelig rapport fra DCE – Nationalt Center for Miljø og Energi nr. 471. Available online at: http://dce2.au.dk/pub/SR471.pdf

Google Scholar

Ewels, P. A., Peltzer, A., Fillinger, S., Patel, H., Alneberg, J., Wilm, A., et al. (2020). The nf-core framework for community-curated bioinformatics pipelines. Nat. Biotechnol. 38, 276–278. doi: 10.1038/s41587-020-0439-x

PubMed Abstract | Crossref Full Text | Google Scholar

EMEP,. (2024) EMEP Status Report 1/2024: “Transboundary particulate matter, photo-oxidants, acidifying and eutrophying components”. Joint MSC-W & CCC & CEIP & CIAM Report. Meteorological Synthesizing Centre-West of EMEP (MSC-W of EMEP) The Norwegian Meteorological Institute..

Google Scholar

Friedl, T. (1996). Evolution of the polyphyletic genus Pleurastrum (Chlorophyta): inferences from nuclear-encoded ribosomal DNA sequences and motile cell ultrastructure. Phycologia 35, 456–469. doi: 10.2216/i0031-8884-35-5-456.1

Crossref Full Text | Google Scholar

Friedl, T., and Rokitta, C. (1997). Species relationships in the lichen alga Trebouxia (Chlorophyta, Trebouxiophyceae): molecular phylogenetic analyses of nuclear-encoded large subunit rRNA gene sequences. Symbiosis 23, 125–148.

Google Scholar

Gärtner, G., and Ingolić, E. (1989). Ein Beitrag zur Kenntnis von Apatococcus lobatus (Chlorophyta, Chaetophorales, Leptosiroideae). Plant Syst. Evol. 164, 133–143. doi: 10.1007/BF00940434

Crossref Full Text | Google Scholar

Geitler, L. (1942). Morphologie, Entwicklungsgeschichte und Systematik neuer bemerkenswerter atmophytischer Algen aus Wien. Flora 136, 1–29. doi: 10.1016/S0367-1615(17)31209-0

Crossref Full Text | Google Scholar

Gruber, A. R., Bernhart, S. H., and Lorenz, R. (2015). “The ViennaRNA web services” in RNA bioinformatics. ed. E. Picardi (New York, NY: Springer New York), 307–326.

Google Scholar

Guiry, M. D., and Guiry, G. M. (2025). Algaebase. World-wide electronic publication, National University of Ireland, Galway. Available online at: http://www.algaebase.org (Accessed April 29, 2025)

Google Scholar

Gurevich, A., Saveliev, V., Vyahhi, N., and Tesler, G. (2013). QUAST: quality assessment tool for genome assemblies. Bioinformatics 29, 1072–1075. doi: 10.1093/bioinformatics/btt086

PubMed Abstract | Crossref Full Text | Google Scholar

Gustavs, L., Eggert, A., Michalik, D., and Karsten, U. (2010). Physiological and biochemical responses of green microalgae from different habitats to osmotic and matric stress. Protoplasma 243, 3–14. doi: 10.1007/s00709-009-0060-9

PubMed Abstract | Crossref Full Text | Google Scholar

Gustavs, L., Görs, M., and Karsten, U. (2011). Polyol patterns in biofilm-forming aeroterrestrial green algae (Trebouxiophyceae, Chlorophyta). J. Phycol. 47, 533–537. doi: 10.1111/j.1529-8817.2011.00979.x

PubMed Abstract | Crossref Full Text | Google Scholar

Gustavs, L., Schumann, R., Karsten, U., and Lorenz, M. (2016). Mixotrophy in the terrestrial green alga Apatococcus lobatus (Trebouxiophyceae, Chlorophyta). J. Phycol. 52, 311–314. doi: 10.1111/jpy.12381

PubMed Abstract | Crossref Full Text | Google Scholar

Hallmann, C., Hoppert, M., Mudimu, O., and Friedl, T. (2016). Biodiversity of green algae covering artificial hard substrate surfaces in a suburban environment: a case study using molecular approaches. J. Phycol. 52, 732–744. doi: 10.1111/jpy.12437

PubMed Abstract | Crossref Full Text | Google Scholar

Hallmann, C., Stannek, L., Fritzlar, D., Hause-Reitner, D., Friedl, T., and Hoppert, M. (2013). Molecular diversity of phototrophic biofilms on building stone. FEMS Microbiol. Ecol. 84, 355–372. doi: 10.1111/1574-6941.12065

Crossref Full Text | Google Scholar

Hamby, K. R., Sims, L., Issel, L., and Zimmer, E. (1988). Direct ribosomal RNA sequencing: optimization of extraction and sequencing methods for work with higher plants. Plant Mol. Biol. Report. 6, 175–192. doi: 10.1007/BF02669591

Crossref Full Text | Google Scholar

Handa, S., and Nakano, T. (1988). Some corticolous algae from Miyajima Island, western Japan. Nov. Hedwig. 46, 165–186.

Google Scholar

Hartmann, A., Glaser, K., Holzinger, A., Ganzera, M., and Karsten, U. (2020). Klebsormidin a and B, two new UV-sunscreen compounds in green microalgal Interfilum and Klebsormidium species (Streptophyta) from terrestrial habitats. Front. Microbiol. 11:499. doi: 10.3389/fmicb.2020.00499

PubMed Abstract | Crossref Full Text | Google Scholar

Hartmann, A., Holzinger, A., Ganzera, M., and Karsten, U. (2016). Prasiolin, a new UV-sunscreen compound in the terrestrial green macroalga Prasiola calophylla (Carmichael ex Greville) Kutzing (Trebouxiophyceae, Chlorophyta). Planta 243, 161–169. doi: 10.1007/s00425-015-2396-z

PubMed Abstract | Crossref Full Text | Google Scholar

Häubner, N., Schumann, R., and Karsten, U. (2006). Aeroterrestrial microalgae growing in biofilms on facades—response to temperature and water stress. Microb. Ecol. 51, 285–293. doi: 10.1007/s00248-006-9016-1

PubMed Abstract | Crossref Full Text | Google Scholar

Helms, G., Friedl, T., Rambold, G., and Mayrhofer, H. (2001). Identification of photobionts from the lichen family Physciaceae using algal-specific ITS rDNA sequencing. Lichenologist 33, 73–86. doi: 10.1006/lich.2000.0298

Crossref Full Text | Google Scholar

Hoang, D. T., Chernomor, O., Von Haeseler, A., Minh, B. Q., and Vinh, L. S. (2017). UFBoot2: improving the ultrafast bootstrap approximation. Mol. Biol. Evol. 35, 518–522. doi: 10.1093/molbev/msx281

PubMed Abstract | Crossref Full Text | Google Scholar

Holzinger, A., and Karsten, U. (2013). Desiccation stress and tolerance in green algae: consequences for ultrastructure, physiological, and molecular mechanisms. Front. Plant Sci. 4:327. doi: 10.3389/fpls.2013.00327

PubMed Abstract | Crossref Full Text | Google Scholar

Holzinger, A., Plag, N., Karsten, U., and Glaser, K. (2023). Terrestrial Trentepohlia sp. (Ulvophyceae) from alpine and coastal collection sites show strong desiccation tolerance and broad light and temperature adaptation. Protoplasma 260, 1539–1553. doi: 10.1007/s00709-023-01866-2

PubMed Abstract | Crossref Full Text | Google Scholar

Hotter, V., Glaser, K., Hartmann, A., Ganzera, M., and Karsten, U. (2018). Polyols and UV-sunscreens in the Prasiola-clade (Trebouxiophyceae, Chlorophyta) as metabolites for stress response and chemotaxonomy. J. Phycol. 54, 264–274. doi: 10.1111/jpy.12619

PubMed Abstract | Crossref Full Text | Google Scholar

Kageyama, H., and Waditee-Sirisattha, R. (2023). Halotolerance mechanisms in salt-tolerant cyanobacteria. Adv. Appl. Microbiol. 124, 55–117. doi: 10.1016/bs.aambs.2023.07.003

PubMed Abstract | Crossref Full Text | Google Scholar

Karsten, U. (2008). “Defense strategies of algae and cyanobacteria against solar UVR” in Algal chemical ecology. ed. C. D. Amsler (Berlin: Springer), 273–296.

Google Scholar

Karsten, U., Lembcke, S., and Schumann, R. (2007). The effects of ultraviolet radiation on photosynthetic performance, growth and sunscreen compounds in aeroterrestrial biofilm algae isolated from building facades. Planta 225, 991–1000. doi: 10.1007/s00425-006-0406-x

PubMed Abstract | Crossref Full Text | Google Scholar

Karsten, U., Peine, M., Gustavs, L., and Schumann, R. (2022). Apatococcus lobatus (Trebouxiophyceae) dominates green algal biofilms on roof tiles in northern Germany. Nova Hedwigia 114, 303–320. doi: 10.1127/nova_hedwigia/2022/0691

Crossref Full Text | Google Scholar

Karsten, U., Thomas, D. N., Weykam, G., Daniel, C., and Kirst, G. O. (1991). A simple and rapid method for extraction and separation of low molecular weight carbohydrates from macroalgae using high-performance liquid chromatography. Plant Physiol. Biochem. 29, 373–378.

Google Scholar

Katoh, K., and Standley, D. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772–780. doi: 10.1093/molbev/mst010

PubMed Abstract | Crossref Full Text | Google Scholar

Keller, A., Schleicher, T., Schultz, J., Müller, T., Dandekar, T., and Wolf, M. (2009). 5.8S-28S rRNA interaction and HMM-based ITS2 annotation. Gene 430, 50–57. doi: 10.1016/j.gene.2008.10.012

PubMed Abstract | Crossref Full Text | Google Scholar

Kinutani, K., Shimamura, M., Friedl, T., and Takeshita, S. (2015). Investigation of the life cycle of Apatococcus lobatus (Chodat) J.B.Petersen (Trebouxiophyceae, Chlorophyta). Hikobia 17, 33–40.

Google Scholar

Kircher, M., Sawyer, S., and Meyer, M. (2012). Double indexing overcomes inaccuracies in multiplex sequencing on the Illumina platform. Nucleic Acids Res. 40:e3. doi: 10.1093/nar/gkr771

PubMed Abstract | Crossref Full Text | Google Scholar

Kranner, I., Beckett, R., Hochman, A., and Nash, T. H. (2008). Desiccation-tolerance in lichens: a review. Bryologist 111, 576–593. doi: 10.1639/0007-2745-111.4.576

Crossref Full Text | Google Scholar

Kützing, F. T. (1845). Phycologia germanica, d.i. Deutschlands Algen in bündigen Beschreibungen. W. Köhne Nordhausen.

Google Scholar

Lemoine, F., Correia, D., Lefort, V., Doppelt-Azeroual, O., Mareuil, F., Cohen-Boulakia, S., et al. (2019). NGPhylogeny.fr: new generation phylogenetic services for non-specialists. Nucleic Acids Res. 47, W260–W265. doi: 10.1093/nar/gkz303

PubMed Abstract | Crossref Full Text | Google Scholar

Letunic, I., and Bork, P. (2024). Interactive tree of life (iTOL) v6: recent updates to the phylogenetic tree display and annotation tool. Nucleic Acids Res. 52, W78–W82. doi: 10.1093/nar/gkae268

PubMed Abstract | Crossref Full Text | Google Scholar

Li, H. (2018). Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics 34, 3094–3100. doi: 10.1093/bioinformatics/bty191

PubMed Abstract | Crossref Full Text | Google Scholar

Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., et al. (2009). The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078–2079. doi: 10.1093/bioinformatics/btp352

PubMed Abstract | Crossref Full Text | Google Scholar

Lokhorst, G. M. (1996). Comparative taxonomic studies on the genus Klebsormidium (Charophyceae) in Europe. Cryptogam. Stud. 5, 1–55.

Google Scholar

Lokhorst, G. M., and Star, W. (1985). Ultrastructure of mitosis and cytokinesis in Klebsormidium mucosum nov. comb., formerly Ulothrix verrucosa (Chlorophyta). J. Phycol. 21, 466–476. doi: 10.1111/j.0022-3646.1985.00466.x

Crossref Full Text | Google Scholar

Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 17, 10–12. doi: 10.14806/ej.17.1.200

Crossref Full Text | Google Scholar

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

PubMed Abstract | Crossref Full Text | Google Scholar

McMurdie, P. J., and Holmes, S. (2014). Waste not, want not: why rarefying microbiome data is inadmissible. PLoS Comput. Biol. 10:e1003531. doi: 10.1371/journal.pcbi.1003531

PubMed Abstract | Crossref Full Text | Google Scholar

Medwed, C., Holzinger, A., Hofer, S., Hartmann, A., Michalik, D., Glaser, K., et al. (2021). Ecophysiological, morphological, and biochemical traits of free-living Diplosphaera chodatii (Trebouxiophyceae) reveal adaptation to harsh environmental conditions. Protoplasma 258, 1187–1199. doi: 10.1007/s00709-021-01620-6

PubMed Abstract | Crossref Full Text | Google Scholar

Mikhailyuk, T. I., Glaser, K., Holzinger, A., and Karsten, U. (2015). Biodiversity of Klebsormidium (Streptophyta) from alpine biological soil crusts (Alps, Tyrol, Austria, and Italy). J. Phycol. 51, 750–767. doi: 10.1111/jpy.12316

PubMed Abstract | Crossref Full Text | Google Scholar

Mikhailyuk, T. I., Sluiman, H. J., Massalski, A., Mudimu, O., Demchenko, E. M., Kondratyuk, S. Y., et al. (2008). New streptophyte green algae from terrestrial habitats and an assessment of the genus Interfilum (Klebsormidiophyceae, Streptophyta). J. Phycol. 44, 1586–1603. doi: 10.1111/j.1529-8817.2008.00606.x

PubMed Abstract | Crossref Full Text | Google Scholar

Minh, B. Q., Schmidt, H. A., Chernomor, O., Schrempf, D., Woodhams, M. D., Von Haeseler, A., et al. (2020). IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Molec. Biol. Evol. 37, 1530–1534. doi: 10.1093/molbev/msaa015

PubMed Abstract | Crossref Full Text | Google Scholar

Nawrocki, E. P. (2009). Structural RNA homology search and alignment using covariance models PhD thesis. Washington, DC: Washington University.

Google Scholar

Nurk, S., Meleshko, D., Korobeynikov, A., and Pevzner, P. A. (2017). MetaSPAdes: a new versatile metagenomic assembler. Genome Res. 27, 824–834. doi: 10.1101/gr.213959.116

PubMed Abstract | Crossref Full Text | Google Scholar

Orfanoudaki, M., Hartmann, A., Kamiya, M., West, J., and Ganzera, M. (2020). Chemotaxonomic study of Bostrychia spp. (Ceramiales, Rhodophyta) based on their mycosporine-like amino acid content. Molecules 25:3273. doi: 10.3390/molecules25143273

PubMed Abstract | Crossref Full Text | Google Scholar

Petersen, J. B. (1915). Studier over danske aërophile alger. D. Kgl. Danske Vidensk. Selsk. Skrifter 7, 269–380.

Google Scholar

Prjibelski, A., Antipov, D., Meleshko, D., Lapidus, A., and Korobeynikov, A. (2020). Using SPAdes De novo assembler. Curr. Protoc. Bioinformatics 70:e102. doi: 10.1002/cpbi.102

PubMed Abstract | Crossref Full Text | Google Scholar

R Core Team (2024). R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing.

Google Scholar

Rindi, F. (2007). “Diversity, distribution and ecology of green algae and cyanobacteria in urban habitats” in Algae and Cyanobacteria in extreme environments. ed. J. Seckbach (Netherlands, Dordrecht: Springer), 619–638.

Google Scholar

Rivers, A., Weber, K., Gardner, T., Liu, S., and Armstrong, S. (2018). ITSxpress: software to rapidly trim internally transcribed spacer sequences with quality scores for marker gene analysis [version 1; peer review: 2 approved]. F1000Res 7:15704. doi: 10.12688/f1000research.15704.1

Crossref Full Text | Google Scholar

Rognes, T., Flouri, T., Nichols, B., Quince, C., and Mahé, F. (2016). VSEARCH: a versatile open source tool for metagenomics. PeerJ 4:e2584. doi: 10.7717/peerj.2584

PubMed Abstract | Crossref Full Text | Google Scholar

Rybalka, N., Blanke, M., Tzvetkova, A., Noll, A., Roos, C., Boy, J., et al. (2023). Unrecognized diversity and distribution of soil algae from maritime Antarctica (Fildes peninsula, King George Island). Front. Microbiol. 14:1118747. doi: 10.3389/fmicb.2023.1118747

PubMed Abstract | Crossref Full Text | Google Scholar

Schliep, K. P. (2011). Phangorn: phylogenetic analysis in R. Bioinformatics 27, 592–593. doi: 10.1093/bioinformatics/btq706

PubMed Abstract | Crossref Full Text | Google Scholar

Shen, W., Le, S., Li, Y., and Hu, F. (2016). SeqKit: a cross-platform and ultrafast toolkit for FASTA/Q file manipulation. PLoS One 11:e0163962. doi: 10.1371/journal.pone.0163962

PubMed Abstract | Crossref Full Text | Google Scholar

Søchting, U. (1997). Epiphyllic cover on spruce needles in Denmark. Ann. Bot. Fenn. 34, 157–164.

Google Scholar

Søchting, U., Jensen, B., and Unger, L. (1992). Epifylfloraen på rødgran : Miljøministeriet, Skov- og Naturstyrelsen and University of Copenhagen.

Google Scholar

Starr, R. C., and Zeikus, J. A. (1993). UTEX - the culture collection of algae at the University of Texas at Austin 1993 list of cultures. J. Phycol. 29, 1–106. doi: 10.1111/j.0022-3646.1993.00001.x

Crossref Full Text | Google Scholar

Straub, D., Blackwell, N., Langarica-Fuentes, A., Peltzer, A., Nahnsen, S., and Kleindienst, S. (2020). Interpretations of environmental microbial community studies are biased by the selected 16S rRNA (gene) amplicon sequencing pipeline. Front.Microbiol. 11:550420. doi: 10.3389/fmicb.2020.550420

PubMed Abstract | Crossref Full Text | Google Scholar

Toju, H., Tanabe, A. S., Yamamoto, S., and Sato, H. (2012). High-coverage ITS primers for the DNA-based identification of ascomycetes and basidiomycetes in environmental samples. PLoS One 7:e40863. doi: 10.1371/journal.pone.0040863

PubMed Abstract | Crossref Full Text | Google Scholar

Vilgalys, R., and Hester, M. (1990). Rapid genetic identification and mapping of enzymatically amplified ribosomal DNA from several Cryptococcus species. J. Bacteriol. 172, 4238–4246. doi: 10.1128/jb.172.8.4238-4246.1990

PubMed Abstract | Crossref Full Text | Google Scholar

Vischer, W. (1960). Reproduktion und systematische stellung einiger rinden- und bodenalgen. Schweiz. Z. Hydrol. 22, 330–349. doi: 10.1007/BF02503280

Crossref Full Text | Google Scholar

Volkmann, M., and Gorbushina, A. A. (2006). A broadly applicable method for extraction and characterization of mycosporines and mycosporine-like amino acids of terrestrial, marine and freshwater origin. FEMS Microbiol. Lett. 255, 286–295. doi: 10.1111/j.1574-6968.2006.00088.x

PubMed Abstract | Crossref Full Text | Google Scholar

White, T. J., Bruns, T., Lee, S., and Taylor, J. W. (1990). “Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics” in PCR protocols: A guide to methods and applications. eds. M. A. Innis, D. H. Gelfand, J. J. Snisky, and T. J. White (New York: Academic Press), 315–322.

Google Scholar

Wickham, H., Averick, M., Bryan, J., Chang, W., Mcgowan, L. D., François, R., et al. (2019). Welcome to the tidyverse. J. Open Source Softw. 4:1686. doi: 10.21105/joss.01686

Crossref Full Text | Google Scholar

Wolf, M., Koetschan, C., and Müller, T. (2014). ITS2, 18S, 16S or any other RNA – simply aligning sequences and their individual secondary structures simultaneously by an automatic approach. Gene 546, 145–149. doi: 10.1016/j.gene.2014.05.065

PubMed Abstract | Crossref Full Text | Google Scholar

Yancey, P. H. (2005). Organic osmolytes as compatible, metabolic and counteracting cytoprotectants in high osmolarity and other stresses. J. Exp. Biol. 208, 2819–2830. doi: 10.1242/jeb.01730

PubMed Abstract | Crossref Full Text | Google Scholar

Zahradníková, M., Andersen, H. L., Tønsberg, T., and Beck, A. (2017). Molecular evidence of Apatococcus, including a. fuscideae sp. nov., as photobiont in the genus. Fuscidea Protist 168, 425–438. doi: 10.1016/j.protis.2017.06.002

PubMed Abstract | Crossref Full Text | Google Scholar

Zeugin, J. A., and Hartley, J. L. (1985). Ethanol precipitation of DNA. Focus 7, 1–2.

Google Scholar

Zuo, Z., Cao, L., Nothia, L.-F., and Mohimani, H. (2021). MS2Planner: improved fragmentation spectra coverage in untargeted mass spectrometry by iterative optimized data acquisition. Bioinformatics 37, i231–i236. doi: 10.1093/bioinformatics/btab279

PubMed Abstract | Crossref Full Text | Google Scholar

Zwerger, M., and Ganzera, M. (2022). Fast and efficient separation of eleven mycosporine-like amino acids by UHPLC-DAD and their quantification in diverse red algae. Mar. Drugs 20:395. doi: 10.3390/md20060395

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: terrestrial epiphytic, amplicon-based metabarcoding, bioindication, metagenomics, Hormidium crenulatum, Klebsormidium, phylogenetic analyses, mycosporine-like amino acids (MAA)

Citation: Søchting U, Friedl T, Moestrup Ø, Grewe F, Sun Y, Çakır YT, Ganzera M, Glaser K, Heesch S, Hammerle F, Nimptsch D, Olberg B and Karsten U (2025) Mass development of a filamentous and likely nitrophilous aerophytic green alga on tree bark: Apatococcus ammoniophilus sp. nov. (Chlorophyta, Trebouxiophyceae). Front. Microbiol. 16:1633308. doi: 10.3389/fmicb.2025.1633308

Received: 23 May 2025; Accepted: 02 July 2025;
Published: 23 July 2025.

Edited by:

Muhammad Zahid Mumtaz, Gansu Agricultural University, China

Reviewed by:

Abdiel Jover Capote, Universidad de Oriente Santiago de Cuba, Cuba
Zain Ul Arifeen Muhammad, Nanjing University, China

Copyright © 2025 Søchting, Friedl, Moestrup, Grewe, Sun, Çakır, Ganzera, Glaser, Heesch, Hammerle, Nimptsch, Olberg and Karsten. 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: Ulrik Søchting, dWxyaWtzQGJpby5rdS5kaw==

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.