- Open Access
Bioprospecting the microbiome of Red Sea Atlantis II brine pool for peptidases and biosynthetic genes with promising antibacterial activity
Microbial Cell Factories volume 21, Article number: 109 (2022)
The search for novel antimicrobial agents is crucial as antibiotic-resistant pathogens continue to emerge, rendering the available antibiotics no longer effective. Likewise, new anti-cancer drugs are needed to combat the emergence of multi-drug resistant tumors. Marine environments are wealthy sources for natural products. Additionally, extreme marine environments are interesting niches to search for bioactive natural compounds. In the current study, a fosmid library of metagenomic DNA isolated from Atlantis II Deep Lower Convective Layer (ATII LCL), was functionally screened for antibacterial activity as well as anticancer effects.
Two clones exhibited antibacterial effects against the marine Bacillus Cc6 strain, namely clones 102-5A and 88-1G and they were further tested against eleven other challenging strains, including six safe relatives of ESKAPE pathogens (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter spp.), a safe relative to Mycobacterium tuberculosis and four resistant clinical isolates. Clone 88-1G resulted in clear zones of inhibition against eight bacterial strains, while clone 102-5A resulted in zones of inhibition against five bacterial strains. The whole cell lysates of clone 88-1G showed 15% inhibition of Mtb ClpP protease -Mycobacterium tuberculosis drug target-, while whole cell lysates of clone 102-5A showed 19% inhibition of Mtb ClpP protease. Whole cell lysates from the selected clones exhibited anticancer effects against MCF-7 breast cancer cells (cell viability at 50% v/v was 46.2% ± 9.9 for 88-1G clone and 38% ± 7 for 102-5A clone), U2OS osteosarcoma cells (cell viability at 50% v/v was 64.6% ± 12.3 for 88-1G clone and 28.3% ± 1.7 for 102-5A clone) and 1BR hTERT human fibroblast cells (cell viability at 50% v/v was 74.4% ± 5.6 for 88-1G clone and 57.6% ± 8.9 for 102-5A clone). Sequencing of 102-5A and 88-1G clones, and further annotation detected putative proteases and putative biosynthetic genes in clones 102-5A and 88-1G, respectively.
The ATII LCL metagenome hosts putative peptidases and biosynthetic genes that confer antibiotic and anti-cancer effects. The tested clones exhibited promising antibacterial activities against safe relative strains to ESKAPE pathogens and Mycobacterium tuberculosis. Thus, searching the microbial dark matter of extreme environments is a promising approach to identify new molecules with pharmaceutical potential use.
Antibiotic resistance is a major global health challenge, which necessitates the search for new antimicrobial agents . Additionally, resistance of cancers to chemotherapy is another challenge that remains to be addressed by new anticancer agents . One possible way to address such challenges is to search for natural products that have novel chemistries . Particularly, searching for specialized metabolites which are produced by a subset of organisms, including microbes, is a valuable approach as it resulted into an array of compounds with antimicrobial and anticancer activities, such as doxorubicin, bleomycin, and salinosporamide .
Specialized metabolites are produced by a subset of organisms, and the enzymes producing them are encoded by biosynthetic gene clusters (BGCs) in the genome of the producing organism [4,5,6]. Currently, mining for BGCs is an expanding field, with the aim to search for antimicrobial agents to cater to the antimicrobial resistance problem, especially the ESKAPE pathogens [7, 8]. Mining for BGCs within genomes and metagenomes, followed by elucidating the chemical structure and further testing of the bioactive compound has unraveled a plethora of bioactive specialized metabolites [9, 10]. Mining for BGCs is facilitated by bioinformatics approaches, including the use of homology and machine learning algorithms [11, 12]. One approach is to exploit DNA isolated from new niches and under-explored environments, followed by functional screening for the activity of interest [13,14,15]. Such metagenomic approaches have enabled researchers to screen the microbial dark matter for new BGCs with antibacterial activities [16,17,18].
Among the marine ecological niches that harbor a huge microbial diversity are environments with extreme conditions [3, 19]. The Red Sea brine pools are areas with higher temperature and higher salinity than the rest of the Red Sea water. There are twenty-five distinct brine pools in the Red Sea [20, 21]. The largest Red Sea brine pool, possessing the highest recorded temperature and with depth of 2194 m, characterized by having large amounts of heavy metals and very low Oxygen levels, is Atlantis II Deep Red Sea brine pool (ATII) [20, 22]. The Lower Convective Layer (LCL) of the brine pool water is characteristic of multiple extreme conditions, as it has a salinity of 270 psu, a pH of 5.3, in addition to the high temperature of 68.2°C [20, 22]. In this work, we used a previously generated fosmid library from the microbial metagenome of ATII LCL [23, 24] and tested two selected clones for antibacterial and anticancer activity, in addition to sequencing and annotating the genes inserted in these two clones.
Earlier, we have screened the metagenomic fosmid library that harbored prokaryotic environmental DNA from the Atlantis II Red Sea brine pool LCL water samples (ATII) [23, 24]. In our earlier study, we focused on two clones that harbored putative orphan biosynthetic gene clusters . We focused this current study on two other clones, namely 88-1G and 102-5A, for antibacterial effects against clinically relevant strains and assessed their cytotoxic activities against different mammalian cell lines.
The aim of this work was to screen the ATII LCL metagenomic fosmid library by a functional assay for anti-bacterial effects. Two positive clones were then selected to assess their antibacterial effects against clinically relevant bacterial strains. Furthermore, cytotoxic effects of the extracts from the two clones were determined against breast cancer and bone cancer cells in addition to non-cancerous fibroblasts. Lastly, the selected clone inserts were sequenced, and their putative protein-coding genes were annotated. The present two clones harbored individual genes which resulted in the observed effects.
Antibacterial effects of clones 88-1G and 102-5A from Red Sea Atlantis II LCL (or Red Sea ATII LCL) fosmid library against clinically relevant strains
In antimicrobial overlay assays, clones 88-1G and 102-5A exhibited zones of inhibition against the challenging strain Bacillus Cc6. The antibacterial overlay assay was also conducted against eleven additional challenging strains (Table 1), including Staphylococcus epidermidis, Erwinia carotovora, Enterococcus raffinosus, Acinetobacter baylyi, Enterobacter aerogenes, Mycobacterium smegmatis, Pseudomonas putida, and the following clinical isolates: MRSA-ZC1, MRSA-ZC2, MRSA-ZC6, and Acinetobacter ZC2. Clone 88-1G showed zones of inhibition against nine strains: Staphylococcus epidermidis, Erwinia carotovora, Enterococcus raffinosus, Enterobacter aerogenes, Mycobacterium smegmatis, Acinetobacter ZC2, MRSA-ZC1, MRSA-ZC2, and Bacillus Cc6 (Additional file 1: Fig. S1). Clone 102-5A showed zones of inhibition against six of the tested strains, namely Staphylococcus epidermidis, Mycobacterium smegmatis, Pseudomonas putida, Acinetobacter ZC2, MRSA-ZC6, and Bacillus Cc6 (Additional file 1: Fig. S2).
Whole cell lysates from 88-1G and 102-5A Red Sea Atlantis II LCL (or Red Sea ATII LCL) fosmid library clones showed inhibition of M. tuberculosis target ClpP Protease
Using the ClpP protease assay to assess the effectiveness of the clone extracts against Mycobacterium tuberculosis, the 88-1G extract showed ~ 15% inhibition of Mtb ClpP1P2 peptidase activity. This inhibition was statistically significant when compared to the negative controls: TRIS (the extraction solvent) and EPI300 control. The 102-5A fosmid extract was slightly more potent and demonstrated an inhibition activity equivalent to 19% that was also statistically significant when compared to the negative controls (Fig. 1).
Differential cytotoxicity of the extracts from 88-1G and 102-5A Red Sea Atlantis II LCL (or Red Sea ATII LCL) fosmid library clones against cancer cell lines
All the extracts were prepared similarly, and the cytotoxic effects were determined for equal concentrations of whole-cell lysates against MCF-7, U2OS cancer cell lines, and 1BR hTERT cells for 48 h. (Fig. 2). Extracts of both 88-1G and 102-5A clones exhibited significant cancer cell cytotoxic effects at 50% v/v concentration. Whole-cell lysate extracts at 50% v/v of clone 102-5A were more cytotoxic to all the cell lines when compared to 88-1G extracts (Fig. 2). For the breast cancer cell line MCF-7 cells, the cell viability was 46.2% ± 9.9 and 38% ± 7 after 48 h exposure to 88-1G and 102-5A clone whole cell lysate extracts, respectively (Fig. 2A). For the osteosarcoma cell line U2OS cells, the cell viability was 64.6% ± 12.3 and 28.3% ± 1.7 after 48 h exposure to 88-1G clone and 102-5A whole cell lysate extracts, respectively (Fig. 2B). Lastly, for the human skin fibroblasts 1BR hTERT cells, the cell viability was 74.4% ± 5.6 and 57.6% ± 8.9 after 48 h exposure to 88-1G and 102-5A whole cell lysate extracts, respectively (Fig. 2C).
Annotated sequences from 88-1G and 102-5A Red Sea Atlantis II LCL (or Red Sea ATII LCL) fosmid library clones are similar to genes from organisms inhabiting similar extreme niches
The fosmid DNA was extracted from both clones (102-5A and 88-1G) and sequenced in an Illumina MiSeq platform. After trimming the vector and E. coli sequences, 102-5A had 8 scaffolds and 88-1G had 19 scaffolds. The assembly metrics are denoted in Table 2.
After the vector and host sequences were trimmed, the number of PEGs (protein-encoding genes) detected by RAST (Rapid Annotations using Subsystems Technology) in each assembly was 21 and 22 in inserts 102-5A and 88-1G, respectively (Table 2). The majority of PEGs detected by RAST coded for hypothetical proteins. This includes 81% and 86% in 102-5A PEGs and 88-1G, respectively (Additional file 1: Tables S1, 2). To further annotate the PEGs, we analyzed each PEG by PSI-BLAST (Position-Specific Iterative Basic Local Alignment Search Tool) to determine the most similar sequences in publicly available databases. The full annotation results for all the detected PEGs are presented in Additional file 1: Tables S1 and 2. PEGs with hits in PSI-BLAST of higher than the threshold E-value 0.005 were depicted as non-significant (Table 3, Additional file 1: Tables S1, Tables S2).
Annotation results of clones 102-5A and 88-1G showed similarity to organisms with significant hits, which are inhabiting similar environments as the Red Sea Atlantis II brine pool. Protein-based phylogeny for clone 102-5A resulted in similarity to organisms such as Candidatus Bathyarchaeota archaeon ex4484_40, an uncultured hydrothermal sediment archaeon from the Guaymas basin in California , candidate division MSBL1 archaeon SCGC-AAA382A20, uncultured archaea from the Red Sea Neurus brine , Thermoanaerobacterium sp. PSU-2, a thermophilic bacterium , and Candidatus Nanosalinarum sp. J07AB56, a halophilic archaeon . MIBiG (Minimum Information about a Biosynthetic Gene cluster) alignments yielded one significant result, which was evident by the alignment of a putative S8 peptidase with the SLI-2138 RiPP (ribosomally synthesized and post-translationally modified peptide) biosynthetic gene cluster [29, 30]. Similarly, clone 88-1G protein-based phylogeny detected significant similarity with: Candidatus Micrarchaeota archaeon CG1_02_49_24, an uncultured archaeon adapted to high carbon dioxide levels , and the haloalkaliphilic archaea Natronococcus occultus  and Natronobacterium texcoconense . MIBiG comparison yielded two significant results, which was evident by the alignment of a putative serine protease AprX with the SLI-2138 biosynthetic gene cluster, and the alignment of a putative glycosyltransferase with a glucosyltransferase within a carotenoid biosynthetic gene cluster .
Annotation of putative biosynthetic genes inserted in clones 102-5A and 88-1G
Clone 88-1G insert harbored the following PEGs (Fig. 3, Table 3): (1) a putative LysR family transcriptional regulator [Hydrocarboniphaga effusa], (2) a putative CopG family transcriptional regulator [Natronococcus occultus], (3) a putative serine protease AprX [Bradyrhizobium erythrophlei] and (4) a putative glycosyltransferase [Candidatus Micrarchaeota archaeon CG1_02_49_24].
Annotation of clone 102-5A insert revealed two putative peptidases and a putative hydrolase, with the best hits being (Fig. 4, Table 3): (1) S9 family peptidase [Nocardia pseudobrasiliensis], (2) peptidase S8 [Streptomyces sp. CB02613] and (3) cobalt-precorrin 5A hydrolase [Hyphomicrobium facile].
Broad spectrum antibacterial activity of clones from Red Sea brine pool Atlantis II (ATII) LCL fosmid library
Clone 88-1G exhibited antibacterial activity against Bacillus Cc6 marine strain, which is a Gram-positive marine bacterial strain (Table 1, Additional file 1: Fig. S1), implying its possible in situ activity, as the metagenomic clones were from a marine environment . Additionally, antibacterial activity was also observed against clinically relevant strains including Staphylococcus epidermidis, Erwinia carotovora, Enterococcus raffinosus and Enterobacter aerogenes, which are safe relatives of ESKAPE pathogens (Table 1, Additional file 1: Fig. S1). Methicillin-resistant Staphylococcus aureus (MRSA) and Acinetobacter (Acinetobacter AB2) clinical isolates were obtained from hospitalized patients in Egypt. Clone 88-1G inhibited the growth of Acinetobacter AB2 strain, MRSA-ZC1, and MRSA-ZC2. Additionally, clone 88-1G resulted in an inhibition zone with the acid-fast bacteria Mycobacterium smegmatis. This result was corroborated with the results of the Mtb ClpP protease assay [30, 31], as the whole-cell lysate of 88-1G resulted in 15% inhibition of the M. tuberculosis target enzyme.
Clone 102-5A similarly exhibited antibacterial activity against Bacillus Cc6 marine strain (Table 1, Additional file 1: Fig. S2) [24, 34], however it showed a broad but different pattern of zones of inhibition than the one observed for 88-1G clone. Antibacterial activity was also observed against Staphylococcus epidermidis and Pseudomonas putida, other safe relatives of ESKAPE pathogens (Table 1, Additional file 1: Fig. S2). Antimicrobial activities against the ESKAPE pathogens are important because of their emerging antimicrobial resistance, particularly in nosocomial settings [32, 33]. Clone 102-5A inhibited the growth of Acinetobacter AB2 strain and MRSA-ZC6. Additionally, it displayed an inhibition zone against Mycobacterium smegmatis (Table 1, Additional file 1: Fig. S2). We confirmed the former results by performing an Mtb ClpP protease assay [35, 36]. The whole-cell lysates of 102-5A resulted in 19% inhibition of the M. tuberculosis target enzyme, which is higher inhibition than that observed for lysates of 88-1G clone.
In conclusion, each of the tested clones exhibited a unique and broad antibacterial activity pattern, including gram-positive, gram-negative strains, and acid-fast bacteria.
Anticancer activity of whole-cell lysates from the Red Sea brine pool Atlantis II (ATII) LCL fosmid library clones
We observed differential cytotoxicity of the two fosmid library clones against the different tested cell lines. For the breast cancer cell line MCF-7, the whole cell lysate extracts of clone 102-5A exhibited a greater cytotoxic effect with cell viability (38% ± 7), when compared to the 88-1G clone, which exhibited a cell viability of 46.2% ± 9.9 (Fig. 2A). 102-5A extracts also exhibited more cytotoxicity against U2OS cells -osteosarcoma cell line- (28.3% ± 1.7) when compared to extracts from the 88-1G (64.6% ± 12.3) (Fig. 2B).
The cell viability of 1BR hTERT cells incubated with clone 88-1G extracts was close to that of the negative control: 74.4% ± 5.6 cell viability and 71.6% ± 5.6 cell viability with 88-1G extracts and buffer control, respectively. On the other hand, 102-5A cell lysates were more lethal to the 1BR hTERT cells (57.6% ± 8.9 viability) (Fig. 2C). The lack of cytotoxicity observed on incubation of 1BR hTERT with 88-1G suggests that 88-1G is a better candidate for selective cancer cytotoxicity.
Putative peptidases and hydrolase detected on clone 102-5A insert of the Red Sea brine pool Atlantis II (ATII) LCL
Clone 102-5A annotation included putative (Fig. 4, Table 3): (1) S9 peptidase, (2) S8 peptidase and (3) a hydrolase. This includes a serine protease of the S9 peptidase family, whose members are not yet fully characterized . A novel putative hydrolase was identified by a similar approach from the metagenomes of microbes associated with marine sponges and was observed to exhibit antibacterial effects . It was hypothesized that hydrolases were produced by eukaryotic hosts to limit the overgrowth of the microbial symbionts . Although the microbial community included in the present study is not a symbiotic community owing to the studied sample, however, it is likely that the detected putative hydrolytic enzymes would have similar antagonistic effects owing to the observed antibacterial activities.
Putative proteases with antibacterial effects were detected in a study using Ralstonia species as the host for heterologous gene expression, rather than E. coli . The detected putative proteases belonged to the metallopeptidase family and the S8A serine protease family . Such antibacterial properties of proteases are likely to act on breaking the amide bond connecting two adjacent present peptidoglycan layers within the cell walls of both gram-positive and gram-negative bacteria, which could explain the broad spectrum of activity observed in this study . Previous studies have reported the antibacterial and anticancer activities of peptidases [15, 38, 39]. It is worth noting that S8 peptidases (subtilases) have industrial significance including the high-alkaline subtilases (family A), which are used in the detergent industry .
The annotated peptidases and hydrolase may play a role in the observed antibacterial and cytotoxic effects, particularly that two putative serine peptidases were detected on 102-5A clone insert; namely an S9 and an S8 peptidase, in addition to a putative hydrolase. Further structure and function analysis studies need to be performed to confirm the hypothesized roles of the annotated peptidases and hydrolase in the antibacterial and cytotoxic effects and to pinpoint the inhibitory role of each protein.
Putative biosynthetic gene and peptidase detected on clone 88-1G insert of the Red Sea brine pool Atlantis II (ATII) LCL
Clone 88-1G fosmid DNA harbored several putative biosynthetic genes linked to the observed activities (Fig. 3, Table 3). Four genes were detected that are likely to contribute to the observed activities: (1) a serine protease, (2) a glycosyltransferase and (3) two putative transcriptional regulators. A serine protease was detected, and as aforementioned, it is likely to function by targeting the synthesis of the peptidoglycan bacterial cell wall .
Furthermore, glycosyltransferases have been reported to play a pivotal role in specialized metabolite biosynthetic pathways [11, 41]. Glycosyltransferase enzymes catalyze the formation of glycosidic bonds and the transfer of sugar moieties to carbohydrate, lipid, small molecule, protein or DNA substrates . Previous studies had identified the role of glycosyltransferases in the synthesis of natural products . Further biochemical studies need to be conducted to elucidate the natural product produced by the putative glycosyltransferase enzyme and its role in the observed antibacterial and/or cytotoxic effects.
Our present study points to the importance of bioprospecting microbial dark matter from extreme environments, particularly in the discovery of antibacterial and anti-cancer agents. Yet the present study use of a mesophilic bacterial strain in the fosmid library construction is likely to limit the expression of extremophilic genes  and therefore further studies should assess the activity in extremophiles.
Additionally, future studies should aim to identify the nature of the antibacterial and/or anticancer active compounds, whether it is a specialized metabolite, or a protein as per the annotations. Peptide active compounds can be challenging in drug delivery . However, several studies are considering AMPs (antimicrobial peptides), RiPPs, bacteriocins and other peptides for anti-bacterial and anti-cancer testing, to increase the array of pharmaceutically-relevant natural products as drug leads [9, 44]. Further structure–function studies would also be interesting in elucidating the cytotoxic and antibacterial cellular mechanisms of the active compound/s in the clones.
Mining microbial metagenomes for specialized metabolites with pharmaceutical use is one way to find new antibacterial and anticancer compounds. Through the approach of mining yet under-explored environments of microbiomes inhabiting harsh marine niches, followed by functional screening, putative peptidases and biosynthetic genes have been identified and further studies are required to fully decipher the specialized metabolites’ structures.
Metagenomic fosmid library screening for antibacterial activity
Water samples from the lower Convective Layer (LCL) of Atlantis II Red Sea brine pool (21° 20.72’ N and 38° 04.59’ E) was previously collected in the 2010 KAUST/WHOI/HCMR expedition . Prokaryotic eDNA was extracted from the 0.1 µm filter as previously described . The ATII LCL fosmid library was previously constructed as per the manufacturer’s protocol and it was composed of 10,656 clones, by using pCC2FOS vector with the CopyControl Fosmid Library Production Kit (Epicenter) . A fresh copy of the aforementioned fosmid library was prepared prior to the downstream assays and further used . The recipient strain for the fosmid was EPI300-T1R E. coli as per the manufacturer’s protocol.
Antimicrobial overlay assay
An antimicrobial overlay assay was used to test for clones with anti-bacterial activity similar to that reported in the literature [34, 46]. E. coli CBAA11 strain producing the antibacterial tambjamine was used as the positive control strain [13, 34]. The challenging strain that was used was Bacillus Cc6 strain which is a marine Bacillus strain associated with the Australian marine sponge Cymbastela concentrica (gift from Prof. Torsten Thomas, University of New South Wales) [13, 34]. Bacillus Cc6 strain was used as a challenging strain because it was isolated from a marine environment and because it was reported in a functional screening assay utilizing the same vector CC2 as used in our study [13, 34]. E. coli CBAA11 strain was chosen as a positive control because it produces the antibacterial tambjamine and was reported to result in clear zones against Bacillus Cc6 [13, 34]. All the plates were supplemented with 0.01% arabinose and 12.5 µg/ml chloramphenicol. All the plates were tested against the clones and included both a positive control (E.coli CBAA11) and a negative control (EPI300).
First the clones were plated on LB/chloramphenicol/arabinose overnight at 37 °C. Chloramphenicol was not added to the other strains, other than Bacillus Cc6 strain. The following day, the plates were incubated overnight at 25 °C and an overnight culture of Bacillus CC6 in LB/chloramphenicol was prepared and incubated at 37 °C with shaking. The following day, the starter culture was used to inoculate 100 ml LB/chloramphenicol of Bacillus CC6 till the OD reached 0.5. Then Bacillus CC6 was poured on the plates using top agar at 1:1000 dilution. The overlaid plates were finally incubated overnight at 25 °C and observed carefully for clear zones in the top layer . Clones exhibiting antibacterial activity were further tested against the eleven challenging strains (Table 1).
Whole-cell extracts were prepared from the clones namely 102-5A and 88-1G, using a method similar to that reported in the literature . Overnight cultures from the E. coli bearing the eDNA from the positive clones -previously supplemented with auto-induction solution and chloramphenicol- were centrifuged at 3,500 rpm for 10 min. Afterwards, the cell pellets were resuspended in 0.01 M Tris–HCl pH 7. The extracts were then sonicated on ice at 20% maximal amplitude for 370 s (with 10 s intervals without sonication). Finally, the extracts were sterile-filtered with 0.2 µm membrane filters (Corning) . Two negative controls were used: (1) cells incubated with only the media, and (2) cells incubated with buffer. The buffer concentration, in the latter, was the same concentration used in the highest extract concentration 50% v/v (Fig. 2).
Cell lines and culture conditions
Three cell lines were used for the cytotoxicity assay, a human breast adenocarcinoma cell line (MCF-7) , an osteosarcoma cell line (U2OS)  (gift from Dr. Andreas Kakarougkas (University of Sussex) and skin fibroblast cells–wild-type and non-cancerous cells- that are immortalized with human telomerase reverse transcriptase (1BR hTERT) [50,51,52]. The cells were cultured in DMEM (Lonza, Germany), supplemented with 10% Fetal Bovine Serum (Lonza, Germany) and 5% Penicillin–Streptomycin (Lonza, Germany). All cells were grown in an incubator adjusted at 37°C and 5% CO2 at a seeding density of 104 cells/well.
Cell viability assay
The initial seeding density was adjusted to 104 cells/well and left overnight to adhere to the bottom of the 96-well plates (Greiner Bio-One, Germany). The old media was discarded the following day and the media containing the different v/v% concentrations of the extracts was added (0, 15, 20 and 50%). The cells were exposed for 48 h to the extracts before assessing the cell viability.
The extracts were incubated for 48 h with MCF-7, U2OS and 1BR hTERT cell lines with different v/v% concentrations (0, 15, 20 and 50%) . Cell viability percentage was assessed by the MTT assay after 48 h exposure to the extracts. Firstly, the media was replaced by 100 µl fresh media containing MTT reagent 3-(4, 5-dimethylthiazolyl-2)-2, 5-diphenyltetrazolium bromide in a concentration of 20 µl MTT (5 mg/ml) (Serva, Germany) for 3 h. Afterwards, the media was discarded and 100 µl DMSO (Sigma-Aldrich, USA) was added to solubilize the purple precipitates.
Several negative controls were used in the experiments. Cells with no cell lysates added were used as the first negative control cells. Cells exposed to the buffer only with no cell lysates were used to test the buffer effect. Additionally, cells exposed to the cell lysates of EPI300 strain were used to test the effect of the strain.
All the negative control cells (A595 control) were supplemented with complete media. The wells used for measuring A595 blank were only supplemented with media and no cells. Absorbance was finally measured at 595 nm (A595), using a SPECTROstar Nano microplate reader (BMG LabTech, Germany). Cell viability was calculated as follows:
Sequencing and bioinformatics
The two clones, 88-1G and 102-5A, were selected for DNA extraction and sequencing . Overnight cultures were supplemented with auto-inducer/chloramphenicol. Fosmid DNA was extracted by QIAprep Spin Miniprep Kit (Qiagen). The fosmid DNA was sequenced using Illumina MiSeq V3 300 bp paired-end read platform (LGC, Germany). After sequencing and quality filtering, the reads were assembled by de novo assembly programs, namely SOAPdenovo2  and CLC Genomics Workbench v 8.0 (Qiagen), respectively. The details of the sequencing and assembly metrics are depicted in Table 2.
Prior to annotation, the vector sequences—pCC2FOS™- were trimmed from the resulting scaffolds. E. coli sequences were also removed from the sequences. NC_010473 DH10B reference sequence was used because the strain used for the fosmid library construction—EPI300™-T1R Phage T1-resistant E. coli- was derived from E. coli DH10B strain. The resulting scaffolds were used to determine the putative PEGs –protein-encoding genes- using RAST platform (analysis conducted March 2017) . Each PEG was further compared to sequences in the publicly available databases by using PSI-BLAST . The PEGs were also compared to the protein sequences curated in the Minimum Information about a Biosynthetic Gene cluster (MIBiG) database by BLASTX . Lastly, PFAM alignments were done for the selected PEGs (Additional file 1: Fig. S4).
M. tuberculosis ClpP Protease Assay
The protein degradation enzyme ClpP protease of M. tuberculosis that is composed of subunits P1 and P2 is a proven drug target against M. tuberculosis. Mtb ClpP assay was used to determine the inhibition of the peptidase activity of Mtb ClpP1P2. The cleavage of the labeled peptide substrate generates fluorescence at 460 nm (excitation at 355 nm). The assay was performed in a F-bottom, black plate, clear bottom 96- well plate with a final reaction volume of 30 µl, as per the manufacturer's instructions (ProFoldin, USA). The assay buffer and Mtb ClP1P2 were mixed. For the negative control, DMSO was used, and for the positive control, the Mtb ClP1P2 was boiled at 95 °C for 10 min (35, 36). Finally, for the extracts under investigation, 0.6 µl of 450 µg/mL of 88-1G, and 102-5A were added and incubated at 37 °C for 20 min. The 100 X labeled substrate was freshly diluted to make the 10 X, added to the assay reaction solution, and set again at 37°C for 60 min. The fluorescence was measured using FLUOstar OPTIMA microplate reader (BMG LabTech, Germany) by exciting at 355 nm and following the increase in fluorescence emission at 460 nm. Percentage inhibition was calculated = F (sample)/F (control) × 100.
Statistical analysis was conducted by one-way ANOVA and Post-hoc Tukey test to calculate the P values for the cytotoxicity results (Fig. 2). R program version 4.1.2 (R Development Core Team 2016) was used to compute the statistical calculations. The P values represent the significant differences between the mean of each condition and the mean of the negative control cells with buffer concentration 50% v/v (* ≤ 0.05, ** ≤ 0.01, *** ≤ 0.001 and **** ≤ 0.0001).
R program was also used to generate Fig. 1, and Kruskal–Wallis one-way analysis of variance (ANOVA) was used to compute the p-values.
Availability of data and materials
The datasets of generated and analysed during the current study, comprising all scaffold sequences and their corresponding RAST annotations, are available with GenBank Accession Numbers BankIt2431428 OK545880-OK545898 for 88-1G, and BankIt2509949 OK545899-OK545906 for 102-5A.
Martens E, Demain AL. The antibiotic resistance crisis, with a focus on the United States. J Antibiotics. 2017;70:520–6.
Galluzzi L, Senovilla L, Vitale I, Michels J, Martins I, Kepp O, et al. Molecular mechanisms of cisplatin resistance. Oncogene. 2012;31:15.
Newman DJ, Cragg GM. Natural products as sources of new drugs over the nearly four decades from 01/1981 to 09/2019. J Nat Prod. 2020;83(3):770–803. https://doi.org/10.1021/acs.jnatprod.9b01285.
Davies J, Ryan KS. Introducing the Parvome: bioactive compounds in the microbial world. ACS Chem Biol. 2011;7(2):252–9. https://doi.org/10.1021/cb200337h.
Davies J. Specialized microbial metabolites: functions and origins. J Antibiotics. 2013;66(7):361–4.
Jensen PR. Natural products and the gene cluster revolution. Trends Microbiol. 2016;24(12):968. /pmc/articles/PMC5123934/ (Accessed 6 Feb 2022)
Long S, Sousa E, Kijjoa A, Pinto MMM. Marine natural products as models to circumvent multidrug resistance. Molecules. 2016;21(7):892.
Caesar LK, Montaser R, Keller NP, Kelleher NL. Metabolomics and genomics in natural products research: complementary tools for targeting new chemical entities. Nat Prod Rep. 2021;38(11):2041–65.
Bauman KD, Butler KS, Moore BS, Chekan JR. Genome mining methods to discover bioactive natural products. Nat Prod Rep. 2021;38(11):2100–29.
Chevrette MG, Gavrilidou A, Mantri S, Selem-Mojica N, Ziemert N, Barona-Gómez F. The confluence of big data and evolutionary genome mining for the discovery of natural products. Nat Product Reports. 2021;38(11):2024–40.
Blin K, Shaw S, Steinke K, Villebro R, Ziemert N, Lee SY, et al. antiSMASH 50: updates to the secondary metabolite genome mining pipeline. Nucleic Acids Res. 2019;47(W1):W81–7.
Hannigan GD, Prihoda D, Palicka A, Soukup J, Klempir O, Rampula L, et al. A deep learning genome-mining strategy for biosynthetic gene cluster prediction. Nucleic Acids Res. 2019; 47(18):e110. /pmc/articles/PMC6765103/
Burke C, Thomas T, Egan S, Kjelleberg S. The use of functional genomics for the identification of a gene cluster encoding for the biosynthesis of an antifungal tambjamine in the marine bacterium Pseudoalteromonas tunicata. Environ Microbiol. 2007;9(3):814–8. https://doi.org/10.1111/j.1462-2920.2006.01177.x.
Chevrette MG, Handelsman J. Needles in haystacks: reevaluating old paradigms for the discovery of bacterial secondary metabolites. Nat Product Reports. 2021;38(11):2083–99.
Iqbal HA, Craig JW, Brady SF. Antibacterial enzymes from the functional screening of metagenomic libraries hosted in Ralstonia metallidurans. FEMS Microbiol Lett. 2014;354(1):19–26.
Donia MS, Ruffner DE, Cao S, Schmidt EW. Accessing the hidden majority of marine natural products through metagenomics. Chembiochem. 2011; 12(8):1230. /pmc/articles/PMC3397154/ (Accessed 6 Feb 2022)
Wilson MC, Piel J. Metagenomic approaches for exploiting uncultivated bacteria as a resource for novel biosynthetic enzymology. Chem Biol. 2013;20(5):636–47.
Robinson SL, Piel J, Sunagawa S. A roadmap for metagenomic enzyme discovery. Nat Product Reports. 2021;38(11):1994–2023.
Milshteyn A, Schneider JS, Brady SF. Mining the metabiome: identifying novel natural products from microbial communities. Chem Biol. 2014;21(9):1211–23.
Qian PY, Wang Y, Lee OO, Lau SCK, Yang J, Lafi FF, et al. Vertical stratification of microbial communities in the Red Sea revealed by 16S rDNA pyrosequencing. ISME J. 2011;5:3.
Antunes A, Ngugi DK, Stingl U. Microbiology of the Red Sea (and other) deep-sea anoxic brine lakes. Environ Microbiol Rep. 2011;3(4):416–33.
Swift SA, Bower AS, Schmitt RW. Vertical, horizontal, and temporal changes in temperature in the Atlantis II and Discovery hot brine pools, Red Sea. Deep Sea Res Part 1 Oceanogr Res Pap. 2012;64:118–28.
Mohamed YM, Ghazy MA, Sayed A, Ouf A, El-Dorry H, Siam R. Isolation and characterization of a heavy metal-resistant, thermophilic esterase from a Red Sea Brine Pool. Sci Rep. 2013;3(1):1–8.
Ziko L, Saqr AHA, Ouf A, Gimpel M, Aziz RK, Neubauer P, et al. Antibacterial and anticancer activities of orphan biosynthetic gene clusters from Atlantis II Red Sea brine pool. Microb Cell Fact. 2019;18(1):1–16. https://doi.org/10.1186/s12934-019-1103-3.
Dombrowski N, Seitz KW, Teske AP, Baker BJ. Genomic insights into potential interdependencies in microbial hydrocarbon and nutrient cycling in hydrothermal sediments. Microbiome. 2017;5(1):106. https://doi.org/10.1186/s40168-017-0322-2.
Mwirichia R, Alam I, Rashid M, Vinu M, Ba-Alawi W, Anthony Kamau A, et al. Metabolic traits of an uncultured archaeal lineage-MSBL1-from brine pools of the Red Sea. Sci Rep. 2016;6(1):1–14.
O-Thong S, Khongkliang P, Mamimin C, Singkhala A, Prasertsan P, Birkeland NK. Draft genome sequence of Thermoanaerobacterium sp. strain PSU-2 isolated from thermophilic hydrogen producing reactor. Genomics Data. 2017;12:49–51.
Narasingarao P, Podell S, Ugalde JA, Brochier-Armanet C, Emerson JB, Brocks JJ, et al. De novo metagenomic assembly reveals abundant novel major lineage of Archaea in hypersaline microbial communities. ISME J. 2012;6(1):81–93.
Kersten RD, Yang YL, Xu Y, Cimermancic P, Nam SJ, Fenical W, et al. A mass spectrometry-guided genome mining approach for natural product peptidogenomics. Nat Chem Biol . 2011;7(11):794. /pmc/articles/PMC3258187/ (Accessed 6 Feb 2022)
Kautsar SA, Blin K, Shaw S, Navarro-Muñoz JC, Terlouw BR, van der Hooft JJJ, et al. MIBiG 20: a repository for biosynthetic gene clusters of known function. Nucleic Acids Res. 2020;48(D1):D454–8.
Probst AJ, Castelle CJ, Singh A, Brown CT, Anantharaman K, Sharon I, et al. Genomic resolution of a cold subsurface aquifer community provides metabolic insights for novel microbes adapted to high CO2 concentrations. Environ Microbiol. 2017;19(2):459–74. https://doi.org/10.1111/1462-2920.13362.
Studdert CA, de Castro RE, Seitz KH, Sánchez JJ. Detection and preliminary characterization of extracellular proteolytic activities of the haloalkaliphilic archaeon Natronococcus occultus. Arch Microbiol. 1997;168(6):532–5.
Ruiz-Romero E, Sánchez-López KB, Coutiño-Coutiño de los MA, González-Pozos S, Bello-López JM, López-Ramírez MP, et al. Natronobacterium texcoconense sp. nov., a haloalkaliphilic archaeon isolated from soil of a former lake. Int J Syst Evolut Microbiol. 2013;63(11):4163–6. https://doi.org/10.1099/ijs.0.053629-0.
Yung PY, Burke C, Lewis M, Kjelleberg S, Thomas T. Novel antibacterial proteins from the microbial communities associated with the sponge Cymbastela concentrica and the green alga Ulva australis. Appl Environ Microbiol. 2011;77(4):1512–5.
Moreno-Cinos C, Goossens K, Salado IG, van der Veken P, de Winter H, Augustyns K. ClpP Protease, a promising antimicrobial target. Int J Mol Sci. 2019; 20(9). /pmc/articles/PMC6540193/ (Accessed 6 Feb 2022)
Nagpal J, Paxman JJ, Zammit JE, Alhuwaider A, Truscott KN, Heras B, et al. Molecular and structural insights into an asymmetric proteolytic complex (ClpP1P2) from Mycobacterium smegmatis. Sci Rep. 2019;9(1):1–12.
Rawlings ND. Identification and prioritization of novel uncharacterized peptidases for biochemical characterization. Database. 2013; 2013. (Accessed 6 Feb 2022)
Firczuk M, Gabrysiak M, Barankiewicz J, Domagala A, Nowis D, Kujawa M, et al. GRP78-targeting subtilase cytotoxin sensitizes cancer cells to photodynamic therapy. Cell Death Dis. 2013;4(7):e741–e741.
Iqbal HA, Feng Z, Brady SF. Biocatalysts and their small molecule products from metagenomic studies. Curr Opin Chem Biol . 2012; 16(1–2):109. /pmc/articles/PMC3335775/ (Accessed 6 Feb 2022)
Niehaus F, Gabor E, Wieland S, Siegert P, Maurer KH, Eck J. Enzymes for the laundry industries: tapping the vast metagenomic pool of alkaline proteases. Microb Biotechnol. 2011; 4(6): 767. /pmc/articles/PMC3815412/ (Accessed 6 Feb 2022)
Skinnider MA, Dejong CA, Rees PN, Johnston CW, Li H, Webster ALH, et al. Genomes to natural products PRediction Informatics for Secondary Metabolomes (PRISM). Nucleic Acids Res. 2015;43(20):9645–62.
Schmid J, Heider D, Wendel NJ, Sperl N, Sieber V. Bacterial Glycosyltransferases: challenges and opportunities of a highly diverse enzyme class toward tailoring natural products. frontiers in microbiology. 2016;7(FEB). /pmc/articles/PMC4757703/ (Accessed 16 Mar 2022)
Lairson LL, Henrissat B, Davies GJ, Withers SG. Glycosyltransferases: structures, functions, and mechanisms. Annu Rev Biochem. 2008;77:521–55.
Bruno BJ, Miller GD, Lim CS. Basics and recent advances in peptide and protein drug delivery. Ther Deliv. 2013;4(11):1443–67.
Rusch DB, Halpern AL, Sutton G, Heidelberg KB, Williamson S, Yooseph S, et al. The Sorcerer II global ocean sampling expedition: Northwest Atlantic through Eastern Tropical Pacific. PLOS Biol. 2007;5(3):e77. https://doi.org/10.1371/journal.pbio.0050077.
Brady SF. Construction of soil environmental DNA cosmid libraries and screening for clones that produce biologically active small molecules. Nat Protoc. 2007;2(5):1297–305.
Felczykowska A, Dydecka A, Bohdanowicz M, Gasior T, Soboň M, Kobos J, et al. The use of fosmid metagenomic libraries in preliminary screening for various biological activities. Microb Cell Fact. 2014;13(1):1–7. https://doi.org/10.1186/s12934-014-0105-4.
Soule HD, Vazquez J, Long A, Albert S, Brennan M. A human cell line from a pleural effusion derived from a breast carcinoma. J Natl Cancer Inst. 1973;51(5):1409–16.
Heldin CH, Johnsson A, Wennergren S, Wernstedt C, Betsholtz C, Westermark B. A human osteosarcoma cell line secretes a growth factor structurally related to a homodimer of PDGF A-chains. Nature. 1986;319(6053):511–4.
Auclair Y, Rouget R, Affar EB, Drobetsky EA. ATR kinase is required for global genomic nucleotide excision repair exclusively during S phase in human cells. Proceed Nat Academy Sci. 2008;105(46):17896–901.
Brunton H, Goodarzi AA, Noon AT, Shrikhande A, Hansen RS, Jeggo PA, et al. Analysis of human syndromes with disordered chromatin reveals the impact of heterochromatin on the efficacy of ATM-dependent G2/M checkpoint arrest. Mol Cell Biol. 2011;31(19):4022–35.
Fujisawa H, Nakajima NI, Sunada S, Lee Y, Hirakawa H, Yajima H, et al. VE-821, an ATR inhibitor, causes radiosensitization in human tumor cells irradiated with high LET radiation. Radiat Oncol. 2015;10(1):1–7. https://doi.org/10.1186/s13014-015-0464-y.
Luo R, Liu B, Xie Y, Li Z, Huang W, Yuan J, et al. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. Gigascience. 2012;1(1):2047–117.
Aziz RK, Bartels D, Best A, DeJongh M, Disz T, Edwards RA, et al. The RAST Server: rapid annotations using subsystems technology. BMC Genomics. 2008;9(1):1–15. https://doi.org/10.1186/1471-2164-9-75.
Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402.
We are grateful to Prof. Torsten Thomas (University of New South Wales) for providing us with the strains for the overlay assay. We would like to thank Prof. Mohamed El-Hadidy (Zewail City for Science and Technology) for providing us with the clinical isolates. We are also thankful to Dr. Nichole Broderick (John Hopkins University) for the safe ESKAPE relative strains. Additionally, we would like to thank Dr. Andreas Kakarougkas (University of Sussex) for providing us with U2OS cells.
This research was funded by the American University in Cairo and partially funded by the Cayman Biomedical Research Institute (CABRI), USA.
Ethics approval and consent to participate
Consent for publication
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Figure S1
. The anti-bacterial overlay assay results of 88-1G fosmid clone, Figure S2. The anti-bacterial overlay assay results of 102-5A fosmid clone, Figure S3. Negative control for the overlay assay. Figure S4. The alignment of the PFAM hits of the selected PEGs in 88_1G and 102_5A clones. Table S1. Annotation of all PEGs after filtering and trimming of the putative orphan biosynthetic gene clusters on 88-1G clone insert, Table S2. Annotation of all PEGs after filtering and trimming of the putative orphan biosynthetic gene clusters on 102-5A clone insert.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Ziko, L., AbdelRaheem, O., Nabil, M. et al. Bioprospecting the microbiome of Red Sea Atlantis II brine pool for peptidases and biosynthetic genes with promising antibacterial activity. Microb Cell Fact 21, 109 (2022). https://doi.org/10.1186/s12934-022-01835-z
- Natural products
- Red Sea
- Specialized metabolites