Next Article in Journal
Roles of the 2-Oxoglutarate-Dependent Dioxygenase Superfamily in the Flavonoid Pathway: A Review of the Functional Diversity of F3H, FNS I, FLS, and LDOX/ANS
Previous Article in Journal
Structural Optimization and Biological Activity of Pyrazole Derivatives: Virtual Computational Analysis, Recovery Assay and 3D Culture Model as Potential Predictive Tools of Effectiveness against Trypanosoma cruzi
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Biosynthesis of Tetrapyrrole Cofactors by Bacterial Community Inhabiting Porphyrine-Containing Shale Rock (Fore-Sudetic Monocline)

1
Department of Geomicrobiology, Institute of Microbiology, Faculty of Biology, University of Warsaw, 02-096 Warsaw, Poland
2
Department of Environmental Microbiology and Biotechnology, Institute of Microbiology, Faculty of Biology, University of Warsaw, 02-096 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Molecules 2021, 26(21), 6746; https://doi.org/10.3390/molecules26216746
Submission received: 10 October 2021 / Revised: 28 October 2021 / Accepted: 4 November 2021 / Published: 8 November 2021

Abstract

:
This study describes for the first time the comprehensive characterization of tetrapyrrole cofactor biosynthetic pathways developed for bacterial community (BC) inhabiting shale rock. Based on the genomic and proteomic metadata, we have detailed the biosynthesis of siroheme, heme, cobalamin, and the major precursor uroporphyrinogen III by a deep BC living on a rock containing sedimentary tetrapyrrole compounds. The obtained results showed the presence of incomplete heme and cobalamin biosynthesis pathways in the studied BC. At the same time, the production of proteins containing these cofactors, such as cytochromes, catalases and sulfite reductase, was observed. The results obtained are crucial for understanding the ecology of bacteria inhabiting shale rock, as well as their metabolism and potential impact on the biogeochemistry of these rocks. Based on the findings, we hypothesize that the bacteria may use primary or modified sedimentary porphyrins and their degradation products as precursors for synthesizing tetrapyrrole cofactors. Experimental testing of this hypothesis is of course necessary, but its evidence would point to an important and unique phenomenon of the tetrapyrrole ring cycle on Earth involving bacteria.

1. Introduction

Many cofactors of key bacterial enzymes are derived from tetrapyrrole and contain a core, complexed metal ion (Fe2+, Co2+, Mg2+, and Ni2+). A well-known example is iron-containing heme (C34H32FeN4O4), which is a prosthetic group in cytochromes, catalases, and peroxidases. Another example is heme-like iron-containing siroheme (C42H36FeN4O16). It is a cofactor at the active site of sulfite reductase, a key enzyme in sulfur assimilation pathway that converts sulfite to sulfide, which is then integrated into the organic compound homocysteine [1]. Another important tetrapyrrole cofactor is cobalt-containing cobalamin (vitamin B12) (C63H88CoN14O14P). Cobalamin serves as a coenzyme in the synthesis of nucleotides and amino acids, degradation of organic compounds (e.g., alcohols and amines), dehalogenation, fermentation, and regulation of gene expression [2,3]. Cobalamin-dependent enzymes include, but are not limited to, ribonucleotide reductase, methyltransferases, and reductive dehalogenases [4,5,6]. Nickel-containing coenzyme F430 (C42H51N6NiO13) is a cofactor of methyl coenzyme M reductase, an enzyme involved in the formation of methane [7].
Microorganisms synthesize tetrapyrroles either by de novo synthesis or by scavenging them from the environment. The biosynthesis of tetrapyrrole is one of the fundamental pathways in living organisms, and the biosynthetic pathways of the abovementioned four porphyrin cofactors are closely related (Figure 1). Uroporphyrinogen III (C40H44N4O16) is a common precursor of tetrapyrrole cofactors, which is converted to precorrin-2 (C42H48N4O16) or protoporphyrin IX (C34H34N4O4) (PPIX). Precorrin-2 in turn acts as a precursor of siroheme and vitamin B12 (while PPIX acts as a direct precursor of heme (protoporphyrin-dependent pathway). Heme can also be synthesized from coproporphyrin (coproporphyrin-dependent path way) and siroheme (siroheme-dependent pathway).
This study aimed to explore the biosynthesis of tetrapyrrole cofactors by a BC that inhabits fossil porphyrin-containing shale rock (Figure 2A–C). These porphyrins found in shale rock are known as geoporphyrins or sedimentary porphyrins. We conducted our research on Kupferschiefer shale rock located in Fore-Sudetic Monocline (SW Poland). It is a highly mineralized sedimentary Lopingian rock rich in fossil organic matter. Previous studies have confirmed the presence of iron, vanadyl, nickel, and cobalt geoporphyrins in Kupferschiefer shale rock (Table 1) [8,9,10,11].
This study is part of a project investigating the effects of bacteria on sedimentary geoporphyrins. Previous studies have shown that bacteria degrades or cleaves aliphatic and aromatic substituents of nickel and vanadyl geoporphyrins, leading to the formation of modified geoporphyrins, and also degrades the tetrapyrrole ring. The end products observed due to the microbial activity on shale rock are mono-, di-, and tripyrrole compounds [12]. Table 1 shows the list of modified geoporphyrins and pyrrole compounds identified in the shale rock inhabited by bacteria. Similar results were reported earlier by Grice et al. [13] who detected, for example, the presence of 1H-pyrrole-2,5-diones.
As part of this research, we investigated eight biosynthetic pathways, including the pathways of uroporphyrinogen III, siroheme, heme, and cobalamin (Figure 1). We analyzed three known biosynthetic pathways (protoporphyrin-, coproporphyrin-, and siroheme-dependent pathways) for heme. Furthermore, we divided the cobalamin biosynthetic pathways into three—aerobic and anaerobic biosynthesis of cob(II)yrinate a,c-diamide and biosynthesis of cobalamin from cob(II)yrinate a,c-diamide. Through metagenomic and metaproteomic analyses of BC, we identified protein-encoding genes (PEGs) and enzymatic proteins expressed as protein sequences (PSs) involved in the biosynthesis of tetrapyrrole cofactors. In addition, we searched for proteins containing the studied cofactors in the BC metaproteome. In the BC examined, we did not detect either PEGs or PSs that are involved in the F430 biosynthetic pathway.

2. Results

2.1. Taxonomic Diversity of BC Inhabiting Black Shale

With a total content of 82%, nearly five phyla of bacteria were found to be predominant in the investigated BC (Figure 2D). The dominant phyla observed were Proteobacteria (60.12%) and Actinobacteria (15.62%). Bacteroidetes, Firmicutes, and Nitrospirae were the other phyla with more than 1% dominance.
The dominant classes of bacteria were Gammaproteobacteria (33.85%), Actinobacteria (15.37%), Alphaproteobacteria (10.66%), and Betaprotobacteria (9.09%) (Figure 2E), while Nitrospira (1.69%) and Bacilli (1.28%) classes were found in smaller numbers. The other identified classes of bacteria accounted for less than 1% of the total.
In the analyzed BC, Pseudomonas was observed as the predominant genus, accounting for up to 20% of all the identified genera (Figure 2F).

2.2. General Characteristics of the Metagenome and Metaproteome of BC

Metagenomic analysis of the studied BC showed that 99.6% of PEGs belonged to the Bacteria domain and the remaining PEGs (0.4%) belonged to Archaea, Eukaryote, and Virus domains (unpresented data). Metaproteomic analysis of BC showed that 100% of PSs belonged to the Bacteria domain (unpresented data). PEG-encoded enzymes involved in metabolism, environmental information processing, cellular processes, and genetic information processing accounted for 4166, 1208, 2901, and 2401 unique reads (URs), respectively (Figure 3A), whereas the metaproteome of BC consisted of 952, 311, 228, and 241 PSs that were involved in the mentioned functional categories, respectively.
To estimate microbial richness at different taxonomic levels and to assess sample saturation for genes involved in porphyrin metabolism as a function of sequencing depth, rarefaction curves were generated (Figure S9). In both cases, the rarefaction curves reached a plateau, suggesting that sufficient sequencing saturation was achieved for taxonomic and genetic evaluation for the analyzed metabolic pathways.
In the studied the BC metagenome analysis showed that 32 PEGs were involved in the major biosynthetic pathways of tetrapyrrole cofactor (Table S2). Among them, the Proteobacteria URs accounted for 91.67%, Actinobacteria URs 7.29%, and Firmicutes URs 1.04% (Figure 3B). In addition, the metaproteomic analysis showed that 66 PSs participated in the biosynthesis of these compounds, of which the PSs of Proteobacteria (83.33%) and Firmicutes (10.61%) were predominant. Acidobacteria, Actinobacteria, and Bacteroidetes together accounted for the remaining 6.06% of PSs (Figure 3C, Table S3). Most of the URs and PSs identified in the studied BC were found to be involved in protoporphyrin-dependent heme metabolism (30.21% URs and 27.27% PSs) (Figure 3B,C).

2.3. Biosynthesis of Cofactors Detected in the Studied BC

2.3.1. Uroporphyrinogen III Biosynthesis

According to the literature data, the biosynthetic pathway of uroporphyrinogen III is depicted in Figure S1 [14,15,16]. Table S1 shows the lists of enzymes involved in this pathway. Figure 4 shows the results of the metagenomic and metaproteomic analyses of the tested BC.
Uroporphyrinogen III biosynthesis occurs with the use of Glycine or L-glutamate is the basic building block utilized in the biosynthesis of the pyrrole rings of uroporphyrinogen III. In the case of glycine substrate pathway, this compound is converted to 5-aminolevulinate by 5-aminolevulinate synthase (Alas) in a single step. L-glutamate is converted to 5-aminolevulinate by glutamyl-tRNA synthetase (Ears), glutamyl-tRNA reductase (HemA), and glutamate-1-semialdehyde 2,1-aminomutase (HemL). In the metagenome of BC, PEGs of all enzymes involved in the biosynthesis of uroporphyrinogen III using L-glutamate as the starting substrate were detected. However, no PSs were detected. The steps involved in the conversion of 5-aminolevulinate to uroporphyrinogen III are common to both pathways. 5-Aminolevulinate is transformed to porphobilinogen by porphobilinogen synthase (HemB), four molecules of which form a four-pyrrole hydroxymethylbilane compound catalyzed by hydroxymethylbilane synthase (HemC). The metagenomic analysis of the study showed the presence of PEGs of both enzymes. In addition, the metaproteomic analysis showed the presence of five and nine PSs that matched these two enzymes, respectively.
In the final reaction of the biosynthetic pathway, the tetrapyrrole ring is closed by uroporphyrinogen III synthase (HemD), resulting in the formation of uroporphyrinogen III. In the studied BC, the PEG of the HemD enzyme was observed, but the matching PSs were not detected.

2.3.2. Siroheme Biosynthesis

According to the literature data, the biosynthetic pathway of siroheme is depicted in Figure S2 [17,18]. Table S1 shows the lists of enzymes involved in this pathway. Figure 4 shows the results of the metagenomic and metaproteomic analyses of the analyzed BC.
The step in which the heme and siroheme pathways are separated is the conversion of uroporphyrinogen III to precorrin 2. The biosynthesis of siroheme is a three-step process. The first step involves the formation of precorrin-2, in which the two methyl groups derived from S-adenosyl-L-methionine on carbons two and seven of uroporphyrinogen III are cleaved. The second step is the formation of sirohydrochlorin, which is a dehydrogenation (NAD+-dependent) process. The third step is the synthesis of siroheme, which involves the chelation of ferrous iron to sirohydrochlorin to form siroheme. The steps that follow are catalyzed by uroporphyrin III C-methyltransferase (CobA) or uroporphyrin III C-methyltransferase/precorrin-2 dehydrogenase/sirohydrochlorin ferrochelatase (CysG). The metagenomic analysis of the study showed the presence of PEGs of all enzymes, and metaproteomic analysis showed eight PSs matching both enzymes.

2.3.3. Heme Biosynthesis

Protoporphyrin-Dependent Heme Biosynthetic Pathway

According to the literature data, the protoporphyrin-dependent heme synthesis pathway is depicted in Figure S3 [19,20]. Table S1 shows the lists of enzymes involved in this metabolic pathway. Figure 4 shows the results of the metagenomic and metaproteomic analyses of the analyzed BC.
The enzyme uroporphyrinogen decarboxylase (HemE) shortens the four substituents derived from propanoic acid to methyl substituents in uroporphyrinogen III, resulting in coproporphyrinogen III. The metagenome of BC analyzed in this study showed the presence of PEG of HemE and four PSs matching the enzyme. Coproporphyrinogen III is converted by enzymes coproporphyrinogen III oxidase (HemF) or oxygen-independent coproporphyrinogen III oxidase (HemN) to protoporphyrinogen IX by cleaving hydrogen from two propyl acid substituents and changing them to ethynyl substituents. The presence of PEGs of both enzymes and one PSs matching HemF and four PSs matching HemN were detected in analyses in the studied BC.
Further conversions of protoporphyrinogen IX to PPIX are catalyzed by three oxidases: protoporphyrinogen/coproporphyrinogen III oxidase (HemY), menaquinone-dependent protoporphyrinogen oxidase (HemG), and protoporphyrinogen IX oxidase (HemJ). In the BC metagenome, only PEGs of the latter enzyme HemJ were detected, whereas no PSs matching the enzyme HemJ were identified. Protoporphyrin/coproporphyrin ferrochelatase (HemH), a PPIX molecule, attaches iron as the central ion, resulting in the formation of protoheme (heme B), which is required in the synthesis of cytochrome B and can transform into different types of heme. The studied BC was found to have the PEG of HemH and five PSs matching this enzyme.
In the next steps, heme B is converted to heme O by heme O synthase (Cox10) and then to heme A by heme A synthase (Cox15). HemA is used as a cofactor of cytochrome A. The study results showed the presence of PEGs of both enzymes; however, four PSs matching only enzyme Cox15 were detected.

Coproporphyrin-Dependent Heme Biosynthetic Pathway

Coproporphyrin-dependent heme synthesis is an alternative heme biosynthetic pathway that uses part of the primary heme pathway [21]. This pathway is shown in Figure S4 based on the literature data [22,23,24]. Table S1 shows the lists of enzymes involved in this pathway. Figure 4 shows the results of the metagenomic and metaproteomic analyses.
In the abovementioned pathway, protoporphyrinogen/coproporphyrinogen III oxidase (HemY) converts coproporphyrinogen III to coproporphyrin III by removing hydrogen atoms from the pyrrole rings. However, the study results showed no PEGs and PSs of HemY in the investigated BC. Similar to the protoporphyrin-dependent heme biosynthetic pathway, protoporphyrin/coproporphyrin ferrochelatase (HemH) converts coproporphyrin III to Fe-coproporphyrin III by adding an iron atom. The results of this study showed PEG of HemH and five PSs of this enzyme in the investigated BC. Further, decarboxylation of Fe-coproporphyrin III substituents by hydrogen peroxide-dependent heme synthase (HemQ) or AdoMet-dependent heme synthase (Ahbd) results in the formation of heme B. The results showed PEG matching only HemQ enzyme in the studied BC, whereas no PSs matching the enzyme HemQ or Ahbd were detected.

Siroheme-Dependent Heme Biosynthetic Pathway

Another lesser known heme biosynthetic pathway is the siroheme pathway [21]. This pathway is depicted in Figure S5 based on the literature data [25,26,27]. Table S1 shows the list of enzymes involved in this pathway. Figure 4 shows the results of the metagenomic and metaproteomic analyses.
The first step in this pathway is the transformation of siroheme to heme B, in which two CO2 groups are separated by siroheme decarboxylase (Ahbab), resulting in the formation of 12,18-didecarboxysiroheme. The newly formed compound can transform to dihydro-heme d1 (the mechanism of this reaction has not been known yet) and then to heme d1 (Nirn) with the help of the enzyme dihydro-heme d1 dehydrogenase. In the investigated BC, only PEG of Ahbab was detected, but no PSs of any of these enzymes were found.
An alternative pathway for the synthesis of heme in this pathway is the conversion of 12,18-didecarboxysiroheme to Fe-coproporphyrin III catalyzed by Fe-coproporphyrin III synthase (Ahbc). In this reaction, two methanoic acid residues, present in the form of substituents, are removed from the porphyrin ring. The results showed PEG and three PSs matching Ahbc in the studied BC.
The final step is the decarboxylation of Fe-coproporphyrin III substituents to heme B by hydrogen peroxide-dependent heme synthase (HemQ) or AdoMet-dependent heme synthase (Ahbd). Metagenomic results of this study showed only PEG of HemQ; however, no PSs matching HemQ or Ahbd enzyme were detected.

2.3.4. Cobalamin Biosynthesis

Cobalamin biosynthesis consists of two steps. Cob(II)yrinate a,c-diamide is synthesized in the first step, which is then converted to cobalamin [28]. The first step may take place under anaerobic or aerobic conditions.

Biosynthesis of Anaerobic and Aerobic Cob(II)yrinate a,c-Diamide

Anaerobic and aerobic synthesis of cob(II)yrinate a,c-diamide is shown in Figures S6 and S7, respectively, based on the literature data [28,29,30,31]. Table S1 shows the list of enzymes involved in these two pathways. Figure 4 shows the results of the BC metagenome and metaproteome analyses.
The results showed that only PEGs of CbiK and CobN were involved in anaerobic and aerobic pathways, respectively, and the metaproteome of BC had two PSs matching CobB-CbiA.

Biosynthesis of Cobalamin from Cob(II)yrinate a,c-Diamide

Biosynthesis of cobalamin from cob(II)yrinate a,c-diamide is shown in Figure S8 based on the literature data [28,32]. Table S1 shows the lists of enzymes involved in this pathway. Figure 4 shows the results of the metagenome and metaproteome analyses.
Cobalamin formation reactions can be carried out by using cob(II)yrinate a,c-diamide or cob(I)yrinate a, c diamide to form adenosyl cobyrinate a, c diamide. In both substrates, Co(III) in the corrinoid ring is reduced to Co(II) by transferring one electron. Furthermore, Co(II) is bound by corrinoid adenosyltransferase, resulting in the displacement of the lower axial substituent by an aromatic residue, as well as the transfer of the deoxyadenosyl group and oxidation of the cobalt atom to the Co(III) state. These reactions are carried out by cob(I)alamin adenosyltransferase (Mmab) in both the substrates. The results of this study showed the presence of PEG and two PSs matching Mmab.
In the next step, adenosylcobyrinic acid a,c-diamide is converted to adenosylcobyric acid by adenosylcobyric acid synthase (cobQ, cbiP), which catalyzes the four-step amidation sequence by forming triamide, tetraamide, and pentaamide as intermediates of adenosylcobyric acid. The PEG of CobQ was detected in this study, but the PS matching this enzyme was not found in the investigated BC.
The compound obtained is converted to adenosyl cobinamide by adenosylcobinamide-phosphate synthase (CbiB, CobD, CobC1). In this reaction, the lower-axis substituent of the NH group is cleaved, resulting in the formation of adenosyl cobinamide. The results of this study showed PEG of CbiB, but PS matching this enzyme was not detected in the BC.
Furthermore, adenosyl cobinamide is converted to adenosyl cobinamide phosphate and then to adenosylcobinamide-GDP by adenosylcobinamide-phosphate guanylyltransferase (CobP, CobU). An orthophosphoric acid residue is attached to the lower substituent in the first reaction, and adenosine is coupled to the attached phosphate in the second reaction. This study showed PEG of CobP and five PSs matching this enzyme.
In the final step, adenosylcobinamide-GDP is converted to cobalamin by adenosylcobinamide-GDP ribazoletransferase (CobS, CobV), which attaches alpha-ribazole to adenosylcobinamide-GDP molecule, resulting in the formation of the final product. This study showed the presence of PEG of CobS, but PS matching the enzyme CobS was not found.
An important part of this metabolic pathway is the formation of alpha-ribazole, which is added to adenosylcobinamide-GDP in the final step of cobalamin biosynthesis. Dimethylbenzimidazole is an important substrate that is converted to 5′-phosphate alpha-ribazole by nicotinate-nucleotide dimethylbenzimidazole phosphoribosyltransferase (CobU, CobT). This reaction attaches to the phospho-alpha-D-ribosyl substrate in place of the hydrogen present at the nitrogen atom. PEG and three PSs matching this enzyme were detected in the analyzed BC.
The resultant product is converted to alpha-ribazole, which is attached to the cortical ring. This reaction is catalyzed by adenosylcobalamin/alpha-ribazole phosphatase (PhpB) by cleaving the phosphoric acid residue from the substituent. The results showed PEG of PhpB, but PSs matching this enzyme were absent.

2.3.5. Transport of Exogenous Porphyrins

In the BC examined, a transport system based on the TonB periplasmic protein with auxiliary heme acquisition protein (HasA), hemoglobin/transferrin/lactoferrin receptor protein (HemR), biopolymer protein transport ExbB, and biopolymer protein transport ExbD was detected (Figure 4). PEGs of TonB, ExbB, and ExbD were identified in the BC metagenome. Furthermore, one PS of tonB, ExbB, and ExbD was identified in the metaproteome (Figure 4).

2.3.6. Tetrapyrrole Cofactor-Containing Enzymes Produced by BC

Table 2 lists the proteins containing tetrapyrrole cofactors in the BC metaproteome. Among these enzymes, (i) heme-containing cytochromes, peroxidases, and catalases; (ii) siroheme-containing sulfite reductases, and (iii) cobalamin-containing methyltrasferases and dehalogenases were detected.

3. Discussion

One of the basic elements of metagenome- and/or metaproteome-based studies of microbial metabolism is the characteristics of the tetrapyrrole cofactor biosynthesis [33,34,35]. However, most of the available literature data describe the biosynthesis of a single tetrapyrrole cofactor, with no detailed analysis of the biosynthesis of the four tetrapyrrole cofactors. Furthermore, a majority of the studies are limited to the analysis of single strain cultivated under laboratory conditions, and no data are available on the synthesis of cofactors by a community of microorganisms in their natural environment. There is also no mention of these metabolic processes with regard to the natural environment in which the microorganisms/communities are found and the potential use of chemical compounds—precursors—present in this environment. In this regard, our study is unique in terms of research material (BC), environment (sedimentary shale rock containing fossil tetrapyrroles), and research methodology. It is also worth emphasizing that the results of our previous geochemical studies of the environment complement the molecular analysis and provide a broader view of the processes under investigation.
The metagenomic and metaproteomic analyses of the tested BCs revealed the existence of tens of thousands of URs and hundreds of PSs involved in the basic metabolism (Figure 3A). Among them, 0.7% URs and 7.6% PSs were involved in the metabolism related to the biosynthesis of tetrapyrrole cofactors, the products of which are essential for many enzymatic reactions. Much higher results were obtained by Feng et al. [36], who employed the same method of analysis for examining the metagenome derived from soil from phosphate-producing chemical plants. Their study showed that approximately 4% of URs were involved in the biosynthesis of tetrapyrrole cofactors. Similar results were obtained by Ghosh and Das [37], who also used the KEGG base to perform a functional investigation of the microbiome from manganese-contaminated mine rocks.
The results of the present study showed that five of the eight metabolic pathways have a set of URs (uroporphyrinogen III, siroheme, protoporphyrin-dependent heme, siroheme-dependent heme, and cobalamin main biosynthetic pathways), including one set of enzymatic proteins (siroheme biosynthetic pathway) (Figure 4 and Figure S5). The remaining pathways lack at least one gene. For example, 66% of PEGs and 33% of PSs were identified in the coproporphyrin-dependent pathway of heme biosynthesis. However, proteins containing all of the examined cofactors were found in BC, and hemoproteins were particularly abundant (Table 2).
When interpreting the above results, one should consider the following questions: (1) do incomplete pathways result from difficulties in DNA isolation and insufficient depth of sequencing? (2) Why do bacteria produce main pathway enzymes if they do not synthesize the basic precursor, i.e., cob(II)yrinate a,c-diamide? (3) What is the role of enzymes that are identified in incomplete pathways? (4) What is the taxonomic origin of the identified PEGs and the relationship between the individual pathways? (5) How can the results obtained be linked to the geochemistry of the environment (shale rock) inhabited by BC?
The tested BC sample can be described as difficult, first of all in the context of DNA isolation. This is evidenced by the amount of sample used for DNA isolation in order to obtain the amount necessary for sequencing. On the other hand, the sequencing results (saturation of the readings) indicate the proper depth of sequencing (Figure S9). At the same time, it should be noted that similar results (incomplete pathways) were obtained for communities of marine bacteria, soil bacteria or even pure bacterial cultures, for which there were no difficulties with DNA isolation and sequencing [35,38,39,40,41].
For example, Doxey et al. [38] examined 430 metagenomes for the presence of genes encoding cobalamin and heme biosynthesis enzymes in marine Thaumarchaeota. No genes coding for: sirohydrochlorin cobaltochelatase, precorrin-6A/cobalt-precorrin-6A reductase, precorrin-3B synthase, precorrin-6A synthase, precorrin-6B C5,15-methyltransferase/cobalt-precorrin-6B C5,C15-methyltransferase, cobaltochelatase CobN participating in the aerobic and anaerobic synthesis of cobyrinate a, c-diamide were detected. The following genes were also not detected: cob(I)alamin adenosyltransferase, adenosylcobinamide kinase/adenosylcobinamide-phosphate guanylyltransferase, nicotinate-nucleotide—dimethylbenzimidazole phosphoribosyltransferase i alpha-ribazole phosphatase participating in the main pathway of cobalamin biosynthesis from cobyrinate a, c-diamide: Despite the lack of genes encoding the necessary enzymes for the biosynthesis of cobalamin, the biosynthesis of this compound takes place. Similar results were obtained for the heme biosynthesis. Only genes encoding enzymes responsible for the biosynthesis of uroporphyrinogen-III and converting heme B to heme H and O were discovered.
In turn, Pan et al. [35] analyzed the metagenomic DNA of Bathyarchaeota found in mangrove sediments. Analysis of metagenomes revealed the presence of only about half of the genes encoding enzymes involved in cobalamin synthesis, including the cobalt chelatase gene, and only some of the genes encoding enzymes responsible for heme synthesis.
Similar results were also obtained by Panek and O′Brian [39], who studied the genomes of pure cultures of microorganisms, in which the genes encoding all enzymes of the heme biosynthetic pathway were also not detected, despite their ability to synthesize this cofactor. Among the examined genomes of bacteria, the gene encoding protoporphyrinogen IX oxidase was most often missing. This applies to strains belonging to α-Proteobacteria, e.g.,: Agrobacterium tumefaciens, Brucella melitensis, Sinorhizobium meliloti, β-Proteobacteria, e.g.,: Neisseria meningitidis (MC58 and Z2491), Ralstonia solanacearum, σ/ε-Proteobacteria, e.g.,: Campylobacter jejuni, Helicobacter pylori (26695 and J99), γ-Proteobacteria, e.g.,: P. aeruginosa, Xylella fastidiosa. The lack of the gene encoding protoporphyrinogen IX oxidase is also detected in heme biosynthesizing microorganisms that do not belong to the Proteobacteria, e.g., Nostoc sp. PCC 7120, Synechocystis PCC 6803. In addition, strains lacking the gene encoding protoporphyrinogen IX oxidase and uroporphyrinogen-III synthase were identified in α-Proteobacteria, e.g., Caulobacter crescentus, Rickettsia conorii, R. prowazekii.
Similar results were obtained by Glanville et al. [40], who in the genome of P. aeruginosa identified only selected genes encoding enzymes of the heme and cobalamin biosynthetic pathway.
According to Lu et al. [41], only about 10% of 155 soil microbial metagenomes contained complete cobalamin biosynthetic pathway. Similar results were obtained by Balabanova et al. [42], who found that the primary producers of cobalamin in the soil microbiome were bacteria of the genera Proteobacteria, Actinobacteria, Firmicutes, and Nitrospirae, which made up a small percentage of the community. The remaining microorganisms were only able to transport cobalamin.
In conclusion, as shown in the above studies, many bacteria do not have the complete heme biosynthetic pathways. The majority of the genomes examined lacked one or more gene-encoding enzymes involved in heme biosynthetic pathways. The same is true for cobalamin biosynthesis.
One explanation for this phenomenon is that there are still undiscovered pathways in the biosynthesis of tetrapyrrole cofactors in bacteria [43]. However, incomplete pathways, including pathways in which microorganisms uptake readymade tetrapyrrole cofactors such as heme [44,45,46] and cobalamin [2,41,42] from the environment, have also been established. These studies have emphasized that the microorganisms transport cofactors from external sources. For example, living cells can discharge cobalamin into the soil or release it via cell lysis. Similarly in the case of heme, microorganisms secrete extracellular hemoproteins, e.g., extracellular peroxidases, which are then taken up by other microorganisms [47,48].
Furthermore, according to the literature data, the addition of certain chemical compounds to microbial medium increases the production of tetrapyrrole cofactors. For example, bacteria have been found to utilize PPIX as a source of the tetrapyrrole ring. These bacteria include, among others, bacteroids [49], Lactococcus lactis [50], and Haemophilus influenzae [51,52]. They have a shorter heme biosynthetic pathway that involves only ferrochelatase reactions, in which an iron atom is incorporated into the core of PPIX, thereby yielding heme. However, how PPIX are transported into the bacterial cell is still unknown, although it is believed that the transport of these compounds takes place via heme transport proteins [45]. Similarly, Jacobs et al. [53] found that adding coproporphyrinogen to the microbial medium increased the heme concentration generated by P. fluorescens, P. denitrificans, and Escherichia coli.
Similar conclusions can be drawn by analyzing cobalamin biosynthesis. In BC, we detected the entire biosynthetic pathway of this key compound from cob(II)yrinate a,c-diamide, with no enzymes involved in the process. According to a study by Lawrence et al. [54], cob(II)yrinate diamide may be extracted from the medium and converted to cobalamin. However, the role of CobB-CbiA in catalyzing the aerobic and anaerobic biosynthetic pathways of cob(II)yrinate a,c-diamide is worth mentioning. Analysis of the metagenome and metaproteome of BC showed that no genes and enzymes preceded or followed the reaction carried out by this enzyme. Two PSs of this enzyme were detected in the metaproteome. Various studies have shown that CobB-CbiA is utilized in the acetylation of proteins in the metabolism of microorganisms [55,56]. Because acetylated proteins are necessary for the normal functioning of metabolism in the cell, this enzyme can therefore be produced independently of cobalamin biosynthesis [57].
The literature data described above lead us to hypothesize whether bacteria can use organic compounds found in sedimentary rocks as tetrapyrrole precursors.
We obtained extremely important information with regard to this hypothesis from Demopoulos et al. [58] and Robinson [59], who showed that monopyrrole, dipyrrole, and tripyrroles compounds can be chemically converted into porphobilinogen and tetrapyrroles. Demopoulos et al. [58] showed the ability of pyrrole, pyrrole-3-acetic acid, 4-acylpyrrole-2-carboxyaldehydes, and pyrrole-2-carboxaldehyde to form porphobilinogen. In addition, Robinson [59] confirmed that uroporphyrinogen can be chemically synthesized from porphobilinogen (monopyrrole compounds) and 3-[5-[(Z)-[(5Z)-4-ethenyl-5-[(4-ethenyl-3-methyl-5-oxopyrrol-2-yl) methylidene]-3-methylpyrrol-2-ylidene]methyl]-2-hydroxy-4-methyl-1H-pyrrol-3-yl] propanoate (tripyrrole compounds).
Taking into account the above information from the literature, it appears that the tested BC may also utilize pyrrole compounds formed as a result of geoporphyrin degradation. As mentioned, the Kupferschiefer shale contains a variety of geoporphyrins including cobalt protoporphyrin (IX) and pyrrole compounds that are formed as a result of bacterial degradation of sedimentary primary geoporphyrins (Table 1).
Finally, it is worth paying attention to the taxonomic origin, which indicates the dominant group of bacteria involved in the biosynthesis of tetrapyrrole cofactors. The results obtained with regard to the biosynthesis of tetrapyrroles revealed that Proteobacteria accounted for 91.67% of URs and 80.7% of PSs (Figure 3B,C). This is most probably related to the dominance of these organisms in BC, constituting more than 60% of the population (Figure 2D). Furthermore, Proteobacteria is one of the major phyla of bacteria characterized by its ability to synthesize cobalamin and iron porphyrins (heme, siroheme) [21,60]. Thus, the results obtained by us are consistent with the data obtained by other research teams.
To summarize, the metagenomic and metaproteomic analyses carried out in this study allowed us to gain a better understanding of the basic metabolism of tetrapyrrole cofactors by BC inhabiting the shale rock. Based on the results obtained, as well as on the literature data, one may consider the possibility of utilizing fossil primary or modified sedimentary porphyrins, as well as their bacterial degradation products, as precursors of tetrapyrrole cofactors. This hypothesis, of course, needs further evaluation, but its evidence would indicate an important and unique phenomenon of the tetrapyrrole ring cycle on Earth involving bacteria.

4. Materials and Methods

4.1. Site and Sample Description

Three BC-inhabited shale rock samples (Figure 2A–C) were aseptically collected at 32 °C and under 1080 hPa from a depth of approx. 800 m below sea level in the Lubin copper mine (KGHM Polska Miedź,Lubin, Poland). The pH of the sample was in the range of 8.2–8.4. Until laboratory processes, the collected samples were stored at −80 °C or −4 °C.

4.2. Isolation of DNA

DNA was isolated from the three BC samples according to the modified procedure of Zhou et al. [61]. Briefly, 100 g of sample was resuspended in 100 mL of DNA extraction buffer (0.1 M Na2EDTA, 0.1 M Tris–HCl, 0.1 M Na2HPO4, 1.5 M NaCl, 1% hexadecyltrimethylammonium bromide, pH 8.0) containing proteinase K and lysozyme (75 µL; 10 mg/mL) and was incubated overnight at 37 °C with horizontal shaking. Then, 20% sodium dodecyl sulfate was added, and the samples were incubated for 4 h at 65 °C and then centrifuged (6000 g, 10 min). The supernatants were harvested, mixed with an equal volume of chloroform/isoamyl alcohol mixture (24:1), and centrifuged again (6000 g, 10 min). After centrifugation, the aqueous phases were collected and precipitated overnight with 0.6 volume of isopropanol at room temperature. Next, the samples were centrifuged once again (16,000 g, 20 min, 4 °C), and the obtained pellets were washed with 70% cold ethanol and dried. Finally, DNA was isolated, resuspended in 50 µl of sterile deionized water, and stored at −80 °C until analysis.

4.3. Sequencing and Analysis of DNA

A barcoded library was prepared using the isolated DNA with Ion Xpress™ Plus Fragment Library Kit (Thermo Fisher, MD, USA) according to the manufacturer’s instructions. The library was then clonally amplified on Ion One Touch 2 system (Thermo Fisher) using Ion PI™ Template OT2 200 Kit v2 (Thermo Fisher) and sequenced on Ion Proton sequencer using Ion PI™ Sequencing 200 Kit v2 (Thermo Fisher) following the manufacturer’s instructions. The reads were demultiplexed with Torrent Suite software.
Qualitative preprocessing of the reads was performed with the Trimmomatic tool (v. 0.38) [62] (SLIDINGWINDOW:4:15, HEADCROP:3, CROP:250, and MINLEN:35). The quality of both pre- and postprocessed reads was assessed using FastQC (v. 0.11.08). Human reads were removed using BMTagger (v 1.1.0) [63] and human genome database GRCh38/hg38. Prior to the taxonomic analysis, artificial duplicate reads were removed using the k-mer approach. The obtained unassembled reads were assigned to the taxa using Kraken2 (v. 2.0.8-beta) [64] and the proteins of prokaryotic, fungal, and viral sequences of the NCBI RefSeq database (ftp://ftp.ncbi.nlm.nih.gov/refseq/, accessed on 7 November 2019. To perform functional analysis the reads were assembled using SPAdes (v.3.11.1) [65] with the “ion torrent” flag in the “careful” mode. Then, the assemblies were annotated using PROKKA (v. 1.13) [66] with default parameters and KofamScan (v. 1.3.0) [67] with KOfam HMM profiles (v. 26.04.2021). Hits having a bitsore of <60 and an e-value of >1e−5 were discarded. Based on the ko numbers and KEGG Orthology classification, the genes were assigned to metabolic pathways. The sequencing reads were realigned back to contigs using BBMap (v. 38.90) [68] with a “vslow” parameter. Duplicates were identified and removed using MarkDuplicates from the Picard tool kit. Only the BC sequences with unique best reads (URs) were taken into account, and duplicates were removed. The coverage for each predicted gene was normalized, expressed as transcripts per million [69], and averaged. Taxonomy for each gene was determined using Kraken2. Rarefaction curves for taxonomic analysis were generated using the krakefaction software (v.0.2.0) (https://github.com/phac-nml/krakefaction) and plotted as the number of unique phyla and genera as a function of sequencing depth. Rarefaction analysis carried out for complex genes that were assigned to the Porphyrin and chlorophyll metabolism pathway (KEGG Pathway: map00860) was performed using the RarefactionAnalyzer tool from the AMRPlusPlus pipeline [70] with 10 iterations per sampling level and 80 gene fraction threshold. The obtained data were plotted as the number of genes as a function of sequencing depth. The results of the metagenomic DNA sequencing are deposited at: http://www.ebi.ac.uk/ena/data/view/PRJEB21520 (ENA-FIRST-PUBLIC: accessed on 9 January 2019).

4.4. Isolation of Proteins

Proteins were extracted from the three BC samples as described by Ram et al. [71] Briefly, 100 g of sample was resuspended in 120 mL of 20 mM Tris–HCl (pH 8), shaken for 3 min, and sonicated on ice for up to 10 times for 1 min, with 1-min pauses in between (Sonics Vibracell; LABOPLUS, ModelCV18head, Sonics & Materials, CT, USA). The unlysed cells and fragments of the cell membrane were removed by adding 100 mL of 0.4 M Na2CO3 (pH 11) and centrifuging the suspension (6000 g, 20 min, 4 °C). The obtained supernatant was filtered through a 0.22-µm filter, and trichloroacetic acid was added at a ratio of 1:10 (v/v). The mixture was incubated overnight at 4 °C to allow for protein precipitation and then centrifuged again (20,000 g, 10 min, 4 °C). After centrifugation, the aqueous phase was discarded, and the pellet containing precipitated protein was resuspended in 0.5 mL of prechilled methanol (4 °C) and centrifuged once again (20,000 g, 10 min, 4 °C). The precipitate obtained was dried and stored at −80 °C. The sample preparation was performed in triplicate.

4.5. Identification of Proteins

Proteins were identified by liquid chromatography coupled to tandem MS (LC–MS-MS/MS) using nanoACQUITY (Waters, MA, USA) LC system and Orbitrap Velos mass spectrometer (Thermo Electron Corp., San Jose, CA, USA).
Identification was carried out at the Environmental Laboratory of Mass Spectrometry, Institute of Biophysics and Biochemistry (Polish Academy of Sciences, Warsaw, Poland). The equipment used for the analysis was sponsored in part by the Center for Preclinical Research and Technology (CePT), a project co-sponsored by the European Regional Development Fund and Innovative Economy, the National Cohesion Strategy of Poland.
Before identification, the proteins were subjected to a standard “in-solution digestion” procedure, during which they were reduced with 50 mM Tris(2-carboxyethyl)phosphine (60 min, 60 °C), alkylated with 200 mM S-methyl methanethiosulfonate (45 min, room temperature), and digested overnight with trypsin (Sequencing Grade Modified Trypsin; Promega V5111). The obtained peptide mixture was applied to an RP-18 precolumn (nanoACQUITY Symmetry® C18; Waters 186003514), with water containing 0.1% trifluoroacetic acid as the mobile phase. Then, the mixture was transferred to a nano-HPLC RP-18 column (nanoACQUITY BEH C18; Waters 186003545) using an acetonitrile (ACN) gradient (5%–35% ACN in 180 min) with 0.05% formic acid (flow rate 250 µL/min). The column outlet was directly coupled to the ion source of the spectrometer, which was working in the regime of data-dependent MS to MS/MS switch. Each analysis was preceded by a blank run to rule out cross-contamination from previous samples.
The acquired raw data were processed by Mascot Distiller followed by Mascot Search (Matrix Science, London, UK; on-site license) against the NCBI protein database. Peptides having a Mascot score exceeding a threshold value corresponding to <5% expectation value, as calculated by the Mascot procedure, were considered as positively identified.
Metaproteomic analysis was performed according to Kanehisa et al. [72] using GhostKOALA automatic annotation and KEGG mapping service.

Supplementary Materials

The following are available online. Figure S1. Biosynthetic pathway of uroporphyrinogen III. Figure S2. Biosynthetic pathway of siroheme. Figure S3. Biosynthetic pathway of protoporphyrin-dependent heme. Figure S4. Biosynthetic pathway of coproporphyrin -dependent heme. Figure S5. Biosynthetic pathway of siroheme-dependent heme. Figure S6. Biosynthetic pathway of anaerobic cob(II)yrinate a,c-diamide. Figure S7. Biosynthetic pathway of aerobic cob(II)yrinate a,c-diamide. Figure S8. Biosynthetic pathway of cobalamin from cob(II)yrinate a,c-diamide. Figure S9. Rarefaction curves representing the relationship between the number of taxa (A) and genes (B) as a function of sequencing depth (proportion of reads sampled). Rarefaction analysis showed that sequencing saturation was achieved in both cases. Table S1. Metabolism of tetrapyrrole cofactors: list of enzymes involved in the biosynthesis of siroheme, heme, and cobalamin [73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93]. Table S2. List of protein-encoding genes detected in the metagenome of bacterial community involved in biosynthetic pathways under study. Table S3. List of protein sequences detected in the metaproteome of bacterial community involved in biosynthetic pathways under study.

Author Contributions

R.S. analyzed data and wrote the manuscript; T.K. performed bioinformatics analyzes and consultation on the description of the results; R.M. invented the idea of research; conceived the idea of the manuscript, analyzed data and wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by a research grant from the National Science Center (Poland) (2012/07/B/NZ8/01904).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in Supplementary Materials.

Acknowledgments

The authors acknowledge KGHM Polska Miedź S.A. for help in collecting the samples for this study.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Not available.

References

  1. Barton, L.L.; Fardeau, M.; Fauque, G.D. Hydrogen sulfide: A toxic gas produced by dissimilatory sulfate and sulfur reduction and consumed by microbial oxidation. In The Metal-Driven Biogeochemistry of Gaseous Compounds in the Environment; Springer: New York, NY, USA, 2014; pp. 237–277. [Google Scholar]
  2. Roth, J.R.; Lawrence, J.G.; Bobik, T.A. Cobalamin (coenzyme B12): Synthesis and biological significance. Annu. Rev. Microbiol. 1996, 50, 137–181. [Google Scholar] [CrossRef] [Green Version]
  3. Romine, M.F.; Rodionov, D.A.; Maezato, Y.; Anderson, L.N.; Nandhikonda, P.; Rodionova, I.A.; Wright, A.T. Elucidation of roles for vitamin B12 in regulation of folate, ubiquinone, and methionine metabolism. Proc. Natl. Acad. Sci. USA 2017, 114, E1205–E1214. [Google Scholar] [CrossRef] [Green Version]
  4. Banerjee, R. Chemistry and Biochemistry of B12; John Wiley & Sons: Hoboken, NJ, USA, 1999. [Google Scholar]
  5. Janssen, D.B.; Oppentocht, J.E.; Poelarends, G.J. Microbial dehalogenation. Curr. Opin. Biotechnol. 2001, 12, 254–258. [Google Scholar] [CrossRef] [Green Version]
  6. Rodionov, D.A.; Vitreschak, A.G.; Mironov, A.A.; Gelfand, M.S. Comparative genomics of the vitamin B12 metabolism and regulation in prokaryotes. J. Biol. Chem. 2003, 278, 41148–41159. [Google Scholar] [CrossRef] [Green Version]
  7. Stephen, W.; Ragdale, G. Chapter 6. Biochemistry of Methyl-Coenzyme M Reductase: The Nickel Metalloenzyme that Catalyzes the Final Step in Synthesis and the First Step in Anaerobic Oxidation of the Greenhouse Gas Methane. Met. Ions Life Sci. 2014, 14, 125–145. [Google Scholar]
  8. Eckardt, C.B.; Wolf, M.; Maxwell, J.R. Iron porphyrins in the Permian Kupferschiefer of the lower Rhine Basin, NW Germany. Org. Geochem. 1989, 14, 659–666. [Google Scholar] [CrossRef]
  9. Wolf, M.; David, P.; Eckardt, C.B.; Hagemann, H.W.; Püttmann, W. Facies and rank of the Permian Kupferschiefer from the Lower Rhine Basin and NW Germany. Int. J. Coal Geol. 1989, 14, 119–136. [Google Scholar] [CrossRef]
  10. Czechowski, F. Metalloporphyrin composition and a model for the early diagenetic mineralization of the Permian Kupferschiefer, SW Poland. In Organic Matter and Mineralisation: Thermal Alteration, Hydrocarbon Generation and Role in Metallogenesis; Springer: Dordrecht, The Netherlands, 2000; pp. 243–259. [Google Scholar]
  11. Szubert, A.; Sadowski, Z.; Gros, C.P.; Barbe, J.M.; Guilard, R. Identification of metalloporphyrins extracted from the copper bearing black shale of Fore Sudetic Monocline (Poland). Miner. Eng. 2006, 19, 1212–1215. [Google Scholar] [CrossRef]
  12. Stasiuk, R.; Matlakowska, R. Postdiagenetic bacterial transformation of nickel and vanadyl sedimentary porphyrins of organic-rich shale rock (Fore-Sudetic Monocline, Poland). Front. Microbiol. 2021. [Google Scholar] [CrossRef]
  13. Grice, K.; Schaeffer, P.; Schwark, L.; Maxwell, J.R. Changes in palaeoenvironmental conditions during deposition of the Permian Kupferschiefer (Lower Rhine Basin, northwest Germany) inferred from molecular and isotopic compositions of biomarker components. Org. Geochem. 1997, 26, 677–690. [Google Scholar] [CrossRef]
  14. Tait, G.H. Aminolaevulinate synthetase of Micrococcus denitrificans. Purification and properties of the enzyme, and the effects of growth conditions on the enzyme activity in cells. Biochemistry 1973, 131, 389–403. [Google Scholar] [CrossRef] [Green Version]
  15. Battersby, A.R. Tetrapyrroles: The pigments of life. Nat. Prod. Rep. 2000, 17, 507–526. [Google Scholar] [CrossRef]
  16. Schauer, S.; Chaturvedi, S.; Randau, L.; Moser, J.; Kitabatake, M.; Lorenz, S.; Jahn, D. Escherichia coli glutamyl-tRNA reductase: Trapping the thioester intermediate. J. Biol. Chem. 2002, 277, 48657–48663. [Google Scholar] [CrossRef] [Green Version]
  17. Schubert, H.L.; Raux, E.; Brindley, A.A.; Leech, H.K.; Wilson, K.S.; Hill, C.P.; Warren, M.J. The structure of Saccharomyces cerevisiae Met8p, a bifunctional dehydrogenase and ferrochelatase. EMBO Rep. 2002, 21, 2068–2075. [Google Scholar] [CrossRef] [Green Version]
  18. Lobo, S.A.; Brindley, A.; Warren, M.J.; Saraiva, L.M. Functional characterization of the early steps of tetrapyrrole biosynthesis and modification in Desulfovibrio vulgaris Hildenborough. Biochemistry 2009, 420, 317–326. [Google Scholar] [CrossRef]
  19. Johansson, P.; Hederstedt, L. Organization of genes for tetrapyrrole biosynthesis in gram-positive bacteria. Microbiology 1999, 145, 529–538. [Google Scholar] [CrossRef] [Green Version]
  20. Layer, G.; Verfu, K.; Mahlitz, E.; Jahn, D. Oxygen-independent coproporphyrinogen-III oxidase HemN from Escherichia coli. J. Biol. Chem 2002, 277, 34136–34142. [Google Scholar] [CrossRef] [Green Version]
  21. Dailey, H.A.; Dailey, T.A.; Gerdes, S.; Jahn, D.; Jahn, M.; O’Brian, M.R.; Warren, M.J. Prokaryotic heme biosynthesis: Multiple pathways to a common essential product. Microbiol. Mol. Biol. Rev. 2017, 81, e00048-16. [Google Scholar] [CrossRef] [Green Version]
  22. Dailey, H.A. Enzymes of heme biosynthesis. J. Biol. Inorg. Chem. 1997, 2, 411–417. [Google Scholar] [CrossRef]
  23. Qin, X.; Sun, L.; Wen, X.; Yang, X.; Tan, Y.; Jin, H.; Shen, Y. Structural insight into unique properties of protoporphyrinogen oxidase from Bacillus subtilis. J. Struct. Biol. 2010, 170, 76–82. [Google Scholar] [CrossRef] [PubMed]
  24. Bali, S.; Lawrence, A.D.; Lobo, S.A.; Saraiva, L.M.; Golding, B.T.; Palmer, D.J.; Warren, M.J. Molecular hijacking of siroheme for the synthesis of heme and d1 heme. Proc. Natl. Acad. Sci. USA 2011, 108, 18260–18265. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Adamczack, J.; Hoffmann, M.; Papke, U.; Haufschildt, K.; Nicke, T.; Bröring, M.; Layer, G. NirN protein from Pseudomonas aeruginosa is a novel electron-bifurcating dehydrogenase catalyzing the last step of heme d1 biosynthesis. J. Biol. Chem. 2014, 289, 30753–30762. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Kühner, M.; Haufschildt, K.; Neumann, A.; Storbeck, S.; Streif, J.; Layer, G. The alternative route to heme in the methanogenic archaeon Methanosarcina barkeri. Archaea 2014, 327–340. [Google Scholar]
  27. Celis, A.I.; Gauss, G.H.; Streit, B.R.; Shisler, K.; Moraski, G.C.; Rodgers, K.R.; DuBois, J.L. Structure-based mechanism for oxidative decarboxylation reactions mediated by amino acids and heme propionates in coproheme decarboxylase (HemQ). J. Am. Chem. Soc. 2017, 139, 1900–1911. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Warren, M.J.; Raux, E.; Schubert, H.L.; Escalante-Semerena, J.C. The biosynthesis of adenosylcobalamin (vitamin B12). Nat. Prod. Rep. 2002, 19, 390–412. [Google Scholar] [CrossRef] [PubMed]
  29. Kajiwara, Y.; Santander, P.J.; Roessner, C.A.; Pérez, L.M.; Scott, A.I. Genetically Engineered Synthesis and Structural Characterization of Cobalt− Precorrin 5A and− 5B, Two New Intermediates on the Anaerobic Pathway to Vitamin B12: Definition of the Roles of the CbiF and CbiG Enzymes. J. Am. Chem. Soc. 2006, 128, 9971–9978. [Google Scholar] [CrossRef] [PubMed]
  30. Moore, S.J.; Biedendieck, R.; Lawrence, A.D.; Deery, E.; Howard, M.J.; Rigby, S.E.; Warren, M.J. Characterization of the enzyme CbiH60 involved in anaerobic ring contraction of the cobalamin (vitamin B12) biosynthetic pathway. J. Biol. Chem. 2013, 288, 297–305. [Google Scholar] [CrossRef] [Green Version]
  31. Lobo, S.A.; Videira, M.A.; Pacheco, I.; Wass, M.N.; Warren, M.J.; Teixeira, M.; Saraiva, L.M. Desulfovibrio vulgaris CbiKP cobaltochelatase: Evolution of a haem binding protein orchestrated by the incorporation of two histidine residues. Environ. Microbiol. 2017, 19, 106–118. [Google Scholar] [CrossRef] [PubMed]
  32. Mera, P.E.; Escalante-Semerena, J.C. Dihydroflavin-driven adenosylation of 4-coordinate Co (II) corrinoids: Are cobalamin reductases enzymes or electron transfer proteins? J. Biol. Chem. 2010, 285, 2911–2917. [Google Scholar] [CrossRef] [Green Version]
  33. Spencer, J.B.; Stolowich, N.J.; Roessner, C.A.; Scott, A.I. The Escherichia coli cysG gene encodes the multifunctional protein, siroheme synthase. FEBS Lett. 1993, 335, 57–60. [Google Scholar] [CrossRef] [Green Version]
  34. Hogle, S.L.; Barbeau, K.A.; Gledhill, M. Heme in the marine environment: From cells to the iron cycle. Metallomics 2014, 6, 1107–1120. [Google Scholar] [CrossRef]
  35. Pan, J.; Zhou, Z.; Béjà, O.; Cai, M.; Yang, Y.; Liu, Y.; Li, M. Genomic and transcriptomic evidence of light-sensing, porphyrin biosynthesis, Calvin-Benson-Bassham cycle, and urea production in Bathyarchaeota. Microbiome 2020, 8, 1–12. [Google Scholar] [CrossRef] [PubMed]
  36. Feng, G.; Xie, T.; Wang, X.; Bai, J.; Tang, L.; Zhao, H.; Zhao, Y. Metagenomic analysis of microbial community and function involved in cd-contaminated soil. BMC Microbiol. 2018, 18, 1–13. [Google Scholar] [CrossRef]
  37. Ghosh, S.; Das, A.P. Metagenomic insights into the microbial diversity in manganese-contaminated mine tailings and their role in biogeochemical cycling of manganese. Sci. Rep. 2018, 8, 1–12. [Google Scholar] [CrossRef]
  38. Doxey, A.C.; Kurtz, D.A.; Lynch, M.D.; Sauder, L.A.; Neufeld, J.D. Aquatic metagenomes implicate Thaumarchaeota in global cobalamin production. ISME J. 2015, 9, 461–471. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Panek, H.; O’Brian, M.R. A whole genome view of prokaryotic haem biosynthesis. Microbiology 2002, 148, 2273–2282. [Google Scholar] [CrossRef] [Green Version]
  40. Glanville, D.G.; Mullineaux-Sanders, C.; Corcoran, C.J.; Burger, B.T.; Imam, S.; Donohue, T.J.; Ulijasz, A.T. A High-Throughput Method for Identifying Novel Genes That Influence Metabolic Pathways Reveals New Iron and Heme Regulation in Pseudomonas aeruginosa. Msystems 2021, 6, e00933-20. [Google Scholar] [CrossRef] [PubMed]
  41. Lu, X.; Heal, K.R.; Ingalls, A.E.; Doxey, A.C.; Neufeld, J.D. Metagenomic and chemical characterization of soil cobalamin production. ISME J. 2020, 14, 53–66. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Balabanova, L.; Averianova, L.; Marchenok, M.; Son, O.; Tekutyeva, L. Microbial and genetic resources for cobalamin (vitamin B12) biosynthesis: From ecosystems to industrial biotechnology. Int. J. Mol. Sci. 2021, 22, 4522. [Google Scholar] [CrossRef]
  43. Boynton, T.O.; Gerdes, S.; Craven, S.H.; Neidle, E.L.; Phillips, J.D.; Dailey, H.A. Discovery of a gene involved in a third bacterial protoporphyrinogen oxidase activity through comparative genomic analysis and functional complementation. Appl. Environ. Microbiol. 2011, 77, 4795–4801. [Google Scholar] [CrossRef] [Green Version]
  44. Cavallaro, G.; Decaria, L.; Rosato, A. Genome-based analysis of heme biosynthesis and uptake in prokaryotic systems. J. Proteome Res. 2008, 7, 4946–4954. [Google Scholar] [CrossRef] [PubMed]
  45. Hopkinson, B.M.; Roe, K.L.; Barbeau, K.A. Heme uptake by Microscilla marina and evidence for heme uptake systems in the genomes of diverse marine bacteria. Appl. Environ. Microbiol. 2008, 74, 6263–6270. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Dent, A.T.; Wilks, A. Contributions of the heme coordinating ligands of the Pseudomonas aeruginosa outer membrane receptor HasR to extracellular heme sensing and transport. J. Biol. Chem. 2020, 295, 10456–10467. [Google Scholar] [CrossRef]
  47. Datta, R.; Kelkar, A.; Baraniya, D.; Molaei, A.; Moulick, A.; Meena, R.S.; Formanek, P. Enzymatic degradation of lignin in soil: A review. Sustainability 2017, 9, 1163. [Google Scholar] [CrossRef] [Green Version]
  48. Wang, M.; Nie, Y.; Wu, X.L. Extracellular heme recycling and sharing across species by novel mycomembrane vesicles of a Gram-positive bacterium. ISME J. 2020, 15, 605–617. [Google Scholar] [CrossRef] [PubMed]
  49. Caldwell, D.R.; White, D.C.; Bryant, M.P.; Doetsch, R.N. Specificity of the heme requirement for growth of Bacteroides ruminicola. J. Bacteriol. Res. 1965, 90, 1645–1654. [Google Scholar] [CrossRef] [Green Version]
  50. Duwat, P.; Sourice, S.; Cesselin, B.; Lamberet, G.; Vido, K.; Gaudu, P. Respiration capacity of the fermenting bacterium Lactococcus lactis and its positive effects on growth and survival. J. Bacteriol. Res. 2001, 183, 4509–4516. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. White, D.C.; Granick, S. Hemin biosynthesis in Haemophilus. J. Bacteriol. Res. 1963, 85, 842–850. [Google Scholar] [CrossRef] [Green Version]
  52. Fournier, C.; Smith, A.; Delepelaire, P. Haem release from haemopexin by HxuA allows Haemophilus influenzae to escape host nutritional immunity. Mol. Microbiol. 2011, 80, 133–148. [Google Scholar] [CrossRef]
  53. Jacobs, N.J.; Jacobs, J.M.; Brent, P. Formation of protoporphyrin from coproporphyrinogen in extracts of various bacteria. J. Bacteriol. Res. 1970, 102, 398–403. [Google Scholar] [CrossRef] [Green Version]
  54. Lawrance, A.D.; Nemoto-Smith, E.; Deery, E.; Baker, J.A.; Schroeder, S.; Brown, D.G.; Warren, M.J. Construction of fluorescent analogs to follow the uptake and distribution of cobalamin (vitamin B12) in bacteria, worms, and plants. Cell Chem. Biol. 2018, 25, 941–951.e6. [Google Scholar] [CrossRef]
  55. Bernal, V.; Castaño-Cerezo, S.; Gallego-Jara, J.; Écija-Conesa, A.; de Diego, T.; Iborra, J.L.; Cánovas, M. Regulation of bacterial physiology by lysine acetylation of proteins. New Biotechnol. 2014, 31, 586–595. [Google Scholar] [CrossRef] [PubMed]
  56. Castaño-Cerezo, S.; Bernal, V.; Post, H.; Fuhrer, T.; Cappadona, S.; Sánchez-Díaz, N.C.; Cánovas, M. Protein acetylation affects acetate metabolism, motility and acid stress response in Escherichia coli. Mol. Syst. Biol. 2014, 10, 762. [Google Scholar] [CrossRef] [PubMed]
  57. Hentchel, K.L.; Escalante-Semerena, J.C. Acylation of biomolecules in prokaryotes: A widespread strategy for the control of biological function and metabolic stress. Microbiol. Mol. Biol. Rev. 2015, 79, 321–346. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Demopoulos, B.J.; Anderson, H.J.; Loader, C.E.; Faber, K. Pyrrole chemistry. XXVI. A synthesis of porphobilinogen from pyrrole. Can. J. Chem. 1983, 61, 2415–2422. [Google Scholar] [CrossRef]
  59. Robinson, R. Structural Relations of Natural Products; Clarendon Press: Oxford, UK, 1955. [Google Scholar]
  60. Heal, K.R.; Qin, W.; Ribalet, F.; Bertagnolli, A.D.; Coyote-Maestas, W.; Hmelo, L.R.; Ingalls, A.E. Two distinct pools of B12 analogs reveal community interdependencies in the ocean. Proc. Natl. Acad. Sci. USA 2017, 114, 364–369. [Google Scholar] [CrossRef] [Green Version]
  61. Zhou, J.; Bruns, M.A.; Tiedje, J.M. DNA recovery from soils of diverse composition. Appl. Environ. Microbiol. 1996, 62, 316–322. [Google Scholar] [CrossRef] [Green Version]
  62. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. J. Bioinform. 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  63. Agarwala, R.; Rotmistrovsky, K. BMTagger: Best Match Tagger for Removing Human Reads from Metagenomics Datasets. 2011. Available online: Ftp://ftp.ncbi.nlm.nih.gov/pub/agarwala/bmtagger/ (accessed on 21 February 2014).
  64. Wood, D.E.; Lu, J.; Langmead, B. mproved metagenomic analysis with Kraken 2. Genome Biol. 2019, 20, 1–13. [Google Scholar] [CrossRef] [Green Version]
  65. Prjibelski, A.; Antipov, D.; Meleshko, D.; Lapidus, A.; Korobeynikov, A. Using SPAdes De Novo Assembler. Curr. Protoc. Bioinform. 2020, 70, e102. [Google Scholar] [CrossRef]
  66. Seemann, T. Prokka: Rapid prokaryotic genome annotation. J. Bioinform. 2014, 30, 2068–2069. [Google Scholar] [CrossRef]
  67. Aramaki, T.; Blanc-Mathieu, R.; Endo, H.; Ohkubo, K.; Kanehisa, M.; Goto, S.; Ogata, H. KofamKOALA: KEGG ortholog assignment based on profile HMM and adaptive score threshold. J. Bioinform. 2020, 36, 2251–2252. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Bushnell, B. BBMap: A Fast, Accurate, Splice-Aware Aligner; Lawrence Berkeley National Lab: Berkeley, CA, USA, 2014. [Google Scholar]
  69. Wagner, G.P.; Kin, K.; Lynch, V.J. Measurement of mRNA abundance using RNA-seq data: RPKM measure is inconsistent among samples. Theory Biosci. 2012, 131, 281–285. [Google Scholar] [CrossRef] [PubMed]
  70. Lakin, S.M.; Dean, C.; Noyes, N.R.; Dettenwanger, A.; Ross, A.S.; Doster, E.; Boucher, C. MEGARes: An antimicrobial resistance database for high throughput sequencing. Nucleic Acids Res. 2017, 45, D574–D580. [Google Scholar] [CrossRef] [PubMed]
  71. Ram, R.J.; VerBerkmoes, N.C.; Thelen, M.P.; Tyson, G.W.; Baker, B.J.; Blake, R.C., II; Shah, M.; Hettich, R.L.; Banfield, J.F. Community proteomics of a natural microbial biofilm. Science 2005, 308, 1915–1920. [Google Scholar] [CrossRef]
  72. Kanehisa, M.; Sato, Y.; Morishima, K. BlastKOALA and GhostKOALA: KEGG tools for functional characterization of genome and metagenome sequences. J. Mol. Biol. 2016, 428, 726–731. [Google Scholar] [CrossRef] [Green Version]
  73. Ravel, J.M.; Wang, S.F.; Heinemeyer, C.; Shive, W. Glutamyl and glutaminyl ribonucleic acid synthetases of Escherichia coli W: Separation, properties, and stimulation of adenosine triphosphate-pyrophosphate exchange by acceptor ribonucleic acid. J. Biol. Chem. 1965, 240, 432–438. [Google Scholar] [CrossRef]
  74. Beale, S.I. Biosynthesis of 5-aminolevulinic acid. In Chlorophylls and Bacteriochlorophylls; Springer: Dordrecht, The Netherlands, 2006; pp. 147–158. [Google Scholar]
  75. Friedman, J.; Lad, L.; Li, H.; Wilks, A.; Poulos, T.L. Structural basis for novel δ-regioselective heme oxygenation in the opportunistic pathogen Pseudomonas aeruginosa. Biochemistry. 2004, 43, 5239–5245. [Google Scholar] [CrossRef]
  76. Tian, B.X.; Erdtman, E.; Eriksson, L.A. Catalytic mechanism of porphobilinogen synthase: The chemical step revisited by QM/MM calculations. J. Phys. Chem. B 2012, 116, 12105–12112. [Google Scholar] [CrossRef]
  77. Battersby, A.R.; Fookes, C.J.; Matcham, G.W.; McDonald, E. Biosynthesis of the pigments of life: Formation of the macrocycle. Nature 1980, 285, 17–21. [Google Scholar] [CrossRef]
  78. Kohno, H.; Furukawa, T.; Yoshinaga, T.; Tokunaga, R.; Taketani, S. Coproporphyrinogen oxidase. Purification, molecular cloning, and induction of mRNA during erythroid differentiation. J. Biol. Chem. 1993, 268, 21359–21363. [Google Scholar] [CrossRef]
  79. Boynton, T.O.; Daugherty, L.E.; Dailey, T.A.; Dailey, H.A. Identification of Escherichia coli HemG as a novel, menadione-dependent flavodoxin with protoporphyrinogen oxidase activity. Biochemistry 2009, 48, 6705–6711. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Möbius, K.; Arias-Cartin, R.; Breckau, D.; Hännig, A.L.; Riedmann, K.; Biedendieck, R.; Schröder, S.; Becher, D.; Magalon, A.; Moser, J.; et al. Heme biosynthesis is coupled to electron transport chains for energy generation. Proc. Natl. Acad. Sci. USA 2010, 107, 10436–10441. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Skotnicová, P.; Sobotka, R.; Shepherd, M.; Hájek, J.; Hrouzek, P.; Tichý, M. The cyanobacterial protoporphyrinogen oxidase HemJ is a new b-type heme protein functionally coupled with coproporphyrinogen III oxidase. J. Biol. Chem. 2018, 293, 12394–12404. [Google Scholar] [CrossRef] [Green Version]
  82. Al-Karadaghi, S.; Hansson, M.; Nikonov, S.; Jönsson, B.; Hederstedt, L. Crystal structure of ferrochelatase: The terminal enzyme in heme biosynthesis. Structure 1997, 5, 1501–1510. [Google Scholar] [CrossRef]
  83. Mogi, T. Over-expression and characterization of Bacillus subtilis heme O synthase. J. Biochem. 2009, 145, 669–675. [Google Scholar] [CrossRef]
  84. Niwa, S.; Takeda, K.; Kosugi, M.; Tsutsumi, E.; Mogi, T.; Miki, K. Crystal structure of heme A synthase from Bacillus subtilis. Proc. Natl. Acad. Sci. USA 2018, 115, 11953–11957. [Google Scholar] [CrossRef] [Green Version]
  85. Spencer, P.; Stolowich, N.J.; Sumner, L.W.; Scott, A.I. Definition of the redox states of cobalt-precorrinoids: Investigation of the substrate and redox specificity of CbiL from Salmonella typhimurium. Biochemistry 1998, 37, 14917–14927. [Google Scholar] [CrossRef]
  86. Debussche, L.; Thibaut, D.; Cameron, B.; Crouzet, J.; Blanche, F. Biosynthesis of the corrin macrocycle of coenzyme B12 in Pseudomonas denitrificans. J. Bacteriol 1993, 175, 7430–7440. [Google Scholar] [CrossRef] [Green Version]
  87. Moore, S.J.; Lawrence, A.D.; Biedendieck, R.; Deery, E.; Frank, S.; Howard, M.J.; Rigby, S.E.; Warren, M.J. Elucidation of the anaerobic pathway for the corrin component of cobalamin (vitamin B12). Proc. Natl. Acad. Sci. USA 2013, 110, 14906–14911. [Google Scholar] [CrossRef] [Green Version]
  88. Fresquet, V.; Williams, L.; Raushel, F.M. Mechanism of cobyrinic acid a, c-diamide synthetase from Salmonella typhimurium LT2. Biochemistry 2004, 43, 10619–10627. [Google Scholar] [CrossRef]
  89. Roessner, C.A.; Warren, M.J.; Santander, P.J.; Atshaves, B.P.; Ozaki, S.I.; Stolowich, N.J.; Iida, K.; Scott, A.I. Expression of 9 Salmonella typhimurium enzymes for cobinamide synthesis Identification of the 11-methyl and 20-methyl transferases of corrin biosynthesis. FEBS Lett. 1992, 301, 73–78. [Google Scholar] [CrossRef] [Green Version]
  90. Deery, E.; Schroeder, S.; Lawrence, A.D.; Taylor, S.L.; Seyedarabi, A.; Waterman, J.; Wilson, K.S.; Brown, D.; Geeves, M.A.; Howard, M.J.; et al. An enzyme-trap approach allows isolation of intermediates in cobalamin biosynthesis. Nat. Chem. Biol. 2012, 8, 933–940. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  91. Roth, J.R.; Lawrence, J.G.; Rubenfield, M.; Kieffer-Higgins, S.; Church, G.M. Characterization of the cobalamin (vitamin B12) biosynthetic genes of Salmonella typhimurium. J. Bacteriol. Res. 1993, 175, 3303–3316. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  92. Cheong, C.G.; Bauer, C.B.; Brushaber, K.R.; Escalante-Semerena, J.C.; Rayment, I. Three-dimensional structure of the L-threonine-O-3-phosphate decarboxylase (CobD) enzyme from Salmonella enterica. Biochemistry 2002, 41, 4798–4808. [Google Scholar] [CrossRef] [PubMed]
  93. Watkins, H.A.; Baker, E.N. Structural and functional characterization of an RNase HI domain from the bifunctional protein Rv2228c from Mycobacterium tuberculosis. J. Bacteriol. Res. 2010, 192, 2878–2886. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Biosynthetic pathways of tetrapyrrole cofactors (protoheme, siroheme, cobalamin and F430) investigated in this study, including the most important precursors (glycine; L-glutamate, uroporphyrinogen III) and intermediates (coproporphyrin III, protoporphyrin IX, precorrin-2, cob(II)yrinate a,c diamide). A detailed description of the pathways is included in the manuscript and in Supplementary Materials (Table S1 and Figures S1–S8). The colour of the arrows indicates whether the pathways have been fully detected (black), partially detected (white) or undetected (grey).
Figure 1. Biosynthetic pathways of tetrapyrrole cofactors (protoheme, siroheme, cobalamin and F430) investigated in this study, including the most important precursors (glycine; L-glutamate, uroporphyrinogen III) and intermediates (coproporphyrin III, protoporphyrin IX, precorrin-2, cob(II)yrinate a,c diamide). A detailed description of the pathways is included in the manuscript and in Supplementary Materials (Table S1 and Figures S1–S8). The colour of the arrows indicates whether the pathways have been fully detected (black), partially detected (white) or undetected (grey).
Molecules 26 06746 g001
Figure 2. Bacterial community inhabiting shale rock: macroscopic micrograph (A) and scanning electron micrographs (B,C). Diversity of dominating (>1%) phyla (D), classes (E), and genera (F) identified in the bacterial community.
Figure 2. Bacterial community inhabiting shale rock: macroscopic micrograph (A) and scanning electron micrographs (B,C). Diversity of dominating (>1%) phyla (D), classes (E), and genera (F) identified in the bacterial community.
Molecules 26 06746 g002
Figure 3. General characteristics of the metagenome and metaproteome of bacterial community. (A) The number of unique reads and protein sequences in the four main functional categories; taxonomic origin of the detected protein-encoding genes (B) and protein sequences (C) involved in biosynthesis of tetrapyrrole cofactors.
Figure 3. General characteristics of the metagenome and metaproteome of bacterial community. (A) The number of unique reads and protein sequences in the four main functional categories; taxonomic origin of the detected protein-encoding genes (B) and protein sequences (C) involved in biosynthesis of tetrapyrrole cofactors.
Molecules 26 06746 g003
Figure 4. Detection of protein-encoding genes in the metagenome and number of protein sequences identified in the metaproteome of bacterial community involved in the biosynthesis of tetrapyrrole cofactors.* List of enzymes involved in anaerobic biosynthesis of cob(II)yrinate diamide detected neither in metagenome nor in metaproteome of bacterial community. ** List of enzymes involved in aerobic biosynthesis of cob(II)yrinate diamide detected neither in metagenome nor in metaproteome of bacterial community.
Figure 4. Detection of protein-encoding genes in the metagenome and number of protein sequences identified in the metaproteome of bacterial community involved in the biosynthesis of tetrapyrrole cofactors.* List of enzymes involved in anaerobic biosynthesis of cob(II)yrinate diamide detected neither in metagenome nor in metaproteome of bacterial community. ** List of enzymes involved in aerobic biosynthesis of cob(II)yrinate diamide detected neither in metagenome nor in metaproteome of bacterial community.
Molecules 26 06746 g004
Table 1. Primary and modified geoporphyrins and pyrrole-containing organic compounds detected in the Kupferschiefer shale rock.
Table 1. Primary and modified geoporphyrins and pyrrole-containing organic compounds detected in the Kupferschiefer shale rock.
Organic CompoundNameReference
Primary
geoporphyrins
Iron aetioporphyrins (octaethyl porphyrins),[8]
Iron cycloalkanoporphyrins
Iron di-cycloalkanoporphyrins
Iron benz-cycloalkanoporphyrins
Vanadyl-cycloalkano-porphyrins [9]
Etio and DPEP iron porphyrins [10]
Etio,DPEP benzo-etio and benzo-DPEP vanadyl porphyrins
Vanadyl porphyrins of series etio/DPEP [11]
Octaethyl nickel porphyrin [12]
Octaethyl vanadyl porphyrin
Meso-tetraphenyl vanadyl porphyrin
Meso-tetraphenyl nickel porphyrin
Protoporphyrin IX cobalt
Modified geoporphyrinsDiphenyl vanadyl porphyrin[12]
Tetraethyl vanadyl porphyrin
Vanadyl porphyrin
Tetraethyl nickel porphyrin
Nickel porphyrin
Organic compounds containing 3 pyrrole rings3-[2-[[3-(2-Carboxyethyl)-5-[(3.4-dimethyl-5-oxopyrrol-2-ylidene)methyl]-4-methyl-1H-pyrrol-2-yl]methylidene]-4-methyl-5-oxopyrrol-3-yl]propanoic acid[12]
3-[(5Z)-5-[[4-Ethenyl-5-[(Z)-(4-ethenyl-3-methyl-5-oxopyrrol-2-ylidene)methyl]-3-methyl-1H-pyrrol-2-yl]methylidene]-4-methyl-2-oxopyrrol-3-yl]propanoate
Organic compounds containing 2 pyrrole rings3,3′-Bipyrrole
2,2′-Bipyrrole
1,1′-Bipyrrole-2,2′,5,5′-tetraone
3,3′,4,4′-Tetramethyl-1H,1′H-2,2′-bipyrrole-5,5′-dicarboxylic acid
Organic compounds containing 1 pyrrole ring1H-Pyrrole
1H-Pyrrole-2-carboxylic acid
2H-Pyrrol-2-one, 5-[[2-[(4-aminophenyl)methylene]-3,4-dimethyll]methylene]-3-ethyl-1,5-dihydro-4-methyl
Indole acetic acid
1H-Indole-2-carboxylic acid
Indole carbaldehyde
1H-pyrrole-2,5-diones [13]
Table 2. Enzymes containing tetrapyrrole cofactors detected in the studied metaproteome of bacterial community.
Table 2. Enzymes containing tetrapyrrole cofactors detected in the studied metaproteome of bacterial community.
EnzymeAccessionScoreSeq(sig)emPAIGenus/Species
Heme-Containing Cytochromes
cbb3-type cytochrome c oxidase subunit iigi|6532510206210.19Dechloromonas agitata
cb-type cytochrome c oxidase ccoo subunitgi|45192095510.12Magnetospirillum magnetotacticum
cytochrome b6gi|49924589712820.13Geobacter sulfurreducens
cytochrome b6 gi|49392452211720.13Legionella drancourtii
cytochrome bd-type quinol oxidase, subunit 1gi|4939753886710.07Desulfovibrio magneticus
cytochrome cgi|5040019957120.1Azospira oryzae
cytochrome c oxidase cbb3-type, subunit iiigi|3309492015610.77Pseudomonas syringa. 1704B
cytochrome c1gi|2685844775810.11Neisseria gonorrhoeae PID18
cytochrome cbb3gi|4996303455510.27Thiobacillus denitrificans
cytochrome d ubiquinol oxidasegi|4981850395510.06Lactobacillus acidipiscis
cytochrome o ubiquinol oxidasegi|6571984194010.12Aeromonas caviae
cytochrome p450gi|5031905892010.1Frankia sp. EuI1c
cytochrome p450gi|6644332994410.07Streptomyces sp. NRRL F-5140
cytochrome soxagi|4996303424410.11Thiobacillus denitrificans
cytochrome soxagi|5173330697410.11Thiobacillus thioparu
multispecies: cytochrome cgi|4949621026810.25Sphingobium sp.
flavocytochrome c sulfide dehydrogenasegi|5190120583910.07Methylotenera sp.
succinate dehydrogenase cytochrome b-556 subunitgi|2546728726610.39Neisseria meningitidis alpha275
thiosulfate reductase cytochrome b subunitgi|4887137806110.15Myxococcus sp.
Heme-Containing Catalases
catalasegi|491324047515110.84Acinetobacter sp.CIP 53.82
catalasegi|39278908610.06Desulfovibrio vulgaris
catalasegi|5204016010.06Haemophilus influenzae
catalasegi|50025165910320.13Pseudomonas stutzeri
catalasegi|5049381316710.09Synechococcus sp. PCC 6312
Heme-Containing Peroxidases
hydroperoxidasegi|6475314744610.04Shewanella marina
hydroperoxidase IIgi|4893756708420.09Pseudomonas stutzeri
hydroperoxidase IIgi|51581522810520.09Pseudomonas stutzeri
peroxidasegi|4911421205810.14Nitrococcus mobilis
Sirroheme-Containing Enzymes
nitrite reductase (NAD(P)H) large subunitgi|1598829756210.04Hydrogenivirga sp. 128-5-R1-1
nitrite reductasegi|4911299945810.08Streptomyces ghanaensis
nitrite reductasegi|233929876510.14uncultured bacterium
sulfite reductasegi|52106509568960.47Thiothrix disciformis
sulfite reductasegi|6550416505310.08Thiothrix lacustris
sulfite reductasegi|488797739207671.18Thiothrix nivea
sulfite reductasegi|48879774025130.27Thiothrix nivea
dissimilatory sulfite reductase alpha subunitgi|305254975210.08uncultured sulfate-reducing bacterium
reverse-type dissimilatory siroheme sulfite reductase subunit Agi|162072844218771.3Thiothrix nivea DSM 5205
reverse-type dissimilatory sulfite reductase (rDSR), alpha subunit (DsrA)gi|3857636984810.07uncultured bacterium 172H5
Cobalamin-Containing Enzymes
methylcrotonoyl-CoA carboxylasegi|6486181953610.06Niabella aurantiaca
methylmalonyl-CoA carboxyltransferasegi|5876411914010.06Skermanella stibiiresistens SB22
methionyl-tRNA synthetasegi|5164100086810.04Erwinia toletana
methionyl-tRNA synthetasegi|5884762333510.04Lactobacillus composti JCM 14202
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Stasiuk, R.; Krucoń, T.; Matlakowska, R. Biosynthesis of Tetrapyrrole Cofactors by Bacterial Community Inhabiting Porphyrine-Containing Shale Rock (Fore-Sudetic Monocline). Molecules 2021, 26, 6746. https://doi.org/10.3390/molecules26216746

AMA Style

Stasiuk R, Krucoń T, Matlakowska R. Biosynthesis of Tetrapyrrole Cofactors by Bacterial Community Inhabiting Porphyrine-Containing Shale Rock (Fore-Sudetic Monocline). Molecules. 2021; 26(21):6746. https://doi.org/10.3390/molecules26216746

Chicago/Turabian Style

Stasiuk, Robert, Tomasz Krucoń, and Renata Matlakowska. 2021. "Biosynthesis of Tetrapyrrole Cofactors by Bacterial Community Inhabiting Porphyrine-Containing Shale Rock (Fore-Sudetic Monocline)" Molecules 26, no. 21: 6746. https://doi.org/10.3390/molecules26216746

Article Metrics

Back to TopTop