Microbiome Fingerprint as Biomarker for Geographical Origin and Heredity in Crocus sativus: A Feasibility Study

Host–microbiome interactions are specific and not random, making them defining entities for the host. The hypothesis proposed by various researchers earlier, that both plants and animals harbor specific inheritable core microbiome, is being augmented in the present study. Additionally, a case for using microbial fingerprint as a biomarker, not only for plant identification but also as a geographical indicator, has been investigated, taking Crocus sativus, saffron, as a study material. Crocus sativus, a monogenetic herb, on account of its male sterility and vegetative propagation, is reported to lack genome based molecular markers. Cormosphere microbiome (microbiome associated with corm) has been compared across three geographical locations, in two continents, to identify the core and unique microbiome, during the vegetative phase of its growth. Microbiome analysis done at phylum and genus level, using next generation sequencing technology, revealed that cormosphere at three locations harbored common phyla. At genus level, 24 genera were found common to all three geographical locations, indicating them to be part of the core microbiome of saffron. However, there were some bacterial genera unique to Kashmir, Kishtwar, and Morocco that can be used to develop microbial markers/geographical indicators for saffron grown in these regions. This is a preliminary study, indicating that the location specific bacterial community can be used to develop microbial barcodes but needs further augmentation with high coverage data from other saffron growing geographical regions.


INTRODUCTION
Molecular markers are important tools for plant genome analysis, crop improvement, and development of barcodes for authentic plant identification (Mishra et al., 2016). The development of molecular markers relies on genetic variation in the plants; however, plants lacking genetic variation due to asexual reproduction harbor unidentified or poor resolution molecular markers. Crocus sativus is one amongst such plants that lack genetic variations. Crocus sativus, commonly known as saffron, is an autumn-flowering perennial sterile triploid plant bearing eight chromosomes per set, i.e., 2n = 3x = 24 (Goldblatt et al., 2006;Ahmad et al., 2021). Due to male sterility, it does not produce viable seeds and propagates vegetatively, thereby lacking any variations at the genomic level. In order to propagate vegetatively, clusters of corms, underground bulb-like, starch-storing organs are dug up, divided, and replanted in the soil (Golmohammadi, 2014;Moshtaghi, 2020). The stigmas of saffron flowers are plucked and dried, to be used as a seasoning and coloring agent in food (Mzabri et al., 2017;Jafari et al., 2020). It is the world's highest priced spice with medicinal and aromatic properties, so referred to as 'golden condiment' (Monika and Neha, 2014;Pandita, 2021;Su et al., 2021). This exotic spice is cultivated worldwide in many countries, among them Iran, Spain, India and Morocco are the major producers (Mykhailenko et al., 2020;Gupta et al., 2021). Saffron is reported to be a monomorphic plant, having no variations at a genomic level across different saffron accessions across the world, as depicted by various molecular markers like RAPD, AFLP, and SSR (Rubio-Moraga et al., 2009;Busconi et al., 2015). However, there are variations in the phenotypic and biochemical characteristics such as the percentage of various metabolites (crocin, picocrocin, and saffranal) reported from different geographical locations (Othman et al., 2020). However, there are no authentic molecular markers that can identify these variations till date (Alavi-Kia et al., 2008;Keify and Beiki, 2012;Mir et al., 2021).
In the last three decades, microbes are established to have a pronounced effect on the biotic or abiotic environments they inhabit (Gupta et al., 2021). No eukaryote lives as an individual but is a meta-organism or holobiont (Guerrero et al., 2013;Vandenkoornhuyse et al., 2015). As established in humans, plants have their own microbiota and studies have shown that the rhizosphere (hotspot for microbial activities) has tremendous microbial density associated with almost 10 11 microbial cells/gram of soil, represented by >30,000 prokaryotic species (Berendsen et al., 2012;Chaparro et al., 2014;Ofek-Lalzar et al., 2014). The sum total of the genomes of microbes associated with any environmental niche is referred to as "microbiome", the term coined by Joshua Lederberg (Lederberg and McCray, 2001).
The rhizosphere of many plants has been reported to contain specific microbiomes that are essential for plant survival and are referred to as the plant's second genome (Berendsen et al., 2012;Ofek-Lalzar et al., 2014). The microbiome and its dynamics have been extensively studied in the rhizosphere of Arabidopsis, Zea mays, Triticum aestivum, Oryza sativa, Vitis vinifera, Crocus sativus, etc. (Ambardar and Vakhlu, 2013;Ambardar et al., 2014Ambardar et al., , 2016Mendes and Raaijmakers, 2015;Kandel et al., 2017). Studies also suggest that though the microbiome associated with any plant is organ specific, it may show temporal variation, primarily during developmental stages, geographical location, and cultivation practices (Peiffer et al., 2013;Chaparro et al., 2014;Edwards et al., 2015;Lin et al., 2020).
Despite temporal and spatial variations, the core microbiome is specific for each plant. Core microbiome may be defined as a congregation of micro-organisms within a host's microbiome that remains unaltered across different growth stages, niches, and geographical locations. Core inheritable microbiome, present in all the organisms, is hypothesized to be inherited as a core genome (Walters et al., 2018). The composition of core microbiome is reported to depend on the host plant, and this hypothesis was tested in Arabidopsis by Lundberg et al. (2012) wherein specific microbial communities common to tested plant varieties were established. The presence of core microbiome has been established in rice (Eyre et al., 2019), sponges (Schmitt et al., 2012), corals (Ainsworth et al., 2015;Hernandez-Agreda et al., 2016, 2017, wheat (Kuzniar et al., 2020), human gut (Qin et al., 2010), human hands (Fierer et al., 2008), and sand beaches (Newton et al., 2013).
The present study was initiated to discover the core microbiome of saffron cormosphere, along with unique genera specific to the saffron cultivated at a particular location. Unique genera can be used to develop molecular markers/geographical indicators. In the present study, core microbiome of saffron cormosphere comprises 24 genera that remain unaltered in three locations across two continents, i.e., Africa and Asia. In addition, unique microbes associated with the corm at each location have also been identified and cataloged, which can be used as molecular markers/geographical indicators.

Sample Collection
Plant samples were collected from three geographical locations, two in Jammu and Kashmir, India, i.e., Kishtwar (33 • 19 ′ 12.00 ′′ N 75 • 46 ′ 12.00 ′′ E) and Kashmir (34.02 •′′ N 74.930 ′′ E), and one from Morocco (Taliouine, 30 • 31 ′ 58.00 ′′ N 7 • 55 ′ 32.00 ′′ W) (Figure 1). Plant samples were collected from Kishtwar during three consecutive years (2011, 2012, and 2013; coded as C11, C12, and C13, respectively), whereas from Kashmir in 2013 and Morocco in 2016. In all the cases, sampling was done during the vegetative phase of saffron life cycle. The sample collection was done as per the protocol developed by Ambardar et al. (2014). The temperature and humidity of the sites at the time of sample collection have been tabulated in Table 1.
In a separate experiment, corm samples and saffron field soil from Kishtwar region [referred to as Kishtwar-L 17 (L for lab)] were collected during the dormant phase of life cycle (August, 2017) and stored at room temperature (20-25 • C) till further use. Corms were sown in the soil (Kishtwar soil) at the start of vegetative phase (October, 2017) and were grown for 2 months. The corms were grown in pots that were incubated in plant growth chamber with day temperature (25 ± 5 • C), night temperature (10 ± 5 • C), and 70% relative humidity and watered after every 5 days with autoclaved distilled water for a period of 2 months. Cormosphere samples were collected after 2 months in vegetative phase.

Whole-Genome Shotgun Sequencing and Bioinformatic Analysis
The library preparation of isolated metagenomic DNA from all the samples was done using Truseq Nano DNA Library preparation kit (Catalog No. 20015964, Illumina, CA, USA) as per the instructions of manufacturer. To ensure maximum yield, a high-fidelity amplification step was performed using sparQ HiFi PCR Master mix (Quanta bio, QIAGEN, Beverly Inc. Catalog No. 95192-050). The quality of amplified libraries was analyzed on Bioanalyzer 2100 (Agilent Technologies, catalog number: G2939BA) and was further sequenced using Illumina NextSeq 500 platform with 2 × 150 base paired-end configuration from Xcleries labs limited, Ahmadabad, Gujarat, India. The raw reads were quality checked using FastQC tool kit (Brown et al., 2017). The low quality reads having phred score less than 30 were filtered and trimmed using software Trimmomatic (Bolger et al., 2014). The de novo assembly of high quality paired end reads for all samples was done using CLC Genomics Workbench 6.0 (CLC bio, Aarhus, Denmark) with the following parameters: minimum contig length: 200, automatic word size, perform scaffolding, mismatch cost: 2, insertion cost: 3, deletion cost: 3, length fraction: 0.5, similarity fraction: 0.8. The assembled contigs were then uploaded on MG-RAST (Metagenomic Rapid Annotation using Subsystem Technology) server version 4.0.3 with default parameters for analyzing taxonomic hit distribution and calculation of Shannon-alpha diversity values (Meyer et al., 2008). Taxonomic classification up to phylum and genus level was represented in the form of heat maps and principal component analysis (PCA) plots using ClustVis tool (Metsalu and Vilo, 2015). Rarefaction curve was generated using PAST software version 4.03 (Hammer et al., 2001). Circos plot was generated using the tool Circos Version 0.63-9 (Krzywinski et al., 2009).

RESULTS
Whole metagenome sequencing data of the six samples (Kishtwar 2011(Kishtwar , 2012(Kishtwar , 2013Kishtwar-L 17, Kashmir 2013;Morocco 2016), collected from three locations during the vegetative period of growth, ranged between 5.3 million reads and 40 million reads. Reads mapped to bacteria and eukaryota ranged from 85.29 to 98.78% and 1 to 14.2% in all the samples, respectively ( Table 2). The dynamics of the microbiome of Kishtwar for three consecutive years (2011, 2012, and 2013) were studied, and the average of the 3 years was taken as reference. Subsequently, the Kishtwar average was compared to the microbiome of saffron cormosphere of Kishtwar-L 17, Kashmir (2013), and Morocco (2016) to discover core and location specific unique microbiome.

Structure and Inheritability of Kishtwar Cormosphere Microbiome for the Three Consecutive Years
At the phylum level, 15 phyla (12 bacterial, 2 fungal, 1 archaeal) were cataloged from the cormosphere of saffron grown in Kishtwar during three consecutive years (2011, 2012, and 2013). The three most dominant phyla were Proteobacteria followed by Actinobacteria and Bacteroidetes. The relative abundance of Proteobacteria, Actinobacteria, and Bacteroidetes varied among the samples, and it was observed that the percentage of Proteobacteria reduced gradually from 2011 to 2013 (54.4 to On the other hand, Actinobacteria and Bacteroidetes increased from 29.5 to 40.9% and 3.6 to 4.36%, respectively (Figure 2). At the genus level, a total 52 genera were cataloged from three Kishtwar samples. Surprisingly, Streptomyces (C11 = 11.5%, C12 = 13.8%, and C13 = 17.7%) belonging to phylum Actinobacteria was the most abundant followed by Mycobacterium, Bradyrhizobium, and Burkholderia all belonging to Proteobacteria (Figure 3). The abundance of three genera Streptomyces, Burkholderia, and Mycobacterium increased from C11 to C13, whereas Bradyrhizobium decreased from 4.58% to 3.22%. Out of total 52 genera cataloged, 42 genera were found common in C11, C12, and C13, which correspond to 81% of total genera. The relative abundance of the 10 genera that were not common (19%) was found to be less than 1% of each genus. Core Bacterial Microbiome of C11, C12, C13 Core bacteriome of the saffron plant grown in Kishtwar comprised of 42 genera that belongs to five different phyla and corresponds to 81% of the total microbiome (52 genera). Based on diversity/relative representation, 50% (21 genera) belonged to Proteobacteria, 31% (13 genera) to Actinobacteria, 7.2% (3 genera) to Bacteriodetes, 7.2% (3 genera) to Verrumicrobia, and 2.3% (1 genus) to phylum Planctomycetes in the core microbiome of Kishtwar. However, on the basis of abundance, 10 most abundant genera were Streptomyces followed by Mycobacterium, Bradyrhizobium, Burkholderia, Mesorhizobium, Conexibacter, Sorangium, Agrobacterium, Rhizobium, and Pseudomonas. The relative abundance of core genera has been tabulated in Table 3.
In addition to 42 common genera, 10 genera were randomly present in one or two of the three samples. Two genera (Stenotrophomonas and Anaeromyxobacter) belonging to phylum Proteobacteria were present only in C11. Three genera such as Bacillus (Firmicutes), Mucilaginibacte (Bacteriodetes), and Brevundimonas (Proteobacteria) were present only in C13. The Shannon α-diversity of C11, C12, and C13 was found to be 443, 395, and 359, respectively, thereby indicating decrease in microbial diversity from 2011 to 2013. However, decrease in microbial diversity may be due to the depletion of few rare genera having less abundance in consecutive years. To conclude, it was observed that 81% of genera were common in the three consecutive years but their relative abundance varied.

Comparative Stability and Inheritability of the Corm Microbiome in Kishtwar Average and Kishtwar-L 17
The Kishtwar microbiome (Kishtwar average of C11, C12, and C13, hereafter referred to as Kishtwar-A) was taken as reference microbiome for comparing the diversity and abundance of microbes in Kishtwar-L 17 (plants grown under controlled laboratory conditions) again to confirm stability and inheritability of the corm microbiome. All 15 phyla present in Kishtwar with 3 year average, referred to as Kishtwar-A, were present in Kashmir-L 17. The three major dominant phyla were the same, i.e., Proteobacteria, Actinobacteria, and Bacteroidetes, but their relative abundance varied between Kishtwar-L 17 and Kishtwar-A. Proteobacteria and Actinobacteria were less abundant, whereas Bacteroidetes were more abundant in Kishtwar-L 17 (Figure 4).

Diversity and Structure of Cormosphere Microbiome of Kashmir and Morocco in Comparison to Reference Kishtwar-A Cormosphere Microbiome
At the phylum level, 15 phyla (12 bacterial, 2 fungal, 1 archaeal) were cataloged from the cormosphere of saffron corms grown in Kashmir and Morocco. The three dominant bacterial phyla among the samples were Proteobacteria, Bacteroidetes, and Actinobacteria. Proteobacteria was most dominant (Kashmir = 48.46%, Morocco = 55.06%) followed by Bacteroidetes (Kashmir = 12.93% and Morocco = 16.92%) and Actinobacteria (Kashmir = 10.98% and Morocco = 16.34%) Proteobacteria (51.23%), followed by Actinobacteria (35.64%), Bacteroidetes (4.06%), and Verrucomicrobia (1.85%). Though the number of phyla was same in Kishtwar-A, Kashmir, and Morocco, heat maps hierarchical clustering clustered Kashmir and Morocco in one clad and Kishtwar-A in different clad. This represented that the co-relation distance and average linkage of Kishtwar-A were different from Kashmir and Morocco (Figure 7). Relative abundance of different phyla of Kishtwar-A, Kashmir, and Morocco is represented in Figure 8.

Comparative Analysis of Core Bacteriome of Three Different Geographical Locations
Of the total 73 genera cataloged, 24 genera were found common in all the three geographical locations, viz., Kishtwar, Kashmir, and Morocco, corresponding to 32.8% similarity in microbial community across geographical regions in Asian and African sub-continent. Core microbiome of saffron across different geographical location comprised of 24 genera; out of 24 common genera, 13 genera (54.2%) belonged to phylum Proteobacteria, 5 genera (20.8) to Actinobacteria, 3 genera (12.5%) to Bacteriodetes, and 3 genera (12.5%) to Verrumicrobia (Figure 10). The core genera with their relative abundance have been tabulated in Table 4.
Notable differences were observed in the cormosphere microbiome composition of Kashmir, Kishtwar, and Morocco as depicted by principal component analysis (PCA) that clusters Kishtwar (C11, C12, and C13) samples in close proximity and Kashmir and Morocco samples that clustered away in different coordinate, thereby indicating differences in diversity/abundance pattern (Figure 11). In addition, alpha diversity indices like Shannon indices were highest in case of Kashmir (507) as compared to Morocco (415) and Kishtwar-A (399). This was further complemented by rarefaction curves wherein Kashmir samples were more enriched as depicted by the higher number of species in Kashmir, in comparison to Kishtwar-A and Morocco (Figure 12).
To conclude, Kishtwar consecutive 3 year samples showed 81% similar genera across the years. The genera that were absent in one or other Kishtwar samples had less than 1% representation in total reads and may be attributed to sequencing bias due to the low relative amount present in whole metagenomic data. This needs to be confirmed by repeated sampling and sequencing across various geographical locations. However, 32.8% of the genera were common in all three locations that can be designated as core microbiome of saffron cormosphere. Bacterial genera that were present in one geographical location and completely absent in the other two locations were observed, thereby representing the specific interactions of these bacterial genera to each geographical location. Unique genera to each place can be developed into biological markers as geographical indicators.

DATA AVAILABILITY STATEMENT
The sequencing data has been submitted to NCBI under the accession number PRJNA705068.

DISCUSSION
The present study was initiated to investigate the possibility of using microbial fingerprints as molecular marker and/or geographical indicator, particularly in monogenetic crops. Microbial biomarkers have been reported to be used for the determination of biogenic chemical concentration (e.g., phospholipid fatty acids; Lupwayi et al., 2017), evaluation of microbial activity (e.g., photosynthetic activity; Corcoll et al., 2011), and capacity (e.g., adaptation; Pesce et al., 2016). Species specific biomarkers are reported to detect the toxicity in the environment caused by a pathogen bacteria or toxic Cyanobacteria (Bonnineau et al., 2012;Guasch et al., 2017). In plants also, microbiome is reported to have significant variation based on different geographical locations, such as reported in the case of vineyards, Amazonian hardwood, and fragrant orchid. The researchers co-relate the variations in the plant microbiome to various factors such as location/genetics/abiotic and biotic factors (Coller et al., 2019;Skaltsas et al., 2019;Lin et al., 2020). However, reports are emerging on the core microbiome that remains unaltered at various developmental stages of the plants or plants from different locations/niches (Chen et al., 2018;Toju et al., 2018;Lin et al., 2020). The significance of molecular markers/barcodes in selecting/identifying different varieties needs no elaboration.
However, some plants are reported to lack high resolution molecular markers on account of sterility. Crocus sativus, saffron, is reported to lack genetic variation on account of vegetative propagation through corms, as it is male sterile plant (Nemati et al., 2019). The present study was initiated to discover if microbiome based markers can be used as geographical indicators in saffron grown at different geographical locations, i.e., Kashmir and Kishtwar in India and Taliouine in Morocco, also to identify core microbiome in saffron across different locations, if any. The microbiome associated with corm, grown at these three locations, was isolated and analyzed at the vegetative stage to keep the development stage as a constant factor.
The bacteriome and mycobiome of Crocus sativus have been reported earlier by our group, but that study was specific to Kashmir for 1 year only . In the present study, microbiome analysis of cormosphere was initiated in saffron grown in Kishtwar during three consecutive years to check the inheritability of cormosphere microbiome. About 81% of the microbial genera were common in the cormosphere of all 3 years, though their abundance varied (Figure 2). The abundance of the genera belonging to phylum Proteobacteria decreased from C11 (54.4%) to C13 (48.3%), whereas an increase in the FIGURE 9 | Relative abundance of top 25 genera cataloged from cormosphere of all the three geographical locations (a total of 73 genera were cataloged from three geographical locations: Kishtwar-A, Kashmir, and Morocco). Streptomyces was the dominant genera in Kishtwar-A followed by Mycobacterim, Bradyrhizobium, and Burkholderia. Pseudomonas was the dominant genera in both Kashmir and Morocco followed by Streptomyces, Chitinophaga, and Sorangium in Kashmir and Caulobacter, Pedobacter, and Flavobacterium in Morocco.
abundance of Actinobacteria and Bacteriodetes was observed from C11 (Actinobacteria−29.5%, Bacteriodetes−3.6%) to C13 (Actinobacteria−40.9%, Bacteriodetes−4.36%). Since saffron is a perennial crop, and farming practices were not changed, the climate was more or less the same and yield also did not show any co-relation (Table 1); the exact reason for this variation in abundance is not clear. Eighty-one percent of the genera, common to all the 3 years, are being proposed as the core microbiome of saffron cormosphere grown in Kishtwar. Core microbiome have been earlier reported in various plants across different geographical location, climatic gradients, genotypes, and seasons, namely, Salvia miltiorrhiza (Chen et al., 2018), angiosperm plants (Fitzpatrick et al., 2018), Coffea arabica (Fulthorpe et al., 2020), Orzya sativa (Eyre et al., 2019), Phaseolus vulgaris (Pérez-Jaramillo et al., 2019), Switchgrass (Panicum virgatum L.) and Miscanthus (Miscanthus x giganteus) (Grady et al., 2019), and Panax ginsen, (Nguyen et al., 2016). On the similar lines, phyllosphere microbiome of Switchgrass was compared for two consecutive years (2016 and 2017) and 60% of leaf communities were found common (Grady et al., 2019), which was comparatively less than the present study (81% of cormosphere community was common in saffron).
In the present study, once the inheritability of the core microbiome in field conditions was established, the same was tested under laboratory conditions wherein the saffron corms and soil were collected from Kishtwar in 2017, and grown in pots in the laboratory. Except for the variation of temperature that was around 25 • C, the rest of the parameters were kept close to natural conditions, such as soil and light intensity. Out of 81% genera common during three consecutive years from Kishtwar, only 65% FIGURE 10 | Abundance of core genera and unique genera in three geographical locations. Central line (ribbons) represents the abundance of each genus and connects genera to samples. The outer three rings are the stacked bars representing the relative contribution of particular genera in every sample.
were retained under laboratory conditions. Eight genera were missing under lab conditions, but additional six genera were present and were completely absent under field conditions. The aim of laboratory cultivation was to check if the unique genera present in Kishtwar samples are retained if corms are dug and grown somewhere else. It was observed that the unique genera were retained as the five unique genera present in Kishtwar 3 years (C11, C12, and C13) were also present in these samples. This further establishes that the unique genera are inherited and can be used as geographical indicators or markers. In addition, all the microbes, representing the core microbiome of saffron corm characterized from different locations, were also present in laboratory grown samples (discussed in the following sections).
In total, 15 phyla (12 bacterial, 2 fungal, 1 archaeal) cataloged were common among three different locations. The most dominant phylum was Proteobacteria followed by Actinobacteria, Bacteroidetes, Streptophyta, Acidobacteria, and Verrucomicrobia (Figure 8). Proteobacteria have also been reported to be the dominant phylum in the rhizosphere of other plants such as Gossypium hirsutum, Artemisia argyi, Ageratum conyzoides, Erigeron annuus, Bidens biternata, Euphorbia hirta, and Viola japonica (Qiao et al., 2017;Lei et al., 2019). Proteobacteria are the main microbial contributors to plant growth promoting rhizobacteria (PGPR), reported so far. They are fast growing and play a role in diverse metabolic activities that enhance plant growth (Philippot et al., 2013;Rampelotto et al., 2013;Schillaci et al., 2019;Mukhtar et al., 2020). Actinobacteria are the gram positive bacteria that are known to promote plant growth and having biocontrol activity (Viaene et al., 2016). Likewise, Bacteriodetes are the pathogen suppressing members and also have a role in plant growth promotion by phosphorus mobilization (Lidbury et al., 2021). This explains the dominance of Proteobacteria, Actinobacteria, and Bacteroidetes in the cormosphere, which are among the most dominant phyla in the soil and maybe getting transferred from the soil, just like in rhizosphere (Miyashita, 2015;Fierer, 2017;Compant et al., 2019). To our knowledge, cormosphere microbiome has not been studied in any other plant, except for saffron reported previously by our group (Ambardar and Vakhlu, 2013;Ambardar, 2014;Ambardar et al., 2014Ambardar et al., , 2016Kour et al., 2018) from Kashmir, India (Chamkhi et al., 2018), and from Taliouine, Morocco. Previous reports by our group have also confirmed the dominance of Proteobacteria in Kashmir saffron rhizosphere and cormosphere, even though the method used was cultivation dependent and 16S targeted metagenome sequencing was different than the methods used in the present study (Ambardar and Vakhlu, 2013;Ambardar et al., 2014) (Ambardar and Vakhlu, 2013;Ambardar et al., 2014;Chen et al., 2018;Kour et al., 2018;Chu et al., 2019;Qessaoui et al., 2019). Pseudomonas isolates have been reported to solubilize phosphate, produce siderophores, ammonia, and indole-3-acetic acid, and colonize the roots of tomato plants (Qessaoui et al., 2019) and also have been reported to induce salt tolerance in Arabidopsis thaliana (Chu et al., 2019). Pseudomonas has also been reported to be dominant in saffron rhizosphere grown in Kashmir, India, analyzed using both cultivation dependent and cultivation independent method (Ambardar and Vakhlu, 2013;Ambardar et al., 2014), and Morocco, analyzed using cultivation dependent approaches (Chamkhi et al., 2018). Streptomyces is the largest genus of Actinobacteria phylum comprising 500 species and has been found to promote plant growth (Lehr et al., 2007(Lehr et al., , 2008Schrey and Tarkka, 2008;Viaene et al., 2016). It is reported to produce secondary metabolites having antibiotic activities and volatile organic compounds (VOCs) that can exert biocontrol activity (Lucas et al., 2012;Viaene et al., 2016).
Streptomyces have been previously reported from the saffron rhizosphere also using cultivation independent approaches . In addition, to Pseudomonas and Streptomyces, 53 have also been reported previously in saffron in Kashmir and Morocco (Ambardar and Vakhlu, 2013;Ambardar et al., 2014Ambardar et al., , 2016Chamkhi et al., 2018;Kour et al., 2018) (Supplementary Material 2). Interestingly, we found that 24 genera were common to all, representing 32.8% of collective bacteriome in the four samples (three field and one laboratory sample). These 24 genera are being proposed to represent the core microbiome of cormosphere of C. sativus across different geographical locations. Core microbiome across different geographical locations has been earlier reported in various plants such as Salvia miltiorrhiza (Chen et al., 2018), Phaseolus vulgaris (Pérez-Jaramillo et al., 2019), Oryza sativa (Eyre et al., 2019), and Gymnadenia conopsea (Lin et al., 2020) and also in Vineyards soil (Coller et al., 2019) and Dalbergia spruceana (Skaltsas et al., 2019). Core microbiome of 21 Salvia miltiorrhiza seeds represented 54% of whole microbiome cataloged from seven different geographic origins, and Pantoea, Pseudomonas, and Sphingomonas were enriched core taxa (Chen et al., 2018). In the case of Phaseolus vulgaris rhizosphere, core microbiome represented 25.9% of total microbiome in native and agricultural soils (Pérez-Jaramillo et al., 2019). Core microbiome of 30 phylogenetically diverse angiosperm plants constitutes 40% of whole microbiome (Fitzpatrick et al., 2018). In rice seeds, the core microbiome was enriched in Rhizobium, Pantoea, Sphingomonas, Methylobacterium, Xanthomonas, Paenibacillus, Alternaria, and Occultifur that have been also reported as PGPB (plant growth promoting bacteria) . Core microbiome in saffron was 32.8% of the total microbiome, which was comparatively less than Salvia miltiorrhiza. Similar to core microbiome in rice and Salvia miltiorrhiza (seed), saffron's core microbiome was also enriched with Pseudomonas, Rhizobium, Sphingomonas, Xanthomonas, and Methylobacterium. In addition, 22 bacterial genera (out of 24) representing the core microbiome of the C. sativus in the present study were also reported previously in Saffron grown in Kashmir, India, by our group, out of which 19 are PGPB in saffron and other plants (Ambardar and Vakhlu, 2013;Kour et al., 2018;Umadevi et al., 2018;Vergani et al., 2019) (Supplementary Material 3). Agrobacterium, Bacillus, Delftia, Luteibacter, Pantoea, Pseudomonas, Rahnella, Rhizobium, and Variovorax from rhizosphere and Agrobacterium, Pseudomonas, Rhizobium, and Variovorax from cormosphere of saffron grown in Morocco using cultivation dependent methods have been reported earlier as well (Chamkhi et al., 2018). Core microbiome enriched with plant growth promoting bacteria indicates that the plant selects the core microbiome based on their function (Supplementary Material 3). To conclude, the saffron cormosphere core microbiome is conserved across two continents, Asia and Africa, as per this study.

Microbial Genera Specific to Geographical Location
Variation in the microbiome of C. sativus across different geographical locations was confirmed by PCA plots, wherein surprisingly Morocco and Kashmir were clustered away from Kishtwar-A, although Kashmir is geographically closer to Kishtwar (Figure 11). These findings were further complemented by Shannon alpha diversity indices, which were maximum in Kashmir (507) compared to Kishtwar-A (354) and Morocco (415). The Shannon alpha diversity index indicates that microbial communities of each sample have significant differences, though they share core microbiome. There are plenty of reports on the effect of temperature on the diversity and abundance of soil microbes, though those studies have been conducted in reference to climate change (Jansson and Hofmockel, 2020).
The effect of geographical location on the microbiome has been discovered in other plants such as barley (Terrazas et al., 2020) and Gymnadenia conopsea (Lin et al., 2020), and significant variation has been reported. Each location in the present study harbored unique bacteria that could be used as a geographical indicator. The cormosphere of Kashmir had five unique bacteria, namely, Marivirga, Planctomyces, Plesiocystis, Candidatus koribacter, and Microscilla, whereas seven bacteria, namely, Rahnella, Chryseobacterium, Bdellovibrio, Algoriphagus, Methylovorus, Geodermatophilus, and Aeromicrobium, were specific to Morocco and five (Rhodococcus, Kribella, Streptosporangium, Saccharopolyspora, and Amycolatopsis) were specific to Kishtwar (Figure 13). These unique genera have also been reported as PGPB in other plants (Supplementary Material 4).
The bacterial genera (8) that have been reported for the first time from saffron underground parts were Marivirga, Microscilla, Starkeya, Actinosynnema, Saccharopolyspora, Bdellovibrio, Algoriphagus, and Methylovorus, whereas 11 bacterial genera (Planctomyces, Plesiocystis, Candidatus koribacter, Chryseobacterium, Rhodococcus, Rahnella FIGURE 13 | Representing the number of common and unique genera identified in three geographical locations. Twenty-four common genera in all the three geographical locations are referred as the core microbiome of saffron. Five genera were unique to Kishtwar, seven to Morocco, and five to Kashmir. Amycolatopsis, Kribbella, Streptosporangium, Geodermatophilus, and Aeromicrobium) have been reported earlier from saffron rhizosphere in our previous studies from Kashmir and Kishtwar in India and Taliouine in Morocco (Ambardar et al., 2021 unpublished data;Chamkhi et al., 2018). However, the bacteria reported from saffron for the first time have been reported from other plants. Marivirga has been reported from Panax ginseng using 16S rRNA sequencing (Nguyen et al., 2016), Microscilla from Tuber melanosporum (Deveau et al., 2016), Starkeya from root nodules of legume plants (Zakhia et al., 2006), Actinosynnema from beech rhizosphere (Colin et al., 2017), Saccharopolyspora from Colocasia rhizosphere (Intra et al., 2019), Rahnella from Oryza sativa (Sun et al., 2020), Chryseobacterium from sand dune plant rhizosphere (Cho et al., 2010), Algoriphagus from mangrove rhizosphere (Song et al., 2020), and Bdellovibrio and Methylovorus from grassland (Macey et al., 2020). Marivirga, Microscilla, Actinosynnema, and Bdellovibrio have been identified using cultivation independent approach (16S rRNA sequencing), whereas Saccharopolyspora, Chryseobacterium, and Algoriphagus have been identified using cultivation based approaches but their PGP activities have not been estimated yet. However, Starkeya and Rahnella have been reported as PGPRs and isolated using cultivation dependent approaches, whereas Methylovorus is reported to be methylotrophs (Zakhia et al., 2006;Cho et al., 2010;Deveau et al., 2016;Nguyen et al., 2016;Colin et al., 2017;Intra et al., 2019;Macey et al., 2020;Song et al., 2020;Sun et al., 2020). The preliminary study suggests that saffron across three geographical locations share common core cormosphere microbiome but have certain unique bacteria genera specific for each location. The unique genera specific to a particular location can be used as the geographical indicators and molecular markers in saffron. The study needs to be complemented further by analyzing the microbiome of saffron grown across various geographical locations, as saffron is cultivated now in both the hemispheres. It also needs to be validated by in-depth sequencing and frequent sampling for a couple of consecutive years.

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

AUTHOR CONTRIBUTIONS
JV, MH, and LS gave plan of study. SS, NB, DT, RZ, and JV designed the experiments. SS, NB, SA, and DT conducted the experiments and analyzed the sequencing data. SS, SA, NB, and SR wrote the manuscript and analyzed the taxonomy data. JV, MH, and LS finalized the manuscript. All authors contributed to the article and approved the submitted version.

ACKNOWLEDGMENTS
Sincere thanks and gratitude to Tahir Ali (Ph.D. Scholar, School of Biotechnology, University of Jammu, Jammu, J&K, India), helping in preparing map showing sample collection sites on the world map and to Mr. Ajmad Malik (Saffron development officer) from Kishtwar (J&K) for providing data on saffron yield. Authors are thankful to UGC for funding the project UGC [42-168/2013 (SR)] and Indo-Austrian project of OeAD GmbH and DST under call number IN2022 to carry out this study.