Molecular fossils within bitumens and kerogens from the ~ 1 Ga Lakhanda Lagerstätte (Siberia, Russia) and their significance for understanding early eukaryote evolution

The emergence and diversification of eukaryotes during the Proterozoic is one of the most fundamental evolutionary developments in Earth’s history. The ca. 1-billion-year-old Lakhanda Lagerstätte (Siberia, Russia) contains a wealth of eukaryotic body fossils and offers an important glimpse into their ecosystem. Seeking to complement the paleontological record of this remarkable lagerstätte, we here explored information encoded within sedimentary organic matter (total organic carbon = 0.01–1.27 wt.%). Major emphasis was placed on sedimentary hydrocarbons preserved within bitumens and kerogens, including molecular fossils (or organic biomarkers) that are specific to bacteria and eukaryotes (i.e. hopanes and regular steranes, respectively). Programmed pyrolysis and molecular organic geochemistry suggest that the organic matter in the analyzed samples is about peak oil window maturity and thus sufficiently well preserved for detailed molecular fossil studies that include hopanes and steranes. Together with petrographic evidence as well as compositional similarities of the bitumens and corresponding kerogens, the consistency of different independent maturity parameters establishes that sedimentary hydrocarbons are indigenous and syngenetic to the host rock. The possible presence of trace amounts of hopanes and absence of steranes in samples that are sufficiently well preserved to retain both types of compounds evidences an environment dominated by anaerobic bacteria with no or very little inputs by eukaryotes. In concert with the paleontological record of the Lakhanda Lagerstätte, our study adds to the view that eukaryotes were present but not significant in Mesoproterozoic ecosystems.


Introduction
The Proterozoic is a critical interval for the early evolution of eukaryotes on Earth. It was during this time that unicellular eukaryotic clades established, sexual reproduction and eukaryophagy evolved, and complex multicellularity developed (e.g., Butterfield 2000Butterfield , 2009Butterfield , 2015Javaux 2007Javaux , 2011Knoll 2011Knoll , 2014Porter 2011Porter , 2004Javaux and Lepot 2018). Not surprisingly, the Proterozoic evolution of eukaryotes has been intensely studied from different perspectives and by using different methodologies. An exhaustive and detailed account on all these endeavors goes far beyond the scope of this paper, and only some of the most important paleontological benchmarks are briefly outlined in the following and illustrated in Fig. 1a.
The oldest undisputed microfossils of eukaryotic origin date back to ca. 1.7-1.5 billion years (in the following 1 3 provided as Giga-annum, Ga) (Fig. 1a). Some of the earliest evidence are acritarchs from the ~ 1.74-1.41 Ga Ruyang Group in North China (Pang et al. 2013;Hu et al. 2014;Lan et al. 2014;Agić et al. 2017). The acritarch-bearing Changzhougou Formation in China was previously dated at 1.80 Ga, which would make it older than the Ruyang Group, but was later dated to ~ 1.67-1.64 Ga (e.g., Lamb et al. 2009;Peng et al. 2009;Miao et al. 2019). This age range is fairly similar to that of other notable acritarch records such as the ~ 1.65 Ga old Mallapunyah Formation in Australia (Javaux 2007) and the ≤ 1.63 Ga old Vindhyan Supergroup in India (Ray et al. 2002;Prasad et al. 2005;Singh and Sharma 2014;Sharma et al. 2016). A further important microfossil record of early eukaryote diversity is the ~ 1.50 Ga old Roper Group in Australia (Javaux et al. 2001;Javaux and Knoll 2017). The 1.0-0.9 Ga old Grassy Bay Formation in Arctic Canada contains the organic-walled microfossil Ourasphaira giraldae, which is interpreted as fungus (Loron et al. 2019a, b).
Possible macrofossil evidence for the presence of eukaryotes comprises Grypania, a carbonaceous fossil that assuredly occurs worldwide between 1.6 and 1.4 Ga (Walter et al. 1990;Kumar 1995) (Fig. 1a) and perhaps even as early as ~ 1.9 Ga (Han and Runnegar 1992;see Schneider et al. 2002 for the age constraint). Early modestly diverse assemblages of carbonaceous compression fossils of macroscopic, multicellular eukaryotes are known from the 1.56 Ga old Gaoyuzhuang Formation in China (Zhu et al. 2016). A particularly important benchmark of early eukaryote evolution is the ~ 1.05 Ga old fossil Bangiomorpha pubescens (previously dated at 1.20 Ga), which is widely accepted Fig. 1 Introduction into early eukaryote evolution and molecular fossils. a Some important benchmarks in early eukaryotic evolution as chronicled by the fossil record (see main text for details and references); b Examples of molecular fossils that are indicative for bacteria (hopanes, left) and eukaryotes (regular steranes, right). R = H or an alkyl group. Such molecular fossils can be identified with gas chromatography-mass spectrometry (GC-MS) and gas chromatogra-phy-tandem mass spectrometry (GC-MS/MS). Particularly important are their mass spectrometric fragmentation patterns: Hopanes show a characteristic main fragment at m/z 191 (corresponding to the A/B ring fragment), while regular steranes exhibit a major fragment at m/z 217 (corresponding to the A/B/C ring fragment). Stratigraphic boundary ages in a from Ogg et al. (2016) as representing a bangiophyte red algae and the first solid evidence of complex multicellularity and sexual reproduction (Butterfield et al. 1990;Butterfield 2000;Gibson et al. 2018). Furthermore, a ~ 1 Ga old benthic chlorophyte (Proterocladus antiquus) has recently been reported from the Nanfen Formation in China (Tang et al. 2020).
Clearly, the micro and macrofossil records provide valuable insights into early eukaryote evolution (Fig. 1a). However, these classic paleontological records are by no means complete since most organisms are not preserved as fossils. Hydrocarbons in ancient sediments can encode information of geobiological value and provide evidence for the presence of organism groups that are not likely to be preserved as fossil (e.g., soft-bodied organisms, microbes). Particularly interesting are molecular fossils (or organic biomarkers), that is, hydrocarbon derivatives of organic biomolecules that retained characteristic structural features and are, therefore, still diagnostic for specific source organisms (e.g., Treibs 1934Treibs , 1936Eglinton 1964;Peters et al. 2005a, b;Briggs and Summons 2014;Schwarzbauer and Jovančićević 2016;Falk and Wolkenstein 2017). Since all types of organisms consist of organic biomolecules, molecular fossils provide an additional, independent source of information that can complement the classic paleontological record. The most established molecular fossils of bacteria and eukaryotes are hopanes and regular (that is, 4-desmethyl) steranes, respectively (Fig. 1b). Both of these are sedimentary hydrocarbons that retained the polycyclic nucleus of the respective precursor biomolecules and can be preserved through geologic time.
As outlined above, eukaryotic body fossils reach back into the Paleoproterozoic (Fig. 1a). In rocks older than ca. 820 million years, however, regular steranes are typically close to, or below, the detection limit (Dutkiewicz et al. 2003;Brocks et al. 2005;Blumenberg et al. 2012;Pawlowska et al. 2013;Flannery and George, 2014;Luo et al. 2015;Suslova et al. 2017;Nguyen et al. 2019;Zhang et al. 2021). This mismatch might reflect a minor relative importance of eukaryotes in early Precambrian ecosystems and/or non-actualistic taphonomic processes in the ancient environments (see Knoll et al. 2007;Knoll 2014;Love and Zumberge, 2021 for detailed reviews). Sedimentary hydrocarbon studies on rocks that contain early eukaryotic body fossils potentially provide important clues that can help to clarify this issue.
The ca. 1.00-1.03 Ga old Lakhanda Lagerstätte in Siberia, Russia (Fig. 2) provides an important window into an early eukaryotic ecosystem. Over many decades of research, the deposits have yielded a wealth of remarkably well-preserved eukaryotic fossils (e.g., Aimonema, Caudosphaera, Jacutianema, Palaeovaucheria, Valeria, Trachyhystrichosphaera, and Germinosphaera) (Jankauskas et al. 1989;Hermann 1990;Rainbird et al. 1998;Semikhatov et al. 2000;Ovchinnikova et al. 2001;Hermann and Podkovyrov 2006, 2010Nagovitsin 2008; German [Hermann] and Podkovyrov 2009;Podkovyrov 2009;Shuvalova et al. 2021b). Numerous new findings, such as potential evidence for eukaryophagy (Shuvalova et al. 2021a), highlight that the archive has not been fully exploited yet. This appears also to be true for the molecular fossil record of this remarkable Lagerstätte, which has been considered to be representative for the Proterozoic and might hold great promise for detailed geobiological studies (Pawlowska et al. 2013).
As mentioned before, molecular fossil studies are a powerful complement to paleontological investigations. However, ancient organic matter records can be compromised by destructive processes as well (e.g., biodegradation, thermal maturation during burial), and might additionally be adulterated by an introduction of non-indigenous hydrocarbons (e.g., through migration and/or contamination) (Sherman et al. 2007;Eigenbrode 2008;Brocks 2011;Duda et al. 2014b;Love and Zumberge 2021). For these reasons, paleoreconstructions that involve molecular fossils must be commenced by a rigorous assessment of organic matter preservation (biodegradation, thermal maturity) and hydrocarbon syngenicity-regardless of the geologic age of the samples. The results must then be reported and documented for scrutiny in any publication arising from the research (Love and Zumberge 2021). So far, however, organic matter maturity and hydrocarbon syngenicity have not been appraised and/ or documented for the Lakhanda Lagerstätte.
Here we present the first systematic assessment of organic matter preservation in the Lakhanda Lagerstätte, involving the detailed analysis of a substantial sample set from bulk to molecular levels. Major emphasis was placed on sedimentary hydrocarbon inventories in the extractable bitumen and-for the first time-non-extractable kerogen phases, including molecular fossils that are specific to eukaryotes (i.e. regular steranes: Fig. 1b). Our integrative approach involved the combined application of various analytical techniques, including programmed temperature pyrolysis, catalytic hydropyrolysis (HyPy), gas chromatography-mass spectrometry (GC-MS) and gas chromatography-tandem mass spectrometry (GC-MS/MS). The study is the first to establish a robust framework for research endeavors on the Lakhanda Lagerstätte that involve organic geochemical analyses. Our findings have important consequences for the taphonomic and paleoecological interpretation of sedimentary hydrocarbons in the Lakhanda Lagerstätte and Proterozoic successions, including molecular fossils.

Sampling and petrography
This study includes the organic-geochemical analysis of 13 samples from the Neryuen Formation of the Lakhanda Group. 12 samples were taken from the Kumakha Subformation at Ytyrynda (Figs. 1, 3a). 1 additional sample (M-01, 1935) was taken from the Nel'kan Subformation at Olemeken (Figs. 1, 3b). It should be noted that the stratigraphic hierarchy has not been firmly established, and that the Kumakha and Nel'kan subformations have been ranked as formations in some recent studies (e.g., Nagovitsin et al. 2008;Shuvalova et al. 2021a, b). Both sections are about 25 km apart (Fig. 2).
Thin sections (vertical cross-sections) were prepared for all samples. Petrographic observation of thin sections (transmitted and reflected light) was conducted with a Keyence VHX-7000 Digital Microscope.

Bulk analyses
Contents of total organic carbon (TOC), sulfur (TS), and nitrogen (TN) were determined with a Hekatech Euro EA elemental analyzer at the Geobiology Department of the University of Göttingen, Germany. TOC was additionally measured with different temperature programmable carbon analyzer instruments in independent labs (Leco RC612: Geobiology Department, University of Göttingen, Germany; Leco SC-632: Applied Petroleum Technology, AS, Oslo, Norway). All values are reported in weight % (wt.%) of bulk sedimentary rock.
Programmed temperature pyrolysis was performed by Applied Petroleum Technology (Oslo, Norway) with a HAWK instrument. The temperature was ramped from 300 °C (held for 3 min) to 650 °C at 25 °C/min. The measurements included quantities of free hydrocarbons (S1, mg HC/g rock) and hydrocarbons yielded from labile kerogen (S2, mg HC/g rock), temperatures at maximum yields of S2 hydrocarbons (T max , °C), as well as the CO 2 generated from organic carbon at higher temperatures up to 650 °C (S3, mg CO 2 /g rock). These data, as well as TOC data obtained by the same lab, were used to calculate the hydrogen index (HI; S2/TOC*100), the oxygen index (OI; S3/TOC*100) and the production index (PI; S1/(S1 + S2)) for all samples.

Analysis of sedimentary hydrocarbons in bitumen and kerogen fractions
Samples were analyzed for sedimentary hydrocarbons using an established methodology (Duda et al. 2014a(Duda et al. , b, 2016(Duda et al. , 2020Reinhardt et al. 2018). All laboratory materials used were heated to 500 °C for 3 h and/or extensively rinsed with acetone. Only distilled solvent was used for extraction and further work.
Exterior parts of the samples were removed with a precleaned precision saw (Buehler IsoMet1000). Interior and exterior parts were then crushed and powdered using a pebble mill (Retsch MM 301). Ground sample material (20-80 g) was sequentially extracted with (i) dichloromethane (DCM), (ii) DCM/n-hexane (1/1; v/v) and n-hexane (20 min ultrasonification, respectively). The obtained total organic extracts (TOEs) were combined, gently concentrated with a cleaned rotary evaporator and N 2 to avoid any loss of low-boiling compounds (Ahmed and George 2004), and desulfurized with activated copper. The TOEs were then fractionated by column chromatography (see below), while the extraction residues (i.e. inorganic matter and macromolecular, non-extractable organic matter) were used for catalytic hydropyrolysis (HyPy).
Molecular fossil studies are usually based on the extractable bitumen fraction. Problematically, however, bitumen is mobile and vulnerable to contamination on geological and anthropogenic scales. Such processes can adversely affect sedimentary hydrocarbon records preserved in ancient rocks. The macromolecular kerogen, in contrast, is immobile and less prone to contamination. Hydrocarbons bound into the kerogen structure are therefore likely syngenetic to the host rock and provide a critical self-consistency check for testing the syngenicity of compounds in the corresponding bitumen (Love et al. 1995(Love et al. , 2008(Love et al. , 2009. Catalytic hydropyrolysis (HyPy) is a state-of-the-art technique for releasing kerogen-bound hydrocarbons with minimal change to structure and stereochemistry, thereby preserving their most diagnostic features (Love et al. 1995(Love et al. , 2008Meredith et al. 2014). The yielded products can then be analyzed by means of gas chromatography-mass spectrometry (GC-MS) and gas chromatography-tandem mass spectrometry (GC-MS/ MS), similarly to the bitumen portions that are conventionally analyzed (see below).
HyPy was conducted following existing protocols (Love et al. 1995Meredith et al. 2014) using an instrument build by Strata Technology Ltd. (Nottingham, UK). Prior to sample runs, the system was heated from ambient temperature up to 500 °C (held for 30 min) to eliminate potential contaminants. For the sample runs, about 500 mg extraction residue was pyrolyzed in the presence of ~ 26 mg sulfided molybdenum catalyst (ammonium dioxydithiomolybdate, [(NH 4 ) 2 MoO 2 S 2 ]) and pre-combusted sea sand. The extraction residues were progressively heated from ambient temperature to 250 °C at 50 °C/min and then to 500 °C at 8 °C/min under 15 MPa pressure and a constant H 2 flow rate of 5 L/min. The HyPy products were trapped downstream on silica powder cooled with dry ice (Meredith et al. 2004). The HyPy pyrolysates were subsequently eluted from the silica powder with DCM, desulfurized overnight with activated copper, and then fractionated by column chromatography (see below).

Gas chromatography-(tandem) mass spectrometry (GC-MS, GC-MS/MS)
The F1 and F2 of all samples were analyzed with GC-MS and highly sensitive GC-MS/MS. The measurements were conducted with a Thermo Scientific Trace 1300 Series GC coupled to a Thermo Scientific Quantum XLS Ultra MS. The instrument was equipped with a capillary column (Phenomenex Zebron ZB-5, 30 m, 0.25 μm film thickness, inner diameter 0.25 mm). The fractions were injected into a spitless injector and transferred to the GC column at 300 °C with helium as carrier gas and a flow rate of 1.5 ml/min. The GC oven temperature was ramped from 80 °C (1 min) to 310 °C with 5 °C/min and held for 20 min. GC-MS was conducted in full scan (mass range of m/z 50-600, scan time = 0.42 s) or GC-MS/MS mode (parent-daughter ion mass transitions for steranes and hopanes). The ion source in the MS was operated in electron impact mode at 250 °C and with 70 eV ionization energy. Argon was used as the collision gas. The extractable aliphatic fraction of a thermally immature black shale of known molecular composition (Posidonia Shale, Lower Toarcian, Holzmaden, southern Germany) containing abundant steranes and hopanes was analyzed in parallel to facilitate the identification of these compounds in the Lakhanda samples.

Petrography and bulk data
The analyzed samples mainly consist of wavy laminated, dark gray to brown argillites that contain abundant pyrite (Fig. 3c, d). The samples do not show any cracks or fissures, and there is no petrographic indication for the migration of bitumen.

Sedimentary hydrocarbons in bitumens and kerogens
Chromatograms of all bitumen and kerogen portions do not exhibit a major hump of an unresolved complex mixture of organic compounds (UCM) and are characterized by similar molecular inventories (Fig. 4). n-Alkanes typically range from C 12 to C 37 and exhibit unimodal distributions with maxima between n-C 15 and n-C 17 . They lack odd or even carbon number preferences as reflected in CPI 24-34 values of ~ 1. Methyl-branched alkanes mostly range from methyl-C 13 to methyl-C 30 and follow the distribution trend of the corresponding n-alkanes (Figs. 4, 5). Phytane is only present in minute traces, while pristane was not found to be present. Eukaryotic steranes and their degradation products (i.e., short-chain steranes) were not detected with GC-MS and GC-MS/MS (Figs. 4, 6), whilst some samples might contain trace amounts of bacterial hopanes (Fig. 7). The reason for the absence of 22,29,30-trisnorhopanes (Ts and Tm) in 370 → 191 GC-MS/MS ion transitions is unclear, but given the overall low abundance of hopanes Ts and Tm might be obscured due to low signal-to-noise ratios.
The aromatic fractions of bitumens and kerogens were only analyzed for the possible presence of aromatic steroids Table 2 Bulk data and programmed temperature pyrolysis parameters Note that some of the data is not considered robust and therefore discounted. Total organic carbon (TOC) contents were measured with different instruments in independent labs (Lab I: Geobiology Department, University of Göttingen, Germany; Lab II: Applied Petroleum Technology, AS, Oslo, Norway; one asterisks: Hekatech Euro EA elemental analyzer; two asterisks: Leco carbon analyzers, see main text for details). TOC and total contents of nitrogen (TN) and sulfur (TS) are all reported in weight % of bulk sedimentary rock. S1 and S2 values are provided in mg HC/g rock, S3 values are given in mg CO 2 /g rock. The hydrogen index (HI) was calculated as S2/TOC*100, the oxygen index (OI) was calculated as S3/TOC*100, and the production index (PI) was calculated as S1/(S1 + S2) (not detected) and maturity-sensitive polyaromatic hydrocarbons (that is, phenanthrene and methylphenanthrenes: Fig. 8).

Data integrity
T max values (440-445 °C, mean = 442 °C: Table 2) correspond to the early oil window, approaching the peak oil generation stage (Peters 1986;Peters and Cassa 1994). The inferred maturity is consistent with PI values (0.06-0.17, mean 0.10) and also in good agreement with low HI and OI values (65-103 mg HC/g TOC and 16-73 mg CO 2 /g TOC, respectively: Table 2), which indicate a relatively low to moderate hydrocarbon generation potential. An analytical bias due to the adsorption of pyrolytic products onto clay minerals-as commonly observed for argillaceous rocks with TOC contents below ca. 0.5 wt.% (Peters 1986)-can likely be excluded since there is no obvious systematic relationship between organic matter contents and programmed temperature pyrolysis parameters. At the same time, sedimentary organic matter is a complex mixture of biochemicals deriving from various organisms, and thus source effects cannot be excluded (see Peters 1986;Peters and Cassa 1994). Taken together, programmed temperature pyrolysis data indicate that the organic matter in both successions has been thermally altered to an extent that principally still allows for the preservation of biologically diagnostic molecular fossils. Thermal maturity of sedimentary organic matter can further be scrutinized on a molecular basis, including both, bitumen and kerogen-based molecular indices. MPI-1 and R C values of bitumen proportions (0.90-1.05 and 0.85-0.95, respectively: Table 1) indicate a peak oil generation stage (Radke and Welte 1983;Radke et al. 1986). This is consistent with MPR values from bitumens (1.31-1.74: Table 1), which in tendency point to the later part of the peak oil window (Radke et al. 1984). Although the MPI-1 was originally  (Radke et al. 1982), values of ≥ 0.7 might also be robust to assess the thermal maturity of Proterozoic marine organic matter (Boreham et al. 1988). Furthermore, we only used aromatic hydrocarbons from the extractable bitumen fractions to evaluate the thermal maturity here, as kerogen-based MPI-1 and MPR ratios are known to be of limited applicability for thermal maturity assessments (Marshall et al. 2007;Pehr et al. 2021). Indeed, the inferred peak-oil window maturity is in good accordance with calculated CPIs (mainly about 1.10) and phytane/n-C 18 and pristane/n-C 17 ratios (<< 1) in the bitumens and pyrolysates (Peters et al. 2005b;ten Haven et al. 1988) (Table 3). At his thermal maturity, molecular fossils such as hopanes and steranes should still be preserved.
At first glance, there seems to be a slight discrepancy between programmed temperature pyrolysis data (early oil window, approaching the peak oil generation stage) and molecular maturity parameters (at or just beyond peak oil generation stage). However, it has been found that programmed temperature pyrolysis data such as T max tend to underestimate actual maturity levels of thermally altered Proterozoic organic matter (Zumberge et al. 2020). Taking this into consideration, the samples analyzed herein are probably rather at the upper end of the maturity range inferred from programmed temperature pyrolysis data (that is, about peak oil window maturity). The data are thus consistent with independent molecular maturity parameters, pointing against an admixture of younger hydrocarbons to the organic pool of the successions. This is further supported by the strikingly similar molecular compositions of bitumens and the corresponding kerogens (Figs. 4,5), as well as by similar hydrocarbon distributions in extracts from interior and exterior portions of all samples. Thus, our results establish that sedimentary hydrocarbons in the Lakhanda Lagerstätte are indigenous and syngenetic to the host rocka crucial prerequisite for geobiological studies that involve molecular fossils such as hopanes and regular steranes.

Meaning for molecular fossil studies on the Lakhanda Lagerstätte
Many of the observed molecular characteristics are quite similar to other mid-Proterozoic shales such as the ~ 1.4 Ga Velkerri Formation in Australia (Dutkiewicz et al. 2003;Flannery and George 2014), the ~ 1.4 Ga Xiamaling Formation in China (Luo et al. 2015;Zhang et al. 2021) and  (Blumenberg et al. 2012). These features include, amongst others, (i) a dominance of low-molecular-weight n-alkanes with maxima at ≤ C 19 , (ii) abundant monomethylalkanes, (iii) low abundances of acyclic isoprenoids such as pristane and phytane (except for Tourist Formation: Blumenberg et al. 2012), and (iv) low contents or even an absence of eukaryotic steranes. It is therefore tempting to speculate that the molecular fossil record of the Lakhanda Lagerstätte might be representative for mid-Proterozoic ecosystems, as proposed earlier based on one sample from the Neryuen Formation (Pawlowska et al. 2013, their "bitumen facies 1").
Fossil lagerstätten are not necessarily representative as they formed under rather unusual conditions (Seilacher 1970). If the molecular fossil record of the Lakhanda Lagerstätte was nevertheless representative for the  (22,29,Ts and Tm), C 29 and C 30 . A thermally immature black shale of known molecular composition (Lower Toarcian, Holzmaden, southern Germany; HM) is provided as reference to verify the retention times of these compounds. In the sample from the Lakhanda Lagerstätte, candidate hopane peaks cooccur with multiple background peaks, which perhaps also explains the seeming absence of Ts and Tm. The chromatograms are scaled to the most abundant peak in each trace mid-Proterozoic, it would allow to reconstruct baseline conditions in early eukaryotic ecosystems. Unfortunately, however, our assessment of organic matter preservation reveals that many of the molecular characteristics observed in the Lakhanda Lagerstätte are due to thermal maturity (about peak oil window maturity). For instance, methyl-branched alkanes (Figs. 4,5) might indeed be sourced by various bacteria (Gelpi et al. 1970;Shiea et al. 1990;Kaneda 1991), but can also form through thermal alteration processes during catagenesis (Kissin 1987). Biosynthetic pathways are highly selective and therefore typically result in particular homologue or isomer preferences. The lack of such preferences and the similar distribution trends of methyl-branched alkanes and n-alkanes in all analyzed samples herein (Figs. 4,5) thus point to a catagenetic origin, which is well in line with independent programmed temperature pyrolysis and molecular maturity parameters.
Another example are acyclic isoprenoids such as pristane and phytane. The detected phytane (Figs. 4,5) likely derived from photoautotrophic bacteria and/or eukaryotes, although archaeal sources cannot be excluded (Peters et al. 2005b). The problem here is that acyclic isoprenoids are known to be considerably affected by thermal degradation (e.g., Peters et al. 2005b;Mißbach et al. 2016). Considering the results from the systematic and rigorous preservation assessment in this study, the low concentrations of phytane in the Neryuen Formation (Figs. 4,5) are indeed likely primarily due to the advanced thermal maturity of organic matter. This means that this feature cannot be used to withdraw solid conclusions on the paleoecology and/or taphonomy of the  Numbers (1,2,3,9) refer to the position of the methyl-substituent. The chromatograms are scaled to the most abundant peak in each trace. Note that kerogen-bound (methyl-) phenanthrene distributions provide no direct information about thermal maturity (Marshall et al. 2007;Pehr et al. 2021) Table 3 Selected hydrocarbon-based ratios for bitumens (B) and HyPy pyrolysates (HyPy). The carbon preference index (CPI) was calculated as 0.5 * [(C 25 − C 33 ) odd /(C 24 − C 32 ) even + (C 25 − C 33 ) odd /(C 26 − C 34 ) even ] (Bray and Evans 1961 Lakhanda Lagerstätte, as done previously for a sample from the Neryuen Formation (Pawlowska et al. 2013). The Lakhanda Lagerstätte contains possible traces of hopanes, but no detectable steranes or aromatic steroids (Figs. 4,6,7). This is in good accordance with findings from other Paleo-to Mesoproterozoic records worldwide where eukaryotic steranes are typically close to, or below, the detection limit (Dutkiewicz et al. 2003;Brocks et al. 2005;Blumenberg et al. 2012;Pawlowska et al. 2013;Flannery and George, 2014;Luo et al. 2015;Suslova et al. 2017;Nguyen et al. 2019;Zhang et al. 2021). Notably, hopanes and steranes have been found in ca. 780-729 million years old rocks from the Chuar Group (Arizona, USA) that are of late oil window maturity (Zumberge et al. 2020). Hence, the thermal maturity of organic matter in the samples analyzed herein (at or just beyond peak oil window) still allows for the preservation of both types of molecular fossils. The possible presence of trace amounts of hopanes and absence of steranes in the Lakhanda Lagerstätte thus requires another explanation.
Bacterial hopanoid synthesis seems to be much more widespread among aerobes, with only a few known exceptions (Rohmer et al. 1984;Thiel et al. 2003;Härtner et al. 2005;Blumenberg et al. 2006). The low abundance of hopanes despite about peak oil window maturity thus may reflects a dominance of anaerobic over aerobic bacteria in the Lakhanda ecosystem. Steranes decompose more rapidly than hopanes during experimental thermal maturation (Mißbach et al. 2016), but it is commonly assumed that both compounds possess similar long-term preservation potential in the geosphere (Love and Zumberge 2021). Accepting this latter assumption, the presence of trace amounts of hopanes and the absence of steranes could suggest a bacterially dominated environment with very little inputs by eukaryotes. Thus, in combination with the paleontological record, the molecular fossils seem to indicate that eukaryotes were present but not significant in the Mesoproterozoic Lakhanda Lagerstätte.
An alternative explanation for the dearth of eukaryotic steranes in Proterozoic samples that preserve hopanes might be taphonomic processes associated with microbial mats (Pawlowska et al. 2013). This interesting hypothesis is based on the premises that microbial mats (i) were much more widespread during the Precambrian and (ii) affected taphonomic processes in a uniform manner. Indeed, microbial mats may have ubiquitously sealed Precambrian sediments, at least in photic zone environments (Walter 1972;Seilacher 1999). However, given the large spatial and temporal extent, these microbial mats were likely microbiologically and biogeochemically divers. Preservation pathways of lipids including hopanoids and steroids are known to be extremely complex, and several studies documented the high variability of the involved multistep processes in modern microbial mats (e.g., Blumenberg et al. 2015;Shen et al. 2018Shen et al. , 2020Lee et al. 2019Lee et al. , 2021. For these reasons, the mat-seal hypothesis appears inadequate to explain the dearth of regular steranes for the Proterozoic molecular fossil record prior to ca. 820 million years.

Conclusions
Our study provides the first thorough investigation and documentation of molecular fossils in the Lakhanda Lagerstätte, including the first systematic and rigorous assessment of organic matter preservation (thermal maturity, hydrocarbon syngeneity). Programmed temperature pyrolysis data and several molecular maturity parameters independently suggest that the organic matter is about peak oil window maturity. Various criteria, as for instance compositional similarities of the bitumens and corresponding kerogens, demonstrate that sedimentary hydrocarbons are indigenous and syngenetic to the host rock. Notably, hopanes appear to be present in traces, while steranes were not detected. Since these molecular fossils should still be preserved at the inferred thermal maturity range, the observed features fit best to an environment dominated by anaerobic bacteria. Our study on molecular fossils in the Lakhanda Lagerstätte thus contributes to the view that eukaryotes were present but not significant in Mesoproterozoic ecosystems, supporting a protracted evolution of this organism group. T Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.