Genomic insights into Dactylosporangium sp. AC04546, a rubber degrader with three latex clearing proteins

With more than 100 rubber-degrading strains being reported, only 9 Lcp proteins isolated from Nocardia, Gordonia, Streptomyces, Rhodococcus, Actinoplanes, and Solimonas have been purified and biochemically characterized. A new strain, Dactylosporangium sp. AC04546 (strain JCM34239), isolated from soil samples collected in Sarawak Forest, was able to grow and utilize natural or synthetic rubber as the sole carbon source. Complete genome of Strain AC04546 was obtained from the hybrid assembly of PacBio Sequel II and Illumina MiSeq. Strain AC04546 has a large circular genome of 13.08 Mb with a G+C content of 72.1%. The genome contains 11,865 protein-coding sequences with 3 latex clearing protein (lcp) genes located on its chromosome. The genetic organization of the lcp gene cluster is similar to two other reported rubber-degrading strains—Actinoplanes sp. OR16 and Streptomyces sp. CFMR 7. All 3 Lcp from strain AC04546 were expressed in Escherichia coli and exhibited degrading activity against natural rubber. The distinctiveness of strain AC04546, along with other characterized rubber-degrading strains, is reported here.

1 Introduction

Rubber manufacturing and production remain strong year after year due to high market demand; the amount of natural rubber produced in 2015 doubled compared to 2000 (Leong et al., 2023). Rubber products are resistant to natural degradation due to the vulcanization process and the presence of additives, antioxidants, and fillers, which take thousands of years to degrade naturally. Over the years, rubber products (e.g., tires) have been taking up landfills, releasing toxic chemicals into the environment, becoming breeding grounds for mosquitoes, and increasing the risk of fire hazards. Recent findings reveal that products that either consist of or contain high levels of rubber had the worst toxicity impact on the microorganisms investigated (Sørensen et al., 2023).

Green technology, such as biodegradation, is important to ensure the success of “responsible consumption and production” under Sustainable Development Goal 12. A safer alternative is needed to discard, process, or convert rubber wastes into high–value-added products. Rubber biodegradation refers to the degradation and devulcanization of rubber products by microorganisms growing on the rubber surface. These microorganisms exert their effect by adhering to the rubber surface or secreting enzymes that degrade the surroundings. Unlike conventional methods of burning or landfilling, rubber biodegradation does not produce harmful and toxic side products. Complex rubber molecules are converted into small oxidized derivatives used as carbon and energy sources by microorganisms (Soares and Steinbüchel, 2022). Both Gram-positive and Gram-negative bacteria have rubber-degrading enzymes, either latex-clearing protein (Lcp) or rubber oxygenase A (RoxA) and B (RoxB). Although structurally different, the enzymes are dioxygenases that catalyze oxidative cleavage of the cis double bonds in the poly (cis-1,4-isoprene) backbone, generating carbonyl and aldehyde end groups (Birke et al., 2017). RoxA and RoxB are key enzymes in rubber degradation by Gram-negative bacteria.

Well-studied rubber-degrading bacteria mostly belong to the Gram-positive group Actinobacteria, known to secrete Lcp. The group includes strains of genera such as Streptomyces, Nocardia, Gordonia, Rhodococcus, and Actinoplanes. We previously reported a newly identified rubber-degrading strain, Dactylosporangium sp. AC04546, which was identified from a group of 940 actinobacteria (Basik et al., 2021). Those strains were screened for their ability to produce clearing zones on latex overlay agar, which is an indication of rubber-degrading enzymes being secreted (Basik et al., 2021).

Dactylosporangium is a genus first discovered in 1967 (Thiemann et al., 1967), and only 17 species have been reported to date, classifying it as a rare genus that is less frequently isolated. Strain AC04546 was isolated from a soil sample collected in a secondary forest near a large Aglaia tree in Sarawak, Malaysia. The tree and its sap may have influenced the evolution of this rubber-degrading strain. Here, we report the complete genome of strain AC04546 compared to the genomes of other related rubber degraders. To date, eight incomplete and five complete Dactylosporangium genomes have been submitted to the GenBank database. The complete genome size ranges from 8.221 to 11.96 Mb. Strain AC04546 has the largest genome at 13.08 Mb, and an in-silico comparison classifies it as a distinct species. We are interested in this strain for its ability to utilize natural and vulcanized rubber as the sole carbon source. Strain AC04546 grew well in the presence of tire samples as the sole carbon and energy source, showing visible changes on the tire surface seen in scanning electron microscope (SEM) images and changes in the Attenuated Total Reflection-Fourier Transform Infrared (ATR-FTIR spectra) within 30 days of incubation (Basik et al., 2021).

Apart from strain AC04546, strains with 3 lcp genes on their chromosomes have been reported for Streptomyces sp. CFMR 7 (Nanthini et al., 2017) and Actinoplanes sp. OR16 (NBRC 114529). Actinoplanes sp. OR16 was isolated from a botanical garden in Japan (Imai et al., 2011), and the lcp genes have been functionally characterized (Gibu et al., 2020). Here, we aim to characterize strain AC04546 with reference to other rubber-degrading strains in which Lcp enzymes have been characterized.

2 Materials and methods

2.1 Growth condition of Dactylosporangium sp. AC04546

Dactylosporangium sp. AC04546 is a mesophilic strain that grows well at 28°C and up to 45°C on yeast malt extract (ISP2) agar. This strain is a slow grower; when cultivated in an ISP2 broth at 28°C, 180 rpm, its exponential growth rate is from day 8 to 13 and reaches a stationary phase on day 14.

2.2 Whole genome sequencing

Cultivation and genomic extraction were conducted using methods previously described (Basik et al., 2021). Genomic DNA was sheared to a size distribution of ~20 kb and used directly for SMRTbell library construction using the SMRTbell^® Express Template Preparation Kit 2.0 (Pacific Biosciences, Menlo Park, CA, USA). Subsequently, digestion of ssDNA ends, DNA damage repair, end repair, and ligation of adapters were performed. The SMRTbell libraries were purified using AMPure^®PB and size-selected to remove hairpin and adapter dimers. DNA fragments without adapters were removed using exonucleases. The size selected SMRTbell libraries were loaded onto the PacBio Sequel II platform according to the manufacturer's instructions (Pacific Biosciences, Menlo Park, CA, USA). Sequencing was performed using the Sequel II Sequencing Kit 2.0 (Pacific Biosciences, Menlo Park, CA, USA).

Subreads generated from the PacBio Sequel II were assembled using Flye v2.9 with default parameters and “–pacbio-raw” for PacBio reads. Quast v3.1 (Gurevich et al., 2013) was used to assess the contiguity of the genome assembly, while BUSCO v4.1.4 (Simão et al., 2015), actinobacteria_class_odb10 profile was used to assess the completeness of the genome. Additionally, short-read data generated from Illumina MiSeq were recruited to further improve the assembly. The short-read sequences were trimmed using Trimmomatic version 0.38 (Bolger et al., 2014) to remove adapters and low-quality (Phred <20) and short reads (<50 bp). Gapfiller v1.10 was used to identify and close the gap in the original assembly with parameters: -m 29 -o 15 -r 0.8 -n 20 -t 10 -T. The gap-filled assembly was polished with the clean short-read data using Pilon v1.24 with default parameters. Circlator v1.5.5 was used to circularize the final assembly with parameters: -b2r_length_cutoff 10,000. The circularized assembly was assessed again using Quast and BUSCO, as mentioned above. A circular plot was plotted using ClicO FS, an online server available at http://103.47.253.210:3000/.

2.3 Phylogenetic analyses

The culture and genome of strain AC04546 were analyzed through methods described in Basik et al. (2021). The results were analyzed by the TYGS (Meier-Kolthoff et al., 2022) for (i) determination of closely related type strains, (ii) pairwise comparison of genome sequences, and (iii) phylogenetic inference.

2.4 Genome structural annotation

The complete genome was structurally annotated wherein protein-coding genes or CDS, tRNAs, and rRNAs were predicted from the genome using Prodigal v2.6.3 with parameters: -c -m; tRNAscan-SE v1.3.1 with default parameters; and rnammer v1.2 (Lagesen et al., 2007) with parameters: -S bac and -multi respectively. After gene prediction, the full repertoire of peptide sequences (≥33 aa) was assessed for its completeness using BUSCO v4.1.4 against the Actinobacteria_class_odb10 profile. The predicted peptide sequences annotated in the genome assembly were searched against RefSeq and SwissProt protein database using DIAMOND binary v2.0.13 (Buchfink et al., 2015) and BLAST+ v2.10.1 (Camacho et al., 2009) with parameters: -max_target_seqs 20 -max_hsps 1 -evalue 1e-5, respectively. The peptide sequences were further functionally annotated through orthology assignment using eggnog-mapper v2.0.0 based on eggNOG 5.0 (Huerta-Cepas et al., 2017, 2019). The peptides were also functionally characterized with protein domains and families using InterProScan v5.44-79.0 (Jones et al., 2014).

2.5 Comparative genome analyses

The newly assembled genome of strain AC04546 (GCF_022095675) and genome sequences of nine additional strains of interest (GCF_024268275, GCF_014647515, GCF_001278095, GCF_020826795, GCF_004001265, GCF_000523235, GCF_000247715, GCF_003104955, and GCF_020885655) were downloaded from NCBI RefSeq database and included in the analysis. All ten genome sequences were assessed for their contiguity and completeness using Quast v3.1 and BUSCO v4.1.4 with parameters as previously described for strain AC04546. Next, genes present or absent from selected genomes of interest were identified using Roary v3.7.0 with parameters: -i 70 -e -n -v -r -p 4. Operons were predicted for all ten genomes using Operon-Mapper, an online web server available at https://biocomputo.ibt.unam.mx/operon_mapper/.

2.6 Data mining

BLASTP v2.10.0+ was used to search for homologs of the reference genes of interest with a general e-value cut-off of 1e-05. Next, a manual assessment of the Blast results was carried out. Integrative Genomics Viewer (IGV) v2.10.2 was used to visualize the rubber-degrading and adjacent genes. Protein-coding genes (CDS) and open reading frame (ORF) on the genome were analyzed and viewed using Rapid Annotation of Microbial Genomes using Subsystems Technology (RAST) (https://rast.nmpdr.org/) and SnapGene Viewer (https://www.snapgene.com/). Molecular weights were also determined using SnapGene Viewer.

Pairwise sequence alignment (PSA) with Matrix EBLOSUM62 using EMBOSS Needle (Madeira et al., 2022) was made for strain AC04546 lcp amino acid sequences (https://www.ebi.ac.uk/jdispatcher/psa/emboss_needle). Lcp 3D structure and identity match were predicted using NCBI BLAST+ Protein similarity search using AlphaFold DB (https://alphafold.ebi.ac.uk/). Multiple sequence alignment of amino acids was viewed using Jalview 2.11 (Martin et al., 2020).

2.7 Expression of lcp genes in E. coli

The primers were designed using In-Fusion Cloning Primer Design Tool v1.0 (https://www.takarabio.com/learning-centers/cloning/primer-design-and-other-tools) with pCold IV as the destination vector and NdeI as the restriction digest. The lcp sequence without twin-arginine translocation (TAT) signal sequence were added as inserts. TAT sequence regions were identified using SignalP 5.0 software (https://services.healthtech.dtu.dk/services/SignalP-5.0/) (Almagro et al., 2019).

The coding regions of lcp1, lcp2, and lcp3 were independently amplified by PCR using lcp1_F and R, lcp2_F and R, and lcp3_F and R primer pairs:

The forward primer was added with an NdeI site to the start codon of the target gene. Amplified fragments were purified using the Fast Gene Plasmid Mini Kit. They were separately cloned into a pCold IV vector, which was subsequently digested with NdeI by using an In-fusion HD Cloning Kit (Takara Bio Inc., Tokyo, Japan).

The resultant plasmids were independently introduced into E. coli Rosetta-gami B(DE3)pLysS and the transformants were grown at 37°C in Luria Bertani (LB) medium containing ampicillin (100 mg/mL), kanamycin (25 mg/mL), chloramphenicol (25 mg/mL), and tetracycline (10 mg/mL). When the absorbance at 600 nm (A600) of the culture reached 0.5, 1 mM isopropyl-β-D-thiogalactopyranoside (IPTG) was added, and the cultures were further incubated at 15°C for 24 h. The resulting cells were harvested and resuspended in 50 mM phosphate buffer (pH 7.4). The cells were disrupted using an ultrasonic disrupter, and the lysates were centrifuged. The supernatant was collected, treated with streptomycin sulfate 10%, and filtered through a 0.45 μm filter.

The supernatant was then purified using the ÄKTA^TM Start Protein Purification System with a HiTrap TALON Crude 1 mL Superflow column (Cytiva, Uppsala, Sweden) previously equilibrated with 50 mM sodium phosphate (pH 7.4) containing 300 mM NaCl. Proteins were allowed to bind for 1 min at 4°C under rotation, followed by washing five times in 5 mL of the same buffer. His-tagged proteins were eluted with 5 mL of 50 mM sodium phosphate (pH 7.4) containing 300 mM NaCl and 150 mM imidazole. The fractions were pooled and concentrated using 30 kDa Amicon (R) Ultra Centrifugal filters.

2.8 Enzyme assay

To test whether the Lcp expressed by recombinant E.coli pCold IV::lcp strains can degrade poly (cis-1,4-isoprene) into aldehydes, crude extracts of Lcp1, 2, and 3 were applied onto a 0.4% deproteinized natural rubber (DPNR)-overlay agar plate and incubated overnight at 30°C. Three mL of Schiff's reagent was applied onto the plate and left for 30 min at room temperature. Destaining was carried out using 5 mL of 10% sodium hydrogen sulfide.

3 Results and discussion

3.1 Whole genome sequencing and assembly

The whole genome of Dactylosporangium sp. AC04546 was sequenced using the PacBio Sequel II platform. The total throughput for the sample was ~115.2 k raw reads or ~1.277 Gb of data, with 98 × genome coverage. A circularized genome comprising one single contig with a total genome size of 13.08 Mb and an average GC content of 72.1% was generated from the hybrid assembly (Genbank submission CP139171). No plasmid was detected. The circularized genome was estimated to achieve 98.3% completeness based on the Actinobacteria class odb10 profile. A total of 11,865 coding sequences (>99 bp), 9 rRNAs, and 83 tRNAs were annotated to the genome (Table 1). The coding sequences (CDS) were categorized into 14% (1,616 CDS) with subsystem categories and 86% (10,740 CDS) without subsystem categories. A subsystem is a set of functional roles that together implement a specific biological process or structural complex (Overbeek et al., 2014). Strain AC04546 has subsystems such as the metabolism of aromatic compounds and sulfur metabolism, which is an added advantage as a rubber degrader should the enzymes be expressed.

Table 1

Table 1. Genome assembly and annotation statistics of strain AC04546.

The coding sequences were functionally annotated to four databases, namely the RefSeq, SwissProt, EggNOG v5.0.0, and InterPro v79.0, with 98.55% of proteins having at least one functional annotation (Table 2).

Table 2

Table 2. Functional annotation statistics of strain AC04546.

3.2 Phylogenetic analyses

Based on its complete genome, strain AC04546 is closely matched to D. sucinum JCM19831, a strain isolated from a soil sample collected from a peat swamp forest soil in Rayong Province, Thailand (Phongsopitanun et al., 2015). There are no reports of D. sucinum for its rubber-degrading activity or the presence of lcp genes. Genome-based phylogeny shows that they are within the same clade (TYGS, https://tygs.dsmz.de) (Figure 1). Bacterial species are defined as similar species when they present more than 95% average nucleotide identity (ANI) and 70% of digital DDH (dDDH) compared to another species in the same genus (Thompson et al., 2013; Volpiano et al., 2021). Digital ddH for strain AC04546 was at only 52.4%, and an ANI value of 91.51% was obtained. Therefore, strain AC04546 is a separate species from Dactylosporangium sucinum.

Figure 1

Figure 1. Tree inferred with FastME 2.1.6.1 (Lefort et al., 2015) from Genome BLAST Distance Phylogeny (GBDP) distances calculated from strain AC04546 and closely related strains genome sequences. The branch lengths are scaled in terms of the GBDP distance formula d5. The numbers above the branches are GBDP pseudo-bootstrap support values > 60% from 100 replications, with an average branch support of 98.6%. The tree was rooted at the midpoint (Farris, 1972).

3.3 Comparative genome analyses

In addition to the current genome sequence, a total of nine additional strains of interest comprising the reference strain D. sucinum, seven strains belonging to type 1 with latex-clearing protein genes (Streptomyces sp. AC04842, Streptomyces sp. CFMR 7, Microtetraspora sp. AC03309, Actinoplanes sp. OR16, Gordonia polyisoprenivorans VH2, Nocardia nova SH22a, and Rhodococcus sp. RDE2), and one strain belonging to type 2 with roxA and roxB genes (Steroidobacter cummioxidans strain:35Y), were included for comparative analysis. The genome sequences downloaded from the NCBI RefSeq database consist of assemblies at the complete, scaffold, and contig levels, with genome sizes ranging from 5.72 to 13.08 Mb and contig or scaffold numbers ranging from 1 to 579. Based on the BUSCO Actinobacteria class odb10 profile, the genome completeness ranged from 75.7 to 100.0% (Table 3). These rubber-degrading strains were isolated from samples collected in Japan, Vietnam, and Brazil. Isolated samples ranged from soil, plant roots, and wastes from rubber processing facilities.

Table 3

Table 3. Sequence metrics and features of ten genomes of interest.

As most strains are taxonomically distant, the gene presence/absence analysis resulted in a total of 62,369 orthologous gene clusters for the ten strains based on a 70% identity cut-off, with 80.47% of them unique to only one strain. Through operon prediction, a total of 2,996–6,605 operons were predicted for the ten strains, which appear to correlate with their genome sizes.

Based on Figure 2, a comparative genome analysis between strain AC04546 and the closely related genome of the Dactylosporangium sucinum genome showed that the genome organization is very different. This is clearly depicted by the same-colored blocks, which indicate the segments, or Locally Collinear Blocks (LCBs), that are conserved among the bacteria compared (linked by the same line). On the other hand, areas that are completely white were not aligned and possibly contained sequence elements specific to a particular genome.

Figure 2

Figure 2. Genome comparison between strain AC04546 and Dactylosporangium sucinum. The alignment display for each genome sequence is organized in a horizontal panel together with the scale showing the sequence coordinates of that genome. Each block outlines the regions of the genome sequence aligned to a portion of another genome. The block above the center line indicates that the aligned region is in the forward orientation relative to the first genome sequence. In contrast, the block below the center line shows that the aligned region is in the inverse orientation. The same-colored blocks indicate the segments, or Locally Collinear Blocks (LCBs), that are conserved among the two bacteria compared.

3.4 Gene mining and visualization

Operons involved in rubber degradation across ten strains were compared by mining the latex clearing protein (lcp) genes, its adjacent genes isoquinoline 1-oxidoreductase subunit alpha (oxiA), isoquinoline 1-oxidoreductase subunit beta (oxiB) and the rubber oxygenase genes (roxA and roxB) (Figure 3). The number of lcp genes detected for all ten strains of interest is consistent with the previous reports. Strain AC04546 (Basik et al., 2021), Streptomyces sp. CFMR 7 (Nanthini et al., 2017), and Actinoplanes sp. OR16 (Gibu et al., 2020) contained three copies of lcp genes, while Streptomyces sp. AC04842 (Basik et al., 2022), Microtetraspora sp. AC03309 (Basik et al., 2022), and Gordonia polyisoprenivorans VH2 (Hiessl et al., 2012) were found to contain two copies of lcp. Nocardia nova SH22a (Luo et al., 2014) and Rhodococcus sp. RDE2 (Gibu et al., 2022) contained only one lcp gene each. On the other hand, the Gram-negative Steroidobacter cummioxidans strain 35Y (Sharma et al., 2018) was found to only contain roxA and roxB but not lcp.

Figure 3

Figure 3. Schematic representation of operon involved in rubber degradation across ten selected strains.

The oxiA and oxiB genes were not detected for strains containing mycolic acid, including Nocardia, Gordonia, and Rhodococcus. Mycolic acid-containing strains produce biofilms where rubber-degrading enzymes are being secreted into; biofilms enable the microbial cells to attach onto hydrophobic surfaces such as rubber (Linos et al., 2000). Clear zone formers, on the other hand, have mycelia which penetrates the rubber surfaces and secrete rubber-degrading enzymes through the mycelial corridors (Röther et al., 2017). There could be a correlation between OxiAB and the presence or absence of clear zones. It has been reported that in Nocardia and Gordonia, which do not generate clear zones, aldehyde dehydrogenase (ALDH) is involved instead of OxiAB. Whereby OxiAB is thought to oxidize oligomers extracellularly because it contains TAT signal sequences, but the ALDH have differences that function intracellularly (Rose et al., 2005; Suzuki et al., 2022). It was suggested that the OxiAB enzyme is required for the clear-zone formation on natural rubber in strain AC04546.

Evolutionary differences of Lcps between strain AC04546 and other Actinomycetes such as Streptomyces sp. AC04842 and Microtetraspora sp. AC03309 was characterized. The amino acid sequences of lcp1, lcp2, and lcp3 for strain AC04546 were closely related to lcp3, lcp1, and lcp2 of Actinoplanes sp. OR16, respectively (Figure 4). This trend was also observed in the estimates of evolutionary divergence between the lcp sequences (Supplementary Table 1).

Figure 4

Figure 4. Phylogenetic distribution of amino acid sequences of Lcp from strain AC04546 (red dot), Streptomyces sp. AC04842, Microtetraspora sp. AC03309 and characterized Lcps from the NCBI database. The Lcp amino acid sequence from Solimonas fluminis (Gram-negative strain) was used as an outgroup. The tree was constructed using maximum likelihood, MEGA11 (Tamura et al., 2021). The bar below the phylogenetic trees represents the scale of sequence divergence.

For strain AC04546, the location and size of nucleotide and amino acid sequences of the related genes are summarized in Table 4. Based on the gene's location, the lcp3, lcp2, and oxiB-oxiA genes are considered part of an operon. The amino acid sequence identity of Lcp1-Lcp2, Lcp2-Lcp3, and Lcp1-Lcp3 is 49.2%, 45.9%, and 60.8%, respectively (Supplementary Table 2). The gene cluster, including lcp3, lcp2, and oxiB-oxiA of strain AC04546, is similar to previously reported strains containing the 3 lcp genes: Streptomyces sp. CFMR 7 (Nanthini et al., 2017) and Actinoplanes sp. strain OR16 (Gibu et al., 2020), where two lcp genes are located adjacent to each other, followed by oxiB-oxiA genes, with the other lcp gene located far apart (Figure 4). The complete genome and lcp genes for strain AC04546 were visualized and are illustrated in Figure 5. The lcp1 gene is located about 3.4 Mb apart from the lcp2, lcp3, and oxiB-oxiA gene, similar to the draft genome reported in Basik et al. (2021). No specific MCE-family genes were detected in the complete genome.

Table 4

Table 4. The location and size of nucleotide and amino acid sequences of the related genes in Dactylosporangium sp. AC04546.

Figure 5

Figure 5. Circular representation of strain AC04546 genome and features. The genome sequence is shown at the outermost ring. The G+C content is indicated in the center of the figure. The second and third outermost rings show the rRNAs and tRNAs, respectively. Locations of latex-clearing protein (lcp) genes on the genome are shown.

The amino acid sequence of oxiA and oxiB strain AC04546 showed 80.5% (100% query coverage, QC) and 79.3% (97% QC) identities with oxiA (ACTI_59660) and oxiB (ACTI_59650) of Actinoplanes sp. OR16, respectively. The amino acid sequence of oxiA and oxiB strain AC04546 showed 74.5% (100% QC) and 77.5% (98% QC) identities with oxiA (AAR26467) and oxiB (AAR26466) of Streptomyces sp. K30 respectively (Rose et al., 2005). Based on the similarity of these sequences, oxiA and oxiB are thought to be involved in the oxidation of oligo isoprene aldehyde in strain AC04546.

Lcp1, Lcp2, and Lcp3 structures were predicted using the AlphaFold Database (Figure 6) (Varadi et al., 2024). The structural and functional analysis of Lcp from Streptomyces sp. K30 identified important residues contributing to the stability of the protein, especially Arg195 and Arg202, while conserved residues Arg164, Tr168, and His198, located close to the haem cofactor, are crucial active site residues (Ilcu et al., 2017). Multiple sequence alignment of amino acid sequences of Lcps from AC04546, OR16, and K30 confirms that Lcps from AC04546 and OR16 have the conserved residues Arg164, Tr168, and His198 and 13-residue-long highly conserved region (Supplementary Figure 1) (Röther et al., 2016). It was suggested that these residues are required for the activity of these Lcps.

Figure 6

Figure 6. Prediction of strain AC04546 Lcp1, Lcp2, and Lcp3 structure (left to right) using AlphaFold Database.

A 1,209 bp gene (located from 4,125,595 bp to 4,126,803 bp) that appears to encode a TetR-type transcriptional regulator is located upstream of the lcp1 gene. It has been reported that TetR-type transcriptional regulators are involved in the regulation of Lcp expression in strain OR16 (Gibu et al., 2020), Streptomyces coelicolor A3(2) (Coenen et al., 2019) and G. polyisoprenivorans VH2 (Oetermann et al., 2019). The genes encoding TetR-type transcriptional regulators were found to be located in the vicinity of the lcp gene. These facts suggest that TetR-type regulators expressed from genes located in the vicinity of lcp1 might be involved in the regulation of the lcp expression in strain AC04546.

3.5 Heterologous expression of Lcp in E. coli

The coding regions of each lcp gene were amplified by PCR and introduced into pCold IV to fuse a 6x histidine tag at the N terminal of the genes. The lcp genes of strain AC04546 were not well expressed in E. coli BL21(DE3); therefore, pCold IV:lcp was transformed into E. coli Rosetta-gami B(DE3)pLysS. Haem-containing proteins, including Lcp have been reported to be brown to red in color (Watcharakul et al., 2016). The dark yellow coloration of the cells of the transformant and supernatant expression of each lcp gene product was observed, but crude Lcp-his could not be purified in this system (Supplementary Figures 2–4). The histidine tag at the N-terminal may cause ineffective binding during Immobilized-metal affinity chromatography (IMAC) purification.

3.6 Enzyme assay

To confirm the enzymatic activity of recombinant Lcp, crude Lcp-his extract of 5 mg/mL was spotted on DPNR-overlay agar and stained with Schiff's reagent. The appearance of purple regions at Lcp1, Lcp2, and Lcp3 positions indicates that aldehyde is produced by the reaction of Lcp from poly(cis-1,4-isoprene). Lcp2 showed more intense purple coloration compared to Lcp1 and Lcp3 (Figure 7). The coloration was not found when the crude extract of E. coli was used with an empty vector. These results strongly suggested that these Lcp can degrade poly(cis-1,4-isoprene) into an oligomer containing an aldehyde group by adding the oxygen molecule.

Figure 7

Figure 7. Detection of strain AC04546 Lcp1, Lcp2, and Lcp3 enzyme activity obtained from heterologous expression of E. coli for using Schiff's reagent for the presence of aldehyde; control plate (left) and test samples plate (right). (A) Application of crude enzyme on filter paper. (B) Application of Schiff's reagent. (C) Control plate (without any samples). (D) After destaining process.

Pairwise sequence alignment of the strain AC04546 and Actinoplanes sp. OR16 lcp amino acid sequences were conducted. The sequences of Lcp1, Lcp2, and Lcp3 from AC04546 showed 62.6%, 80.1%, and 76.4% identities with those of Lcp3, Lcp1, and Lcp2 from Actinoplanes sp. OR16, respectively. It has been reported that Lcp1 has the highest activity to degrade poly(cis-1,4-isoprene) in Actinoplanes sp. OR16 (Gibu et al., 2020). Based on the sequence similarity, it was suggested that these Lcp enzymes from strain AC04546 have the activity to degrade poly(cis-1,4-isoprene), and Lcp2 has the highest activity in strain AC04546 (see Section 3.4).

4 Conclusion

Dactylosporangium sp. AC04546, a newly discovered species, has three Lcp enzymes that are expressed. Lcp proteins of this strain did not match other known Lcp in the database, and AlphaFold prediction of all three structures of Lcp1, Lcp2, and Lcp3 differ from each other. Differences in degradation capability were also observed during the enzyme assay. Moreover, the expression levels of the three Lcp were different, with Lcp2 being more highly expressed. This is the first report of strain AC04546 and its Lcp enzymes compared to rubber-degrading enzymes from other known species. The complete genome of this strain is a useful reference for future characterization works.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI - PRJNA896744.

Author contributions

AB: Writing—original draft. NG: Investigation, Supervision, Writing—review & editing. YK: Investigation, Writing—review & editing, Supervision. S-MN: Formal analysis, Investigation, Methodology, Visualization, Writing—review & editing. TY: Funding acquisition, Project administration, Writing—review & editing. KS: Conceptualization, Project administration, Supervision, Writing—review & editing. DK: Funding acquisition, Methodology, Project administration, Supervision, Writing—review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This project was funded by the Sarawak State Government and part of this work was also supported by JST-Mirai Program Grant Number JPMJMI19E6, Japan.

Acknowledgments

We would like to acknowledge the Sarawak Biodiversity Centre and Nagaoka University of Technology of Japan for providing their support.

Conflict of interest

S-MN was employed by Codon Genomics S/B.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2024.1378082/full#supplementary-material

References

Almagro, A. J. J., Tsirigos, K. D., Sønderby, C. K., Petersen, T. N., Winther, O., Brunak, S., et al. (2019). SignalP 5.0 improves signal peptide predictions using deep neural networks. Nat. Biotechnol. 37, 420–423. doi: 10.1038/s41587-019-0036-z

PubMed Abstract | Crossref Full Text | Google Scholar

Basik, A. A., Nanthini, J., Yeo, T. C., and Sudesh, K. (2021). Rubber degrading strains: Microtetraspora and Dactylosporangium. Polymers 13:3524. doi: 10.3390/polym13203524

PubMed Abstract | Crossref Full Text | Google Scholar

Basik, A. A., Trakunjae, C., Yeo, T. C., and Sudesh, K. (2022). Streptomyces sp. AC04842: genomic insights and functional expression of its latex clearing protein genes (lcp1 and lcp2) when cultivated with natural and vulcanized rubber as the sole carbon source. Front. Microbiol. 13:854427. doi: 10.3389/fmicb.2022.854427

PubMed Abstract | Crossref Full Text | Google Scholar

Birke, J., Röther, W., and Jendrossek, D. (2017). RoxB is a novel type of rubber oxygenase that combines properties of rubber oxygenase RoxA and latex clearing protein (Lcp). Appl. Environ. Microbiol. 83, e00721–e00717. doi: 10.1128/AEM.00721-17

PubMed Abstract | Crossref Full Text | Google Scholar

Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for illumina sequence data. Bioinformatics 30, 2114–2120. doi: 10.1093/bioinformatics/btu170

PubMed Abstract | Crossref Full Text | Google Scholar

Buchfink, B., Xie, C., and Huson, D. H. (2015). Fast and sensitive protein alignment using DIAMOND. Nat. Methods 12, 59–60. doi: 10.1038/nmeth.3176

PubMed Abstract | Crossref Full Text | Google Scholar

Camacho, C., Coulouris, G., Avagyan, V., Ma, N., Papadopoulos, J., Bealer, K., et al. (2009). BLAST+: architecture and applications. BMC Bioinf. 10, 1–9. doi: 10.1186/1471-2105-10-421

PubMed Abstract | Crossref Full Text | Google Scholar

Coenen, A., Oetermann, S., and Steinbüchel, A. (2019). Identification of LcpRB A3(2) a novel regulator of lcp expression in streptomyces coelicolor A3 (2). AppI. Microbiol. Biotechnol. 103, 5715–5726. doi: 10.1007/s00253-019-09896-8

PubMed Abstract | Crossref Full Text | Google Scholar

Farris, J. S. (1972). Estimating phylogenetic trees from distance matrices. Am. Natur. 106, 645–668. doi: 10.1086/282802

Crossref Full Text | Google Scholar

Gibu, N., Arata, T., Kuboki, S., Linh, D. V., Fukuda, M., Steinbüchel, A., et al. (2020). Characterization of the genes responsible for rubber degradation in Actinoplanes sp. strain OR16. Appl. Microbiol. Biotechnol. 104, 7367–7376. doi: 10.1007/s00253-020-10700-1

PubMed Abstract | Crossref Full Text | Google Scholar

Gibu, N., Linh, D. V., Suzuki, N., Ngan, N. T. T., Fukuda, M., Anh, T. K., et al. (2022). Identification and transcriptional analysis of poly(cis-1, 4-isoprene) degradation gene in Rhodococcus sp. strain RDE2. J. Biosci. Bioeng. 133, 452–458. doi: 10.1016/j.jbiosc.2022.01.013

PubMed Abstract | Crossref Full Text | Google Scholar

Gurevich, A., Saveliev, V., Vyahhi, N., and Tesler, G. (2013). QUAST: quality assessment tool for genome assemblies. Bioinformatics 29, 1072–1075. doi: 10.1093/bioinformatics/btt086

PubMed Abstract | Crossref Full Text | Google Scholar

Hiessl, S., Schuldes, J., Thürmer, A., Halbsguth, T., Bröker, D., Angelov, A., et al. (2012). Involvement of two latex-clearing proteins during rubber degradation and insights into the subsequent degradation pathway revealed by the genome sequence of Gordonia polyisoprenivorans strain VH2. Appl. Environ. Microbiol. 78, 2874–2887. doi: 10.1128/AEM.07969-11

PubMed Abstract | Crossref Full Text | Google Scholar

Huerta-Cepas, J., Forslund, K., Coelho, L. P., Szklarczyk, D., Jensen, L. J., von Mering, C., et al. (2017). Fast genome-wide functional annotation through orthology assignment by eggNOG-Mapper. Molec. Biol. Evol. 34, 2115–2122. doi: 10.1093/molbev/msx148

PubMed Abstract | Crossref Full Text | Google Scholar

Huerta-Cepas, J., Szklarczyk, D., Heller, D., Hernández-Plaza, A., Forslund, S. K., Cook, H., et al. (2019). eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology rresource based on 5090 organisms and 2502 viruses. Nucl. Acids Res 47, D309–D314. doi: 10.1093/nar/gky1085

PubMed Abstract | Crossref Full Text | Google Scholar

Ilcu, L., Röther, W., Birke, J., Brausemann, A., Einsle, O., and Jendrossek, D. (2017). Structural and functional analysis of latex clearing protein (Lcp) provides insight into the enzymatic ceavage of rubber. Sci. Rep. 7:6179. doi: 10.1038/s41598-017-05268-2

PubMed Abstract | Crossref Full Text | Google Scholar

Imai, S., Ichikawa, K., Muramatsu, Y., Kasai, D., Masai, E., and Fukuda, M. (2011). Isolation and characterization of Streptomyces, Actinoplanes, and Methylibium strains that are involved in degradation of natural rubber and synthetic poly (cis-1, 4-isoprene). Enzyme Micr. Technol. 49, 526–531. doi: 10.1016/j.enzmictec.2011.05.014

PubMed Abstract | Crossref Full Text | Google Scholar

Jones, P., Binns, D., Chang, H. Y., Fraser, M., Li, W., McAnulla, C., et al. (2014). InterProScan 5: genome-scale protein function classification. Bioinformatics 30, 1236–1240. doi: 10.1093/bioinformatics/btu031

PubMed Abstract | Crossref Full Text | Google Scholar

Lagesen, K., Hallin, P., Rødland, E. A., Stærfeldt, H. H., Rognes, T., and Ussery, D. W. (2007). RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucl. Acids Res. 35, 3100–3108. doi: 10.1093/nar/gkm160

PubMed Abstract | Crossref Full Text | Google Scholar

Lefort, V., Desper, R., and Gascuel, O. (2015). FastME 2.0: a comprehensive, accurate, and fast distance-based phylogeny inference program. Molec. Biol. Evol. 32, 2798–2800. doi: 10.1093/molbev/msv150

PubMed Abstract | Crossref Full Text | Google Scholar

Leong, S. Y., Lee, S. Y., Koh, T. Y., and Ang, D. T. C. (2023). 4R of rubber waste management: current and outlook. J. Mater. Cycles Waste Manag. 25, 37–51. doi: 10.1007/s10163-022-01554-y

PubMed Abstract | Crossref Full Text | Google Scholar

Linos, A., Berekaa, M. M., Reichelt, R., Keller, U., Schmitt, J., Flemming, H. C., et al. (2000). Biodegradation of cis-1, 4-polyisoprene rubbers by distinct actinomycetes: microbial strategies and detailed surface analysis. Appl. Environ. Microbiol. 66, 1639–1645. doi: 10.1128/AEM.66.4.1639-1645.2000

PubMed Abstract | Crossref Full Text | Google Scholar

Luo, Q., Hiessl, S., Poehlein, A., Daniel, R., and Steinbüchel, A. (2014). Insights into the microbial degradation of rubber and gutta-percha by analysis of the complete genome of Nocardia nova SH22a. Appl. Environ. Microbiol. 80, 3895–3907. doi: 10.1128/AEM.00473-14

PubMed Abstract | Crossref Full Text | Google Scholar

Madeira, F., Pearce, M., Tivey, A. R., Basutkar, P., Lee, J., Edbali, O., et al. (2022). Search and sequence analysis tools services from EMBL-EBI in 2022. Nucl. Acids Res. 50, W276–W279. doi: 10.1093/nar/gkac240

PubMed Abstract | Crossref Full Text | Google Scholar

Martin, D., Procter, J., Soares, B., Waterhouse, A., Shehata, S., Giang, N., et al. (2020). Jalview 2.11. Manual and Introductory Tutorial.

PubMed Abstract | Google Scholar

Meier-Kolthoff, J. P., Carbasse, J. S., Peinado-Olarte, R. L., and Göker, M. (2022). TYGS and LPSN: a database tandem for fast and reliable genome-based classification and nomenclature of prokaryotes. Nucl. Acids Res. 50, D801–D807. doi: 10.1093/nar/gkab902

PubMed Abstract | Crossref Full Text | Google Scholar

Nanthini, J., Ong, S. Y., and Sudesh, K. (2017). Identification of Three Homologous Latex-clearing protein (lcp) genes from the Genome of Streptomyces sp. strain CFMR 7. Gene 628, 146–155. doi: 10.1016/j.gene.2017.07.039

PubMed Abstract | Crossref Full Text | Google Scholar

Oetermann, S., Jongsma, R., Coenen, A., Keller, J., and Steinbüchel, A. (2019). LcpR VH2—regulating the expression of latex-clearing proteins in Gordonia polyisoprenivorans VH2. Microbiology 165, 343–354. doi: 10.1099/mic.0.000755

PubMed Abstract | Crossref Full Text | Google Scholar

Overbeek, R., Olson, R., Pusch, G. D., Olsen, G. J., Davis, J. J., Disz, T., et al. (2014). The SEED and the rapid annotation of microbial genomes using subsystems technology (RAST). Nucl. Acid Symp. Ser. 42, D206–D214. doi: 10.1093/nar/gkt1226

PubMed Abstract | Crossref Full Text | Google Scholar

Phongsopitanun, W., Kudo, T., Ohkuma, M., Suwanborirux, K., and Tanasupawat, S. (2015). Dactylosporangium sucinum sp. nov., isolated from thai peat swamp forest soil. J. Antib. 68, 379–384. doi: 10.1038/ja.2014.170

PubMed Abstract | Crossref Full Text | Google Scholar

Rose, K., Tenberge, K. B., and Steinbüchel, A. (2005). Identification and characterization of genes from Streptomyces sp. strain K30 responsible for clear zone formation on natural rubber latex and poly(cis-1,4-isoprene) rubber degradation. Biomacromolecules 6, 180–188. doi: 10.1021/bm0496110

PubMed Abstract | Crossref Full Text | Google Scholar

Röther, W., Austen, S., Birke, J., and Jendrossek, D. (2016). Cleavage of rubber by the latex clearing protein (Lcp) of Streptomyces sp. strain K30: molecular insights. Appl. Environ. Microbiol. 82, 6593–6602. doi: 10.1128/AEM.02176-16

PubMed Abstract | Crossref Full Text | Google Scholar

Röther, W., Birke, J., Grond, S., Beltran, J. M., and Jendrossek, D. (2017). Production of functionalized oligo-isoprenoids by enzymatic cleavage of rubber. Micr. Biotechnol. 10, 1426–1433. doi: 10.1111/1751-7915.12748

PubMed Abstract | Crossref Full Text | Google Scholar

Sharma, V., Siedenburg, G., Birke, J., Mobeen, F., Jendrossek, D., and Prakash, T. (2018). Metabolic and taxonomic insights into the gram-negative natural rubber degrading bacterium Steroidobacter cummioxidans sp. nov., strain 35Y. PLoS ONE 13:e0197448. doi: 10.1371/journal.pone.0197448

PubMed Abstract | Crossref Full Text | Google Scholar

Simão, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V., and Zdobnov, E. M. (2015). BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210–3212. doi: 10.1093/bioinformatics/btv351

PubMed Abstract | Crossref Full Text | Google Scholar

Soares, F. A., and Steinbüchel, A. (2022). Natural rubber degradation products: fine chemicals and reuse of rubber waste. Eur. Poly. J. 165:111001. doi: 10.1016/j.eurpolymj.2022.111001

PubMed Abstract | Crossref Full Text | Google Scholar

Sørensen, L., Gomes, T., Igartua, A., Lyngstad, I. L., Almeida, A. C., Wagner, M., et al. (2023). Organic chemicals associated with rubber are more toxic to marine algae and bacteria than those of thermoplastics. J. Hazard. Mater. 458:131810. doi: 10.1016/j.jhazmat.2023.131810

PubMed Abstract | Crossref Full Text | Google Scholar

Suzuki, N., Suda, D., Ngan, N. T. T., Gibu, N., Huong, N. L., Anh, T. K., et al. (2022). Characterization of latex-clearing protein and aldehyde dehydrogenases involved in the utilization of poly(cis-1, 4-isoprene) by Nocardia farcinica NBRC 15532. Microorganisms 10:2324. doi: 10.3390/microorganisms10122324

PubMed Abstract | Crossref Full Text | Google Scholar

Tamura, K., Stecher, G., and Kumar, S. (2021). MEGA11: molecular evolutionary genetics analysis version 11. Molec. Biol. Evol. 38, 3022–3027. doi: 10.1093/molbev/msab120

PubMed Abstract | Crossref Full Text | Google Scholar

Thiemann, J. E., Pagani, H., and Beretta, G. (1967). A new genus of the actinoplanaceae: Dactylosporangium, gen. nov. Arch. Mikrobiol. 58, 42–52. doi: 10.1007/BF00691167

PubMed Abstract | Crossref Full Text | Google Scholar

Thompson, C.C., Chimetto, L., Edwards, R.A., Swings, J., Stackebrandt, E., and Thompson, F.L. (2013). Microbial genomic taxonomy. BMC Gen. 14:913. doi: 10.1186/1471-2164-14-913

PubMed Abstract | Crossref Full Text | Google Scholar

Varadi, M., Bertoni, D., Magana, P., Paramval, U., Pidruchna, I., Radhakrishnan, M., et al. (2024). AlphaFold protein structure database in 2024: providing structure coverage for over 214 million protein sequences. Nucl. Acids Res. 52, D368–D375. doi: 10.1093/nar/gkad1011

PubMed Abstract | Crossref Full Text | Google Scholar

Volpiano, C. G., Sant'Anna, F. H., Ambrosini, A., de São José, J. F., Beneduzi, A., Whitman, W. B., et al. (2021). Genomic metrics applied to Rhizobiales (Hyphomicrobiales): species reclassification, identification of unauthentic genomes and false type strains. Front. Microbiol. 12:614957. doi: 10.3389/fmicb.2021.614957

PubMed Abstract | Crossref Full Text | Google Scholar

Watcharakul, S., Röther, W., Birke, J., Umsakul, K., Hodgson, B., and Jendrossek, D. (2016). Biochemical and spectroscopic characterization of purified latex clearing protein (lcp) from newly isolated rubber degrading Rhodococcus rhodochrous strain RPK1 reveals novel properties of Lcp. BMC Microbiol. 16, 1–13. doi: 10.1186/s12866-016-0703-x

PubMed Abstract | Crossref Full Text | Google Scholar