Metagenomic Analysis of Lichen Associated Bacterial Community Profiling In Roccella Montagnei

doi:10.21203/rs.3.rs-915860/v1

PDF

Research Article

Metagenomic Analysis of Lichen Associated Bacterial Community Profiling In Roccella Montagnei

https://doi.org/10.21203/rs.3.rs-915860/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 22 Dec, 2021

Read the published version in Archives of Microbiology →

You are reading this latest preprint version

A lichen is a composite organism formed of algae or cyanobacteria that live in a mutually advantageous symbiotic relationship with the filaments (hyphae) of fungus. Three lichen samples were obtained from diverse sites at Kuppanasamy temple, Pollachi, a terrestrial habitat located in Coimbatore, Nithiravilai, Nagercoil, and Ramarpatham, Vedaranyam, both coastal habitats located in Kanyakumari and Nagapattinam districts of Tamil Nadu, for this study. Amplification and sequencing of 16S rRNA V3 - V4 regions were used for metagenomic study. Aside from the NGS data, distinct types of lichen microbiome profiles were clearly shown. The bacterial diversity in the lichen genera of Roccella montagnei growing in coastal and terrestrial environments was further investigated using common and unique operational taxonomic units (OTUs) and the QIIME pipeline (1.9.1). By using similarity clustering, the heat map analysis depicts the abundance information of chosen OTUs as well as the similarity and difference between OTUs and lichen samples. Using multiple methods, the Alpha and Beta diversity analysis revealed that there were differences in all of the samples. However, UPGMA tree inference of comparable bacterial community in coastal habitat lichen samples compared to terrestrial habitat validates their evolutionary lineage. As a result, the bacterial population associated with corticolous lichen is dependent on geographic locations, growth substrate, and climatic circumstances of similar lichen genera produced in different habitats and tree substrates.

Molecular Biology

Lichen

Roccella montagnei

microbiome

16S rRNA V3-V4 regions

OTUs

QIIME

Lichens are symbiotic organisms, which were formed by the mycobiont and photobiont partner. Lichens have recently been documented as advanced ecosystems that contain a distinct range of microbial lineages from all around the tree, rock and soil of life (Arnold et al. 2009; Hodkinson & Lutzoni 2009; Cardinale et al. 2012; Bates et al. 2011; Hodkinson et al. 2012). Lichens have high tolerance to ecological stresses such as surviving in extremes of temperature, low nutrient or water availability and UV light intensities and therefore can be found in a wide range of habitats, where other organisms cannot survive (Beckett et al. 2008; Boddy et al. 2010).

The organic (process, biological process) process of a lichen structure begins when it comes into contact with a suitable photobiont. In most situations, the mycobiont dictates the appearance of a lichen structure, however there are few exceptions where the photobiont controls the environment. However, the photobiont has an impact on lichen ontogeny, and the lichen structure is not fully grown until each partner contributes (Budel & Scheidegger 2008). In their early stages, the partners create a tangled cell mass, eventually creating the lichen structure, which is largely made up of fungal mycelia, as well as the majority of biomass made up of mycobionts (Gilbert 2000). The lichen structure is generally a visibly overlaid structure that includes the cortex, affiliated protista layer, medulla, and rhizinae in certain circumstances (Jahns 1973). The dependent partners make contact at intervals in the protista layer, but the relationship between the protista and the fungal hyphae is different. The top cortex is made up of densely packed blood cells that cover the lichen structure's highest point. This lichen component layer can be hundreds of micrometres thick, yet it can also appear as a thin coating (Budel & Scheidegger 2008). Apart from offering protection against high radiation such as ultraviolet (UV) (Solhaug et al. 2003), actinic radiation, and actinic ray light, the upper cortex also has numerous roles such as altering energy budgets and defending against lichenivorous animals (Budel & Scheidegger 2008). The protista layer is found in the medulla's higher reaches and gets more light for photosynthesis by the photobiont (Gilbert 2000).

Though the fungus hyphae surround the protoctist cells, they're gently interlinking and leave enough area for optimum intensity. All an equivalent, the photobionts square measure rarely set deep within the lichen structure and there aren't any unattached protoctist within the lichen structure (Jahns 1973). The medullary layer makes up an oversized a neighborhood of the inner lichen structure in mostly leafy or fruticose lichens; this layer consists of long and loosely packed fungus hyphae, forming a soft layer with high internal airspace, facilitating gas exchange at intervals the lichen structure (Honegger 1993). In some lichens, notably fruticose lichens, supporting tissue is typically designed at intervals the medulla, making irregular forms that might be discovered in genera like genus Cladonia and Usnea (Budel & Scheidegger 2008).

The bacterial population found in lichens was first reported in the early twentieth century (Uphof 1925; Henkel and Yuzhakova 1936; Iskina 1938). Various bacterial species were identified to be linked with lichens in early investigations based on culture-dependent methods, including classes such as Gammaproteobacteria, Alphaproteobacteria and Firmicutes (Iskina 1938; Panosyan and Nikogosyan 1966; Henkel and Plotnikova 1973). The study of lichen associated bacterial community were analysed from three lichen samples such as Cladonia arbuscula, Lecanora polytropa and Umbilicaria cylindrical (Grube etal 2009)

Metagenomic analysis involves within the application of bioinformatics tools to review the genetic material from environmental, uncultured microorganisms. The Next generation sequencing of 16S rRNA lets the evaluation of bacterial diversity and distinguish thousands of organisms. There are three primary phases in metagenomic data processing: 1) assembly, 2) annotation, and 3) statistical analysis (Thomas et al. 2012).

Lichen is a symbiotic association of fungi and algae occurs and distributed in tree, rock and soil habitat. The symbiotic association of lichen establishment in tree, rock, soil is influenced by several bacteria, actinobacteria and cyanobacteria community. The type, kind of bacterial community in the lichen growing substrate like tree, rock, soil is unique and distinct in geographical locations and substrate chemical composition. Nevertheless, the report on corticolous lichen associated bacterial community analysis in Roccella montagnei in Tamil Nadu, India is very scanty. Keeping above points, and to accomplish the knowledge gap, Therefore, we are concerned to invent the bacterial profiling of Roccella montagnei corticolous lichens, which are grown under different and same tree substrate in distinct geographical locations by Next Generation Sequencing approach of Illumina Sequencing by 16S rRNA amplicon analysis.

2.1 Lichen sample collection

Corticolous lichen Roccella montagnei samples (Fig. 1) were collected from various locations in Kuppannasamy Temple, Pollachi (Latitude: 10.65° N, Longitude: 77.00° E) a Terrestrial habitat located in Coimbatore, Nithiravilai, Nagercoil (Latitude: 8.32° N, Longitude: 77.34° E) and Ramarpatham, Vedaranyam (Latitude: 10.28° N, Longitude: 79.82° E) both coastal habitats located in Kanyakumari and Nagapattinam districts of Tamil Nadu, India (Fig. 2). Lichen samples were excised using a collection kit and labelled KR, PR, and RR lichen sample from Kanyakumari, Pollachi, and Ramarpatham respectively.

2.2 Lichen sample preparation

The lichens were thoroughly washed with distilled water, surface disinfected by immersion in absolute ethanol for 5-10 seconds, followed by immersion in hydrogen peroxide (4%) for 90 seconds, and then rinsed three times for 60 seconds each with sterile demineralized water for surface sterilization.

2.3 Extraction of lichen metagenomic DNA

Following the modified procedures of Silva et al. (2014) and Hassan et al. (2018) total metagenomic DNA was extracted from the lichen samples. 2 g of the lichen samples were transferred to a sterile 15 mL centrifuge tube after crushing them into a paste form with a mortar and pestle using lysis buffer. Heat shock method was administrated by keeping the samples at 95℃ water bath for 1 hour and incubated in -80℃ freezer for 1 hour. Immediately thawing samples by placing it back in 95℃ water bath. Absolute ethanol was added double the sample volume, vortexed the samples, were centrifuged at 5000 rpm for 10 min. Then supernatant was discarded from each samples and add 2 ml of lysis buffer, 10µl lysozyme, 10µl proteinase K and 10µl RNase and vortexed for 15 minutes. The set up was incubated for 1.5 hours at 60°C with intermittent vortexing and a buffer-saturated phenol: chloroform: isoamyl alcohol (1:1:1) was transferred to the mixture at an equal volume. The aqueous layer was transferred to a sterile centrifuge tube after 20 minutes of centrifugation at 8000 rpm at 4°C, and an equal amount of chloroform: isoamyl alcohol (1:1) was added to the aqueous layer for the elimination of phenol traces and repeated twice. The supernatant was transferred to a sterile centrifuge tube with an equivalent amount of isopropanol and incubated on ice for 20 minutes after centrifugation. The supernatant was removed after centrifugation, and the pellet was rinsed with 100 µL of 80% ethanol and centrifuged for 5 minutes in 4°C at 8000 rpm. The supernatant was removed, and the DNA pellets were dried completely at room temperature. The pellet was resuspended in 50 µL of TE buffer and kept at -20°C until needed. A total of 9 lichen metagenomics DNA (3 samples from each location) was extracted from Roccella montagnei and mixed together and used for analysis.

2.4 16S rRNA V3/V4 gene amplification by NGS:

The following procedures were followed for the creation of next-generation sequencing libraries and Illumina MiSeq sequencing. A Qubit 2.0 Fluorometer was used to measure DNA samples (Invitrogen, Carlsbad, CA, USA). A MetaVxTM Library Preparation kit was used to produce amplicons from 30–50 ng of DNA (GENEWIZ, Inc., South Plainfield, NJ, USA). The hypervariable V3 and V4 regions of bacterial 16S rRNA were chosen for amplicon sequencing, followed by taxonomic analysis. GENEWIZ created a set of proprietary primers aiming at the relatively conserved areas of the V3 and V4 hypervariable regions of the 16S rRNA of bacteria.

2.5 Sequencing and Library preparation

The Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA) was used to confirm library quality control, and the Qubit 2.0 Fluorometer was used to quantify it. Following the manufacturer's recommendations, DNA libraries were multiplexed and put into an Illumina MiSeq device (Illumina, San Diego, CA, USA). The MiSeq control software package (MCS) integrated in the MiSeq instrument was used to sequence using a 2300 paired end (PE) configuration; image analysis and base calling were accomplished using the MiSeq control software package (MCS). The metagenome sequence of Lichen associated bacterial community in Roccella montagnei was submitted in NCBI with the accession number of bioproject PRJNA741221.

3.1 Sampling and Process of Lichen Roccella montagnei

The Roccella montagnei was collected from two different habitats of sampling locations like coastal and terrestrial regions. In the present study, we have examined the bacterial community profiling of the Roccella montagnei collected from the different habitat, corticolous lichen grown under different and same tree substrate in distinct geographical locations (Table S1). A total of 3 Roccella montagnei lichen samples were used for metagenomic DNA extraction. The lichen samples were abbreviated as PR lichen from terrestrial habitat and RR, KR lichen from coastal habitat.

3.2 Lichen Microbiome profiling

QIIME pipeline (1.9.1) was used to analyses the bacterial diversity in the lichen genera of Roccella montagnei obtained from the coastal and terrestrial habitats. In the 16S rRNA dataset 2.1 lakhs sequences number were obtained having a combined length of 390 - 470 base pairs (bp) (Table S2).

3.3 Common and unique operational taxonomic units (OTUs)

16S rRNA V3 – V4 gene profiles were generated from the lichen samples to compare and evaluate the diversity and abundance of bacterial association in the lichen Roccella montagnei. The circles with different colors in the Venn diagram (Fig. 3), each represents a different sample of corticolous lichen grown under same tree substrate in distinct geographical locations and numbers within them represents the number OTUs which are either common or unique to the samples, the number inside the overlapping part denotes the common OTUs whereas the in the non-overlapping ones are unique to the samples.

A sum total of 424 OTUs were observed from three samples, out of which 205 OTUs were in common with all three lichen samples, 35 OTUs are unique to KR lichen sample from coastal habitat, 34 for PR sample from terrestrial habitat and about 42 for RR lichen from coastal habitat. It is also observed that the highest OTUs were shared between KR and RR lichen sample which are both from coastal habitats; on the other hand, PR and RR lichen share the lowest of 20 OTUs.

The heatmap analysis illustrates the abundance of bacterial community information of selected OTU as well as the similarity and difference across OTUs and lichen samples by similarity clustering. The figure showed the top 30 OTUs with the highest abundance of bacteriome.

3.4 Taxonomic analysis of bacterial community abundance on Roccella montagnei

The bacterial community abundance was analyzed in Roccella montagnei for all taxonomic levels and the results were observed as follows, their respective stacked bar plots.

The most abundant phyla in the 3 lichen Roccella montagnei samples, which were in coastal and terrestrial habitat were (i) Proteobacteria with 95.78% present in Nagercoil lichen sample, 96.54% in Pollachi lichen sample, and a 96.16% in lichen sample collected from Vedaranyam; (ii) Bacteroidetes having a presence of 1.36% in Nagercoil lichen sample, 1.87% in Pollachi lichen sample and 1.77% in Vedaranyam sample. The other phyla are in decreasing order as follow: Actinobacteria, Acidobacteria, Firmicutes, Unclassified phyla, Cyanobacteria, Fusobacteria, Gemmatimonadetes, Tenericutes, WPS-2, Chloroflexi, Verrucomicrobia (Fig. 4 (a)). There are 4 common and significantly abundant classes in all 3 corticoles lichen Roccella montagnei grown under coastal and terrestrial habitat in the tree bark of Cocos nucifera and Maba buxifolia. (i) Gammaproteobacteria: 65.68% in Nagercoil lichen in Cocos nucifera (Coastal), 58.92% in Pollachi lichen from Cocos nucifera (Terrestrial), 56.85% in Vedaranyam lichen from Maba buxifolia (Coastal). (ii) Alphaproteobacteria: 25.71% in Nagercoil sample, 32.43% in Pollachi sample, 34.73% in Vedaranyam sample; (iii) Betaproteobacteria: 4.36% in Nagercoil sample, 4.96% in Pollachi sample, 4.53% in Vedaranyam sample; (iv) Flavobacteriia: 1.14% in Nagercoil sample, 1.81% in Pollachi sample, 1.55% in Vedaranyam sample. The other classes are in decreasing order as follows: Actinobacteria, Acidobacteriia, Unclassified classes, Bacilli, Saprospirae, Deltaproteobacteria, Chloroplast, Clostridia, Oscillatoriophycideae, Fusobacteriia, Gemmatimonadetes, Mollicutes, Solibacteres, Synechococcophycideae, Rubrobacteria, Anaerolineae, Sphingobacteriia, Acidimicrobiia, Spartobacteria (Fig. 4 (b)). There are about 9 significantly abundant and orders common to all the 3 corticolous lichen Roccella montagnei grown under coastal and terrestrial habitat in the tree bark of Cocos nucifera and Maba buxifolia. (i) Pseudomonadales: 42.04% in Nagercoil sample, 34.96% in Pollachi sample, 33.61% in Vedaranyam sample, (ii) Rhizobiales: 15.99% in Nagercoil sample, 27.89% in Pollachi sample, 29.1% in Vedaranyam sample, (iii) Xanthomonadales: 10.98% in Nagercoil sample, 17.48% in Pollachi sample, 8.82% in Vedaranyam sample, (iv) Enterobacteriales: 12.46% in Nagercoil sample, 6.27% in Pollachi sample, 12.85% in Vedaranyam sample, (v) Burkholderiales: 4.3% in Nagercoil sample, 4.9% in Pollachi sample, 4.48% in Vedaranyam sample, (vi) Rhodospirillales: 5.73% in Nagercoil sample, 0.74% in Pollachi sample, 1.77% in Vedaranyam sample, (vii) Caulobacterales: 2.4% in Nagercoil sample, 2.3% in Pollachi sample, 1.74% in Vedaranyam sample, (viii) Sphingomonadales: 1.31% in Nagercoil sample, 1.27% in Pollachi sample, 1.98% in Vedaranyam sample, (ix) Flavobacteriales: 1.14% in Nagercoil sample, 1.81% in Pollachi sample, 1.55% in Vedaranyam sample (Fig. 5 (c)). In family level there are about 11 common important families abundant in the 3 lichen corticolous lichen Roccella montagnei grown under coastal and terrestrial habitat in the tree bark of Cocos nucifera and Maba buxifolia. (i) Pseudomonadaceae: 32.9% in Nagercoil sample, 24.75% in Pollachi sample, 26.21% in Vedaranyam sample, (ii) Brucellaceae: 7.92% in Nagercoil sample, 23.41% in Pollachi sample, 18.14% in Vedaranyam sample, (iii) Xanthomonadaceae: 10.89% in Nagercoil sample, 17.48% in Pollachi sample, 8.67% in Vedaranyam sample, (iv) Enterobacteriaceae: 12.46% in Nagercoil sample, 6.27% in Pollachi sample, 12.85% in Vedaranyam sample, (v) Moraxellaceae: 9.14% in Nagercoil sample, 10.21% in Pollachi sample, 7.4% in Vedaranyam sample, (vi) Alcaligenaceae: 3.41% in Nagercoil sample, 4.4% in Pollachi sample, 3.67% in Vedaranyam sample, (vii) Methylocystaceae: 1.99% in Nagercoil sample, 2.28% in Pollachi sample, 5.98% in Vedaranyam sample, (viii) Rhizobiaceae: 4.03% in Nagercoil sample, 0.91% in Pollachi sample, 3.27% in Vedaranyam sample, (ix) Acetobacteraceae: 5.72% in Nagercoil sample, 0.74% in Pollachi sample, 1.75% in Vedaranyam sample, (x) Caulobacteraceae: 2.4% in Nagercoil sample, 2.3% in Pollachi sample, 1.74% in Vedaranyam sample, (xi) Sphingomonadaceae: 1.31% in Nagercoil sample, 1.27% in Pollachi sample, 1.97% in Vedaranyam sample (Fig. 5 (d)). Eight commonly noteworthy genera are present in all 3 corticolous lichen Roccella montagnei grown under coastal and terrestrial habitat in the tree bark of Cocos nucifera and Maba buxifolia Along the genera level, one genus is signified as unclassified. (i) Unclassified: 32.34% in Nagercoil sample, 13.35% in Pollachi sample, 26.79% in Vedaranyam sample, (ii) Pseudomonas: 23.75% in Nagercoil sample, 24.42% in Pollachi sample, 23.06% in Vedaranyam sample,(iii) Ochrobactrum: 7.84% in Nagercoil sample, 21.27% in Pollachi sample, 18.09% in Vedaranyam sample, (iv) Stenotrophomonas: 10.15% in Nagercoil sample, 17.47% in Pollachi sample, 8.05% in Vedaranyam sample, (v) Acinetobacter: 9.12% in Nagercoil sample, 10.19% in Pollachi sample, 7.4% in Vedaranyam sample, (vi) Achromobacter: 3.39% in Nagercoil sample, 4.37% in Pollachi sample, 3.63% in Vedaranyam sample, (vii) Agrobacterium: 3.83% in Nagercoil sample, 0.79% in Pollachi sample, 3.14% in Vedaranyam sample, (viii) Brevundimonas: 0.7% in Nagercoil sample, 2.05% in Pollachi sample, 0.52% in Vedaranyam sample (Fig. 6 (e)).

3.5 Alpha Diversity analysis of the lichen microbiome communities

Alpha diversity is mainly used to reflect the diversity of each lichen sample, which estimates the number of species in the microbial community as well as the abundance and diversity of species in environmental communities through a series of statistical indices. In which are summarized in terrestrial and coastal habitat lichens in distinct geographic locations of tree as growing substrate of Cocos nucifera and Maba buxifolia in Nagercoil, Vedaranyam, Pollachi of Tamil Nadu, India. The indices for community richness calculation include: Chao and ACE index. The Shannon and Simpson index for community diversity calculations. The Goods coverage indices to measure sequencing depth (Coverage). Nagercoil lichen had the highest Shannon index (4.75) which indicates abundance of the lichen associated bacterial community present in a significantly maximum in the trees of Cocos nucifera in coastal habitat in comparison with other two lichen samples i.e., Pollachi sample (4.12) and Vedaranyam sample (4.515).

The Simpson’s index shows that quantify the diversity of species richness and evenness is almost the same in Nagercoil lichen from Cocos nucifera in terrestrial locations (0.923) and Vedaranyam lichen from Maba buxifolia in Coastal locations (0.909) and both the samples were higher than the Pollachi Cocos nucifera in terrestrial location (0.884).

The Fig. 7 (a) shows the rarefaction curve i.e., it tries to determine all of the bacterial diversity in a sample at a particular sequencing depth by sub sampling the sequence to a minimum and plotting it against the observed OTUs which enables the comparison among the different lichen samples in of different and same tree substrate in distinct geographical locations. The curves represent the bacterial diversity as a function of sequencing depth. The X axis is the number of valid sequences extracted, and the Y axis is the number of OTUs (Observed OTUs). Each sample is represented by one curve with a unique color. The number of OTUs increases with the increase of extracted sequence count until reaching a plateau, which indicates the number of detected OTUs will not increase with the number of extracted sequences and reflects the reasonable sequence depth.

Rank-abundance curve is used to analyze bacterial diversity. Rank-abundance curve reflects both bacterial species abundance and species uniformity in lichen. The abundance of species is reflected by the length of the curve on the X axis. The more extended on the X axis, the more abundant the species is. Species uniformity is reflected by the shape of the curve. The smoother the curves, the higher the species uniformity in lichen Roccella montagnei.

Each curve in the figure above corresponds to an individual lichen Roccella montagnei samples from three locations. The X axis is the relative abundance of the OTU in descending order. The Y axis is the relative abundance of the OTU. '100' on the X axis indicates the OTU in the sample is ranked as the 100th abundant in descending order, and the corresponding value on the Y axis is the percentage of the sequence count in the OTU (the number of sequences of the OTU divided by the total number of sequences) (Fig. 7 (b)).

3.6 Beta diversity of the lichen microbiome communities:

Beta diversity analysis is used to determine how different the lichen bacterial community in coastal samples is from that in terrestrial habitats, as well as to disclose elements of microbial ecology that aren't apparent simply looking at the sample composition. Beta diversity analysis is frequently divided into two categories. They can be quantitative (e.g., Bray-Curtis or weighted UniFrac) or qualitative (e.g., Bray-Curtis or weighted UniFrac) (Goodrich et al. 2014). The NMDS, PCA, and PCoA figures were all created using beta diversity distance matrices.

3.6.1 PCA and PCoA plot analysis:

PCA analysis (Principal Component Analysis): The difference and distance between samples can be reflected by the analysis of the gene functional distribution of different samples of Roccella montagnei. The first, second, and third main components are represented by PC1, PC2, and PC3, respectively (Fig. 8 (a)). The percentage after the main component, which reflects how much information the primary component can extract from the original data, represents the contribution rate of this component to sample difference. The distance between coastal and terrestrial lichen demonstrates the similarity of the distribution of functional categories in the sample.

The PCoA (Principal Co-ordinates Analysis) is a method for analysing and visualising data differences and similarities. Illustrations from the same group are identified by the similar colour and form. The PCoA plot for the first and second principal coordinates is shown in the PC1 vs PC2 image (Fig. 8 (b)), with the first and second axes representing the first and second principal coordinates, respectively. The contribution of the related coordinate to the sample variations is shown by the percentage number in the axis label, which illustrates how much information is recovered from the original data. The distance between the sample points demonstrates the similarity of the bacterial population in the lichen.

According to the principle component analysis (PCA) and Principal Co-ordinates Analysis (PCoA), the lichen samples from Nagercoil, Pollachi, and Vedaranyam were differed significantly. Mostly, the Nagercoil lichen from Cocos nucifera in coastal habitat and Vedaranyam lichen from Maba buxifolia in coastal locations, indicated similarity than the other sample comparisons.

3.6.2 UPGMA Tree

Cluster analysis uses evolutionary data from sample sequencing to determine if samples in a well-known environment are significantly different from an evolutionary branch in bacterial communities. In UPGMA tree each branch represents a sample and difference of colours will represent different groups present. The results (Fig. 9), infers the all the 3-lichen bacterial community belongs to the same group and the evolutionary microbial lineage between Nagercoil lichen from Cocos nucifera from coastal habitat and Vedaranyam lichen from Maba buxifolia from Coastal habitat are similar and the highly different from the Pollachi lichen from Cocos nucifera from terrestrial habitat.

The composition of bacterial microbiota has been well studied in lichens, nevertheless there are only limited reports in respect to comparing the bacteriome of the same lichen species from different types of locations in Tamil Nadu, India and Rocella montagnei is a species whose bacteriome had not been explored previously. In this study we have collected Rocella montagnei from three different locations that were near coastal and terrestrial habitats and metagenomic analysis and their bacteriobiome records through NGS- Illumina sequencing platform. The predominance of proteobacteria, bacteroidetes, actinobacteria, acidobacteria and verrucomicrobia in phylum level revealed and consistent with previous reports reports (Hodkinson et al. 2012; Leiva et al. 2021; Printzen et al. 2012; Sierra et al. 2020). In the lichen genera such as Cora, Hypotrachyna, Cladonia, Usnea, Sticta, Stereocaulon and Peltigera the presence of these members of this phylum are considered to play important roles such as symbioses by supplying nutrients, iron and phosphate mobilization, and fixing nitrogen (Grube et al. 2015; Sigurbjornsdottir and Vilhelmsson 2016). The dominance of Gammaproteobacteria over Alphaproteobacteria at class level was found in our results and this finding is supported by similar discovery reported by Sierra et al. (2020) while most of the reports majorly showed the dominance of alphaproteobacterial which is ubiquitous in all the lichen of different habitats (Almendras et al. 2018; Aschenbrenner et al. 2014; Cernava and Müller et al. 2015; Hodkinson et al. 2012; Leiva et al. 2021; Printzen et al. 2012). The predominance of Pseudomonadales of Gammaproteobacteria indicates a different result from the previous studies that revealed typically the dominance of Rhizobiales, which is only the second most dominant order in our study, other orders such as Rhodospirillales, Sphigomonadales were observed in reports of Almendras et al. (2018); Aschenbrenner et al. (2014); Cernava et al. (2017); Hodkinson et al. (2012); Leiva et al. (2021); Printzen et al. (2012); Sierra et al. (2020) these orders help in fixing nitrogen in the lichens. In the family level Psedumondaceae are highly dominant followed by Brucellaceae followed by Aschenbrenner et al. (2014); Hodkinson et al. (2012); Printzen et al. (2012); Sierra et al. (2020). Xanthomonadaceae and others like Sphingomonadaceae and Acetobacteraceae that are mostly similar to the reports. The genus level bacteriome exhibits the major dominance of unclassified genera on the other hand the classified genera were dominated by Pseudomonas (Aschenbrenner et al. 2017; Cernava and Aschenbrenner et al. 2015; Cernava and Müller et al. 2015; Lee et al. 2014), Ochrobactrum which is surprisingly never before found in the lichen bacteriome and is followed by Stenotrophomonas (Cernava and Aschenbrenner, et al. 2015; Cernava and Muller et al. 2015). The occurrence and presence of Stenotrophomonas indicates that the lichen can survives drought and high levels of salinity by producing protecting compounds such as spermidine and different osmolytes, which is relatable in our study as two of the samples where from coastal habitat, Sphingomonas (Aschenbrenner et al. 2017; Graham et al. 2018; Lee et al. 2014) that helps the lichen grow by producing growth promoting hormones, Acinetobacter (Park et al. 2016) and others. The alpha and Beta diversity analysis are done to identify the similarity of the microbial population of the lichen form the different habitats and all the analysis has led to the conclusion of similarity of coastal lichens and then distinctively different from the terrestrial lichen as proven by the UPGMA tree at an evolutionary lineage level in the microbial communities.

The metagenomic analysis of lichen Roccella montagnei data that has been generated, the findings understand that a total 424 OTUs were obtained from three lichens. Among the total OTU Vedaranyam lichen coastal location sample showed the greatest number of unique OTUs (42) and lowest was Pollachi lichen from Cocos nucifera in terrestrial habitat. However, the majority OTUs (57) were shared between Nagercoil lichen from Cocos nucifera (Coastal); Vedaranyam lichen from Maba buxifolia in Coastal locations than by Pollachi sample the others. The taxonomic analysis of microbial community abundance on Roccella montagnei showed the abundance of Proteobacteria, Bacteroidetes and Actinobacterias in Phylum level; Gammaproteobacteria, Alphaproteobacteria and Betaproteobacteria in Class level; Pseudomonadales, Rhizobiales and Xanthomonadales in Family level; Unclassified generas, Pseudomonas, Ochrobactrum and Stenotrophomonas in the Genus level; and finally Unclassified species, diminuta, marcescens and meningoseptica in Species level of taxonomy. All the Alpha and Beta diversity analysis using various tools have showed variations present in all sample but it has led to the conclusion of similarity of Nagercoil and Vedharanyam coastal lichen samples, the UPGMA tree confirms their evolutionary lineage. The result reveals that, the bacterial community variation in lichen is due to the geographic origin, the biotic and abiotic factors surrounding the environment, the substrate in which they grow, the microhabitat conditions around the lichen sample.

The data were analysed in terms of the lichen-associated bacterial community's stability under shifting abiotic settings. They add to our knowledge of a fascinating lifeform by demonstrating how the bacterial component of the holobiont maintains its integrity in the face of adverse environmental conditions. Finally, this study data will be a contribution of bacterial community profile that were uncultivable and their diversity, similarity and dissimilarity has been revealed to this society.

Acknowledgements

The first author thanks the Research Foundation of Bharathidasan University for the fellowship (URF) (021505/URF/DIR-RES/2019 Date: 22.10.2019). DST-FIST (Department of Science and Technology-Fund for improvement of S & T infrastructure), New Delhi, is express gratitude for the instrumentation facility (DST Sanction Order No.: SR/FIST/LSI-013/ 2012/Dt.13.08.2012). We acknowledge the financial assistance of the Department of Science and Technology, Government of India, under DST-Promotion of University Research and Scientific Excellence (PURSE) scheme - Phase II, Rashtriya Uchchatar Shiksha Abhiyan (RUSA)-2. O.

Almendras K, Garcia J, Caru M, Orlando J (2018) Nitrogen-fixing bacteria associated with Peltigera Cyanolichens and Cladonia Chlorolichens. Molecules 23(12):3077
Arnold AE, Miadlikowska J, Higgins KL, Sarvate SD, Gugger P, Way A, Hofstetter V, Kauff F, Lutzoni F (2009) A phylogenetic estimation of trophic transition networks for ascomycetous fungi: are lichens cradles of symbiotrophic fungal diversification? Syst Biol 58:283–297
Aschenbrenner IA, Cardinale M, Berg G, Grube M (2014) Microbial cargo: Do bacteria on symbiotic propagules reinforce the microbiome of lichens? Environ Microbiol 16(12):3743–3752
Aschenbrenner IA, Cernava T, Erlacher A, Berg G, Grube M (2017) Differential sharing and distinct co-occurrence networks among spatially close bacterial microbiota of bark, mosses and lichens. Mol Ecol 26(10):2826–2838
Bates ST, Cropsey GWG, Caporaso JG, Knight R, Fierer N (2011) Bacterial communities associated with the lichen symbiosis. Appl Environ Microbiol 77:1309–1314
Beckett RP, Kranner I, Minibayeva FV (2008) Stress physiology and the symbiosis. In Lichen Biology, Edited by T. H. Nash. Cambridge: Cambridge University Press 2nd edition, pp. 134– 151
Boddy L, Dyer PS, Helfer S (2010) Plant pests and perfect partners. In from Another Kingdom, Edited by L. Boddy & M. Coleman. Edinbugh: Royal Botanic Gardens pp. 52–65
Budel B, Scheidegger DC (2008) Thallus morphology and anatomy. In: Nash THIII (ed) Lichen biology, 2nd edn. Cambridge University Press, Cambridge, pp 40–68
Cardinale M, Steinova J, Rabensteiner J, Berg G, Grube M (2012) Age, sun and substrate: triggers of bacterial communities in lichens. Environmental Microbiology Reports 4:23–28
Cernava T, Aschenbrenner IA, Grube M, Liebminger S, Berg G (2015) A novel assay for the detection of bioactive volatiles evaluated by screening of lichen-associated bacteria. In Frontiers in Microbiology (Vol. 6, Issue MAY)
Cernava T, Muller H, Aschenbrenner IA, Grube M, Berg G (2015) Analyzing the antagonistic potential of the lichen microbiome against pathogens by bridging metagenomic with culture studies. In Frontiers in Microbiology (Vol. 6, Issue JUN)
Gilbert O (2000) Lichens. Harper Collins Publishers, New York (NY)
Goodrich JK, Di Rienzi SC, Poole AC, Koren O, Walters WA, Caporaso JG, Ley RE (2014) Conducting a Microbiome Study. Cell 158(2):250–262
Graham LE, Trest MT, Will-Wolf S, Miicke NS, Atonio LM, Piotrowski MJ, Knack JJ (2018) Microscopic and metagenomic analyses of Peltigera ponojensis (Peltigerales, ascomycota). Int J Plant Sci 179(3):241–255
Grube M, Cardinale M, de Castro JV, Muller H, Berg G (2009) Species-specific structural and functional diversity of bacterial communities in lichen symbioses. ISME J 3(9):1105–1115
Hassan M, Essam T, Megahed S (2018) Illumina sequencing and assessment of new cost-efficient protocol for metagenomic-DNA extraction from environmental water samples. Braz JMicrobiol 49:S1–S8
Henkel PA, Plotnikova TT (1973) Nitrogen-fixing bacteria in lichens. Izv Akad Nauk Ser Biol 1973: 807–813
Henkel PA, Yuzhakova LA (1936) Nitrogen-fixing bacteria in lichens. Izv Biol Inst Permsk Gos Univ 10:9–10
Hodkinson BP, Lutzoni F (2009) A microbiotic survey of lichen-associated bacteria reveals a new lineage from the Rhizobiales. Symbiosis 49:163–180
Hodkinson BP, Gottel NR, Schadt CW, Lutzoni F (2012) Photoautotrophic symbiont and geography are major factors affecting highly structured and diverse bacterial communities in the lichen microbiome. Environ Microbiol 14:147–161
Honegger R (1993) Tansley Review No. 60. Developmental Biology of Lichens. New Phytol 125(4):659–677
Iskina RY (1938) On nitrogen fixing bacteria in lichens. Isv Biol Inst Permsk 11:133–139
sol Anatomy, morphology and development. London, England: Academic Press, In Ahmadjian V (ed), Lichens. pp. 3–57
Lee YM, Kim EH, Lee HK, Hong SG (2014) Biodiversity and physiological characteristics of Antarctic and Arctic lichens-associated bacteria. World Journal of Microbiology Biotechnology 30(10):2711–2721
Leiva D, Fernandez-Mendoza F, Acevedo J, Caru M, Grube M, Orlando J (2021)
The Bacterial Community of the Foliose Macro-lichen Peltigera frigida Is More than a Mere Extension of the Microbiota of the Subjacent Substrate. Microbial Ecology, 81(4), 965–976
Panosyan AK, Nikogosyan VG (1966) The presence of Azotobacter in lichens. Akad Nauk Armian SSR Biol Zhurn Armen 19:3–11
Park CH, Kim KM, Kim OS, Jeong G, Hong SG (2016) Bacterial communities in Antarctic lichens. Antarct Sci 28(6):455–461
Printzen C, Fernandez-Mendoza F, Muggia L, Berg G, Grube M (2012) Alphaproteobacterial communities in geographically distant populations of the lichen Cetraria aculeata. FEMS Microbiol Ecol 82(2):316–325
Sierra MA, Danko DC, Sandoval TA, Pishchany G, Moncada B, Kolter R, Mason CE, Zambrano MM (2020) The Microbiomes of Seven Lichen Genera Reveal Host Specificity, a Reduced Core Community and Potential as Source of Antimicrobials. Front Microbiol 11(March):1–12,398
Sigurbjornsdottir MA, Vilhelmsson O (2016) Selective isolation of potentially phosphate-mobilizing, biosurfactant-producing and biodegradative bacteria associated with a sub-Arctic, terricolous lichen, Peltigera membranacea. FEMS Microbiol Ecol 92:fiw090
Silva MAL, d Medeiros Z, Soares CRP, da Silva ED, Miranda-Filho DB, deMelo FL (2014) A comparison of four DNA extraction protocols for the analysis of urine from patients with visceral leishmaniasis. Rev Soc Bras Med Trop 47(2):193–197
Solhaug KA, Gauslaa Y, Nybakken L, Bilger W (2003) UV-induction of sun-screening pigments in lichens. New Phytol 158(1):91–100
Thomas T, Gilbert J, Meyer F (2012) Metagenomics-a guide from sampling to data analysis. Microb Inform Exp 2(1):3–8
Uphof JT (1925) Purple bacteria as symbionts of a lichen. Science 61(1568):67–67

Supplementaryfiles.docx

PDF

Journal Publication

published 22 Dec, 2021

Read the published version in Archives of Microbiology →

First submitted to journal
04 Oct, 2021
Editorial decision: Minor revisions
19 Sep, 2021

You are reading this latest preprint version

Metagenomic Analysis of Lichen Associated Bacterial Community Profiling In Roccella Montagnei

Status:

Journal Publication

Version 1

Abstract

Figures

1 Introduction

2. Materials And Methods

2.1 Lichen sample collection

2.2 Lichen sample preparation

2.3 Extraction of lichen metagenomic DNA

2.4 16S rRNA V3/V4 gene amplification by NGS:

2.5 Sequencing and Library preparation

3 Results

3.1 Sampling and Process of Lichen Roccella montagnei

3.2 Lichen Microbiome profiling

3.3 Common and unique operational taxonomic units (OTUs)

3.4 Taxonomic analysis of bacterial community abundance on Roccella montagnei

3.5 Alpha Diversity analysis of the lichen microbiome communities

3.6 Beta diversity of the lichen microbiome communities:

3.6.1 PCA and PCoA plot analysis:

3.6.2 UPGMA Tree

4 Discussion

5 Conclusions

Declarations

References

Supplementary Files

Status:

Journal Publication

Version 1