**Data archiving statement**

Data for this study are available in the Dryad Digital Repository: <http://dx.doi.org/10.5061/dryad.540q6>

The 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*. [2008](#pala12238-bib-0010){ref-type="ref"}; Li *et al*. [1984](#pala12238-bib-0033){ref-type="ref"}, [1997](#pala12238-bib-0034){ref-type="ref"}; Zhang *et al*. [2010](#pala12238-bib-0077){ref-type="ref"}; Carney *et al*. [2012](#pala12238-bib-0009){ref-type="ref"}; Lindgren *et al*. [2011](#pala12238-bib-0036){ref-type="ref"}), the evolutionary history of visual signalling (Clarke *et al*. [2008](#pala12238-bib-0010){ref-type="ref"}; Li *et al*. [1984](#pala12238-bib-0033){ref-type="ref"}, [1997](#pala12238-bib-0034){ref-type="ref"}; Zhang *et al*. [2010](#pala12238-bib-0077){ref-type="ref"}; Carney *et al*. [2012](#pala12238-bib-0009){ref-type="ref"}) and of physiology (Li *et al*. [1997](#pala12238-bib-0034){ref-type="ref"}) and the fidelity of the fossil record (Glass *et al*. [2014](#pala12238-bib-0018){ref-type="ref"}, [2012](#pala12238-bib-0019){ref-type="ref"}; Lindgren *et al*. [2011](#pala12238-bib-0036){ref-type="ref"}). 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*. [1989](#pala12238-bib-0073){ref-type="ref"}; Barden *et al*. [2011](#pala12238-bib-0003){ref-type="ref"}; Knight *et al*. [2011](#pala12238-bib-0029){ref-type="ref"}; Carney *et al*. [2012](#pala12238-bib-0009){ref-type="ref"}), birds (Clarke *et al*. [2008](#pala12238-bib-0010){ref-type="ref"}; Zhang *et al*. [2010](#pala12238-bib-0077){ref-type="ref"}; Wogelius *et al*. [2008](#pala12238-bib-0076){ref-type="ref"}), dinosaurs (Clarke *et al*. [2008](#pala12238-bib-0010){ref-type="ref"}; Li *et al*. [1984](#pala12238-bib-0033){ref-type="ref"}, [1997](#pala12238-bib-0034){ref-type="ref"}; Zhang *et al*. [2010](#pala12238-bib-0077){ref-type="ref"}), fish (Lindgren *et al*. [1998](#pala12238-bib-0035){ref-type="ref"}), marine reptiles (Lindgren *et al*. [2011](#pala12238-bib-0036){ref-type="ref"}), mammals (Li *et al*. [1997](#pala12238-bib-0034){ref-type="ref"}) and squid (Glass *et al*. [2014](#pala12238-bib-0018){ref-type="ref"}, [2012](#pala12238-bib-0019){ref-type="ref"}). These fossil microstructures are superficially similar to some bacteria in terms of their morphology and size (Vinther *et al*. [2010](#pala12238-bib-0077){ref-type="ref"}; Zhang *et al*. [2010](#pala12238-bib-0077){ref-type="ref"}; Knight *et al*. [2011](#pala12238-bib-0029){ref-type="ref"}) 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](#pala12238-bib-0004){ref-type="ref"}; Moyer *et al*. [2010](#pala12238-bib-0050){ref-type="ref"}). Certain trace elements (e.g. Cu, Zn, Ca, and S) have been proposed as biomarkers for melanin (Wogelius *et al*. [2008](#pala12238-bib-0076){ref-type="ref"}; Manning *et al*. [2000b](#pala12238-bib-0038){ref-type="ref"}). Fossil melanosomes can also retain organic geochemical evidence for the melanin molecule itself (Glass *et al*. [2014](#pala12238-bib-0018){ref-type="ref"}, [2012](#pala12238-bib-0019){ref-type="ref"}; Lindgren *et al*. [1998](#pala12238-bib-0035){ref-type="ref"}, [2011](#pala12238-bib-0036){ref-type="ref"}). 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](#pala12238-bib-0045){ref-type="ref"}) and other taxa (Colleary *et al*. [2010](#pala12238-bib-0011){ref-type="ref"}) 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*. [1998](#pala12238-bib-0035){ref-type="ref"}). Analyses of melanosomes preserved in Jurassic squid indicate diagenetic cross‐linking (Glass *et al*. [2014](#pala12238-bib-0018){ref-type="ref"}). Previous studies have referred to the high recalcitrance of melanin (Glass *et al*. [2014](#pala12238-bib-0018){ref-type="ref"}): it is a highly cross‐linked inert polymer of dihydroxyindole and dihydroxyindole carboxylic acid (McGraw [2010](#pala12238-bib-0040){ref-type="ref"}) that resists microbial degradation (Goldstein *et al*. [2013](#pala12238-bib-0020){ref-type="ref"}), oxidation and hydrolysis by acids and bases (Riley [2013b](#pala12238-bib-0055){ref-type="ref"}). 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](#pala12238-bib-0008){ref-type="ref"}; Stankiewicz *et al*. [2012](#pala12238-bib-0064){ref-type="ref"}; Gupta *et al*. [1996](#pala12238-bib-0024){ref-type="ref"}, [2011](#pala12238-bib-0025){ref-type="ref"}; Manning *et al*. [2000a](#pala12238-bib-0037){ref-type="ref"}, [2002](#pala12238-bib-0039){ref-type="ref"}; Dutta *et al*. [2006](#pala12238-bib-0015){ref-type="ref"}; Gupta & Briggs [1998](#pala12238-bib-0023){ref-type="ref"}; Barden *et al*. [2011](#pala12238-bib-0003){ref-type="ref"}; Edwards *et al*. [2011](#pala12238-bib-0017){ref-type="ref"}). This process generates an aliphatic component similar to Type I and II kerogens (Briggs [1999](#pala12238-bib-0008){ref-type="ref"}) 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](#pala12238-bib-0003){ref-type="ref"}, [2014](#pala12238-bib-0004){ref-type="ref"}; Glass *et al*. [2012](#pala12238-bib-0019){ref-type="ref"}), 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*. [2015](#pala12238-bib-0012){ref-type="ref"}; Schouten *et al*. [2014](#pala12238-bib-0057){ref-type="ref"}; Sinninghe Damsté *et al*. [1992](#pala12238-bib-0060){ref-type="ref"}, [*b*](#pala12238-bib-0061){ref-type="ref"}, [2004](#pala12238-bib-0063){ref-type="ref"}; Kok *et al*. [2003](#pala12238-bib-0030){ref-type="ref"}, [*b*](#pala12238-bib-0031){ref-type="ref"}; Werne *et al*. [1984](#pala12238-bib-0075){ref-type="ref"}; van Dongen *et al*. [1998a](#pala12238-bib-0068){ref-type="ref"}, [*b*](#pala12238-bib-0069){ref-type="ref"}). 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*. [2003](#pala12238-bib-0030){ref-type="ref"}; Werne *et al*. [1984](#pala12238-bib-0075){ref-type="ref"}; Sinninghe Damsté *et al*. [2004](#pala12238-bib-0063){ref-type="ref"}). Further, *in situ* polymerization via sulfurization, also known as natural vulcanization (Tegelaar *et al*. [1994](#pala12238-bib-0065){ref-type="ref"}), can enhance the preservation potential of organic materials on geological timescales (Sinninghe Damsté *et al*. [1992](#pala12238-bib-0060){ref-type="ref"}; Melendez *et al*. [2006](#pala12238-bib-0047){ref-type="ref"}, [*b*](#pala12238-bib-0048){ref-type="ref"}). It has been invoked to explain the preservation of very labile organic components, such as carbohydrates (van Kaam‐Peters *et al*. [2000](#pala12238-bib-0072){ref-type="ref"}; Sinninghe Damsté *et al*. [2006](#pala12238-bib-0061){ref-type="ref"}; Kok *et al*. [2003](#pala12238-bib-0030){ref-type="ref"}; van Dongen *et al*. [1998a](#pala12238-bib-0068){ref-type="ref"}, [*b*](#pala12238-bib-0069){ref-type="ref"}, [1999](#pala12238-bib-0070){ref-type="ref"}), as well as soft tissues in fossils as sulfur‐rich organic remains (McNamara *et al*. [2012](#pala12238-bib-0041){ref-type="ref"}, [2014](#pala12238-bib-0043){ref-type="ref"}). Intriguingly, many examples of melanosome‐bearing soft tissues in fossils contain organic sulfur compounds (OSCs; Glass *et al*. [2014](#pala12238-bib-0018){ref-type="ref"}, [2012](#pala12238-bib-0019){ref-type="ref"}; Lindgren *et al*. [2011](#pala12238-bib-0036){ref-type="ref"}) or are associated with elevated levels of sulfur (Barden *et al*. [2011](#pala12238-bib-0003){ref-type="ref"}, [2014](#pala12238-bib-0004){ref-type="ref"}; Lindgren *et al*. [1998](#pala12238-bib-0035){ref-type="ref"}; Pinheiro *et al*. [2013](#pala12238-bib-0052){ref-type="ref"}) , although the latter does not preclude preservation of certain S‐bearing amino acids such as cysteine (Barden *et al*. [2011](#pala12238-bib-0003){ref-type="ref"}, [2014](#pala12238-bib-0004){ref-type="ref"}) or of S‐bearing phaeomelanins (rather than eumelanin, which does not contain S) (Lindgren *et al*. [2011](#pala12238-bib-0036){ref-type="ref"}). 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*. [2011](#pala12238-bib-0036){ref-type="ref"}), although this has been questioned (Colleary *et al*. [2010](#pala12238-bib-0011){ref-type="ref"}). 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*. [2010](#pala12238-bib-0011){ref-type="ref"}). 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ätte of Libros, NE Spain (Fig. [1](#pala12238-fig-0001){ref-type="fig"}A). The soft tissues of the frogs include dark brown to black layers of carbonaceous, spherical to ovoid microstructures previously interpreted as fossil bacteria (Fig. [1](#pala12238-fig-0001){ref-type="fig"}A--B; McNamara *et al*. [2012](#pala12238-bib-0041){ref-type="ref"}, [2012](#pala12238-bib-0042){ref-type="ref"}). 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*. [2012](#pala12238-bib-0041){ref-type="ref"}). 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*. [2014](#pala12238-bib-0018){ref-type="ref"}, [2012](#pala12238-bib-0019){ref-type="ref"}; Lindgren *et al*. [1998](#pala12238-bib-0035){ref-type="ref"}, [2011](#pala12238-bib-0036){ref-type="ref"}). 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.

![Melanosome‐like microbodies in the Libros frogs. A, MNCN 63663; inset, detail of thin, dark brown carbonaceous film defining soft tissues in area indicated; asterisk indicates region of sediment analysed in Fig. [4](#pala12238-fig-0004){ref-type="fig"}E. B--H, scanning (B--C, H) and transmission (D--G) electron micrographs showing details of melanosome‐like microbodies in the brown layer. B--C, densely packed melanosome‐like microbodies (B) with detail of surface texture (C). D--F, unstained (D--E) and stained (F) TEM sections of microbodies showing uniform electron density; note internal vacuoles in microbodies in D. G--H, unstained sections of microbodies immediately adjacent to the phosphatized skin showing electron‐dense margin of calcium phosphate (G) and nanocrystalline surface texture (H). Scale bars represent: 50 mm (A); 1 mm (inset in A); 5 μm (B, F); 500 nm C, E--F); 1 μm (D); 2 μm (G).](PALA-59-337-g001){#pala12238-fig-0001}

Geological background {#pala12238-sec-0002}
=====================

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ón *et al*. [1992](#pala12238-bib-0002){ref-type="ref"}); low oxygen isotope values from primary gypsum deposits and abundant gammacerane in the laminated mudstones indicate persistent high sulfide levels and intense bacterial sulfate reduction in the monimolimnion (de las Heras *et al*. [2003a](#pala12238-bib-0013){ref-type="ref"}; Ortí *et al*. [2012](#pala12238-bib-0051){ref-type="ref"}). The kerogen from the laminated mudstones contains abundant OSCs indicative of extensive sulfurization (del Río *et al*. [2003b](#pala12238-bib-0014){ref-type="ref"}) and there is no evidence for the presence of benthic microbial mats (McNamara *et al*. [2009](#pala12238-bib-0044){ref-type="ref"}).

Material and method {#pala12238-sec-0003}
===================

Fossil specimens {#pala12238-sec-0004}
----------------

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 μm wide (Fig. [1](#pala12238-fig-0001){ref-type="fig"}A; McNamara *et al*. [2012](#pala12238-bib-0041){ref-type="ref"}, [2012](#pala12238-bib-0042){ref-type="ref"}). The context, morphology and ultrastructure of the fossil microbodies have been described previously (McNamara *et al*. [2012](#pala12238-bib-0042){ref-type="ref"}). In brief, the microbodies are restricted to the outline of the soft tissues and most are located internal to the phosphatized skin where they envelop the bones and stomach contents (McNamara *et al*. [2012](#pala12238-bib-0042){ref-type="ref"}). The microbodies typically exhibit granular to botryoidal surface nanotextures (Fig. [1](#pala12238-fig-0001){ref-type="fig"}B--C) and are usually uniformly electron‐dense in unstained (Fig. [1](#pala12238-fig-0001){ref-type="fig"}D--E) and stained (Fig. [1](#pala12238-fig-0001){ref-type="fig"}F) TEM sections; electron‐lucent inclusions may be apparent (Fig. [1](#pala12238-fig-0001){ref-type="fig"}D--E). Microbodies located immediately adjacent to the phosphatized skin (i.e. within *c*. 5 μm) exhibit a more electron‐dense margin interpreted as a diagenetic precipitate of calcium phosphate (Fig. [1](#pala12238-fig-0001){ref-type="fig"}G), 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. [1](#pala12238-fig-0001){ref-type="fig"}H) and are rich in Ca and P in addition to C.

###  {#pala12238-sec-0005}

#### Institutional abbreviation {#pala12238-sec-0006}

MNCN, Museo Nacional de Ciencias Naturales, Madrid.

Electron microscopy {#pala12238-sec-0007}
-------------------

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*. ([2012](#pala12238-bib-0041){ref-type="ref"}). 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 μm diameter.

Py‐GCMS {#pala12238-sec-0008}
-------

Samples of the fossil tissues and host sediment were collected and stored in aluminium foil and crushed in a mortar. A 10 μL 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 μm 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 {#pala12238-sec-0009}
----

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 {#pala12238-sec-0010}
--------

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 μs pulses). The instrument was operated with a mass resolution *m*/Δ*m c*. 1000 and a lateral resolution *c*. 1 μm. Images comprising 256 × 256 pixels of total ion signal or selected *m*/*z* values were acquired from fossil frog samples over an area ranging from 200 × 200 μm^2^ to 1000 × 1000 μm^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.

Results {#pala12238-sec-0011}
=======

Py‐GCMS {#pala12238-sec-0012}
-------

The pyrogram of the fossil microbodies is dominated by propenoic acid butyl ester (compound 'c' in Fig [2](#pala12238-fig-0002){ref-type="fig"}A) and a range of organic sulfur compounds (OSCs), including dimethyl disulfide (compound 'a' in Fig. [2](#pala12238-fig-0002){ref-type="fig"}A), methylsulfonic acid methyl ester (compound 'b' in Fig. [2](#pala12238-fig-0002){ref-type="fig"}A) and a series of short chain alkylated thiophenes (Fig. [2](#pala12238-fig-0002){ref-type="fig"}A; Table [1](#pala12238-tbl-0001){ref-type="table-wrap"}). 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. [2](#pala12238-fig-0002){ref-type="fig"}B; Table [1](#pala12238-tbl-0001){ref-type="table-wrap"}). Hopanes and OSCs, predominantly short chain alkylated thiophenes, are also present, but are relatively minor components.

![Partial TMAH‐assisted Py‐GCMS total ion current chromatograms of Libros microstructures (A) and matrix (B); insets (not to scale) show the *m/z* 97 + 98 + 111 + 112 + 125 + 126 + 139 + 140 + 153 + 154 chromatograms of dotted areas revealing the presence and distribution of alkylated thiophene moieties. \*, alkylated thiophene; +, fatty acid methyl ester; --, *n*‐alkane/*n*‐alkene doublets; \#, benzene derivative; ○, phenol derivative; ▲, benzoic acid derivative; ![](PALA-59-337-g005.jpg "image"), hopane moieties; a, dimethyl disulfide; b, methylsulfonic acid methyl ester; c, propenoic acid butyl ester; d, succinic acid dimethylester; e, methyl succinic acid dimethyl ester; f, methyl indolinone; g, phthalate; h, 2,5‐dimethoxy benzoic acid methyl ester (or isomer); i, contaminant; j, methyl phenylsulfide; k, 3,5‐dimethyl‐2‐(methylsulfanyl) thiophene. c~x~ indicates the carbon chain length and numbers refer to the alkylated thiophenes listed in Table [1](#pala12238-tbl-0001){ref-type="table-wrap"}.](PALA-59-337-g002){#pala12238-fig-0002}

###### 

Alkylated thiophenes identified in the TMAH‐assisted pyrolysates of Libros microstructures and sedimentary matrix

  ---- -------------------------------------------------------
  1    2‐methylthiophene
  2    3‐methylthiophene
  3    2‐ethylthiophene
  4    2,5‐dimethylthiophene
  5    2,4‐dimethylthiophene
  6    2,3‐dimethylthiophene
  7    3,4‐dimethylthiophene
  8    2‐propylthiophene
  9    2‐ethyl‐5‐methylthiophene
  10   2‐ethyl‐4‐methylthiophene
  11   Ethylmethylthiophene
  12   2,3,5‐trimethylthiophene
  13   2,3,4‐trimethylthiophene
  14   C~8~H~12~S[†](#pala12238-note-0002){ref-type="fn"}
  15   C~6~H~6~O~2~S[‡](#pala12238-note-0002){ref-type="fn"}
  16   C~7~H~8~O~2~S[§](#pala12238-note-0002){ref-type="fn"}
  ---- -------------------------------------------------------

Numbers refer to Fig. [2](#pala12238-fig-0002){ref-type="fig"}; †Mixture of 2‐methyl‐5‐propylthiophene, 2,5‐diethylthiophene, 2‐butylthiophene, 2‐ethyl‐3,5‐dimethylthiophene, ethyldimethylthiophene and/or 5‐ethyl‐2,3‐dimethylthiophene; ‡Mixture of methyl‐2‐thiophene carboxylate and methyl‐3‐thiophene carboxylate; §Mixture of methyl‐5‐methyl‐2‐thiophene carboxylate and methyl‐3‐methyl‐2‐thiophene carboxylate.

John Wiley & Sons, Ltd

FTIR {#pala12238-sec-0013}
----

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](#pala12238-fig-0003){ref-type="fig"}, Table [2](#pala12238-tbl-0002){ref-type="table-wrap"}). The FTIR spectrum of the latter is dominated by three major absorption bands (Hong & Simon [2006](#pala12238-bib-0031){ref-type="ref"}; Centeno *et al*. [2008](#pala12238-bib-0010){ref-type="ref"}; Glass *et al*. [2014](#pala12238-bib-0018){ref-type="ref"}). 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 in‐plane bending vibrations, respectively. The minor band at 775 cm^−1^ may reflect bending modes of aromatic C‐H or C=C bonds.

![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](#pala12238-tbl-0002){ref-type="table-wrap"}.](PALA-59-337-g003){#pala12238-fig-0003}

###### 

Infrared peak assignments for the fossil microbodies and a melanin standard (*Sepia officinalis*)

  Notation in Fig. [3](#pala12238-fig-0003){ref-type="fig"}   Bond              Wavenumber (cm^−1^)   Mode of vibration   Melanin standard   Fossil microbodies
  ----------------------------------------------------------- ----------------- --------------------- ------------------- ------------------ --------------------
  1                                                           O‐H               3200--3207            Stretch             X                  X
  2                                                           C=O               1606--1620            Stretch             X                  X
  3                                                           O‐H ± N‐H         1290--1315            Bend                X                  X
  4                                                           C‐H (aromatic)    775--780              Bend                X                  X
  \*                                                          C‐H (aliphatic)   2920                  Stretch             --                 X
  \*\*                                                        C=O (ketone)      1713                  Stretch             --                 X
  †                                                           C‐H (aliphatic)   1428                  Bend                --                 X
  ‡                                                           C‐S ± C‐O         1040                  Stretch             --                 X
  ◦                                                           C=S               1152                  Stretch             --                 X
  ●                                                           C‐S               598                   --                  --                 X

Numbers and symbols refer to Fig. [3](#pala12238-fig-0003){ref-type="fig"}.

John Wiley & Sons, Ltd

These absorption bands also dominate FTIR spectra of the fossil microbodies (Fig. [3](#pala12238-fig-0003){ref-type="fig"}; Table [2](#pala12238-tbl-0002){ref-type="table-wrap"}). 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](#pala12238-bib-0006){ref-type="ref"}); 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*. [2014](#pala12238-bib-0018){ref-type="ref"}); 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*. [2014](#pala12238-bib-0018){ref-type="ref"}).

ToF‐SIMS {#pala12238-sec-0014}
--------

The samples analysed using ToF‐SIMS comprise exclusively ovoid microstructures (Fig. [4](#pala12238-fig-0004){ref-type="fig"}A--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. [4](#pala12238-fig-0004){ref-type="fig"}D). Key peaks in the spectrum for the melanin standard at *m/z* 50, 66, 73, 74, 97, 98, 121 and 122 reflect the nitrogen‐rich chemical structure of melanin (Lindgren *et al*. [1998](#pala12238-bib-0035){ref-type="ref"}). These are also among the most prominent peaks in the fossil microbodies (Fig. [4](#pala12238-fig-0004){ref-type="fig"}D). 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. [4](#pala12238-fig-0004){ref-type="fig"}E). Almost all other major peaks in the spectrum from the fossil soft tissues can be tentatively assigned to fragment ions consistent with the structure of nitrogen‐rich melanins (Table [3](#pala12238-tbl-0003){ref-type="table-wrap"}). Further, ion images showing the spatial signal intensity distribution for key melanin peaks (*m/z* 50 and 66) superimposed on SEM images demonstrate a clear correlation between the presence of the fossil microbodies and the peaks for melanin (Fig. [4](#pala12238-fig-0004){ref-type="fig"}F--I). With the exception of the peaks at *m/z* 80 and 97 (see above), 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. [4](#pala12238-fig-0004){ref-type="fig"}D--E). The dominant peaks in the sediment spectrum correspond to carbon chains of various length, some possibly bearing O or N; proteinaceous moieties have been previously identified in the organic matter from the Libros oil shales (del Río *et al*. [2003b](#pala12238-bib-0014){ref-type="ref"}).

![Melanin in Miocene frogs from Libros. A--C, light micrograph (A) and scanning electron micrographs (B--C) of freshly fractured sample of soft tissues from MNCN 63663. B, detail of area indicated by box in A. C, detail of area indicated by \* in B, 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. [1](#pala12238-fig-0001){ref-type="fig"}A. 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 μm (A); 300 μm (B); 4 μm (C); 5 μm (H--I).](PALA-59-337-g004){#pala12238-fig-0004}

###### 

Tentative assignments and *m/z* ratios for peaks in negative ToF‐SIMS spectra of melanin and microstructures in the Libros frogs

  Tentative assignment     Theoretical mass (u)   Melanin standard   Fossil frog   Sediment
  ------------------------ ---------------------- ------------------ ------------- ----------
  C~4~                     48.00                  48.00              48.00         47.98
  C~4~H                    49.01                  49.01              49.02         49.00
  C~3~N                    50.00                  50.00              50.03         50.00
  C~3~HN                   51.01                  51.01              51.03         51.02
  C~3~H~2~N                52.02                  52.02              52.02         52.01
  C~5~                     60.00                  60.00              60.01         59.99
  C~5~H                    61.01                  61.01              61.00         61.00
  C~4~N/C~5~H~2~           62.00/62.02            62.01              62.02         62.00
  C~4~HN                   63.01                  63.01              63.02         63.00
  C~4~H~2~N/CH~4~O~3~      64.02/64.02            64.01              63.98         63.96
  C~4~HO/H~3~NO~3~         65.00/65.01            65.01              65.02         65.00
  C~3~NO                   66.00                  66.00              66.02         65.99
  C~6~                     72.00                  72.00              72.00         72.00
  C~6~H/C~2~H~3~NO~2~      73.01/73.02            73.01              73.02         73.01
  C~5~N/C~2~H~2~O~3~       74.00/74.00            74.00              74.02         74.01
  SO~3~                    79.96                  79.98              79.97         79.96
  HSO~3~                   80.97                  80.99              80.98         80.98
  C~7~                     84.00                  84.00              84.00         84.00
  C~7~H                    85.01                  85.01              85.01         85.02
  C~6~N/C~7~H~2~           86.00/86.02            86.01              86.01         86.00
  C~6~HN                   87.01                  87.01              87.02         87.02
  C~6~H~2~N/C~3~H~4~O~3~   88.02/88.02            88.02              88.02         87.99
  C~6~HO/C~2~H~3~NO~3~     89.00/89.01            89.01              89.02         88.99
  C~5~NO                   90.00                  90.00              90.01         90.00
  C~7~H~7~                 91.06                  91.02              91.02         91.01
  C~8~                     96.00                  96.02              95.96         95.97
  HSO~4~/C~4~H~3~NO~2~     96.96/97.02            97.01              96.97         96.98
  C~7~N/C~4~H~2~O~3~       98.00/98.00            98.00              98.00         98.00
  C~9~                     108.00                 108.00             108.00        108.00
  C~9~H                    109.01                 109.01             109.00        108.99
  C~8~N/C~9~H~2~           110.00/110.02          110.01             110.01        110.01
  C~8~HN                   111.01                 111.01             111.02        111.02
  C~8~H~2~N/C~5~H~4~O~3~   112.02/112.02          112.02             112.00        111.98
  C~8~HO/C~4~H~3~NO~3~     113.00/113.01          113.02             112.99        112.99
  C~7~NO                   114.00                 114.02             114.01        114.00
  C~10~                    120.00                 120.00             119.99        112.00
  C~10~H/C~6~H~3~NO~2~     121.01/121.02          121.01             121.00        120.98
  C~9~N                    122.00                 122.00             121.99        121.98
  C~12~H/C~8~H~3~NO~2~     145.01/145.02          145.01             144.99        144.96

Grey shading highlights key N‐bearing molecular ions likely to be derived from melanin.

John Wiley & Sons, Ltd

Discussion {#pala12238-sec-0015}
==========

The nature of the fossil microbodies {#pala12238-sec-0016}
------------------------------------

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 (McNamara *et al*. [2012](#pala12238-bib-0042){ref-type="ref"}). As others have noted (Zhang *et al*. [2010](#pala12238-bib-0077){ref-type="ref"}), the size and geometry of the microbodies from the fossil amphibians cannot be considered to be definitive. The structures range between *c*. 0.4 μm and 1.0 μm in width (Fig. [1](#pala12238-fig-0001){ref-type="fig"}B; McNamara *et al*. [2012](#pala12238-bib-0042){ref-type="ref"}); this falls within the range of structures identified as melanosomes in other fossil vertebrates (Li *et al*. [1984](#pala12238-bib-0033){ref-type="ref"}; Zhang *et al*. [2010](#pala12238-bib-0077){ref-type="ref"}; Lindgren *et al*. [1998](#pala12238-bib-0035){ref-type="ref"}) but alone does not exclude their being bacteria. The fossil microbodies envelop the bones and the stomach contents (McNamara *et al*. [2012](#pala12238-bib-0042){ref-type="ref"}), i.e. they occur on both ventral and dorsal surfaces of the body, and most are located internal to the layer of phosphatized skin. This distribution is consistent with their being melanosomes (Zhang *et al*. [2010](#pala12238-bib-0077){ref-type="ref"}), but not an external 'microbial overgrowth' *sensu* Moyer *et al*. ([2010](#pala12238-bib-0050){ref-type="ref"}). Further, there is no evidence in the host sediments for benthic microbial mats (McNamara *et al*. [2009](#pala12238-bib-0044){ref-type="ref"}). The seemingly connected microbodies resembling bacterial cells in the process of fission could simply represent closely juxtaposed melanosomes. Indeed, the dense packing of microbodies in the Libros amphibians and of melanosomes reported from other fossil taxa (Clarke *et al*. [2008](#pala12238-bib-0010){ref-type="ref"}; Vinther *et al*. [2010](#pala12238-bib-0077){ref-type="ref"}; Zhang *et al*. [2010](#pala12238-bib-0077){ref-type="ref"}; Barden *et al*. [2011](#pala12238-bib-0003){ref-type="ref"}; Glass *et al*. [2014](#pala12238-bib-0018){ref-type="ref"}; Lindgren *et al*. [1998](#pala12238-bib-0035){ref-type="ref"}, [2011](#pala12238-bib-0036){ref-type="ref"}) is notable. This could originate via degradation and collapse of more labile soft tissues that enveloped the melanosomes *in vivo*, thus forming a layer of melanosomes that does not retain the original spatial arrangement and density of melanosomes in the tissue (and may even derive from different tissues; M.E. McNamara, J.S. Kaye, M.J. Benton, P.J. Orr & N.P. Lockyer, unpub. data). Notably in the case of the adult frogs from Libros, the microbodies external to the phosphatized skin show a markedly different distribution on the dorsal and ventral surfaces. Microbodies are typically densely packed on the dorsal surface of the specimen, but sparse on the ventral surface. This distribution is consistent with the distribution of integumentary melanosomes in the integument of many extant frogs, which generates countershading (Beebee & Griffiths [2000](#pala12238-bib-0005){ref-type="ref"}).

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. [1](#pala12238-fig-0001){ref-type="fig"}E); electron density is enhanced, but remains uniform among microbodies, in stained sections (Fig. [1](#pala12238-fig-0001){ref-type="fig"}F). In stained TEM sections, melanosomes in extant vertebrates are also uniformly electron‐dense (Prum [2013a](#pala12238-bib-0054){ref-type="ref"}; Hellström *et al*. [2006](#pala12238-bib-0026){ref-type="ref"}; Moyer *et al*. [2010](#pala12238-bib-0050){ref-type="ref"}), but bacteria exhibit a marginal, electron‐dense cell wall surrounding a relatively electron‐lucent cell lumen (Mesnage *et al*. [2009](#pala12238-bib-0049){ref-type="ref"}; Saar‐Dover *et al*. [1998](#pala12238-bib-0056){ref-type="ref"}). A similar ultrastructure (i.e. an electron‐dense margin surrounding an electron‐lucent core) is present in some of the microbodies from the Libros amphibians (Fig. [1](#pala12238-fig-0001){ref-type="fig"}G). This feature is only evident, however, in microbodies located immediately adjacent to (i.e. within *c*. 5 μm of) the phosphatized skin. The electron‐dense margin of these microbodies has been interpreted as a diagenetic precipitate of calcium phosphate, with the phosphate ions derived from the decaying skin (which is rich in calcium phosphate *in vivo*; McNamara *et al*. [2012](#pala12238-bib-0042){ref-type="ref"}). These microbodies exhibit a nanocrystalline texture on their surfaces (Fig. [1](#pala12238-fig-0001){ref-type="fig"}H) (McNamara *et al*. [2012](#pala12238-bib-0042){ref-type="ref"}) and are rich in Ca and P, in contrast to the other, uniformly electron‐dense, microbodies, which are carbonaceous in composition (McNamara *et al*. [2012](#pala12238-bib-0041){ref-type="ref"}). Many of the microbodies from the Libros amphibians exhibit electron‐lucent inclusions in TEM images (Fig. [1](#pala12238-fig-0001){ref-type="fig"}E) that are strikingly similar to the internal vacuoles of melanosomes in extant vertebrate tissues (Akazaki *et al*. [2014](#pala12238-bib-0001){ref-type="ref"}; Moyer *et al*. [2010](#pala12238-bib-0050){ref-type="ref"}), but are rare in bacteria. Further, the carbonaceous fossil microbodies from Libros typically exhibit subtle granular to botryoidal surface nanotextures (Fig. [1](#pala12238-fig-0001){ref-type="fig"}C) that resemble those of modern melanosomes (Moyer *et al*. [2010](#pala12238-bib-0050){ref-type="ref"}).

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 conditions (Edwards *et al*. [2010](#pala12238-bib-0016){ref-type="ref"}). 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](#pala12238-fig-0003){ref-type="fig"}, [4](#pala12238-fig-0004){ref-type="fig"}D); ToF‐SIMS data also demonstrate that key peaks for melanin are localized to the fossil microbodies (Fig. [4](#pala12238-fig-0004){ref-type="fig"}F--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) {#pala12238-sec-0017}
----------------------------------------------------------

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*. [2011](#pala12238-bib-0036){ref-type="ref"}) 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*. [2011](#pala12238-bib-0036){ref-type="ref"}). 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*. [2006](#pala12238-bib-0028){ref-type="ref"}; Prota [1992](#pala12238-bib-0060){ref-type="ref"}) 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*. [2012](#pala12238-bib-0041){ref-type="ref"}) 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*. [2003b](#pala12238-bib-0014){ref-type="ref"}) and confirms that sulfurization of organic matter was widespread in the Libros sediments. The extended series 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](#pala12238-bib-0028){ref-type="ref"}). Sulfurization of organic matter is characteristic of environments where inputs of bacterially generated sulfide exceed the availability of reactive iron (Sinninghe Damsté *et al*. [1992](#pala12238-bib-0060){ref-type="ref"}). 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é *et al*. [1992](#pala12238-bib-0060){ref-type="ref"}; Melendez *et al*. [2006](#pala12238-bib-0047){ref-type="ref"}, [*b*](#pala12238-bib-0048){ref-type="ref"}). 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é *et al*. [2012](#pala12238-bib-0059){ref-type="ref"}, [1992](#pala12238-bib-0060){ref-type="ref"}, [1997](#pala12238-bib-0062){ref-type="ref"}, [2004](#pala12238-bib-0063){ref-type="ref"}; Grice *et al*. [1992](#pala12238-bib-0022){ref-type="ref"}, [2004](#pala12238-bib-0021){ref-type="ref"}; van Kaam‐Peters *et al*. [2000](#pala12238-bib-0072){ref-type="ref"}; Kok *et al*. [2003](#pala12238-bib-0030){ref-type="ref"}, [*b*](#pala12238-bib-0031){ref-type="ref"}; Werne *et al*. [1984](#pala12238-bib-0075){ref-type="ref"}; Kolonic *et al*. [2013](#pala12238-bib-0032){ref-type="ref"}; van Dongen *et al*. [1999](#pala12238-bib-0070){ref-type="ref"}; Jaraula *et al*. [2009](#pala12238-bib-0027){ref-type="ref"}). 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*. [2003b](#pala12238-bib-0014){ref-type="ref"}; McNamara *et al*. [2012](#pala12238-bib-0041){ref-type="ref"}). 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 [2013b](#pala12238-bib-0055){ref-type="ref"}). 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*. [2011](#pala12238-bib-0025){ref-type="ref"}; Gupta & Briggs [1998](#pala12238-bib-0023){ref-type="ref"}). Sulfurization, however, is known to affect biomolecules of different recalcitrance, including carbohydrates, lipids, and structural biopolymers, including melanin (Jimbow *et al*. [2006](#pala12238-bib-0028){ref-type="ref"}; Melendez *et al*. [2006](#pala12238-bib-0047){ref-type="ref"}, [*b*](#pala12238-bib-0048){ref-type="ref"}). 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é *et al*. [2006](#pala12238-bib-0061){ref-type="ref"}; Kok *et al*. [2006](#pala12238-bib-0031){ref-type="ref"}; van Dongen *et al*. [1998a](#pala12238-bib-0068){ref-type="ref"}) and kerogen from the Kimmeridge Clay Formation (van Kaam‐Peters & Sinninghe Damsté [2007](#pala12238-bib-0071){ref-type="ref"}, van Kaam‐Peters *et al*. [2000](#pala12238-bib-0072){ref-type="ref"}; van Dongen *et al*. [1999](#pala12238-bib-0070){ref-type="ref"}). 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*. [1993](#pala12238-bib-0067){ref-type="ref"}; Brassell *et al*. [1986](#pala12238-bib-0007){ref-type="ref"}; Sinninghe Damsté & de Leeuw [2003](#pala12238-bib-0058){ref-type="ref"}; Sinninghe Damsté *et al*. [1992](#pala12238-bib-0060){ref-type="ref"}, [*b*](#pala12238-bib-0061){ref-type="ref"}). Preservation via sulfurization is therefore likely to enhance the inherent resistance (Goldstein *et al*. [2013](#pala12238-bib-0020){ref-type="ref"}) of melanins to microbial degradation.

Conclusions {#pala12238-sec-0018}
===========

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.

We thank Begoña Sanchez and Patricia Perez for access to specimens and sampling permission, Stuart Kearns for assistance with SEM analyses, Matthew James for assistance with Py‐GCMS analyses and Cormac O\'Connell and Dave Cottell for assistance with TEM analyses. The research was funded by Enterprise Ireland Basic Research Grant SC/2002/138 to PJO, a Marie Curie Career Integration Grant awarded to MMN and NERC Grant NE/J023426/1 to BEvD.

MMN designed the research and undertook SEM and TEM analyses. BVD and IDB undertook GC‐MS and Py‐GCMS analyses. NPL undertook ToF‐SIMS analyses. MMN wrote the paper with input from all other authors.
