On the use of high-throughput sequencing for the study of cyanobacterial diversity in Antarctic aquatic mats

Spatiotemporal diversity and community structure of cyanobacteria and associated bacteria in the large shallow subtropical Lake Okeechobee (Florida, United States)

Lake Okeechobee is a large eutrophic, shallow, subtropical lake in south Florida, United States. Due to decades of nutrient loading and phosphorus rich sediments, the lake is eutrophic and frequently experiences cyanobacterial harmful algal blooms (cyanoHABs). In the past, surveys of the phytoplankton community structure in the lake have been conducted by morphological studies, whereas molecular based studies have been seldom employed. With increased frequency of cyanoHABs in Lake Okeechobee (e.g., 2016 and 2018 Microcystis-dominated blooms), it is imperative to determine the diversity of cyanobacterial taxa that exist within the lake and the limnological parameters that drive bloom-forming genera. A spatiotemporal study of the lake was conducted over the course of 1 year to characterize the (cyano)bacterial community structure, using 16S rRNA metabarcoding, with coincident collection of limnological parameters (e.g., nutrients, water temperature, major ions), and cyanotoxins. The objectives of this study were to elucidate spatiotemporal trends of community structure, identify drivers of community structure, and examine cyanobacteria-bacterial relationships within the lake. Results indicated that cyanobacterial communities within the lake were significantly different between the wet and dry season, but not between periods of nitrogen limitation and co-nutrient limitation. Throughout the year, the lake was primarily dominated by the picocyanobacterium Cyanobium. The bloom-forming genera Cuspidothrix, Dolichospermum, Microcystis, and Raphidiopsis were highly abundant throughout the lake and had disparate nutrient requirements and niches within the lake. Anatoxin-a, microcystins, and nodularins were detected throughout the lake across both seasons. There were no correlated (cyano)bacteria shared between the common bloom-forming cyanobacteria Dolichospermum, Microcystis, and Raphidiopsis. This study is the first of its kind to use molecular based methods to assess the cyanobacterial community structure within the lake. These data greatly improve our understanding of the cyanobacterial community structure within the lake and the physiochemical parameters which may drive the bloom-forming taxa within Lake Okeechobee.

Novel diversity of polar Cyanobacteria revealed by genome-resolved metagenomics

Benthic microbial mats dominated by Cyanobacteria are important features of polar lakes. Although culture-independent studies have provided important insights into the diversity of polar Cyanobacteria, only a handful of genomes have been sequenced to date. Here, we applied a genome-resolved metagenomics approach to data obtained from Arctic, sub-Antarctic and Antarctic microbial mats. We recovered 37 metagenome-assembled genomes (MAGs) of Cyanobacteria representing 17 distinct species, most of which are only distantly related to genomes that have been sequenced so far. These include (i) lineages that are common in polar microbial mats such as the filamentous taxa Pseudanabaena, Leptolyngbya, Microcoleus/Tychonema and Phormidium; (ii) the less common taxa Crinalium and Chamaesiphon; (iii) an enigmatic Chroococcales lineage only distantly related to Microcystis; and (iv) an early branching lineage in the order Gloeobacterales that is distributed across the cold biosphere, for which we propose the name Candidatus Sivonenia alaskensis. Our results show that genome-resolved metagenomics is a powerful tool for expanding our understanding of the diversity of Cyanobacteria, especially in understudied remote and extreme environments.

Cyanoseq: A database of cyanobacterial 16S rRNA gene sequences with curated taxonomy

Cyanobacteria are photosynthetic bacteria that occupy various habitats across the globe, playing critical roles in many of Earth's biogeochemical cycles both in both aquatic and terrestrial systems. Despite their well-known significance, their taxonomy remains problematic and is the subject of much research. Taxonomic issues of Cyanobacteria have consequently led to inaccurate curation within known reference databases, ultimately leading to problematic taxonomic assignment during diversity studies. Recent advances in sequencing technologies have increased our ability to characterize and understand microbial communities, leading to the generation of thousands of sequences that require taxonomic assignment. We herein propose CyanoSeq (https://zenodo.org/record/7569105), a database of cyanobacterial 16S rRNA gene sequences with curated taxonomy. The taxonomy of CyanoSeq is based on the current state of cyanobacterial taxonomy, with ranks from the domain to genus level. Files are provided for use with common naive Bayes taxonomic classifiers, such as those included in DADA2 or the QIIME2 platform. Additionally, FASTA files are provided for creation of de novo phylogenetic trees with (near) full-length 16S rRNA gene sequences to determine the phylogenetic relationship of cyanobacterial strains and/or ASV/OTUs. The database currently consists of 5410 cyanobacterial 16S rRNA gene sequences along with 123 Chloroplast, Bacterial, and Vampirovibrionia (formally Melainabacteria) sequences.

Metagenomic mapping of cyanobacteria and potential cyanotoxin producing taxa in large rivers of the United States

Cyanobacteria and cyanotoxin producing cyanobacterial blooms are a trending focus of current research. Many studies focus on bloom events in lentic environments such as lakes or ponds. Comparatively few studies have explored lotic environments and fewer still have examined the cyanobacterial communities and potential cyanotoxin producers during ambient, non-bloom conditions. Here we used a metagenomics-based approach to profile non-bloom microbial communities and cyanobacteria in 12 major U.S. rivers at multiple time points during the summer months of 2019. Our data show that U.S. rivers possess microbial communities that are taxonomically rich, yet largely consistent across geographic location and time. Within these communities, cyanobacteria often comprise significant portions and frequently include multiple species with known cyanotoxin producing strains. We further characterized these potential cyanotoxin producing taxa by deep sequencing amplicons of the microcystin E (mcyE) gene. We found that rivers containing the highest levels of potential cyanotoxin producing cyanobacteria consistently possess taxa with the genetic potential for cyanotoxin production and that, among these taxa, the predominant genus of origin for the mcyE gene is Microcystis. Combined, these data provide a unique perspective on cyanobacteria and potential cyanotoxin producing taxa that exist in large rivers across the U.S. and can be used to better understand the ambient conditions that may precede bloom events in lotic freshwater ecosystems.

Population structure of an Antarctic aquatic cyanobacterium

BACKGROUND

Ace Lake is a marine-derived, stratified lake in the Vestfold Hills of East Antarctica with an upper oxic and lower anoxic zone. Cyanobacteria are known to reside throughout the water column. A Synechococcus-like species becomes the most abundant member in the upper sunlit waters during summer while persisting annually even in the absence of sunlight and at depth in the anoxic zone. Here, we analysed ~ 300 Gb of Ace Lake metagenome data including 59 Synechococcus-like metagenome-assembled genomes (MAGs) to determine depth-related variation in cyanobacterial population structure. Metagenome data were also analysed to investigate viruses associated with this cyanobacterium and the host's capacity to defend against or evade viruses.

RESULTS

A single Synechococcus-like species was found to exist in Ace Lake, Candidatus Regnicoccus frigidus sp. nov., consisting of one phylotype more abundant in the oxic zone and a second phylotype prevalent in the oxic-anoxic interface and surrounding depths. An important aspect of genomic variation pertained to nitrogen utilisation, with the capacity to perform cyanide assimilation and asparagine synthesis reflecting the depth distribution of available sources of nitrogen. Both specialist (host specific) and generalist (broad host range) viruses were identified with a predicted ability to infect Ca. Regnicoccus frigidus. Host-virus interactions were characterised by a depth-dependent distribution of virus type (e.g. highest abundance of specialist viruses in the oxic zone) and host phylotype capacity to defend against (e.g. restriction-modification, retron and BREX systems) and evade viruses (cell surface proteins and cell wall biosynthesis and modification enzymes).

CONCLUSION

In Ace Lake, specific environmental factors such as the seasonal availability of sunlight affects microbial abundances and the associated processes that the microbial community performs. Here, we find that the population structure for Ca. Regnicoccus frigidus has evolved differently to the other dominant phototroph in the lake, Candidatus Chlorobium antarcticum. The geography (i.e. Antarctica), limnology (e.g. stratification) and abiotic (e.g. sunlight) and biotic (e.g. microbial interactions) factors determine the types of niches that develop in the lake. While the lake community has become increasingly well studied, metagenome-based studies are revealing that niche adaptation can take many paths; these paths need to be determined in order to make reasonable predictions about the consequences of future ecosystem perturbations. Video Abstract.

A Summer of Cyanobacterial Blooms in Belgian Waterbodies: Microcystin Quantification and Molecular Characterizations

In the context of increasing occurrences of toxic cyanobacterial blooms worldwide, their monitoring in Belgium is currently performed by regional environmental agencies (in two of three regions) using different protocols and is restricted to some selected recreational ponds and lakes. Therefore, a global assessment based on the comparison of existing datasets is not possible. For this study, 79 water samples from a monitoring of five lakes in Wallonia and occasional blooms in Flanders and Brussels, including a canal, were analyzed. A Liquid Chromatography with tandem mass spectrometry (LC-MS/MS) method allowed to detect and quantify eight microcystin congeners. The mcyE gene was detected using PCR, while dominant cyanobacterial species were identified using 16S RNA amplification and direct sequencing. The cyanobacterial diversity for two water samples was characterized with amplicon sequencing. Microcystins were detected above limit of quantification (LOQ) in 68 water samples, and the World Health Organization (WHO) recommended guideline value for microcystins in recreational water (24 µg L⁻¹) was surpassed in 18 samples. The microcystin concentrations ranged from 0.11 µg L⁻¹ to 2798.81 µg L⁻¹ total microcystin. For 45 samples, the dominance of the genera Microcystis sp., Dolichospermum sp., Aphanizomenon sp., Cyanobium/Synechococcus sp., Planktothrix sp., Romeria sp., Cyanodictyon sp., and Phormidium sp. was shown. Moreover, the mcyE gene was detected in 75.71% of all the water samples.

Temporal Patterns of Bacterial and Viral Communities during Algae Blooms of a Reservoir in Macau

Compositions of microbial communities associated with blooms of algae in a storage reservoir in Macau, China were investigated between 2013 and 2016. Algae were enumerated by visible light microscopy. Profiles of organisms in water were examined by 16S rRNA sequences and viral metagenomics, based on next generation sequencing. Results of 16S rRNA sequencing indicated that majority of the identified organisms were bacteria closely related to Proteobacteria, Cyanobacteria, Verrucomicrobia, Bacteroidetes, and Actinobacteria. Metagenomics sequences demonstrated that the dominant virus was Phycodnavirus, accounting for 70% of the total population. Patterns of relative numbers of bacteria in the microbial community and their temporal changes were determined through alpha diversity indices, principal coordinates analysis (PCoA), relative abundance, and visualized by Venn diagrams. Ways in which the bacterial and viral communities are influenced by various water-related variables were elucidated based on redundancy analysis (RDA). Relationships of the relative numbers of bacteria with trophic status in a reservoir used for drinking water in Macau, provided insight into associations of Phycodnavirus and Proteobacteria with changes in blooms of algae.

Diversity Assessment of Toxic Cyanobacterial Blooms during Oxidation

Fresh-water sources of drinking water are experiencing toxic cyanobacterial blooms more frequently. Chemical oxidation is a common approach to treat cyanobacteria and their toxins. This study systematically investigates the bacterial/cyanobacterial community following chemical oxidation (Cl₂, KMnO₄, O₃, H₂O₂) using high throughput sequencing. Raw water results from high throughput sequencing show that Proteobacteria, Actinobacteria, Cyanobacteria and Bacteroidetes were the most abundant phyla. Dolichospermum, Synechococcus, Microcystis and Nostoc were the most dominant genera. In terms of species, Dolichospermum sp.90 and Microcystis aeruginosa were the most abundant species at the beginning and end of the sampling, respectively. A comparison between the results of high throughput sequencing and taxonomic cell counts highlighted the robustness of high throughput sequencing to thoroughly reveal a wide diversity of bacterial and cyanobacterial communities. Principal component analysis of the oxidation samples results showed a progressive shift in the composition of bacterial/cyanobacterial communities following soft-chlorination with increasing common exposure units (CTs) (0–3.8 mg·min/L). Close cyanobacterial community composition (Dolichospermum dominant genus) was observed following low chlorine and mid-KMnO₄ (287.7 mg·min/L) exposure. Our results showed that some toxin producing species may persist after oxidation whether they were dominant species or not. Relative persistence of Dolichospermum sp.90 was observed following soft-chlorination (0.2–0.6 mg/L) and permanganate (5 mg/L) oxidation with increasing oxidant exposure. Pre-oxidation using H₂O₂ (10 mg/L and one day contact time) caused a clear decrease in the relative abundance of all the taxa and some species including the toxin producing taxa. These observations suggest selectivity of H₂O₂ to provide an efficient barrier against toxin producing cyanobacteria entering a water treatment plant.

Toxic benthic freshwater cyanobacterial proliferations: Challenges and solutions for enhancing knowledge and improving monitoring and mitigation

1. This review summarises knowledge on the ecology, toxin production, and impacts of toxic freshwater benthic cyanobacterial proliferations. It documents monitoring, management, and sampling strategies, and explores mitigation options. 2. Toxic proliferations of freshwater benthic cyanobacteria (taxa that grow attached to substrates) occur in streams, rivers, lakes, and thermal and meltwater ponds, and have been reported in 19 countries. Anatoxin- and microcystin-containing mats are most commonly reported (eight and 10 countries, respectively). 3. Studies exploring factors that promote toxic benthic cyanobacterial proliferations are limited to a few species and habitats. There is a hierarchy of importance in environmental and biological factors that regulate proliferations with variables such as flow (rivers), fine sediment deposition, nutrients, associated microbes, and grazing identified as key drivers. Regulating factors differ among colonisation, expansion, and dispersal phases. 4. New -omics-based approaches are providing novel insights into the physiological attributes of benthic cyanobacteria and the role of associated microorganisms in facilitating their proliferation. 5. Proliferations are commonly comprised of both toxic and non-toxic strains, and the relative proportion of these is the key factor contributing to the overall toxin content of each mat. 6. While these events are becoming more commonly reported globally, we currently lack standardised approaches to detect, monitor, and manage this emerging health issue. To solve these critical gaps, global collaborations are needed to facilitate the rapid transfer of knowledge and promote the development of standardised techniques that can be applied to diverse habitats and species, and ultimately lead to improved management.

Latitudinal gradient of cyanobacterial diversity in tidal flats

Latitudinal diversity gradients are well-known for plants and animals, but only recently similar patterns have been described for some specific microbial communities in distinct habitats. Although microbial diversity is well-investigated worldwide, most of the studies are spatially too restricted to allow general statements about global diversity patterns. Additionally, methodological differences make it hard and often impossible to compare several studies. This study investigated the cyanobacterial diversity in tidal flats along geographical and ecological gradients based on high-throughput sequencing of 16S rRNA gene fragments (Illumina MiSeq) and environmental data on a large spatial scale from the subtropics to the Arctic. Latitude and strongly correlated environmental parameters (e.g. temperature) were identified as important drivers of cyanobacterial diversity on global scale resulting in a latitudinal diversity gradient similar to that known from plants and animals. Other non-correlated parameters (e.g. grain size) were shown to be more important on local scales, although no consistent pattern occurred across different locations. Among a total number of 989 operational taxonomic units (OTUs) only one cosmopolitan (classified as Coleofasciculus chthonoplastes), but many location-specific and putative endemic ones (78%) were detected. High proportions of rare members of the community (up to 86%) were found in all samples. Phylogenetic beta diversity was shown to be influenced by the developmental stage of the mat community becoming increasingly similar with increasing stabilization.

Unmasking the identity of toxigenic cyanobacteria driving a multi-toxin bloom by high-throughput sequencing of cyanotoxins genes and 16S rRNA metabarcoding

Cyanobacterial harmful algal blooms (CyanoHABs) are complex communities that include coexisting toxic and non-toxic strains only distinguishable by genetic methods. This study shows a water-management oriented use of next generation sequencing (NGS) to specifically pinpoint toxigenic cyanobacteria within a bloom simultaneously containing three of the most widespread cyanotoxins (the hepatotoxins microcystins, MCs; and the neurotoxins anatoxin-a, ATX, and saxitoxins, STXs). The 2013 summer bloom in Rosarito reservoir (Spain) comprised 33 cyanobacterial OTUs based on 16S rRNA metabarcoding, 7 of which accounted for as much as 96.6% of the community. Cyanotoxins and their respective biosynthesis genes were concurrently present throughout the entire bloom event including: MCs and mcyE gene; ATX and anaF gene; and STXs and sxtI gene. NGS applied to amplicons of cyanotoxin-biosynthesis genes unveiled 6 toxigenic OTUs, comprising 3 involved in MCs production (Planktothrix agardhii and 2 Microcystis spp.), 2 in ATX production (Cuspidothrix issatschenkoi and Phormidium/Tychonema spp.) and 1 in STXs production (Aphanizomenon gracile). These toxigenic taxa were also present in 16S rRNA OTUs list and their relative abundance was positively correlated with the respective toxin concentrations. Our results point at MC-producing P. agardhii and ATX-producing C. issatschenkoi as the main contributors to the moderate toxin concentrations observed, and suggest that their distribution in Southern Europe is broader than previously thought. Our findings also stress the need for monitoring low-abundance cyanobacteria (<1% relative abundance) in cyanotoxicity studies, and provide novel data on the presence of picocyanobacteria and potentially ATX-producing benthic taxa (e.g., Phormidium) in deep thermally-stratified water bodies. This study showcases a straightforward use of amplicon metagenomics of cyanotoxin biosynthesis genes in a multi-toxin bloom thus illustrating the broad applicability of NGS for water management in risk-oriented monitoring of CyanoHABs.

Marked Succession of Cyanobacterial Communities Following Glacier Retreat in the High Arctic

Cyanobacteria are important colonizers of recently deglaciated proglacial soil but an in-depth investigation of cyanobacterial succession following glacier retreat has not yet been carried out. Here, we report on the successional trajectories of cyanobacterial communities in biological soil crusts (BSCs) along a 100-year deglaciation gradient in three glacier forefields in central Svalbard, High Arctic. Distance from the glacier terminus was used as a proxy for soil age (years since deglaciation), and cyanobacterial abundance and community composition were evaluated by epifluorescence microscopy and pyrosequencing of partial 16S rRNA gene sequences, respectively. Succession was characterized by a decrease in phylotype richness and a marked shift in community structure, resulting in a clear separation between early (10-20 years since deglaciation), mid (30-50 years), and late (80-100 years) communities. Changes in cyanobacterial community structure were mainly connected with soil age and associated shifts in soil chemical composition (mainly moisture, SOC, SMN, K, and Na concentrations). Phylotypes associated with early communities were related either to potentially novel lineages (< 97.5% similar to sequences currently available in GenBank) or lineages predominantly restricted to polar and alpine biotopes, suggesting that the initial colonization of proglacial soil is accomplished by cyanobacteria transported from nearby glacial environments. Late communities, on the other hand, included more widely distributed genotypes, which appear to establish only after the microenvironment has been modified by the pioneering taxa.

Diversity and Cyclical Seasonal Transitions in the Bacterial Community in a Large and Deep Perialpine Lake

High-throughput sequencing (HTS) was used to analyze the seasonal variations in the bacterioplankton community composition (BCC) in the euphotic layer of a large and deep lake south of the Alps (Lake Garda). The BCC was analyzed throughout two annual cycles by monthly samplings using the amplification and sequencing of the V3-V4 hypervariable region of the 16S rRNA gene by the MiSeq Illumina platform. The dominant and most diverse bacterioplankton phyla were among the more frequently reported in freshwater ecosystems, including the Proteobacteria, Cyanobacteria, Bacteroidetes, Verrucomicrobia, Actinobacteria, and Planctomycetes. As a distinctive feature, the development of the BCC showed a cyclical temporal pattern in the two analyzed years and throughout the euphotic layer. The recurring temporal development was controlled by the strong seasonality in water temperature and thermal stratification, and by cyclical temporal changes in nutrients and, possibly, by the remarkable annual cyclical development of cyanobacteria and eukaryotic phytoplankton hosting bacterioplankton that characterizes Lake Garda. Further downstream analyses of operational taxonomic units associated to cyanobacteria allowed confirming the presence of the most abundant taxa previously identified by microscopy and/or phylogenetic analyses, as well as the presence of other small Synechococcales/Chroococcales and rare Nostocales never identified so far in the deep lakes south of the Alps. The implications of the high diversity and strong seasonality are relevant, opening perspectives for the definition of common and discriminating patterns characterizing the temporal and spatial distribution in the BCC, and for the application of the new sequencing technologies in the monitoring of water quality in large and deep lakes.

Genomic variation and biogeography of Antarctic haloarchaea

BACKGROUND

The genomes of halophilic archaea (haloarchaea) often comprise multiple replicons. Genomic variation in haloarchaea has been linked to viral infection pressure and, in the case of Antarctic communities, can be caused by intergenera gene exchange. To expand understanding of genome variation and biogeography of Antarctic haloarchaea, here we assessed genomic variation between two strains of Halorubrum lacusprofundi that were isolated from Antarctic hypersaline lakes from different regions (Vestfold Hills and Rauer Islands). To assess variation in haloarchaeal populations, including the presence of genomic islands, metagenomes from six hypersaline Antarctic lakes were characterised.

RESULTS

The sequence of the largest replicon of each Hrr. lacusprofundi strain (primary replicon) was highly conserved, while each of the strains' two smaller replicons (secondary replicons) were highly variable. Intergenera gene exchange was identified, including the sharing of a type I-B CRISPR system. Evaluation of infectivity of an Antarctic halovirus provided experimental evidence for the differential susceptibility of the strains, bolstering inferences that strain variation is important for modulating interactions with viruses. A relationship was found between genomic structuring and the location of variation within replicons and genomic islands, demonstrating that the way in which haloarchaea accommodate genomic variability relates to replicon structuring. Metagenome read and contig mapping and clustering and scaling analyses demonstrated biogeographical patterning of variation consistent with environment and distance effects. The metagenome data also demonstrated that specific haloarchaeal species dominated the hypersaline systems indicating they are endemic to Antarctica.

CONCLUSION

The study describes how genomic variation manifests in Antarctic-lake haloarchaeal communities and provides the basis for future assessments of Antarctic regional and global biogeography of haloarchaea.

Toxic Cyanobacteria in Svalbard: Chemical Diversity of Microcystins Detected Using a Liquid Chromatography Mass Spectrometry Precursor Ion Screening Method

Cyanobacteria synthesize a large variety of secondary metabolites including toxins. Microcystins (MCs) with hepato- and neurotoxic potential are well studied in bloom-forming planktonic species of temperate and tropical regions. Cyanobacterial biofilms thriving in the polar regions have recently emerged as a rich source for cyanobacterial secondary metabolites including previously undescribed congeners of microcystin. However, detection and detailed identification of these compounds is difficult due to unusual sample matrices and structural congeners produced. We here report a time-efficient liquid chromatography-mass spectrometry (LC-MS) precursor ion screening method that facilitates microcystin detection and identification. We applied this method to detect six different MC congeners in 8 out of 26 microbial mat samples of the Svalbard Archipelago in the Arctic. The congeners, of which [Asp³, ADMAdda⁵, Dhb⁷] MC-LR was most abundant, were similar to those reported in other polar habitats. Microcystins were also determined using an Adda-specific enzyme-linked immunosorbent assay (Adda-ELISA). Nostoc sp. was identified as a putative toxin producer using molecular methods that targeted 16S rRNA genes and genes involved in microcystin production. The mcy genes detected showed highest similarities to other Arctic or Antarctic sequences. The LC-MS precursor ion screening method could be useful for microcystin detection in unusual matrices such as benthic biofilms or lichen.

Allelopathy-mediated competition in microbial mats from Antarctic lakes

Microbial mats are vertically stratified communities that host a complex consortium of microorganisms, dominated by cyanobacteria, which compete for available nutrients and environmental niches, within these extreme habitats. The Antarctic Dry Valleys near McMurdo Sound include a series of lakes within the drainage basin that are bisected by glacial traverses. These lakes are traditionally independent, but recent increases in glacial melting have allowed two lakes (Chad and Hoare) to become connected by a meltwater stream. Microbial mats were collected from these lakes, and cultured under identical conditions at the McMurdo Station laboratory. Replicate pairings of the microbial mats exhibited consistent patterns of growth inhibition indicative of competitive dominance. Natural products were extracted from the microbial mats, and a disk diffusion assay was utilized to show that allelochemical compounds mediate competitive interactions. Both microscopy and 16S rRNA sequencing show that these mats contain significant populations of cyanobacteria known to produce allelochemicals. Two compounds were isolated from these microbial mats that might be important in the chemical ecology of these psychrophiles. In other disk:mat pairings, including extract versus mat of origin, the allelochemicals exhibited no effect. Taken together, these results indicate that Antarctic lake microbial mats can compete via allelopathy.

Unexpected associated microalgal diversity in the lichen Ramalina farinacea is uncovered by pyrosequencing analyses

The current literature reveals that the intrathalline coexistence of multiple microalgal taxa in lichens is more common than previously thought, and additional complexity is supported by the coexistence of bacteria and basidiomycete yeasts in lichen thalli. This replaces the old paradigm that lichen symbiosis occurs between a fungus and a single photobiont. The lichen Ramalina farinacea has proven to be a suitable model to study the multiplicity of microalgae in lichen thalli due to the constant coexistence of Trebouxia sp. TR9 and T. jamesii in long-distance populations. To date, studies involving phycobiont diversity within entire thalli are based on Sanger sequencing, but this method seems to underestimate the diversity. Here, we aim to analyze both the microalgal diversity and its community structure in a single thallus of the lichen R. farinacea by applying a 454 pyrosequencing approach coupled with a careful ad hoc-performed protocol for lichen sample processing prior to DNA extraction. To ascertain the reliability of the pyrosequencing results and the applied bioinformatics pipeline results, the thalli were divided into three sections (apical, middle and basal zones), and a mock community sample was used. The developed methodology allowed 40448 filtered algal reads to be obtained from a single lichen thallus, which encompassed 31 OTUs representative of different microalgae genera. In addition to corroborating the coexistence of the two Trebouxia sp. TR9 and T. jamesii taxa in the same thallus, this study showed a much higher microalgal diversity associated with the lichen. Along the thallus ramifications, we also detected variations in phycobiont distribution that might correlate with different microenvironmental conditions. These results highlight R. farinacea as a suitable material for studying microalgal diversity and further strengthen the concept of lichens as multispecies microecosystems. Future analyses will be relevant to ecophysiological and evolutionary studies to understand the roles of the multiple photobionts in lichen symbioses.