Disentangling the drivers of functional complexity at the metagenomic level in Shark Bay microbial mat microbiomes

The functional metagenomic potential of Shark Bay microbial mats was examined for the first time at a millimeter scale, employing shotgun sequencing of communities via the Illumina NextSeq 500 platform in conjunction with defined chemical analyses. A detailed functional metagenomic profile has elucidated key pathways and facilitated inference of critical microbial interactions. In addition, 87 medium-to-high-quality metagenome-assembled genomes (MAG) were assembled, including potentially novel bins under the deep-branching archaeal Asgard group (Thorarchaetoa and Lokiarchaeota). A range of pathways involved in carbon, nitrogen, sulfur, and phosphorus cycles were identified in mat metagenomes, with the Wood–Ljungdahl pathway over-represented and inferred as a major carbon fixation mode. The top five sets of genes were affiliated to sulfate assimilation (cysNC cysNCD, sat), methanogenesis (hdrABC), Wood–Ljungdahl pathways (cooS, coxSML), phosphate transport (pstB), and copper efflux (copA). Polyhydroxyalkanoate (PHA) synthase genes were over-represented at the surface, with PHA serving as a potential storage of fixed carbon. Sulfur metabolism genes were highly represented, in particular complete sets of genes responsible for both assimilatory and dissimilatory sulfate reduction. Pathways of environmental adaptation (UV, hypersalinity, oxidative stress, and heavy metal resistance) were also delineated, as well as putative viral defensive mechanisms (core genes of the CRISPR, BREX, and DISARM systems). This study provides new metagenome-based models of how biogeochemical cycles and adaptive responses may be partitioned in the microbial mats of Shark Bay.