Microbial communities in sandy beaches from the three domains of life differ by microhabitat and intertidal location

The microbial communities of sandy beaches are poorly described despite the biogeochemical importance and ubiquity of these ecosystems. Using metabarcoding of the 16S and 18S rRNA genes, we investigated the diversity, microhabitats (with or between sand grains) and intertidal distributions of microorganisms (including meiofauna) from pristine sandy beaches in British Columbia, Canada, and hypothesized that abiotic variations due to microhabitat or intertidal gradients influence the distribution of microorganisms on local scales. Bacterial, archaeal and protistan communities of the sand were clearly distinct from interstitial communities, and from planktonic communities of the overlying seawater, which correlated with differences in function and lifestyle (e.g., sulphur reduction and gliding motility). In contrast, meiofaunal communities could not be distinguished by sample type, suggesting that they are more frequently mobilized between these microhabitats. Across intertidal zones, high intertidal, mid intertidal and low intertidal/swash communities were distinct and correlated with moisture, organic carbon and phosphate content, implying that the distribution of microorganisms is influenced by intertidal abiotic gradients. However, few taxa at the genus or species level individually contributed to this zonation pattern; rather, a unique combination of multiple microbial taxa was probably responsible. Although significant differences in microbial community composition on sandy beaches can be attributed to microhabitat and intertidal gradients, further investigations are needed to assess community assembly processes, the consistency of these distributions, and the functions of the majority of the microorganisms observed in the sand and their effects on the biogeochemistry and ecology of sandy beaches.