A 16S rRNA Gene-Based Metabarcoding of Phosphate-Rich Deposits in Muierilor Cave, South-Western Carpathians

Muierilor Cave is one of Romania’s most important show caves, with paleontological and archeological deposits. Recently, a new chamber was discovered in the cave, with unique yellow calcite crystals, fine-grained crusts, and black sediments. The deposits in this chamber were related to a leaking process from the upper level that contains fossil bones and a large pile of guano. Samples were taken from the new chamber and another passage to investigate the relationship between the substrate and microbial community. Chemical, mineralogical, and whole community 16S rRNA gene-based metabarcoding analyses were undertaken, and the base of the guano deposit was radiocarbon dated. Our study indicated bacteria linked to the presence of high phosphate concentration, most likely due to the nature of the substrate (hydroxyapatite). Bacteria involved in Fe, Mn, or N cycles were also found, as these elements are commonly identified in high concentrations in guano. Since no bat colonies or fossil bones were present in the new chamber, a high concentration of these elements could be sourced by organic deposits inside the cave (guano and fossil bones) even after hundreds of years of their deposition and in areas far from both deposits. Metabarcoding of the analyzed samples found that ∼0.7% of the identified bacteria are unknown to science, and ∼47% were not previously reported in caves or guano. Moreover, most of the identified human-related bacteria were not reported in caves or guano before, and some are known for their pathogenic potential. Therefore, continuous monitoring of air and floor microbiology should be considered in show caves with organic deposits containing bacteria that can threaten human health. The high number of unidentified taxa in a small sector of Muierilor Cave indicates the limited knowledge of the bacterial diversity in caves that can have potential applications in human health and biotechnology.

Muierilor Cave is one of Romania's most important show caves, with paleontological and archeological deposits. Recently, a new chamber was discovered in the cave, with unique yellow calcite crystals, fine-grained crusts, and black sediments. The deposits in this chamber were related to a leaking process from the upper level that contains fossil bones and a large pile of guano. Samples were taken from the new chamber and another passage to investigate the relationship between the substrate and microbial community. Chemical, mineralogical, and whole community 16S rRNA genebased metabarcoding analyses were undertaken, and the base of the guano deposit was radiocarbon dated. Our study indicated bacteria linked to the presence of high phosphate concentration, most likely due to the nature of the substrate (hydroxyapatite). Bacteria involved in Fe, Mn, or N cycles were also found, as these elements are commonly identified in high concentrations in guano. Since no bat colonies or fossil bones were present in the new chamber, a high concentration of these elements could be sourced by organic deposits inside the cave (guano and fossil bones) even after hundreds of years of their deposition and in areas far from both deposits. Metabarcoding of the analyzed samples found that ∼0.7% of the identified bacteria are unknown to science, and ∼47% were not previously reported in caves or guano. Moreover, most of the identified human-related bacteria were not reported in caves or guano before, and some are known for their pathogenic potential. Therefore, continuous monitoring of air and floor microbiology should be considered in show caves with organic deposits containing bacteria that can threaten human health. The high number of unidentified taxa in a small sector of Muierilor Cave indicates the limited knowledge of the bacterial diversity in caves that can have potential applications in human health and biotechnology.
The opening of the first show cave in the early 17th century (Vilenica Cave, Slovenia) was just the beginning of the development of speleo-tourism (Cigna, 2019). Since then, more than 600 caves have been visited by tourists worldwide, receiving more than 25 million visitors per year (Cigna and Burri, 2000;Novas et al., 2017;Tièar et al., 2018). With the growing number of tourists, proper management of the caves' natural heritage is required against the impacts affecting cave habitats and microclimate (Pulido-Bosch et al., 1997;Hoyos et al., 1998;Northup and Lavoie, 2001;Mann et al., 2002;Paksuz and Özkan, 2012;Bercea et al., 2018Bercea et al., , 2019Debata, 2020;Constantin et al., 2021). However, the biological resources of caves can be reservoirs for pathogens and may be regarded as potential environments for pathogen transmission (Dean, 1957;Jurado et al., 2010;Igreja, 2011). Consequently, visitors must be informed of the threats, particularly where the management of show caves is rudimentary (Constantin et al., 2021).
The increasing interest in caves also encompasses the potential biotechnological applications in biological and medical sciences of the cave microbiomes (Jiang et al., 2015;Rangseekaew and Pathom-Aree, 2019;Zada et al., 2021). Already known bacteria reveal new characteristics, and the potential of pathogenicity has been found in some bacteria (Jurado et al., 2010;Moldovan et al., 2020). Alongside these discoveries, unreached sides of metabolism in organisms are starting to be discovered and understood.
Herein, we present the results of the studies on the microbiome of different deposit types in a protected sector of Muierilor Cave (Romania) to unravel the diversity of microorganisms related to the substrate chemistry and mineralogy and identify their possible sources. Muierilor Cave is one of the most important caves from paleontological, biological, and archeological points of view . With more than 130,000 visitors each year, it is also one of the most visited show caves in the country (Burghele et al., 2018;Constantin et al., 2021;Mirea et al., 2021).

Site Description and Sampling
Muierilor Cave (45 • 11 31.78 N and 23 • 45 14.07 E) is developed in limestone and located in Baia de Fier, south-western Romania, and is one of the most-visited show caves in the country (Figure 1). The cave is situated on the right slope of Galbenul Gorges at 650 m.a.s.l and is developed on four distinct levels. The second one has been a tourist and scientific attraction since the late 19th century. It is the first cave to be fitted with electric lights from Romania since 1963. The cave is inhabited by four species of hibernating bats, with the most abundant being Miniopterus schreibersii and Rhinolophus ferrumequinum, each forming colonies of around 1,000 individuals, especially along the touristic path.
Muierilor Cave is also known for the rich presence of hydroxyapatite [Ca 5 (PO 4 ) 3 (OH)] on rocks and boulders in different sectors of the cave. Hydroxyapatite is a naturally occurring mineral found in the bones and teeth of mammals, which is white when pure and yellow, green, or brown in nature. The lower level is the Scientific Reserve, from where the studied samples were collected, specifically from the Yellow Chamber.
The Yellow Chamber was resurveyed in 2019 during the paleontological excavation in the Scientific Reserve . It was named after the yellow calcite crystals that cover the rock surface in some of its parts. The Yellow Chamber is situated roughly underneath the Guano Chamber (which owes its name to a big pile of guano) located along the Touristic Passage (Figure 1). The Touristic Passage is also known for the abundance of fossil remains (belonging mainly to Ursus spelaeus sensu lato) . Part of the Yellow Chamber floor is covered with a black deposit of sandy clay with very sparse drips of water coming from the upper level, while two big boulders in its middle are covered with white and black crusts at a few tens of centimeters apart. The fragile crusts cover a small surface of the boulders, and at present, no water is dripping in that area. None of the other chambers or passages of the cave have such features, and they have not been seen in other Romanian caves. The ceiling of the Yellow Chamber is about 14 m at the highest point. The walls are covered with seemingly an organic deposit brought by the percolating water from the upper level.
The Yellow Chamber is a side passage with little ventilation. The geomorphological features develop on secondary faults (perpendicular to the main N-S fault) and connect with the upper levels through shafts and fissures. This feature allows the percolating water to precipitate in the lower levels. The temperature at the entrance of the Yellow Chamber (registered with a permanent data logger) was about 10 • C, and the humidity was >92%. Spot measurements of CO 2 were around 560 ppm (Constantin et al., 2021).
Sediment samples were collected in sterile Falcon tubes for microbiological analysis. A sample of 10 g of sediments was taken directly from the Yellow Chamber's black, sandy clay floor in a Falcon tube (PM11), and about 2 g of each crust was sampled with a sterilized scalpel (PMW = white crust, PMB = black crust). An additional sample of 10 g of brown sandy clay sediments was taken directly in a sterile Falcon tube from a different passage of the Scientific Reserve (PM7). The samples were transported for further laboratory analysis in an icebox and kept in the freezer at -60 • C until extraction.
For the chemical and mineralogical analysis, 10 g of sediments from PM7 and PM11 and 2 g of both crusts were collected in clean plastic bags and transported to the laboratory for analysis.

Dating of Guano Deposits
A guano sample was collected with a small corer from the bottom of the guano pile in the touristic passage for radiocarbon dating. The AMS 14 C dating was done commercially at the Poznan Radiocarbon Laboratory (Poland), following the procedures described in Mirea et al. (2021). The calibration of AMS 14 C ages was performed using the program OxCal ver. 4.4.2 (Bronk Ramsey, 2009, 2020) against the INTCAL13 radiocarbon calibration curve (Reimer et al., 2020).

Chemical and Mineralogical Analysis
The pH was measured using a Seven Excellence Multiparameter (Mettler Toledo, Greifensee, Switzerland) in a 1/5 (m/v) solid to water suspension. N and C were measured using a Flash 2000 CHNS/O analyzer (Thermo Fisher Scientific, Waltham, MA, United States) by combustion of 2 mg of sample. For metal, S, and P analysis, 1 g of dried and grounded sample (<250 µm) was digested in 28 ml of 3/1 (v/v) HCl (37%)/HNO 3 (62%) mixture. The concentration of Na, Mg, K, Ca, Al, Fe, P, and S was measured by inductively coupled plasma optical emission spectrometry using a 5300 Optima DV (Perkin Elmer, Waltham, MA, United States) spectrometer. The concentration of trace elements was determined by inductively coupled plasma mass spectrometry (ICP-MS) using an Elan DRC II (Perkin Elmer, Waltham, MA, United States) spectrometer.
To preserve crust deposits as well as possible, the mineralogical analysis was performed only on sediment samples. Sampling for mineralogical analysis could lead to the destruction of crusts and subsequently could limit the possibility of further studies. Samples were analyzed using powdered X-ray diffraction method with Rigaku Ultima IV diffractometer in parallel beam geometry equipped with CuKα radiation (wavelength 1.5406 Å). The XRD patterns were collected in the 2θ range between 5 and 80 with a speed of 2 • /min and a step size of 0.02 • . Rigaku PDXL Software connected to the ICDD database was used for phase identification. The quantitative determination was made using the RIR (reference intensity ratio) methodology.

DNA Extraction and Sequencing
DNA from sediments (PM7 and PM11) and crusts (PMW and PMB) were extracted according to the following protocol. Sediments were extracted in duplicates, while crusts were extracted in triplicates.
Cells were disrupted using FastPrep-24TM (MP Biomedicals), and genomic DNA was extracted in duplicates from each sample using the commercial kit Quick-DNA Fecal/Soil Miniprep kit (Zymo Research, Irvine, CA, United States), according to the manufacturer's instructions. DNA was quantified using SpectraMax QuickDrop (Molecular Devices, San Jose, CA, United States) and was further used as a template to investigate the composition of microbial communities in the samples. MiSeq 16S V3-V4 Metagenome Sequencing was performed by a commercial company (Macrogen, Amsterdam, Netherlands). PCR of the V3-V4 hypervariable regions of the bacterial and archeal SSU rRNA gene was performed using primers 341F (50-CCTACGGGNGGCWGCAG-3) and 805R (50-GACTACHVGGGTATCTAATCC-30), according to Illumina's 16S amplicon-based metagenomics sequencing protocol.

Metabarcoding Analysis
Sequencing primers from both forward and reverse reads and reads containing any N characters were removed using Cutadapt v2.9 (Martin, 2011). Only reads with a minimum length of 250 nt and a maximum of 301 nt were kept for further analysis. Pairedend reads for the 10 samples were processed using the DADA2 package (Callahan et al., 2016a) implemented in R by adapting existing pipelines [Callahan et al., 2016b;DADA2 Pipeline Tutorial (1.16)]. DADA2 allows accurate differentiation between sequencing errors and true biological variation, thus avoiding the use of operational taxonomic units (OTUs). Instead, DADA2 infers amplicon sequence variants (ASVs). Following primer removal, paired reads were loaded into the DADA2 pipeline and trimmed (forward reads 3 truncated at 280 nt, reverse reads 3 truncated at 250 nt), filtered (max. 2 errors per read, minimum length after trimming = 200 nt), and finally merged with a minimum required overlap of 50 nt. Chimeras were removed from merged pairs (1,576 bimeras out of 7,439; 95.8% pass rate). Following filtration and chimera removal, a total of 413,539 merged reads were retained (min per sample = 29,706, max per sample = 55,460, average = 41,353). Taxonomic classification of curated ASVs was performed using the SILVA 138.1 database (Pruesse et al., 2007). From the obtained ASVs, a mean value for triplicates (crusts) and duplicates (sediments) was used in the further analysis. Sequence data generated during this study are available in the European Nucleotide Archive (ENA) as part of project PRJEB51350 with accession numbers ERR9118894-ERR9118903.

Statistical Analysis
Differences in community composition and statistics were computed using the phyloseq (McMurdie and Holmes, 2013) package in R. Bray distance-based hierarchical clustering and alpha-diversity were estimated after merging replicates by calculating their mean values. Counts and relative abundances were further generated at the genus, class, family, and phylum levels after taxonomic agglomeration using the tax_glom function provided by phyloseq (McMurdie and Holmes, 2013). Taxa merged at each taxonomic rank were considered for abundance estimations only if they were assigned a minimum of 10 counts in at least one averaged sample. Domain Bacteria were considered for further analysis of microbial composition. Shannon and Simpson's diversity indices provided information about the composition of samples by considering both the number of species and the abundance of each species. Shannon's index gives a better description of a sample's diversity by considering the species richness and rare species, while Simpson's index considers evenness and common species. The Chao1 index estimates diversity from abundance data and gives more weight to low-abundance species.
Venn analysis was performed using Web software developed by Heberle et al. (2015). It was used to indicate the distribution of the genera abundances between the different samples. TBtools were used for heatmap representations (Chen et al., 2020).
We used multidimensional scaling (MDS) to represent the sediment samples in a bi-dimensional space. The MDS was done from a similarity matrix between the chemical characteristics of the samples to the coordinates of these in bi-dimensional space for easy visualization. The analysis was done in XLSTAT 2021.4.1 (Addinsoft, New York, NY, United States).

Chemical and Mineralogical Analysis
The chemical composition of the samples is presented in Table 1. The pH of both sediments and crusts was slightly alkaline, ranging from 8.3 to 7.4. The sediment sample collected in the Scientific Reserve (PM7) has a slightly higher pH (8.3) than the sediments collected in the Yellow Chamber (PM11, pH = 7.4) and the crusts (pH = 7.7). The carbon content was low (<3%) in the PM11 sediment and the crusts and much higher in PM7 sediment (11.6%), while the N content was low in all the samples (<1%). The MDS that uses the chemical composition of samples (Supplementary Figure 1) shows the similarity between the crusts separated from the sediment samples. Spatially and thus chemically, PM11 is the nearest to the crusts.
Calcium is the main element present in all the samples, accounting for 30-40% of the sample. The sediment sample collected in the Yellow Chamber (PM11) had higher major and trace element concentrations than the sediment in PM7. Except for Al, Fe, K, and Ba, the sediment in the Yellow Chamber had a lower concentration of metals than the crusts. The high P content in the samples from Yellow Chamber (about 6% in crusts and 2.5% in PM11 sediment) suggests the presence of minerals containing P, most probably hydroxyapatite (Supplementary Figure 2). The PM7 sediment sample has a Ca content similar to that of the crusts, but a much lower P content. The white crust (PMW) has one order of magnitude higher contents of S, Na, Mg, K, Al, and Fe and slightly higher contents of trace elements, suggesting that aluminosilicates constitute these secondary minerals. The high amount of these elements is also reflected in the different abundance of microbial composition in our samples.
X-ray diffraction patterns show that the mineralogical composition of sediments collected from the cave contains mostly primary minerals of detrital origin (Supplementary  Figure 2). The mineralogical association is relatively similar, with both samples containing quartz, muscovite, and albite in various concentrations. In addition to the mentioned silicates, the PM11 sample collected from the Yellow Chamber includes hydroxyapatite.

Guano Deposit Age
The guano sample collected from the bottom of the deposit in the Guano Chamber ( Figure 1C) yielded a 14 C age of 1,315 ± 30 years BP, which corresponds to a calibrated age between AD 654 and 775 (95.4% probability).
Few differences were identified when comparing the presence and absence of phyla, with some present only in sediments or crusts and others present only in one of the samples. Acetothermia, Campylobacterota, and Spirochaetota were only identified in PMW, while Fibrobacterota and Hydrogenedentes were observed only in PM7. In addition, Fusobacteriota was only found in crust samples, while Zixibacteria and TX1A-33 were found in sediment samples.
The relative abundance of genera was different among the samples and even between the samples of the same type. The most abundant genera in each sample ( Figure 5) were Nitrospira, Flavobacterium, and Pseudarthrobacter in PM7; wb1-P19, Nitrospira, NA ASV82, and Pseudarthrobacter in PM11; Aquicella, Delftia, and Coxiella in PMW; and wb1-P19, Pseudarthrobacter, and Massilia in PMB. The abundance of the total unassigned genera was significantly high (PM7, 31%; PM11, 31%; PMW, 29%; and PMB, 25%), and representatives were found to be among the firsts (Figure 5).
Further analysis indicated the presence of human-related bacteria mainly in Yellow Chamber samples and some presenting potential of being pathogenic (Supplementary Table 1), such as Actinomyces, Prevotella, and Bacteroides. Almost half of the identified genera were not mentioned to be present in the caves or guano before, with crust samples having the higher abundance (38% in PMW and 25% in PMB).
We performed a Venn analysis to better visualize the differences at the genus level (Figure 6). Of the 378 assigned genera, 90 were found in all four samples. The number of genera identified in PM7 (219) was lower than that identified in PM11 (228). The number of genera identified in PMW (248) was higher than in PMB (192). The crust sample PMW (72) had the highest number of unique genera, followed by the sediment samples PM7 (26) and PM11 (24), while PMB (3) had the lowest number.
The calculated alpha diversity indices (Shannon, Chao1, Simpson's, InvSimpson; Table 2) showed that the bacterial diversity differed between the samples, with a higher diversity in the sediment samples than in the crust samples. Richness estimator, Chao1, showed high values in the sediment samples.

Chemical and Mineralogical Considerations
The four analyzed samples had a high content of Ca originating from other minerals in this sample, most probably from the limestone rock or calcite dissolution. Among the four samples, only the ones collected from the Yellow Chamber (crusts and PM11) had high P content, which can be attributed to the presence of hydroxyapatite as one of their main constituents. The higher contents of S, Na, Mg, K, Al, Fe, and trace elements in the white crust suggest the presence of aluminosilicates in these secondary minerals. Detrital minerals (quartz, muscovite, and albite) are related to the weathering and transport of soils and geological formations that outcrop outside the cave (Bosch and White, 2004;Kurečić et al., 2021) by running or percolating waters. Hydroxyapatite can be formed in the cave environment by interacting with phosphorus ions and the carbonate component (Audra et al., 2019). The source of phosphate in caves is usually represented either by a fossil bone deposit or a guano accumulation (Hill, 1999;Queffelec et al., 2018). Elements, such as C, P, N, S, K, Na, Cl, Ca, Mg, Fe, Al, Zn, and Ba, were usually identified in guano deposits (Miko et al., 2001;Wurster et al., 2015;Misra et al., 2019), and high concentration of these elements was also found in both the crusts of Muierilor Cave. Considering that no bat colonies or fossil bones are present in the Scientific Reserve, the mechanisms involved in the precipitation of hydroxyapatite in the Yellow Chamber could be related to the presence of organic matter in the above level (the Touristic Passage). The possible source of organic matter in the Touristic Passage is represented by a recent guano deposit (less than c. 1,400 years BP) and fossil bone remains with ages ranging between ∼20,000 and ∼50,000 ka . The fossil bone deposits may represent the main source of PO 4 , but the leaching of guano should also be considered. The actual source of phosphate for the Scientific Reserve is challenging to assess, and further studies are required to understand the mechanisms responsible for hydroxyapatite precipitation.

Microbial Diversity of Sediments vs. Crusts
The microbiological surveys of caves indicate that the most abundant phyla usually found in caves are Proteobacteria and Actinobacteriota in different proportions, followed by Firmicutes, Bacteroidota, and Chloroflexi (Porca et al., 2012;Kieraite-Aleksandrova et al., 2015;Oliveira et al., 2017;Dhami et al., 2018;Zhu et al., 2019;Tok et al., 2021), a fact which  is consistent with our results. Proteobacteria was the most abundant phylum in all the samples and is characterized by a vast metabolic diversity (Kersters et al., 2006;Dong et al., 2020;Tok et al., 2021). Actinobacteria was the second most abundant phylum in PM7, PM11, and PMB samples. Members of this phylum are a key community in soil, where they play an important role in the decomposition of organic matter (Madigan et al., 2018;Hazarika and Thakur, 2020;Scheublin et al., 2020). The members of Acidobacteriota (second most abundant in PMW and third in PM11) are also known for their involvement in nitrogen assimilation and metabolism of iron, a finding that coincides with the high concentration of N and Fe in our samples (Kielak et al., 2016;Dedysh and Damsté, 2018). The representatives of Firmicutes (third most abundant in PM7 and PMB) are commonly found in soils, particularly in the rhizosphere (Kumar et al., 2012). Their presence in caves could be justified by the highly resistant endospores, which are generally unaffected by the environmental stress factors, such as high/low temperatures, unfavorable pH, UV, and lack of water or nutrients (Parkes and Sass, 2009). Bacterial communities found on the bat's skin represented the members of the families found in high abundance in our samples, such as members of Micrococcaceae, Flavobacteriaceae, Comamonadaceae, and Planococcaceae (González-Quiñonez et al., 2014;Lemieux-Labonté et al., 2017). In our samples, families found in high abundance were those previously mentioned in caves or related to bats. For instance, members of Micrococcaceae and Planococcaceae were sampled from cave walls (Laiz et al., 1999;Armetta et al., 2022), members of Flavobacteriaceae from cave water (Shabarova and Pernthaler, 2009), members of Nitrospiraceae from speleothem surface (Baskar et al., 2016), and members of Nitrosococcaceae from sediment and water (Zhu et al., 2019).
A slight correlation could be drawn when commenting on the relationship between chemistry and microbial composition, and the genera identified mirror the chemical composition of the samples. The Yellow Chamber samples (crusts and PM11) were higher in P, and the identified genera involved in the P cycle (Chen et al., 2006;Sarikhani et al., 2019) showed a higher abundance in these samples than those found in PM7. Such genera are Chryseobacterium, Enterobacter (PM11, PMW, and PMB), Delftia, Stenotrophomonas (PM11 and PMW), and Pantoea (PMW). Higher amounts of Fe, N, and Mn were also found in the samples collected from the Yellow Chamber when compared to the PM7 sediments. Some genera, such as Rheinheimera (PMW) and Ralstonia (PM11 and PMW), are known to be involved in the iron cycle (Swanner et al., 2011;Schröder et al., 2016), while Aphanizomenon NIES81 (PMW), Herbaspirillum, Methylocella, and Ralstonia (PM11 and PMW) are involved in the N cycle (Ploug et al., 2010;Dalsing et al., 2015;Dedysh and Dunfield, 2016;Waller et al., 2021). Bacteria involved in the Mn cycle (Yang et al., 2013;Bai et al., 2021;Lee et al., 2021) were identified as Escherichia/Shigella, Halomonas, Microbacterium (PMW), and Cupriavidus (PMW and PMB). Elements, such as N, Fe, and P, are usually found in guano (Miko et al., 2001;Wurster et al., 2015;Misra et al., 2019), P is also associated with bones (Audra et al., 2019), and Mn deposits are usually biogenic in caves (Northup and Lavoie, 2001).
Since we analyzed the chemical differences between PM7 and the Yellow Chamber samples, we could also consider this grouping when discussing human-and bat/guano-related genera. Most bat/guano-and human-related genera were identified in the Yellow Chamber samples (Supplementary Table 1). PM7 was collected from a passage without guano or other organic accumulation.
Bats are a significant reservoir of pathogens and a well-known vector for disease transmission (Ogórek et al., 2018;Agustin et al., 2019). Genera found in crusts, such as Fusobacterium and Rothia, were previously related to bats, while Microbacterium was previously related to guano. Some genera identified in our samples, most of them in crusts, are also found in the guano of other caves in Europe, which are habitats for the same bat species. These are Enterococcus, Acinetobacter, Pseudomonas, Paenibacillus, Bacillus, Staphylococcus, Rahnella, Micrococcus, Enterobacter, Lysinibacillus, and Sphingobacerium (Tomova et al., 2013;Vandžurová et al., 2013;Veikkolainen et al., 2014;Vengust et al., 2018;Wolkers-Rooijackers et al., 2019;Dimkić et al., 2020;Gerbáčová et al., 2020). Their abundance did not dominate the microbial composition in our samples, but their presence could indicate a possible connection to the upper Guano Chamber.
The presence of a significant number of genera never reported in a cave before (Supplementary Table 1) suggests the limited knowledge regarding the bacterial diversity of caves. Almost half of the genera identified in our samples were never reported in a cave or were never related to bats/guano. Some genera are known members of human microbiota: Abiotrophia (Rasic et al., 2020), Gemella (García López and Martín-Galiano, 2020), Granulicatella (Shailaja et al., 2013), Treponema (Seña et al., 2015), and Haemophilus (Johnson, 2018). Genera, such as Bacillus, Escherichia, and Staphylococcus, which are known as "human indicator bacteria, " are usually found in high concentrations after extended cave visits (Lavoie and Northup, 2006;Bercea et al., 2018;Mudgil et al., 2018) and were identified in our samples in high abundance. Recent studies in Muierilor Cave compared the microbial composition in touristic and nontouristic sites. Based on an air survey and human exposed surfaces (Bercea et al., 2018(Bercea et al., , 2019, the touristic sector of Muierilor Cave fits in the high-risk class (500-2,000 CFU/m 3 , according to the European Commission's report on Biological Particles in Indoor Environments and Commission of the European Communities, 1994). Moldovan et al. (2020) concluded a clear difference in the microbial composition of water habitats between touristic and non-touristic sites, indicating the strong impact of tourism. These studies found a high abundance of members of Bacillus and Staphylococcus, which are also abundant in our samples.
Muierilor Cave, as a touristic cave, may be considered as a reservoir of novel and allochthonous species, since animals and humans transit and use it for different purposes. Moreover, visitors play an important role in the spread of bacteria and continuously influence the composition of the microbial community.

CONCLUSION
Most of the identified ASVs in four samples, two sediments and two crusts, collected from Muierilor Cave belong to bacteria. The crusts and the nearby sediment sample in a small chamber (the Yellow Chamber) of the cave are phosphate-rich deposits, with abundant bacteria involved in the P cycle. Moreover, bacteria involved in Fe, Mn, and N cycles were found, as these elements are commonly identified in high concentrations in guano. Since no bat colonies or fossil bones were present in the chamber, a high concentration of these elements could be sourced from organic deposits (fossil bones and guano) located at the upper levels of the cave. The high diversity of human-related genera identified only in the crusts and sediment samples collected from the Yellow Chamber, such as Capnocytophaga, Anaerococcus, Abiotrophia, Actinomyces, Alloprevotella, and Eikenella, could only be related to the upper touristic sector of the cave. Tourism in caves with bat colonies and guano accumulation should benefit from continuous monitoring of the air and floor microbiomes.
Approximately 47% of non-identified bacteria in a small sector of the complex subterranean system of Muierilor Cave indicates the limited knowledge of the bacterial diversity in caves, and these bacteria have high potential in human health and biotechnology applications.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: ENA -PRJEB51350.

AUTHOR CONTRIBUTIONS
OM and CH designed the research. OM, RN-B, CH, IM, and LF collected the field data. RN-B, CH, and PB contributed to the extraction, bioinformatic analysis, and interpretation of the molecular data. EL, OC, LF, VF, and IA contributed to the chemical and mineralogical analyses. SC and IM contributed to the geological context and radiocarbon data interpretation. CH, OM, PB, EL, IM, and LF wrote the first draft. All the authors corrected and approved the final version.

ACKNOWLEDGMENTS
We are grateful to the editor and the reviewers for their valuable suggestions that improved the manuscript clarity considerably.