Data from: Selective yet opportunistic assemblies of symbionts in a lichen community, Cladia asiatica, in Taiwan
Data files
Mar 23, 2026 version files 482.98 KB
-
Kuan_Cladia_asiatica.zip
479.41 KB
-
README.md
3.58 KB
Abstract
Lichens are ecological communities composed primarily of fungal (mycobiont) and algal or cyanobacterial (photobiont) partners. However, the evolutionary dynamics underlying these multi-partner associations remain poorly understood. Using Cladia asiatica lichens from Taiwan, we investigated the phylogeography and co-evolutionary relationships among mycobionts and their putative photobionts. Our results revealed that fungal hosts exhibit strong global-scale phylogenetic clustering among regional lineages, supporting the recent taxonomic treatment of an Asian endemic clade, but showing little fine-scale phylogeographic structure within Taiwan. In contrast, green algae include multiple evolutionarily distinct lineages coexisting in Taiwan and display statistically significant spatial associations that may imply that some geographic or environmental filtering factors may have determined their distribution. Furthermore, weak cophylogenetic congruence between hosts and photobionts indicates that C. asiatica lichens are locally assembled rather than dispersing as integrated evolutionary units. We argue that each symbiotic partner follows independent evolutionary and biogeographic trajectories. Importantly, algal diversity appears to shape the ecological breadth of lichens, implying that shifts in algal distributions may constrain lichen persistence under climate change. Our findings reveal the need to integrate photobiont diversity into evolutionary frameworks to better understand how lichens function as bioindicators and how their biodiversity may respond to global environmental change.
Dataset DOI: 10.5061/dryad.2jm63xt4b
Description of the data and file structure
The associated sequence alignment data and the reconstructed phylogenies of the study.
The files also include the geographic coordinates of the sampling localities and the associated bioclim variables retrieved from worldclim data base
Files and variables
File: Kuan_Cladia_asiatica.zip
Description: A compressed file including folders that contain different sets of data.
After decompressed:
BEAST trees
-MCC_trees
--Taiwan folder: including four folders, each containing a MCC tree output from BEAST analysis (.tree file) and an input file (.xml format) for BEAST analysis. 16s beast_taiwan: result from 16s cyanobacterial dataset; ITS1A beast_taiwan: associated result and input for green algal ITS1 dataset; ITS1F beast_taiwan: associated result and input for fungal ITS1 dataset; rbcL beast_taiwan: associated result and input for green algal rbcL dataset.
--Global folder: including four folders, each containing a MCC tree output from BEAST analysis (.tree file) and an input file (.xml format) for BEAST analysis. 16s beast: result from 16s cyanobacterial dataset; ITS1A beast: associated result and input for green algal ITS1 dataset; ITS1F beast: associated result and input for fungal ITS1 dataset; rbcL beast: associated result and input for green algal rbcL dataset.
PhyML trees
-PhyML_outputs
--Taiwan folder: including the maximum likelihood trees (.txt files) and the smart model selection results (.csv files) from the analyses. 16s phyml_taiwan: results from 16s cyanobacterial dataset; its1a phyml_taiwan: associated results for green algal ITS1 dataset; ITS1F phyml_taiwan: associated results for fungal ITS1 dataset; rbcL phyml_taiwan: associated results for green algal rbcL dataset.
For the smart model section, the results are in the CSV files:
- Model: The nucleotide substitution model tested (e.g., HKY and GTR), which describes the rates of substitution between nucleotides.
- Decoration: Additional model parameters indicating rate heterogeneity and base frequency assumptions (e.g.,+G for gamma-distribution rate variation among sites, +I for proportion of invariant sites).
- K: The number of free parameters in the model.
- Llk: The log-likelihood of the model given the sequence alignment data; higher values indicate a better fit.
- AIC: Akaike information criterion, calculated from the log-likelihood and number of parameters; lower values indicate a better fit.
- BIC: Bayesian information criterion, calculated from the log-likelihood and number of parameters; lower values indicate a better fit.
sequence alignments
The sequence alignments for the four loci used in the study are in fasta format. There are two subfolders, one includes the sequence alignments from Taiwanese samples only (Taiwan_samples), while the other one includes global sampling (Global_samples).
geo coordinates
The geographic coordinate for each sample included in each sequence alignment. The different columns denote the longitude and latitude, as well as elevation (if available) associated with each sampled individual (taxa). The coordinates data can be used to retrieve climatic variables from the WorldClim data based as described in the study.
Access information
Other publicly accessible locations of the data:
- GenBank
