Metagenomic Insights Into the Diversity of Halophilic Microorganisms Indigenous to the Karak Salt Mine, Pakistan

Hypersaline regions are terrestrial analogs of the Earth’s primitive ecosystem and extraterrestrial environment. The salt range in Pakistan is considered among a few of the ancient salt deposits in the subcontinent. Karak salt mine is situated at the Northwest end in Pakistan. Despite the fact that halophiles initiated the formation of terrestrial ecosystems, their products and identities remain hidden. Some preliminary studies limited to culture-dependent isolations have been reported. Characterizing the microbiome that spans over centuries of ecosystem development is crucial, given their role in shaping landscape succession and biogeochemical cycles. Here, we used metagenomics techniques to explore the microbial diversity of the Karak salt mine. We used 16S rRNA Illumina amplicon sequencing to characterize the halophilic communities entrapped in Karak mine. The results were interpreted using Illumina Basespace, QIIME, and Cytoscape. Cultures were isolated at 16–25% salinity. Metagenomics data was consistent with our preliminary culturing data, indicating remarkable species to strain-level diversity of unique halophiles. A total of 107,099 (brine) and 122,679 (salt) reads were obtained. 16S rRNA based sequencing revealed a microbiome with bacteria (66% brine and 72% salt) dominated by Bacteroidetes and Proteobacteria with a strikingly high abundance of Archaea (18% brine and 13% salt). Alpha diversity has higher values in salt than in the brine. The study of the halophiles in the Karak salt mine provides clues for species contributing to the maintenance of biogeochemical cycles of the ecosystem. This is the first report of a metagenomic study of any hypersaline region of Pakistan.


INTRODUCTION
It has only recently been accepted that microorganisms commonly inhabit evaporite formations. Microbes in halite crystals serve as excellent models to study evolutionary patterns of ancient microbes because signs such as initial crystallization, surface sterilization, in situ observations of the halites and efficient reproducibility of experiments can authenticate the results (DasSarma and DasSarma, 2017). Determining the species colonizing salt crystals, aids in understanding evolutionary patterns and factors influencing microbial survival in paleoenvironments (Jaakkola et al., 2016). Pakistan possesses all three kinds of salt (rock, sea, lake) available on Earth. The hypersaline environments are distributed at the entire length of Pakistan in the form of rock salt, soda lakes, and salt lakes. The salt range in Pakistan represent the most ancient salt deposits on the Asian subcontinent and reveals a series of rock formations ranging from the pre-Cambrian era to the most current geological periods (Cremo, 2001;Jehangiri et al., 2015;Leena et al., 2018). Such exposures are seldom observed and therefore, are remarkable geological phenomena. In the Province of Khyber Pakhtunkhwa (KPK), reserves of salt deposits are found in the mountain range in Karak district at the Northwest end of the Khewra salt range. The salt deposits, also called Kohat rock salt, outcrops the length of 12 km and width of 1/2 km. These mines have unique stratigraphic constraints and fossils, dating its origin millions of years ago. A Pakistan based geological survey assigned the age of Karak salt to the early Eocene period (Roohi et al., 2014;Khattak et al., 2016). This suggests that the microbes entrapped in their salt crystals may have preserved ancient DNA sequences representing the Eocene life. Despite the fact that halophiles may have initiated the formation of terrestrial ecosystems, their products and identities remain enigmatic. A few preliminary studies limited to culture-dependent bacterial isolations and halotolerant enzymes producing bacteria from Karak salt mines have been conducted so far (Roohi et al., 2014;Shah et al., 2017Shah et al., , 2018. With access to next-generation sequencing, metagenomic approaches have been effectively used to analyze microbial diversity in extreme conditions. They have proven very useful in the reconstruction of the draft and even closed genome sequences of uncultured microorganisms and have recently been used to explore microbial population and composition in diverse hypersaline environments (Mesbah et al., 2007;Ghai et al., 2011;Fernández et al., 2016;Vavourakis et al., 2016). Metagenomic analyses shed light on the interplay between microbial communities and their environment representing ecological diversification. Furthermore, metagenomic studies also allow us to understand the interactions between environmental and evolutionary processes and the role of microbes in the process of nutrient cycling. In this report, we provide the first taxonomic snapshot of microbial communities from the Karak salt mines using metagenomic analysis (Illumina 16S rRNA amplicon sequencing).

Site Description
The Karak salt mines are located at 32 • 39 0 N, 73 • 1 0 E at an altitude of 678 m/2227 ft. The region consists of rounded dry hills. The rock salts are colored in the shades of light gray to dark gray (Figure 1). Rainfall is meager and intermittent in the area. Also, the environment is dynamic. The region is very hot in summer and very cold in winter (Khan et al., 2011). The geological horizon predicts the age of these rocks as early Eocene (Sharif et al., 2007). The raw salt from mine was recorded to have 0.12% moisture content and 7.98 (mg/Kg) water-insoluble impurities.
The composition of salt is fine to medium crystalline, friable and loose, with a slightly sugary texture. The feature distinguishing its halite from the rest of the salt range crystals is its porphyritic nature (Sharif et al., 2007;Nafees et al., 2013;Roohi et al., 2014).

Sampling and Physiochemical Analysis
Salt and brine samples were collected from the mine under conditions preventing the ingress of exogenous microbial contamination. Brine dripping from the efflorescence of the mine wall and from saturated brine pools was collected in sterilized sampling bottles. The salt crystals were scraped from the mine walls using a sterilized spatula and stored in sterilized bottles and stored in 4 • C until further analysis. All samples were sealed and transported to the laboratory (Quaid-i-Azam University, Islamabad) within a few hours there to at 4 • C. GPS location coordinates were recorded from each sampling location. The pH and temperature of the brine were determined at the sampling site using pH paper, and later measured under laboratory conditions by a pH meter (Sartorius professional meter PP15). The salinity was measured using a hand-held refractometer and total dissolved solids were also calculated. The chemical analysis of major ions (Sulfates, Chlorides, Phosphates, Bicarbonates) in the brine was performed with Spectroquant (Pharo 100 Spectroquant R Merck, Germany) using standard protocols for the examination of water and wastewater by (APHA) American Public Health Association.

Reagents and Chemicals
Reagents and chemicals used in this study were of analytical grade and obtained from BDH (Dorset, United Kingdom), Sigma-Aldrich (St. Louis, MA, United States), E Merck (Darmstadt, Germany), Difco Laboratories (Detroit, United States), AVACADO Research Chemicals Ltd (London, United Kingdom).

Isolation of Halophilic Microorganisms
For our preliminary analysis of detecting culturable halophiles in samples, we used the standard dilution-plating method to count and isolate cultivable microorganisms in salt and brine samples. A series of dilutions (using sterile 1.5 M NaCl solution) was prepared for each sample from the stock suspension, 1 g salt, and l mL brine in 99 mL of 1.5 M NaCl solution. Aliquots, 0.1 mL, from various dilutions (10 −1 to 10 −4 ) were spread on the surfaces of plates containing the solid Nutrient Broth medium (NA) composed of (in grams per liter) (pH 7.0, peptone 5, yeast extract 3, and NaCl (4, 16, and 25%). Solid media were amended with 2% agar. We used the NA medium for the isolation of halophilic bacteria. For haloarchaea, the mineral constituents of the specific medium CM+ have the following composition in (g/L); NaCl: 250, Magnesium Sulfate: 20, (Tri) Sodium Citrate: 3, Bacteriological peptone: 10, Trace metals: 0.1. To solidify the media, 2% agar was added, and the medium pH was adjusted to 7.2. The plates were sealed and incubated at 37 • C for 10 days to 2 weeks. Three replicate plates were prepared for each dilution. The total colony numbers were counted, and colony-forming units (CFU) per gram of sample were calculated. Three replicate plates of dilutions containing countable numbers of colonies (100-200) were pooled. Identical colonies (colony shapes, sizes, colors, margins, texture, etc. and cell shapes and sizes as well as the cell Gram stain reaction) were counted, their proportions of the total CFU were calculated, and five replicate representatives were sub-cultured, purified and maintained for further study. For the molecular characterization of the isolates, the genomic DNA was extracted using the ethanol precipitation method (Delbès et al., 2000) with slight modification (Leena et al., 2018). The genomic DNA is confirmed by gel electrophoresis by running it on 1% gel in 0.5 × TBE buffer. Gel images were observed by Bio-Rad gel documentation system. DNA. The 16S rRNA genes in the total genomic DNAs were amplified by polymerase chain reaction (PCR) using the universal primer combination 8F: 5 AGAGTTTGATCCTGGCTCAG and 1492R: 5 AAGTCGTAACAAGGTAACC for bacteria and 4F: 5 TCCG GTTGATCCTGCC and 1492R: 5 AAGTCGTAACAAGGTAA CC for Archaea. Approximately, 5-10 ng of NanoDrop quantified DNA was used for amplifying the entire 16S rRNA with universal primers. The PCR reaction was carried out with PCR mixtures (25 µL) containing 25 pmol of each primer, 1 µL of 100% DMSO, 0.25 µL of 10 mM dNTPs, 0.25 µL of Taq polymerase, 2.5 µL of 10x buffer and 2 µL of the DNA template. The PCR program consisted of an initial denaturation step at 94 • C for 5 min; 35 cycles of denaturation at 94 • C for 30 s, primer annealing at 55 • C for 45 s, and extension at 72 • C for 2 min, followed by a final 5 min elongation step at 72 • C. Partial sequencing of the 16S rRNA genes was performed using the 3130 XL Genetic Analyzer (BASLab Baltimore, MD, United States). Sequences were subjected to a basic local alignment search tool analysis (BLAST) with the National Center for Biotechnology Information (NCBI) GenBank database. Sequenced 16S rRNA genes data was assembled manually after checking for quality assurance utilizing BioEdit software version 7.3.5. (Hall et al., 2011).

DNA Extraction and Sample Preparation
Environmental DNA was extracted from brine and salt samples. The brine samples were filtered through 0.45 µm mixed cellulose ester (MCE) membranes. The membrane is used for DNA extractions. DNA extraction protocol from salt and filtered brine were carried out by bead beating using the PowerLyzer R PowerSoil R DNA Isolation kit (DNeasy PowerLyzer PowerSoil, QIAGEN Inc., Germantown, MD) with some minor modifications. The alternate lysis protocol that is suggested for hard to lyse cells was used. DNA was quantified using NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific Inc.). Quality was additionally checked by electrophoresis on an Agarose gel (1%), and DNA was stored at -80 o C. Throughout the DNA extraction, preparation, and sequencing processes, a positive control, and negative control, without DNA were used to monitor for contamination. Approximately, 5-10 ng of NanoDrop quantified DNA was amplified using standard Universal primers provided by Illumina targeting the V3-V4 region of the 16S rRNA gene. The PCR products were then subsequently used for Miseq Illumina sequencing according to the manufacturer's instructions.

Illumina Sequencing
To profile the prokaryotic community, 16S rRNA amplicon sequencing libraries were prepared according to the manufacturer's instructions (Illumina, San Diego, CA, United States). Briefly, the hypervariable V3-V4 region of the 16S rRNA gene was amplified using the primer pair designed against the surrounding conserved regions (Klindworth et al., 2013). Purified PCR products were pooled in equimolar amounts, and sequenced using the Illumina MiSeq platform (250 bp paired-end reads). The DNA library was sequenced on an Illumina MiSeq platform at the Institute of Marine and Environmental Technology (IMET) using a MiSeq Reagent Kit 300 v2 (Illumina, Baltimore, MD, United States).

Sequence Analysis
Raw reads were preprocessed using CLC Genomic Workbench (version 10.1.1) (Qiagen). Adapter sequences were removed and read pairs were quality trimmed (qual. limit = 0.05; ambiguous nucleotide max. = 2; min sequence length = 100 bp) and merged (mismatch cost = 2; gap cost = 3; max. unaligned = 0; min score = 8). Bioinformatics analyses were performed using Quantitative Insights into Microbial Ecology (QIIME), version 1.33, and Illumina Basespace. Open reference OTU picking was performed against the GreenGenes13_8 database (similarity threshold = 0.97; OTU picking method = Uclust; minimum cluster size = 2, however, to generate the OTU network minimum count was kept 10, while OTU picking. OTU networks were generated using the make_otu_network.py script. For this analysis, OTUs were re-picked according to the above method, however, the minimum cluster size was increased to 10 in order to reduce the number of nodes. The network was visualized using Cytoscape v3.2.1 using an edge-weighted spring-embedded layout. QIIME pipeline was used to align pair-end reads, assign taxonomy, generating taxonomic plots, and phylogenetic tree. Taxonomies were summarized at multiple levels (L2-L6) using summarize_taxa.py script. Representative sequences were aligned with the PyNast algorithm against the Greengenes core set 1 . For comparative analysis, the data is rarefied to an even depth using single_rarefaction.py script in QIIME. Using alpha_diversity.py script, several diversity and richness metrics were calculated for each community such as observed OTUs including Shannon index for diversity and Chao1 and Simpson richness estimators. The rarefaction plots (Figure 3) were computed at 97% similarity, as part of QIIME's alpha diversity pipeline. OTU networks were constructed from the output BIOM tables and visualized in Cytoscape (version 3.3.0). Multivariate biome charts for relative abundance were constructed using Krona tools (version: 2.6.1).

Analysis of Samples Isolated From Karak Brine
The pH of the brine was 7.14 and the temperature was 28 • C at the time of sampling. The salinity of brine measured by a refractometer at 32%, and the Total dissolved solids (TDS) were 35%. Chloride values from brine samples from the Karak mine were estimated as 2002.2 mg/l and sulfates and bicarbonates were measured as 770 and 333 mg/l, respectively.

Preliminary Data of Isolation and Characterization of Extremely Halophilic Archaea and Bacteria
A total of six halophilic strains were isolated from brine and soil samples, out of which five extreme halophilic bacteria KPS3A, KPS8A, SWM2A, KK1A, and KK2A were selected that showed growth above 16% NaCl concentration and one Archaea was isolated at 25% NaCl, designated KKW1. These selected isolates were then characterized morphologically, and their cellular morphology was observed using phase contrast microscopy. BLASTn results determined the isolates had 16S rRNA gene sequences most similar to Oceanobacillus sp., Oceanobacillus oncorhynchi subsp. incaldanensis, Brevibacterium linens, Salicola marasensis, Chromohalobacter salexigens, and Haloferax alexandrinus strains. The sequences were submitted in the NCBI GeneBank with the following accession numbers KPS3A-MT406254, KPS8A-MT406253, SWM2A-MT406257, KK1A-MT406258, KK2A-MT406255, and KKW1-MT406256, respectively, where each designated specie is followed by its respective accession number (Table 1).

Culture-Independent Studies
16S rRNA reads generated from both brine and salt datasets were used to build a representation of the high-level phylogenetic groups in each dataset using QIIME. Both datasets display a remarkably diverse collection of taxonomic clusters. The outcomes of all methods of classification were largely in agreement and showed exciting merits.

Statistics
The brine (KKB) and salt (KKS) samples from the Karak salt mines were subjected to high-throughput sequencing analyses. Sequences acquired after the quality filtering were trimmed, and only sequences of high quality were preserved for further analysis. A total of 229,778 high-quality prokaryotic reads were obtained. After the rarefaction of data to an even depth, the analysis presented a total of 5491 operational taxonomic units (OTUs) using ≥97% sequence identity. Interestingly, salt sample KKS (3911 OTUs), corresponding to brine sample KKB (2234) contained a higher number of distinct OTUs, and a total of 635 shared OTUs were detected between both samples. KKSsalt contains the highest number of unclassified OTUs (27%). We recovered 46 classes, 77 orders, 106 families, and 156 genera across both sampling sites, from reads categorized at each taxonomic level. For both sample datasets, over 80% of all quality reads were constantly categorized at the Phylum level. However, about 16.3% of reads remained unclassified (Figure 2). Alpha diversity was measured using metrics such a Good's coverage, Shannon index, Simpson, and nonparametric Chao1 estimator ( Table 2). The results indicated that between the two samples, KKS-salt had the highest diversity as compared to KKS-brine as demonstrated by Shannon indices of 7.11 and 5.55, respectively. Due to the elimination of singletons, the number of observed OTUs was observed to be very close to the measured richness values using the indices mentioned above. The results of Good's coverage (99% coverage), and rarefaction curves (Figure 3) suggest an almost complete sampling. A comparison of the rarefaction curves with the number of OTUs determined by the Chao1 richness estimator and other diversity indices revealed that the samples have significant microbial diversity.

Community Structure
Comparative analysis of both the brine (KKB) and salt (KKS) datasets at the Phylum level showed that 27% reads from the KKS-salt sample and 3% of the reads are from KKB-brine sample remain could not be assigned to any Phylum. For both the datasets, the microbial community comprised of 60% of bacteria and 22% of archaea. A total of 28 bacterial phyla were identified in both the Karak salt and brine samples. The top ten abundant phyla contributed up to  (Figure 4). Euryarchaeota contributed 100% to the total archaeal diversity. At the genus level, 73 and 95% of the 16S rRNAs reads were assigned to KKB-brine and KKS-salt datasets, respectively. The archaea were dominated by numerous clades of Halobacteria that dwell in salt lakes and salterns and can propagate on or within salt crystals. Many archaeal OTUs between the two samples are shared. However, a significant number of unique archaea were detected in the sample with high salinity (KKS-salt) (Figure 4) Figure S1). Phylum Proteobacteria were ubiquitous across both samples. FIGURE 3 | Rarefaction curves computed through QIIME depicting distribution and diversity of microbial communities in salt (blue) and brine (red) samples using Chao1 richness estimator and observed OTUs at 97% sequence similarity.  (Figures 5, 6).

DISCUSSION
To fully comprehend the role that microorganisms play in Earth's Subsurface environment, it is essential to study the life entrapped in the closed environment of the salt mine. Microbes in salt crystals have implications for the preservation of life in the fluid trapped in the salt crystals at hypersalinity (Adamski et al., 2006;Pasteris et al., 2006). Salt crystals from multiple environments such as salt mines; salt plains, rock salts, and food-grade salts have been shown to harbor viable, ancient microbes (Vreeland et al., 1998;Caton et al., 2009;Henriet et al., 2014;Robinson et al., 2015). However, to the best of our knowledge, metagenomic community profiles of halophiles from these salt mines have not been reported. Hence, in this study, we present for the first time the novel insights into the community structure of indigenous microbiota of Karak salt mines by using metagenomics, which provides a snapshot of microbial community diversity, distribution, and functionality. The salt crystals of Karak salt mine are reported as 98% pure (halite). The salinity of the brine pools running within the mine runs up to 32%, hence reaching close to the higher range of salinity of brines from Andean salares which had also been investigated for prokaryotic diversity (Maturrano et al., 2006;Haferburg et al., 2017). In this study, isolations were performed at near-neutral pH and temperature ranging from 37 • C in accordance with the conditions measured at the investigated site. Our results on cultivable halophiles from Karak salt mine coincide with a few studies conducted over the years to isolate halophiles from the rock salts starting from isolations from the Permian salt mine, Kansas (Reiser and Tasch, 1960). Recent isolations of halophiles were reported from Great Salt Plains of Oklahoma (Caton et al., 2009) Death Valley California (Mormile et al., 2003;Ram et al., 2016) China, Bolivia (DasSarma et al., 2019 and Anatolia, Turkey (Yildiz et al., 2012;Xiao et al., 2013). In Pakistan, only a few studies regarding bacterial isolations from saline environments have been conducted, reporting the dominance of genus Bacillus in the rock salt along with Planococcus, Halomonas, Brevibacterium, Kocuria, Salinivibrio, Salinicoccus, and Kushneria (Muyzer et al., 1993;Akhtar et al., 2008;Roohi et al., 2014). Our preliminary isolation study was focused on the isolation of extreme halophiles (both archaea and bacteria) so that in addition to earlier reports of isolation, which comprise of mostly, slight to moderate halophiles, our report contributes in capturing the extreme halophilic culturable microbial community of the rock salt deposits of Karak mine. Isolates Chromohalobacter (KK2) Haloferax (KKW1) and Oceanobacillus sp. (KPS3A and KPS8A) were also detected in the metagenomic data. In this investigation, isolation yielded several species such as Haloferax and Salicola that probably belong to the unique biodiversity that is not adequately present to be identified by metagenomics but may be opportunists and effortlessly grow under the applied cultivation conditions. Many genera that are detected by metagenomics, may not be isolated soon due to specie's yet unknown growth requirements (Bolhuis et al., 2004;Naghoni et al., 2017). To attain a better comparative evaluation between cultivation approaches and metagenomics, a multitude of isolates would be required to grow under variable conditions. The metagenomic technique serves as a powerful tool to study microbial diversity directly from the environment and also allows the exploration of novel microorganisms and their adaptation mechanisms in relatively unexplored extreme habitats (Behzad et al., 2016). Miseq-Illumina analysis showed the dominance of Proteobacteria and Bacteroidetes with a strikingly high abundance of Archaea. Members of these taxa were formerly verified in other hypersaline and tropical soda lake ecosystems such as Salar de Uyuni, Lake Nakuru, Lake Bogoria, Kae Magadi as well as from the Lonar Lake ecosystem (Duckworth et al., 1996;Margesin and Schinner, 2001;Vargas et al., 2004;Kambura et al., 2016;Paul et al., 2016). The Karak salt mine has a metabolically diverse halotolerant bacterial population. Members of Proteobacteria phylum were found to be the principal inhabitants of the Karak mine, covering almost half of the bacterial OTUs (48.3%) within the data. Within Phylum Proteobacteria, members of the Altermonadaceae, Rhodospirillaceae, Ectothiorhodospira, Oceanospirillaceae, and Thiohalorhabdaceae were found to be prevalent in our taxonomic review and had also been independently verified in such ecosystems using culture-based methods (Chandrasekhar et al., 2002;Paul et al., 2016). These bacteria are known for their metabolic diversity and ubiquitous nature, which allows them to survive in extreme environments by utilizing various carbon compounds to preserve both aerobic and anaerobic lifestyles (Janssen, 2006). OTUs belonging to the Phylum Bacteroidetes amount up to ∼18% of the total reads. The decreased number of reads associated with anoxygenic photosynthetic (sulfur and non-sulfur purple, and green bacteria) in Karak samples at high salinity was consistent with other reports. However, further study is needed to determine the impact of these groups of Archaea to the primary productivity in the salt mines samples. Investigation of the average community composition of the brine and salt samples showed that the most greatly represented genus was Marinobacter which contributed up to 23.6% of the total data followed by Salinibacter (10%), Halorubrum (3%), Haloarcula (2.5%), and Halorhabdus (2.1%). Within Archaea, clear predominance of the Class, Halobacteria was observed in Karak mine samples. A massive number of the archaeal OTUs FIGURE 4 | OTU Network demonstrating the number of unique and shared OTUs in the brine and salt samples of Karak mine. The OTU table was generated using 97% cut-off sequence value and mc-10. The taxonomic clustering was displayed in the form of OTU network using the edge-weighted, spring embedded layout in Cytoscape. Node color denotes consensus taxonomy (phylum). Edges connecting OTU and sample indicate the presence of the OTU in that sample.
All of the archaeal sequences fell into the Class Halobacteria. In the present study, Halorubrum was found to be the most abundant genus in brine, and Haloarcula was the most abundant genus in salt samples suggesting their significant roles in their respective habitat. Abundant strains of the genera Halorubrum and Haloarcula in the halite crystals have also been observed in other culture-based studies of salt mines in China and foodgrade salt crystals (Gibtan et al., 2017). Halorhabdus is the most extreme halophilic member of the family Haloarculaceae and is present in much more abundance in salt (3.3%) than in brine (0.8%). Other archaeal genera such as Halobaculum (0.6%) and Halotinus (0.6%) were also found in relatively more abundance in the salt sample. Members of genus Halonotius are reported to be ubiquitously present in highly saline environments (Durán-Viseras et al., 2019). However, Haloquadratum, which is one of the most highly represented archaeal genera in salterns (Welsh, 2000;Oren, 2012;Dillon et al., 2013) appeared to be absent in the Karak mine samples. Haloquadratum was also reported as an insignificant member in commercial salts suggesting its role limited to hypersaline pond systems and salterns. It has previously been suggested that widely used 16S rRNA-targeted primer sets may not be suitable to indicate the presence of Nanohaloarchea and so using archaeal specific primers could probably provide a more in-depth diversity analysis (Narasingarao et al., 2012). The diversity of microbes within the sediment and water samples was distinct concerning community composition. The OTU network generated from the16S rRNA-based metagenomic analysis of the assembled reads verified the overlap in the diversity between the brine and salt samples. The overlap in sequences mostly consists of Archaea, Bacteriodetes, and Proteobacteria. Candidate Div. OP1, Firmicutes, and Planctomycetes were exclusively noted in the salt samples, while Gemmatimonadetes were found chiefly in the brine samples. The diverse patterns of microbial community assemblage across salt and water are likely to be caused by variances in salinity, nutrient availability, and redox state in this ecosystem.
The most highly represented bacterial OTUs belonged to genera Marinobacter. It was found to be the most abundant genus (35%) in both salt crystals and brine samples. Its presence in rock salts has also been described in the salt mines of China, Romania, and Turkey, as well as in the food-grade salts (Carpa et al., 2014;Henriet et al., 2014;Cınar and Mutlu, 2016;Haferburg et al., 2017). Members belonging to Marinobacter family were found in diverse environments such as saline soils, petroleum field brines, surface seawater, ballast water (Arctic), marine snow, hydrothermal plumes, hot springs at coasts, deep seawater, volcanic basalts and beneath the iron ranges of Minnesota (Singer et al., 2011;Handley and Lloyd, 2013;Bonis and Gralnick, 2015). Subsequently, the strains show diverse potential to survive poly extreme conditions such as heavy metal, antibiotic  resistance, phage resistance, Fe (II) oxidation, and electrogenic Perchlorate-Reduction allowing for active exploitation of an unstable environment (Stepanov et al., 2016). Marinobacter strains have substantial metabolic versatility, which makes them attractive for biotechnological exploitation. The dominance of Marinobacter in the Karak salt mines sheds light on its significant contribution to the maintenance of the ecosystem. Hypersaline environments are considered as the terrestrial analogs of the extra-terrestrial environment. Given the discovery of perchlorate among the salts detected by the Phoenix Lander on Mars, an abundance of Marinobacter in Karak salt mines, its perchlorate reducing ability and its ability to survive under poly-extreme conditions implies that presence of perchlorate does not preclude the prospect of halophilic life on Mars. This argument is also supported by the abundance of other perchlorate-resistant halophiles in the salt samples such as [Haloarcula, Haloferax, Halobacterium, Halorubrum (Laye and DasSarma, 2018) and the bacterium Halomonas elongata] (Oren, 2014). Hypersaline habitats have been known as a sanctuary to preserve biological material on time periods exceeding that of non-saline systems. Thus, the presence of paleo-sulfate lakes on Mars is of particular significance for astrobiology, as these might have preserved signs of previous microbial inhabitation on Mars (Carrier et al., 2020). Further biological characterization of the Karak mines is essential to define the limits to microbial growth in these hypersaline environments. With a relative abundance of 10%, the genus Salinibacter represents the second most abundant extremely halophilic member thriving in the brine sample. The genus is distributed globally in hypersaline environments (Oren, 2013).

Taxa With Known Ecological Roles
In addition to salt tolerance, a number of geochemically relevant species were detected in our data. A low abundance of anoxygenic phototrophs (Chloroflexi, Chlorobi, Rhodobacterales, and Chromatiales) was detected. The low abundance of phototrophs implies that in the Karak mine, particularly in salt samples, carbon is fixed using alternative pathways. Perhaps little known and relatively uncharacterized phylogenetic groups, exclusive to salt crystals such as the Planctomycetes, Firmicutes, Acetothermia, and Euryarchaeota, are fulfilling this significant role. Acetothermia belonging to Candidate div OP1 are currently uncultivated but has been identified in a variety of environments using culture-independent methods. One report revealed genes encoding for the (acetyl-CoA) pathway for CO 2 fixation in Candidatus "Acetothermum autotrophicum." The presence of genes for CO 2 fixation in Acetothermia suggests that they could contribute to primary production through chemolithoautotrophic acetogenesis (Farias et al., 2017). The sulfur cycle seemed more active in the salt samples as compared to brine with the presence of Thiohaldorhabdus, Halorhabdus, and Halorhodospira involved in sulfur cycling. Members of family Desulfohalobiaceae (Deltaproteobacteria) dealing with the reductive part of the sulfur cycle were present in less abundance. These are reported to be more active at lower salinity (Lane, 1991;Jackson et al., 2014). Among Gammaproteobacteria, Thiohalorhabdus and members from order Chromatiales such as Halorhodospira grow by autotrophically utilizing energy derived from the oxidation of elemental sulfur and reduced inorganic sulfur compounds. The presence of both sulfur-reducing and sulfur-oxidizing prokaryotes indicates the presence of an active sulfur cycle, especially in the salt samples. The presence of the nitrogen cycle was also detected with the presence of OTUs affiliated with efficient ammoniaoxidizing bacteria (Beta-proteobacteria) Nitrosomonas. These chemolithotrophic ammonia-oxidizing bacteria AOB commonly belong to the Beta-and Gammaproteobacteria (Purkhold et al., 2000). Phosphate solubilizing archaea, Haloarcula, Halobacterium, and Halolamina were observed in our dataset. These organisms had implicated involvement in the maintenance of Phosphorous cycle in Hypersaline environments (Yadav et al., 2015). Marinobacter representing the major genera in Karak samples was also reported to have incredible metabolic potential and had shown capabilities of anaerobic respiration using dimethylsulfoxide (DMSO) or nitrate (Yau et al., 2013).

Diversity Indices
Alpha diversity was measured using metrics such a Good's coverage, Shannon index and specie richness were calculated with a nonparametric Chao1 estimator ( Table 2). With a number of 3956 OTUs the salt samples of Karak mine (KKS) displayed the more diverse habitat whereas the brine KKB with 2215 OTUs showed a somewhat lower diversity. The diversity metrics analysis indicated that between the two samples, KKS (salt) had the highest diversity as compared to KKS (brine). The results of Good's coverage (99% coverage), and rarefaction curves (Figure 3) suggest an almost complete sampling. Perhaps, the higher abundance of different minerals in halite salt as compared to water samples acts as an abiotic regulator on microbial community composition (Paul et al., 2016). The statistics revealed a high proportion of co-shared phylotypes (635) between the salt and brine samples, signifying the presence of a well-adapted and stable halophilic core community distribution. However, unique OTUs in salt samples such as Acetothmeria and Planctomycetes highlights their specific role in the survival and functioning of microbial community entrapped within salt crystals.

CONCLUSION
The current examination provides first-ever insight into the composition of halophilic microbial communities from an uncharacterized ecosystem of Karak salt mine, Pakistan. The significant existence of halophiles advocates the capacity of microbes to grow and divide in this extreme environment. The distribution of archaeal and bacterial communities in the Karak salt mines is diverse, with the members of Phyla Euryarchaeota, Bacteroidetes, Proteobacteria, Firmicutes, and Class Acetothermia as top 5 phyla. Representatives of the genera Marinobacter and Salinibacter were the principal inhabitants in samples suggesting their crucial role in this environment. In addition to the adaptation to hypersalinity, the perchlorateresistant ability of bacteria (Marinobacter, Halomonas) and archaea (Haloarcula, Haloferax) dominating the salt mine, make this environment and isolates the best model for astrobiological studies. Although the uncharacterized OTUs covered a relatively large proportion of the overall population in our analysis, the number of novel genera spotted revealed significant diversity, suggesting that a more detailed examination needs to be performed to characterize them. Future research with multiple sampling points, e.g., halite crystals from multiple-layered sediments and at different time periods within the mine will increase our understanding of the ecology of extreme halophiles. Our future prospects include a detailed analysis of culturable diversity for a more comprehensive snapshot of microbial diversity. The metagenome analysis cannot calculate protein synthesis and gene expression. Therefore, in order to investigate the functional complexity in microbial ecosystems, a metatranscriptomic method is essential. Analogous halites on Mars would signify an excellent setting to investigate for conserved organic material and associated biomarkers, which might be concentrated through evaporative processes and subsequently conserve life that could persist over geologic time.

AUTHOR CONTRIBUTIONS
LC drafted the manuscript, performed all laboratory experiments, did downstream analysis with the results, and performed bioinformatics analysis of the raw RNA-seq data using the transcriptomic pipeline. FH and SD provided experimental guidance and manuscript feedback and editing. RM, MA, and WP provided guidance and critical feedback and editing of the manuscript. All authors contributed to the article and approved the submitted version.