Gypsum Cave Biofilm Communities are Strongly Influenced by Bat- And Arthropod-Related Fungi
Caves are populated by a wide variety of microorganisms, including bacteria, fungi, algae, protozoa, and more, which thrive in both pristine and show caves. These microorganisms are found to be associated with biofilms that cover cave walls and sediments. Most caves present conspicuous biofilms, observable with the naked eye, in the form of round colonies, patches, or extensive microbial mats of different colors, generally white, yellow, and grey . Furthermore, show caves subjected to artificial lighting exhibit wide phototrophic biofilms mainly composed of cyanobacteria and microalgae . The interest in the study of biofilms lies in their contribution to the physical and chemical deterioration of rocks and speleothems, which is of particular importance for the conservation of rock art that house many caves . The case of Altamira Cave, Spain, is paradigmatic with respect to the colonization of its walls by different colored biofilms . Jurado et al.  studied both biofilms and walls in this cave without visible microbial colonizations, detecting only bacterial DNA and no fungal DNA. The absence of fungi was further established by scanning electron microscopy . This was later confirmed by González-Riancho Fernández . This absence was attributed to the production of bioactive compounds with antifungal properties by bacteria that form the biofilms. In particular, many strains of Crossiella, Streptomyces, Micromonospora, Nocardia, Pseudonocardia, Lechevalieria, Kibdelosporangium, Amycolatopsis, Saccharothrix, among others, were isolated from these biofilms . For example, recently, two Crossiella strains isolated from Altamira white biofilms inhibited the growth of other bacteria and fungi. Their genomes showed an exclusive combination of gene clusters involved in the synthesis of lanthipeptides, lasso peptides, sactipeptides, furans, polyketide synthases, and non-ribosomal peptide synthetases . Further studies corroborated the absence of fungi within cave biofilms. Martin-Pozas et al.  investigated yellow biofilms in Pindal Cave, Spain, applying specific CARD-FISH probes targeting Ascomycota and Basidiomycota, with negative results. Additionally, sequence analysis using fungal primers did not detect any fungal sequences, concluding that fungi were not present in the yellow biofilms. However, other authors found fungi in cave biofilms from different countries . It should be noted that fungi are common in cave air, as well as in excrements, guano, and animal cadavers frequently found in caves and sediments . Many cave fungi have been related to arthropods and bats , while others appear to be specific to cave ecosystems . Although considerable attention has been directed towards fungi in limestone and volcanic caves , studies on gypsum caves are scarce, if not lacking. Cunningham et al.  found in Lechuguilla, a gypsum cave in New Mexico, patches of fungi on deposits, and aspergilli associated with Fe-, Mn- and S-rich encrustations. Cacchio et al.  investigated the involvement of bacteria in the formation of speleothems in the gypsum cave Grave Grubbo, Italy. Mihajlovski et al.  studied the fungal diversity in gypsum efflorescences in a French cave and described the abundant presence of Isaria fumosorosea and Engyodontium album , two well-known entomopathogens . Some gypsum karsts are important in Western Europe because they are excellent examples (geoheritage sites) of emerged outcrops of selenitic gypsum, which formerly corresponded to different peripheral basins around the Mediterranean with a diachronous precipitation of evaporite deposits during the Messinian age. These karstified gypsum outcrops are mainly located in the Apennines and Sicily, Italy , while in Spain the Gypsum Karst of Sorbas in Almeria is notable . The geology, geomorphology, hydrogeology, and microclimatology of this last karst have been extensively studied in the last two decades ; however, the microbial communities of their gypsum caves are unknown. In the gypsum karst of Sorbas, Covadura Cave harbored an abundant bat population, although this has decreased considerably in recent years, with only isolated individuals currently detected in the galleries sampled in this study. In C3 Cave there is an absence of roosting bats. Altamira and Pindal caves, located in northern Spain, do not have bats in their confined and close environments. Nonetheless, all the four caves show abundant white and yellow biofilms composed of bacteria, but no fungi were retrieved from the Altamira and Pindal caves , and unpublished data. Here, we study the fungal diversity in the biofilms that cover the walls of two gypsum caves to determine whether bat microbiomes indeed have influence on the wall biofilms. 

The Gypsum Karst of Sorbas The Gypsum Karst of Sorbas, is located in the province of Almeria, southeast Spain (Fig. ). This is formed by significant gypsiferous Messinian evaporates within a 120 m thick cyclic sequence formed by alternating gypsum and marly units, and many caves amounting to 100 km of passages, including the six most relevant caves for their extension and representativeness in the local karst hydrogeology (Covadura Cave, C3 Cave, Gypsum Cave, Water Cave, Treasure Cave and GEP Complex) . The Covadura Cave is 120 m deep and 275 m long in the studied area, although the system has more than 4 km of galleries. The C3 Cave, at some 200 m from the Covadura Cave, is only about 3 m deep and has a length of 150 m . A diversity of gypsum speleothems, such as stalactites, coralloids, gypsum crusts, etc. have been described in these caves, as well as colored biofilms (Fig. ). The Covadura Cave harbored a bat population of more than 5,000 bats, composed of the species Rhinolophus ferrumequinum, Rhinolophus hipposideros , and Miniopterus schreibersii . Rhinolophus euryale has occasionally been found . Access to the caves is forbidden from October to March to protect the bats during the hibernation period. In Covadura Cave there has been some activation of the watercourse and partial flooding washed away guano; however, the area of wall biofilms has not been affected and may have retained previously existing fungal populations. However, through visits and sampling carried out in recent years, the disappearance of guano in the cave has been confirmed. In Cave C3, located about 200 m from Covadura Cave, there have never been bat colonies. The climate in the Sorbas area is semi-arid Mediterranean, with a mean annual temperature and rainfall of 17 °C and 274 mm, respectively. Estimated annual potential evapotranspiration is nearly five times the mean annual precipitation, and about 80% annual rainfall during low-frequency rainstorm events, usually in the autumn. The vegetation comprises scattered shrubs and dwarf shrubs, as well as some crass perennial plants are also common. The cover of perennial vegetation is below 40 % at both study sites and is mainly made up of tussock grass. Plant interspaces are usually colonized by well-developed biocrusts dominated by lichens. Preparation and Characterization of Samples by Scanning Electron Microscopy (SEM) Biofilms were collected and transported at 4°C without separating them from the mineral support on which they were growing. Sample preparation for SEM analysis was carried out as described by Martin-Pozas . Biofilm structures were examined under high and low vacuum conditions using an environmental scanning electron microscope Inspect (FEI, USA). Sampling Method In 2010 and 2022, nine white and yellow biofilm samples on the gypsum walls were collected in Covadura Cave, while only two biofilm samples from the C3 Cave were sampled in 2022 (Supplementary Figure ). Biofilm samples were collected using sterile scalpels and preserved in Lifeguard preservation solution (Qiagen, Hilden, Germany) until their arrival in the laboratory, where they were stored at -80°C. DNA Extraction, Sequencing and Data Processing Genomic DNA extraction was performed using the FastDNA SPIN Kit for Soil (MP Biomedicals, Illkirch, France). DNA concentrations were quantified using a Qubit 2.0 fluorometer (Invitrogen, Carlsbad, CA, USA). High-throughput sequencing of the extracted DNA was carried out by FISABIO Institute (Valencia). The intergenic transcribed spacers (ITS) of fungal organisms were amplified according to Toju et al. . The primer sequences used were ITS3_KYO2 (GATGAAGAACGYAGYRAA) and ITS4_KYO1 (TCCTCCGCTTWTTGWTWTGC), using Illumina MiSeq and 2 × 300 paired end sequencing. Sequence data have been performed using the Qiime2 pipeline . Metataxonomy analysis was performed using some of the qiime2 plugins. Denoising, joining of the paired-ends, and depletion of the chimera was performed starting from data from the paired ends using DADA2 pipeline . Taxonomic affiliations were assigned using the Naive Bayesian classifier integrated in the quiime2 plugins. The database used for taxonomic assignment was UNITE version 10.0 . All computations and statistics were performed within the R Statistics environment . Taxonomic data obtained in QIIME2 were imported into R and analyzed with the phyloseq package. From these taxonomic data, ecological roles were assigned using the FUNGuildR tool and the FUNGuild database . To assess ecological alpha diversity within each sample, we used non-phylogenetic diversity indices such as Chao1, Shannon, and Simpson. These common metrics were calculated using raw species-level data with R packages phyloseq and vegan. The plots and heatmap were built using the R packages ggplot 2 and pheatmap. Deposit of Sequences The raw reads were deposited into the NCBI Sequence Read Archive (SRA) database under the project number PRJNA1103123. 

The Gypsum Karst of Sorbas, is located in the province of Almeria, southeast Spain (Fig. ). This is formed by significant gypsiferous Messinian evaporates within a 120 m thick cyclic sequence formed by alternating gypsum and marly units, and many caves amounting to 100 km of passages, including the six most relevant caves for their extension and representativeness in the local karst hydrogeology (Covadura Cave, C3 Cave, Gypsum Cave, Water Cave, Treasure Cave and GEP Complex) . The Covadura Cave is 120 m deep and 275 m long in the studied area, although the system has more than 4 km of galleries. The C3 Cave, at some 200 m from the Covadura Cave, is only about 3 m deep and has a length of 150 m . A diversity of gypsum speleothems, such as stalactites, coralloids, gypsum crusts, etc. have been described in these caves, as well as colored biofilms (Fig. ). The Covadura Cave harbored a bat population of more than 5,000 bats, composed of the species Rhinolophus ferrumequinum, Rhinolophus hipposideros , and Miniopterus schreibersii . Rhinolophus euryale has occasionally been found . Access to the caves is forbidden from October to March to protect the bats during the hibernation period. In Covadura Cave there has been some activation of the watercourse and partial flooding washed away guano; however, the area of wall biofilms has not been affected and may have retained previously existing fungal populations. However, through visits and sampling carried out in recent years, the disappearance of guano in the cave has been confirmed. In Cave C3, located about 200 m from Covadura Cave, there have never been bat colonies. The climate in the Sorbas area is semi-arid Mediterranean, with a mean annual temperature and rainfall of 17 °C and 274 mm, respectively. Estimated annual potential evapotranspiration is nearly five times the mean annual precipitation, and about 80% annual rainfall during low-frequency rainstorm events, usually in the autumn. The vegetation comprises scattered shrubs and dwarf shrubs, as well as some crass perennial plants are also common. The cover of perennial vegetation is below 40 % at both study sites and is mainly made up of tussock grass. Plant interspaces are usually colonized by well-developed biocrusts dominated by lichens. 

Biofilms were collected and transported at 4°C without separating them from the mineral support on which they were growing. Sample preparation for SEM analysis was carried out as described by Martin-Pozas . Biofilm structures were examined under high and low vacuum conditions using an environmental scanning electron microscope Inspect (FEI, USA). 

In 2010 and 2022, nine white and yellow biofilm samples on the gypsum walls were collected in Covadura Cave, while only two biofilm samples from the C3 Cave were sampled in 2022 (Supplementary Figure ). Biofilm samples were collected using sterile scalpels and preserved in Lifeguard preservation solution (Qiagen, Hilden, Germany) until their arrival in the laboratory, where they were stored at -80°C. 

Genomic DNA extraction was performed using the FastDNA SPIN Kit for Soil (MP Biomedicals, Illkirch, France). DNA concentrations were quantified using a Qubit 2.0 fluorometer (Invitrogen, Carlsbad, CA, USA). High-throughput sequencing of the extracted DNA was carried out by FISABIO Institute (Valencia). The intergenic transcribed spacers (ITS) of fungal organisms were amplified according to Toju et al. . The primer sequences used were ITS3_KYO2 (GATGAAGAACGYAGYRAA) and ITS4_KYO1 (TCCTCCGCTTWTTGWTWTGC), using Illumina MiSeq and 2 × 300 paired end sequencing. Sequence data have been performed using the Qiime2 pipeline . Metataxonomy analysis was performed using some of the qiime2 plugins. Denoising, joining of the paired-ends, and depletion of the chimera was performed starting from data from the paired ends using DADA2 pipeline . Taxonomic affiliations were assigned using the Naive Bayesian classifier integrated in the quiime2 plugins. The database used for taxonomic assignment was UNITE version 10.0 . All computations and statistics were performed within the R Statistics environment . Taxonomic data obtained in QIIME2 were imported into R and analyzed with the phyloseq package. From these taxonomic data, ecological roles were assigned using the FUNGuildR tool and the FUNGuild database . To assess ecological alpha diversity within each sample, we used non-phylogenetic diversity indices such as Chao1, Shannon, and Simpson. These common metrics were calculated using raw species-level data with R packages phyloseq and vegan. The plots and heatmap were built using the R packages ggplot 2 and pheatmap. 

The raw reads were deposited into the NCBI Sequence Read Archive (SRA) database under the project number PRJNA1103123. 

Biofilms are mainly composed of bacteria, but fungal spores were also evident (Fig. ). SEM images of these biofilms were similar to those obtained in the biofilms from Castañar Cave, with a high abundance of bacteria and fungi . The bacterial composition of the biofilms will be discussed in a further work, but we can anticipate that white and yellow biofilms in these gypsum caves are mainly composed of bacteria with a relatively high abundance of two main phyla: Actinobacteriota and Pseudomonadota , followed by Acidobacteriota and Planctomycetota . In the study of fungi, a total of 2,515,260 high-quality sequence reads were obtained, ranging from 8,652 to 309,474 from the white and yellow biofilms (Supplementary Table ). In total, 422 amplicon sequence variants (ASVs) were assigned based on 99% sequence similarity. According to diversity indices, samples with the highest Chao1, Shannon, and Simpson values are found in the white and yellow biofilms from Covadura Cave, indicating greater species diversity and more uniform distribution of relative species abundance compared to samples from Cave C3 (Supplementary Table and Figure ). In summary, Covadura Cave exhibits greater species richness and diversity than C3 Cave, and specifically, the yellow biofilms exhibit higher species diversity, particularly those from 2010. In the white biofilms collected in the Covadura Cave, Ascomycota (relative abundance 62.2% – 99.1%) and Basidiomycota (0.4% – 37.7%) were the most abundant phyla. However, yellow biofilms showed a different and variable pattern, with higher abundances of Mortierellomycota in three samples (19.2% – 97.4%), Ascomycota in four samples (39.8% – 97.0%) and Basidiomycota in the other two samples (11.0% – 20.4%). Additionally, a yellow biofilm sample presented 19.0% relative abundance of a phylum at an unclassified taxonomic rank. Negligible or no relative abundances were observed for other phyla, including unidentified sequences (Fig. ). In general, a decrease in Ascomycota relative abundances was detected in white biofilms in 2022 (75.6 – 62.2%) compared to levels recorded in 2010 (77.3 – 99.1%), although no consistent pattern was observed in yellow biofilms, with variable relative abundances between samples. This lack of consistency was also evident in the abundances of Basidiomycota and Mortierellomycota . In C3 Cave, the fungal community in the white biofilm was predominantly composed of Ascomycota , representing 95.9% of the relative abundance, followed by Mortierellomycota (3.8%). The yellow biofilm exhibited a reverse trend, with Mortierellomycota comprising 94.8% of the community and Ascomycota representing around 5.0%. In all samples, the most abundant classes were Sordariomycetes, Mortierellomycetes , and Dothideomycetes. Sordariomycetes were represented in all samples from the two caves, although with variable abundance depending on the samples (1.8% – 95.4%) (Fig. ). Mortierellomycetes were relatively abundant in the yellow biofilms of Covadura Cave (0.2% – 97.4%), and in C3 Cave (94.8%), but only appeared in a sample of white biofilms from Covadura Cave (1.6%) and in C3 (3.8%). Dothideomycetes were more abundant in the yellow biofilms collected in the Covadura Cave in 2010 (9.6% – 75.5%) than in 2022 (0% – 11.1%) as well as for the white biofilms (6.0% – 57.2% in 2010 vs 0.3% – 0.5% in 2022). An unidentified Basidiomycota taxonomic class presented relative abundances of 21.7% and 37.1% in the two white biofilms from Covadura Cave collected in 2022. Eurotiomycetes and Agaricomycetes were relatively abundant in only one sample of Covadura Cave white biofilms collected in 2010 (16.1% and 13.0%, respectively) . In the white biofilms samples collected in 2022 in the Covadura Cave, Eurotiomycetes were absent and Agaricomycetes showed low abundances (0.4 – 0.6%). In the yellow biofilms of Covadura Cave Eurotiomycetes reached low abundances in 2022 (0% – 1.3%) and in 2010 (0.3% – 2.4%) and Agaricomycetes slightly higher in 2022 (0% – 8.6%) than in 2010 (1.4% – 3.5%). In the C3 Cave, the abundances of these two classes were negligible. Relative abundances for other classes were missing or negligible, except for the 2010 white ( Tremellomycetes 9.8%) and yellow biofilm samples ( Leotiomycetes 5.8%, Saccharomycetes 5.7% and Malasseziomycetes 9.2%). More than 100 fungal genera were identified, with an additional forty-six genera remaining unidentified in the Covadura and C3 caves. Significant differences were observed between white and yellow biofilms on the abundance of the genera (Fig. ). In Covadura Cave, several genera presented relative abundances above 10% in different biofilms, but the distribution between white and yellow biofilms was different. The 2010 white biofilms were characterized by the abundant occurrence (>10%) of Engyodontium, Toxicocladosporium, Lecanicillium, Neocamarosporium and Trechispora , and in the 2022 white biofilms of Lecanicillium, Lasionectria and an unidentified Basidiomycota. The yellow biofilms showed a high abundance of Lecanicillium, Cladosporium, Podila, and Mortierella in the 2010 samples, and of Lecanicillium and Podila in the 2022 samples. Other genera also well represented (5-10%) included Alternaria, Malbranchea and Chrysosporium in white biofilms and Malassezia and Hirsutella in yellow biofilms. In the C3 Cave, only Lecanicillium and Verticillium were important (>10%) in the white biofilms and Mortierella in the yellow ones. Metacordyceps (7.2%) was identified in the white biofilms. The most abundant and common genus in all samples was Lecanicillium , an insect pathogen, both in Covadura white and yellow biofilms, as well as in C3 white and yellow biofilms. Lecanicillium antillanum was relatively abundant (23.2%) in a yellow biofilm from 2010. Other well-known entomopathogenic fungal genera are Engyodontium (=Beauveria), Hirsutella, Simplicillium , Metacordyceps , etc., also found in Covadura biofilms. Engyodontium was only abundant (74.1%) in a white biofilm collected in 2010, while Hirsutella (5.1%) and Simplicillium (2.6%) were found in yellow biofilms of the same year. The Mortierellaceae include the genus Mortierella , Linnemannia, and Podila among others. M. alpina reached a relative abundance of 94.6% in the yellow biofilm of the C3 Cave and 17.0% in a 2010 yellow biofilm from the Covadura Cave. M. elongata was identified in a 2010 yellow biofilm (4.9%). The genus Podila was established in 2020 to accommodate some Mortierella species that are difficult to identify by ITS sequences . Podila was represented in yellow biofilms with 97.3% relative abundance in a yellow biofilm from 2022 and 27.5% in another from 2010. The genus Cladosporium was relatively abundant in three 2010 yellow Covadura biofilms (5.9 –64.6%), less abundant in 2010 white biofilms (4.5 – 7.9%), and largely absent in most 2022 samples. Its abundance in the C3 Cave was also insignificant (0.0 – 0.4%), Cladosporium halotolerans was only identified in a yellow biofilm from 2010 (1.0%). The species Verticillium leptobactrum (= Leptobacillium leptobactrum ), a recognized entomopathogenic fungus, was identified in the white biofilm of C3 Cave (10.0% relative abundance) and 0.1% in a few Covadura and C3 biofilms. Metacordyceps reached 7.2% in the white biofilm of C3 Cave and 0.1% in a yellow Covadura biofilm. Abundant genera, but only present in a unique sample of a Covadura 2010 white biofilm, were Toxicocladosporium , and Neocamarosporium . These two genera and a few others such as Catenulostroma, Libertasomyces, Xylodon, Symmetrospora, and Thelephora , were recovered in low abundances in the biofilms. A detailed classification of fungal traits is presented in Fig. . There, it can be observed that the abundances of saprotrophs and pathotrophs were higher than 60% of the total relative abundances in C3, while the abundances of saprotrophs were also high in Covadura Cave, but the pathotrophs decreased to 40%. Other important traits were wood saprotroph and animal pathogen. 

Many studies reported the presence of fungi in caves, particularly in phototrophic biofilms. Predominantly, these fungi belong to the phyla Ascomycota, Zygomycota and Basidiomycota . The dominant fungal phyla in the Covadura and C3 gypsum caves were Ascomycota, Mortierellomycota , and Basidiomycota. This is in agreement with data from other authors in karstic caves . The most abundant classes in the two caves were Sordariomycetes, Mortierellomycetes , and Dothideomycetes. Similar abundances of these three classes have been reported in a few caves from different countries and rock surfaces . The Sordariomycetes are represented by the genera Lecanicillium , Engyodontium, and Lasionectria, with abundances over 10% in at least one of the samples. The genus Lecanicillium comprises insect pathogens and also parasitizes other arthropods, fungi, and plants . Different species of Lecanicillium were reported in Spanish and other European caves with an abundance of arthropods . However, only one species of Lecanicillium was identified, L. antillanum in the 2010 Covadura white and yellow biofilms and the C3 yellow biofilm. The caves in the Gypsum Karst of Sorbas harbor a wide diversity of arthropods, along with bats. Indeed, many authors have recorded the occurrence of numerous arthropod taxa in these caves, including members of the orders Araneida, Anactinotrichida , Opilionida , Pseudoescorpionida, Isopoda, Quilopoda, Collembola, Psocoptera, Thysanoptera, Thysanura, Hemiptera, Coleoptera, Siphonaptera, Diptera, Diplura, Orthoptera , Trichoptera , among others . Most of these orders and others have been reported in bat guano in different geographical locations . Therefore, the abundant occurrence of Lecanicillium is probably linked to cave arthropods and bat guano. The Mortierellomycetes , (Mortierella , Linnemannia, Podila , etc.), are saprobic organisms with a widespread distribution in a wide range of habitats . Among all Mortierella species, the most commonly recovered in caves were M. alpina and M. elongata (= Linnemannia elongata ) . Many species of the genus Mortierella have been isolated from bats, dung samples collected in caves and mines, wings and bat carcasses, in addition to soils where they are abundant . Podila species are frequent in soils, compost and dung . In addition to Mortierella, a large number of genera found in the biofilms of Covadura Cave have previously been recovered from bat guano. These include Lecanicillium, Cladosporium, Alternaria, Aureobasidium, Botrytis, Simplicillium, Debaryomyces, Talaromyces, Chrysosporium , Apiotrichum , and Cutaneotrichosporon , while only Lecanicillium, Cladosporium, Alternaria , and Aureobasidium were found in C3 Cave . Furthermore, most of these genera and others ( Filobasidium, Malassezia, Vishniacozyma, Wallemia, Comoclasthris, Cephalotrichum, Lasionectria, Solicozyma, Acremonium, etc.) were isolated from bat mycobiomes . The Dothideomycetes are represented by the genera Cladosporium , Toxicocladosporium , and Neocamarosporium with significant abundances. Cladosporium halotolerans, the most frequently isolated indoor species , was found in the two gypsum caves, and previously in other Spanish caves . Other species of Cladosporium are ubiquitous plant pathogens and were also isolated from caves . Toxicocladosporium was detected in biofilms from stones , and plants , as well as Neocamarosporium . Important ecological traits observed in the genera found in the Covadura and C3 biofilms were the chitinolytic and keratinolytic activities. Chitinolytic activity is related to the occurrence of insectivorous bats and chitin remains in the guano. Chitinolytic fungal genera were identified in both caves but with a higher diversity in Covadura ( Cladosporium, Alternaria, Verticillium, Aureobasidium, Acremonium, Trichoderma, Talaromyces, Botrytis, Stachybotrys , and Mortierella than in C3 ( Verticillium, Aureobasidium, Acremonium, Mortierella ). Keratinolysis is common in members of the order Onygenales ( Malbranchea, Chrysosporium ) and Eurotiales ( Cladosporium , Talaromyces ). Another fungal group of interest in Covadura Cave is the class Tremellomycetes that includes the basidiomycetous yeasts Filobasidium, Apiotrichum, Vishniacozyma, Cutaneotrichosporon, and Dioszegia . These genera have been associated with bats . In the C3 Cave only Vishniacozyma (= Cryptococcus ) was identified. Members of this genus are common in soils and plants . Of no less interest is the increased abundance of a few genera ( Lecanicillium, Lasionectria, unidentified Basidiomycota ) in the Covadura 2022 sampling of white biofilms and of Lecanicillium and Podila in the yellow biofilms. The decrease in abundance in 2022 of a few genera that were abundant in the 2010 yellow biofilms ( Engyodontium, Mortierella ) is noteworthy. These changes can be attributed to fluctuations in arthropod/bat populations and/or severe droughts in recent years. When comparing the Covadura and C3 cave biofilms with those of the Altamira and Pindal limestone caves, the strongest difference was the absence of fungi in the limestone biofilms . In the limestone caves no bat roosts were found. The entrances to the gypsum caves are open, allowing the entry and roosting of bats. Therefore, the high number of fungal genera identified in the bat microbiome by a few authors , and their finding in cave biofilms, are not surprising. Furthermore, Hathaway et al.  reported that bat microbiomes have an influence on the cave wall microbiome, and Kokurewicz et al.  stated that the number of bats was a key factor determining the occurrence of fungal spores at the hibernation site. There is a general consensus on the influence of arthropods, bats and guano on cave fungi. Yoder et al.  reported the presence of entomopathogenic fungi associated to cave spiders. Vanderwolf et al.  studied the fungi in three caves, one populated by bats and the other two without bats. They concluded that caves with active bat hibernacula supported bat-related fungi, while caves without bats, but abundant arthropods, contained entomopathogenic fungal genera. The conclusions of this and other papers  were confirmed in this study on biofilms from gypsum caves. In summary, the fungi recovered from the different white and yellow biofilms in Covadura Cave showed a wide diversity, with differences depending on their location and year of sampling, and were highly influenced by the presence of bats, guano, and arthropods that thrive in the guano. In the C3 Cave the occurrence of bat-related fungi was low, probably due to the absence of roosting bats, with the cave environment dominated by soil- and arthropod-related fungi. 

Supplementary Material 1: (PDF 433 kb) 

