Fossilized anaerobic and possibly methanogenesis-fueling fungi identified deep within the Siljan impact structure, Sweden

Recent discoveries of extant and fossilized communities indicate that eukaryotes, including fungi, inhabit energy-poor and anoxic environments deep within the fractured igneous crust. This subterranean biosphere may constitute the largest fungal habitat on our planet, but knowledge of abyssal fungi and their syntrophic interactions with prokaryotes and their concomitant metabolisms is scarce. Here we report findings of fossilized, chitin-bearing fungal hyphae at ~540 m depth in fractured bedrock of the Siljan impact structure, the largest crater in Europe. Strong 13C-enrichment of calcite precipitated with and on the fungi suggests formation following methanogenesis, and that the anaerobic fungi decomposed dispersed organic matter producing for example H2 that may have fueled autotrophic methanogens. An Eocene age determined for the calcite infers the first timing constraint of fossilized fungi in the continental igneous crust. Fungi may be widespread decomposers of organic matter and overlooked providers of H2 to autotrophs in the vast rock-hosted deep biosphere. Fossilized anaerobic fungi are identified in the largest impact structure in Europe and may have fuelled methane production in energy-poor anoxic environments, according to microscale observations and isotopic analyses

T he deep biosphere is the largest microbial habitat on Earth by volume [1][2][3][4] . The subsurface encompasses diverse environments and ecosystems, including severely energy limited deep biosphere ecosystems, inhabited by ancient evolutionary microbial lineages, adapted to subsurface conditions 1,5,6 . Knowledge of these lifeforms and processes has important implications for life in extreme environments and evolution on Earth, as well as significance for the quest for life on other planets and icy moons 7 . Estimates of the total biomass in the continental deep biosphere range from a couple of percent to almost 20% of Earth's total 4,8 . Research on active communities in the deep biosphere has so far focused on prokaryotes. Recent studies have, however, shown that microorganisms of all three domains of life are abundant and active in deep aquifers in Precambrian crystalline rocks 9 , which make up the largest volumes of the continental crust 10 . Eukaryotes have been detected at several kilometers depth 11 , but the general trend is increased proportion of prokaryotes with declination 12 . Of the eukaryotes, organisms of the kingdom Fungi have been detected in several studies of deep present-day fracture waters 13,14 . Deep remnants of ancient fungi have been found in fractured igneous rock in Sweden 15,16 . Fungi are evidently inhabiting the deep igneous biosphere, but knowledge about their ecological role, and how long they have dwelled in this environment is far from resolved. From other settings, ancient remains of fungi or fungus-like structures have been detected in cavities and/or beddings of Neoproterozoic sedimentary rocks 17,18 and in vesicular basalt of Devonian 19 and Paleoproterozoic age 20 . The most convincing fungal fossils are Phanerozoic in age 21 .
The shift to anoxic conditions generally occurs within the upper tens of meters in continental igneous aquifers 22 suggesting that the majority of deep biosphere microorganisms, including fungi, have adapted to an anaerobic lifestyle. Most studies of obligate anaerobic fungi, and hence most well-established metabolic models, are of the phylum Neocallimastigomycota living in the rumens of herbivores. Neocallimastigomycota possess redox organelles called hydrogenosomes instead of mitochondria that produce H 2 (as well as CO 2 , lactate, and formate) as a metabolic waste product [23][24][25] . In the rumen, fungi consort symbiotically with prokaryotes, such as methanogens 26,27 , and symbiotic relationships between gut fungi and methanogens in relation to changing diets have been reported 28 . For the deep biosphere, anaerobic fungi have been shown to have distinct survival strategies to adapt to anaerobic ecosystems (e.g., ethanol fermentation and production of amino acids by Schizophyllum commune 20R-7-F01 in sediments at 2 km below the seafloor 29 ) and a similar type of syntrophic consortium as in rumens has been proposed, where fungi produce H 2 that can be used by indigenous autotrophic prokaryotes, such as sulfate reducers 16,30,31 . The hypothesis includes degradation of prokaryotic biofilms by heterotrophic fungi and production of H 2 as a waste product that then fuels the autotrophic communities, marking a potentially key ecological role for fungi in the deep biosphere 31 . When anaerobic fungi are enriched together with methanogens, they are able to degrade lignocellulose and produce more methane than bacteria and methanogens do 32 . This suggests that consortia of anaerobic fungi and methanogens in the deep subsurface hold the potential for methane production, not only from lignocellulose, but also from other sources of organic matter that may be more recalcitrant to degradation by bacteria and methanogens. This opens up a unique niche and new role for subsurface, fungalmethanogen consortia to fill, but in situ evidence of such consortia has not yet been reported form the deep subsurface.
Here, we explore the late Devonian (Frasnian) Siljan impact structure in Sweden for fossil fungi in anoxic environments and use imaging techniques, staining and biomarkers to identify and describe the fossilized fungal remains. Stable isotopes and organic molecular remains of coeval secondary minerals are used to link the fungi to prokaryotic metabolisms, such as methanogenesis. Microbial methane has formed in this impact structure at least since the Late Cretaceous 33 , making it an optimal site for exploration of potential ancient fungi-methanogen relationships. Our findings suggest a link between anaerobic fungi and methanogens at several hundred meters depth in the impact crater, in a putative process where fungi fueled methanogenesis via H 2 production during degradation of organic matter infiltrated from black shales 33 . This complements H 2 from other sources such as abiotic and prokaryotic fermentation, and manifests a potentially widespread and long-term ecological influence of fungi for methane formation in the vast deep biosphere environment. This may be of great importance, since methane is a very potent greenhouse gas 34 . We also establish the first timing constraints of anaerobic fungi in the continental deep biosphere, by linking it to previously reported U-Pb dating of coeval calcite 33 and, in addition, detect a tetrahymanol biodegradation product in the form of 30nor-gammacerane in the fungal remnants.

Results
The eroded Siljan crater has an age of 380.9 ± 4.6 Ma 35 and has down-faulted Devonian and Silurian sedimentary rocks with thicknesses of up to 400-500 m 36,37 in a circular crater rim depression (Fig. 1). Cored boreholes of up to 700 m depth penetrate the rim zone 38 . Open fractures in seven of these drill cores representing both the sedimentary units and the underlying Paleoproterozoic felsic host rock were screened for the presence of fossilized microorganisms. A methane-dominated gas was encountered at several depths during borehole drilling, both in the sedimentary units and in the deeper basement 33 . Bituminous material and seep oil of shale origin are scattered in the pore space and fractures of both sedimentary rocks 39,40 and the deeper igneous rocks 33 . A highly fractured and porous core section in granite-rhyolite at depths between 534 and 542 m in borehole CC1 contained scattered filamentous carbonaceous material resembling fungal hyphae. Five fractures in this section were sampled for detailed analysis. This rock section is >125 m below the contact to the Paleozoic sedimentary units and represents a sub-vertically dipping fracture network structure. The fractures are partly sealed by mineral infillings, but abundant open cavities with euhedral crystals remain, providing open pore space for colonization by microorganisms (Fig. 1c, d). The open cavities were exposed in the laboratory by opening the fractures physically with a chisel. The exposed fracture surfaces contain carbonaceous material, euhedral crystals of calcite, pyrite, and clay minerals grown on older euhedral quartz that lines the host rock contact (Fig. 1e). The morphology of the calcite crystals reveals more than one generation; an older core and an overgrowth (Calcite-1 and Calcite-OG, Fig. 1e).
The carbonaceous material occurs as thin films resembling extracellular polymeric substance (EPS) biofilms on the fracture surface, preferentially on quartz ( Fig. 2a) but also on the euhedral calcite crystals (Fig. 1e). From this presumed biofilm, filamentous structures protrude into the open cavity as single filaments or as intertwined mycelium-like networks (Fig. 2b, c) that are intergrown with the outermost parts of the calcite overgrowths. The filaments have diameters of~10 ± 5 µm, and extend as undulating features of more than a millimeter (Fig. 3a) from the onset at the biofilm (Fig. 3b). At contact with mineral surfaces the filaments widen in a tapered fashion similar to so called hyphal bridges known from fungal hyphae 38 (Fig. 3b). The filaments have septa at regular distances of around 50-60 µm (Fig. 3c), and branch frequently (Fig. 2c). Mineralization to dominantly clay minerals is a common feature and is partial (Fig. 2d) to complete (Fig. 2e). Fluorescence microscopy confirmed that the filaments were positively stained with Calcofluor White, which indicates the presence of chitin in some parts of the filaments (Fig. 4f). Parts of the filaments did not show a positive reaction, presumably due to degradation and mineralization.
The outermost growth zone of the euhedral calcite crystals is intergrown and partly overgrown with filaments (Figs. 1e, and 4a). Synchrotron radiation X-ray tomographic microscopy (SRXTM) shows that mineralized filaments exist within the euhedral calcite crystals, and that the calcite overgrowth has completely covered the suggested organisms in places (Fig. 4b). Environmental Scanning Electron Microscopy (ESEM) reveals that the individual filaments feature <5 µm sized precipitates of clay minerals and euhedral calcite crystals on the surface of the filaments (Fig. 4c). Carbonaceous material appears as small filamentous-like ridges with widths of~1 µm and length of a few tens of micrometers, on the surfaces of the main filaments (Fig. 4c). These smaller filamentous-like features appear to form networks that transform into carbonaceous biofilm-like material interpreted as remnants of EPS, which contains calcite mineralization.
Microcrystalline aggregates of pyrite and late-stage calcite occur within the carbonaceous mycelium-like network, and are physically separated from the fracture wall and older euhedral pyrite and calcite (Calcite-1) coating (Fig. 4d, e and Supplementary Fig. 1). Spheres with diameters ranging between 5 and 15 µm occur on the tip of the hyphae (Fig. 3d) and as individual features on the fracture surface ( Supplementary Fig. 2). Spheres within the smaller diameter range occur as connected pairs which resemble budding cell division (binary division). Orbs within the larger diameter range are usually singular. The morphology of the spheres and filaments is intact in the SEM and their chemical composition is close to identical ( Supplementary Fig. 3).
The organic extracts of the fracture minerals including parts of the fungal biofilm contained predominantly highly degraded organic matter in an unresolved complex mixture, in common with bitumen-rich coatings in the fracture network 33 . Hydrocarbons C 16 -C 35 detected with gas chromatography-mass spectrometry (GC-MS) show a bimodal distribution with maxima at C 19 and C 25 . Intriguingly, 30nor-gammacerane is detected, along with a number of hopanoids ranging from C 27 to C 35, but no steranes (Fig. 5a). This 30nor-gammacerane is most concentrated in the organic extract of permineralized hyphae together with elemental sulfur (Fig. 5b). In contrast to the mineral-dominated sample, the permineralized hyphae sample also contains fatty acids between C 14 and C 22 , dominated by C 16 and C 18 (full chromatogram not shown). Minor but distinct peaks of unsaturated C 16:1 and C 18:1 occur. Odd-chain fatty acids, especially C 15; C 15i and C 15ai , and C 17 ; C 17i and C 17ai are also detected.
Secondary ion mass spectrometry (SIMS) microanalysis transects through calcite crystal cross-sections (Fig. 6) show a shift from relatively isotopically light δ 13 C values in the inner part of the calcite crystals (Calcite-1) to isotopically heavy in the outermost overgrowth that embeds the filaments (Calcite-OG). The δ 13 C calcite values are as high as +8.1 ± 0.4‰ (n SIMS = 14, median +6.0‰ in Calcite-OG, Supplementary Data 1). Pyrite occurring within the mycelium-like network (Fig. 4d, e) has isotopically light and relatively homogenous δ 34 S values of −40 ± 1‰ (n SIMS = 37) as presented previously 33 .

Discussion
The size and septation of the filaments as well as their occurrence in mycelium-like networks with frequent branching is in agreement with a fungal, rather than a prokaryotic, interpretation 41 . The mineralization to clay minerals is also a common fungal feature 42,43 and the chemistry derived by energy dispersive X-ray spectroscopy (EDS, Supplementary Fig. 4) resembles previous findings of smectite/montmorillonite-dominated clay-mineralization of fungal hyphae 42 . Such mineralization of organic matter prevents further degradation and so explains the exceptional preservation of the mycelial network 44 . The detection of chitin further strengthens a fungal interpretation and excludes a prokaryotic affinity. Chitin can be resistant to degradation and diagenesis for long periods of time, and has previously been detected in fossilized fungi from the oceanic crust of Miocene age 45 , in 2 Ma old fungal fossils from a paleo-hydrothermal system 46 , as well as in fungal fossils of an alleged Neoproterozoic age in dolomitic shale of the Mbuji-Mayi Supergroup, Democratic Republic of Congo 17 . The behavior of hyphae originating from a basal film corresponds to previously described fossilized hyphae of the deep biosphere 16,21 . The hyphal size and morphology, including repetitive septa suggests an affinity to Dikarya (ascomycetes or basidiomycetes), but a conclusive b Geological map of Siljan with location and simplified geology of cored borehole CC1 shown (no detailed mapping of the sedimentary units is available), including the fractured section at 534-542 m sampled for detailed studies. The map, excluding the drill core log, is modified from 33 and is a electronically re-drawn in Adobe Illustrator after a map originally presented in the ref. 87 . c Drill core section with partly open fractures and cavities in rhyolitic rock. d Exposed surface of the walls of mechanically opened cavities from c. e Back-scattered ESEM-image of minerals and carbonaceous and mineralized filaments (arrows "Hypha" and "Mineralized hypha") on the surface of the exposed fracture. Minerals are quartz (Qz), pyrite, clay minerals, and calcite. The calcite crystals show an older core ("Calcite-1") and a younger overgrowth ("Calcite-OG", borders marked in one of the crystals by a gray line).
affinity of the fungal fossils is challenging to determine. The spheres of small diameter exhibit likely cell division, and the single spheres of larger diameter correspond well with budding yeast and subsequent growth of single yeast cells. Yeasts have previously been reported from~400 m depth in an anaerobic crystalline bedrock aquifer at Äspö, Sweden 13 . We suggest the spheres in our material represent fossil yeast-like organisms, although the interpretation is uncertain as it is merely based on morphology. Yeasts are facultative anaerobes and can survive under anoxic conditions throughout their lifecycle 47 . We posit that the contemporaneous occurrence of filaments and spheres, combined with their identical composition and preservation, indicates mold-like and yeast-like growth stages of a dimorphic fungus, rather than two different fungal taxa. Dimorphic fungi can exist as a yeast-like or a filamentous phase and can switch between these two during their lifecycle 48 . We hypothesize that the ability to grow and switch between a multicellular hyphal stage and a unicellular yeast confers a selective advantage for survival and colonization in the oligotrophic deep subsurface biome. Hyphal growth facilitates the exploration for nutrients and carbon, whilst the yeast form enables rapid reproduction and widespread dissemination of cells via fluids in rock fracture systems.
The findings of the fungal fossils at depths well below the redox front 22,49 , together with the associated unaltered sulfide minerals, indicate fungal growth under strictly anoxic conditions. Facultative anaerobic fungi are known from all major fungal divisions and are frequently found in extreme environments including the deep subsurface 14 . Obligate anaerobic fungi, however, are largely from the phylum Neocallimastigomycota and are best known from the rumen of mammalian herbivores but have been reported from the gut and coelomic fluid of the coastal sediment dwelling sea urchin Echinocardium cordatum 50 . They also occur in the guts of the algae-grazing marine iguana Amblyrhynchus cristatus 51 . They have been reported from environmental sites such as landfill soils 52 as well as lacustrine 53 , estuarine 54 , and marine sediments 55 . Thus, a growing number of observations indicate that Neocallimastigomycota exist not only as endosymbionts in ruminants but also in marine animals. It still remains unclear however, whether the findings of free-living Neocallimastigomycota in anoxic sediments and soils, which are based on molecular data alone, are due to contamination, as no free-living example organisms have been enriched or isolated yet. The growing literature of extant and fossil fungi from strictly anoxic environments indicates a gap in our understanding of anaerobic fungi across the fungal taxonomic tree, and that the abundance and diversity of obligate anaerobic fungi are far from well understood.
The detection of 30nor-gammacerane exclusively in the hyphae and not in the other carbonaceous materials is interesting considering the isolated and strictly anoxic nature of the samples in which the fossil fungi are found. Nor-gammacerane is the decomposition product of tetrahymanol, a key molecule in anaerobic metabolism in various micro-eukaryotes including fungi 56,57 . Even though the precursor of tetrahymanol may be unclear, this lipid is only produced under anaerobic conditions, and thus in the geological record, gammacerane is a useful ecosystem biomarker and proxy for anoxic settings. Sterols play essential roles in cell physiology, especially in cell membranes, in most eukaryotes. Sterols requires molecular oxygen at multiple steps of biosynthesis. Even though the enzymes responsible can function at extremely low levels of oxygen, eukaryotes need additional strategies to survive in the absence of oxygen 58 . Possible strategies for anaerobic eukaryotes are either to acquire sterols from the environment, or use alternative molecules. Sterols, or the diagenetically degraded forms, steranes, can be stable for long periods of time as the primary degradation pathways are oxygen dependent 59 . Steranes can be detected at substantial quantities in various environments throughout the geologic record, and therefore act as available substrates for scavenged by other eukaryotes 60 . Most eukaryotes have transporters for acquisition of sterols from the environment. The yeast Saccharomyces cerevisiae, responsible for production of wine or beer is an example. Oxygen depletion prior to alcohol fermentation makes the yeasts unable to produce sterols. Instead they acquire sterols from the environment using Aus1 and Pdr11 transporters, which have homologs in many other fungi. Presence of either transporter is enough for growth during fermentation, but loss of both is lethal 61 . In isolated and strictly anoxic environments, eukaryotes can neither produce sterols nor scavenge off other steroid-bearing eukaryotes. Instead, they use a sterol-like molecule that can be produced in the absence of oxygen-tetrahymanol or other hopanoids-as an alternative 62 . Neocalligmastigomycota produce tetrahymanol as a characteristic membranereinforcing lipid instead of sterols as an adaptation to oxygen-free environments 58,63 . The rumen fungus Piromyces communis provides such an example, and lacks both sterols and homologs of the transporter Aus1 in its genome, but has tetrahymanol in the cell membranes 63 . Piromyces, as well as some other eukaryotes living in anoxic environments, can synthesize this molecule owing to a laterally transferred gene from bacteria 64 . This suggests that such fungi may have unique strategies for survival in environments permanently lacking oxygen 29 . Discoveries of actively anaerobic and facultatively anaerobic fungi in the deep biosphere, belonging mainly to the Ascomycota phylum 14 , suggest that this adaptation may be more widespread, but the tetrahymanol pathway has not been confirmed among deep biosphere fungi previously.
The traces of bitumen have infiltrated down into the crystalline basement fractures 33 from the organic-rich Upper Ordovician Fjäcka shale existing in the sedimentary successions of the crater rim 40 . The bitumen has undergone substantial thermal maturation and biodegradation 39,40 . Seep-oil and solid bitumen in  c Close-up ESEM-image of a filament (arrow "Hypha"), featuring tiny precipitates of clay minerals and euhedral calcite crystals (<5 µm in size, composition in Supplementary  Fig. 4) and filamentous carbonaceous material on the surface of the main filament (arrow "EPS", carbonaceous composition documented in Supplementary Fig. 5). The EPS grades from carbonaceous to completely mineralized. d Back-scattered ESEM-image of microcrystalline pyrite aggregates on mineralized mycelium together with organically preserved filaments (arrows "Hypha"). e SRXTM surface-rendering of the mycelium (orange) together with fine-grained pyrite (bright yellow) and calcite (dark purple) on the older calcite fracture coating in the cavity. f Staining of a filament. The filament parts are either showing positive reaction for chitin, or no reaction due to mineralization. The hypha was detached from the fracture surface and background is glass vial and staining fluid. CFW Calcofluor White.
limestone fractures have carbon numbers of steranes and diasteranes, which might indicate predominantly terrigenous higher plant input deposited in oxidizing environments 40,65 . Investigations of bitumen in deeper fractures of the crystalline basement, however, indicate an absence of steranes 33,65 . Because of the molecular structure it is likely that steranes in the seep oils were incorporated in locally fixed kerogen structures and did not migrate through the fracture network. As the fungus lived in an anoxic part deep within the fracture network, where sterols could neither be readily acquired from the environment nor produced, it would have been necessary to find alternative ways to build sterol-like molecules to maintain cell wall stability. Although  hopanoids synthesized via oxygen-independent pathways could also substitute position and function of the sterols in the membrane, the only known pathway for sterol substitution in anaerobic fungi is via tetrahymanol. The latter is in line with our findings of 30nor-gammacerane (Fig. 5). 30nor-gammacerane may thus be an important biomarker for anaerobic fungi in the rock record.
The fungal hyphae are succeeded by, and intergrown with, calcite overgrowths (Figs. 4 and 6). The outermost growth zone of the euhedral calcite crystals has been dated with high spatial resolution U-Pb geochronology to an age of 39.2 ± 1.4 Ma 33 . This marks the age of the fungi embedded in this calcite growth zone. However, partly organically preserved hyphae in the cavities are likely of younger age. Sampling of the drill core immediately after drilling assures that the fungi are not modern contaminants. The ancient nature of the fungal hyphae is supported by several lines of evidence, including 1) the partly to completely mineralized and degraded nature of the hyphae, as well as absence of any recent functionalized molecules in the extracted organic matter (e.g., intact lipids, or sugars and amino acids common in EPS), 2) presence of sulfides on the hyphae confirming formation under anaerobic conditions and not at oxic, atmospheric conditions after core retrieval, and 3) no collapsing of the filaments or spheres under vacuum in the SEM.
Relatively isotopically light, 12 C-rich, methane is produced during microbial methanogenesis, and because of this isotopic discrimination against 13 C, the residual CO 2 becomes 13 C-rich 66 . Incorporation of the residual carbon in authigenic carbonate minerals inherits the 13 C-rich signature and can therefore be used as a marker for methanogenesis 67 . The significantly 13 C-enriched calcite observed in the fungus-related calcite overgrowth (δ 13 C calcite values as heavy as +8.1‰, Fig. 6) is therefore proposed to reflect precipitation following microbial methanogenesis in situ. This is in agreement with previous widespread findings of isotopically heavy calcite in this fracture system, with δ 13 C calcite values up to +21.5‰ 33 . The older calcite (Calcite-1) with lighter δ 13 C calcite values is, however, not linked to methanogenesis based on the C isotope values (Fig. 6). Instead, onset of methanogenesis occurred in temporal relation with the appearance of fungi. The significantly 13 C enriched calcite values point to methanogenesis through the carbonate reduction pathway 35,68 , the dominant terminal process in petroleum biodegradation in the subsurface 69 . The in situ formation of microbial methane through the carbonate reduction pathway is further supported by isotopically light methane (δ 13 C CH4 : −64 ± 2‰) and isotopically heavy carbon dioxide (δ 13 C CO2 : +5-9‰) sampled in boreholes nearby 33 . As this is an active gas prospecting site with ongoing production well tests, no sampling for active microbial communities in the water has been possible, and hence the presence of active methanogenic communities has not been possible to determine.
The supply of the electron donor H 2 needed for microbial methanogenesis through CO 2 reduction may originate from fermentation of hydrocarbon related organic matter, in addition to geological sources (hydrolysis) in a first step and utilization of H 2 by autotrophic methanogens in a second step (CO 2 + 4H 2 = CH 4 + 2H 2 O) 69 . Based on the spatial link between 13 C-rich calcite and fungal hyphae, we propose that anaerobic fungi that produce H 2 as a metabolic waste product [23][24][25] have facilitated autotrophic methanogenesis through degradation of organic matter, in addition to H 2 originating from prokaryotic fermentation degradation. The significance of fungal derived H 2 versus prokaryotic derived H 2 for methanogenesis cannot be assessed in this fossil system, but the presence of 13 C-rich calcite on the hyphae suggests a close connection between the fungi and methanogenesis. The available data do not suggest the presence of prokaryotic fermenters, but the possibility cannot be ruled out since they are known from subsurface settings 70,71 . Culture-dependent evidence of fungi-methanogen relationship has, so far, only been established in gut biomes and not in environmental samples 72,73 . However, fungi have been largely neglected in microbial ecology, not least in the deep biosphere, and anaerobic fungi as possible key players in anoxic systems is a relatively new concept 30,31 . Isolation and culture-dependent investigations of deep biosphere biota will hopefully be more frequently used in the future and shed more light on subsurface ecosystems. An in vitro approach using "synthetic" co-cultures of various candidate filamentous fungi and methanogens is of particular interest, as well as the enrichment, isolation, and co-culturing of native subsurface fungal-prokaryote community members. The organic matter may be remnants of bitumen that the fungi can degrade, and/or prokaryotic metabolites and mats 16 , the main constituent of dissolved organic matter at large depth within this type of setting 74 . It has been shown that soil-derived fungi are able to produce methane through biodegradation in relationship with methanogens 75 , but also without 76 . In the latter case, methane was produced by saprotrophic aerobic fungi, but so far, no anaerobic fungi are known to produce methane.
Microbial sulfate reduction (MSR) can also be involved, both in the oxidative degradation of organic matter of sedimentary origin, but also by autotrophic utilization of H 2 73 . Pyrite formed by MSR is typically strongly depleted in 34 S 77,78 . The very low δ 34 S pyrite values (−40 ± 1‰, V-CDT) 33 of pyrite precipitated in the mineralized mycelium (Fig. 4d, e) are thus proposed to reflect MSR, in similarity with findings of fungal hyphae with 34 Sdepleted pyrite at 740 m depth at Laxemar, Sweden 16 . The sulfate reducers were either feeding on the H 2 or other fungal waste products, or scavenging the fungal biomass as a direct energy source. The sulfate reducers and methanogens were probably not completely coeval as the former outcompete methanogens for H 2 and other substrates when sulfate concentrations are elevated 79 . Exhaustion of sulfate by MSR, or infiltration of fresh water may have resulted in temporal fluctuation of conditions in favor of methanogenesis.
The smaller filamentous EPS structures on the larger fungal hyphae are closely connected to calcite with methanogenesis signals ( 13 C-rich). EPS formation is common in mixed fungal-bacterial biofilms 80 , but in the case presented here, neither body fossils resembling archaea, nor diagnostic biomarkers of methanogens remain in the fossilized biofilm on the fracture surface. This is in line with previous observations of ancient deep igneous-rock hosted microbial communities 16 .
The upper crystalline continental crust environment makes up one of the largest, but yet least explored, deep biosphere habitats on Earth. Considering the vastness of this biome, it is reasonable to assume that these microbial processes affect the energy cycles of our planet, but to what degree remains largely elusive. Of particular interest is the in situ microbial formation of methane, which if/when released to the atmosphere is a potent greenhouse gas 34 . The extent, continuity and physicochemical prerequisites for methane accumulation require more attention in order to assess the significance of this underexplored greenhouse gas source on a global scale. Autotrophic methanogenesis requires a supply of H 2 as electron donor. In specific environments, such as serpentinized systems and at subduction zones, the flux of abiotic H 2 may be large and can therefore provide enough energy for autotrophic communities such as methanogens 81 . In granitoid-hosted aquifers of Precambrian cratons, several different processes have been proposed for the origin of H 2 , including abiotic radiolysis 82 that may build up in the waters over the long residence times recorded for deep brine fluids at kilometers depths 83 . Igneous-rock hosted groundwaters from the Fennoscandian shield show scattered elevated H 2 concentrations, of up to~0.2 mL/L 84 , that are neither correlated with depth nor with residence times. The residence times of these waters, in the order of a couple of thousand years 85 , are too short for build-up of significant H 2 concentrations by radiolysis (cf. 82 ). Instead, our results suggest that anaerobic fungi may be a neglected source of H 2 , that fuels autotrophic methanogens in the deep continental subsurface and complement potential prokaryotic degradation pathways of organic matter involved in H 2 formation, such as fermentation. The heterotrophic fungal metabolism can be sustained either by scavenging of descended organic matter 33,41 , or feeding off of indigenous organic prokaryotic metabolites suspended in the water 74 or available in microbial mats 16 , and we propose future co-culturing studies to gain more knowledge about the syntrophic relationships of fungi and methanogens in these deep environments. Our results also suggest that impact craters enable deep colonization by eukaryotes in the pore space and connected conductivity that result from deep crustal fracturing during meteorite bombardments, partly by introducing surficial organics to great depth. Impact events may thus have enabled intermittent eukaryotic propagation into new ecological niches in the deep biosphere throughout Earth's history.

Methods
Materials. Open fractures from a drill core CC-1 at Mora, Sweden ( Fig. 1) were sampled for carbonaceous filaments and coeval secondary fracture-coating minerals. Samples are from 534 to 542 m vertical depth from the ground-surface, in igneous rock fractures. Calcite mainly occurs as euhedral crystals of short c-axis type (Fig. 1e). Polished cross sections reveal growth zonation in several calcite crystals (Fig. 6a). The paragenesis includes pyrite, clay minerals, harmotome, apophyllite, sphalerite, galena, and quartz but most of these minerals are related to the oldest growth phases of calcite (exceptions are pyrite and clay minerals). Carbonaceous material resembling fragments of solid bitumen are present. The borehole was drilled at the crater rim zone through an upper layer of depressed Paleozoic sedimentary rocks into the deeper, fractured Proterozoic crystalline bedrock. The interface between the rock units is at 406 m in the studied borehole. The mineralogy and appearance of the fracture coatings were examined using a Hitachi S-3400N Scanning Electron Microscope (SEM) equipped with an integrated EDS system and an Environmental SEM (ESEM) with a FEI QUANTA FEG 650 (Oxford Instruments, UK) EDS and an Oxford T-Max 80 detector. The analyses were performed in low vacuum to minimize surficial charging effects. This enables the use of uncoated samples and, thus, EDS analyses of the C content. The acceleration voltage was 20 or 15 kV depending on the nature of the sample, and the instrument was calibrated with a cobalt standard. Peak and element analyses were done using INCA Suite 4.11 software and normalized to 100 wt%. Element mapping was done using Aztec software. The calcite crystals were handpicked under the microscope for stable isotope analysis and carbonaceous filaments were hand-picked for staining and biomarker analysis. Larger fragments of the mycelium-like structures and mineral coatings were detached from the drill core sample for SRXTM.

SIMS.
Calcite crystals related to hyphae were mounted in epoxy, polished to expose cross-sections and examined with SEM to trace zonation and impurities prior to SIMS analysis. SIMS-analysis (10 μm lateral beam dimension, 1-2 μm depth dimension) of carbon isotopes in calcite was performed on a Cameca IMS1280 ion microprobe following the analytical settings and tuning reported previously 86 . Analyses were performed in automated sequences, with each analysis comprising a 70 s pre-sputter to remove the gold coating over a rastered 15 × 15 µm area, centering of the secondary beam in the field aperture to correct for small variations in surface relief and data acquisition in sixteen four second integration cycles. The magnetic field was locked at the beginning of the session using an NMR field sensor. Secondary ion signals were detected simultaneously using a Faraday detector/Electron Multiplier combination with mass resolution 2500 on the 12 C peak and 4000 on the 13 C peak to resolve it from 12 C 1 H. Results are reported as per mil (‰) δ 13 C based on the Pee Dee Belemnite (V-PDB) reference value. Data were normalized for instrumental mass fractionation (IMF) using matrix matched reference materials mounted together with the sample mounts and analyzed after every sixth sample analysis. The calcite reference material, S0161 from a granulite facies marble in the Adirondack Mountains, was kindly provided by R.A. Stern (Univ. of Alberta). The values used for IMF correction were determined by conventional stable isotope mass spectrometry at Stockholm University on ten separate pieces, yielding δ 13 C = −0.22 ± 0.11‰V-PDB (1 std. dev.). Precision was ±0.3-0.5‰. Values of the reference material measurements are listed together with the samples in Supplementary Data 1.
GC-MS. The fungi sample and a sample with calcite and fungi were gently acetone flushed to remove surface contaminations and then ground with an agate pistil. The sample powders (30 mg) were extracted in Teflon-capped glass vials using 2 mL of pre-distilled dichloromethane/methanol (ultrasonication, 15 min, 40°C) and the supernatant was decanted after centrifuging. Extraction was repeated twice, with dichloromethane and hexane as solvents. Once evaporated, the combined extracts and re-dissolution in pure dichloromethane, the solvents were dried with N 2 . The total organic extract was derivatized by addition of 50 μL BSTFA (50°C, overnight). After the sample was dried with N 2 , it was mobilized with 20 µL n-hexane and stored frozen at −18°C. The sample remnants were dissolved and demineralized by adding 5 mL of TMCS/methanol (1 + 9) for 12 h and then derivatized for 90 min (80°C). After cooling, the samples were mixed with n-hexane and the decanted supernatants collected separately. This procedure was repeated three times. The combined supernatants were dried with nitrogen and remobilized with 50 µL nhexane. Each sample extract (1 µL) was analyzed via on-column injection into a Varian CP-3800 GC/1200-quadrupole MS (70 keV) equipped with a fused silica column (Phenomenex ZB-5; 30 m length, 0.32 µm inner diameter, 0.25 μm film thickness). The GC oven programming was from 80°C (held 3 min) to 325°C (held 40 min) at 6°C min −1 . Carrier gas was He, at 1.4 mL min −1 . Compounds were assigned by comparison with spectral databases and published mass spectral data.
SRXTM. The tomographic measurements were carried out on the TOMCAT beamline at the Swiss Light Source, Paul Scherrer Institute, Villigen, Switzerland. Projections were acquired equiangularly over 180°, post-processed online and rearranged into flat-corrected and darkfield-corrected sinograms. The beam energy used for the five aliquots was respectively 20 (n = 2), 23 (2), and 25 (1) keV, for maximum absorption contrast. Specimens were scanned with objectives 10× (n = 2, exposure time 300 ms) and 20× (n = 3, exposure time 1300 ms) at two different sessions. During the scanning process a LuAg:Ce 20 µm scintillator was used. Visualization has been done using Avizo 9.5.0 (FEI Company). With the 20× lens used, the resulting voxel size was 0.325 μm and, respectively, 0.65 µm with 10×.
Staining. Samples were treated with sterile stainless steel forceps and not touched by ungloved hands prior to staining to reduce introduction of fluorescence particles. A few single fungal hyphae were detached from the fracture surface under the microscope, and investigated under fluorescence microscopy using a Olympus BX51 microscope connected to an external Lumen Dynamics X-CITE series 120 Q, prior to staining to exclude autoflourescence. The samples were stained with Calcofluor White (ordered from VWR) for ten minutes and investigated with fluorescence light illumination under the microscope.

Data availability
All relevant data are included in the Supplementary Information of this article and stored publically available in the archive at Swedish National Data Service (doris.snd.gu.se) Received: 15 September 2020; Accepted: 19 January 2021;