- Institute of Process Engineering in Life Sciences, II, Technical Biology, Karlsruhe Institute of Technology, Karlsruhe, Germany
After major mass extinction events, ancient plants and terrestrial vertebrates were faced with various challenges, especially ultraviolet (UV) light. These stresses probably resulted in changes in the biosynthetic pathways, which employed the MIO (3,5-dihydro-5-methylidene-4H-imidazole-4-one)-dependent enzymes (ammonia-lyase and aminomutase), leading to enhanced accumulation of metabolites for defense against UV radiation, pathogens, and microorganisms. Up to now, the origin and evolution of genes from this superfamily have not been extensively studied. In this report, we perform an analysis of the phylogenetic relations between the members of the aromatic amino acid MIO-dependent enzymes (AAM), which demonstrate that they most probably have a common evolutionary origin from ancient bacteria. In early soil environments, numerous bacterial species with tyrosine ammonia-lyase genes (TAL; EC 4.3.1.23) developed tyrosine aminomutase (TAM; EC 5.4.3.6) activity as a side reaction for competing with their neighbors in the community. These genes also evolved into other TAL-like enzymes, such as histidine ammonia-lyase (HAL, EC 4.3.1.3) and phenylalanine ammonia-lyase (PAL; EC 4.3.1.24), in different bacterial species for metabolite production and accumulation for adaptation to adverse terrestrial environmental conditions. On the other hand, the existence of phenylalanine aminomutase (PAM; EC 5.4.3.10) and phenylalanine/tyrosine ammonia-lyase (PTAL; EC 4.3.1.25) strongly indicates the horizontal gene transfer (HGT) between bacteria, fungi, and plants in symbiotic association after acquiring the PAL gene from their ancestor.
Introduction
MIO Prosthetic Group
Usually, the side chains of proteogenic amino acids act as nucleophiles in enzymatic catalysis. Due to the lack of strongly electrophilic groups, enzymes use metal ions and organic molecules to assist in electrophilic catalysis (Bischoff and Schlüter, 2012; Punekar, 2018). Besides cofactors from the environment, posttranslational modifications/conversions (PTMs) of amino acid side chains that provide a strongly electrophilic center are of the same importance in enzymes. (Müller, 2018). MIO (3,5-dihydro-5-methylidene-4H-imidazole-4-one) is such a catalytic moiety for the elimination of ammonia from arylalanine amino acids, which belong to arylalanine ammonia-lyase and aminomutase. This highly electrophilic moiety is spontaneously folded by an inner amino acid triad (Ala/Ser/Cys/Thr)-Ser-Gly. The glycine amide lone pair attacks the π* orbital of carbonyl in the amino acid, which is located at two positions preceding glycine. This nucleophilic attack is electronically unfavorable; therefore, mechanical compression from neighboring residues and connections with internal water molecules play a vital role by promoting backbone cyclization (Baedeker and Schulz, 2002; Sánchez-Murcia et al., 2016). The formation mechanism of MIO is similar to the chromophore in green (GFP) and red fluorescent protein (RFP) (Figure 1) (Reid and Flynn, 1997; Baedeker and Schulz, 2002).
FIGURE 1. The proposed mechanism for the MIO in HAL and GFP formation by posttranslational modification.
The MIO prosthetic group performs a nucleophilic attack on the substrate arylalanine amino acid, forming a covalent amino-MIO intermediate to enhance its acidity. In the most accepted E1cB mechanism of ammonia elimination in AAM, the intermediate binds the MIO at the benzylic position of the amino acid and is deprotonated by the enzymatic base. Then, the intermediate yields a carbanion intermediate and releases ammonia in the subsequent step. In the alternative mechanism (Friedel–Crafts), this reaction occurs at the aryl side chain of the substrate (Figure 2) (Punekar, 2018).
FIGURE 2. Two proposals for ammonia elimination by R. toruloides PAL. The other residues help to orient and stabilize the carboxyl group of arylalanine amino acids.
Structures of MIO-dependent Enzymes
According to the crystal structure, the members of AAM are homotetramers, comprising four identical active sites that build among the residues of three interlocking monomeric subunits in a nose-to-tail way separately. Each monomer is composed of a rigid central core domain, a globular N-terminal domain, and an elongated C-terminal domain. In the N-terminal region, the MIO prosthetic group is located in a narrow tunnel, which is created by an inner and an outer loop. These two loops cap the tunnel and enclose the active site from solvents in HAL. However, they are more mobile in eukaryotic PAL, which restricts the access of the substrate and influences the mutase and lyase activity in MIO-dependent enzymes. The catalytic tyrosine residue (Tyr110 in Rhodosporidium toruloides PAL) on the inner loop is highly conserved in HAL, TAL, and PAL, which is necessary for catalytic activity since the mutants of this residue to alanine and phenylalanine are completely inactive (Röther et al., 2002; Wybenga et al., 2014). The residues in the carboxyl binding pocket promote the MIO prosthetic group formation and also interact with the substrate, representing essential binding sites in the active center (Wybenga et al., 2014), whereas the hydrophobic binding pocket provides sufficient space for the aromatic ring of the substrates and modulates the substrate specificity (Nagy et al., 2019).
Phenylpropanoid Pathway
Replying to abiotic environmental stress, plants activate the phenylpropanoid pathway to accumulate phenolic secondary metabolites. As a major component of phenylpropanoids, flavonoids provide significant protective effects to plants in response to various unfavorable conditions (drought, heavy metals, salinity, and UV radiations) (Sharma et al., 2019). Photoprotection is the most important functional role of flavonoids. The epidermal flavonoids reduce protein and DNA damage by preventing dimerization of thymine, adsorbing radiations, and scavenging the reactive oxygen species (Kootstra, 1994; Treutter, 2006). Other specific metabolites, monolignols, confer tolerance to plant cell walls against chilling stress. The phenylpropanoid accumulation is regulated by the gene expression of corresponding enzymes in the biosynthesis pathway (Sharma et al., 2019).
PAL is the first regulatory enzyme that transforms l-phenylalanine into trans-cinnamic acid, controlling the carbon flux from the shikimate pathway to phenylpropanoid metabolism. The subsequent enzyme is cinnamic acid 4-hydroxylase (C4H; EC 1.14.13.11), which reduces trans-cinnamic acid, leading to the formation of p-coumaric acid. 4-coumaric acid-CoA ligase (4CL; EC 6.2.1.12) catalyzes the ATP-dependent formation of the p-coumaroyl CoA, which serves as the branch point in phenylpropanoid biosynthesis. The conversions by these initial three enzymes are necessary and form the main pathway skeleton in the higher plants. In general, l-phenylalanine is starting immediately from the shikimate pathway, and in certain monocot species, p-coumaric acid may directly be produced from l-tyrosine through PTAL with bifunctional activity, bypassing the hydroxylation by C4H (Figure 3). In comparison to dicotyledonous plants, these monocots form the cell walls with higher proportions of syringyl (S)-rich lignins, more esterified coumaric acid as well as flavonoid tricin (Barros et al., 2016). This preferential composition suggests the existence of alternative routes with PTAL in the lignin biosynthetic pathway. The methylation of caffeic acid in this parallel pathway employs the caffeate/5-hydroxyferulate 3-O-methyltransferase (COMT), encouraging an efficient conversion step of p-coumaric acid to caffeic acid through the cytosolic coumarate 3-hydroxylase (C3H) (Figure 3). This alternative pathway with PTAL and C3H in monocots bypasses the synthesis of H and G lignins by membrane-bound cytochrome P450, leading to the enhancing of (S)-lignin content efficiently (Barros et al., 2019).
FIGURE 3. The phenylpropanoid pathway in plants for the production of flavonoids and monolignols under abiotic stress conditions.
Histidine Pathway
Histidine catabolism is initiated by HAL, which is a universal enzyme to form trans-urocanic acid and release ammonia in the bacterial pathway. trans-Urocanic acid is transformed into 4-imidazole-5-propionate by urocanate hydratase (EC 4.2.1.49). Subsequently, the ring is cleaved by imidazolonepropionase (EC 3.5.2.7). Depending on the organism species, l-glutamate is generated by various enzymes from N-formimino-l-glutamate. In mammals, the N-formimino-l-glutamate is hydrolyzed by tetrahydrofolate (THF)-dependent glutamate formiminotransferase (EC 2.1.2.5). In bacteria, some genera (Bacillus, Klebsiella, and Salmonella) eliminate the formimino group in one step by formininoglutamase (EC 3.5.3.8), whereas the other genera, such as Pseudomonas, coregulate the hutF and hutG genes and express formimidoylglutamate deiminase (EC 3.5.3.13) and N-formylglutamate deformylase (EC 3.5.1.68) to yield the l-glutamate (Figure 4) (Bender, 2012; Kohlmeier, 2015).
HAL is crucial for growth in children, primarily in the skin and the liver. For example, one disease, histidinemia, results from a deficiency of HAL (Brosnan and Brosnan, 2020). The trans-urocanic acids, as UV-absorbing compounds, accumulate in the stratum corneum of the skin because of the absence of urocanase to eliminate them. Under UV light, trans-urocanic acid is isomerized into cis-urocanic acid until the quantity equation of these two isomers. The latter form probably initiates the immunoregulation under UV exposure (Norval, 2001).
Results
Phylogenetic Analysis
To explore the evolutionary history of AAM, a phylogeny of 268 available biochemical-characterized protein sequences from the SwissProt was reconstructed (Bairoch et al., 2004). As displayed in Figure 5 and Supplementary Figure S1, AMM sequences from various phyla are grouped in two well-supported clusters (bootstrap support: 96%). One group contained all PAL sequences across different phyla and one eukaryotic HAL from Dictyostelium discoideum. The second cluster comprised all bacterial TALs/TAMs and the rest of the HALs. Overall, the phylogeny suggests that the arylalanine amino acids ammonia-lyase and aminomutase likely share common ancestry.
FIGURE 5. Phylogenetic tree of arylalanine amino acid ammonia-lyase and aminomutase homologs retrieved from the Swissprot database. Characterized ammonia-lyases are shown in the square while characterized aminomutases are shown in the circle. The color coding corresponds to the substrates of the enzymes: histidine (yellow), phenylalanine (green), tyrosine (purple). The source organisms in the tree are color-coded in the outer ring according to their origin.
In the first group, the eukaryotic PALs represent a well-supported (100%) major phylogenetic branch with two groups in fungi and plants, which acquired the corresponding genes from a bacterial ancestor. The closest relatives of these sister groups are bacterial PAL from Anabaena variabilis and Nostoc punctiforme (91%). Both of these species are known to occur in diazotrophic symbiotic associations with fungi and plants that support the assumption of quite early genetic transfer of PALs from the bacterial ancestor to fungi and plants through an ancient symbiosis (Meeks et al., 2001). In contrast to PALs in bacteria and fungi, the phylogeny of the plant PALs diversified into three major sub-branches, including the acrogymnospermae, the liliopsida, and the eudicotyledons, which were described as the largest source of PALs. Among the plant PAL clusters, the acrogymnosperm clade was located between fungi and angiosperm and split into two mainly well-supported branches (100%) with monofunctional PAL from Pinus taeda in one branch as well as the PAL from genus Taxus (Taxus canadensis and Taxus wallichiana var. chinensis), which exhibits PAM activity as a side reaction, in the other branch. Furthermore, D. discoideum, as an obvious outsider among the eukaryotic HAL with the strong support of its position (100%), formed the deepest branching lineage in this large group and linked to the PALs clusters, indicating a probable single event of HGT from a prokaryotic ancestor in early age.
The second large group was further clearly separated into two branches with great statistical support (96%). One branch included all bacterial TALs/TAMs, which were found within the deep branching stable lineage, indicating that they were closer to the common ancestor than bacterial HALs. In the other branch, the genus Streptomyces (Streptomyces avermitilis, Streptomyces coelicolor, and Streptomyces griseus) developed a distinct bacterial HAL clade that divided off at the branching point from genus Thermoplasma (Thermoplasma acidophilum and Thermoplasma volcanium). Although this clade was located at a poorly supported position (44%), it still suggests that the HGT of the HAL gene could have occurred between the domains of Bacteria and Archaea. The phylogenetic tree shows the bacterial HAL mainly distributed into two large branches with all Proteobacteria in one branch and the Terrabacteria, Fusobacteria, and Spirochaetes in another branch. In the second branch, the eukaryotic HALs arose at an unsupported position, which was close to the bacterial clade of genus Deinococcus.
Sequence Alignment
Usually, PALs and TALs have strict specificity for their natural substrates. Among them, the monofunctional TALs occur mainly in the bacterial kingdom though they are relatively rare in the archaea and eukaryotes. However, partial ammonia-lyases in fungi (R. toruloides) and monocots (Zea mays) are bifunctional with similar efficiency to phenylalanine and tyrosine. It is intriguing to assume the appearance of these PTALs by genetic transformation directly from bacterial TALs. Alternatively, as the bifunctional enzymes emerge in random species, this neofunctionalization event might be explained by mutations of residues in the enzymatic active center.
To explain the appearance of this bifunctional enzyme in plants and fungi, the sequence alignment was created with all enzymes with TAL activity and PALs in monocots from the phylogenetic tree (Figure 5). The carboxylate group of arylalanine amino acids is bound in a network of hydrogen bonds with the highly conserved residues of Asn247, Asn380, and Arg341 (in ZmPTAL) to stabilize the carbanion intermediate (Calabrese et al., 2004; Wu et al., 2012). The residues that surround the aromatic ring of amino acids (His123, Leu124, Tyr440, Lys443, Ile447, and Asp471 in ZmPTAL) are usually the key contributor to substrate selectivity (Supplementary Figure S2) (Baedeker and Schulz, 2002; Ritter and Schulz, 2004). Comparison of these residues from resulting multiple sequence alignments suggests that the RtPTAL and ZmPTAL shared a higher degree of similarity in residue groupings with monofunctional PALs in monocots than bacterial TAL/TAMs. It supports the hypothesis that the bifunctional PTAL might originate in monocots by mutation and then were transferred through HGT to fungi.
Based on mutagenesis research and sequence alignment analysis, it is apparent that residues 137 and 138 strongly influence substrate selectivity in the AAM family. The His89Phe (His123 in ZmPTAL) mutation introduced the hydrophobic residue in TAL from Rhodobacter sphaeroides, resulting in complete substrate selectivity switch from tyrosine to phenylalanine (Watts et al., 2006).
Discussion
The Origin of AAMs and Mass Extinction Events
Interestingly, members of the AAM family usually are reported to play key roles in the production of radioprotective intermediates in both animals and plants. For instance, urocanic acid, the deaminated product from histidine by HALs, is a major epidermal chromophore that provides protection against UV-induced immunosuppression (Noonan et al., 1992). Similarly, in the plant kingdom, PALs are involved in the biosynthesis of flavonoids, which serve as DNA-protective metabolites against UV damage (Kootstra, 1994). Before the formation of the stratospheric ozone layer, ancient organisms suffered under higher UV exposure, which limited the possibility of land colonization (Hessen, 2008; Voosen, 2020). Thus, the ocean surface provided the primary protective shield for early marine organisms by adsorbing the most of the sun’s harmful ultraviolet radiation (Whitehead et al., 2000). The presence of the stratospheric ozone layer in the earlier Cambrian (around 600 million years ago) is thought to have influenced the divergence of multicellular animals, which led to the Cambrian explosion (Falkowski, 2011). The ozone layer served as a protective shield and promoted the colonization of land by higher plants and arthropods. The latter invaded and colonized the land during the Ordovician period, whereas the land plants appeared later (Rota-Stabelli et al., 2013). During land colonization, these organisms were exposed to a higher level of UV radiation than under water (Averof et al., 1995; Rubinstein et al., 2010), which presumably favored the evolution of the AAM gene in the common ancestors of early land dwellers. During this evolution of the related gene, the ancestors of vascular plants and terrestrial animals may have obtained or stimulated the ability to prevent damage from UV radiation.
The phylogenetic relations demonstrate that eukaryotic HAL probably evolved from prokaryotic ancestors earlier than plant PAL in multiple independent events, consistent with molecular time trees (Rota-Stabelli et al., 2013). Both arthropod and green algae are considered to be earlier colonizers of the terrestrial environment. Curiously, available genomic data indicate that the HUT pathway is missing in most lower metazoans and the arthropod, whereas the PALs are absent in red and green algae based on the available genomic data (Emiliani et al., 2009; Bender, 2012). The waxy cuticle of green algae and the exoskeleton or chitinous cuticles of arthropods might be considered to be various strategies to escape excessive UV exposure on the land. It can also be argued that AAMs did not evolve immediately after land colonization. The major extinction events were synchronous with volcanism. The volcanic gases resulted in ozone layer depletion, thereby elevating UV radiation on land (Lindström et al., 2019). Sculpture malformation in plant spores, spore tetrads, and pollen indicates the biological stress and evolutionary pressure on plants by increased UV radiation during the extinction interval (Benca et al., 2018; Marshall et al., 2020). Under high UV intensities, the frequency of this malformation increased in plant spores and pollen because of the DNA damage before the formation of the protective wall layer (Fields et al., 2020). Most immediately, ancient organisms, especially the ones that lived on land and in shallow water, were exposed to a high level of UV radiation during this period. As a strong mutagenic agent, UV irradiation might favor the appearance of AAM genes or at least the natural selection of organisms with AAM. Notably, the absence of the pal gene in terrestrial vertebrates and the disappearance of the HUT pathway in fungi and plants may also suggest a huge elimination of the old dominant species during the mass extinction event.
Microbial Coevolution in the Early Terrestrial Ecosystem
Ancient bacteria faced the major challenge of limited natural resources and space and conflict with their neighbors within a community by the production of small antibiotic compounds. As the precursor of various antibiotics, β-tyrosine was generated through bacteria with TAM to compete against others. Based on the available genomic information, the representative bacterial species with TAL/TAM activity could be isolated from soil, for example, Cupriavidus, Streptomyces, and Myxococcus, which suggests that this gene arose in an early soil environment.
Except for the RsTAL from Rhodobacter sphaeroides, the other TALs exhibit tyrosine aminomutase activity with different preferred enantioselectivity, matching the configuration of their β-tyrosine-containing secondary metabolites: (S)-TAM from Streptomyces globisporus involves the synthesis of (S)-3-chloro-5-hydroxy-β-tyrosine in the antitumor antibiotic C-1027; In Chondromyces crocatus, TAM converts (S)-α-tyrosine to (R)-β-tyrosine for production of the cytotoxic chondramides (Rachid et al., 2007). Furthermore, the arrangement of enantioselectivity within the bacterial TAM clades reflects an independent evolutionary relationship, indicating that these traits of TAM may have evolved convergently (Krug and Müller, 2009). According to the most deeply branching line formed by RsTAL in the phylogenetic tree, it could be assumed that the ancestors of bacterial TAMs may have evolved from ancient bacterial TAL for the production of the chemical inhibitor and to impair nearby competitive organisms.
The hutH gene, encoding the cytosolic enzyme HAL, is considered an ancient and basal gene, participating in a core metabolism for the degradation pathway of l-histidine. It distributes broadly in bacteria and promotes them to utilize histidine as carbon and nitrogen sources (Fuchs and Kane, 1985). The loss of the HUT pathway in bacteria only leads to their lack of sufficient capacity for the utilization of histidine, but this is not lethal. Based on the phylogeny, the most coherent hypothesis for the HAL origin would be the HAL-like enzyme gene transfer from ancient TAL-bacteria to Streptomyces bacteria and/or thermophilic or halophilic archaea at quite early ages. Among the other archaea as well as among lower eukaryotes, the appearance of HAL is irregular and rare, whereas all vertebrates, especially mammals, require the hutH gene to metabolize histidine and suffer the disease caused by a deficiency of HAL (Ghadimi et al., 1962; Taylor et al., 1991). Of the protozoa HAL, it has so far only been assured in D. discoideum, which is identified as two homologs: the one placed in the eukaryotic HAL clade and the other one situated close to the cyanobacteria clade. It is noteworthy that the HUT pathway is absent from cyanobacteria. In cyanobacteria, the enzyme to catalyze phenylalanine conversion was the other member of AAM: PAL.
The bacterial PALs were rooted at the deepest position among the PAL clades, which is a clear indication of the original PAL emergence in bacteria. According to the evolutionary analysis of Glomeromycotina, it is proposed that the ancestral fungi have been found in symbiotic association with green algae or cyanobacteria before terrestrial colonization (Prasad et al., 2020). Nevertheless, neither HAL nor PAL orthologs in red and green algae were identified in the available genomic data. Besides this, the existing genomic information of bacteria and fungi with the pal gene, even the nearby amoebae with HAL activity, have a relatively high dependency on the soil environment, which might indicate that the gene transfer arose during coevolution between bacterial, fungi, and plants in the early terrestrial environment (Emiliani et al., 2009).
In most plant species, PAL is encoded by a multigene family, containing up to five members that express differently in numerous tissues or reply to various environmental stress conditions. For example, in Arabidopsis thaliana, pal3 is expressed at a relatively lower level than the other three genes (pal1, pal2, and pal4) in stem tissue, whereas pal3 seems to be mainly expressed in leaves, the most UV exposed plant tissue. Among these three genes, pal1 and pal2, sharing the common promoter elements, mainly exist in roots and stems in biosynthesis processes related to abiotic imports, whereas the gene pal4 was expressed also in seeds (Raes et al., 2003). This gene duplication could be traced back to gymnosperms: P. taeda owns five pal genes, representing different tissues with various levels of expression (Bagal et al., 2012). It seems that plants received the pal genes from the bacterial ancestor and expanded mainly through gene duplication. Moreover, Shang et al. report that seven pal genes are tandemly arranged with four pseudogenes in two chromosomes from Cucumis sativus and point out that these pal genes were duplicated more recently as the splitting of cucumber from the other dicots (Shang et al., 2012).
Structurally, the quaternary sizes of the prokaryotic HALs (∼500 amino acids) and eukaryotic PALs (∼710 amino acids) from the AAM superfamily are diverse because of the presence of an additional C-terminal multihelix domain in the latter enzyme (Figure 6) (Allwood et al., 1999; Ritter and Schulz, 2004; Pilbák et al., 2006). These extended domains in plants are implicated to destabilize the enzyme for rapid regulation of the phenylpropanoid biosynthesis and adaptation of the varying environmental stresses. The phosphorylation site of PAL from Phaseolus vulgaris is determined as Thr545, whereas the most accessible cleavage sites by trypsin and chymotrypsin from R. toruloides PAL are identified as Arg123 and Tyr110, which are located in the extended regions to decrease the lifetime of PAL (Gámez et al., 2005).
FIGURE 6. (A) Overlap X-ray crystal structures of PpHAL from Pseudomonas putida (PDB entry 1GKM, blue), RtPAL from R. toruloides (PDB entry 1T6J, yellow), and RsTAL from Rhodobacter sphaeroides (PDB entry 2O6Y, pink) (B) The inner (blue) and outer (red) loops in RsTAL from two chains. The MIO is colored in red.
PTAL and PAM
PTAL, a specifical ammonia-lyase in monocots (especially the Poaceae family), is associated with the biosynthesis of stress-induced syringyl-rich lignins (Barros and Dixon, 2020). Among eight pal genes from the Brachypodium distachyon, only BdPTAL1 was expressed as the enzyme with additional TAL activity (Barros et al., 2016). The monocots with PTAL activity clustered with each other and were distinct from dicots. It can be assumed that this gene acquired an inverted substrate selectivity toward tyrosine after the early duplication event in monocots. The pattern of distribution of fungal and plant PTALs in the phylogenetic tree (Figure 4) and sequence alignment results (Supplementary Figure S2) might suggest that this gene originated in monocots and then transferred through HGT to fungi.
In other instances, yew species produce defensive metabolites, paclitaxel derived from β-phenylalanine, to protect themselves against widespread wood-degrading fungi (Malik et al., 2011). As long-lived species, yew trees are reasonably susceptible to fungal infection because they can form lateral buds on the old branches and the stem, which lead to a crash of barks (Thomas and Polwart, 2003). In these bark-cracking tissues, paclitaxel is mainly accumulated to hinder the pathogen attack (Talbot, 2015). The rate-limiting step of the side chain assembly process in the taxol biosynthesis is catalyzed by PAM that presumably evolved from the plant PAL ancestors in gymnosperms (Walker et al., 2004; Steele et al., 2005). Compared with the yew species, the other gymnosperms, such as P. taeda, employ the monofunctional PAL without any aminomutase activity. According to phylogenetic analysis, it is tough to conclude whether the PAL from P. taeda lost the catalysis ability as PAM during evolution in the early ages or the PAM gene in Taxus species comes from the ancient bacteria that obtained the PAM-like aminomutase. The taxol-producing endophytic fungi, including the Penicillium species, which tend to have PAL with aminomutase activity, could be isolated from the genus Taxus (Soliman and Raizada, 2013). Moreover, epiphytic and pathogenic fungi from the other hosts as well as saprophytic fungi are reported as taxol producers, especially the pathogenic fungi Pestalotiopsis malicola, which is isolated from soil, proposing another possible evolutionary scenario that points to a probable PAM origin from fungi (Bi et al., 2011). Based on this hypothesis, the old gymnosperms might obtain genes through HGT against fungal infection.
The dynamic simulation indicates that the inner loop is more conformationally flexible at higher temperatures or with more hydrophilic residues (A77T, I79S, C89T, and L97G in Taxus chinesis PAM), affording lyase features to aminomutase. Furthermore, the reaction temperature also influences the distance between phenolic O-atom of catalytic essential residue tyrosine and exocyclic methylene C-atom of the MIO prosthetic group in PcPAL from Pelargonium crispum, whereas HALs are highly tolerant to temperatures due to their stable inner loop (Pilbák et al., 2006; Heberling et al., 2015; Attanayake et al., 2018). Thus, the inner loop flexibility and environmental temperature are the determinants that distinguish the mutase vs. lyase activities. It is noteworthy that the growth rates of various fungal pathogens, like taxol producer P. islandicum, are significantly enhanced by the combination of warm temperatures and high humidity (Mannaa and Kim, 2018). It may indicate that the PAM emergence in plants is related to their protecting themselves against the pathogen fungi.
Conclusion
The AAMs share a similar structure with highly electrophilic MIO and a common catalytic mechanism. Among them, eukaryotic PALs and HALs are involved in various metabolic and catabolic pathways to form several protection compounds under UV radiation. The corresponding genes most probably originated because of the fluctuating UV intensity. Our results indicate that the bacterial TALs have been developed and can produce antibacterial compounds under water-limited environments in the soils. Under UV exposure at an early age, some ancient bacteria were able to metabolize histidine by HAL or accepted phenylalanine as substrates by PAL, respectively. During further evolution, land plants and fungi obtained the PAL gene from bacteria through an early symbiosis while the terrestrial vertebrates inherited the HAL gene from their bacterial ancestor.
To conflict with the other organisms and withstand the pathogen infraction, the PAL gene was expressed in certain plants or fungi species as bifunctional PTAL or PAM. These corresponding genes were transferred through HGT between these two phyla. It may indicate that a symbiotic association involving bacteria, fungi, plants, and amoebae occurred in ancient terrestrial environments. Unlike PAM and TAM, there is no obvious genome evidence of the existence of aminomutase to produce β-histidine. It cannot be excluded that histidine aminomutase did not provide great advantages and, thus, disappeared during later evolution.
This review provides a phylogenetic framework for further evolutionary research of AAM, which has multiple meanings for studies of engineering in metabolic pathways and enzymatic biotechnology. Furthermore, more genomic data of new species and phylogenetic lines of this enzyme superfamily as well as substrate specificity data from well-investigated enzymes are also needed to find the missing puzzle pieces and integrate the evolutionary history.
Methods
Phylogenetic Analysis
To explore the phylogenetic history of arylalanine amino acid AAM, 268 available sequences were identified from SwissProt (Bairoch et al., 2004). These sequences were aligned using MAFFT version 7.3 (Katoh and Standley, 2013). After deleting the gaps within the resulting alignments by Gblocks (Talavera and Castresana, 2007), the 100% identified sequences were manually trimmed. The maximum likelihood trees were performed using the IQ-TREE version 1.4 (Minh et al., 2020) with the best-fit substitution model (LG + G) inferred by PhyML SMS servers (Lefort et al., 2017). Branch supports were assessed with 1000 bootstrap replicates and rooted using a minimal ancestor deviation (MAD) approach (Tria et al., 2017). The resulting phylogenetic tree was visualized with annotations by the online tool “interactive Tree of Life” (iTol) (Letunic and Bork, 2019).
Sequence Analysis
The sequences of all TALs and certain PALs from monocots were identified from the SwissPort in UniProt knowledgebase (UniProtKB) and aligned with MAFFT version 7.3. Abbreviations and accession numbers from Swissport are RsTAL, Rhodobacter sphaeroides TAL (Q3IWB0); SgTAM, Streptomyces globisporus TAM (Q8GMG0); CmTAL, Cupriavidus metallidurans TAL (Q1LRV9); CcTAM, Chondromyces crocatus TAM (Q0VZ68); MfTAM, Myxococcus fulvus TAM (B8ZV93); RtPTAL, R. toruloides PTAL (P11544); ZmPTAL, Zea mays PTAL (Q8VXG7); BfPAL, Bromheadia finlaysoniana PAL (Q42609); NpPAL, Narcissus pseudonarcissus PAL (A0A2H5AIY6); OsPAL, Oryza sativa subsp. Japonica (P14717); TaPAL, Triticum aestivum PAL (Q43210).
Data Availability Statement
The data sets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author Contributions
Methodology, HA and FP; phylogenetic analysis, HA and FP; data analysis, FP; writing the manuscript, FP; revision of the manuscript and supervision, JR, UE and HA. All authors have read and agreed to the published version of the manuscript.
Funding
This work is supported by China Scholarship Council and is part of the bioeconomy graduate program BBW ForWerts of the German state of Baden-Württemberg. We gratefully thank the Open Access Publishing Fund of Karlsruhe Institute of Technology.
Conflict of Interest
The 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/fgene.2022.851738/full#supplementary-material
Supplementary Figure S1 | Bootstrap tree of arylalanine amino acid ammonia-lyase and aminomutase homologs retrieved from Swissprot.
Supplementary Figure S2 | Partial sequence alignment of eukaryotic PTALs, prokaryotic TAL/TAMs, and certain PAL in monocots. Shown is the location of the residues in the active center that are suggested to represent major roles for substrate selectivity. The specific residues in various domains are highlighted in different colors.
Abbreviations
4CL, 4-coumaric acid-CoA ligase; AAM, Aromatic amino acid MIO-dependent enzyme; C3H, Cytosolic coumarate 3-hydroxylase; C4H, Cinnamic acid 4-hydroxylase; COMT, 3-O-methyltransferase; GFP, Green fluorescent protein; HAL, Histidine ammonia-lyase; HGT, Horizontal gene transfer; HUT, Histidine pathway; MIO, 3,5-dihydro-5-methylidene-4H-imidazole-4-one; PAL, Phenylalanine ammonia-lyase; PAM, Phenylalanine aminomutase; PTAL, Phenylalanine/tyrosine ammonia-lyase; PTM, Posttranslational modification; RFP, Red fluorescent protein; TAL, Tyrosine ammonia-lyase; TAM, Tyrosine aminomutase; UV, ultraviolet.
References
Allwood, E. G., Davies, D. R., Gerrish, C., Ellis, B. E., and Bolwell, G. P. (1999). Phosphorylation of Phenylalanine Ammonia-Lyase: Evidence for a Novel Protein Kinase and Identification of the Phosphorylated Residue. FEBS Lett. 457, 47–52. doi:10.1016/S0014-5793(99)00998-9
Attanayake, G., Walter, T., and Walker, K. D. (2018). Understanding Which Residues of the Active Site and Loop Structure of a Tyrosine Aminomutase Define its Mutase and Lyase Activities. Biochemistry 57, 3503–3514. doi:10.1021/acs.biochem.8b00269
Averof, M., Crc, W., and Akam, M. (1995). Insect-crustacean Relationships: Insights from Comparative Developmental and Molecular Studies. Phil. Trans. R. Soc. Lond. B 347, 293–303. doi:10.1098/rstb.1995.0028
Baedeker, M., and Schulz, G. E. (2002). Autocatalytic Peptide Cyclization during Chain Folding of Histidine Ammonia-Lyase. Structure 10, 61–67. doi:10.1016/S0969-2126(01)00692-X
Bagal, U. R., Leebens-Mack, J. H., Lorenz, W. W., and Dean, J. F. (2012). The Phenylalanine Ammonia Lyase (PAL) Gene Family Shows a Gymnosperm-specific Lineage. BMC Genomics 13, 1–9. doi:10.1186/1471-2164-13-s3-s1
Bairoch, A., Boeckmann, B., Ferro, S., and Gasteiger, E. (2004). Swiss-prot: Juggling between Evolution and Stability. Brief. Bioinform. 5, 39–55. doi:10.1093/bib/5.1.39
Barros, J., and Dixon, R. A. (2020). Plant Phenylalanine/Tyrosine Ammonia-Lyases. Trends Plant Sci. 25, 66–79. doi:10.1016/j.tplants.2019.09.011
Barros, J., Escamilla-Trevino, L., Song, L., Rao, X., Serrani-Yarce, J. C., Palacios, M. D., et al. (2019). 4-Coumarate 3-hydroxylase in the Lignin Biosynthesis Pathway Is a Cytosolic Ascorbate Peroxidase. Nat. Commun. 10, 1–11. doi:10.1038/s41467-019-10082-7
Barros, J., Serrani-Yarce, J. C., Chen, F., Baxter, D., Venables, B. J., and Dixon, R. A. (2016). Role of Bifunctional Ammonia-Lyase in Grass Cell wall Biosynthesis. Nat. Plants 2, 1–9. doi:10.1038/NPLANTS.2016.50
Benca, J. P., Duijnstee, I. A. P., and Looy, C. V. (2018). UV-B-induced forest Sterility: Implications of Ozone Shield Failure in Earth's Largest Extinction. Sci. Adv. 4 (2), e1700618. doi:10.1126/sciadv.1700618
Bender, R. A. (2012). Regulation of the Histidine Utilization (Hut) System in Bacteria. Microbiol. Mol. Biol. Rev. 76, 565–584. doi:10.1128/mmbr.00014-12
Bi, J., Ji, Y., Pan, J., Yu, Y., Chen, H., and Zhu, X. (2011). A New Taxol-Producing Fungus (Pestalotiopsis Malicola) and Evidence for Taxol as a Transient Product in the Culture. Afr. J. Biotechnol. 10, 6647–6654. doi:10.5897/AJB11.170
Bischoff, R., and Schlüter, H. (2012). Amino Acids: Chemistry, Functionality and Selected Non-enzymatic post-translational Modifications. J. Proteomics 75, 2275–2296. doi:10.1016/j.jprot.2012.01.041
Brosnan, M. E., and Brosnan, J. T. (2020). Histidine Metabolism and Function. J. Nutr. 150, 2570S–2575S. doi:10.1093/jn/nxaa079
Calabrese, J. C., Jordan, D. B., Boodhoo, A., Sariaslani, S., and Vannelli, T. (2004). Crystal Structure of Phenylalanine Ammonia Lyase: Multiple Helix Dipoles Implicated in Catalysis,. Biochemistry 43, 11403–11416. doi:10.1021/bi04905310.1021/bi049053+
Emiliani, G., Fondi, M., Fani, R., and Gribaldo, S. (2009). A Horizontal Gene Transfer at the Origin of Phenylpropanoid Metabolism: A Key Adaptation of Plants to Land. Biol. Direct 4, 1–12. doi:10.1186/1745-6150-4-7
Falkowski, P. G. (2011). The Biological and Geological Contingencies for the Rise of Oxygen on Earth. Photosynth. Res. 107, 7–10. doi:10.1007/s11120-010-9602-4
Fields, B. D., Melott, A. L., Ellis, J., Ertel, A. F., Fry, B. J., Lieberman, B. S., et al. (2020). Supernova Triggers for End-Devonian Extinctions. Proc. Natl. Acad. Sci. USA 117 (35), 21008–21010. doi:10.1073/pnas.2013774117
Fuchs, R. L., and Kane, J. F. (1985). In Vivo synthesis of Histidine by a Cloned Histidine Ammonia-Lyase in Escherichia coli. J. Bacteriol. 162, 98–101. doi:10.1128/jb.162.1.98-101.1985
Gámez, A., Sarkissian, C. N., Wang, L., Kim, W., Straub, M., Patch, M. G., et al. (2005). Development of Pegylated Forms of Recombinant Rhodosporidium Toruloides Phenylalanine Ammonia-Lyase for the Treatment of Classical Phenylketonuria. Mol. Ther. 11, 986–989. doi:10.1016/j.ymthe.2005.02.013
Ghadimi, H., Partington, M. W., and Hunter, A. (1962). Inborn Error of Histidine Metabolism. Pediatrics 29, 714–728. doi:10.1542/peds.29.5.714
Heberling, M. M., Masman, M. F., Bartsch, S., Wybenga, G. G., Dijkstra, B. W., Marrink, S. J., et al. (2015). Ironing Out Their Differences: Dissecting the Structural Determinants of a Phenylalanine Aminomutase and Ammonia Lyase. ACS Chem. Biol. 10, 989–997. doi:10.1021/cb500794h
Katoh, K., and Standley, D. M. (2013). MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability. Mol. Biol. Evol. 30, 772–780. doi:10.1093/molbev/mst010
Kohlmeier, M. (2015). “Chapter 8-Amino Acids and Nitrogen Compounds,” in Nutrient Metabolism. Editor M. Kohlmeier. Second Edition (San Diego: Academic Press), 265–477. doi:10.1016/B978-0-12-387784-0.00008-0
Kootstra, A. (1994). Protection from UV-B-Induced DNA Damage by Flavonoids. Plant Mol. Biol. 26, 771–774. doi:10.1007/BF00013762
Krug, D., and Müller, R. (2009). Discovery of Additional Members of the Tyrosine Aminomutase Enzyme Family and the Mutational Analysis of CmdF. ChemBioChem 10, 741–750. doi:10.1002/cbic.200800748
Lefort, V., Longueville, J. E., and Gascuel, O. (2017). SMS: Smart Model Selection in PhyML. Mol. Biol. Evol. 34, 2422–2424. doi:10.1093/molbev/msx149
Letunic, I., and Bork, P. (2019). Interactive Tree of Life (iTOL) V4: Recent Updates and New Developments. Nucleic Acids Res. 47, 256–259. doi:10.1093/nar/gkz239
Lindström, S., Sanei, H., Van De Schootbrugge, B., Pedersen, G. K., Lesher, C. E., Tegner, C., et al. (2019). Volcanic Mercury and Mutagenesis in Land Plants during the End-Triassic Mass Extinction. Sci. Adv. 5 (10), eaaw4018. doi:10.1126/sciadv.aaw4018
Malik, S., Cusidó, R. M., Mirjalili, M. H., Moyano, E., Palazón, J., and Bonfill, M. (2011). Production of the Anticancer Drug Taxol in Taxus Baccata Suspension Cultures: A Review. Process. Biochem. 46, 23–34. doi:10.1016/j.procbio.2010.09.004
Mannaa, M., and Kim, K. D. (2018). Effect of Temperature and Relative Humidity on Growth of Aspergillus and Penicillium Spp. And Biocontrol Activity of Pseudomonas Protegens AS15 against Aflatoxigenic Aspergillus Flavus in Stored rice Grains. Mycobiology 46, 287–295. doi:10.1080/12298093.2018.1505247
Marshall, J. E., Lakin, J., Troth, I., and Wallace-Johnson, S. M. (2020). UV-B Radiation Was the Devonian-Carboniferous Boundary Terrestrial Extinction Kill Mechanism. Sci. Adv. 6 (22), eaba0768. doi:10.1126/sciadv.aba0768
Meeks, J. C., Elhai, J., Thiel, T., Potts, M., Larimer, F., Lamerdin, J., et al. (2001). An Overview of the Genome of Nostoc Punctiforme, a Multicellular, Symbiotic Cyanobacterium. Photosynth. Res. 70, 85–106. doi:10.1023/a:1013840025518
Minh, B. Q., Schmidt, H. A., Chernomor, O., Schrempf, D., Woodhams, M. D., Von Haeseler, A., et al. (2020). IQ-TREE 2: New Models and Efficient Methods for Phylogenetic Inference in the Genomic Era. Mol. Biol. Evol. 37, 1530–1534. doi:10.1093/molbev/msaa015
Müller, M. M. (2018). Post-translational Modifications of Protein Backbones: Unique Functions, Mechanisms, and Challenges. Biochemistry 57, 177–185. doi:10.1021/acs.biochem.7b00861
Nagy, E. Z. A., Tork, S. D., Lang, P. A., Filip, A., Irimie, F. D., Poppe, L., et al. (2019). Mapping the Hydrophobic Substrate Binding Site of Phenylalanine Ammonia-Lyase from Petroselinum crispum. ACS Catal. 9, 8825–8834. doi:10.1021/acscatal.9b02108
Noonan, F. P., and De Fabo, E. C. (1992). Immunosuppression by Ultraviolet B Radiation: Initiation by Urocanic Acid. Immunol. Today 13, 250–254. doi:10.1016/0167-5699(92)90005-R
Norval, M. (2001). “Chapter 5-Effects of Solar Radiation on the Human Immune System,” in Sun Protection in Man Comprehensive Series in Photosciences. Editor P. U. Giacomoni (Elsevier), 91–113. doi:10.1016/S1568-461X(01)80040-5
Pilbák, S., Tomin, A., Rétey, J., and Poppe, L. (2006). The Essential Tyrosine-Containing Loop Conformation and the Role of the C-Terminal Multi-helix Region in Eukaryotic Phenylalanine Ammonia-Lyases. FEBS J. 273, 1004–1019. doi:10.1111/j.1742-4658.2006.05127.x
Prasad, R., Chhabra, S., Gill, S. S., Singh, P. K., and Tuteja, N. (2020). The Microbial Symbionts: Potential for Crop Improvement in Changing Environments. Adv. Crop Improv. Tech., 233–240. doi:10.1016/b978-0-12-818581-0.00014-0
Rachid, S., Krug, D., Weissman, K. J., and Müller, R. (2007). Biosynthesis of (R)-β-tyrosine and its Incorporation into the Highly Cytotoxic Chondramides Produced by Chondromyces Crocatus. J. Biol. Chem. 282, 21810–21817. doi:10.1074/jbc.M703439200
Raes, J., Rohde, A., Christensen, J. H., Van De Peer, Y., and Boerjan, W. (2003). Genome-Wide Characterization of the Lignification Toolbox in Arabidopsis. Plant Physiol. 133, 1051–1071. doi:10.1104/pp.103.026484
Reid, B. G., and Flynn, G. C. (1997). Chromophore Formation in green Fluorescent Protein. Biochemistry 36, 6786–6791. doi:10.1021/bi970281w
Ritter, H., and Schulz, G. E. (2004). Structural Basis for the Entrance into the Phenylpropanoid Metabolism Catalyzed by Phenylalanine Ammonia-Lyase. Plant Cell 16, 3426–3436. doi:10.1105/tpc.104.025288
Rota-Stabelli, O., Daley, A. C., and Pisani, D. (2013). Molecular Timetrees Reveal a Cambrian Colonization of Land and a New Scenario for Ecdysozoan Evolution. Curr. Biol. 23, 392–398. doi:10.1016/j.cub.2013.01.026
Röther, D., Poppe, L., Morlock, G., Viergutz, S., and Rétey, J. (2002). An Active Site Homology Model of Phenylalanine Ammonia-Lyase from Petroselinum crispum. Eur. J. Biochem. 269, 3065–3075. doi:10.1046/j.1432-1033.2002.02984.x
Rubinstein, C. V., Gerrienne, P., de la Puente, G. S., Astini, R. A., and Steemans, P. (2010). Early Middle Ordovician Evidence for Land Plants in Argentina (Eastern Gondwana). New Phytol. 188, 365–369. doi:10.1111/j.1469-8137.2010.03433.x
Sánchez-Murcia, P. A., Bueren-Calabuig, J. A., Camacho-Artacho, M., Cortés-Cabrera, Á., and Gago, F. (2016). Stepwise Simulation of 3,5-Dihydro-5-Methylidene-4h-Imidazol-4-One (MIO) Biogenesis in Histidine Ammonia-Lyase. Biochemistry 55, 5854–5864. doi:10.1021/acs.biochem.6b00744
Shang, Q. M., Li, L., and Dong, C. J. (2012). Multiple Tandem Duplication of the Phenylalanine Ammonia-Lyase Genes in Cucumis Sativus L. Planta 236, 1093–1105. doi:10.1007/s00425-012-1659-1
Sharma, A., Shahzad, B., Rehman, A., Bhardwaj, R., Landi, M., and Zheng, B. (2019). Response of Phenylpropanoid Pathway and the Role of Polyphenols in Plants under Abiotic Stress. Molecules 24, 1–22. doi:10.3390/molecules24132452
Soliman, S. S. M., and Raizada, M. N. (2013). Interactions between Co-habitating Fungi Elicit Synthesis of Taxol from an Endophytic Fungus in Host Taxus Plants. Front. Microbiol. 4, 1–14. doi:10.3389/fmicb.2013.00003
Steele, C. L., Chen, Y., Dougherty, B. A., Li, W., Hofstead, S., Lam, K. S., et al. (2005). Purification, Cloning, and Functional Expression of Phenylalanine Aminomutase: The First Committed Step in Taxol Side-Chain Biosynthesis. Arch. Biochem. Biophys. 438, 1–10. doi:10.1016/j.abb.2005.04.012
Talavera, G., and Castresana, J. (2007). Improvement of Phylogenies after Removing Divergent and Ambiguously Aligned Blocks from Protein Sequence Alignments. Syst. Biol. 56, 564–577. doi:10.1080/10635150701472164
Talbot, N. J. (2015). Plant Immunity: A Little Help from Fungal Friends. Curr. Biol. 25, R1074–R1076. doi:10.1016/j.cub.2015.09.068
Taylor, R. G., Levy, H. L., and McInnes, R. R. (1991). Histidase and Histidinemia. Clinical and Molecular Considerations. Mol. Biol. Med. 8, 101–116. Available at: http://europepmc.org/abstract/MED/1943682.
Thomas, P. A., and Polwart, A. (2003). Taxus Baccata L. J. Ecol. 91, 489–524. doi:10.1046/j.1365-2745.2003.00783.x
Treutter, D. (2006). Significance of Flavonoids in Plant Resistance: A Review. Environ. Chem. Lett. 4, 147–157. doi:10.1007/s10311-006-0068-8
Tria, F. D. K., Landan, G., and Dagan, T. (2017). Phylogenetic Rooting Using Minimal Ancestor Deviation. Nat. Ecol. Evol. 1, 193. doi:10.1038/s41559-017-0193
Voosen, P. (2020). UV Radiation Blamed in Ancient Mass Extinction. Science 368, 926. doi:10.1126/science.368.6494.926
Walker, K. D., Klettke, K., Akiyama, T., and Croteau, R. (2004). Cloning, Heterologous Expression, and Characterization of a Phenylalanine Aminomutase Involved in Taxol Biosynthesis. J. Biol. Chem. 279, 53947–53954. doi:10.1074/jbc.M411215200
Watts, K. T., Mijts, B. N., Lee, P. C., Manning, A. J., and Schmidt-Dannert, C. (2006). Discovery of a Substrate Selectivity Switch in Tyrosine Ammonia-Lyase, a Member of the Aromatic Amino Acid Lyase Family. Chem. Biol. 13, 1317–1326. doi:10.1016/j.chembiol.2006.10.008
Whitehead, R. F., de Mora, S. J., and Demers, S. (2000). Enhanced UV Radiation–A New Problem for the marine Environment. Eff. UV Radiat. Mar. Environ. 10, 1–34. doi:10.1017/cbo9780511535444.002
Wu, B., Szymański, W., Wybenga, G. G., Heberling, M. M., Bartsch, S., de Wildeman, S., et al. (2012). Mechanism-Inspired Engineering of Phenylalanine Aminomutase for Enhanced β-Regioselective Asymmetric Amination of Cinnamates. Angew. Chemie 124, 497–501. doi:10.1002/ange.201106372
Keywords: MIO-dependent enzymes, a mass extinction event, microbial coevolution, rooted phylogenetic tree, minimal ancestor deviation
Citation: Peng F, Engel U, Aliyu H and Rudat J (2022) Origin and Evolution of Enzymes with MIO Prosthetic Group: Microbial Coevolution After the Mass Extinction Event. Front. Genet. 13:851738. doi: 10.3389/fgene.2022.851738
Received: 10 January 2022; Accepted: 16 February 2022;
Published: 29 March 2022.
Edited by:
Sunil Kumar Raghav, Institute of Life Sciences (ILS), IndiaReviewed by:
Kirk Matthew Schnorr, Novozymes, DenmarkRodrigo Ligabue-Braun, Federal University of Health Sciences of Porto Alegre, Brazil
Copyright © 2022 Peng, Engel, Aliyu and Rudat. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jens Rudat, amVucy5ydWRhdEBraXQuZWR1