Unique Chemical Parameters and Microbial Activity Lead to Increased Archaeological Preservation within Vindolanda, UK.

Unique Chemical Parameters and Microbial Activity Lead to Increased Archaeological Preservation within Vindolanda, UK. | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Unique Chemical Parameters and Microbial Activity Lead to Increased Archaeological Preservation within Vindolanda, UK. C. H Orr, R Williams, H Halldórsdóttir, A Birley, E Greene, A Nelson, and 2 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-244364/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 04 Aug, 2021 Read the published version in Scientific Reports → Version 1 posted 8 You are reading this latest preprint version Abstract Waterlogged burial conditions impact upon artefact preservation. One major determinant of preservation is presence and behaviour of microorganisms, however, unravelling the mechanisms is challenging. In this study, we analysed elemental composition, bacterial diversity and community structure from excavation trenches at the Roman Site of Vindolanda, Northumberland, UK, using pXRF and 16s rRNA gene amplicon sequencing. Excavation trenches provide layers which can be classed as phosphorus sinks and which are strongly linked to the abundance of vivanite. The results indicated that microbial communities were dominant in Firmicutes, Bacteroidetes and Proteobacteria at a phylum level. Samples which also had visible vivianite presence showed that there were marked increases in Methylophilus . Methylophilus might be associated with favourable preservation in these anaerobic conditions. More research is needed to clearly link the presence of Methylophilus with vivianite production. The study emphasises the need for further integration of chemical and microbiome approaches, especially in good preservation areas, to explore microbial and chemical degradation mechanisms. General Microbiology Molecular Biology Vindolanda Anaerobic Amplicon Sequencing Vivianite Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Introduction Vindolanda, a Roman auxiliary fort situated south of Hadrian’s Wall near Bardon Mill, Northumberland UK, is a site known for its excellent preservation of leather and wooden artefacts 1 . This is typified by the writing tablets showing the correspondence of Julius Verecundus, Prefect of the First Cohort of Tungrians, dated to between 1st and 2nd centuries; an extensive assemblage of leather artefacts; and the discovery of the oldest Roman boxing gloves c.a. 120 AD. The fort is situated upon the Alston formation, which is calcium carbonate sedimentary bedrock, overlaid by loamy and clayey soils. However, this alone does not explain the exceptional level of artefact preservation 1 . Between Roman occupation periods, wooden and stone buildings were destroyed, sealed with thick layers of clays and then re-built upon, forming layers in which oxygen was excluded from the decomposing material underneath. This lack of oxygen led to the formation of anaerobic layers which are ideal preservation environments 2 and created areas of waterlogging on top of dense clay layers. Waterlogged environments have unique and complex hydrological chemistry, hence they are often too complex to determine individual indicators of good presevation 3 . A general lack of understanding of waterlogged environments has led English Heritage to adopt a reburial policy for in-situ preservation as is typified by the Rose and Globe Theatres in London 4 . Preservation indicators in waterlogged environments can include and be influenced by i) redox reactions, such as Fe 3+ reducing to Fe 2+ , and formations of chemical complexes like vivianite through combination of Fe 3+ and phosphate 2 , ii) microbiological species, such as decreasing abundance of Proteobacteria, Actinobacteria and Firmicutes , iii) hydrological indicators, such as acidity (pH values), conductivity and oxidation-reduction potential, iv) the decomposition material of surrounding highly organic matter, including those rich in phosphates 3 . Redox reactions are one of the most important and visually impactful processes present at Vindolanda. Chemically vivianite (Fe 3 (PO 4 ) 2 .8H 2 O) can form, when there is an availability of iron and phosphorus in the soil, lack of oxygen, low sulphur content and low pH values 2 , 5 , 6 . The phosphorus released from the decaying organic matter, such as floor material, animal faeces and plants is retained, through sorption to iron oxides and as such vivianite acts as a phosphorus sink. Upon excavation and exposure to oxygen, Fe 3+ is converted to Fe 2+ , which produces the unique blue colouration of vivianite 2 . Interestingly, vivianite is associated with other instances of good preservation including that of human remains: For example human DNA of increased quality can be extracted from samples located within the Brisbane burial grounds from underneath vivianite crusts 5 . The formation of vivianite is also thought to be microbially mediated with studies showing a range of microorganisms capable of producing the substance in anaerobic culture when exposed to the correct conditions 6 – 8 . Thus, despite vivianite acting as a visible preservation indicator, the exact chemical and microbiological formation, production and generation mechanisms are not completely understood. Microorganisms play a key role in decomposition, with above ground artefacts degraded by a range of bacteria, fungi and insects relatively quickly 9 . The speed of this degradation is impacted by oxygen availability, sample age and material which in turn impact upon the microbial community capable of carrying out decay 10 . Upon burial, particularly under anaerobic conditions, many fungal species (for example, soft rot fungi) typically involved in degradation cannot act. This means that the majority of degradation in anaerobic conditions is carried out by bacteria 11 . There is also evidence in marine sediments that key microbial communities, linked to anoxic environments can be central to increased preservation 12 . Ultimately, the mechanisms responsible for the excellent preservation observed at Vindolanda are unknown and will be explored within this study. Our assumption is that the exceptional levels of preservation seen at Vindolanda are brought about by a combination of environmental factors specifically linked to the unusual anaerobic and high iron conditions observed, through an interaction between local geochemistry, hydrology and past human activity at the site (clay sealing). Such conditions in turn impact on the microbial community within the soil, producing vivianite and limiting growth of degrading organisms. The initial aim of this study was to characterise the microbial community within different occupation layers of archaeological soil at Vindolanda. Bacterial communities within the different occupation layers which exhibit various degrees of preservation were compared with non-occupation environmental soil layers. The overall aim was then to correlate changes in environmental variables within the site and begin to unpick whether the microorganisms, environmental conditions or a combination of the two may be responsible for the good preservation observed at Vindolanda. As the factors influencing preservation and decomposition are poorly understood it is hoped that this study will help to shape current thinking of factors responsible for improved preservation at archaeological sites, to allow strengthening of in-situ preservation techniques. To our knowledge, this is the first reported study of a Roman burial environment using both chemistry and microbiological data, obtained using next generation sequencing techniques. Results Environmental variables pH values were consistent between the soil locations (depth1p (d1p) = 6.33 ± 0.06, depth2p (d2p) = 6.35 ± 0.20, depth3i (d3i) = 5.56 ± 0.88, depth4t (d4t) = 6.40 ± 0.16, depth5c (c1) = 6.07 ± 0.29, depth6c (c2) = 6.58 ± 0.47). Kruskal-Wallis tests showed no significant difference in pH value (χ2 = 46.38, df = 45, p = 0.42; supplementary date Fig. 1) between soil locations and significant differences in moisture content (χ2 = 53, df = 17, p < 0.001; supplementary date Fig. 1). Nemenyi tests showed that d1p and d2p had significantly more moisture than c1 and c2, but all other comparisons were not significantly different. Elemental analysis obtained from pXRF analysis showed significant changes in elemental composition between the differing locations. D1p and d2p, the sites of best preservation, displayed different elemental composition with respect to phosphorus, sulphur and iron (P < 0.001). Pairwise comparisons using Dwass-Steele-Critchlow-Fligner (DSCF) showed that the interaction between phosphorus, sulphur and iron in d1p and d2p was not significantly different to each other, but were significantly different to d3i, d4t, c1 and c2. Elemental composition for major and minor elements are provided in Fig. 1. Phosphorus has a range of Max = 1.094; Min = 0.512; Range = 0.582; Mean = 0.813; Median = 0.806 in d1p and d2p. Microbial community profile There is a shift in the microbial community profile between the preservation samples at d1p and d2p and interface layer d3i when compared with other locations (Figure 2). Spearman’s rank correlation analysis showed significant positive correlations with depth and Acidobacteria (ρ = 0.825, P = 2.52 x10 -5 ), Actinobacteria (ρ = 0.887, P = 9,189 x10 -7 ), Chloroflexi (ρ = 0.724, P = 6.78 x10 -4 ) and Planctomyces (ρ = 0.699, P = 1.24 x10 -3 ) and negative correlations with Firmicutes (ρ = -0.937, P = 9.76 x10 -9 ), Bacteroidetes (ρ = -0.687, P = 1.65 x10 -3 ) and Spirochaetes (ρ = -0.901, P = 3.347 x10 -7 ). There was no significant correlation linking pH values to relative abundance of any phyla. At OTU and genera level there were marked differences between the soil locations clearly separating samples from occupation layers with good preservation and interface samples from others (Figure 3). Specifically, identifiable dominating genera within preservation and interface layers were Sulfuricurvum , Flavobacterium and Methylophilus , with Ktedonobacter , Geobacter and Methylobacter dominating in non-preservation layers. Although the composition of the community changes this was not reflected in a shift in diversity with both Shannon diversity and Evenness measures showing little difference between the locations: d1p H’ 5.724 ± 0.14, Evenness 0.747 ± 0.001; d2p H’ 6.340 ± 0.12, Evenness 0.773 ± 0.006; d3i H’ 5.286 ± 0.10, Evenness 0.688 ± 0.011, d4t H’ 6.741 ± 0.03, Evenness 0.825 ± 0.000; c1 H’ 5.936 ± 0.09, Evenness 0.7778 ± 0.002; c2 H’ 6.120 ± 0.09, Evenness 0.779 ± 0.009. Microorganisms indicative of preservation layers LEfSe analysis of the microbial community revealed specific community structure associated with preservation layers characterised by Firmicutes, Epsilonproteobacteria, BRC1 and Chloroflexi with no phyla indicative of non-preservation layers due to varied communities. Tentative organisms which may be indicative of preservation layers are identified as: the Firmicute Turicibacter ; Alphaproteobacteria Filomicrobium and Sphingobium ; Betaproteobacteria Pusillimonas and Vogesella ; Epsilonproteobacteria Arcobacter; and Gammaproteobacteria Thiovirga (Figure 4) Composition of samples that contain vivianite. When compared with soils within 10 cm, soils with vivianite formation showed comparable Shannon diversity (visible vivianite = 3.211 ± 1.65; no visible vivianite = 3.636 ± 1.97) and evenness (visible vivianite = 0.598 ± 0.21; no visible vivianite = 0.618 ± 0.23) at phyla level. Adjacent soils profiles were not dissimilar to the preservation layers shown in Figure 2 (Figure 5). However, sections with vivianite were dominated by a single genus (50.3% of the total OTU within one of the samples) Methylophilus compared with sections without (P=0.0162) (Figure 5). Discussion Elemental composition within the soil layers Overall, there were significant shifts in elemental composition between soil samples at Vindolanda specifically with increases in phosphorus, sulphur and iron within d1p and d2p preservation locations. Non-preservation locations d3i, d4t, c1 and c2 showed consistent elemental composition regardless of the different trench locations. Whilst preservation of artefacts is poor in d3i and d4t there is evidence of Roman occupation. However, there is no such evidence in c1 and c2 suggesting they sit outside of the settlement. Interestingly, there was little observable difference within pH values between any of the samples. Iron and sulphate reduction, denitrification and nitrate respiration are associated with anaerobic environments 13 , with dissimilatory Fe 3+ reduction being the most abundant anaerobic metabolism 14 . Anaerobic soil conditions are typically characterised by the accumulation of iron cations (Fe 2+ and Fe 3+ ) as soil-available Fe is reduced. In aerobic conditions oxidation of Fe 2+ to Fe seems to correlate to a speeding up of microbial metabolism and a decomposition of organic matter both within field and pot experiments 15 . However, pXRF analysis used in the current study cannot differentiate between the different oxidative states of iron; it does however, suggest that there is an increase in the amount of iron present within these locations not just a change in redox state. Phosphorus release is associated with increased iron reduction within anaerobic soil environments with the dynamics of phosphorus release and Fe 3+ reduction closely linked 16 . However, there is also a relationship between the cycling of iron and phosphorus with the concentration of sulphate found within waterlogged environments. Within such conditions sulphate reduction results in increases in phosphorus release 17 . The iron and phosphorus concentrations are higher at Vindolanda within locations showing preservation and could, therefore suggest why increased vivianite are observed here 2 . Ratios of less than 1.5 sulphur:iron are needed for vivianite formation within marine sediments 18 ; all soils used within this study met these conditions. However, an excess of phosphorus is also needed which is only seen in d1p and d2p. Certain soil types are more likely to result in good archaeological preservation than others with anaerobic soils rich in organic material being more likely to show reduced biological oxidation of wood, leather and soft tissue 19 , partly due to the microorganism which flourish in such environments 10 . Reduced forms of phosphorus, sulphur and iron are more likely to be found within anaerobic soils compared with aerobic. Furthermore, the changes in moisture conditions at Vindolanda are linked to the periodic release of agricultural leachate which would be expected to increase soil degradation but may allow reaction with organic material within the soil, thus facilitating a natural ‘tanning’ process 19 . The Microbial community within the Vindolanda soil layers Preservation and interface locations d1p, d2p and d3i are distinct from other samples in terms of composition of phyla and genera. Shifts in bacterial phyla are correlated to location of sample but not pH values or elemental composition. Preservation locations are dominated by Firmicutes, Bacteroidetes and Proteobacteria at a phylum level with LEfSe analysis suggesting that Firmicutes, Epsilonproteobacteria, BRC1 and Chloroflexi may be indicative of preservation. Non-preservation and control locations were typified by increased abundance of Acidobacteria, Actinobacteria and Planctomycetes . Shannon diversity remained high within preservation layers and was comparable to the surrounding soil. Maintenance of diversity and evenness within good archaeological preserving environments has been shown previously in No. 1 Wangshanqiao Chu Tomb in Jingzhou, China 20 . Previous studies have shown that at the phylum level soils within the top approximately 20 cm are likely to be dominated by Acidobacteria, Actinobacteria and Alphaproteobacteria with Firmicutes and Spirochaetes being more abundant below 50cm 21 . Acidobacteria and Actinobacteria are increased within d4t, c1 and c2 compared with d1p, d2p and d3i. Actinobacteria are widely recognised as being crucial to aerobic degradation of organic matter particularly that which is cellulose based 22 , 23 . Acidobacteria are a diverse phylum which are mostly heterotrophic in metabolism and can degrade a wide range of carbohydrate sources 24 . Members of the Acidobacteria group 1–4 have been linked with degradation of cellulose based material 25 and are seen in significant abundance within layers 4–6 of this study. Interestingly they are greatly reduced within the preservation layers despite the pH values remaining consistent. This is probably due to the anaerobic conditions present within these layers as the majority of Acidobacteria are aerobic 24 . Firmicutes, Bacteroidetes and Proteobacteria were significantly increased within the preservation layers. However, this is most likely a result of the anaerobic conditions found rather than these phyla not being decomposers. Proteobacteria, Firmicutes and Bacteroidetes are the most dominant phyla in solid waste decomposition within anaerobic phases 26 with Firmicutes and Bacteroidetes proposed as the main microorganisms involved in cellulose degradation 25 . There are very few studies assessing microbiology within archaeological sites. However, those available correspond with our findings of consistent microbial diversity between layers which show preservation and those which do not 20 , 27 . Siles (2018) also found varied soil profiles that were common to soil environments and recorded no dominant ‘rare’ microbial taxa as clear indicators of preservation 27 . At the genus level preservation locations showed increased abundance of Sulficurvum , Flavobacterium and Methylophilus with additional genera Turicibacter , Filomicrobium , Sphingobium , Pusillimonas , Vogesella , Arcobacter and Thiovirga being identified as indicative of preservation by LEfSe. These microorganisms are not typically indicative of depth or oxygen availability as changes at phylum level were. Instead Sulficurvum , Flavobacterium and Methylophilus are all capable of degrading polycyclic aromatic hydrocarbons (PAHs) and are found to be encouraged by increased concentrations of iron and sulphur 28 , 29 . Flavobacterium and Methylophilus are aerobic, whilst Sulfuricurvum are facultative anaerobes which are major contributors to the sulphur cycle given their abilities to oxidise sulphur and are observed in the presence of excess methane 30 – 32 . This coupled with the changes in soil chemistry involved suggests an unusual soil environment within the preservation layers at Vindolanda. Microbial community with vivianite Interestingly, further investigation of samples with visible vivianite showed that these soils contained marked increases in Methylophilus compared with adjacent soils. Methylophilus are methylotrophic Betaproteobacteria commonly occurring in activated sludge, mud, and aquatic habitats, and are capable of utilizing methanol or methylamine as a sole source of carbon and energy 33 . Members of the genus are thought to be strictly aerobic and have been extracted from various soils 32 . Usually Methylophilus would be expected to be found in low abundance within soil environments. However, previous studies have found increased amounts of CO 2 within soil samples can lead to proliferation of Methylophilus presumably due to an increased abundance of methane 34 . Methylophilus is capable of solubilizing inorganic phosphorus 35 and produces siderophores and indole acetic acid 32 . Studies have shown that Methylophilus species are capable of reducing iron species such as ferrihydrite within conditions of increased methanol with evidence that the ferrihydrite was transformed into crystalline magnetite, a mixed Fe 3+ oxide 35 . Furthermore, other microorganisms capable of producing magnetite can also produce vivianite in phosphate rich environments 36 . Studies in Brazil have shown that, with the exception of environments that are high in sulphur, phosphorus and iron the production of iron species such as magnetite and haematite are relatively common while the production of vivianite is more restricted 37 . These conditions are similar to those we have identified within the preservation layers at Vindolanda. Additionally, analysis within the South China Sea showed that concentrations of iron-bound phosphorus and organic-phosphorus were influenced by sulfate-driven and metal-driven anaerobic oxidation of methane resulting in the release of phosphates and Fe 2+ and a subsequent significantly increased formation of vivianite 38 . Conclusion Our findings taken together suggest that the soil at Vindolanda contains distinctive chemical and microbial profiles associated with layers that show increased preservation of organic/inorganic artefacts. These shifts are not simply the result of changes in depth and oxygen availability but suggest that the specific conditions potentially allow the microbial production of vivianite which may lead to enhanced preservation of artefacts. The role of the occupation layers as phosphorus sinks is strongly linked to the abundance of vivanite. More research is needed to clearly link the presence of Methylophilus with vivianite production. Methods Sampling and archaeological context Soil was sampled from two trenches within the Vindolanda site (54.9911°N 2.3608°W). The external surface was scraped back by approximately 5cm. Samples were then taken in triplicate at four equal depths (d1p, d2p, d3i and d4t), each representing a separate archaeological context (see supplementary Fig. 2). The top 40cm of the trench had was exposed for approximately 3 months prior to the sampling taking place down to the interface level of d4t. The lower 1m section was excavated to the level of d1p four weeks prior to the sampling taking place. D1p and d2p correlate to the deepest layers and are associated with the best quality preservation. D3i and D4t are associated with lower quality preservation while d3i forming an interface (i) layer. Soil was sampled in triplicate (A-B) from a control (c) trench showing no evidence occupation or archaeological preservation, at two depths (c1-c2) (see supplementary Fig. 3). In addition, further samples were taken from a rampart section, as detailed in Taylor et al 2019. The rampart sample was sampled in triplicate from the vivianite visible layer (blue) as well as triplicate samples within 10 cm which showed no visible vivianite formation. Moisture Triplicate soils were weighed before and after drying in an oven overnight at 105°C to determine the moisture content of the original sample using Eq. 1. Equation 1: 𝑤𝑒𝑡 𝑤𝑒𝑖𝑔ℎ𝑡 (𝑔) − 𝑑𝑟𝑦 𝑤𝑒𝑖𝑔ℎ𝑡 (𝑔) / 𝑑𝑟𝑦 𝑤𝑒𝑖𝑔ℎ𝑡 (𝑔) × 100 = 𝑚𝑜𝑖𝑠𝑡𝑢𝑟𝑒 𝑐𝑜𝑛𝑡𝑒𝑛𝑡 % pH Dried soils (10g) were placed into plastic cups with distilled water added (20ml) for a 1:2 proportion. The solution was stirred with a glass rod and settled for 30 minutes. pH values was analysed using an Advanced Jenway 3510 pH meter calibrated with commercial buffers of pH 4, 7 and 10. The electrode probe was rinsed between each sample. pXRF Soils were analysed following the method provided by Williams (2020) 39 . The dried soil samples were homogenised with mortar and pestle for 140 seconds and sieved to 2mm. Soil were loaded into XRF sample cups (SPEX CertiPrep™ 3529) and covered with 5 µm polypropylene thin-film (SPEX™ SamplePrep 3520 window film). A Thermo Niton™ XL3t GOLDD + pXRF with an Ag anode (6-50kV, 0-200 µA max X-ray tube) was used to analyse the soil. The pXRF was warmed up, system checked against the internal 1¼ Cr–½ Mo coupon, and scanned against a blank standard and NIST 2709a standard (y = 0.9674x -0.0089, r 2 = 0.9998). The NIST 2709a is a San Joaquin soil, used because it was certified for most elements of interest within the single calibration sample. The pXRF was periodically reset and system checked to account for drift. Soil samples were analysed using Mining mode (fundamental parameters), with 30-second scans for the main filter (50 kV, ≤ 50 µA), low filter (20 kV, ≤ 100 µA) and high filter (50 kV, ≤ 40 µA), and a 60-second scan for the light filter (6 kV, ≤ 200 µA). DNA extraction DNA was extracted from 0.25 g of each soil sample using the Qiagen DNeasy PowerSoil Kit following manufacturer’s instructions. DNA quality and yield were measured using NanoDrop 1000 spectrophotometer (V3.8.1, ThermoFisher) and agarose gel electrophoresis. DNA was stored at -80°C prior to further analysis. Next Generation Sequencing Bacterial community DNA sequencing was made (NU-OMICs, Northumbria University, Newcastle Upon Tyne, U.K.) by targeting the V4 region of the 16S rRNA gene according to Kozich (2013) 40 . The raw sequencing reads were processed in FASTQ format and analysed with Mothur software package (version 1.36.1) (University of Michigan, U.S.A.). UCHIME was used to quality check and filter the FASTA formatted sequences for chimeras. These were aligned to the SILVA reference, and taxonomic identification of the reads was assessed by assigning sequences to operational taxonomic units (OTUs) using Ribosomal Database Project (RDP) classification. PCR negative controls were run and sequenced in parallel to the experimental samples. The OTUs recorded in negative controls and samples were excluded from further analysis. Non-bacterial sequences (Archaea/Eukaryota/Unclassified) were discarded. OTUs less than 3% were classified as rare taxa with these and the unclassified OTUs omitted from the plots. Statistical analysis Data were analysed in R 4.0.2 41 using the car, DescTools and PMCMR packages 42–44 . Normality and variance failed the Q-Q plot, Shapiro-Wilk’s and Levene’s assumption tests. Data were analysed using Kruskal-Wallis rank sum tests, followed by Nemenyi tests with Chi-Square distribution to identify significant differences between burial layers. Tukey box plots of the median and quartiles were produced using the ggplot2 package 45 to show the change in raw elemental concentration with increasing burial depth. The tidypaleo package in tidyverse 46 was used to produce stratigraphic plots of the moisture content and pH values against burial depth. Operational taxonomic unit (OUT) tables and metadata were analysed using R version 4.0.2. Samples were faceted to rarefy, relative abundance data was compiled and most abundant phyla and genera were sought. Within the readr package 47 significant differences were identified using ANOVA and TukeyHSD tests and Spearman Correlation was used to determine significant correlation between variables. Figures were created using gplots 48 . Taxa indicative of preservation were identified, and histograms generated following Linear Discriminant Analysis Effect Size (LEfSe) analysis completed in Galaxy 49 . Declarations Acknowledgements The authors thank Dr Andrew Birley and the Vindolanda Trust for providing access to the archaeological samples. There are no conflicts of interest to declare. The data that support the findings of this study are available in the additional supporting information. Author contributions GT and CO wrote the first draft of the manuscript and participated in sample collection and data analysis. RW and HH prepared and analysed the pXRF samples. AN organized and conducted the Next Generation Sequencing. AB and EG was involved in the fieldwork and sample collection. TKR sample preparation for Next generation sequencing. References Birley, R. Vindolanda: A Roman frontier post on Hadrian's Wall (Thames and Hudson, London, 1977). Taylor, G., Hefford, R. J. W., Birley, A. & Huntley, J. P. Identifying the “blue substance” at the Roman site of Vindolanda, Northumberland. Journal of Archaeological Science: Reports. 24 , 582–587 (2019). Holden, J. et al. Hydrological controls of in situ preservation of waterlogged archaeological deposits. Earth-Science Reviews. 78 (1–2), 59–83 (2006). Caple, C. Reburial of waterlogged wood, the problems and potential of this conservation technique. International Biodeterioration and Biodegradation. 34 (1), 61–72 (1994). McGowan, G. & Prangnell, J. The significance of vivianite in archaeological settings. Geoarchaeology. 21 (1), 93–111 (2006). Sánchez-Román, M., Puente-Sánchez, F., Parro, V. & Amils, R. Nucleation of Fe-rich phosphates and carbonates on microbial cells and exopolymeric substances. Frontiers in Microbiology. 2015;6(SEP). Miot, J. et al. Transformation of vivianite by anaerobic nitrate-reducing iron-oxidizing bacteria. Geobiology. 7 (3), 373–384 (2009). Volkland, H. P. et al. Bacterial phosphating of mild (unalloyed) steel. Applied and Environmental Microbiology. 66 (10), 4389–4395 (2000). Blanchette, R. A., Nilsson, T., Daniel, G. & Abad, A. Biological degradation of wood. American Chemical Society. 1990:141–174. Björdal, C. G., Nilsson, T. & Daniel, G. Microbial decay of waterlogged archaeological wood found in Sweden. Applicable to archaeology and conservation. International Biodeterioration and Biodegradation. 43 (1–2), 63–73 (1999). Björdal, C. G. Microbial degradation of waterlogged archaeological wood. Journal of Cultural Heritage. 13 (3 SUPPL.), S118–S122 (2012). Jessen, G. L. et al. Hypoxia causes preservation of labile organic matter and changes seafloor microbial community composition (Black Sea). Sci Adv. 2017;3(2). Tiedje, J. M., Sexstone, A. J., Parkin, T. B. & Revsbech, N. P. Anaerobic processes in soil. Plant and Soil. 76 (1–3), 197–212 (1984). Keiluweit, M., Wanzek, T., Kleber, M., Nico, P. & Fendorf, S. Anaerobic microsites have an unaccounted role in soil carbon stabilization. Nat Commun. 2017;8(1). Hall, S. J. & Silver, W. L. Iron oxidation stimulates organic matter decomposition in humid tropical forest soils. Global Change Biol. 19 (9), 2804–2813 (2013). Forsmann, D. M. & Kjaergaard, C. Phosphorus release from anaerobic peat soils during convective discharge - Effect of soil Fe: P molar ratio and preferential flow. Geoderma. 223–225 (1), 21–32 (2014). Lehtoranta, J., Ekholm, P., Wahlström, S., Tallberg, P. & Uusitalo, R. Labile organic carbon regulates phosphorus release from eroded soil transported into anaerobic coastal systems. Ambio. 44 (2), 263–273 (2015). Rothe, M., Kleeberg, A. & Hupfer, M. The occurrence, identification and environmental relevance of vivianite in waterlogged soils and aquatic sediments. Earth-Science Reviews. 158 , 51–64 (2016). Kibblewhite, M., Tóth, G. & Hermann, T. Predicting the preservation of cultural artefacts and buried materials in soil. Science of the Total Environment. 529 , 249–263 (2015). Xu, J., Wei, Y., Jia, H., Xiao, L. & Gong, D. A new perspective on studying burial environment before archaeological excavation: Analyzing bacterial community distribution by high-throughput sequencing. Scientific Reports. 2017;7. Douterelo, I., Goulder, R. & Lillie, M. Soil microbial community response to land-management and depth, related to the degradation of organic matter in English wetlands: Implications for the in situ preservation of archaeological remains. Appl Soil Ecol. 44 (3), 219–227 (2010). Wang, C. et al. Metagenomic analysis of microbial consortia enriched from compost: New insights into the role of Actinobacteria in lignocellulose decomposition. Biotechnol Biofuels. 2016;9(1). Pankratov, T. A., Dedysh, S. N. & Zavarzin, G. A. The leading role of actinobacteria in aerobic cellulose degradation in Sphagnum peat bogs. Dokl Biol Sci. 410 (1), 428–430 (2006). Kielak, A. M., Barreto, C. C., Kowalchuk, G. A., van Veen, J. A. & Kuramae, E. E. The ecology of Acidobacteria: Moving beyond genes and genomes. Frontiers in Microbiology. 2016;7(MAY). Eichorst, S. A., Kuske, C. R. & Schmidt, T. M. Influence of plant polymers on the distribution and cultivation of bacteria in the phylum acidobacteria. Applied and Environmental Microbiology. 77 (2), 586–596 (2011). Yang, S. & Song, L. Succession of bacterial community structure and metabolic function during solid waste decomposition. Bioresource Technology. 2019;291. Siles, J. A., Öhlinger, B., Cajthaml, T., Kistler, E. & Margesin, R. Characterization of soil bacterial, archaeal and fungal communities inhabiting archaeological human-impacted layers at Monte Iato settlement (Sicily, Italy). Scientific Reports. 2018;8(1). Sun, W. et al. Profiling microbial community in a watershed heavily contaminated by an active antimony (Sb) mine in Southwest China. Science of the Total Environment. 550 , 297–308 (2016). Haaijer, S. C. M. et al. Bacteria associated with iron seeps in a sulfur-rich, neutral pH, freshwater ecosystem. ISME Journal. 2 (12), 1231–1242 (2008). Li, X. et al. Metagenome-assembled genome sequence of sulfuricurvum sp. strain IAE1, isolated from a 4-chlorophenol-degrading consortium. Micro Res Ann. 2019;8(31). Hubert, C. R. J. et al. Massive dominance of Epsilonproteobacteria in formation waters from a Canadian oil sands reservoir containing severely biodegraded oil. Environmental Microbiology. 14 (2), 387–404 (2012). Madhaiyan, M., Poonguzhali, S., Kwon, S. W. & Sa, T. M. Methylophilus rhizosphaerae sp. nov., a restricted facultative methylotroph isolated from rice rhizosphere soil. Int J Syst Evol Microbiol. 59 (11), 2904–2908 (2009). Doronina, N., Kaparullina, E. & Trotsenko, Y.. The Family Methylophilaceae. In: Rosenberg E, DeLong EF, Lory S, Stackebrandt E, Thompson F, eds. The Prokaryotes: Alphaproteobacteria and Betaproteobacteria. Berlin, Heidelberg:Springer Berlin Heidelberg; 2014:869–880. Ma, J. et al. Short-term effects of CO2 leakage on the soil bacterial community in a simulated gas leakage scenario. PeerJ. 2017 , 11 (2017). Yang, Y. et al. Extracellular electron transfer of Methylophilus methylotrophs. Process Biochem. 94 , 313–318 (2020). Veeramani, H. et al. Products of abiotic U(VI) reduction by biogenic magnetite and vivianite. Geochimica et Cosmochimica Acta. 75 (9), 2512–2528 (2011). Volpato, J. Relation between geochemistry and minerals (magnetite, hematite, pyrite and vivianite) in contamineted estuarine sediments, Santos, Sao Paulo, Brazil., Universidade Estadual de Campinas; 2015. Wu, D. et al. Zone of metal-driven anaerobic oxidation of methane is an important sink for phosphorus in the Taixinan Basin, South China Sea. Marine Geology. 427 , 1–10 (2020). Williams, R., Taylor, G. & Orr, C. pXRF method development for elemental analysis of archaeological soil. Archaeometry. 62 (6), 1145–1163 (2020). Kozich, J., Westcott, S. L., Baxter, N. T., Highlander, S. K. & Schloss, P. D. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the miseq illumina sequencing platform. Applied and Environmental Microbiology. 79 (17), 5112–5120 (2013). R: A language and environment for statistical computing [computer program]. Vienna:R Foundation for Statistical Computing; 2020. Signorell, A. DescTools: Tools for descriptive statistics. 2016. Pohlert, T. The pairwise multiple comparison of mean ranks package (PMCMR). 2015. Fox, J. & Weisberg, S. An {R} companion to applied regression. 2011. Wickham, H. ggplot2: Elegant graphics for data analysis. 2016. Wickham, H. et al. Welcome to the Tidyverse. The journal of open source software. 4 (43), 1686 (2019). Wickham, H., Hester, J. & Francois, R. readr: read rectangular Text Data. R package version 1.3.1.. 2018; https://CRAN.R-project.org/package=readr . Warnes, G. et al. L, Gplots:Various R Programming Tools for Plotting Data. R package 3.03.. 2020; https://CRAN.R-project.org/package=gplots . Segata, N. et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12(6). Additional Declarations No competing interests reported. Supplementary Files TaylorSupplentaryFigureFINAL.docx Cite Share Download PDF Status: Published Journal Publication published 04 Aug, 2021 Read the published version in Scientific Reports → Version 1 posted Editorial decision: Major revision 08 Apr, 2021 Reviews received at journal 29 Mar, 2021 Reviewers agreed at journal 19 Mar, 2021 Reviewers invited by journal 08 Mar, 2021 Editor assigned by journal 03 Mar, 2021 Editor invited by journal 03 Mar, 2021 Submission checks completed at journal 03 Mar, 2021 First submitted to journal 15 Feb, 2021 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-244364","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":14656036,"identity":"212c772f-6aa8-4bc1-95b6-1d8fa7a5284a","order_by":0,"name":"C. H Orr","email":"","orcid":"","institution":"Teesside University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"C.","middleName":"H","lastName":"Orr","suffix":""},{"id":14656037,"identity":"81ff5dc4-32fa-4948-abbe-4b2d7ba51920","order_by":1,"name":"R Williams","email":"","orcid":"","institution":"Teesside University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"R","middleName":"","lastName":"Williams","suffix":""},{"id":14656038,"identity":"8cdc0ac7-798d-471c-b6e3-b64993f21b6f","order_by":2,"name":"H Halldórsdóttir","email":"","orcid":"","institution":"Teesside University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"H","middleName":"","lastName":"Halldórsdóttir","suffix":""},{"id":14656039,"identity":"9f49b635-7ea2-4425-a0ca-4cac3fd7da21","order_by":3,"name":"A Birley","email":"","orcid":"","institution":"Vindolanda Trust","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"A","middleName":"","lastName":"Birley","suffix":""},{"id":14656040,"identity":"9bce9408-2487-4dd5-ad93-6d19e1124280","order_by":4,"name":"E Greene","email":"","orcid":"","institution":"University of Western Ontario","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"E","middleName":"","lastName":"Greene","suffix":""},{"id":14656041,"identity":"d4451133-5a7d-4bb3-b77e-3aa078e7440f","order_by":5,"name":"A Nelson","email":"","orcid":"","institution":"Northumbria University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"A","middleName":"","lastName":"Nelson","suffix":""},{"id":14656042,"identity":"026545a4-77e7-4547-b0c8-98450d147885","order_by":6,"name":"T.K Ralebitso-Senior","email":"","orcid":"","institution":"Liverpool John Moores University","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"T.K","middleName":"","lastName":"Ralebitso-Senior","suffix":""},{"id":14656043,"identity":"0415c38a-fb0a-4ce7-ad6e-7a61fa106e96","order_by":7,"name":"G Taylor","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA2ElEQVRIie2RMQrCQBBFvy6szcK2EYu9woRAQBDxKAFhLbUULIwE4iU8h3VCIF7BTkSwshELIzauNnZZ7QT3wRQD//EHBnA4fpYJIN8bt8WZGQLa6dcKlZ8qaiDzi6DeONgl+b5CX8HTUa3iZ4x1BOnuuuTDQGDox57O6pWYwSgFhaUIO+bMCN4otins/lSCVF5vFeZ2RYHxVwtxwSFQGMVyGDEWdlekySt1YNyNn4pjVN+yXBy2p2mPZFLsz9V0pmRLU31LATTFe/3gkSoGGpUt5XA4HP/NA+NKNrV+Fm2MAAAAAElFTkSuQmCC","orcid":"","institution":"Teesside University","correspondingAuthor":true,"submittingAuthor":false,"prefix":"","firstName":"G","middleName":"","lastName":"Taylor","suffix":""}],"badges":[],"createdAt":"2021-02-15 22:29:02","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-244364/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-244364/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1038/s41598-021-94853-7","type":"published","date":"2021-08-04T15:04:57+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":6799457,"identity":"7407a1fe-e4c6-4654-9586-96a46b7c5b87","added_by":"auto","created_at":"2021-03-10 15:01:16","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":134088,"visible":true,"origin":"","legend":"Elemental analysis for minor and major elements from sample d1p, d2p, d3i, d4t, c1 and c2.","description":"","filename":"floatimage1.png","url":"https://assets-eu.researchsquare.com/files/rs-244364/v1/4f3fb096e52184ea17e244d7.png"},{"id":6799020,"identity":"5304e6ec-376f-4e54-829e-5d2f628183dc","added_by":"auto","created_at":"2021-03-10 14:58:09","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":187670,"visible":true,"origin":"","legend":"Relative abundance of bacterial phyla across six locations of the study. Replicates are averaged. The ten most abundant phyla are shown with the remainder shown as ‘other’ ","description":"","filename":"2.png","url":"https://assets-eu.researchsquare.com/files/rs-244364/v1/1659e774b64e76cf14389106.png"},{"id":6799453,"identity":"a8e17d0d-baa9-41cd-a9ca-a32fd6c7ca3d","added_by":"auto","created_at":"2021-03-10 15:01:15","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":426503,"visible":true,"origin":"","legend":"Heatmap showing the distribution of OTU within each replicate sample. Replicate samples are shown with –a, -b and –c. Data is simplified so that OTU are only depicted if they contribute more than 1% of a single sample. Dark red = 1, pale yellow = 0. Within OTU names ‘un’ = unclassified. ","description":"","filename":"3.png","url":"https://assets-eu.researchsquare.com/files/rs-244364/v1/62cebdca29328622e793e95a.png"},{"id":6799454,"identity":"2234cca3-f5d7-4c4f-b584-d473ce3a328b","added_by":"auto","created_at":"2021-03-10 15:01:15","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":307647,"visible":true,"origin":"","legend":"Histograms generated following LEfSe analysis showing relative abundance of genera identified as potentially indicative of preservation locations. Class shows if samples are within a preservation location. P values are listed in corner of each figure. NS = non-significant. Tukey’s test within each graph is represented by letters. Layers which share letters are not significantly different from each other. ","description":"","filename":"4.png","url":"https://assets-eu.researchsquare.com/files/rs-244364/v1/b438bcdf22545acd0b925385.png"},{"id":6799024,"identity":"c3c9523c-8e01-459b-9002-0b2202fef918","added_by":"auto","created_at":"2021-03-10 14:58:09","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":175883,"visible":true,"origin":"","legend":"Relative abundance of bacterial phyla (A) and genera (B) within samples showing visible vivianite (Viv) or no visible vivianite (No Viv). Replicates are averaged. The ten most abundant phyla/genera are shown with the remainder shown as ‘other’","description":"","filename":"5.png","url":"https://assets-eu.researchsquare.com/files/rs-244364/v1/e40fb1d8b565d7e6916784a2.png"},{"id":13676888,"identity":"23eea342-451a-4f0d-b2ae-a3934cf50afc","added_by":"auto","created_at":"2021-09-17 11:32:45","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":1370537,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-244364/v1/41d7b3dc-b268-4745-88c1-0519eb83bd19.pdf"},{"id":6799452,"identity":"82a48256-c541-48c4-b930-9a03786d03db","added_by":"auto","created_at":"2021-03-10 15:01:15","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":7781488,"visible":true,"origin":"","legend":"","description":"","filename":"TaylorSupplentaryFigureFINAL.docx","url":"https://assets-eu.researchsquare.com/files/rs-244364/v1/298ca424324fb4df42a57c3b.docx"}],"financialInterests":"No competing interests reported.","formattedTitle":"Unique Chemical Parameters and Microbial Activity Lead to Increased Archaeological Preservation within Vindolanda, UK.","fulltext":[{"header":"Introduction","content":" \u003cp\u003eVindolanda, a Roman auxiliary fort situated south of Hadrian\u0026rsquo;s Wall near Bardon Mill, Northumberland UK, is a site known for its excellent preservation of leather and wooden artefacts\u003csup\u003e\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e\u003c/sup\u003e. This is typified by the writing tablets showing the correspondence of Julius Verecundus, Prefect of the First Cohort of Tungrians, dated to between 1st and 2nd centuries; an extensive assemblage of leather artefacts; and the discovery of the oldest Roman boxing gloves c.a. 120 AD. The fort is situated upon the Alston formation, which is calcium carbonate sedimentary bedrock, overlaid by loamy and clayey soils. However, this alone does not explain the exceptional level of artefact preservation\u003csup\u003e\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e\u003c/sup\u003e. Between Roman occupation periods, wooden and stone buildings were destroyed, sealed with thick layers of clays and then re-built upon, forming layers in which oxygen was excluded from the decomposing material underneath. This lack of oxygen led to the formation of anaerobic layers which are ideal preservation environments\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e and created areas of waterlogging on top of dense clay layers.\u003c/p\u003e \u003cp\u003eWaterlogged environments have unique and complex hydrological chemistry, hence they are often too complex to determine individual indicators of good presevation\u003csup\u003e\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e\u003c/sup\u003e. A general lack of understanding of waterlogged environments has led English Heritage to adopt a reburial policy for \u003cem\u003ein-situ\u003c/em\u003e preservation as is typified by the Rose and Globe Theatres in London\u003csup\u003e\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e\u003c/sup\u003e. Preservation indicators in waterlogged environments can include and be influenced by i) redox reactions, such as Fe\u003csup\u003e3+\u003c/sup\u003e reducing to Fe\u003csup\u003e2+\u003c/sup\u003e, and formations of chemical complexes like vivianite through combination of Fe\u003csup\u003e3+\u003c/sup\u003e and phosphate\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e, ii) microbiological species, such as decreasing abundance of \u003cem\u003eProteobacteria, Actinobacteria\u003c/em\u003e and \u003cem\u003eFirmicutes\u003c/em\u003e, iii) hydrological indicators, such as acidity (pH values), conductivity and oxidation-reduction potential, iv) the decomposition material of surrounding highly organic matter, including those rich in phosphates\u003csup\u003e\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eRedox reactions are one of the most important and visually impactful processes present at Vindolanda. Chemically vivianite (Fe\u003csub\u003e3\u003c/sub\u003e(PO\u003csub\u003e4\u003c/sub\u003e)\u003csub\u003e2\u003c/sub\u003e.8H\u003csub\u003e2\u003c/sub\u003eO) can form, when there is an availability of iron and phosphorus in the soil, lack of oxygen, low sulphur content and low pH values\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e,\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e,\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e\u003c/sup\u003e. The phosphorus released from the decaying organic matter, such as floor material, animal faeces and plants is retained, through sorption to iron oxides and as such vivianite acts as a phosphorus sink. Upon excavation and exposure to oxygen, Fe\u003csup\u003e3+\u003c/sup\u003e is converted to Fe\u003csup\u003e2+\u003c/sup\u003e, which produces the unique blue colouration of vivianite\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e. Interestingly, vivianite is associated with other instances of good preservation including that of human remains: For example human DNA of increased quality can be extracted from samples located within the Brisbane burial grounds from underneath vivianite crusts\u003csup\u003e\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e\u003c/sup\u003e. The formation of vivianite is also thought to be microbially mediated with studies showing a range of microorganisms capable of producing the substance in anaerobic culture when exposed to the correct conditions\u003csup\u003e\u003cspan additionalcitationids=\"CR7\" citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e\u003c/sup\u003e. Thus, despite vivianite acting as a visible preservation indicator, the exact chemical and microbiological formation, production and generation mechanisms are not completely understood.\u003c/p\u003e \u003cp\u003eMicroorganisms play a key role in decomposition, with above ground artefacts degraded by a range of bacteria, fungi and insects relatively quickly\u003csup\u003e\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e\u003c/sup\u003e. The speed of this degradation is impacted by oxygen availability, sample age and material which in turn impact upon the microbial community capable of carrying out decay\u003csup\u003e\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e\u003c/sup\u003e. Upon burial, particularly under anaerobic conditions, many fungal species (for example, soft rot fungi) typically involved in degradation cannot act. This means that the majority of degradation in anaerobic conditions is carried out by bacteria\u003csup\u003e\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e\u003c/sup\u003e. There is also evidence in marine sediments that key microbial communities, linked to anoxic environments can be central to increased preservation\u003csup\u003e\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e\u003c/sup\u003e. Ultimately, the mechanisms responsible for the excellent preservation observed at Vindolanda are unknown and will be explored within this study. Our assumption is that the exceptional levels of preservation seen at Vindolanda are brought about by a combination of environmental factors specifically linked to the unusual anaerobic and high iron conditions observed, through an interaction between local geochemistry, hydrology and past human activity at the site (clay sealing). Such conditions in turn impact on the microbial community within the soil, producing vivianite and limiting growth of degrading organisms.\u003c/p\u003e \u003cp\u003eThe initial aim of this study was to characterise the microbial community within different occupation layers of archaeological soil at Vindolanda. Bacterial communities within the different occupation layers which exhibit various degrees of preservation were compared with non-occupation environmental soil layers. The overall aim was then to correlate changes in environmental variables within the site and begin to unpick whether the microorganisms, environmental conditions or a combination of the two may be responsible for the good preservation observed at Vindolanda. As the factors influencing preservation and decomposition are poorly understood it is hoped that this study will help to shape current thinking of factors responsible for improved preservation at archaeological sites, to allow strengthening of \u003cem\u003ein-situ\u003c/em\u003e preservation techniques. To our knowledge, this is the first reported study of a Roman burial environment using both chemistry and microbiological data, obtained using next generation sequencing techniques.\u003c/p\u003e "},{"header":"Results","content":"\u003cdiv\u003e\n\u003cp\u003e\u003cstrong\u003eEnvironmental variables\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003epH values were consistent between the soil locations (depth1p (d1p)\u0026thinsp;=\u0026thinsp;6.33\u0026thinsp;\u0026plusmn;\u0026thinsp;0.06, depth2p (d2p)\u0026thinsp;=\u0026thinsp;6.35\u0026thinsp;\u0026plusmn;\u0026thinsp;0.20, depth3i (d3i)\u0026thinsp;=\u0026thinsp;5.56\u0026thinsp;\u0026plusmn;\u0026thinsp;0.88, depth4t (d4t)\u0026thinsp;=\u0026thinsp;6.40\u0026thinsp;\u0026plusmn;\u0026thinsp;0.16, depth5c (c1)\u0026thinsp;=\u0026thinsp;6.07\u0026thinsp;\u0026plusmn;\u0026thinsp;0.29, depth6c (c2)\u0026thinsp;=\u0026thinsp;6.58\u0026thinsp;\u0026plusmn;\u0026thinsp;0.47). Kruskal-Wallis tests showed no significant difference in pH value (\u0026chi;2\u0026thinsp;=\u0026thinsp;46.38, df\u0026thinsp;=\u0026thinsp;45, p\u0026thinsp;=\u0026thinsp;0.42; supplementary date Fig.\u0026nbsp;1) between soil locations and significant differences in moisture content (\u0026chi;2\u0026thinsp;=\u0026thinsp;53, df\u0026thinsp;=\u0026thinsp;17, p\u0026thinsp;\u0026lt;\u0026thinsp;0.001; supplementary date Fig.\u0026nbsp;1). Nemenyi tests showed that d1p and d2p had significantly more moisture than c1 and c2, but all other comparisons were not significantly different.\u003c/p\u003e\n\u003cp\u003eElemental analysis obtained from pXRF analysis showed significant changes in elemental composition between the differing locations. D1p and d2p, the sites of best preservation, displayed different elemental composition with respect to phosphorus, sulphur and iron (P\u0026thinsp;\u0026lt;\u0026thinsp;0.001). Pairwise comparisons using Dwass-Steele-Critchlow-Fligner (DSCF) showed that the interaction between phosphorus, sulphur and iron in d1p and d2p was not significantly different to each other, but were significantly different to d3i, d4t, c1 and c2. Elemental composition for major and minor elements are provided in Fig.\u0026nbsp;1. Phosphorus has a range of Max\u0026thinsp;=\u0026thinsp;1.094; Min\u0026thinsp;=\u0026thinsp;0.512; Range\u0026thinsp;=\u0026thinsp;0.582; Mean\u0026thinsp;=\u0026thinsp;0.813; Median\u0026thinsp;=\u0026thinsp;0.806 in d1p and d2p.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eMicrobial community profile\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThere is a shift in the microbial community profile between the preservation samples at d1p and d2p and interface layer d3i when compared with other locations (Figure 2). Spearman\u0026rsquo;s rank correlation analysis showed significant positive correlations with depth and\u003cem\u003e Acidobacteria\u003c/em\u003e (\u0026rho; = 0.825, P = 2.52 x10\u003csup\u003e-5\u003c/sup\u003e), \u003cem\u003eActinobacteria\u003c/em\u003e (\u0026rho; = 0.887, P = 9,189 x10\u003csup\u003e-7\u003c/sup\u003e), \u003cem\u003eChloroflexi\u003c/em\u003e (\u0026rho; = 0.724, P = 6.78 x10\u003csup\u003e-4\u003c/sup\u003e) and \u003cem\u003ePlanctomyces \u003c/em\u003e(\u0026rho; = 0.699, P = 1.24 x10\u003csup\u003e-3\u003c/sup\u003e) and negative correlations with \u003cem\u003eFirmicutes\u003c/em\u003e (\u0026rho; = -0.937, P = 9.76 x10\u003csup\u003e-9\u003c/sup\u003e),\u003cem\u003e Bacteroidetes\u003c/em\u003e (\u0026rho; = -0.687, P = 1.65 x10\u003csup\u003e-3\u003c/sup\u003e) and \u003cem\u003eSpirochaetes\u003c/em\u003e (\u0026rho; = -0.901, P = 3.347 x10\u003csup\u003e-7\u003c/sup\u003e). There was no significant correlation linking pH values to relative abundance of any phyla.\u003c/p\u003e\n\u003cp\u003eAt OTU and genera level there were marked differences between the soil locations clearly separating samples from occupation layers with good preservation and interface samples from others (Figure 3). Specifically, identifiable dominating genera within preservation and interface layers were \u003cem\u003eSulfuricurvum\u003c/em\u003e, \u003cem\u003eFlavobacterium\u003c/em\u003e and \u003cem\u003eMethylophilus\u003c/em\u003e, with \u003cem\u003eKtedonobacter\u003c/em\u003e, \u003cem\u003eGeobacter\u003c/em\u003e and \u003cem\u003eMethylobacter\u003c/em\u003e dominating in non-preservation layers. Although the composition of the community changes this was not reflected in a shift in diversity with both Shannon diversity and Evenness measures showing little difference between the locations: d1p H\u0026rsquo; 5.724 \u0026plusmn; 0.14, Evenness 0.747 \u0026plusmn; 0.001; d2p H\u0026rsquo; 6.340 \u0026plusmn; 0.12, Evenness 0.773 \u0026plusmn; 0.006; d3i H\u0026rsquo; 5.286 \u0026plusmn; 0.10, Evenness 0.688 \u0026plusmn; 0.011, d4t H\u0026rsquo; 6.741 \u0026plusmn; 0.03, Evenness 0.825 \u0026plusmn; 0.000; c1 H\u0026rsquo; 5.936 \u0026plusmn; 0.09, Evenness 0.7778 \u0026plusmn; 0.002; c2 H\u0026rsquo; 6.120 \u0026plusmn; 0.09, Evenness 0.779 \u0026plusmn; 0.009.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eMicroorganisms indicative of preservation layers\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eLEfSe analysis of the microbial community revealed specific community structure associated with preservation layers characterised by \u003cem\u003eFirmicutes, Epsilonproteobacteria, BRC1\u003c/em\u003e and \u003cem\u003eChloroflexi\u003c/em\u003e with no phyla indicative of non-preservation layers due to varied communities. Tentative organisms which may be indicative of preservation layers are identified as: the Firmicute \u003cem\u003eTuricibacter\u003c/em\u003e; \u003cem\u003eAlphaproteobacteria\u003c/em\u003e\u003cem\u003eFilomicrobium\u003c/em\u003e and \u003cem\u003eSphingobium\u003c/em\u003e; Betaproteobacteria \u003cem\u003ePusillimonas\u003c/em\u003e and \u003cem\u003eVogesella\u003c/em\u003e; Epsilonproteobacteria \u003cem\u003eArcobacter;\u003c/em\u003e and Gammaproteobacteria \u003cem\u003eThiovirga\u003c/em\u003e (Figure 4)\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eComposition of samples that contain vivianite.\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWhen compared with soils within 10 cm, soils with vivianite formation showed comparable Shannon diversity (visible vivianite = 3.211 \u0026plusmn; 1.65; no visible vivianite = 3.636 \u0026plusmn; 1.97) and evenness (visible vivianite = 0.598 \u0026plusmn; 0.21; no visible vivianite = 0.618 \u0026plusmn; 0.23) at phyla level. Adjacent soils profiles were not dissimilar to the preservation layers shown in Figure 2 (Figure 5). However, sections with vivianite were dominated by a single genus (50.3% of the total OTU within one of the samples) \u003cem\u003eMethylophilus\u003c/em\u003e compared with sections without (P=0.0162) (Figure 5).\u003c/p\u003e\n\u003c/div\u003e"},{"header":"Discussion","content":" \u003cdiv id=\"Sec8\" class=\"Section2\"\u003e \u003ch2\u003eElemental composition within the soil layers\u003c/h2\u003e \u003cp\u003eOverall, there were significant shifts in elemental composition between soil samples at Vindolanda specifically with increases in phosphorus, sulphur and iron within d1p and d2p preservation locations. Non-preservation locations d3i, d4t, c1 and c2 showed consistent elemental composition regardless of the different trench locations. Whilst preservation of artefacts is poor in d3i and d4t there is evidence of Roman occupation. However, there is no such evidence in c1 and c2 suggesting they sit outside of the settlement. Interestingly, there was little observable difference within pH values between any of the samples.\u003c/p\u003e \u003cp\u003eIron and sulphate reduction, denitrification and nitrate respiration are associated with anaerobic environments\u003csup\u003e\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e\u003c/sup\u003e, with dissimilatory Fe\u003csup\u003e3+\u003c/sup\u003e reduction being the most abundant anaerobic metabolism\u003csup\u003e\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e\u003c/sup\u003e. Anaerobic soil conditions are typically characterised by the accumulation of iron cations (Fe\u003csup\u003e2+\u003c/sup\u003e and Fe\u003csup\u003e3+\u003c/sup\u003e) as soil-available Fe is reduced. In aerobic conditions oxidation of Fe\u003csup\u003e2+\u003c/sup\u003e to Fe seems to correlate to a speeding up of microbial metabolism and a decomposition of organic matter both within field and pot experiments\u003csup\u003e\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e\u003c/sup\u003e. However, pXRF analysis used in the current study cannot differentiate between the different oxidative states of iron; it does however, suggest that there is an increase in the amount of iron present within these locations not just a change in redox state.\u003c/p\u003e \u003cp\u003ePhosphorus release is associated with increased iron reduction within anaerobic soil environments with the dynamics of phosphorus release and Fe\u003csup\u003e3+\u003c/sup\u003e reduction closely linked\u003csup\u003e\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e\u003c/sup\u003e. However, there is also a relationship between the cycling of iron and phosphorus with the concentration of sulphate found within waterlogged environments. Within such conditions sulphate reduction results in increases in phosphorus release\u003csup\u003e\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e\u003c/sup\u003e. The iron and phosphorus concentrations are higher at Vindolanda within locations showing preservation and could, therefore suggest why increased vivianite are observed here\u003csup\u003e\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e. Ratios of less than 1.5 sulphur:iron are needed for vivianite formation within marine sediments\u003csup\u003e\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e\u003c/sup\u003e; all soils used within this study met these conditions. However, an excess of phosphorus is also needed which is only seen in d1p and d2p.\u003c/p\u003e \u003cp\u003eCertain soil types are more likely to result in good archaeological preservation than others with anaerobic soils rich in organic material being more likely to show reduced biological oxidation of wood, leather and soft tissue\u003csup\u003e\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e\u003c/sup\u003e, partly due to the microorganism which flourish in such environments\u003csup\u003e\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e\u003c/sup\u003e. Reduced forms of phosphorus, sulphur and iron are more likely to be found within anaerobic soils compared with aerobic. Furthermore, the changes in moisture conditions at Vindolanda are linked to the periodic release of agricultural leachate which would be expected to increase soil degradation but may allow reaction with organic material within the soil, thus facilitating a natural \u0026lsquo;tanning\u0026rsquo; process\u003csup\u003e\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cdiv id=\"Sec9\" class=\"Section3\"\u003e \u003ch2\u003eThe Microbial community within the Vindolanda soil layers\u003c/h2\u003e \u003cp\u003ePreservation and interface locations d1p, d2p and d3i are distinct from other samples in terms of composition of phyla and genera. Shifts in bacterial phyla are correlated to location of sample but not pH values or elemental composition. Preservation locations are dominated by \u003cem\u003eFirmicutes, Bacteroidetes\u003c/em\u003e and \u003cem\u003eProteobacteria\u003c/em\u003e at a phylum level with LEfSe analysis suggesting that \u003cem\u003eFirmicutes, Epsilonproteobacteria, BRC1\u003c/em\u003e and \u003cem\u003eChloroflexi\u003c/em\u003e may be indicative of preservation. Non-preservation and control locations were typified by increased abundance of \u003cem\u003eAcidobacteria, Actinobacteria\u003c/em\u003e and \u003cem\u003ePlanctomycetes\u003c/em\u003e. Shannon diversity remained high within preservation layers and was comparable to the surrounding soil. Maintenance of diversity and evenness within good archaeological preserving environments has been shown previously in No. 1 Wangshanqiao Chu Tomb in Jingzhou, China\u003csup\u003e\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003ePrevious studies have shown that at the phylum level soils within the top approximately 20 cm are likely to be dominated by \u003cem\u003eAcidobacteria, Actinobacteria\u003c/em\u003e and \u003cem\u003eAlphaproteobacteria\u003c/em\u003e with \u003cem\u003eFirmicutes\u003c/em\u003e and \u003cem\u003eSpirochaetes\u003c/em\u003e being more abundant below 50cm \u003csup\u003e\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u003c/sup\u003e. \u003cem\u003eAcidobacteria\u003c/em\u003e and \u003cem\u003eActinobacteria\u003c/em\u003e are increased within d4t, c1 and c2 compared with d1p, d2p and d3i. \u003cem\u003eActinobacteria\u003c/em\u003e are widely recognised as being crucial to aerobic degradation of organic matter particularly that which is cellulose based\u003csup\u003e\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e,\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e\u003c/sup\u003e. \u003cem\u003eAcidobacteria\u003c/em\u003e are a diverse phylum which are mostly heterotrophic in metabolism and can degrade a wide range of carbohydrate sources\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e. Members of the \u003cem\u003eAcidobacteria\u003c/em\u003e group 1\u0026ndash;4 have been linked with degradation of cellulose based material\u003csup\u003e\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e\u003c/sup\u003e and are seen in significant abundance within layers 4\u0026ndash;6 of this study. Interestingly they are greatly reduced within the preservation layers despite the pH values remaining consistent. This is probably due to the anaerobic conditions present within these layers as the majority of \u003cem\u003eAcidobacteria\u003c/em\u003e are aerobic\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003e \u003cem\u003eFirmicutes, Bacteroidetes\u003c/em\u003e and \u003cem\u003eProteobacteria\u003c/em\u003e were significantly increased within the preservation layers. However, this is most likely a result of the anaerobic conditions found rather than these phyla not being decomposers. \u003cem\u003eProteobacteria, Firmicutes\u003c/em\u003e and \u003cem\u003eBacteroidetes\u003c/em\u003e are the most dominant phyla in solid waste decomposition within anaerobic phases\u003csup\u003e\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e\u003c/sup\u003e with \u003cem\u003eFirmicutes\u003c/em\u003e and \u003cem\u003eBacteroidetes\u003c/em\u003e proposed as the main microorganisms involved in cellulose degradation \u003csup\u003e\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e\u003c/sup\u003e. There are very few studies assessing microbiology within archaeological sites. However, those available correspond with our findings of consistent microbial diversity between layers which show preservation and those which do not \u003csup\u003e\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e,\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e\u003c/sup\u003e. Siles (2018) also found varied soil profiles that were common to soil environments and recorded no dominant \u0026lsquo;rare\u0026rsquo; microbial taxa as clear indicators of preservation\u003csup\u003e\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eAt the genus level preservation locations showed increased abundance of \u003cem\u003eSulficurvum\u003c/em\u003e, \u003cem\u003eFlavobacterium\u003c/em\u003e and \u003cem\u003eMethylophilus\u003c/em\u003e with additional genera \u003cem\u003eTuricibacter\u003c/em\u003e, \u003cem\u003eFilomicrobium\u003c/em\u003e, \u003cem\u003eSphingobium\u003c/em\u003e, \u003cem\u003ePusillimonas\u003c/em\u003e, \u003cem\u003eVogesella\u003c/em\u003e, \u003cem\u003eArcobacter\u003c/em\u003e and \u003cem\u003eThiovirga\u003c/em\u003e being identified as indicative of preservation by LEfSe. These microorganisms are not typically indicative of depth or oxygen availability as changes at phylum level were. Instead \u003cem\u003eSulficurvum\u003c/em\u003e, \u003cem\u003eFlavobacterium\u003c/em\u003e and \u003cem\u003eMethylophilus\u003c/em\u003e are all capable of degrading polycyclic aromatic hydrocarbons (PAHs) and are found to be encouraged by increased concentrations of iron and sulphur\u003csup\u003e\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e,\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e\u003c/sup\u003e. \u003cem\u003eFlavobacterium\u003c/em\u003e and \u003cem\u003eMethylophilus\u003c/em\u003e are aerobic, whilst \u003cem\u003eSulfuricurvum\u003c/em\u003e are facultative anaerobes which are major contributors to the sulphur cycle given their abilities to oxidise sulphur and are observed in the presence of excess methane\u003csup\u003e\u003cspan additionalcitationids=\"CR31\" citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e\u003c/sup\u003e. This coupled with the changes in soil chemistry involved suggests an unusual soil environment within the preservation layers at Vindolanda.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec10\" class=\"Section3\"\u003e \u003ch2\u003eMicrobial community with vivianite\u003c/h2\u003e \u003cp\u003eInterestingly, further investigation of samples with visible vivianite showed that these soils contained marked increases in \u003cem\u003eMethylophilus\u003c/em\u003e compared with adjacent soils. \u003cem\u003eMethylophilus\u003c/em\u003e are methylotrophic \u003cem\u003eBetaproteobacteria\u003c/em\u003e commonly occurring in activated sludge, mud, and aquatic habitats, and are capable of utilizing methanol or methylamine as a sole source of carbon and energy\u003csup\u003e\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e\u003c/sup\u003e. Members of the genus are thought to be strictly aerobic and have been extracted from various soils\u003csup\u003e\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e\u003c/sup\u003e. Usually \u003cem\u003eMethylophilus\u003c/em\u003e would be expected to be found in low abundance within soil environments. However, previous studies have found increased amounts of CO\u003csub\u003e2\u003c/sub\u003e within soil samples can lead to proliferation of \u003cem\u003eMethylophilus\u003c/em\u003e presumably due to an increased abundance of methane\u003csup\u003e\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003e \u003cem\u003eMethylophilus\u003c/em\u003e is capable of solubilizing inorganic phosphorus\u003csup\u003e\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e\u003c/sup\u003e and produces siderophores and indole acetic acid\u003csup\u003e\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e\u003c/sup\u003e. Studies have shown that \u003cem\u003eMethylophilus\u003c/em\u003e species are capable of reducing iron species such as ferrihydrite within conditions of increased methanol with evidence that the ferrihydrite was transformed into crystalline magnetite, a mixed Fe\u003csup\u003e3+\u003c/sup\u003e oxide\u003csup\u003e35\u003c/sup\u003e. Furthermore, other microorganisms capable of producing magnetite can also produce vivianite in phosphate rich environments\u003csup\u003e\u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eStudies in Brazil have shown that, with the exception of environments that are high in sulphur, phosphorus and iron the production of iron species such as magnetite and haematite are relatively common while the production of vivianite is more restricted\u003csup\u003e\u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e\u003c/sup\u003e. These conditions are similar to those we have identified within the preservation layers at Vindolanda. Additionally, analysis within the South China Sea showed that concentrations of iron-bound phosphorus and organic-phosphorus were influenced by sulfate-driven and metal-driven anaerobic oxidation of methane resulting in the release of phosphates and Fe\u003csup\u003e2+\u003c/sup\u003e and a subsequent significantly increased formation of vivianite \u003csup\u003e\u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e "},{"header":"Conclusion","content":" \u003cp\u003eOur findings taken together suggest that the soil at Vindolanda contains distinctive chemical and microbial profiles associated with layers that show increased preservation of organic/inorganic artefacts. These shifts are not simply the result of changes in depth and oxygen availability but suggest that the specific conditions potentially allow the microbial production of vivianite which may lead to enhanced preservation of artefacts. The role of the occupation layers as phosphorus sinks is strongly linked to the abundance of vivanite. More research is needed to clearly link the presence of \u003cem\u003eMethylophilus\u003c/em\u003e with vivianite production.\u003c/p\u003e"},{"header":"Methods","content":"\u003cp\u003e\u003cstrong\u003eSampling and archaeological context\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eSoil was sampled from two trenches within the Vindolanda site (54.9911\u0026deg;N 2.3608\u0026deg;W). The external surface was scraped back by approximately 5cm. Samples were then taken in triplicate at four equal depths (d1p, d2p, d3i and d4t), each representing a separate archaeological context (see supplementary Fig.\u0026nbsp;2). The top 40cm of the trench had was exposed for approximately 3 months prior to the sampling taking place down to the interface level of d4t. The lower 1m section was excavated to the level of d1p four weeks prior to the sampling taking place. D1p and d2p correlate to the deepest layers and are associated with the best quality preservation. D3i and D4t are associated with lower quality preservation while d3i forming an interface (i) layer. Soil was sampled in triplicate (A-B) from a control (c) trench showing no evidence occupation or archaeological preservation, at two depths (c1-c2) (see supplementary Fig.\u0026nbsp;3).\u003c/p\u003e\n\u003cp\u003eIn addition, further samples were taken from a rampart section, as detailed in Taylor \u003cem\u003eet al\u003c/em\u003e 2019. The rampart sample was sampled in triplicate from the vivianite visible layer (blue) as well as triplicate samples within 10 cm which showed no visible vivianite formation.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eMoisture\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTriplicate soils were weighed before and after drying in an oven overnight at 105\u0026deg;C to determine the moisture content of the original sample using Eq.\u0026nbsp;1.\u003c/p\u003e\n\u003cp\u003eEquation 1: 𝑤𝑒𝑡 𝑤𝑒𝑖𝑔ℎ𝑡 (𝑔) \u0026minus; 𝑑𝑟𝑦 𝑤𝑒𝑖𝑔ℎ𝑡 (𝑔) / 𝑑𝑟𝑦 𝑤𝑒𝑖𝑔ℎ𝑡 (𝑔) \u0026times; 100 = 𝑚𝑜𝑖𝑠𝑡𝑢𝑟𝑒 𝑐𝑜𝑛𝑡𝑒𝑛𝑡 %\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003epH\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eDried soils (10g) were placed into plastic cups with distilled water added (20ml) for a 1:2 proportion. The solution was stirred with a glass rod and settled for 30 minutes. pH values was analysed using an Advanced Jenway 3510 pH meter calibrated with commercial buffers of pH 4, 7 and 10. The electrode probe was rinsed between each sample.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003epXRF\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eSoils were analysed following the method provided by Williams (2020)\u003csup\u003e39\u003c/sup\u003e. The dried soil samples were homogenised with mortar and pestle for 140 seconds and sieved to 2mm. Soil were loaded into XRF sample cups (SPEX CertiPrep\u0026trade; 3529) and covered with 5 \u0026micro;m polypropylene thin-film (SPEX\u0026trade; SamplePrep 3520 window film). A Thermo Niton\u0026trade; XL3t GOLDD\u0026thinsp;+\u0026thinsp;pXRF with an Ag anode (6-50kV, 0-200 \u0026micro;A max X-ray tube) was used to analyse the soil. The pXRF was warmed up, system checked against the internal 1\u0026frac14; Cr\u0026ndash;\u0026frac12; Mo coupon, and scanned against a blank standard and NIST 2709a standard (y\u0026thinsp;=\u0026thinsp;0.9674x -0.0089, \u003cem\u003er\u003c/em\u003e\u003csup\u003e2\u003c/sup\u003e\u0026thinsp;=\u0026thinsp;0.9998). The NIST 2709a is a San Joaquin soil, used because it was certified for most elements of interest within the single calibration sample. The pXRF was periodically reset and system checked to account for drift. Soil samples were analysed using Mining mode (fundamental parameters), with 30-second scans for the main filter (50 kV, \u0026le;\u0026thinsp;50 \u0026micro;A), low filter (20 kV, \u0026le;\u0026thinsp;100 \u0026micro;A) and high filter (50 kV, \u0026le;\u0026thinsp;40 \u0026micro;A), and a 60-second scan for the light filter (6 kV, \u0026le;\u0026thinsp;200 \u0026micro;A).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eDNA extraction\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eDNA was extracted from 0.25 g of each soil sample using the Qiagen DNeasy PowerSoil Kit following manufacturer\u0026rsquo;s instructions. DNA quality and yield were measured using NanoDrop 1000 spectrophotometer (V3.8.1, ThermoFisher) and agarose gel electrophoresis. DNA was stored at -80\u0026deg;C prior to further analysis.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eNext Generation Sequencing\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eBacterial community DNA sequencing was made (NU-OMICs, Northumbria University, Newcastle Upon Tyne, U.K.) by targeting the V4 region of the 16S rRNA gene according to Kozich (2013)\u003csup\u003e40\u003c/sup\u003e. The raw sequencing reads were processed in FASTQ format and analysed with Mothur software package (version 1.36.1) (University of Michigan, U.S.A.). UCHIME was used to quality check and filter the FASTA formatted sequences for chimeras. These were aligned to the SILVA reference, and taxonomic identification of the reads was assessed by assigning sequences to operational taxonomic units (OTUs) using Ribosomal Database Project (RDP) classification. PCR negative controls were run and sequenced in parallel to the experimental samples. The OTUs recorded in negative controls and samples were excluded from further analysis. Non-bacterial sequences (Archaea/Eukaryota/Unclassified) were discarded. OTUs less than 3% were classified as rare taxa with these and the unclassified OTUs omitted from the plots.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eStatistical analysis\u003c/strong\u003e\u003c/p\u003e\n\u003cdiv\u003e\n\u003cp\u003eData were analysed in R 4.0.2\u003csup\u003e41\u003c/sup\u003e using the car, DescTools and PMCMR packages\u003csup\u003e42\u0026ndash;44\u003c/sup\u003e. Normality and variance failed the Q-Q plot, Shapiro-Wilk\u0026rsquo;s and Levene\u0026rsquo;s assumption tests. Data were analysed using Kruskal-Wallis rank sum tests, followed by Nemenyi tests with Chi-Square distribution to identify significant differences between burial layers. Tukey box plots of the median and quartiles were produced using the ggplot2 package\u003csup\u003e45\u003c/sup\u003e to show the change in raw elemental concentration with increasing burial depth. The tidypaleo package in tidyverse\u003csup\u003e46\u003c/sup\u003e was used to produce stratigraphic plots of the moisture content and pH values against burial depth.\u003c/p\u003e\n\u003cp\u003eOperational taxonomic unit (OUT) tables and metadata were analysed using R version 4.0.2. Samples were faceted to rarefy, relative abundance data was compiled and most abundant phyla and genera were sought. Within the readr package\u003csup\u003e47\u003c/sup\u003e significant differences were identified using ANOVA and TukeyHSD tests and Spearman Correlation was used to determine significant correlation between variables. Figures were created using gplots\u003csup\u003e48\u003c/sup\u003e. Taxa indicative of preservation were identified, and histograms generated following Linear Discriminant Analysis Effect Size (LEfSe) analysis completed in Galaxy\u003csup\u003e49\u003c/sup\u003e.\u003c/p\u003e\n\u003c/div\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eAcknowledgements\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors thank Dr Andrew Birley and the Vindolanda Trust for providing access to the archaeological samples. There are no conflicts of interest to declare. The data that support the findings of this study are available in the additional supporting information.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthor contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eGT and CO wrote the first draft of the manuscript and participated in sample collection and data analysis. RW and HH prepared and analysed the pXRF samples. AN organized and conducted the Next Generation Sequencing. AB and EG was involved in the fieldwork and sample collection. TKR sample preparation for Next generation sequencing.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eBirley, R. \u003cem\u003eVindolanda: A Roman frontier post on Hadrian's Wall\u003c/em\u003e (Thames and Hudson, London, 1977).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTaylor, G., Hefford, R. J. W., Birley, A. \u0026amp; Huntley, J. P. Identifying the \u0026ldquo;blue substance\u0026rdquo; at the Roman site of Vindolanda, Northumberland. \u003cem\u003eJournal of Archaeological Science: Reports.\u003c/em\u003e \u003cb\u003e24\u003c/b\u003e, 582\u0026ndash;587 (2019).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHolden, J. \u003cem\u003eet al.\u003c/em\u003e Hydrological controls of in situ preservation of waterlogged archaeological deposits. \u003cem\u003eEarth-Science Reviews.\u003c/em\u003e \u003cb\u003e78\u003c/b\u003e (1\u0026ndash;2), 59\u0026ndash;83 (2006).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCaple, C. Reburial of waterlogged wood, the problems and potential of this conservation technique. \u003cem\u003eInternational Biodeterioration and Biodegradation.\u003c/em\u003e \u003cb\u003e34\u003c/b\u003e (1), 61\u0026ndash;72 (1994).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMcGowan, G. \u0026amp; Prangnell, J. The significance of vivianite in archaeological settings. \u003cem\u003eGeoarchaeology.\u003c/em\u003e \u003cb\u003e21\u003c/b\u003e (1), 93\u0026ndash;111 (2006).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eS\u0026aacute;nchez-Rom\u0026aacute;n, M., Puente-S\u0026aacute;nchez, F., Parro, V. \u0026amp; Amils, R. Nucleation of Fe-rich phosphates and carbonates on microbial cells and exopolymeric substances. \u003cem\u003eFrontiers in Microbiology.\u003c/em\u003e 2015;6(SEP).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMiot, J. \u003cem\u003eet al.\u003c/em\u003e Transformation of vivianite by anaerobic nitrate-reducing iron-oxidizing bacteria. \u003cem\u003eGeobiology.\u003c/em\u003e \u003cb\u003e7\u003c/b\u003e (3), 373\u0026ndash;384 (2009).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eVolkland, H. P. \u003cem\u003eet al.\u003c/em\u003e Bacterial phosphating of mild (unalloyed) steel. \u003cem\u003eApplied and Environmental Microbiology.\u003c/em\u003e \u003cb\u003e66\u003c/b\u003e (10), 4389\u0026ndash;4395 (2000).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBlanchette, R. A., Nilsson, T., Daniel, G. \u0026amp; Abad, A. Biological degradation of wood.\u003cem\u003eAmerican Chemical Society.\u003c/em\u003e1990:141\u0026ndash;174.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBj\u0026ouml;rdal, C. G., Nilsson, T. \u0026amp; Daniel, G. Microbial decay of waterlogged archaeological wood found in Sweden. Applicable to archaeology and conservation. \u003cem\u003eInternational Biodeterioration and Biodegradation.\u003c/em\u003e \u003cb\u003e43\u003c/b\u003e (1\u0026ndash;2), 63\u0026ndash;73 (1999).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBj\u0026ouml;rdal, C. G. Microbial degradation of waterlogged archaeological wood. \u003cem\u003eJournal of Cultural Heritage.\u003c/em\u003e \u003cb\u003e13\u003c/b\u003e (3 SUPPL.), S118\u0026ndash;S122 (2012).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJessen, G. L. \u003cem\u003eet al.\u003c/em\u003e Hypoxia causes preservation of labile organic matter and changes seafloor microbial community composition (Black Sea).\u003cem\u003eSci Adv.\u003c/em\u003e2017;3(2).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eTiedje, J. M., Sexstone, A. J., Parkin, T. B. \u0026amp; Revsbech, N. P. Anaerobic processes in soil. \u003cem\u003ePlant and Soil.\u003c/em\u003e \u003cb\u003e76\u003c/b\u003e (1\u0026ndash;3), 197\u0026ndash;212 (1984).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKeiluweit, M., Wanzek, T., Kleber, M., Nico, P. \u0026amp; Fendorf, S. Anaerobic microsites have an unaccounted role in soil carbon stabilization.\u003cem\u003eNat Commun.\u003c/em\u003e2017;8(1).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHall, S. J. \u0026amp; Silver, W. L. Iron oxidation stimulates organic matter decomposition in humid tropical forest soils. \u003cem\u003eGlobal Change Biol.\u003c/em\u003e \u003cb\u003e19\u003c/b\u003e (9), 2804\u0026ndash;2813 (2013).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eForsmann, D. M. \u0026amp; Kjaergaard, C. Phosphorus release from anaerobic peat soils during convective discharge - Effect of soil Fe: P molar ratio and preferential flow. \u003cem\u003eGeoderma.\u003c/em\u003e \u003cb\u003e223\u0026ndash;225\u003c/b\u003e (1), 21\u0026ndash;32 (2014).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLehtoranta, J., Ekholm, P., Wahlstr\u0026ouml;m, S., Tallberg, P. \u0026amp; Uusitalo, R. Labile organic carbon regulates phosphorus release from eroded soil transported into anaerobic coastal systems. \u003cem\u003eAmbio.\u003c/em\u003e \u003cb\u003e44\u003c/b\u003e (2), 263\u0026ndash;273 (2015).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRothe, M., Kleeberg, A. \u0026amp; Hupfer, M. The occurrence, identification and environmental relevance of vivianite in waterlogged soils and aquatic sediments. \u003cem\u003eEarth-Science Reviews.\u003c/em\u003e \u003cb\u003e158\u003c/b\u003e, 51\u0026ndash;64 (2016).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKibblewhite, M., T\u0026oacute;th, G. \u0026amp; Hermann, T. Predicting the preservation of cultural artefacts and buried materials in soil. \u003cem\u003eScience of the Total Environment.\u003c/em\u003e \u003cb\u003e529\u003c/b\u003e, 249\u0026ndash;263 (2015).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eXu, J., Wei, Y., Jia, H., Xiao, L. \u0026amp; Gong, D. A new perspective on studying burial environment before archaeological excavation: Analyzing bacterial community distribution by high-throughput sequencing.\u003cem\u003eScientific Reports.\u003c/em\u003e2017;7.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDouterelo, I., Goulder, R. \u0026amp; Lillie, M. Soil microbial community response to land-management and depth, related to the degradation of organic matter in English wetlands: Implications for the in situ preservation of archaeological remains. \u003cem\u003eAppl Soil Ecol.\u003c/em\u003e \u003cb\u003e44\u003c/b\u003e (3), 219\u0026ndash;227 (2010).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWang, C. \u003cem\u003eet al.\u003c/em\u003e Metagenomic analysis of microbial consortia enriched from compost: New insights into the role of Actinobacteria in lignocellulose decomposition.\u003cem\u003eBiotechnol Biofuels.\u003c/em\u003e2016;9(1).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePankratov, T. A., Dedysh, S. N. \u0026amp; Zavarzin, G. A. The leading role of actinobacteria in aerobic cellulose degradation in Sphagnum peat bogs. \u003cem\u003eDokl Biol Sci.\u003c/em\u003e \u003cb\u003e410\u003c/b\u003e (1), 428\u0026ndash;430 (2006).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKielak, A. M., Barreto, C. C., Kowalchuk, G. A., van Veen, J. A. \u0026amp; Kuramae, E. E. The ecology of Acidobacteria: Moving beyond genes and genomes. \u003cem\u003eFrontiers in Microbiology.\u003c/em\u003e 2016;7(MAY).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eEichorst, S. A., Kuske, C. R. \u0026amp; Schmidt, T. M. Influence of plant polymers on the distribution and cultivation of bacteria in the phylum acidobacteria. \u003cem\u003eApplied and Environmental Microbiology.\u003c/em\u003e \u003cb\u003e77\u003c/b\u003e (2), 586\u0026ndash;596 (2011).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYang, S. \u0026amp; Song, L. Succession of bacterial community structure and metabolic function during solid waste decomposition.\u003cem\u003eBioresource Technology.\u003c/em\u003e2019;291.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSiles, J. A., \u0026Ouml;hlinger, B., Cajthaml, T., Kistler, E. \u0026amp; Margesin, R. Characterization of soil bacterial, archaeal and fungal communities inhabiting archaeological human-impacted layers at Monte Iato settlement (Sicily, Italy).\u003cem\u003eScientific Reports.\u003c/em\u003e2018;8(1).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSun, W. \u003cem\u003eet al.\u003c/em\u003e Profiling microbial community in a watershed heavily contaminated by an active antimony (Sb) mine in Southwest China. \u003cem\u003eScience of the Total Environment.\u003c/em\u003e \u003cb\u003e550\u003c/b\u003e, 297\u0026ndash;308 (2016).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHaaijer, S. C. M. \u003cem\u003eet al.\u003c/em\u003e Bacteria associated with iron seeps in a sulfur-rich, neutral pH, freshwater ecosystem. \u003cem\u003eISME Journal.\u003c/em\u003e \u003cb\u003e2\u003c/b\u003e (12), 1231\u0026ndash;1242 (2008).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi, X. \u003cem\u003eet al.\u003c/em\u003e Metagenome-assembled genome sequence of sulfuricurvum sp. strain IAE1, isolated from a 4-chlorophenol-degrading consortium.\u003cem\u003eMicro Res Ann.\u003c/em\u003e2019;8(31).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHubert, C. R. J. \u003cem\u003eet al.\u003c/em\u003e Massive dominance of Epsilonproteobacteria in formation waters from a Canadian oil sands reservoir containing severely biodegraded oil. \u003cem\u003eEnvironmental Microbiology.\u003c/em\u003e \u003cb\u003e14\u003c/b\u003e (2), 387\u0026ndash;404 (2012).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMadhaiyan, M., Poonguzhali, S., Kwon, S. W. \u0026amp; Sa, T. M. Methylophilus rhizosphaerae sp. nov., a restricted facultative methylotroph isolated from rice rhizosphere soil. \u003cem\u003eInt J Syst Evol Microbiol.\u003c/em\u003e \u003cb\u003e59\u003c/b\u003e (11), 2904\u0026ndash;2908 (2009).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDoronina, N., Kaparullina, E. \u0026amp; Trotsenko, Y.. The Family Methylophilaceae. In: Rosenberg E, DeLong EF, Lory S, Stackebrandt E, Thompson F, eds. The Prokaryotes: Alphaproteobacteria and Betaproteobacteria. Berlin, Heidelberg:Springer Berlin Heidelberg; 2014:869\u0026ndash;880.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMa, J. \u003cem\u003eet al.\u003c/em\u003e Short-term effects of CO2 leakage on the soil bacterial community in a simulated gas leakage scenario. \u003cem\u003ePeerJ.\u003c/em\u003e \u003cb\u003e2017\u003c/b\u003e, 11 (2017).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYang, Y. \u003cem\u003eet al.\u003c/em\u003e Extracellular electron transfer of Methylophilus methylotrophs. \u003cem\u003eProcess Biochem.\u003c/em\u003e \u003cb\u003e94\u003c/b\u003e, 313\u0026ndash;318 (2020).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eVeeramani, H. \u003cem\u003eet al.\u003c/em\u003e Products of abiotic U(VI) reduction by biogenic magnetite and vivianite. \u003cem\u003eGeochimica et Cosmochimica Acta.\u003c/em\u003e \u003cb\u003e75\u003c/b\u003e (9), 2512\u0026ndash;2528 (2011).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eVolpato, J. Relation between geochemistry and minerals (magnetite, hematite, pyrite and vivianite) in contamineted estuarine sediments, Santos, Sao Paulo, Brazil., Universidade Estadual de Campinas; 2015.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWu, D. \u003cem\u003eet al.\u003c/em\u003e Zone of metal-driven anaerobic oxidation of methane is an important sink for phosphorus in the Taixinan Basin, South China Sea. \u003cem\u003eMarine Geology.\u003c/em\u003e \u003cb\u003e427\u003c/b\u003e, 1\u0026ndash;10 (2020).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWilliams, R., Taylor, G. \u0026amp; Orr, C. pXRF method development for elemental analysis of archaeological soil. \u003cem\u003eArchaeometry.\u003c/em\u003e \u003cb\u003e62\u003c/b\u003e (6), 1145\u0026ndash;1163 (2020).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKozich, J., Westcott, S. L., Baxter, N. T., Highlander, S. K. \u0026amp; Schloss, P. D. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the miseq illumina sequencing platform. \u003cem\u003eApplied and Environmental Microbiology.\u003c/em\u003e \u003cb\u003e79\u003c/b\u003e (17), 5112\u0026ndash;5120 (2013).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003e\u003cem\u003eR: A language and environment for statistical computing\u003c/em\u003e [computer program]. Vienna:R Foundation for Statistical Computing; 2020.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSignorell, A. DescTools: Tools for descriptive statistics. 2016.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePohlert, T. The pairwise multiple comparison of mean ranks package (PMCMR). 2015.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFox, J. \u0026amp; Weisberg, S. An {R} companion to applied regression. 2011.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWickham, H. ggplot2: Elegant graphics for data analysis. 2016.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWickham, H. \u003cem\u003eet al.\u003c/em\u003e Welcome to the Tidyverse. \u003cem\u003eThe journal of open source software.\u003c/em\u003e \u003cb\u003e4\u003c/b\u003e (43), 1686 (2019).\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWickham, H., Hester, J. \u0026amp; Francois, R. readr: read rectangular Text Data. R package version 1.3.1.. 2018; \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://CRAN.R-project.org/package=readr\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWarnes, G. \u003cem\u003eet al.\u003c/em\u003e L, Gplots:Various R Programming Tools for Plotting Data. R package 3.03.. 2020; \u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://CRAN.R-project.org/package=gplots\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSegata, N. \u003cem\u003eet al.\u003c/em\u003e Metagenomic biomarker discovery and explanation.\u003cem\u003eGenome Biol.\u003c/em\u003e2011;12(6).\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"scientific-reports","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"scirep","sideBox":"Learn more about [Scientific Reports](http://www.nature.com/srep/)","snPcode":"","submissionUrl":"","title":"Scientific Reports","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Scientific Reports","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"Vindolanda, Anaerobic, Amplicon Sequencing, Vivianite","lastPublishedDoi":"10.21203/rs.3.rs-244364/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-244364/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eWaterlogged burial conditions impact upon artefact preservation. One major determinant of preservation is presence and behaviour of microorganisms, however, unravelling the mechanisms is challenging. In this study, we analysed elemental composition, bacterial diversity and community structure from excavation trenches at the Roman Site of Vindolanda, Northumberland, UK, using pXRF and 16s rRNA gene amplicon sequencing. Excavation trenches provide layers which can be classed as phosphorus sinks and which are strongly linked to the abundance of vivanite. The results indicated that microbial communities were dominant in \u003cem\u003eFirmicutes, Bacteroidetes\u003c/em\u003e and \u003cem\u003eProteobacteria\u003c/em\u003e at a phylum level. Samples which also had visible vivianite presence showed that there were marked increases in \u003cem\u003eMethylophilus\u003c/em\u003e. \u003cem\u003eMethylophilus\u003c/em\u003e might be associated with favourable preservation in these anaerobic conditions. More research is needed to clearly link the presence of \u003cem\u003eMethylophilus\u003c/em\u003e with vivianite production. The study emphasises the need for further integration of chemical and microbiome approaches, especially in good preservation areas, to explore microbial and chemical degradation mechanisms.\u003c/p\u003e","manuscriptTitle":"Unique Chemical Parameters and Microbial Activity Lead to Increased Archaeological Preservation within Vindolanda, UK.","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2021-03-10 14:58:07","doi":"10.21203/rs.3.rs-244364/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Major revision","date":"2021-04-08T07:50:30+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2021-03-30T02:41:02+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"b71933d5-01e4-48df-804e-589596497c06","date":"2021-03-19T11:55:10+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2021-03-08T19:20:27+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2021-03-03T18:41:22+00:00","index":"","fulltext":""},{"type":"editorInvited","content":"","date":"2021-03-03T15:30:00+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2021-03-03T15:26:12+00:00","index":"","fulltext":""},{"type":"submitted","content":"Scientific Reports","date":"2021-02-15T22:18:23+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"scientific-reports","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"scirep","sideBox":"Learn more about [Scientific Reports](http://www.nature.com/srep/)","snPcode":"","submissionUrl":"","title":"Scientific Reports","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Scientific Reports","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"0acb3f7b-1684-4a3f-bf96-e0e13428475d","owner":[],"postedDate":"March 10th, 2021","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[{"id":2872074,"name":"General Microbiology"},{"id":2872075,"name":"Molecular Biology"}],"tags":[],"updatedAt":"2021-08-22T15:19:40+00:00","versionOfRecord":{"articleIdentity":"rs-244364","link":"https://doi.org/10.1038/s41598-021-94853-7","journal":{"identity":"scientific-reports","isVorOnly":false,"title":"Scientific Reports"},"publishedOn":"2021-08-04 15:04:57","publishedOnDateReadable":"August 4th, 2021"},"versionCreatedAt":"2021-03-10 14:58:07","video":"","vorDoi":"10.1038/s41598-021-94853-7","vorDoiUrl":"https://doi.org/10.1038/s41598-021-94853-7","workflowStages":[]},"version":"v1","identity":"rs-244364","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-244364","identity":"rs-244364","version":["v1"]},"buildId":"cBFmMYwuxLRRLfASyISRj","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}