Fig. 2

A Beta diversity patterns for hot spring photosynthetic biofilms (N = 395) indicated statistically supported biogeographic regions. Beta diversity was estimated using Principal Coordinate Analysis of unweighted UniFRAC distances (Confidence interval 99%), where NT = North Thailand, CT = Central Thailand, ST = South Thailand, NM = North Malaysia, SM = South Malaysia, SG = Singapore. Symbol shapes denote hot spring location codes (see Supplementary Table S1). B A Differential heat tree matrix showing: hierarchal taxonomic coverage represented as nodes, number of ASVs assigned to each taxon as node size/width, and significant differences in pairwise comparisons of median proportions of reads between biogeographic regions as colored nodes determined using a Wilcox rank-sum test. Bacterial and archaeal lineages that were relatively enriched in biofilms are depicted as colored regions in the rows (green branches and nodes) and columns (brown branches and nodes). The labelled gray tree on the lower left is a key for the smaller unlabeled trees in the matrix. Additionally, a heat tree showing all levels of taxonomic coverage with all nodes labeled is shown in Supplementary Fig. S6. C) Functional profiling of hot spring photosynthetic biofilms using FAPROTAX revealed all biofilm samples were dominated by oxygenic photoautotrophy and aerobic chemoheterotrophy