Microbial Communities in Vermiculation Deposits from an Alpine Cave

Morgana Cave is located in Val di Scerscen, Central Italian Alps. The cave opens at an altitude of 2,600 m a.s.l. close to the retreating glacier Vedretto di Scerscen, and its entrance was discovered 30 years ago hidden underneath the glacier. A characteristic of this cave is the occurrence of vermiculation deposits on the walls and ceiling. In general, the composition of the microbial communities in cave vermiculations is relatively unknown and rarely investigated. Here we present the data of a geomicrobiological study of vermiculations from an Alpine cave subjected to extreme climate conditions. The microbial communities were dominated by 13 main phyla of Bacteria, and contained a negligible percentage (<1%) of Archaea. The two major bacterial classes were Gammaproteobacteria and Betaproteobacteria, whose metabolic traits were mainly associated with the nitrogen cycle. In addition, psychrophilic and methanotrophic bacterial groups were identified. The occurrence of a large number of uncultured members, at the lowest taxonomic ranks, indicated the presence of still unexplored microbial taxa in the vermiculations.


INTRODUCTION
Global warming is causing a rapid glacial retreat. The total glacier-covered surface area in the Alps decreased from an estimated 4,500 km 2 around 1850 to around 2,200 km 2 in 2001, 1,800 km 2 in 2010 (Haeberli et al., 2013) and 1,650 km 2 in 2015 (Sommer et al., 2020).
The glacier Vedretta di Scerscen is located in Val di Scerscen, Central Italian Alps (Northern Lombardy), in the scenic setting of Valmalenco ( Figures 1A,G). The glacier retreat has uncovered a few caves, the most important of which, opening at about 2,600 m a.s.l., are: Veronica Cave (638 m long), Morgana Cave (770 m), and Tana dei Marsooi (77 m). These caves represent very interesting systems from both pure geological-geomorphological and microbiological standpoints (De Waele et al., 2014).
Morgana Cave was hidden underneath the rapidly retreating glacier until the end of the last century ( Figure 1B). In the Morgana Cave area, snow and ice remain until the beginning of June. At the sampling time (July 17th, 2017) the entrance of the cave was open, but there was still snow around and above the cave. The cave had very low temperatures, slightly above 0°C, with freezing occurring close to the entrance, and rather stable conditions occurring deeper into the cave (reaching 2°C 50 m inward).
Morgana Cave is characterized by phreatic morphologies (Figures 1C,D) and contains glacial sediments ( Figure 1E). The cave probably formed before the last ice age. The entrance area of Morgana cave and many smaller nearby caves were all apparently cut both by valley slope retreat and ice scouring. The caves are characterized by shallow vadose entrenchments ( Figure 1D) which record a later phase of free-running meltwaters in the cave passages. Air circulation triggered by the temperature difference between the cave interior and exterior induces condensationcorrosion processes, which are actually shaping the cave wall and ceiling morphologies ( Figure 1F), with smooth reliefs, notches and cupola (D'Angeli et al., 2018).
Cave vermiculations are thin, irregular and discontinuous deposits of incoherent particles. Generally, they exhibit several kinds of morphologies (i.e. dots, dendritic, hieroglyphic), colors (red, brown, grey or white) and dimensions (De Waele et al., 2014;Addesso et al., 2020). According to Bini et al. (1978) vermiculations are mainly composed of clays and are the Opportunities to obtain samples are rare, and the samples are small and hard to collect due to their mineral composition. The small sample size sometimes makes it difficult to perform a complete microbiological analysis.
Vermiculations, being aggregates of weathered minerals and clays, represent an interesting and suitable niche for microorganisms. In addition, organic compounds are largely adsorbed onto the clays. Clays also increase water retention with respect to the host rock, facilitating colonization and thriving of microorganisms. Therefore, vermiculation deposits are ideal sites for studying the diversity of cold-adapted microbial communities in Morgana Cave.
Here we report the microbial community of different vermiculation deposits (Figure 2) collected on the walls of Morgana Cave, in an attempt to expand the understanding on the geomicrobiology of Alpine caves.
The samples, stored in sterile containers, were maintained at 5°C until laboratory arrival and those used for field emission scanning electron microscopy (FESEM) were fixed in situ. The samples for molecular biology were preserved in RNAlater and then conserved at −80°C in the laboratory.
Vermiculations showed different colors (whitish, greyish, brownish); the whitish are dominated by dolomite clasts, whereas clay minerals have been solely observed in the greyish and brownish deposits.

Field Emission Scanning Electron Microscopy
The morphology of the different vermiculations was studied by field emission scanning electron microscopy (FESEM) combined with energy dispersive X-ray spectroscopy (EDS). Previously, samples were dried at 50°C for 24 h and directly mounted on sample stubs, sputter coated with gold, and subsequently examined in a FEI Teneo FESEM equipped with an EDS detector. FESEM examinations were operated using the secondary electron (SE) detection mode (Everhart-Thomley) with an acceleration voltage of 5 kV for ultra-high resolution images and 10 kV for elemental microanalysis.

Nucleic Acid Extraction
Genomic DNA and RNA were extracted from 300 mg of samples using the DNeasy PowerBiofilm DNA and RNA isolation kits according to manufacturer's protocols (Qiagen) and quantified using a Qubit 2.0 fluorometer (Invitrogen). DNA and RNA were stored at −80°C until their use.

Next-Generation Sequencing and Data Analysis
The extracted DNA of all samples (except M1) was analyzed by next generation sequencing (NGS). We investigated the bacterial V3 and V4 regions of the 16S SSU rRNA gene using Illumina MiSeq and 2 × 250 paired end sequencing, according to the Illumina 16S metagenomic library preparation protocol used by Macrogen (Seoul, Korea).
The primer sequences used in this study were 341F (CCTACGGGNGGCWGCAG) and 805R (GACTACHVGGGTATCTAATCC). Takahashi et al. (2014) showed that this prokaryotic universal primer for simultaneous analysis had a coverage rate of 98% for Bacteria and 94.6% for Archaea. Raw data were processed in QIIME 1.9.1 (Caporaso et al., 2010). Quality control and trimming were performed using FASTQC (http://www. bioinformatics.babraham.ac.uk/projects/fastqc/) and Trimmomatic (0.36 version), respectively. Paired end reads were assembled using PEAR (Zhang et al., 2014). Primers removal was carried out using a PERL script developed by the Environmental Microbiology group (IRNAS-CSIC). Operational Taxonomic Units (OTUs) were clustered at 97% cutoff using UCLUST (Edgar, 2010). The SILVA database (version 132) was used to assign the taxonomic classification of each 16S rRNA gene sequence with a threshold of 80% (Quast et al., 2013). The alpha_diversity. py command was employed for alpha diversity measurements.
Phylogenetic trees were constructed using the maximumlikelihood (Felsenstein, 1981) method in MEGA X (Kumar et al., 2018). A bootstrap analysis of 1000 re-samplings was used to evaluate the trees' robustness. Frontiers in Earth Science | www.frontiersin.org December 2020 | Volume 8 | Article 586248

Deposit of Sequences
Raw sequences were deposited in the NCBI SRA database under accession numbers ERS3506978-ERS3506982. Clone sequences (≈360 bp for nifH gene, 650 bp form archaeal amoA gene and 490 bp for bacterial amoA gene) were deposited at GenBank under the accession numbers LR596555-LR5965594, LR596553, and LR778119-LR778132.

Field Emission Scanning Electron Microscopy
The study of the samples by field emission scanning electron microscopy (FESEM) revealed the presence of filamentous microorganisms in the vermiculations of Morgana Cave ( Figure 3). Most of the microbial features observed in samples M1, M3, M4, M5, and M6 resemble the enigmatic reticulated filaments reported in caves worldwide (Melim et al., 2008;Miller et al., 2012), which are long tubular filaments with peculiar surface ornamentation (Miller et al., 2012;Melim et al., 2015). In general, these long filaments (<100 µm) were found embedded in a slimy matrix of extracellular polymeric substances ( Figures 3A,D,F). These filaments depict mineralized empty sheaths ( Figure 3B), which can be attributed to the lysis of the cells and disappearance of living microorganisms.
The reticulated quadrangular-shaped filament depicted in Figure 3C is another type of intriguing structure, which, to our knowledge, has not been observed in other caves. In addition, flat filaments were found, which are probably mineral based on their morphology ( Figure 3E).
The analysis of all the samples with energy dispersive X-ray spectroscopy (EDS) revealed that the bulk samples M1, M2, M4, and M6 contain high amounts of Ca and Mg (Supplementary
Chao1, Shannon and Simpson indices were calculated to estimate alpha diversity ( Table 1). Regardless of variation in  Figure 4) revealed the presence of 137,298 OTUs that are unique to vermiculation deposits. Most of these unique taxa were found in sample M5 (33,646), followed by M3 (23,836) and M2 (27,320). The common microbial core comprised 6,829 distinct OTUs. The largest OTUs number (30,691) was shared between M4 and M5 with the highest bacterial diversity, whilst the lowest common prokaryotic OTUs were found between M2 and M3 (21,767) and M5 and M6 (21,628). Table 2 shows that the microbial communities of the samples from Morgana Cave were almost totally composed of Bacteria, ranging from 99.33 (in M2) up to 99.79% (in M5). Archaea exhibited a low percentage ranging from 0.65 (in M2) to 0.19% (in M5). Unclassified members at domain level ranged from 0.05 to 0.02%. This distribution of taxa is similar to those reported by Itcus et al. (2018) from a perennial ice cave, where the range was 94.90-99.03% for Bacteria and 0-1.03% for Archaea.

Archaeal Phyla
Four Archaea phyla were found in the vermiculations of Morgana Cave, Diapherotrites, Euryarchaeota, Nanoarchaeota and Thaumarchaeota ( Figure 5), although with very low relative abundances of OTUs. The class Thermoplasmata (Euryarchaeota phylum) ranged between 0.17 and 0.04%, Woesearchaeia (Nanoarchaeota phylum) between 0.37 and 0.08%, and Nitrososphaeria (Thaumarchaeota phylum) between 0.16 and 0.06%. The class Thermoplasmata was mainly present in M2, M3, and M6, and includes the order Methanomassiliicoccales with methanogenic archaea which accounted for 9.7% of all the Archaea present in the vermiculations. In addition, we observed an uncultured order of Thermoplasmata accounting for 17.3%.
The class Woesearchaeia was the most abundant (50%), considering all Archaea classes, particularly in M2 and M4. The class Nitrososphaeria was mainly present in M2, with the order Nitrosopumilales, accounting for 17.4% of all the Archaea. This last order includes Nitrosopumilus maritimus, an ammoniaoxidizing archaea (Pester et al., 2011).
The class Betaproteobacteria was composed of sequences of uncultured members of families from which only Nitrosomonadaceae and Burkholderiaceae were identified ( Figure 7; Table 4; Supplementary Table S1).
In Morgana Cave Planctomycetacia members were affiliated with the family Gemmataceae and other uncultured families and genera, with relatively low relative abundances (Supplementary Table S1).
Actinobacteria were mainly composed by the classes Acidimicrobiia and Rubrobacteria ( Figure 6). In the class Acidimicrobiia we observed sequences of the uncultured actinobacterium IMCC26256 (Figure 7), clearly separated from other Acidimicrobiia species (Hu et al., 2018). The class Rubrobacteria with the genus Gaiella is well represented in the different vermiculations (Supplementary Table S1).
Ignavibacteria is a recently established class of Bacteroidetes and was denoted by clone BVS26 in the vermiculations of Morgana Cave (mainly in M2 and M4).  Table S1).
Of the 26 existing subdivisions in the phylum Acidobacteria only five were found in Morgana Cave ( Figure 6): Acidobacteriia, Blastocatellia, Holophagae and subdivisions 17 and 6. Within the Blastocatellales the families Blastocatellaceae and Pyrinomonadaceae contained the uncultured genera JGI and RB41, respectively (Supplementary Table S1). In general, the relative abundances of the different subdivisions were low, and in most cases rarely exceeded 1%.
Regarding the distribution in classes in the different vermiculations ( Figure 6; Supplementary Table S1) vermiculation M2 presented the highest relative abundances in Beta-and Gammaproteobacteria. M3 showed a high diversity with the majority of relative abundances in 12 different classes, corresponding to eight phyla. In this vermiculation, Beta-and Gammaproteobacteria were also highly significant.
Vermiculation M4 revealed high relative abundances with respect to other vermiculations only in two classes of Chloroflexi and one of Bacteroidetes. M5 showed high relative abundances in 2 classes of Bacteroidetes, and one each of Acidobacteria, Planctomycetes, Zixibacteria and Verrucomicrobia. M6 presented high relative abundances in two classes of Chloroflexi, and one of Planctomycetes and Actinobacteria.

Functional Genes
To determine the involvement of the vermiculation communities in the nitrogen cycle, we surveyed the presence of functional genes involved in nitrogen fixation (nifH) and ammonia oxidation to nitrite (amoA). We observed nifH presence in the vermiculations M1, M2, M5, and M6. The phylogenetic analysis placed the sequences among nitrogen-fixing species (Supplementary Figure S2).
In addition, we detected the presence of functional genes involved in ammonia oxidation (amoA) in M1, M2, and M6 corresponding to Archaea. Regarding the bacteria, we identified genes related to Nitrosospira in M2 and M3. Phylogenetic   Figure S3A) indicated that the closest relatives of bacterial amoA gene sequences was the genus Nitrosospira, while the genus Nitrosopumilus was the closest to archaeal sequences (Supplementary Figure S3B).

DISCUSSION
Microbial communities typically show a few dominant phyla and simultaneously a high number of phyla with a quite low abundance, in most cases <0.1%. Sogin et al. (2006) referred to this as the "rare biosphere", microorganisms that cannot be identified without deep-sequencing NGS approaches. Only one third of the 35 less abundant phyla retrieved from the vermiculations of Morgana Cave have been also found in other caves (Ortiz et al., 2013;De Mandal et al., 2016;Hershey and Barton, 2018). For Pascoal et al. (2020) the rare biosphere constitutes an important "genomic reservoir" or "pool of diversity" and likely plays fundamental roles in ecosystem functioning. Frontiers in Earth Science | www.frontiersin.org December 2020 | Volume 8 | Article 586248 8 Vermiculations form an ecological niche suitable for microorganisms due to the presence of different mineral particles and water retention. Samples from vermiculations M3 and M5 were mainly composed of clay minerals (Figure 3; Supplementary Figure S1), and M5 showed greater diversity than dolomite-rich vermiculations. Organic compounds and bacteria are generally adsorbed onto clay minerals (Bingham and Cotrufo, 2015;Cuadros, 2017), which might explain the abundance of taxa in the clay-containing vermiculations with respect to dolomite-rich vermiculations.
In both types of vermiculations, EPS and reticulated filaments were evident. These enigmatic filaments have been found in microbial mats from caves worldwide, such as limestone caves and basaltic lava tubes (Melim et al., 2008;Jones, 2009Jones, , 2010Jones, , 2011Miller et al., 2014;Riquelme et al., 2015a), as well as in a granite spring water tunnel in Portugal (Miller et al., 2012). They were characterized as long microbial filaments with hexagonal and diamond-shaped chambers resembling honeycombed structures, but attempts to decipher their taxonomic affiliation have failed thus far.  Identities were based on nucleotide sequence similarities with their closest relatives in NCBI databases using BLAST hits.
Frontiers in Earth Science | www.frontiersin.org December 2020 | Volume 8 | Article 586248 Melim et al. (2008) reported long reticulated filaments (up to 75 μm long and 0.5 μm in diameter) in calcite speleothems from limestone caves in New Mexico, and in lava tubes from Cape Verde Islands. (Jones 2009(Jones , 2010(Jones , 2011 described the occurrence of reticulated filaments in caves from the British West Indies, and suggested they are calcified filaments. Contrasting with those filaments described by Melim et al. (2008) and (Jones 2009(Jones , 2010, rich in carbon and calcium, respectively, Miller et al. (2012) found reticulated filaments rich in Mn oxides. These filaments were composed of an inner cellular structure coated by an exterior sheath of manganese oxides (Miller et al., 2012). Melim et al. (2015) reported living reticulated filaments associated with clay minerals in a limestone cave in Germany. They were rich in carbon, without mineralized sheaths. Most of these filaments were coated with clays, obscuring wall ornamentation, as also observed in samples from Morgana Cave. In fact, this association with clays could explain their presence in the vermiculations M3 and M5. Shedding light on the origin and nature of these intricate reticulated filaments is a still a challenge for geomicrobiologists.
Archaea are ubiquitous in terrestrial and marine ecosystems (Eme et al., 2017), but the richness of this domain in the vermiculations of Morgana Cave was scarce. Tebo et al. (2015) did not find any archaeal rRNA gene sequences in the volcanic ice caves of Mt. Erebus, Antarctica. Itcus et al. (2018) reported the lack of Archaea in surface ice blocks from a perennial ice cave and negligible percentages in ice cores. In addition, Icelandic and Antarctic lakes were dominated by Bacteria with no detectable or few archaeal members (Achberger et al., 2017).
In Morgana Cave a few groups of microorganisms could be involved in the nitrogen cycle. Ward et al. (2019) recognized four clades of ammonia oxidizers, which included members of the Thaumarchaeota (Archaea), the Nitrococcaceae family (Gammaproteobacteria), the Nitrosomonadaceae family (Betaproteobacteria) and some members of Nitrospirae.
The order Nitrosopumilales (Thaumarchaeota) contains the families Nitrosopumilaceae and Candidatus Nitrosotenuaceae, which includes Nitrosopumilus maritimus and Candidatus Nitrosotenuis cloacae, respectively (Qin et al., 2017). Ammoniaoxidizing archaea were detected in the vermiculations through the identification of the gene encoding monooxygenae subunit A (amoA), however the relative abundance is very low as shown in Table 2. Among the bacteria, amoA genes related to Nitrosospira (Betaproteobacteria) were also found.
The high proportion of bacteria over archaea found in the vermiculations is consistent with previous reports indicating that in terrestrial ecosystems nitrate stimulates the activity of ammonia-oxidizing bacteria, but not that of archaea (Meinhardt et al., 2018). In addition, Zhang et al. (2019a) reported that bacterial amoA gene abundance was significantly higher than that of archaea due to the activity of cold adapted Nitrosospira in rivers from the Qinghai-Tibet Plateau. In this context, Zhang et al. (2019b) showed that soil microbial communities were adapted to high nitrogen levels, strengthening the function of ammonia oxidation in order to deplete the excess of nitrogen.
Among Bacteria, the phylum Proteobacteria and the classes Gammaproteobacteria and Betaproteobacteria were the majority in the vermiculations, with percentages in all the samples exceeding 14%.
The gammaproteobacterial order wb1-P19 was the most abundant group in the vermiculations. Holmes et al. (2001) proposed that wb1-P19 clustered phylogenetically with sulfur or nitrite-oxidizing autotrophic bacteria. Unfortunately, records on this lineage were very scarce, with only three references in the literature (Holmes et al., 2001;Schabereiter-Gurtner et al., 2002;Zhu et al., 2019). These three papers reported this order in caves. Remarkably, Zhu et al. (2019) found that in Chinese caves, in addition to wb1-P19, uncultured Rokubacteriales and Gaiellales were also present, similar to our results in the vermiculations. Table 4 shows the taxonomic affiliation of the representative sequences of OTUs classified as Gammaproteobacteria. In this case, the closest related sequences corresponded to uncultured bacteria from volcanic caves in the Terceira and La Palma islands, located in the Macaronesia region (Riquelme et al., 2015b;Gonzalez-Pimentel et al., 2018). This is striking due to the different environmental and climate (subtropical) conditions of the region and the mineral composition of the volcanic rocks (basalts) (Morgana is an Alpine cave developed on dolomitic marble).
Unfortunately, little information on metabolic traits has been reported for these uncultured volcanic cave bacteria. Gonzalez-Pimentel et al. (2018) stated that sulfur and nitrogen were associated with the chemical composition of basaltic rocks in lava tubes, and sulfur and nitrate oxidizers or reducers, such as bacteria from the Proteobacteria phylum, were very abundant in La Palma cave. Also, this phylum has similar abundances in Terceira cave (Riquelme et al., 2015b). Interestingly, Hathaway et al. (2014a) studied the diversity of ammonia oxidation (amoA) and nitrogen fixation (nifH) genes in lava caves from the Azores. They found that Nitrosospira-like sequences dominated the ammonia-oxidizing bacteria community, and the nitrogen fixation community was dominated by Klebsiella pneumoniae-like sequences. Note that Nitrosospira sequences were also present in vermiculations from Morgana Cave.
Apart from the close relationship between analyzed representative sequences and uncultured bacteria from the lava caves, secondary affiliations are indicated in Table 4. Some sequences from Morgana Cave were closely related to the gammaproteobacterial genera Methylococcus, Methylosarcina, Thioprofundum, Halochromatium, Lamprobacter, Thiocapsa, and Steroidobacter, with identity percentages ranging between 88.39 and 94.62%. According to Yarza et al. (2014), who established taxonomic boundaries for higher taxonomic levels for taxa of bacteria, the Morgana Cave uncultured bacterium close to Methylococcus and Methylosarcina must be included in the family Methylococcaceae. This family comprises methaneoxidizing bacteria that carry out the first step of nitrification, the oxidation of ammonia to nitrite, through the intermediate hydroxylamine (Hazeu et al., 1980;Wise et al., 2001;Sutka et al., 2003;Poret-Peterson et al., 2008).
The uncultured bacterium close to Thioprofundum must be placed within this genus. Thioprofundum hispidum is an Frontiers in Earth Science | www.frontiersin.org December 2020 | Volume 8 | Article 586248 obligately chemolithoautotrophic sulfur-oxidizing bacterium that bases its metabolism on sulfur oxidation and carbon dioxide fixation (Mori and Suzuki, 2014). The uncultured bacteria related to Halochromatium, Lamprobacter, and Thiocapsa should be included in the Chromatiaceae, a family of phototrophic purple sulfur bacteria growing under anoxic conditions using sulfide as photosynthetic electron donor, which is oxidized to sulfate via intermediate accumulation of globules of elemental sulfur (Imhoff, 2014).
The uncultured bacterium related to Steroidobacter, a genus classified under incertae sedis, is close to the family Solimonadaceae. Steroidobacter anaerobically metabolizes steroids with nitrate reduction to dinitrogen monoxide, and further to dinitrogen (Fahrbach et al., 2008) A Steroidobacter OTU was also found in Californian lava caves (Lavoie et al., 2017).
The class Betaproteobacteria was placed with lower relative abundance, represented by Nitrosomonadales and some uncultured members. The order Nitrosomonadales is composed of the genera IS-44 and MND1 of the family Nitrosomonadaceae. All the cultivated representatives of this family are lithoautotrophic ammonia oxidizers (Prosser et al., 2014). IS-44 is rarely cited in the literature and, among the scarce reports, was found to be a member of the microbial community of composts (Neher et al., 2013). MND1 was found in ferromanganese nodules, soils and caves and may be diverse in functions in a wide range of environmental conditions Spain et al., 2009). Table 4 shows the complete taxonomic affiliation of the representative sequences of OTUs classified as Betaproteobacteria. In this class, also OTUs from Morgana Cave were closely related (97-99%) with uncultured bacteria from the volcanic caves. In addition, OTUs affiliated to the order TRA3-20 were closely related to uncultured bacteria from a United States subsurface basaltic site (Lin et al., 2012) and from an Alpine cave in Austria (Reitschuler et al., 2015). However, other more distant affiliations, within the family Nitrosomonadaceae, were Nitrosospira lacus, a psychrotolerant ammonia-oxidizing bacterium (Urakawa et al., 2015). The finding of bacterial ammonia oxidizer communities at elevated altitude on Mt. Everest suggests that this group is sufficiently adapted to survive in Alpine areas (Zhang et al., 2009). Nitrosospira sp. and Nitrosovibrio tenuis are also ammoniaoxidizing bacteria (Prosser et al., 2014).
An uncultured bacterium within the family Burkholderiaceae is related to sequences of Piscinibacter aquaticus, Piscinibacterium candidicorallinum, Rhizobacter sp. and Methylibium sp., genera of uncertain position placed as unclassified Burkholderiales. Most strains from these genera are able to perform nitrate reduction (Sheu et al., 2016).
Within the genus MND1, family Nitrosomonadaceae, an uncultured bacterium closely related to another from Terceira lava cave, and more distant from Thiobacter subterraneus, a thermophilic, sulfur-oxidizing bacterium (Hirayama et al., 2005) and Cupriavidus sp., a genus including nitrate reducing strains (Vandamme and Coenye, 2004) was included.
Therefore, the scope of the metabolic diversity within the two main Proteobacteria classes, Gamma-and Betaproteobacteria, in the Morgana Cave vermiculations remains largely unknown, although there is great evidence of possible participations in the carbon, nitrogen and sulfur cycles.
In a lower level of abundance, two orders were found in the vermiculations: the actinobacterial Gaiellales, and Nitrospirales. Gaiellales comprises the genus Gaiella with only one known species: G. occulta, isolated from a deep aquifer with poorly mineralized waters. Distinctive characteristics are the presence of internally branched iso-fatty acids found in no other bacterium and the utilization of organic compounds (Albuquerque et al., 2011). The genome was studied by Severino et al. (2019) identifying genes of nitrate reduction. This genus was also found in volcanic caves (Riquelme et al., 2015a;Gonzalez-Pimentel et al., 2018).
The genus Nitrospira comprises ammonia-oxidizing bacteria which are relatively common in caves (Saiz-Jimenez, 2015). Its metabolic versatility, including the participation in different nitrogen cycling processes, is a key factor for the ubiquity and the high diversity of Nitrospira in caves (Tomczyk-Żak and Zielenkiewicz, 2016). Nitrospira was previously found in volcanic caves from the Azores, New Mexico, Hawai'i, and California (Northup et al., 2011;Riquelme et al., 2015b;Lavoie et al., 2017).
The phyla Rokubacteria, Gemmatimonadetes, Latescibacteria, Patescibacteria, Verrucomicrobia, and Zixibacteria had the lowest abundances in the vermiculations. They are commonly represented by uncultured bacteria that can only be identified in deep-sequencing NGS (Sogin et al., 2006). Most of these phyla are part of a rare cave biosphere formed by uncultured bacteria present in low numbers (Hershey and Barton, 2018).
Nitrogen fixation (nifH) genes have been found in Azorean and Brazilian caves (Hathaway et al., 2014a;Marques et al., 2018), suggesting this metabolic trait is widespread in the subsurface environment. Nearly all known diazotrophs contain a minimum of six conserved genes: nifH, nifD, nifK, nifE, nifN, and nifB, but for practical reasons in environmental studies only the nifH gene is commonly used as an indicator of nitrogen fixation activity. Although Dos Santos et al. (2012) stated that false positives can be obtained if only nifH is used, the data here reported and the phylogenetic tree of Supplementary Figure S2 place the sequences retrieved from the vermiculations among wellknown nitrogen-fixing species, thus suggesting that potential nitrogen-fixing bacteria are present in the cave.
In addition, Table 4 shows that some sequences from the vermiculations were related to the betaproteobacterial nitrogenfixing Azoarcus, and the gammaproteobacterial Thiocapsa. Furthermore, the sequence similarities of nifH genes with the genera Thiocystis, Chlorobium, and Aquaspirillum support the belief that uncultured bacteria have metabolic traits related to nitrogen-fixation in the vermiculations. Uncultured members of the Rhizobiales, Bulkholderiales, Methylococcales, Dehalococcoidia, and Euryarchaeota were also observed in the Frontiers in Earth Science | www.frontiersin.org December 2020 | Volume 8 | Article 586248 vermiculations and many species of these groups were reported as diazotrophs by Dos Santos et al. (2012). In Alpine caves it is important to investigate how microbes acquire energy and obtain organic carbon. In this context, one of the major impacts on oligotrophic caves is organic carbon, which enriches the ecosystem, brought into the cave from outside.
Atmospheric transport and deposition of carbon and nitrogen is an important dispersal process on the earth's surface. Several papers have reported the presence of organic carbon in the ice and snow of the Alps (Gröllert et al., 1997;Legrand et al., 2007;Preunkert and Legrand, 2013). This organic carbon derives from anthropogenic and biogenic sources, and is mainly composed of mono-and dicarboxylic acids, aliphatic alcohols, monoterpenes, phenols and humic-like substances Legrand et al., 2013;Müller-Tautges et al., 2016;Materic et al., 2019). These compounds constitute suitable carbon sources for microorganisms.
The presence of organic carbon, nitrogen and sulfur in the Alps suggests that its content in the meltwaters from the glacier Vedretto di Scerscen has likely enhanced the activity of cave microorganisms involved in the geochemical cycles, and particularly of the most copiotrophic groups, such as Betaproteobacteria and Gammaproteobacteria involved in the nitrogen cycle ( Figure 6; Table 4).
Chemolithoautotrophic bacteria provide a way to evaluate whether microbes are capable of producing their own energy and organic carbon to sustain an ecosystem. A possible metabolic pathway in vermiculations, in addition to nitrogen and sulfur, is methanotrophy. This is supported by the presence of the Rokubacteria class NC10, able of anaerobic oxidation of methane coupled to denitrification, as well as members of the Methanomassiliicoccales order which includes methanogenic archaea (Borrel et al., 2014). Uncultured members of the Methylococcales were present in the vermiculations suggesting that gammaproteobacterial methanotrophs may occur in this cave. Rokubacteria abundance ranged from 2 to 5%, Methylococcales from 0.3 to 2.3%, while Methanomassiliicoccales abundance was lower than 0.1% in all the vermiculations. These data indicate that bacteria metabolizing organic matter brought into the cave from external sources probably contribute more to the ecosystem than methanotrophic bacteria.
On the other hand, Alpine ecosystems are covered by snow and ice for most of the year, and ice and glacier caves are characterized by low temperatures below 0°C and inhabited by extremophiles. In these ecosystems, cave microorganisms have to cope with severe environmental constraints. Bowman (2017) studied the genomes of 257 bacterial and archaeal genera related to cold ecosystems, from which 59 genera could be considered highly cold adapted and 137 genera included psychrotolerant species. From the list only Flavobacterium was present in the vermiculations analyzed in this work. Members of the Flavobacteria group colonized cold springs in the Qinghai-Tibetan Plateau (Li et al., 2012) and represented 12% of the isolates from high-altitude ice cores in the Tibetan Plateau (Liu et al., 2019).
Flavobacterium, Cytophaga species and Ignavibacteria contain antifreeze DUF3494-type proteins, which are present only in cold-adapted taxa (Bowman, 2017). The finding in the vermiculations of uncultured members of the classes Cytophagia and Ignavibacteria, in addition to the genus Terrimonas, from which novel species were isolated from Arctic environments (Kim et al., 2017), suggests that psychrophilic bacteria may be present in the vermiculations. Psychrophilic ammonia-oxidizing bacteria were also reported in cold environments (Zhang et al., 2009;Urakawa et al., 2015).
In spite of the extensive investigations on the microbiology of several glacial environments, very little is known on the diversity, activity and environmental role of microorganisms colonizing Alpine caves. In previous studies of ice caves, the presence of unclassified bacteria was common (Itcus et al., 2018;Purcarea, 2018). This also occurs in the vermiculations of Morgana Cave, where most of the retrieved phyla contain a large number of uncultured representatives, with unknown metabolic traits. As a matter of fact, many taxa common in subsurface soils are represented by high proportions of novel and undescribed microbial lineages, less amenable to currently known cultivation techniques (Brewer et al., 2019).
It is noteworthy that some of the Beta-and Gammaproteobacteria OTUs from the Alpine Morgana Cave (Table 4) were closely related to sequences from two subtropical volcanic caves in the Macaronesia region (Azores and Canary islands). Hathaway et al. (2014b) suggested that the distribution of lava cave organisms is not cosmopolitan, as they found substantial differences in the microbial diversity of lava caves from the Azores and Hawai'i. In fact, at the OTU level, there was only 5% overlap between the two archipelagoes. Nevertheless, the similarities of OTUs between our high Alpine and two subtropical volcanic caves are remarkable. Porca et al. (2012) studied microbial communities from geographically distinct limestone caves and reported the existence of a common core of microorganisms that were morphologically similar and related in phylogeny. These authors suggested that, in nutrient-limited environments, in order to efficiently exploit nutrients, cave microorganisms could develop complex cooperative and mutualistic associations which could explain the broad distribution of bacterial phyla in caves.
Conversely, it was clear that the comparison of the microbial communities from Morgana Cave vermiculations with those of the Frasassi and Fetida sulfuric acid caves revealed striking differences, due to the extremely different environments (alkaline in Morgana and highly acid in sulfuric caves). In fact, Fetida Cave vermiculations showed high abundance of unclassified Betaproteobacteria and sulfur-oxidizing Hydrogenophilales (including Sulfuriferula), Acidiferrobacterales (including Sulfurifustis), sulfur-reducing Desulfurellales, and ammoniaoxidizing Planctomycetes and Nitrospirae (D' Angeli et al., 2019b).
The Frasassi Cave vermiculation community was diverse, including as major lineages Betaproteobacteria (composed of Frontiers in Earth Science | www.frontiersin.org December 2020 | Volume 8 | Article 586248 unclassified bacteria, some related to sulfur-oxidizing bacteria), Gammaproteobacteria (composed of relatives of Piscirickettsiaceae, Nevskiaceae, Acidothiobacillaceae, and unclassified bacteria), Acidobacteria, Nitrospirae, and Planctomycetes . The data from Frasassi and Fetida sulfuric acid caves suggests that pH and sulfur are important drivers in the composition of their vermiculation microbial communities. On the contrary, vermiculations from the carbonate Morgana Cave have different communities, indicating that the microbial composition of vermiculations largely depends on environmental conditions and other factors. To summarize, the microbial groups identified in the vermiculations of this Alpine cave differed significantly from those retrieved from sulfuric acid caves, suggesting that mineral composition and pH play a major role in the colonization and establishment of communities.
This study provides new insights into vermiculation microbial communities from Morgana Cave which are composed of cosmopolitan bacteria in addition to a few niche-specific members. It is remarkable that this Alpine cave, supposed to be oligotrophic in nature due to the altitude and its remoteness, revealed that the 2 major bacterial classes Gammaproteobacteria and Betaproteobacteria, are composed of copiotrophic members whose metabolic traits are mainly associated with the nitrogen cycle. Evidences of the presence of methane-and sulfur-oxidizing bacteria are also found. This highlights the importance of atmospheric transport and deposition of anthropogenic pollutants which enhance the colonization and activity of microorganisms in pristine environments.

CONCLUSION
The microbial communities from Morgana Cave vermiculations were composed of different major evolutionary lineages. The vermiculations were dominated by the phyla Proteobacteria, Actinobacteria, Acidobacteria, Chloroflexi, Planctomycetes, and Bacteroidetes. The most abundant classes were Gamma-and Betaproteobacteria, followed by Planctomycetacia, Acidimicrobiia, Rubrobacteria, Ignavibacteria, NC10 from the phylum Rokubacteria, and Blastocatellia. The microorganisms linked to the nitrogen cycle are noteworthy, and their abundance must be related to atmospheric deposition of different forms of nitrogen coming from glacier meltwaters and rain/snow melt infiltrating into the cave. In addition, psychrophilic and methanotrophic bacterial groups were identified. The high number of phyla and candidate phyla (48) found in this oligotrophic high mountain cave is remarkable. The high percentage of uncultured bacteria present in the vermiculations, particularly at the lowest taxonomic ranks (only four genera were identified) indicates that the microbial diversity of Alpine caves deserves further research.

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 in the article/Supplementary Material.

ETHICS STATEMENT
Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

AUTHOR CONTRIBUTIONS
VJ performed the molecular biology, VJ, BH, and JG-P contributed to bioinformatics and statistical data analyses and AM to FESEM. ID'A and PT participated to cave sampling and described the sampling locations and the caves. JDW and CS-J led the research. JDW participated in discussions and corrected the manuscript. CS-J wrote the manuscript. All of the authors reviewed the manuscript and approved the final version.