Next Article in Journal / Special Issue
Isolation, Characterization, and Antibacterial Activity of Hard-to-Culture Actinobacteria from Cave Moonmilk Deposits
Previous Article in Journal
Geographic Variation in Antibiotic Consumption—Is It Due to Doctors’ Prescribing or Patients’ Consulting?
Previous Article in Special Issue
Novel Aspects of Polynucleotide Phosphorylase Function in Streptomyces
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Antibiotics
Volume 7
Issue 2
10.3390/antibiotics7020027
antibiotics-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High-Throughput Sequencing Analysis of the Actinobacterial Spatial Diversity in Moonmilk Deposits

by
Marta Maciejewska
1,
Magdalena Całusińska
2,
Luc Cornet
3,
Delphine Adam
1,
Igor S. Pessi
1,
Sandrine Malchair
4,
Philippe Delfosse
2,
Denis Baurain
3,
Hazel A. Barton
5,
Monique Carnol
4 and
Sébastien Rigali
1,*
1
InBioS—Centre for Protein Engineering, Institut de Chimie B6a, University of Liège, B-4000 Liège, Belgium
2
Environmental Research and Innovation Department, Luxembourg Institute of Science and Technology, Belvaux, Luxembourg
3
InBioS—PhytoSYSTEMS, Eukaryotic Phylogenomics, University of Liège, B-4000 Liège, Belgium
4
InBioS—Plant and Microbial Ecology, Botany B22, University of Liège, B-4000 Liège, Belgium
5
Department of Biology, University of Akron, Akron, OH 44325, USA
*
Author to whom correspondence should be addressed.
Antibiotics 2018, 7(2), 27; https://doi.org/10.3390/antibiotics7020027
Submission received: 12 February 2018 / Revised: 15 March 2018 / Accepted: 16 March 2018 / Published: 21 March 2018
(This article belongs to the Special Issue Actinomycetes: The Antibiotics Producers)
Browse Figures
Versions Notes

Abstract

:
Moonmilk are cave carbonate deposits that host a rich microbiome, including antibiotic-producing Actinobacteria, making these speleothems appealing for bioprospecting. Here, we investigated the taxonomic profile of the actinobacterial community of three moonmilk deposits of the cave “Grotte des Collemboles” via high-throughput sequencing of 16S rRNA amplicons. Actinobacteria was the most common phylum after Proteobacteria, ranging from 9% to 23% of the total bacterial population. Next to actinobacterial operational taxonomic units (OTUs) attributed to uncultured organisms at the genus level (~44%), we identified 47 actinobacterial genera with Rhodoccocus (4 OTUs, 17%) and Pseudonocardia (9 OTUs, ~16%) as the most abundant in terms of the absolute number of sequences. Streptomycetes presented the highest diversity (19 OTUs, 3%), with most of the OTUs unlinked to the culturable Streptomyces strains that were previously isolated from the same deposits. Furthermore, 43% of the OTUs were shared between the three studied collection points, while 34% were exclusive to one deposit, indicating that distinct speleothems host their own population, despite their nearby localization. This important spatial diversity suggests that prospecting within different moonmilk deposits should result in the isolation of unique and novel Actinobacteria. These speleothems also host a wide range of non-streptomycetes antibiotic-producing genera, and should therefore be subjected to methodologies for isolating rare Actinobacteria.

1. Introduction

Molecular approaches evaluating microbial communities in caves have revealed a level of diversity greater than initially expected [1]. Microorganisms have been found to inhabit virtually all subterranean niches, including cave walls, ceilings, speleothems, soils, sediments, pools, and aquifers [2]. Cave bacteria often represent novel taxonomic groups [3,4,5,6,7], which are frequently more closely related to other cave-derived bacterial lineages than to the microbiota of other environments [8,9,10].
Among cave speleothems, moonmilk draws a particular scientific attention due to its distinctive crystal morphology. The origins of various moonmilk crystalline habits, including monocrystalline rods, polycrystalline chains, and nanofibers, are tentatively attributed to the moonmilk indigenous microbial population [11]. Among a moonmilk microbiome comprising Archaea, Fungi, and Bacteria [9,10,12,13,14,15,16,17,18,19], the indigenous filamentous Fungi [20] and Actinobacteria [11,21] are believed to mediate moonmilk genesis with cell surfaces promoting CaCO3 deposition [11,20,21]. Actinobacteria were additionally reported to be metabolically capable of inducing favorable conditions for CaCO3 precipitation, or even directly precipitating carbonate minerals [12,21]. Members of the phylum Actinobacteria are routinely found in this speleothem [9,10,12,13,14,18,19], as well as in the other subterranean deposits within limestone caves [3,8,22,23], volcanic caves [24,25,26], and ice caves [27]. The broad distribution of Actinobacteria in the subsurface systems stimulates investigation in order to understand the factors driving their existence in mainly inorganic and highly oligotrophic environments, and the processes that enable them to mediate speleogenesis. The successful adaptation of Actinobacteria to a wide range of environments could probably be a consequence of their broad-spectrum metabolism, which includes prolific secreted hydrolytic systems that are capable of generating nutrient sources from various substrates, along with their extraordinary faculty to produce specialized metabolites (metal chelators, antimicrobials, hormones, etc.) [28].
As recently reported, moonmilk Actinobacteria represent novel microorganisms, which is a discovery that opens great avenues for the bioprospecting of novel drugs [6,10,18]. Rooney et al. (2010) [13] showed that spatially separated moonmilk speleothems in Ballynamintra Cave are inhabited by taxonomically distinct fungal and bacterial communities. Instead, in our attempt to isolate moonmilk-dwelling Actinobacteria for assessing their potential for participating in the genesis of these speleothems [21] and producing antimicrobial compounds [10], we only recovered members of the genus Streptomyces. Such a dominance of streptomycetes was rather unexpected, according to other moonmilk microbial diversity studies performed through culture-dependent [10,12,13,18] and culture-independent approaches using clone libraries [9], denaturing gradient gel electrophoresis (DGGE) fingerprinting [14,16,17], automated ribosomal intergenic spacer analysis (ARISA) [13], and, more recently, high-throughput sequencing (HTS) [19]. The actinobacterial genera identified in those studies included Rhodococcus, Pseudonocardia, Propionibacterium, Nocardia, Amycolatopsis, Saccharothrix, Geodermatophilus, Mycobacterium, Aeromicrobium, Kribella, Nocardioides, Actinomycetospora, Nonomuraea, Euzebya, Rubrobacter, and Arthrobacter, in addition to Streptomyces. Nonetheless, the diversity of the moonmilk actinobacterial microbiome still remains largely unknown and, beyond evaluating “what and how much have we missed in our culture-dependent bioprospecting approach” [10], a major important question that arises is: to what extent are moonmilk-dwelling Actinobacteria different between the moonmilk deposits within a single cave, or in different caves?
In this work, we carried out a comparative (HTS) of 16S small subunit (SSU) rRNA gene from DNA extracted from spatially separated moonmilk deposits within the same cave, “Grotte des Collemboles” (Springtails’ Cave) in Comblain-au-Pont, Belgium (Figure S1), in order to draw a detailed taxonomic picture of the intra-phylum diversity. Identifying the presence of rare Actinobacteria and unveiling to which degree they exhibit a spatial variability would help determining whether it is worth prospecting from different moonmilk deposits to isolate unique and novel natural compound producers.

2. Results

2.1. Actinobacterial Abundance within the Whole Moonmilk Bacterial Microbiome

Libraries spanning the V4–V6 variable regions of the 16S rRNA gene using universal bacterial primers were used to assess the proportion of Actinobacteria in comparison to the whole bacterial community of three moonmilk deposits of the cave “Grotte des Collemboles” (Table S1a). The observed bacterial communities differed in species richness, evenness, and diversity between the three sampling points (Table 1). Phylotype richness (total number of operational taxonomic units (OTUs) per site) was the highest in COL4 (1863 OTUs), followed by COL1 and COL3, with 1332 and 1161 OTUs, respectively (Table 1, Figure 1a). Across the three sampling points, we found a total of 2301 different OTUs, amongst which 710 (31%) were common to all of the deposits (Figure 1a). Interestingly, pairwise comparison revealed highly similar percentages (~31.7 ± 0.53%) of shared bacterial OTUs between moonmilk deposits (Table 2, Figure 1a). A total of 956 OTUs (42%) were found to be exclusive to one sampling site, with COL4 having the highest number of unique bacterial phylotypes (584 OTUs), along with the most diverse bacterial population, as reflected by the highest diversity indices (Table 1, Figure 1a).
Bacterial OTUs were grouped into 21 phyla and 18 candidate phyla (Table S2, Figure 2). Actinobacteria represented 9%, 23%, and 10% of the total bacterial population in COL1, COL3, and COL4, respectively (Table S2, Figure 2a). In terms of abundance, they were the most common phylum after Proteobacteria, which accounted for 52%, 34%, and 30% of the total community in COL1, COL3, and COL4, respectively (Table S2, Figure 2a). The other major phyla of the moonmilk microbiome included Acidobacteria, Nitrospirae, Chloroflexi, Gemmatimonadetes, Planctomycetes, Latescibacteria, Verrucomicrobia, Zixibacteria, Armatimonadetes, Bacteroidetes, and Parcubacteria (Table S2, Figure 2a). Together, these phyla constituted 93.4%, 94.7%, and 91.5% of the total community in COL1, COL3, and COL4, respectively (Table S2, Figure 2a). The remaining phyla (with a relative abundance of <1%) were pooled as ’other’ (Figure 2a), and included most of the candidate divisions identified in this study (Figure 2b). Sequences that could not be affiliated to any bacterial phylum accounted for 4%, 3%, and 5% of the sequences in COL1, COL3, and COL4, respectively (Table S2). Some fraction of the moonmilk microbial diversity still remains to be discovered for all of the three sampling sites, as the rarefaction curves did not reach a plateau (Figure S2a).

2.2. Actinobacterial Diversity in Moonmilk Deposits

Evaluation of the actinobacterial profile was performed with libraries spanning the V6–V7 variable regions of 16S rRNA gene, and using modified Actinobacteria-specific primers (Table S1a). The specificity of the primers was confirmed by the detection of only 1%, 0.2%, and 2% of non-actinobacterial sequences in COL1, COL3, and COL4, respectively (Figure S3). In contrast to the bacterial dataset, the diversity of Actinobacteria appeared to be exhaustively sampled with the phylum-specific primers (Figure S2b).
The diversity indices for Actinobacteria showed the same trends as the ones observed for the whole Bacteria domain, i.e., evenness and diversity were the highest in COL4, followed by COL3, and COL1 (Table 1). Phylotype richness was the highest in COL4 with 211 OTUs, followed by COL1 and COL3 with 150 OTUs and 147 OTUs, respectively (Figure 1b and Table 1). Among the 243 different OTUs, 105 OTUs (43%) were found in all three of the studied moonmilk deposits (Figure 1b). Hence, the moonmilk-associated actinobacterial community appeared to be more conservative than the moonmilk-associated bacterial population (31%, Figure 1a). If we also include OTUs shared between at least two sampling points, the level of conservation rises to 66% of OTUs for Actinobacteria, and 58% for Bacteria. Still, 34% of the 243 OTUs (14, 15, and 54 OTUs in COL1, COL3, and COL4, respectively) remained specific to a moonmilk deposit, despite the close localization of collection points within the studied cave (Figure 1b). COL4 was characterized not only with the highest number of unique phylotypes (54 OTUs) (Figure 1b), but also with the most diverse population, as revealed by diversity indices (Table 1). As observed for the bacterial dataset, pairwise comparisons showed highly similar percentages (~36.4% ± 0.41%) of shared actinobacterial OTUs between moonmilk deposits (Table 2).
A taxonomic assignment of actinobacterial OTUs revealed the presence of two major classes—Acidimicrobiia and Actinobacteria, next to the low-abundant Thermoleophilia class (Table S3, Figure 3a). Acidimicrobiia was represented by one single order, the Acidimicrobiales, which dominated sample COL4, constituting 55.3% of the population (Table S3, Figure 3a). The Acidimicrobiales order consisted of two families, i.e., Acidimicrobiaceae and Iamiaceae (Table S3). The Actinobacteria class was represented by 15 orders, with Corynebacteriales dominating in COL1, and Pseudonocardiales in COL3 and COL4 (Table S3, Figure 3b). The most abundant families among the Actinobacteria class were Pseudonocardiaceae, Nocardiaceae, and Streptomycetaceae (Table S3). The proportion of unclassified and uncultured sequences at the family level ranged from 9% in COL3, to 25% in COL1, and 53% in COL4 (Table S3).
Among 28 families, a total of 47 genera were identified across the investigated samples (Table S3 and Table 3), with 35 genera identified for the first time in moonmilk (Table 3). COL1 was dominated by Rhodococcus (38.37%), while uncultured and unclassified Actinobacteria were the most abundant in COL3 and COL4 (Table 3). When only known genera were taken into account, Pseudonocardia prevailed in those samples, accounting for 20% and 18% of the population in COL3 and COL4, respectively (Table 3). Other genera, which constituted more than 1% of the population in at least one moonmilk deposit, included Streptomyces, Arthrobacter, Sporichthya, Planotetraspora, Nocardia, Mycobacterium, and Frankia (Table 3). While accounting in average for only 3% of the actinobacterial community, streptomycetes displayed the highest diversity, with 19 OTUs identified across the three moonmilk deposits (Table 3).
Some taxa showed important differences in their relative abundance between investigated samples, particularly Rhodococcus, which was approximately four and 14 times more abundant in COL1 than in COL3 and COL4, respectively (Table 3). The Streptomyces genus represented only 0.8% of the population in COL3, while it was detected at the level of 5.3% and 3% in the COL1 and COL4, respectively (Table 3). An important discrepancy in the relative abundance between speleothems was also observed for the genera Planotetraspora, Mycobacterium, and Frankia, whereas some taxa (e.g., Pseudoclavibacter, Lentzea, Propionibacterium) were exclusively found in a single sampling site (Table 3).

2.3. Analysis of the Most Abundant Actinobacterial OTUs

In order to obtain more information about the most dominant moonmilk-dwelling Actinobacteria, a detailed analysis was conducted for the 41 most abundant OTUs (~17% of all of the OTUs) accounting together for 90% (413,739 out of 456,878) of the sequences obtained via our HTS approach (Table 4). Out of the subset of 41 OTUs, 16 phylotypes belonged to the class Acidimicrobiia, with most of them being uncultured at the family level, and the remaining 25 OTUs belonged to the class Actinobacteria (Table 4). In the latter case, all of the OTUs were associated with major families previously identified in moonmilk deposits, including Pseudonocardiaceae, Propionibacteriaceae, Micrococcaceae, Nocardiaceae, Streptomycetaceae, and Streptosporangiaceae (Table 4). Only 16 OTUs could be classified at the genus level and were affiliated to genera Rhodococcus, Pseudonocardia, Arthrobacter, Sporichthya, Streptomyces, Planotetraspora, and Nocardia (Table 4).
Taking into account the spatial differences in terms of the most abundant taxa across the cave, COL1 was highly dominated by OTU1, affiliated to the genus Rhodococcus, and accounting for 38% of the total population in this speleothem (Table 4). This phylotype highly outnumbered other two Rhodococcus OTUs detected in COL1 (Table 3), which together constituted only 0.1% (data not shown). The predominant phylotypes identified in speleothems COL3 and COL4 were OTU2, representing an unclassified Pseudonocardiaceae in COL3 (29%), and OTU4, representing uncultured bacterium from Acidimicrobiia class in COL4 (11%) (Table 4). Among the known genera, Rhodococcus (OTU1, 10%) was also prevailing in COL3, while Pseudonocardia (OTU262, 6%) was found to be the most abundant in COL4 (Table 4).
In total, 40 out of 41 OTUs were present in all the three studied moonmilk deposits, often with an extreme variation in terms of their relative abundance across the different collection points. This is well demonstrated by OTU2 (Pseudonocardiaceae, unclassified at the genus level), which largely dominated the actinobacterial community in COL3 (29%), while only representing 0.3% of the actinobacterial microbiome in COL1 (Table 4).

2.4. Comparison of Moonmilk Streptomyces OTUs and Streptomyces Strains Isolated via the Culture-Dependent Approach

The true diversity of microbial communities is known to be strongly biased by cultivation-based methods in comparison to molecular techniques; therefore, we wanted to assess how much of the Streptomyces moonmilk-dwelling community we managed to isolate in our previous bioprospection work [10]. For this purpose, we compared the 16S rRNA sequences of the 19 Streptomyces OTUs retrieved from the HTS approach with the sequences of the 31 previously isolated Streptomyces phylotypes (MM strains), which were trimmed to the corresponding V6–V7 variable regions of HTS amplicons. Figure 4 presents the phylogenetic tree generated by maximum likelihood with all the 252 nt 16S rRNA sequences from the Streptomyces phylotypes (MM strains) and OTUs. The identity threshold for clustering sequences in the same branch of the tree was fixed to 97%, i.e., the same threshold as the one used to define OTUs in our HTS approach (see methods for details). As deduced from the generated phylogenetic tree, the 31 isolated Streptomyces strains matched with only five of the 19 Streptomyces OTUs, suggesting that the isolated strains represent a minor fraction of the Streptomyces species dwelling in the moonmilk deposits of the studied cave. Expectedly, Figure 4 further shows that we isolated Streptomyces species that are associated with the most abundant Streptomyces OTUs, e.g., OTU15, OTU21, OTU30, and OTU99 (Table 4), which together represent 79% of the Streptomyces sequences retrieved by our HTS approach. Moreover, 21 out of the 31 phylotype strains (68%) clustered together with OTU21 (Figure 4). Finally, two Streptomyces isolates, i.e., MM24 and MM106, did not cluster with any of the identified Streptomyces OTUs (Figure 4).

3. Discussion

3.1. New Insights into Moonmilk Bacterial Diversity Revealed by High-Throughput Sequencing

Previous investigations on the moonmilk microbiome revealed a very diverse microbial community in these deposits [9,13,14,15,16,17,19]. The high-throughput sequencing approach used in this work complemented previous findings by providing an in-depth picture of the bacterial population, together with a detailed taxonomic fingerprint of the phylum Actinobacteria.
Comparison of the bacterial diversity in moonmilk between earlier investigations and the present work is limited to some extent by the differences in experimental procedures, such as DNA isolation and PCR-based approaches, and the sensitivity of the sequencing techniques. Nonetheless, the profile of the major taxonomic groups found in this work is consistent with that observed for the moonmilk communities in the caves “Grotta della Foos” and “Bus della Genziana” in Italy, which were obtained from 16S rRNA clone libraries [9]. All of the phyla detected in the above-mentioned caves, including Bacteroidetes, Acidobacteria, Chloroflexi, Planctomycetes, Verrucomicrobia, Actinobacteria, Firmicutes, Nitrospirae, Chlorobi, Proteobacteria, and WS3 (now Latescibacteria), were also identified in the cave “Grotte des Collemboles”, although their relative abundance varied between the studies. While Proteobacteria were found to be the most abundant phylum in both cases, the second most abundant population identified in Italian caves was the phylum Bacteroidetes, which constituted a minor part of the bacterial community in the present study. The Actinobacteria population was found to be an important part of the moonmilk microbiome in the “Grotte des Collemboles” (from 9% to 23%), but instead represented only a minor fraction (<2%) of the bacterial population in the two Italian caves investigated by Engel et al. (2013). Very recently, a study by Dhami et al. has reported the moonmilk microbiome profile in the Australian “Lake Cave” using an HTS approach [19], as in this work. The presence of Proteobacteria, Actinobacteria, Acidobacteria, Chloroflexi, Nitrospirae, Gemmatimonadetes, Firmicutes, and Bacteroidetes were detected in the moonmilk deposit of the “Lake Cave”, similarly to the “Grotte des Collemboles”. However, many of the low-abundance taxa identified in the Belgian cave were not reported, possibly because their phylogenetic profiles were based on different regions of 16S rRNA gene—V3/V4 for the “Lake Cave”, and V6–V7 for the “Grotte des Collemboles”. Interestingly, unlike in Italian and Belgian caves, the “Lake Cave” moonmilk deposit was strongly dominated by Actinobacteria, which were more than twice as abundant as the Proteobacteria [19]. The highly sensitive HTS amplicon sequencing approach employed in this work revealed the presence of 26 phyla within the moonmilk microbiome that had not been previously described in this speleothem. These included Zixibacteria (formerly RBG-1), Armatimonadetes (formerly OP10), and Parcubacteria (formerly OD1) among the main phyla of moonmilk microbiome (Figure 2a), which have been previously reported from other subterranean environments [24,29,30,31,32,33], and 23 low-abundant taxa that were found below the level of 1%, and included many candidate divisions (Figure 2b).
This new study uncovered a surprisingly diverse Actinobacteria taxonomic profile that demonstrates the limitations of our previous cultivation-based screening, in which only the Streptomyces species could be isolated from the three moonmilk deposits [10]. Here, a total of 47 actinobacterial genera from 28 families were identified across the investigated samples. Beyond the previously reported members of the Actinomycetales family—including Nocardia and Rhodococcus (Nocardiaceae) [15,18], Pseudonocardia, Amycolatopsis and Saccharothrix (not identified in our study) (Pseudonocardiaceae) [12,14,19], Propionibacterium (Propionibacteriaceae) [14], Streptomyces (Streptomycetaceae) [10,12,18,19], Arthrobacter (Micrococcaceae) [13], Mycobacterium (Mycobacteriaceae) [19], Nocardioides, Aeromicrobium, and Kribbella (Nocardioidaceae) [19], and Geodermatophilus (Geodermatophilaceae) [19]—35 other genera were identified in the moonmilk deposits of the “Grotte des Collemboles”. The population of each investigated sample was also found to include representatives of the Acidimicrobiia class, which were not previously reported in moonmilk. Their presence in all the three sampling sites, with an abundance up to 55% in COL4, and the dominance of the unclassified Acidimicrobiia phylotype (OTU 4) within the community of COL4, suggest that the chemical composition of the investigated moonmilk would be particularly suitable for the development of the representatives of this class of Actinobacteria, of which the ecology and metabolism are still largely unknown.

3.2. Moonmilk Deposits as Appealing Source of Novel Producers of Bioactive Compounds

Extreme environmental niches have recently become the main targets for intense bioprospecting, as they are expected to host diverse yet-unknown microorganisms, which could offer unexplored chemical diversity. While Streptomyces are reported as the most prolific “antibiotic makers”, advances in the cultivation and characterization of rare Actinobacteria revealed similarly promising capabilities for the production of bioactive natural compounds [34,35,36]. The results obtained in this work suggest a significant biodiversity of the moonmilk-dwelling actinobacterial population, with a wide spectrum of rare genera. Next to Streptomyces, other members of Actinobacteria with valuable secondary metabolism were detected at a high proportion, such as Pseudonocardia, Amycolatopsis, Streptosporangium, Nocardia, Nocardioides, and Rhodococcus. Such findings clearly prompt to apply appropriate selective cultivation methods to isolate rare Actinobacteria from moonmilk deposits.
Moreover, particular importance should be also focused on Acidimicrobiia, which constituted an important part of the community in the studied deposits. Members of this class are a recently identified taxonomic unit [37] that is considered to represent an early-branching lineage within the phylum [38]. Due to their phylogenetic isolation and novelty, they are likely to hide a yet-uncovered valuable bioactive arsenal.
The great potential of moonmilk as a source of diverse and metabolically beneficial Actinobacteria is illustrated by the comparison of Streptomyces isolated in our previous study and the Streptomyces OTUs identified in this work (Figure 4). Most Streptomyces OTUs are phylogenetically distinct from culturable representatives (Figure 4), indicating that a great number of species still remain to be isolated. On the other hand, our culture-dependent study identified Streptomyces strains (MM24 and MM106, Figure 4) that were not associated with OTUs deduced from the HTS approach, confirming that both strategies are complementary, and should be used in parallel for microbial diversity assessment [39,40]. In addition, next to the identification methods themselves, our data suggests that the diversity level can be also biased by the identity threshold that is used for OTU definition. The tree revealed that a single OTU (OTU21, Figure 4) clustered together with most of the phylotypes deduced from MLSA (multilocus sequence analysis), each most likely representing a distinct species [10]. This indicates that the 97% sequence homology threshold applied to the comparative analysis of the V6–V7 regions of the 16S RNA gene largely underestimated the number of Streptomyces species dwelling in a studied environmental niche.

4. Materials and Methods

4.1. Site description and Sampling

The cave “Grotte des Collemboles” (Springtails’ Cave), located in Comblain-au-Pont (GPS coordinates 50°28′41″ N, 5°36′35″ E), Belgium (Figure S1, Maciejewska et al. 2017 for full description), was formed in Visean limestone and has the shape of a 70-m long meander. White to brown–orange moonmilk deposits are found on the walls in the first narrow chamber located at the entrance of the cave, as well as in the narrow passages leading deeper into the cave (Figure S1). Moonmilk samples used for total DNA extractions were aseptically collected in January 2012 from three spatially separated locations along about a 20-m transect in the cave. Soft moonmilk speleothems were scratched with sterile scalpels into sterile Falcon tubes from the wall in the first chamber, adjacent to the cave entrance (COL4), and from the walls in a narrow passage after the first chamber (COL1, COL3) (Figure S1). COL4 was located approximately 6 m from COL1, and 20 m from COL3 (Figure S1). Samples were immediately transferred to the laboratory, freeze-dried on a VirTis Benchtop SLC Lyophilizer (SP Scientific, Warminster, PA, USA), and stored at −20 °C.

4.2. Total DNA Extraction and 16S rRNA Gene Amplicon High-Throughput Sequencing

The metagenetic approach applied in this work was performed on DNA extracted from three moonmilk deposits (COL1, COL3, and COL4) originating from the “Grotte des Collemboles”. Environmental genomic DNA isolation was carried out from 200 mg of the freeze-dried moonmilk samples COL1, COL3, and COL4 (Figure S1), using the PowerClean Soil DNA kit (MoBio, Carlsbad, CA, USA), according to manufacturer’s instructions. The integrity of purified DNA was assessed by agarose gel electrophoresis (1% w/v), and the dsDNA concentration was evaluated by Qubit fluorometer (Invitrogen, Carlsbad, CA, USA).
The 16S rRNA gene amplicon libraries were generated using bacterial (S-D-Bact-0517-a-S-17/S-D-Bact-1061-a-A-17 spanning V4–V6 region [41]) and actinobacterial (Com2xf/Ac1186r, spanning V6–V7 region [42]) specific primer pairs. The Illumina platform-compatible dual index paired-end approach was designed as previously described [43] (detailed description provided in the Table S1a). Each forward and reverse primer consisted of an Illumina-compatible forward/reverse primer overhang attached to the 5′ end. Additionally, a heterogeneity spacer of four degenerate nucleotides (Ns) was added to the forward primer, between the primer overhang and the locus-specific sequence. The Illumina barcodes and sequencing adapters were added during the subsequent cycle-limited amplification step using Nextera XT Index kit (Illumina, San Diego, CA, USA). Triplicated PCR reactions were performed for each sample in 25 µL of volume containing 2.5 µL of total DNA, 5 µL of each primer (1 µM), and 12.5 µL of 2× Q5 High-Fidelity Master Mix (New England Biolabs, Ipswich, MA, USA). Amplification conditions for each set of primers are listed in Table S1b. The triplicated amplicons were visualized on 3% agarose gel, pulled, purified with Agencourt AMPure XP beads (Beckman Coulter, Brea, CA, USA), and quantified with the Qubit HS dsDNA assay kit (Invitrogen, Carlsbad, CA,) before being processed for index ligation, using the Nextera XT Index kit (Illumina, San Diego, CA, USA). The PCR amplifications were performed with the same enzyme and cycling conditions as described above [43], with the total number of cycles reduced to eight, and an annealing temperature of 55 °C. The resulting amplicons were purified with the Agencourt AMPure XP magnetic beads (Beckman Coulter, CA, USA), quantified, and pooled in equimolar concentrations. The library concentration was quantified by qPCR using a Kappa SYBR FAST kit (Kapa Biosystems, Wilmington, MA, USA), and subsequently, the library was normalized to 4 nM, denaturated, and diluted to the final concentration of 8 p.m. The resulting pool was mixed with the PhiX control and subjected to 2 × 300 bp paired-end sequencing on Illumina MiSeq platform (Illumina, San Diego, CA, USA). Raw sequences were deposited in the NCBI Sequence Read Archive (SRA) database under the Bioproject PRJNA428798 with accession numbers SRX3540524–SRX3540529.

4.3. 16S rRNA Amplicon Analysis

16S rRNA amplicon analysis was based for both Bacteria and Actinobacteria on forward reads only, owing to the poor quality of reverse reads. Quality trimming (prohibiting mismatches and ambiguities, ensuring a minimum quality score of 20 and removing the four degenerate nucleotides from the 5′ end) was carried out using CLC Genomic Workbench (Qiagen, Hilden, Germany). USEARCH [44] was applied for length trimming (minimum length = 240 nt) and dereplication. Operational taxonomic units (OTUs) for both bacterial and actinobacterial datasets were defined using a 97% identity threshold on 16S rRNA sequences. OTUs were clustered using the UPARSE algorithm [45], and their taxonomic position was assigned by MOTHUR [46] with SILVA v128 database [47]. OTUs were further classified using BLASTN [48] analyses against a local mirror of NCBI nt database (downloaded on 9 August 2017), through manual and automatic analyses. For the automatic approach, a last common ancestor (LCA) classification was performed with a custom parser mimicking the MEGAN algorithm [49], which we developed for analyses of genome contamination (Cornet et al., 2017, under review). A maximum number of 100 hits per OTU were taken into account. To consider a BLASTN hit, the E-value threshold was set at 1e-15, the minimum identity threshold was set at 95.5%, the minimum bit score was set at 200, and the bit score percentage threshold was set at 99% of the best hit. These thresholds were defined through preliminary analyses (data not shown). When the BLASTN hits are too numerous, the MEGAN-like algorithm frequently yields high-ranking LCAs (e.g., Bacteria) that are not informative in practice. In order to minimize this effect, we decided to skip uncultured/unclassified hits whenever other, more informative, hits also passed the thresholds. Moreover, when computing LCAs, we only considered the most frequent taxa, provided that they represented ≥95% of the (up to 100) accumulated BLASTN hits, so as to avoid uninformative classifications due to a few (possibly aberrant) outliers.
Normalized OTU abundance data was used to calculate α-diversity and β-diversity estimators using MOTHUR [46]. Community richness, evenness, diversity, and differential OTU abundance between samples were calculated using sobs, the Simpson index, the inverse Simpson index and Venn diagrams, respectively.
The 19 OTUs identified as Streptomyces were combined to 31 sequences (16S rRNA region V6–V7) from previously isolated Streptomyces phylotypes (MM strains) and dereplicated with the UCLUST algorithm [44] using an identity threshold of 97%. This yielded 21 clusters, to which we added the homologous region of Corynebacterium diphtheriae JCM-1310 as an outgroup. A multiple sequence alignment was built with MUSCLE [50] (default parameters), and then analyzed with PhyML [51] under a K80 + Γ4 model. Due to the limited amount of phylogenetic signal (short sequences from very related organisms), the resolution of the tree was low (bootstrap proportions <50 for nearly all nodes; data not shown).

5. Conclusions

Before the advent of metagenomics, bioprospecting was carried out blindly, with poor knowledge on the real potential of an ecological niche mined for novel organisms, enzymes, or bioactive compounds. The results of the metagenetic study presented here confirmed that different moonmilk deposits host their own indigenous microbial population, and thus each individual speleothem can be a source of a great biodiversity. Consequently, the observed important differences in the spatial diversity of Actinobacteria imply that bioprospecting within different moonmilk deposits—from different caves or within the same cave—could result in the isolation of unique and novel natural compound producers. Our study also revealed how many and which actinobacterial genera have been missed in our first attempt to isolate antibiotic producers. We now know that the Streptomyces strains of our collection isolated from the moonmilk deposits of the cave ‘Grotte des Collemboles’ [10] are just the tip of the iceberg. These results prompted us to apply a series of ‘tips and tricks’ to isolate other Streptomyces and representatives of other antibiotic-producing Actinobacteria that are present in different proportions in each moonmilk deposit. The results of our adapted protocols for the isolation of rare Actinobacteria are presented in the article ‘Isolation, Characterization, and Antibacterial Activity of Hard-to-Culture Actinobacteria from Cave Moonmilk Deposits’, which is published in the same special issue [52].

Supplementary Materials

The following are available online at https://www.mdpi.com/2079-6382/7/2/27/s1, Figure S1: Localization of the “Grotte des Collemboles” (Springtails’ Cave) together with the cave map and visualization of the moonmilk deposit sampling points, Figure S2: Rarefaction curves of OTUs clustered at 97% sequence identity across the three moonmilk-sampling points for Bacteria (a) and Actinobacteria (b), Figure S3: Taxonomic profile of bacterial phyla generated with Actinobacteria-specific primers. Note the high specificity of Actinobacteria primers, Table S1: Details of the PCR primers used for community profiling of moonmilk samples (a) and—PCR conditions used for 16S rRNA amplification from moonmilk samples (b), Table S2: Relative abundance (%) of bacterial phyla identified in the three moonmilk deposits in “Grotte des Collemboles”. Low-abundant taxa with relative abundance <1% are marked in red, Table S3: Relative abundance (%) of the phylum Actinobacteria at different taxonomic levels identified in the three moonmilk deposits in the “Grotte des Collemboles”.

Acknowledgments

The authors are grateful to Luc Willems for the introduction to the subject and help with sampling. MM, DA, and LC work was supported by a Research Foundation for Industry and Agriculture (FRIA) grant. MC and PD were supported by the Luxembourg National Research Fund (FNR CORE 2011 project GASPOP; C11/SR/1280949: Influence of the Reactor Design and the Operational Parameters on the Dynamics of the Microbial Consortia Involved in the Biomethanation Process). Computational resources (“durandal” grid computer) were funded by three grants from the University of Liège, “Fonds spéciaux pour la recherche”, “Crédit de démarrage 2012” (SFRD-12/03 and SFRD-12/04) and “Crédit classique 2014” (C-14/73) and by a grant from the F.R.S.-FNRS “Crédit de recherche 2014” (CDR J.0080.15). This work is supported in part by the Belgian program of Interuniversity Attraction Poles initiated by the Federal Office for Scientific Technical and Cultural Affairs (PAI No. P7/44). SR is a Research Associate at Belgian Fund for Scientific Research (F.R.S-FNRS). The authors declare no conflict of interest. We dedicate the work to the memory of Leonard Maculewicz (1936–2017) who always supported our work with great enthusiasm.

Author Contributions

M.M., Ma.C., S.M., Mo.C., P.D., and S.R. designed and performed experiments. Bioinformatic analyses were performed by M.M., Mo.C., Ma.C., L.C., D.B., and S.R. Data were analyzed by all authors. The manuscript was written and/or corrected by all authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Engel, A.S. Microbial diversity of cave ecosystems. In Geomicrobiology: Molecular and Environmental Perspective; Springer: Berlin, Germany, 2010; ISBN 9789048192038. [Google Scholar]
  2. Barton, H.A.; Northup, D.E. Geomicrobiology in cave environments: Past, current and future perspectives. J. Cave Karst Stud. 2007, 69, 163–178. [Google Scholar]
  3. Groth, I.; Vettermann, R.; Schuetze, B.; Schumann, P.; Saiz-Jimenez, C. Actinomycetes in Karstic caves of northern Spain (Altamira and Tito Bustillo). J. Microbiol. Methods 1999, 36, 115–122. [Google Scholar] [CrossRef]
  4. Jurado, V.; Kroppenstedt, R.M.; Saiz-Jimenez, C.; Klenk, H.P.; Mouniée, D.; Laiz, L.; Couble, A.; Pötter, G.; Boiron, P.; Rodríguez-Nava, V. Hoyosella altamirensis gen. nov., sp. nov., a new member of the order Actinomycetales isolated from a cave biofilm. Int. J. Syst. Evol. Microbiol. 2009, 59, 3105–3110. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Jurado, V.; Groth, I.; Gonzalez, J.M.; Laiz, L.; Saiz-Jimenez, C. Agromyces salentinus sp. nov. and Agromyces neolithicus sp. nov. Int. J. Syst. Evol. Microbiol. 2005, 55, 153–157. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Maciejewska, M.; Pessi, I.S.; Arguelles-Arias, A.; Noirfalise, P.; Luis, G.; Ongena, M.; Barton, H.; Carnol, M.; Rigali, S. Streptomyces lunaelactis sp. nov., a novel ferroverdin A-producing Streptomyces species isolated from a moonmilk speleothem. Antonie van Leeuwenhoek, Int. J. Gen. Mol. Microbiol. 2015, 107, 519–531. [Google Scholar] [CrossRef] [PubMed]
  7. Lee, S.D.; Kang, S.O.; Hah, Y.C. Hongia gen, nov., a new genus of the order Actinomycetales. Int. J. Syst. Evol. Microbiol. 2000, 50, 191–199. [Google Scholar] [CrossRef] [PubMed]
  8. Ortiz, M.; Neilson, J.W.; Nelson, W.M.; Legatzki, A.; Byrne, A.; Yu, Y.; Wing, R.A.; Soderlund, C.A.; Pryor, B.M.; Pierson, L.S.; et al. Profiling Bacterial Diversity and Taxonomic Composition on Speleothem Surfaces in Kartchner Caverns, AZ. Microb. Ecol. 2013, 65, 371–383. [Google Scholar] [CrossRef] [PubMed]
  9. Engel, A.S.; Paoletti, M.G.; Beggio, M.; Dorigo, L.; Pamio, A.; Gomiero, T.; Furlan, C.; Brilli, M.; Dreon, A.L.; Bertoni, R.; et al. Comparative microbial community composition from secondary carbonate (moonmilk) deposits: Implications for the Cansiliella servadeii cave hygropetric food web. Int. J. Speleol. 2013, 42, 181–192. [Google Scholar] [CrossRef]
  10. Maciejewska, M.; Adam, D.; Martinet, L.; Naômé, A.; Calusinska, M.; Smargiasso, N.; De Pauw, E.; Barton, H.; Carnol, M.; Hanikenne, M.; et al. A Phenotypic and Genotypic Analysis of the Antimicrobial Potential of Cultivable Streptomyces isolated from Cave Moonmilk Deposits. Front. Microbiol. 2016. [Google Scholar] [CrossRef] [PubMed]
  11. Cañaveras, J.C.; Cuezva, S.; Sanchez-Moral, S.; Lario, J.; Laiz, L.; Gonzalez, J.M.; Saiz-Jimenez, C. On the origin of fiber calcite crystals in moonmilk deposits. Naturwissenschaften 2006, 93, 27–32. [Google Scholar] [CrossRef] [PubMed]
  12. Cañaveras, J.C.; Hoyos Gómez, M.; Sánchez-Moral, S.; Sanz Rubio, E.; Bedoya, J.; Hoyos, V.; Groth, I.; Schumann, P.; Laiz Trobajo, L. Microbial Communities Associated with Hydromagnesite and Needle-Fiber Aragonite Deposits in a Karstic Cave (Altamira, Northern Spain). Geomicrobiol. J. 1999, 16, 9–25. [Google Scholar]
  13. Rooney, D.C.; Hutchens, E.; Clipson, N.; Baldini, J.; McDermott, F. Microbial Community Diversity of Moonmilk Deposits at Ballynamintra Cave, Co. Waterford, Ireland. Microb. Ecol. 2010, 60, 753–761. [Google Scholar] [CrossRef] [PubMed]
  14. Portillo, M.C.; Gonzalez, J.M. Moonmilk Deposits Originate from Specific Bacterial Communities in Altamira Cave (Spain). Microb. Ecol. 2011, 61, 182–189. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Reitschuler, C.; Lins, P.; Schwarzenauer, T.; Spotl, C.; Wagner, A.O.; Illmer, P. New Undescribed Lineages of Non-extremophilic Archaea Form a Homogeneous and Dominant Element within Alpine Moonmilk Microbiomes. Geomicrobiol. J. 2015, 32. [Google Scholar] [CrossRef]
  16. Reitschuler, C.; Spötl, C.; Hofmann, K.; Wagner, A.O.; Illmer, P. Archaeal Distribution in Moonmilk Deposits from Alpine Caves and Their Ecophysiological Potential. Microb. Ecol. 2016, 71, 686–699. [Google Scholar] [CrossRef] [PubMed]
  17. Reitschuler, C.; Lins, P.; Wagner, A.O.; Illmer, P. Cultivation of moonmilk-born non-extremophilic Thaum and Euryarchaeota in mixed culture. Anaerobe 2014, 29, 73–79. [Google Scholar] [CrossRef] [PubMed]
  18. Axenov-Gibanov, D.V.; Voytsekhovskaya, I.V.; Tokovenko, B.T.; Protasov, E.S.; Gamaiunov, S.V.; Rebets, Y.V.; Luzhetskyy, A.N.; Timofeyev, M.A. Actinobacteria isolated from an underground lake and moonmilk speleothem from the biggest conglomeratic karstic cave in Siberia as sources of novel biologically active compounds. PLoS ONE 2016, 11. [Google Scholar] [CrossRef] [PubMed]
  19. Dhami, N.K.; Mukherjee, A.; Watkin, E. Characterisation of Mineralogical-Mechanical-Microbial properties of calcitic speleothems and the in vitro biomineralization potential of associated microbial communities. Front. Microbiol. 2018, 9, 40. [Google Scholar] [CrossRef] [PubMed]
  20. Bindschedler, S.; Milliere, L.; Cailleau, G.; Job, D.; Verrecchia, E.P. Calcitic nanofibres in soils and caves: A putative fungal contribution to carbonatogenesis. Geol. Soc. Lond. Spec. Publ. 2010, 336, 225–238. [Google Scholar] [CrossRef]
  21. Maciejewska, M.; Adam, D.; Naômé, A.; Martinet, L.; Tenconi, E.; Calusinska, M.; Delfosse, P.; Hanikenne, M.; Baurain, D.; Compère, P.; et al. Assessment of the potential role of Streptomyces in cave moonmilk formation. Front. Microbiol. 2017, 8. [Google Scholar] [CrossRef] [PubMed]
  22. Nimaichand, S.; Devi, A.M.; Tamreihao, K.; Ningthoujam, D.S.; Li, W.J. Actinobacterial diversity in limestone deposit sites in Hundung, Manipur (India) and their antimicrobial activities. Front. Microbiol. 2015, 6. [Google Scholar] [CrossRef] [PubMed]
  23. Wu, Y.; Tan, L.; Liu, W.; Wang, B.; Wang, J.; Cai, Y.; Lin, X. Profiling bacterial diversity in a limestone cave of the western Loess Plateau of China. Front. Microbiol. 2015, 6. [Google Scholar] [CrossRef] [PubMed]
  24. Northup, D.E.; Melim, L.A.; Spilde, M.N.; Hathaway, J.J.M.; Garcia, M.G.; Moya, M.; Stone, F.D.; Boston, P.J.; Dapkevicius, M.L.N.E.; Riquelme, C. Lava Cave Microbial Communities within Mats and Secondary Mineral Deposits: Implications for Life Detection on Other Planets. Astrobiology 2011, 11, 601–618. [Google Scholar] [CrossRef] [PubMed]
  25. Riquelme, C.; Hathaway, J.J.M.; Dapkevicius, M.d.L.N.E.; Miller, A.Z.; Kooser, A.; Northup, D.E.; Jurado, V.; Fernandez, O.; Saiz-Jimenez, C.; Cheeptham, N. Actinobacterial diversity in volcanic caves and associated geomicrobiological interactions. Front. Microbiol. 2015, 6. [Google Scholar] [CrossRef] [PubMed]
  26. Cheeptham, N.; Sadoway, T.; Rule, D.; Watson, K.; Moote, P.; Soliman, L.C.; Azad, N.; Donkor, K.K.; Horne, D. Cure from the cave: Volcanic cave actinomycetes and their potential in drug discovery. Int. J. Speleol. 2013, 42, 35–47. [Google Scholar] [CrossRef]
  27. Tebo, B.M.; Davis, R.E.; Anitori, R.P.; Connell, L.B.; Schiffman, P.; Staudigel, H. Microbial communities in dark oligotrophic volcanic ice cave ecosystems of Mt. Erebus, Antarctica. Front. Microbiol. 2015, 6. [Google Scholar] [CrossRef] [PubMed]
  28. Chater, K.F.; Biró, S.; Lee, K.J.; Palmer, T.; Schrempf, H. The complex extracellular biology of Streptomyces. FEMS Microbiol. Rev. 2010, 34, 171–198. [Google Scholar] [CrossRef] [PubMed]
  29. Castelle, C.J.; Hug, L.A.; Wrighton, K.C.; Thomas, B.C.; Williams, K.H.; Wu, D.; Tringe, S.G.; Singer, S.W.; Eisen, J.A.; Banfield, J.F. Extraordinary phylogenetic diversity and metabolic versatility in aquifer sediment. Nat. Commun. 2013, 4. [Google Scholar] [CrossRef] [PubMed]
  30. Schabereiter-Gurtner, C.; Saiz-Jimenez, C.; Pinar, G.; Lubitz, W.; Rolleke, S. Phylogenetic diversity of bacteria associated with Paleolithic paintings and surrounding rock walls in two Spanish caves (Llonin and La Garma). FEMS Microbiol. Ecol. 2004, 47, 235–247. [Google Scholar] [CrossRef]
  31. Pedersen, K.; Arlinger, J.; Ekendahl, S.; Hallbeck, L. 16S rRNA gene diversity of attached and unattached bacteria in boreholes along the access tunnel to the Äspö hard rock laboratory, Sweden. FEMS Microbiol. Ecol. 1996, 19, 249–262. [Google Scholar] [CrossRef]
  32. Zhou, J.; Gu, Y.; Zou, C.; Mo, M. Phylogenetic diversity of bacteria in an earth-cave in Guizhou Province, Southwest of China. J. Microbiol. 2007, 45, 105–112. [Google Scholar] [PubMed]
  33. De Mandal, S.; Chatterjee, R.; Kumar, N.S. Dominant bacterial phyla in caves and their predicted functional roles in C and N cycle. BMC Microbiol. 2017, 17, 90. [Google Scholar] [CrossRef] [PubMed]
  34. Tiwari, K.; Gupta, R.K. Rare actinomycetes: A potential storehouse for novel antibiotics. Crit. Rev. Biotechnol. 2012, 32, 108–132. [Google Scholar] [CrossRef] [PubMed]
  35. Choi, S.S.; Kim, H.J.; Lee, H.S.; Kim, P.; Kim, E.S. Genome mining of rare actinomycetes and cryptic pathway awakening. Process Biochem. 2015, 50, 1184–1193. [Google Scholar] [CrossRef]
  36. Rigali, S.; Anderssen, S.; Naômé, A.; van Wezel, G.P. Cracking the regulatory code of biosynthetic gene clusters as a strategy for natural product discovery. Biochem. Pharmacol. 2018. [Google Scholar] [CrossRef] [PubMed]
  37. Norris, P. Class II. Acidimicrobiia class. nov. In Bergey’s Manual of Systematic Bacteriology; Whitman, W., Goodfellow, M., Kämpfer, P., Busse, H.-J., Trujillo, M., Ludwig, W., Suzuki, K.-I., Parte, A., Eds.; Springer: New York, NY, USA, 2012; Volume 5, p. 1968. [Google Scholar]
  38. Bull, A.T. Actinobacteria of the extremobiosphere. In Extremophiles Handbook; Springer: Tokyo, Japan, 2011; pp. 1203–1240. ISBN 9784431538981. [Google Scholar]
  39. Stefani, F.O.P.; Bell, T.H.; Marchand, C.; de la Providencia, I.E.; El Yassimi, A.; St-Arnaud, M.; Hijri, M. Culture-Dependent and -Independent Methods Capture Different Microbial Community Fractions in Hydrocarbon-Contaminated Soils. PLoS ONE 2015, 10, e0128272. [Google Scholar] [CrossRef] [PubMed]
  40. Vaz-Moreira, I.; Egas, C.; Nunes, O.C.; Manaia, C.M. Culture-dependent and culture-independent diversity surveys target different bacteria: A case study in a freshwater sample. Antonie van Leeuwenhoek, Int. J. Gen. Mol. Microbiol. 2011, 100, 245–257. [Google Scholar] [CrossRef] [PubMed]
  41. Klindworth, A.; Pruesse, E.; Schweer, T.; Peplies, J.; Quast, C.; Horn, M.; Glöckner, F.O. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2013, 41. [Google Scholar] [CrossRef] [PubMed]
  42. Schäfer, J.; Jäckel, U.; Kämpfer, P. Development of a new PCR primer system for selective amplification of Actinobacteria. FEMS Microbiol. Lett. 2010, 311, 103–112. [Google Scholar] [CrossRef] [PubMed]
  43. Goux, X.; Calusinska, M.; Fossépré, M.; Benizri, E.; Delfosse, P. Start-up phase of an anaerobic full-scale farm reactor—Appearance of mesophilic anaerobic conditions and establishment of the methanogenic microbial community. Bioresour. Technol. 2016, 212, 217–226. [Google Scholar] [CrossRef] [PubMed]
  44. Edgar, R.C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 2010, 26, 2460–2461. [Google Scholar] [CrossRef] [PubMed]
  45. Edgar, R.C. UPARSE: Highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 2013, 10, 996–998. [Google Scholar] [CrossRef] [PubMed]
  46. Schloss, P.D.; Westcott, S.L.; Ryabin, T.; Hall, J.R.; Hartmann, M.; Hollister, E.B.; Lesniewski, R.A.; Oakley, B.B.; Parks, D.H.; Robinson, C.J.; et al. Introducing mothur: Open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 2009, 75, 7537–7541. [Google Scholar] [CrossRef] [PubMed]
  47. Pruesse, E.; Peplies, J.; Glöckner, F.O. SINA: Accurate high-throughput multiple sequence alignment of ribosomal RNA genes. Bioinformatics 2012, 28, 1823–1829. [Google Scholar] [CrossRef] [PubMed]
  48. Boratyn, G.M.; Schäffer, A.A.; Agarwala, R.; Altschul, S.F.; Lipman, D.J.; Madden, T.L. Domain enhanced lookup time accelerated BLAST. Biol. Direct 2012, 7. [Google Scholar] [CrossRef] [PubMed]
  49. Huson, D.H.; Auch, A.F.; Qi, J.; Schuster, S.C. MEGAN analysis of metagenomic data. Genome Res. 2007, 17, 377–386. [Google Scholar] [CrossRef] [PubMed]
  50. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [PubMed]
  51. Guindon, S.; Dufayard, J.F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O. New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0. Syst. Biol. 2010, 59, 307–321. [Google Scholar] [CrossRef] [PubMed]
  52. Adam, D.; Maciejewska, M.; Naômé, A.; Martinet, L.; Coppieters, W.; Karim, L.; Baurain, D.; Rigali, S. Isolation, Characterization, and Antibacterial Activity of Hard-to-Culture Actinobacteria from Cave Moonmilk Deposits. Antibiotics 2018, in press. [Google Scholar] [CrossRef]
Antibiotics 07 00027 g001 550
Figure 1. Venn diagrams showing the numbers of shared and unique bacterial (a) and actinobacterial (b) OTUs between the three moonmilk sampling points (COL1, COL3, COL4).
Figure 1. Venn diagrams showing the numbers of shared and unique bacterial (a) and actinobacterial (b) OTUs between the three moonmilk sampling points (COL1, COL3, COL4).
Antibiotics 07 00027 g001
Antibiotics 07 00027 g002 550
Figure 2. Taxonomic profiles of the moonmilk-associated microbiome at the phylum level across the three moonmilk sampling points (COL1, COL3, COL4). The main phyla of the microbiome are presented on the left (a), while the pattern of low-abundance taxa, named as ‘other’ (with a relative abundance of <1%) is displayed on the right (b).
Figure 2. Taxonomic profiles of the moonmilk-associated microbiome at the phylum level across the three moonmilk sampling points (COL1, COL3, COL4). The main phyla of the microbiome are presented on the left (a), while the pattern of low-abundance taxa, named as ‘other’ (with a relative abundance of <1%) is displayed on the right (b).
Antibiotics 07 00027 g002
Antibiotics 07 00027 g003 550
Figure 3. Taxonomic profiles of moonmilk-associated Actinobacteria at different taxonomic levels—(a) class; (b) order; (c) family—observed across the three moonmilk-sampling points (COL1, COL3, COL4). ‘Other’ includes orders and families with a relative abundance of <1%.
Figure 3. Taxonomic profiles of moonmilk-associated Actinobacteria at different taxonomic levels—(a) class; (b) order; (c) family—observed across the three moonmilk-sampling points (COL1, COL3, COL4). ‘Other’ includes orders and families with a relative abundance of <1%.
Antibiotics 07 00027 g003
Antibiotics 07 00027 g004 550
Figure 4. Phylogenetic relationships between culturable and non-culturable Streptomyces originating from moonmilk of “Grotte des Collemboles”. The tree was inferred by maximum likelihood. Scale bar is in substitution per site. Numbers between brackets reflect the predicted mean abundance of Streptomyces OTUs in the studied deposits based on the percentage of sequences retrieved from the HTS analysis. Streptomyces phylotypes isolated in our previous bioprospection study (MM strains) are marked in blue.
Figure 4. Phylogenetic relationships between culturable and non-culturable Streptomyces originating from moonmilk of “Grotte des Collemboles”. The tree was inferred by maximum likelihood. Scale bar is in substitution per site. Numbers between brackets reflect the predicted mean abundance of Streptomyces OTUs in the studied deposits based on the percentage of sequences retrieved from the HTS analysis. Streptomyces phylotypes isolated in our previous bioprospection study (MM strains) are marked in blue.
Antibiotics 07 00027 g004
Table
Table 1. Richness, specificity, diversity, and evenness of the bacterial and actinobacterial communities in the three moonmilk deposits of the “Grotte des Collemboles”. OTUs: operational taxonomic units.
Table 1. Richness, specificity, diversity, and evenness of the bacterial and actinobacterial communities in the three moonmilk deposits of the “Grotte des Collemboles”. OTUs: operational taxonomic units.
Target groupSiteTotal OTUs (Richness)Unique OTUs (Specificity)Inverse Simpson Index (Diversity)Simpson Index (Evenness)
BacteriaCOL11332238 (17.9%)13.230.01
COL31161134 (11.6%)58.940.05
COL41863584 (31.3%)155.310.08
ActinobacteriaCOL115014 (9.3%)6.210.04
COL314715 (10.2%)7.740.05
COL421154 (25.6%)24.130.11
Table
Table 2. Pairwise comparisons of shared OTUs between the moonmilk deposits.
Table 2. Pairwise comparisons of shared OTUs between the moonmilk deposits.
Target groupCOL1 and COL3COL1 and COL4COL4 and COL3
Bacteria776/2493 (31.1%)1028/3195 (32.2%)961/3024 (31.8%)
Actinobacteria108/297 (36.4%)133/361 (36.9%)129/358 (36.0%)
Table
Table 3. Actinobacterial genera pattern in moonmilk deposits of the “Grotte des Collemboles” based on 16S rRNA amplicon libraries.
Table 3. Actinobacterial genera pattern in moonmilk deposits of the “Grotte des Collemboles” based on 16S rRNA amplicon libraries.
GenusCOL1COL3COL4TOTAL
Seq.%OTUsSeq.%OTUsSeq.%OTUsSeq.Av. %Diff. OTUs
Uncultured46,34632.386558,90334.806392,44463.9890197,69343.796
Rhodococcus54,92038.37316,6369.83341022.84375,65817.04
Pseudonocardia †12,9139.02733,64019.87825,54517.68972,09815.59
Unclassified57624.031852,90831.262086686.002667,33813.832
Streptomyces †76285.331012920.761143473.011813,2673.019
Arthrobacter43233.02432141.90426491.83510,1862.35
Sporichthya *7710.5415560.33225151.74238420.92
Planotetraspora *27691.9312940.1711130.08131760.71
Nocardia10230.7121570.09212120.84223920.52
Mycobacterium †12990.9121270.0824400.30318660.43
Frankia *10670.752770.0511260.09312700.33
Luedemannella *5550.3922120.1322760.19210430.22
Longispora *3710.2613740.2211100.0818550.21
Agromyces *3100.2221360.0823550.2528010.22
Actinoplanes *5910.412850.0521110.0827870.22
Nakamurella *4160.2911140.0711010.0716310.11
Nocardioides †3600.254520.0331580.1185700.19
Geodermatophilus †780.051700.0423740.2625220.12
Catellatospora *950.073700.0421410.1043060.074
Kribbella †1200.081360.0211350.0922910.072
Kocuria *2610.181230.011---2840.12
Actinomyces *2470.17510.001160.00412540.065
Corynebacterium *1510.111320.023100.0121930.045
Rhizocola *---270.0211450.1011720.061
Microbacterium *1080.081480.03150.00311610.041
Iamia *1070.071---310.0211380.051
Pseudoclavibacter *1380.101------1380.11
Lentzea *------1110.0811110.081
Aeromicrobium †730.051---280.0221010.042
Amycolatopsis860.061---50.0031910.032
Cryptosporangium *---710.041100.011810.021
Glycomyces *350.021---450.031800.031
Streptosporangium *610.041---190.011800.031
Smaragdicoccus *430.031---340.021770.031
Propionibacterium570.041-----0570.041
Kineosporia *440.031---80.011520.021
Jatrophihabitans *---210.011240.021450.011
Promicromonospora *---220.011210.011430.012
Millisia *---220.01160.0041280.011
Rothia *20.0011130.01170.0051220.0051
Tessaracoccus *---170.011---170.011
Acidothermus *------160.011160.011
Marmoricola *------140.012140.012
Dermacoccus *------110.0081110.011
Ponticoccus *---80.0051---80.0051
Stackebrandtia *------80.006180.011
Umezawaea *------20.001120.0011
Actinospica *------10.001110.0011
Propionimicrobium *------10.001110.0011
For each taxon, the number of obtained sequences (Seq.) and their relative abundance (%), together with the number of OTUs, are given. The total number of sequences, average relative abundance, and total number of different OTUs obtained per genus are shown in the last three columns. Taxa marked with an asterisk (*) were reported for the first time in moonmilk deposits in this studyTaxa marked with a cross (†) were detected in moonmilk deposits in this work, and in the high-throughput sequencing (HTS)-based study of Dhami et al. [19]. Taxa underlined represent the ones that were also detected in other moonmilk microbial diversity studies [12,13,14,15,18]. Cases filled in grey highlight the most abundant genera in each studied sampling point. Abbreviations: Seq.—number of sequences identified.
Table
Table 4. The relative abundance (%) and taxonomy assignment of the most abundant actinobacterial OTUs found across moonmilk samples within the “Grotte des Collemboles”.
Table 4. The relative abundance (%) and taxonomy assignment of the most abundant actinobacterial OTUs found across moonmilk samples within the “Grotte des Collemboles”.
OTUCOL1COL3COL4Av. %ClassFamilyGenus
OTU138.289.682.7416.90ActinobacteriaNocardiaceaeRhodococcus
OTU20.3128.892.1310.44ActinobacteriaPseudonocardiaceaeunclassified
OTU80.7513.960.585.10ActinobacteriaPseudonocardiaceaeuncultured
OTU43.351.1311.475.32Acidimicrobiiauncultureduncultured
OTU30.527.804.874.40ActinobacteriaPseudonocardiaceaePseudonocardia
OTU2623.454.125.894.49ActinobacteriaPseudonocardiaceaePseudonocardia
OTU67.531.003.984.17Acidimicrobiiauncultureduncultured
OTU122.976.611.823.80ActinobacteriaPseudonocardiaceaeuncultured
OTU53.531.705.873.70Acidimicrobiiauncultureduncultured
OTU130.464.913.653.01ActinobacteriaPseudonocardiaceaePseudonocardia
OTU980.991.076.943.00Acidimicrobiiauncultureduncultured
OTU2031.730.315.922.66Acidimicrobiiauncultureduncultured
OTU4323.681.401.662.25ActinobacteriaPseudonocardiaceaePseudonocardia
OTU1420.721.912.351.66ActinobacteriaPseudonocardiaceaeunclassified
OTU701.703.311.67ActinobacteriaPseudonocardiaceaeuncultured
OTU192.060.521.651.41ActinobacteriaMicrococcaceaeArthrobacter
OTU1900.460.312.841.20Acidimicrobiiauncultureduncultured
OTU100.810.291.680.93AcidimicrobiiaAcidimicrobiaceaeuncultured
OTU90.240.152.390.93Acidimicrobiiauncultureduncultured
OTU2510.770.651.120.85ActinobacteriaPseudonocardiaceaePseudonocardia
OTU140.540.311.730.86ActinobacteriaSporichthyaceaeSporichthya
OTU302.110.130.310.85ActinobacteriaStreptomycetaceaeStreptomyces
OTU3600.400.171.680.75Acidimicrobiiauncultureduncultured
OTU241.930.170.080.73ActinobacteriaStreptosporangiaceaePlanotetraspora
OTU110.380.031.720.71Acidimicrobiiauncultureduncultured
OTU200.710.470.810.67AcidimicrobiiaIamiaceaeuncultured
OTU470.750.260.980.66Acidimicrobiiauncultureduncultured
OTU150.300.171.480.65ActinobacteriaStreptomycetaceaeStreptomyces
OTU1920.011.650.010.56ActinobacteriaPseudonocardiaceaeuncultured
OTU160.240.511.110.62Acidimicrobiiauncultureduncultured
OTU500.880.540.370.60ActinobacteriaPseudonocardiaceaeuncultured
OTU221.610.160.060.61ActinobacteriaPropionibacteriaceaeunclassified
OTU230.271.190.100.52ActinobacteriaMicrococcaceaeArthrobacter
OTU540.130.920.460.50ActinobacteriaPseudonocardiaceaePseudonocardia
OTU180.480.121.040.54Acidimicrobiiauncultureduncultured
OTU250.560.090.830.49ActinobacteriaNocardiaceaeNocardia
OTU210.800.190.340.45ActinobacteriaStreptomycetaceaeStreptomyces
OTU991.050.090.170.44ActinobacteriaStreptomycetaceaeStreptomyces
OTU360.220.040.960.41Acidimicrobiiauncultureduncultured
OTU440.800.140.260.40Acidimicrobiiauncultureduncultured
OTU320.220.030.920.39Actinobacteriaunclassifiedunclassified

Share and Cite

MDPI and ACS Style

Maciejewska, M.; Całusińska, M.; Cornet, L.; Adam, D.; Pessi, I.S.; Malchair, S.; Delfosse, P.; Baurain, D.; Barton, H.A.; Carnol, M.; et al. High-Throughput Sequencing Analysis of the Actinobacterial Spatial Diversity in Moonmilk Deposits. Antibiotics 2018, 7, 27. https://doi.org/10.3390/antibiotics7020027

AMA Style

Maciejewska M, Całusińska M, Cornet L, Adam D, Pessi IS, Malchair S, Delfosse P, Baurain D, Barton HA, Carnol M, et al. High-Throughput Sequencing Analysis of the Actinobacterial Spatial Diversity in Moonmilk Deposits. Antibiotics. 2018; 7(2):27. https://doi.org/10.3390/antibiotics7020027

Chicago/Turabian Style

Maciejewska, Marta, Magdalena Całusińska, Luc Cornet, Delphine Adam, Igor S. Pessi, Sandrine Malchair, Philippe Delfosse, Denis Baurain, Hazel A. Barton, Monique Carnol, and et al. 2018. "High-Throughput Sequencing Analysis of the Actinobacterial Spatial Diversity in Moonmilk Deposits" Antibiotics 7, no. 2: 27. https://doi.org/10.3390/antibiotics7020027

APA Style

Maciejewska, M., Całusińska, M., Cornet, L., Adam, D., Pessi, I. S., Malchair, S., Delfosse, P., Baurain, D., Barton, H. A., Carnol, M., & Rigali, S. (2018). High-Throughput Sequencing Analysis of the Actinobacterial Spatial Diversity in Moonmilk Deposits. Antibiotics, 7(2), 27. https://doi.org/10.3390/antibiotics7020027

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop