Fossilization of melanosomes via sulfurization

Abstract Fossil melanin granules (melanosomes) are an important resource for inferring the evolutionary history of colour and its functions in animals. The taphonomy of melanin and melanosomes, however, is incompletely understood. In particular, the chemical processes responsible for melanosome preservation have not been investigated. As a result, the origins of sulfur‐bearing compounds in fossil melanosomes are difficult to resolve. This has implications for interpretations of original colour in fossils based on potential sulfur‐rich phaeomelanosomes. Here we use pyrolysis gas chromatography mass spectrometry (Py‐GCMS), fourier transform infrared spectroscopy (FTIR) and time of flight secondary ion mass spectrometry (ToF‐SIMS) to assess the mode of preservation of fossil microstructures, confirmed as melanosomes based on the presence of melanin, preserved in frogs from the Late Miocene Libros biota (NE Spain). Our results reveal a high abundance of organosulfur compounds and non‐sulfurized fatty acid methyl esters in both the fossil tissues and host sediment; chemical signatures in the fossil tissues are inconsistent with preservation of phaeomelanin. Our results reflect preservation via the diagenetic incorporation of sulfur, i.e. sulfurization (natural vulcanization), and other polymerization processes. Organosulfur compounds and/or elevated concentrations of sulfur have been reported from melanosomes preserved in various invertebrate and vertebrate fossils and depositional settings, suggesting that preservation through sulfurization is likely to be widespread. Future studies of sulfur‐rich fossil melanosomes require that the geochemistry of the host sediment is tested for evidence of sulfurization in order to constrain interpretations of potential phaeomelanosomes and thus of original integumentary colour in fossils.

Abstract: Fossil melanin granules (melanosomes) are an important resource for inferring the evolutionary history of colour and its functions in animals. The taphonomy of melanin and melanosomes, however, is incompletely understood. In particular, the chemical processes responsible for melanosome preservation have not been investigated. As a result, the origins of sulfur-bearing compounds in fossil melanosomes are difficult to resolve. This has implications for interpretations of original colour in fossils based on potential sulfur-rich phaeomelanosomes. Here we use pyrolysis gas chromatography mass spectrometry (Py-GCMS), fourier transform infrared spectroscopy (FTIR) and time of flight secondary ion mass spectrometry (ToF-SIMS) to assess the mode of preservation of fossil microstructures, confirmed as melanosomes based on the presence of melanin, preserved in frogs from the Late Miocene Libros biota (NE Spain). Our results reveal a high abundance of organosulfur compounds and non-sulfurized fatty acid methyl esters in both the fossil tissues and host sediment; chemical signatures in the fossil tissues are inconsistent with preservation of phaeomelanin. Our results reflect preservation via the diagenetic incorporation of sulfur, i.e. sulfurization (natural vulcanization), and other polymerization processes. Organosulfur compounds and/or elevated concentrations of sulfur have been reported from melanosomes preserved in various invertebrate and vertebrate fossils and depositional settings, suggesting that preservation through sulfurization is likely to be widespread. Future studies of sulfur-rich fossil melanosomes require that the geochemistry of the host sediment is tested for evidence of sulfurization in order to constrain interpretations of potential phaeomelanosomes and thus of original integumentary colour in fossils.
T H E study of fossil colour is an emerging field in palaeontology and evolutionary biology. Much research to date has focused on the preservation of morphological and chemical evidence of melanin in fossils and its implications for reconstructions of the original colouration of ancient animals (Clarke et al. 2010;Li et al. 2010Li et al. , 2012Zhang et al. 2010;Carney et al. 2012;Lindgren et al. 2014), the evolutionary history of visual signalling (Clarke et al. 2010;Li et al. 2010Li et al. , 2012Zhang et al. 2010;Carney et al. 2012) and of physiology (Li et al. 2012) and the fidelity of the fossil record (Glass et al. 2012(Glass et al. , 2013Lindgren et al. 2014). Micron-sized granules of melanin, i.e. melanosomes, have been reported from various vertebrate and invertebrate fossils ranging in age from the Devonian to the Eocene, i.e. feathers (Vinther et al. 2008;Barden et al. 2011;Knight et al. 2011;Carney et al. 2012 (Glass et al. 2012(Glass et al. , 2013. These fossil microstructures are superficially similar to some bacteria in terms of their morphology and size (Vinther et al. 2008;Zhang et al. 2010;Knight et al. 2011) but can be most plausibly interpreted as fossil melanosomes based on details of their arrangement, location, precise context within soft tissues, and in particular their geochemistry (for further discussion see Barden et al. 2014;Moyer et al. 2014). Certain trace elements (e.g. Cu, Zn, Ca, and S) have been proposed as biomarkers for melanin (Wogelius et al. 2011;Manning et al. 2013). Fossil melanosomes can also retain organic geochemical evidence for the melanin molecule itself (Glass et al. 2012(Glass et al. , 2013Lindgren et al. 2012Lindgren et al. , 2014. Somewhat surprisingly, despite such intensive study, current understanding of the taphonomy of melanin and melanosomes is limited. In particular, no mechanism has been proposed to explain the widespread preservation of melanin and melanosomes in fossils. Maturation experiments using feathers from extant birds (McNamara et al. 2013) and other taxa (Colleary et al. 2015) have shown that melanosome size is altered during diagenesis. Chemical analyses of melanosomes preserved in fossil fish have suggested that melanin is preferentially preserved while other original chemical components of the organelles, such as proteins and lipids, are almost completely degraded (Lindgren et al. 2012 (Riley 1997). It is unclear, however, whether such intrinsic chemical properties of the melanin molecule alone are responsible for widespread preservation of melanin and melanosomes in fossils, or whether other geochemical processes contribute to preservation.
There is increasing evidence that the organic preservation of many non-biomineralized tissues, such as the cuticle of fossil arthropods and plants and the periderm of graptolites, reflects in situ polymerization of authochthonous lipids during diagenesis (Briggs 1999;Stankiewicz et al. 2000;Gupta et al. 2006Gupta et al. , 2009Manning et al. 2009Manning et al. , 2014Dutta et al. 2010;Gupta & Briggs 2011;Barden et al. 2011;Edwards et al. 2014). This process generates an aliphatic component similar to Type I and II kerogens (Briggs 1999) that can be identified in fossils, typically as a series of alkane/alkene doublets. Indeed, such aliphatic compounds are present in melanosome-bearing soft tissues from various fossils (Barden et al. 2011(Barden et al. , 2014Glass et al. 2013), suggesting that polymerization processes are also involved in the preservation of melanin and melanosomes. In situ polymerization of organic molecules may be promoted via the incorporation of S (in the form of sulfides or polysulfides) from the surrounding environment (de Graaf et al. 1992;Schouten et al. 1994;Sinninghe Damst e et al. 1998a, b, 2007Kok et al. 2000a, b;Werne et al. 2000;van Dongen et al. 2003a, b). This process creates resistant macromolecules comprising carbon chains cross-linked by (poly-)sulfide bridges and can be completed during very early diagenesis (i.e. within 10 ka; Kok et al. 2000a;Werne et al. 2000;Sinninghe Damst e et al. 2007). Further, in situ polymerization via sulfurization, also known as natural vulcanization (Tegelaar et al. 1989), can enhance the preservation potential of organic materials on geological timescales (Sinninghe Damst e et al. 1998a;Melendez et al. 2013a, b). It has been invoked to explain the preservation of very labile organic components, such as carbohydrates (van Kaam-Peters et al. 1998;Sinninghe Damst e et al. 1998b;Kok et al. 2000a;van Dongen et al. 2003avan Dongen et al. , b, 2006, as well as soft tissues in fossils as sulfur-rich organic remains (McNamara et al. 2006(McNamara et al. , 2010. Intriguingly, many examples of melanosome-bearing soft tissues in fossils contain organic sulfur compounds (OSCs; Glass et al. 2012Glass et al. , 2013Lindgren et al. 2014) or are associated with elevated levels of sulfur (Barden et al. 2011(Barden et al. , 2014Lindgren et al. 2012;Pinheiro et al. 2012) , although the latter does not preclude preservation of certain S-bearing amino acids such as cysteine (Barden et al. 2011(Barden et al. , 2014 or of S-bearing phaeomelanins (rather than eumelanin, which does not contain S) (Lindgren et al. 2014). It is therefore possible that preservation of melanosomes in many fossils is promoted by diagenetic sulfurization. Indeed, the possibility of diagenetic incorporation of sulfur into melanosomes has been posited (Lindgren et al. 2014), although this has been questioned (Colleary et al. 2015). The latter study of diverse fossil taxa reported clustering of time of flight secondary ion mass spectrometry (ToF-SIMS) chemical data on the basis of taxonomy, rather than depositional setting; these data were interpreted as evidence that diagenetic sulfurization of eumelanin is 'negligible'. However, data derived from ToF-SIMS analysis cannot differentiate unequivocally between preserved phaeomelanin and diagenetically altered eumelanin (Colleary et al. 2015). This uncertainty regarding the nature of S-bearing moieties in fossil melanosomes limits our understanding of melanosome taphonomy and, critically, our ability to identify fossil phaeomelanin, and thus has important implications for inferences of original colour in fossils. Complementary geochemical analyses of both fossil melanosomes and associated sediments are therefore required in addition to ToF-SIMS analyses to determine the origins of S-bearing moieties in fossil melanosomes.
We investigated sulfur-rich, organically preserved melanosome-like microstructures preserved in fossil frogs, and the host sediment, from the Late Miocene Lagerst€ atte of Libros, NE Spain (Fig. 1A). The soft tissues of the frogs include dark brown to black layers of carbonaceous, spherical to ovoid microstructures previously interpreted as fossil bacteria ( Fig. 1A-B; McNamara et al. 2006McNamara et al. , 2009. The high concentration of sulfur (c. 10% S) in the fossil microstructures has been attributed to sulfurization, but this has not yet been confirmed using organic geochemical techniques (McNamara et al. 2006). Geochemical data have provided critical evidence that similar microstructures associated with the soft tissue outlines of other vertebrate fossils represent preserved melanosomes and not fossil bacteria (Glass et al. 2012(Glass et al. , 2013Lindgren et al. 2012Lindgren et al. , 2014. In light of this, we reassessed our initial interpretations of the nature of the Libros microstructures and tested previous hypotheses that sulfurization contributed to their preservation. Here, we couple detailed analysis of the morphology and spatial context of the microstructures with new geochemical data from pyrolysis gas chromatography mass spectrometry (Py-GCMS), fourier transform infrared spectroscopy (FTIR) and ToF-SIMS. Our new results provide strong evidence that the Libros microstructures are not fossil bacteria but melanosomes and that sulfurization played an important role in their preservation.

GEOLOGICAL BACKGROUND
The Libros frogs are hosted within organic-rich laminated mudstones (oil shales) of the Libros Gypsum Unit (Vallesian) which outcrops c. 25 km SE of Teruel in NE Spain. The unit was deposited in deep-water zones of a stratified freshwater lake (Anad on et al.

Fossil specimens
Most specimens (n = 73) are preserved as well-articulated skeletons with soft tissues defined by a combination of pale-toned phosphatized skin and dark-toned layers (one thick layer internal to the phosphatized skin, and one thin layer, external) of densely packed carbonaceous microbodies c. 0.4-1.0 lm wide ( Fig. 1A; McNamara et al. 2006McNamara et al. , 2009. The context, morphology and ultrastructure of the fossil microbodies have been described previously Microbodies located immediately adjacent to the phosphatized skin (i.e. within c. 5 lm) exhibit a more electron-dense margin interpreted as a diagenetic precipitate of calcium phosphate (Fig. 1G), the phosphate ions for which were derived from the decaying skin (which is rich in calcium phosphate in vivo). These microbodies exhibit a nanocrystalline texture on their surfaces (Fig. 1H) and are rich in Ca and P in addition to C.

Electron microscopy
Small (1-2 mm 2 ) samples of the dark layers were dissected using sterile tools and prepared for scanning-and transmission-electron microscopy (SEM and TEM) as in McNamara et al. (2006). Samples for SEM were sputter coated with gold and examined using a Hitachi S-3500N variable pressure SEM at an accelerating voltage of 15 kV. For TEM, unstained and stained ultrathin (80-90 nm thick) microtome sections were placed on Cu grids and examined using a JEOL 2000TEMSCAN operating at 80 kV and using an objective aperture of 10 lm diameter.

Py-GCMS
Samples of the fossil tissues and host sediment were collected and stored in aluminium foil and crushed in a mortar. A 10 lL measure of tetramethylammonium hydroxide (TMAH) solution was added to each sample to allow thermochemolysis. Py-GCMS was performed using a Chemical Data Systems (CDS) 5200 series pyroprobe pyrolysis unit attached to an Agilent 7890A gas chromatograph (GC) fitted with an HP-5MS fused column (J&W Scientific 5% diphenyl-dimethylpolysiloxane; 30 m, 0.25 mm i.d.; 0.25 lm film thickness) and an Agilent 5975C MSD single quadrupole mass spectrometer operated in electron ionization (EI) mode (EI source temperature 230°C, MS quadrupole 150°C and interface 280°C) with helium as a carrier gas (1 mL/min). Samples were pyrolysed in a quartz tube at 600°C for 20 s, transferred to the GC using a pyrolysis transfer line (350°C) and injected onto the GC using a split ratio of 5:1; the injector port temperature was maintained at 350°C. The oven was programmed to heat at 4°C/min from 40°C (held for 4 min) to 300°C (held for 15 min). The mass spectrometric detector (MSD) mass range scanned from was m/z 60-600 at one scan per second and there was a solvent delay of 4 min. Compounds were identified using the National Institute of Standards and Technology (NIST) database and by comparison with spectra from the literature.

FTIR
Samples of the fossil tissues and a natural melanin standard (Sepia officinalis) were compressed onto the surface of individual pellets of indium. IR spectra of the sample were collected using a Perkin Elmer Spectrum Two spectrometer (wavenumber range 500-4000 cm À1 ) in attenuated total reflectance mode with a resolution of 4 cm À1 . Absorption bands were identified with reference to the NIST reference library and published literature.

ToF-SIMS
ToF-SIMS analysis was performed on a BioToF-SIMS instrument (Kore Technology Ltd, Ely, UK) using a 20 keV Au 3 + primary ion beam (Ionoptika Ltd, Chandlers Ford, UK). The samples were mounted in indium for analysis. Samples of fossil soft tissues and associated sediment were cleaved immediately prior to analysis using a sterile scalpel. The sample of melanin standard was prepared as a powder crushed in indium. The 500 pA primary ion beam was pulsed (120 ns) and scanned over the sample and data acquired in negative ion mode. Secondary ions were extracted into the dual-stage reflectron mass analyser using a delayed extraction field of 2500 V/cm, and detected with a post-acceleration voltage of 7 kV. Charge compensation was provided by a 50 eV electron floodgun (100 nA, 100 ls pulses). The instrument was operated with a mass resolution m/Dm c. 1000 and a lateral resolution c. 1 lm. Images comprising 256 9 256 pixels of total ion signal or selected m/z values were acquired from fossil frog samples over an area ranging from 200 9 200 lm 2 to 1000 9 1000 lm 2 . All analyses were conducted under static SIMS conditions with a primary ion fluence <10 12 ion cm À2 to ensure ion beam damage did not contribute to the measured signal. Optical images of the analysed areas were recorded to allow subsequent SEM imaging of the corresponding regions.

Py-GCMS
The pyrogram of the fossil microbodies is dominated by propenoic acid butyl ester (compound 'c' in Fig 2A) and a range of organic sulfur compounds (OSCs), including dimethyl disulfide (compound 'a' in Fig. 2A), methylsulfonic acid methyl ester (compound 'b' in Fig. 2A) and a series of short chain alkylated thiophenes ( Fig. 2A; Table 1). Also present are an extended series of fatty acid methyl esters that range up to C 19 and have maxima at C 9 and/or C 16 , and derivatives of benzene, phenol and benzoic acid. Notably, a series of n-alkane-/ alkene doublets and bacterially-derived biomarkers, such as hopanoid pyrolysis products, are absent. In contrast, analyses of the host sediment produced a pyrogram dominated by extended series of n-alkane-/alkene doublets and fatty acid methyl esters that range up to C 29 and have a maximum at C 24 ( Fig. 2B; Table 1). Hopanes and OSCs, predominantly short chain alkylated thiophenes, are also present, but are relatively minor components.

FTIR
The measured FTIR spectra reveal a close similarity between the preserved organic functional groups present in the fossil microbodies and the melanin standard (Fig. 3, Table 2). The FTIR spectrum of the latter is dom-inated by three major absorption bands (Hong & Simon 2006;Centeno et al. 2008;Glass et al. 2012). The broad peak centred at 3200 cm À1 corresponds to absorption due to the stretching mode of the OH bond. The relatively sharp bands at 1606 cm À1 and 1290 cm À1 correspond to the carbonyl C=O stretching vibration, probably in indole quinone, and a combination of OH and NH inplane bending vibrations, respectively. The minor band at 775 cm À1 may reflect bending modes of aromatic C-H or C=C bonds.
These absorption bands also dominate FTIR spectra of the fossil microbodies ( Fig. 3; Table 2). In addition, the fossil spectra also show a major band at 1152 cm À1 and minor bands at 2920 cm À1 , 1713 cm À1 , 1428 cm À1 , 1040 cm À1 and 598 cm À1 . Several of these bands can be attributed to the presence of sulfur: the major band at 1152 cm À1 corresponds to the C=S stretching vibration and the minor bands at 1040 cm À1 and 598 cm À1 may represent stretching vibrations in CS functional groups (Bloxham et al. 2002); CO stretching may also contribute to the band at 1040 cm À1 . The minor band at 2920 cm À1 is due to stretching vibrations in CH x functional groups, probably lipids (Glass et al. 2012); this is consistent with the minor band at 1428 cm À1 , which reflects the bending mode of CH bonds in CH 2 or CH 3 groups. The band at 1713 cm À1 is attributed to the carbonyl ketone stretch and may indicate lower levels of indole quinone in the fossil sample than in the melanin standard (Glass et al.

ToF-SIMS
The samples analysed using ToF-SIMS comprise exclusively ovoid microstructures ( Fig. 4A-C). There is a very close agreement between the spectra from the fossil microbodies and the melanin standard, in terms of both the position and relative intensity of peaks (Fig. 4D). Key peaks in the spectrum for the melanin standard at m/z 50, 66, 73, 74, 97, 98, 121 and 122 reflect the nitrogenrich chemical structure of melanin (Lindgren et al. 2012). These are also among the most prominent peaks in the fossil microbodies (Fig. 4D). Additional peaks at m/z 80 and m/z 97 in the fossil sample correspond to SO 3 À and HSO 4 À and also occur in the host sediment (Fig. 4E).
F I G . 3 . FTIR absorption spectra for the fossil microbodies in the Libros frogs and in a melanin standard (Sepia officinalis). Details of the absorption bands and associated symbols are in Table 2.
ToF-SIMS spectra from the sediment and fossil tissues are notably dissimilar; none of the key peaks in the melanin standard are prominent in the sediment (Fig. 4D

DISCUSSION
The nature of the fossil microbodies The microbodies associated with the soft tissues of the Libros frogs were originally interpreted as fossilized bacteria on the basis of their size, geometry, restriction to the soft tissues of the specimens, and the presence of conjoined microbodies suggestive of cell mitosis ( In addition to these data on anatomical and sedimentary context, the microbodies from the Libros amphibians exhibit key ultrastructural features that support their being melanosomes. In unstained and stained TEM sections, the fossil microbodies are uniformly electron-dense (Fig. 1E); electron density is enhanced, but remains uniform among microbodies, in stained sections (Fig. 1F) is present in some of the microbodies from the Libros amphibians (Fig. 1G). This feature is only evident, however, in microbodies located immediately adjacent to (i.e. within c. 5 lm of) the phosphatized skin. The electrondense margin of these microbodies has been interpreted as a diagenetic precipitate of calcium phosphate, with the T A B L E 2 . Infrared peak assignments for the fossil microbodies and a melanin standard (Sepia officinalis). Fig. 3 Bond Wavenumber (cm À1 ) Mode of vibration Melanin standard Fossil microbodies

Notation in
Numbers and symbols refer to Our new chemical data provide strong evidence that the preserved microbodies in the Libros frogs are not fossil bacteria. Neither Py-GCMS nor ToF-SIMS analyses reveal evidence of bacterially derived compounds. This is not a taphonomic artefact: bacterially-derived hopane moieties are present in the host sediment. Hopanoid producing microbes are predominantly aerobic; abundant hopane moieties in other fossil vertebrate soft tissues has been interpreted as evidence of extended decay in aerobic showing densely packed, spheroidal to ovoid microbodies. D, negative ion ToF-SIMS spectra for the region of the fossil tissue sample highlighted in yellow in F, and for a synthetic melanin standard; filled circles indicate ions used for ToF-SIMS images in G and I. E, negative ion ToF-SIMS spectrum for a sample of sediment from the area indicated by * in Fig. 1A. F, H, ToF-SIMS images of total ion counts; F, sample shown in A; H, a selected region of the counterpart of the sample. G, I, ion images showing the distribution of ions derived from melanin (blue: m/z 50; red: m/z 66) and organosulfur compounds (yellow: m/z 80; green: m/z 97) superimposed onto scanning electron micrographs; G, is the sample in A; I, shows a detailed view of densely packed microbodies from the area indicated by a box in H; the signal is strongest in the lower part of the sample as this region was orientated towards the detector. Scale bars represent: 500 lm (A); 300 lm (B); 4 lm (C); 5 lm (H-I).
conditions (Edwards et al. 2011). The lack of hopane moieties in the Libros frog tissues could indicate predominantly anaerobic decay before and after deposition. The presence of hopanes in the host sediment could reflect slower rates of settling of organic particles through the water column and thus more prolonged decay in aerobic waters. Further, FTIR and ToF-SIMS analyses reveal a strong spectral agreement between the fossil tissues and the melanin standard (Figs 3, 4D); ToF-SIMS data also demonstrate that key peaks for melanin are localized to the fossil microbodies ( Fig. 4F-I). These data confirm that the microbodies contain high concentrations of melanin and/or its derivatives. Both Py-GCMS and ToF-SIMS analyses reveal striking differences between the chemical composition of the fossil microbodies in the Libros frogs and that of the sedimentary matrix. In particular, the microbodies and sediment differ markedly in the abundance of n-alkane-/alkene doublets and in the fatty acid methyl ester distribution patterns and maxima evident in the pyrograms, and in peak position and intensity in the ToF-SIMS spectra. These data strongly indicate that the melanin-based, and other, organic compounds in the microbodies are indigenous and are not derived from the host sediment. Collectively, the context, ultrastructure and chemical composition of the fossil microbodies in the Libros amphibians are strong evidence that the fossil microbodies in the Libros frogs are fossil melanosomes.

Preservation through sulfurization (natural vulcanization)
Py-GCMS and ToF-SIMS analyses confirm the presence of diverse OSCs in the Libros melanosomes and the sedimentary matrix; this is supported by FTIR spectra, which indicate the presence of various sulfur moieties not present in the eumelanin standard. These OSCs are unlikely to reflect the presence of S-bearing phaeomelanins as typical pyrolysis products of phaeomelanin, for example, S-bearing compounds such as benzothiazine and benzothiazole (Lindgren et al. 2014) were not detected. Further, ToF-SIMS spectra for phaeomelanin show prominent peaks for S-bearing molecular fragments at m/z 57, 58, 81, 82, 105, 106, 118 and 134 (Lindgren et al. 2014). With the exception of the peak at m/z 81 (which we have attributed to HSO 3 ), none of these peaks are prominent in our fossil samples. Small amounts of sulfur can occur within natural eumelanins, but their origins are unclear (Jimbow et al. 1984;Prota 1992) and they are unlikely to be responsible for the high abundance and diversity of OSCs in the Libros melanosomes. Many of the OSCs in the Libros melanosomes are also present in the host sediment, although in lower abundance, where a melanin origin is highly improbable. Thus the OSCs in the Libros melanosomes most likely indicate the presence of sulfurized organic matter. This confirms previous suggestions (McNamara et al. 2006) that the carbonaceous soft tissues in the Libros amphibians are preserved via sulfurization. The presence of abundant sulfurized organic matter in the sedimentary matrix is consistent with previous reports of OSCs in the sedimentary kerogen (del R ıo et al. 2004) and confirms that sulfurization of organic matter was widespread in the Libros sediments. The extended series T A B L E 3 . Tentative assignments and m/z ratios for peaks in negative ToF-SIMS spectra of melanin and microstructures in the Libros frogs. of fatty acid methyl esters also present in the fossil soft tissues and sedimentary matrix indicates that part of the organic matter is preserved via polymerization processes that did not involve the incorporation of sulfur. S-mediated reduction processes are critical to the preservation of organic matter in sediments (Hebting et al. 2006). Sulfurization of organic matter is characteristic of environments where inputs of bacterially generated sulfide exceed the availability of reactive iron (Sinninghe Damst e et al. 1998a). Low concentrations of reactive iron are, however, not required for sulfurization: certain organic molecules are predisposed to sulfurization and can be sulfurized very early during diagenesis, even within the water column (Sinninghe Damst e et al. 1998a;Melendez et al. 2013a, b). Sulfurization may therefore occur under a broad spectrum of depositional conditions; unsurprisingly, OSCs have been reported from various lacustrine and marine settings from the Permian to the Recent (Sinninghe Damst e et al. 1993, 1998a, 1999, 2007Grice et al. 1996Grice et al. , 1998van Kaam-Peters et al. 1998;Kok et al. 2000a, b;Werne et al. 2000;Kolonic et al. 2002;van Dongen et al. 2006;Jaraula et al. 2013). Sulfurization of organic molecules in the sedimentary kerogen and fossil melanosomes from Libros is attributed to deposition beneath anoxic, sulfidic bottom waters poor in reactive iron (del R ıo et al. 2004;McNamara et al. 2006). Such conditions are likely to characterize other fossil deposits in which melanosomes are preserved. Indeed, the presence of elevated levels of sulfur and, especially, of OSCs in fossil melanosomes associated with other vertebrate and invertebrate fossils (see above) strongly suggests preservation via sulfurization. These fossils are hosted within diverse sediments, including freshwater lacustrine varves (Xiagou Formation, China), volcanic maar oil shales (Oligocene of Enspel, Germany), marine diatomite (Fur Formation, Denmark) and marine oil shales (Cretaceous Boquillas Formation; Jurassic of Dorset, UK). Geochemical analyses of the host sediments at these and other localities offer a test for the possibility that sulfurization occurred and thus constrain the confidence with which potential phaeomelanosomes can be identified in fossils.
Melanin is a large and relatively inert biopolymer (Riley 1997). Other biopolymers, e.g. cellulose, lignin and chitin, are not usually the primary contributors to the preservation of organic soft tissues in fossils as the process of in situ polymerization involves primarily endogenous lipids (Gupta et al. 2009;Gupta & Briggs 2011). Sulfurization, however, is known to affect biomolecules of different recalcitrance, including carbohydrates, lipids, and structural biopolymers, including melanin (Jimbow et al. 1984;Melendez et al. 2013a, b). In pyrograms of the Libros melanosomes, the high abundance and distribution pattern of short-chain alkylated (C 0 -C 5 ) thiophenes resemble those of sulfurized monosaccharides (carbohydrates; Sinninghe Damst e et al. 1998b;Kok et al. 2000b;van Dongen et al. 2003a) and kerogen from the Kimmeridge Clay Formation (van Kaam-Peters & Sinninghe Damst e 1997, van Kaam-Peters et al. 1998van Dongen et al. 2006). Diagenetic sulfurization is known to generate a spectrum of organic compounds (i.e. low molecular weight OSCs and macromolecular S-bearing aggregates) that are resistant to bacterial degradation and remineralization and are thus preserved in the fossil record (Valisolalao et al. 1984;Brassell et al. 1986;Sinninghe Damst e & de Leeuw 1990;Sinninghe Damst e et al. 1998a, b). Preservation via sulfurization is therefore likely to enhance the inherent resistance (Goldstein et al. 2004) of melanins to microbial degradation.

CONCLUSIONS
Although the anatomical and sedimentological context of the microbodies associated with the soft tissues of the Libros frogs is not conclusive, their geochemical composition (dominated by melanin and distinct from the host sediment) confirms that the fossil microbodies represent fossilized melanosomes. The high abundance of OSCs in the fossil melanosomes confirms that the latter are preserved via diagenetic incorporation of sulfur, i.e. preservation through sulfurization. This process is likely to enhance the inherent recalcitrance and thus preservation potential of melanin. Given that sulfurization can occur in various marine and freshwater depositional settings, and that many examples of fossil melanosomes are associated with sulfur and/or organosulfur compounds, this mode of preservation is likely to contribute widely to preservation of melanosomes in the fossil record. Future studies of sulfur-rich fossil melanosomes will require the host sediment to be tested for geochemical evidence of sulfurization in order to constrain interpretations of potential phaeomelanosomes and thus of original integumentary colour in fossils.