Investigating variability in microbial community composition in replicate environmental DNA samples down lake sediment cores

Lake sediments are natural archives that accumulate information on biological communities and their surrounding catchments. Paleolimnology has traditionally focussed on identifying fossilized organisms to reconstruct past environments. In the last decade, the application of molecular methodologies has increased in paleolimnological studies, but further research investigating factors such as sample heterogeneity and DNA degradation are required. In the present study we investigated bacterial community heterogeneity (16S rRNA metabarcoding) within depth slices (1-cm width). Sediment cores were collected from three lakes with differing sediment compositions. Samples were collected from a variety of depths which represent a period of time of approximately 1,200 years. Triplicate samples were collected from each depth slice and bacterial 16S rRNA metabarcoding was undertaken on each sample. Accumulation curves demonstrated that except for the deepest (oldest) slices, the combination of three replicate samples were insufficient to characterise the entire bacterial diversity. However, shared Amplicon Sequence Variants (ASVs) accounted for the majority of the reads in each depth slice (max. shared proportional read abundance 96%, 86%, 65% in the three lakes). Replicates within a depth slice generally clustered together in the Non-metric multidimensional scaling analysis. There was high community dissimilarity in older sediment in one of the cores, which was likely due to the laminae in the sediment core not being horizontal. Given that most paleolimnology studies explore broad scale shifts in community structure rather than seeking to identify rare species, this study demonstrates that a single sample is adequate to characterise shifts in dominant bacterial ASVs.