Abstract
Copper-containing amine oxidases (CuAOs) catalyze the terminal oxidation of polyamines (PAs), producing ammonium, an aminoaldehyde, and hydrogen peroxide (H2O2). Plant CuAOs are induced by stress-related hormones such as methyl-jasmonate (MeJA), abscisic acid (ABA), and salicylic acid (SA). Mammalian copper-containing amine oxidases (CAOs), encoded by four genes (AOC1-4) that catalyze the oxidation of primary amines to aldehydes, regulate various biological processes and are linked to diseases like inflammatory conditions and histamine intolerance. To understand the evolutionary history and functional divergence of CuAOs, we conducted phylogenetic and expression analyses of CuAOs in plants and animals. In this study, the copper amine oxidase (CuAO) genes were identified by HMMER and BLASTP, and verified by CDD/HMM/SMART. Multiple sequence alignment was performed using Muscle5, and the phylogenetic tree was constructed by IQ-TREE2. The syntenic relationship was analyzed by MCScanX and CIRCOS. Meanwhile, the expression data of Arabidopsis thaliana and human and other species were integrated for analysis. Here, 950 and 264 CuAO orthologues were identified in 188 plant and 79 animal genomes. Phylogenetic analyses indicate that CuAO originated in the common ancestor before the divergence of plants and animals. The copy numbers of CuAOs vary significantly across plant species, whereas they remain relatively stable in animal species, generally maintaining 3-4 copies per species. During the evolutionary process, plant CuAOs formed two clades (I and II), while animal CuAOs formed three clades (CAO-like, AOC1, AOC2-4). Interestingly, plant clade I CuAOs lacks the active site motif T/S-X1-X2-N-Y-D. The further differentiation of plant clade II CuAOs is related to the preference for X1 and X2 active sites. CAO-like and AOC1 are monophyletic branches. Mammalian AOC2-4 is separated from non-mammalian AOC2-4, and the differentiation of mammalian AOC3 and AOC4 occurs in a species-specific manner. Our study provides a comprehensive understanding of the evolutionary trajectory of the CuAO gene family in plants and animals at the genome-wide level. These findings lay a crucial foundation for future research to conduct in-depth functional characterization.
Introduction
PAs are a class of small aliphatic nitrogen-containing compounds that are widely present in animals, plants, bacteria, and fungi (). As important regulators of cell growth and differentiation, PAs are involved in basic cellular processes such as DNA replication and transcription, RNA modification, protein synthesis, ion channel regulation, free radical scavenging, cell cycle regulation, signal transduction, and programmed cell death (; ). In plants, PAs play crucial roles in flower development and differentiation, leaf development and aging, fruit ripening, and various stress responses (). In mammalian cells, dysregulation of PA metabolism is associated with cancer development, making it an important target for anti-cancer therapy (; ; ). The three most common PAs found in plants and animals are putrescine (Put), spermidine (Spd), and spermine (Spm) (). Their oxidative degradation is mainly catalyzed by copper amine oxidase (CuAO) and polyamine oxidase (PAO). CuAOs mainly catalyze the oxidation of putrescine to produce NH4+, 4-aminobutyrate, and H2O2, with 4-aminobutyrate further converting to γ-aminobutyric acid (). CuAOs not only play roles in plants, but also are widely distributed in mammalian tissues, and their activity changes are related to various diseases, such as histamine intolerance and rheumatism (; ; ).
CuAOs are widely occurred enzymes that catalyze the oxidative deamination of primary amines and PAs such as diamines and histamine (). Their molecular form is a dimer, with subunit molecular weights ranging from 70 to 95 kDa (; ). The enzymatic activity of CuAOs depends on two cofactors: covalently bound 2,4,5-trihydroxyphenylalanine quinone (TPQ) and a copper ion involved in the catalytic cycle (). TPQ is produced by the modification of endogenous tyrosine residues, and copper ions play a key role in the biogenesis and catalytic reactions of TPQ-mediated substrate oxidation (). Mammalian CuAOs have a highly conserved T/S-X1-X2-N-Y-D motif at the active site, where tyrosine is post-translationally modified to TPQ (). CuAOs primarily catalyze the terminal oxidation of diamines (such as Put) and, in some cases, the oxidation of higher PAs (such as Spd), thus playing a crucial role in the early steps of PA oxidative metabolism and regulating specific subcellular compartments of PA levels (). CuAOs catalytic reactions involve three main processes: (i) TPQ biosynthesis, (ii) amine substrate oxidation to reduce TPQ, and (iii) molecular oxygen reduction via reduced TPQ (). All of these reactions require the participation of molecular oxygen ().
CuAOs have various biological functions in plants and animals. In Arabidopsis thaliana, there are 10 CuAO genes, of which AtCuAOε1 (AT4G12270) and AtCuAOε2 (AT4G12280) lack essential active residues and are considered to be fragment copies of AtCuAOδ (AT4G12290) (; ; ). AtCuAOα3 (AT1G31710) and AtCuAOζ (AT2G42490) play roles in regulating the stability of PA in peroxisomes (). AtCuAOγ1 (AT1G62810) primarily involves NO production mediated by PAs and/or abscisic acid (ABA), while AtCuAOα2 (AT1G31690) regulates arginase activity to affect the availability of arginine and thus participates in NO production, revealing a new regulatory pathway for NO production in plants (; ). AtCuAOδ plays a significant role in ABA-induced stomatal closure (), while AtCuAOβ (AT4G14940) is involved in the stomatal regulation process under mechanical damage conditions (). Furthermore, AtCuAOβ is involved in both MeJA-induced early root protoxylem differentiation and the response to wounding (; ). The expression of AtCuAOζ in guard cells indicates that it is involved in ABA-mediated stomatal opening control (). AtCuAOβ, AtCuAOγ1, and AtCuAOγ2 are located in the apoplast (; ). Additionally, AtCuAOα2 and AtCuAOα3 are located in peroxisomes (). To date, the only vacuolar located CuAO isoform is AtCuAOδ ().There are 8 CuAO genes in Solanum lycopersicum genome, among which 6 CuAOs (SlCuAO1, 2, 3, 4, 6, 7) are specifically expressed in the roots, SlCuAO5 is specifically expressed in the fruits, and SlCuAO8 is specifically expressed in the flowers (). This indicates that the expression of amine oxidases is diverse in different tissues, and it is particularly prominent in the roots (). In humans, there are four genes that encode CuAO: AOC1 (diamine oxidase), AOC2 (retinal-specific amine oxidase), AOC3 (vascular adhesion protein-1, VAP-1), and AOC4 (pseudogene) (). AOC1 is mainly expressed in the kidney, placenta, intestine, thymus, and seminal vesicle, and is the main enzyme for metabolizing histamine (; ). Its activity is related to the risk of pregnancy (). AOC2 was initially cloned from the retina and is also expressed in adipose tissue, playing a role in the process of adipocyte differentiation (; ; ). The substrates of AOC2 include 2-phenylethylamine, tryptamine and p-tyramine, but it does not oxidize histamine (). AOC3 is highly expressed in adipocytes, smooth muscle cells, and endothelial cells, and is abundantly present in the lungs, aorta, liver, and ileum ().
While the roles of CuAOs are slowly being clarified, studies on these enzymes remain confined to a limited number of model organisms, including humans and A. thaliana. Currently, our understanding of the origin and evolution of CuAOs is sparse, and a comprehensive investigation into CuAOs across both plant and animal kingdoms is still missing. In this study, we explored the evolutionary path of CuAOs by examining the complete genomic sequences from various plants and animals. As far as we are aware, this represents the most thorough phylogenetic analysis of CuAOs conducted so far. This research offers an in-depth view of the evolution of CuAOs, examining their origins, evolutionary processes, and functional diversity, while also establishing a strong basis for future studies on functional resolution and molecular evolution.
Results
Identification and distribution of CuAO genes
In order to comprehensively identify CuAO in plants and animals, we selected a diverse range of 188 plant species and 79 animal species, including rhodophytes, chlorophytes, charophytes, bryophytes (hornworts, mosses, liverworts), ferns, lycophytes, gymnosperms, basal angiosperms, angiosperms, invertebrates and vertebrates (Table 1, Supplementary Table S1). Through homology search and domain prediction, a total of 1214 candidate protein sequences of CuAO were retrieved, including 950 plant sequences and 264 animal sequences (Table 1, Supplementary Table S2). CuAO copy numbers are conserved among different animal lineages, while they vary significantly among different plant lineages (Table 1). The average count of CuAOs was 2 in rhodophytes, 1.6 in chlorophytes, 1.8 in charophytes, 1.29 in hornworts, 3.38 in mosses, 2.59 in liverworts, 3 in ferns, 3.14 in lycophytes, 4.07 in gymnosperms, 2 in basal angiosperms, 2.36 in magnoliids, 5.57 in monocots, and 7.53 in eudicots. The copy numbers of CuAO were lower in algal species, higher in land plants, and significantly elevated in seed plants. This change may be closely related to the gradual increase in biological complexity or the whole genome duplication (WGD) events experienced during plant evolution. Interestingly, most animal lineages have a relatively stable number of CuAOs, typically ranging from 3 to 4.
Table 1
| Taxonomy | Number of species | Number of CuAO | Average number of CuAO per species |
|---|---|---|---|
| Rhodophytes | 2 | 4 | 2.00 |
| Chlorophytes | 5 | 8 | 1.60 |
| Charophytes | 5 | 9 | 1.80 |
| hornworts | 7 | 9 | 1.29 |
| mosses | 8 | 27 | 3.38 |
| liverworts | 17 | 44 | 2.59 |
| Lycophytes | 7 | 22 | 3.14 |
| Ferns | 1 | 3 | 3.00 |
| Gymnosperms | 15 | 61 | 4.07 |
| Basal angiosperms | 3 | 6 | 2.00 |
| Magnoliids | 14 | 33 | 2.36 |
| Monocots | 30 | 167 | 5.57 |
| Eudicots | 74 | 557 | 7.53 |
| Brassicaceae | 36 | 273 | 7.58 |
| invertebrates | 15 | 53 | 3.53 |
| vertebrates | 64 | 211 | 3.30 |
The number of CuAO in plant species and animal species.
Evolutionary origin of CuAOs in plants and animals
To further explore the evolutionary origin of CuAOs, we selected representative species and constructed phylogenetic trees with IQ-TREE2 tool. The phylogenetic tree topologies distinctly categorized CuAOs from plants and animals into two separate clades, suggesting that the origin of CuAOs predates the divergence between plants and animals (Figure 1). Additionally, the phylogenetic tree shows that CuAOs underwent a distinct expansion in plants, forming two branches.
Figure 1
Phylogenetic classification of plant CuAOs
An IQ tree was constructed for CuAO proteins from plant lineages including rhodophytes, chlorophytes, charophytes, bryophytes (hornworts, mosses, liverworts), ferns, lycophytes, gymnosperms, basal angiosperms, angiosperms (Figure 2, Supplementary Figure S1). The phylogenetic tree shows that CuAO expanded and split into two branches in the most recent common ancestor of plants: clade I and clade II. Clade I is also known as the ζ branch and is a monophyletic branch, indicating that its members share a common ancestor. They evolved independently during radiation of plants. Clade II expanded after the division of vascular plants into the lycophytes and ferns, forming the subbranches of γ + δ and α + β. The differentiation of γ and δ occurred in the common ancestor of the mesangiosperms, while the differentiation of α and β occurred in the common ancestor of the eudicots and magnoliids. Clade I lost the green algae, and Clade II lost hornwort.
Figure 2
Phylogenetic classification of the animal CuAOs
To better understand the evolutionary relationships of animal CuAOs, we constructed an IQ tree using the CuAO sequences of 79 animals (15 invertebrates and 64 vertebrates) (Figure 3). In invertebrates, the CuAO from Coelenterata, Mollusca, and Leptocardii formed highly supported monophyletic clades. Vertebrates include mammals, birds, reptiles, amphibians and fish. The vertebrate CuAO can be divided into two main branches, AOC1 and AOC2-4. In AOC1, mammals and non-mammals form a monophyletic branch, while in AOC2-4, mammals and non-mammals are separated. Non-mammalian AOC2 and AOC3 form their own subbranches, while mammalian AOC2-4 further divide into the subbranches of AOC2 and AOC3/4. Interestingly, the differentiation of AOC3 and AOC4 in mammals occurs in species-specific ways, which may be due to different evolutionary pressures. AOC1 and AOC2 are present in all vertebrates, AOC3 is present in amphibians, reptiles, and mammals, and AOC4 is unique to mammals. CAO-like is present in invertebrates and reptiles.
Figure 3
Phylogenetic classification of CuAOs in Brassicaceae
Certain species within the Brassicaceae family are widely cultivated and possess considerable economic importance. To further elucidate the systematic relationships within the Brassicaceae CuAO family, we constructed a phylogenetic tree using 36 Brassicaceae plants and three sister groups (Carica papaya, Tarenaya hassleriana, Cleome violacea) with a total of 291 CuAOs (Figure 4). The five subbranches of the Brassicaceae are all monophyletic branches, and the α branch has a significantly higher number of CuAOs than the other subbranches (Figure 4A).
Figure 4
Previous studies have shown that the variable residues X1 and X2 in the T/S-X1-X2-N-Y-D active site motif of proteases follow the substrate preference of mammalian CuAOs, which helps distinguish their subfamilies (; ; ). This leads us to believe that the residues at positions X1 and X2 may help distinguishing between the Brassicaceae CuAO subfamilies. The results show that the Brassicaceae CuAOζ branch lacks the T/S-X1-X2-N-Y-D active site motif in its entirety, which may affect the function of the enzyme or its role in metabolic pathways, suggesting that this branch has different biological functions from other branches with active sites (Figure 4B). CuAOγ and CuAOδ tend to have the first amino acid in the active site motif more biased towards serine (S), while CuAOβ and CuAOα are more biased towards threonine (T). The active site motifs of CuAOγ and CuAOδ are highly conserved, while those of CuAOβ and CuAOα differ. In CuAOβ, the residue at position X1 is always leucine (L), while in CuAOα, the residue is more inclined towards valine (V). Therefore, the specific bias of residue X1 corresponds with the phylogenetic classification of the Brassicaceae CuAOα and CuAOβ.
Expansion of CuAOs during plant and animal evolution
In most species, tandem and segmental duplications serve as essential forces driving evolution (). To explore the expansion of CuAO genes in plants and animals, we conducted a synteny analysis (Figure 5, Supplementary Figure S2, Supplementary Table S3). 10, 12, 2, 4, 3, 2, 1, 16, 7, and 2 pairs of segmental duplication genes were identified in B. napus, G. max, M. domestica, M. acuminata, P. somniferum, P. vulgaris, P. patens, T. aestivum, P. virgatum, and X. laevis, respectively. Furthermore, tandem duplication resulted in 4, 3, 2, 2, 4, 2, 2, 2, 5, 2, 6, 2, 2, 6, 4, 5, 6, 3, 2, 6, 2, 5, 2, 2, and 2 additional CuAO genes in A. comosus, A. thaliana, C. papaya, G. biloba, G. max, M. domestica, O. sativa, P. somniferum, P. vulgaris, S. tuberosum, T. cacao, V. vinifera, P. virgatum, B. indicus, B. taurus, B. gargarizans, C. hircus, G. gorilla, H. sapiens, M. musculus, N. parkeri, N. scutatus, R. bivittatum, X. laevis, and X. tropicalis, respectively. Various mechanisms contribute to the expansion of CuAO genes in different plant species. Segmental duplication primarily drives this expansion in G. max, P. somniferum, and P. virgatum, while tandem duplication is the main mechanism for P. vulgaris. Segmental duplication occurs only in B. napus, M. acuminata, P. patens, and T. aestivum, whereas tandem duplication is exclusive to A. thaliana, A. comosus, C. papaya, G. biloba, O. sativa, S. tuberosum, T. cacao, and V. vinifera. Interestingly, in animal species, the vast majority only have tandem duplication. This emphasizes that different plant species adopt diverse evolutionary strategies to expand the CuAO genes, while animal species tend to have tandem duplication.
Figure 5
The analysis of Ka/Ks ratios (the rate of non-synonymous to synonymous substitutions) was performed to evaluate selection pressure on gene evolution (Zhang et al., 2024). A Ka/Ks ratio above 1 indicates positive selection, while a ratio below 1 suggests purifying or negative selection (Zhang et al., 2024). These results revealed that the Ka/Ks ratios for all CuAO paralogs were below 1, suggesting that they underwent purifying selection throughout their evolution (Supplementary Table S4).
To delve deeper into the evolutionary relationships of the CuAO gene, inter-species synteny analyses were conducted on 15 plant genomes and 12 animal genomes (Figure 6, Supplementary Table S5). The CuAO homologous genes of individual plants demonstrated one-to-one collinear relationships between C. papaya and A. thaliana, R. sativus and L. sativa, O. sativa and M. acuminata, M. acuminata and D. cayenensis (Figure 6A). In addition, either one-to-many or many-to-one homozygosity was identified between A. trichopoda and P. somniferum, P. somniferum and C. papaya, A. thaliana and B. napus, B. napus and B. oleracea, B. oleracea and B. rapa, B. rapa and R. sativus, L. sativa and A. chinensis, A. chinensis and G. max, G. max and T. aestivum, T. aestivum and O. sativa (Figure 6A). Interestingly, in vertebrates, AOC1 shows interspecies synteny only with AOC1, while AOC2, AOC3, AOC4, and AOCX show interspecies synteny among themselves (Figure 6B).
Figure 6
Expression patterns of CuAOs in key evolutionary lineages of plants and animals
To investigate the expression of CuAO genes, we performed a gene expression analysis of CuAOs in eight representative species from plants and animals: A. thaliana, B. napus, B. oleracea, G. max, O. sativa, H. sapiens, M. musculus, and P. troglodytes (Figure 7). In A. thaliana, we examined 23 different tissues and developmental stages (Figure 7A). AtCuAOδ, AtCuAOζ, AtCuAOγ1, AtCuAOα3, AtCuAOβ, and AtCuAOγ2 are widely expressed in various tissues, with AtCuAOζ showing significantly higher expression than the other CuAO genes. Interestingly, in various tissues, AtCuAOγ1 is generally expressed at higher levels than AtCuAOγ2, while AtCuAOα3 shows significantly higher expression than both AtCuAOα2 and AtCuAOα1. Notably, AtCuAOα1 is highly expressed specifically in siliques, whereas AtCuAOα2 is predominantly expressed in rosette leaves. In B. oleracea, BolCuAOδ2, BolCuAOζ2, and BolCuAOζ1 are widely expressed, with BolCuAOδ2 showing higher expression levels (Figure 7B). In leaves, siliques, and callus, the expression of BolCuAOζ2 exceeds that of BolCuAOζ1, while in flowers and buds, BolCuAOζ1 is more highly expressed than BolCuAOζ2. Additionally, BolCuAOα is specifically highly expressed in buds. B. napus and G. max have higher number of CuAOs, but their expression levels are generally low across all detected tissues (Figures 7C, D). BnaCuAOδ1 and BnaCuAOδ2 exhibit highly similar expression patterns, showing significantly higher specific expression than other genes in filaments, petals, sepals, leaves, seeds, and stem peels (Figure 7C). GmaCuAOα4 is specifically expressed at high levels in root hairs, roots, and nodules, while GmaCuAOα6 is specifically expressed at high levels in pods (Figure 7D). In rice, OsaCuAOζ2 is specifically highly expressed in the glumes, anthers, stems, and roots, while OsaCuAOδ is specifically highly expressed in the stems (Figure 7E). In summary, duplicated CuAO genes in Brassicaceae exhibit differing expression across various tissues, indicating functional divergence among these paralogs. In soybean and rice, their specific expression suggests roles in different growth stages or plant tissues, potentially involving development, growth, and nutrient metabolism.
Figure 7
Most animal lineages have a relatively few CuAOs, usually 3 to 4. AOC2 and AOC3 are widely expressed in most tissues (Figures 7F-H). In humans and chimpanzees, AOC3 is expressed significantly more than AOC2 (Figures 7F, G), while in mice, MumAOC2 is expressed significantly more than MumAOC3 (Figure 7H). The expression of HsaAOC1 is tissue-specific, with high specific expression in the rectum, adult mammalian kidney, and metanephros cortex in humans and in the adult mammalian kidney in chimpanzees (Figure 7F-G). Additionally, PtrAOC4 is specifically highly expressed in the liver of chimpanzees.
Discussion
This study conducted phylogenetic classification and examination of the characteristics of thousands of CuAOs in hundreds of plant and animal species, deepening our understanding of CuAOs in plants and animals, particularly those in angiosperms and vertebrates. The phylogenetic perspective aids in uncovering the molecular and biological roles of different CuAO proteins (Figure 8).
Figure 8
Phylogenetic relationship of CuAOs in plant and animal
Since genome-based identification of CuAO proteins might overlook certain proteins due to the quality of genome assembly and annotation accuracy, we incorporated numerous species within each phylogenetic group to mitigate the effects of any potentially missing proteins. Our search for the CuAO proteins revealed that the CuAO gene family is widely distributed across both plant and animal lineages, with varying copy numbers in different plant species and a more conservative distribution in animal species, with approximately 3-4 copies per species (Table 1). Our study suggests that the origin of CuAOs predates the divergence of plants and animals, implying that the bioammonia degradation pathway was established in the common ancestor of plants and animals (Figure 1) (; ). After the divergence of plant and animal species, CuAOs underwent significant expansion in plants, forming two branches (Clade I and Clade II) (Figure 2). It is worth noting that, consistent with previous research results, we identified a CuAO gene located in the ζ branch in Spirodela polyrhiza (). Interestingly, in the ζ branch, most monocots have multiple copies of the CuAO gene, while a few monocots (S. polyrhiza, A. americanus, K. littledalei, Z. mays, B. distachyon) contain only a single copy (Supplementary Figure S1). Previous studies have shown that mammalian CuAOs can be classified based on the two residues X1 and X2 in the active site motif T/S-X1-X2-N-Y-D (). Our study shows that Clade I is a monophyletic branch that has lost the T/S-X1-X2-N-Y-D active site motif, which may affect the function of the enzyme or its role in the metabolic pathway, suggesting that Clade I is functionally diverse from Clade II (Figures 2, 4). In previous studies, the plant CuAOs were divided into three subfamilies (clades I-III), clade I consists of the α and β subgroups, clade II consists of the γ and δ subgroups, and clade III consists of the ζ subgroup (). In our study, Clade II further expanded during evolution and formed four major branches: α, β, γ, and δ. CuAOγ and CuAOδ tend to have serine (S) as the first amino acid in the active site motif, while CuAOβ and CuAOα are more inclined towards threonine (T) (Figure 4B). The active site motifs of CuAOγ and CuAOδ are highly conserved, but those of CuAOβ and CuAOα differ. In CuAOβ, the residue at position X1 is always leucine (L), while in CuAOα, it is more likely to be valine (V). Thus, the specific bias of residue X1 correlates with the phylogenetic classification of CuAOα and CuAOβ (Figure 4).
In the phylogenetic analysis of vertebrate CuAOs, both mammalian and non-mammalian AOC1 proteins are grouped in the same clade, showing a robust bootstrap support of 100% (Figure 3). In contrast to AOC1s, non-mammalian AOC2-4 proteins deviate from their mammalian equivalents with a bootstrap value of 100% as well (Figure 3). Non-mammalian AOC2 and AOC3 proteins form distinct sub-branches, while mammalian AOC2-4 proteins are further classified into separate sub-branches for AOC2 and a combined group for AOC3/AOC4, both presenting high bootstrap values. These results are in line with previous studies (). Interestingly, the differentiation between mammalian AOC3 and AOC4 proteins occurs in a species-specific manner, perhaps due to different evolutionary pressures. Furthermore, unlike previous studies, our research showed that cartilaginous fish have not only AOC1 but also AOC2 ().
Insights into CuAO evolution and functional diversification
To understand the evolutionary implications of the phylogenetic architecture of the CuAO family, it is optimal to cross-reference the phylogenetic relationships of CuAOs and their functions. Repeated events frequently cause the expansion and functional diversification of gene families (). Generally, CuAOs within the same clade have similar and specific functions. CuAOδ is highly expressed in flowers, and this is the case in A. thaliana, B. oleracea, B. napus, G. max and O. sativa (Figure 7A-E). Studies have shown that Malus domestica CuAOδ (MdAO2) contributes to the formation of flower fragance (Zarei et al., 2015). Furthermore, different members of the same clade might have functional specialization or redundancy. For example, the expression of AtCuAOα3 is significantly higher than that of AtCuAOα1 and AtCuAOα2, and the expression of AtCuAOγ1 is much higher than that of AtCuAOγ2; while the expression patterns of BolCuAOζ1 and BolCuAOζ2 are alike, and so are those of BnaCuAOδ1 and BnaCuAOδ2 (Figure 7A-E). SlCuAO1 (SlyCuAOβ3), SlCuAO2 (SlyCuAOζ1), SlCuAO3 (SlyCuAOζ2), SlCuAO4 (SlyCuAOζ3), SlCuAO6 (SlyCuAOβ2) and SlCuAO7 (SlyCuAOβ1) in S. lycopersicum are mainly specifically expressed in roots, while SlCuAO5 (SlyCuAOδ) is specifically expressed in fruits, and SlCuAO8 (SlyCuAOα) is mainly expressed in flowers (). The CuAO gene family in tomato shows obvious tissue-specific expression: the β and ζ subtypes are mainly expressed in roots, the δ subtype is enriched in fruits, and the α subtype is highly expressed in flowers. This branch-specific expression pattern reflects the functional differentiation of the CuAO gene family.
Compared with A. thaliana CuAO, human, mouse, and chimpanzee CuAO isoforms show significant functional differences. For example, in humans and chimpanzees, the expression level of AOC3 is significantly higher than that of AOC2 and AOC1, while in mice, the expression level of AOC2 is significantly higher than that of AOC3 and AOC1 (Figure 7F-H). The expression of AOC1 is tissue-specific and is highly expressed in the kidneys (). AOC1 gene encodes diamine oxidase (), and human AOC1 is thought to have in vivo functions in cell proliferation, inflammation, allergic reactions, and ischemia (). AOC2 gene encodes retinal-specific amine oxidase, and human AOC2 is thought to play a role in hereditary retinal diseases (). AOC3 gene encodes vascular adhesion protein-1, and human AOC3 expression is upregulated during inflammation and plays an important role in the migration of white blood cells from the blood into tissues (; ; ). These difference suggest different functions of these three isoforms.
Materials and methods
Data sources and sequence acquisition
We selected 267 species and obtained their genome and proteome sequences from public databases like NCBI (http://www.ncbi.nlm.nih.gov/), Ensembl Plants (http://plants.ensembl.org/index.html), and BRAD (http://brassicadb.cn/#/) (Supplementary Table S1). To identify CuAO genes in animals and plants, we first performed an HMMER search (E-value = 1e-10) using the Hidden Markov Model profile of the CuAO domain (PF01179) in local databases (; ). The HMMER parameters are set as follows: hmmsearch -E 1e-10 PF01179.hmm protein.fa. We also used BLASTP with A. thaliana CuAO amino acid sequences against the protein database, setting an E-value threshold of less than 1e-6. The putative CuAOs were further validated using online tools CDD, HMM, and SMART (; ).
Multiple sequence alignment (MSA) and phylogenetic analysis
Multiple sequence alignments of CuAO were performed using Muscle5 (). Phylogenetic trees were generated from the full-length protein sequences of CuAO. The maximum likelihood (ML) tree was constructed with IQ-TREE2 using parameters ‘-m MFP -B 1000’ and 1000 bootstrap replicates (). The consensus tree topology was visualized using iTOL ().
Synteny analysis
To identify homologous pairs in a specific species, we used the all-to-all BLASTP method. Syntenic blocks were inferred using MCScanX with default parameters, including an E-value threshold of 1e-10 and a minimum of 5 BLAST hits (). The resulting gene arrangement was visualized with CIRCOS, linking potential duplicated genes to show their connections (). The duplicated CuAO gene pairs were subjected to TBtools software for calculation of synonymous (Ka) and non-synonymous (Ks) substitution rates ().
Expression analysis of CuAO genes
Expression data for ten species were downloaded from public databases, including A. thaliana (The Arabidopsis Information Resource), B. napus (https://yanglab.hzau.edu.cn/BnIR), B. oleracea (http://bogdb.com/), G. max (https://www.soybase.org/), O. sativa (http://rice.uga.edu/index.shtml), H. sapiens (https://www.bgee.org/), M. musculus (https://www.bgee.org/), and P. troglodytes (https://www.bgee.org/) (Supplementary Table S6). Heat map was established using the R package “pheatmap” (https://github.com/raivokolde/pheatmap).
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author/s.
Author contributions
ZZ: Funding acquisition, Writing – original draft. TX: Software, Visualization, Writing – original draft. KL: Writing – original draft. KH: Writing – original draft. SW: Software, Visualization, Writing – review & editing. LW: Software, Visualization, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the grants from National Natural Science Foundation of China (32170351), the General Research Projects of Zhejiang Provincial Department of Education (Y202351039), Huzhou Science and Technology Plan Project (2023GZ44), and Research Program of Huzhou College (2023HXKM09).
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.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
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/fpls.2025.1544527/full#supplementary-material
Supplementary Table 1List of complete predicted proteomes databases used in this study.
Supplementary Table 2The identified CuAO protein.
Supplementary Table 3Segmental duplication and tandem gene.
Supplementary Table 4The Ka/Ks analysis of CuAO genes.
Supplementary Table 5Segmental duplication gene pairs of CuAOs in interspecies.
Supplementary Table 6Source of expression data.
Supplementary Figure 1Phylogenetic relationship of representative CuAO genes from plants. The phylogenetic tree was built using IQ-TREE2. Distinct branch colors represent various plant lineages.
Supplementary Figure 2Intraspecies syntenic relationships of CuAO genes in representative plants. The syntenic paralogs of CuAO genes are connected by black lines. The asterisk indicates tandem duplicate pairs.
References
1
BoomsmaF.BhaggoeU. M.van der HouwenA. M.van den MeirackerA. H. (2003). Plasma semicarbazide-sensitive amine oxidase in human (patho)physiology. Biochim. Biophys. Acta1647, 48–54. doi: 10.1016/S1570-9639(03)00047-5
2
BourS.DaviaudD.GresS.LefortC.PrévotD.ZorzanoA.et al. (2007). Adipogenesis-related increase of semicarbazide-sensitive amine oxidase and monoamine oxidase in human adipocytes. Biochimie89, 916–925. doi: 10.1016/j.biochi.2007.02.013
3
BrazeauB. J.JohnsonB. J.WilmotC. M. (2004). Copper-containing amine oxidases. Biogenesis and catalysis; a structural perspective. Arch. Biochem. Biophys.428, 22–31. doi: 10.1016/j.abb.2004.03.034
4
ChenH.WangT.HeX.CaiX.LinR.LiangJ.et al. (2022). BRAD V3.0: an upgraded Brassicaceae database. Nucleic Acids Res.50, D1432–d1441. doi: 10.1093/nar/gkab1057
5
ChenC.WuY.LiJ.WangX.ZengZ.XuJ.et al. (2023). TBtools-II: A “one for all, all for one” bioinformatics platform for biological big-data mining. Mol. Plant16, 1733–1742. doi: 10.1016/j.molp.2023.09.010
6
DunkelP.GelainA.BarloccoD.HaiderN.GyiresK.SperlághB.et al. (2008). Semicarbazide-sensitive amine oxidase/vascular adhesion protein 1: recent developments concerning substrates and inhibitors of a promising therapeutic target. Curr. Med. Chem.15, 1827–1839. doi: 10.2174/092986708785133022
7
EdgarR. C. (2022). Muscle5: High-accuracy alignment ensembles enable unbiased assessments of sequence homology and phylogeny. Nat. Commun.13, 6968. doi: 10.1038/s41467-022-34630-w
8
ElmoreB. O.BollingerJ. A.DooleyD. M. (2002). Human kidney diamine oxidase: heterologous expression, purification, and characterization. J. Biol. Inorganic. Chem.7, 565–579. doi: 10.1007/s00775-001-0331-1
9
FinneyJ.MoonH. J.RonnebaumT.LantzM.MureM. (2014). Human copper-dependent amine oxidases. Arch. Biochem. Biophys.546, 19–32. doi: 10.1016/j.abb.2013.12.022
10
FortesA. M.Agudelo-RomeroP.PimentelD.AlkanN. (2019). Transcriptional modulation of polyamine metabolism in fruit species under abiotic and biotic stress. Front. Plant Sci.10, 816. doi: 10.3389/fpls.2019.00816
11
FraudentaliI.GhugeS. A.CarucciA.TavladorakiP.AngeliniR.ConaA.et al. (2019). The copper amine oxidase atCuAOδ Participates in abscisic acid-induced stomatal closure in arabidopsis. Plants (Basel Switzerland).8, 1–15. doi: 10.3390/plants8060183
12
FraudentaliI.GhugeS. A.CarucciA.TavladorakiP.AngeliniR.Rodrigues-PousadaR. A.et al. (2020a). Developmental, hormone- and stress-modulated expression profiles of four members of the Arabidopsis copper-amine oxidase gene family. Plant Physiol. Biochem.: PPB147, 141–160. doi: 10.1016/j.plaphy.2019.11.037
13
FraudentaliI.PedalinoC.D’IncàR.TavladorakiP.AngeliniR.ConaA. (2023). Distinct role of AtCuAOβ- and RBOHD-driven H(2)O(2) production in wound-induced local and systemic leaf-to-leaf and root-to-leaf stomatal closure. Front. Plant Sci.14, 1154431. doi: 10.3389/fpls.2023.1154431
14
FraudentaliI.Rodrigues-PousadaR. A.AngeliniR.GhugeS. A.ConaA. (2021). Plant copper amine oxidases: key players in hormone signaling leading to stress-induced phenotypic plasticity. Int. J. Mol. Sci.22, 1–15. doi: 10.3390/ijms22105136
15
FraudentaliI.Rodrigues-PousadaR. A.TavladorakiP.AngeliniR.ConaA. (2020b). Leaf-wounding long-distance signaling targets atCuAOβ Leading to root phenotypic plasticity. Plants (Basel Switzerland).9, 1–15. doi: 10.3390/plants9020249
16
GhugeS. A.CarucciA.Rodrigues-PousadaR. A.TisiA.FranchiS.TavladorakiP.et al. (2015a). The apoplastic copper AMINE OXIDASE1 mediates jasmonic acid-induced protoxylem differentiation in arabidopsis roots. Plant Physiol.168, 690–707. doi: 10.1104/pp.15.00121
17
GhugeS. A.CarucciA.Rodrigues-PousadaR. A.TisiA.FranchiS.TavladorakiP.et al. (2015b). The MeJA-inducible copper amine oxidase AtAO1 is expressed in xylem tissue and guard cells. Plant Signaling Behav.10, e1073872. doi: 10.1080/15592324.2015.1073872
18
GroßF.RudolfE. E.ThieleB.DurnerJ.AstierJ. (2017). Copper amine oxidase 8 regulates arginine-dependent nitric oxide production in Arabidopsis thaliana. J. Exp. Bot.68, 2149–2162. doi: 10.1093/jxb/erx105
19
HeniquezA.MeissonnierG.VisentinV.PrévotD.CarpénéC. (2003). High expression of semicarbazide-sensitive amine oxidase genes AOC2 and AOC3, but not the diamine oxidase gene AOC1 in human adipocytes. Inflammation Res.52 Suppl 1, S74–S75. doi: 10.1007/s000110300061
20
ImamuraY.KubotaR.WangY.AsakawaS.KudohJ.MashimaY.et al. (1997). Human retina-specific amine oxidase (RAO): cDNA cloning, tissue expression, and chromosomal mapping. Genomics40, 277–283. doi: 10.1006/geno.1996.4570
21
JalkanenS.KarikoskiM.MercierN.KoskinenK.HenttinenT.ElimaK.et al. (2007). The oxidase activity of vascular adhesion protein-1 (VAP-1) induces endothelial E- and P-selectins and leukocyte binding. Blood110, 1864–1870. doi: 10.1182/blood-2007-01-069674
22
KaitaniemiS.ElovaaraH.GrönK.KidronH.LiukkonenJ.SalminenT.et al. (2009). The unique substrate specificity of human AOC2, a semicarbazide-sensitive amine oxidase. Cell. Mol. Life Sci.: CMLS.66, 2743–2757. doi: 10.1007/s00018-009-0076-5
23
KiviE.ElimaK.AaltoK.NymalmY.AuvinenK.KoivunenE.et al. (2009). Human Siglec-10 can bind to vascular adhesion protein-1 and serves as its substrate. Blood114, 5385–5392. doi: 10.1182/blood-2009-04-219253
24
KrzywinskiM.ScheinJ.BirolI.ConnorsJ.GascoyneR.HorsmanD.et al. (2009). Circos: an information aesthetic for comparative genomics. Genome Res.19, 1639–1645. doi: 10.1101/gr.092759.109
25
KurkijärviR.AdamsD. H.LeinoR.MöttönenT.JalkanenS.SalmiM. (1998). Circulating form of human vascular adhesion protein-1 (VAP-1): increased serum levels in inflammatory liver diseases. J. Immunol. (Baltimore. Md.: 1950).161, 1549–1157.
26
LetunicI.BorkP. (2024). Interactive Tree of Life (iTOL) v6: recent updates to the phylogenetic tree display and annotation tool. Nucleic Acids Res. 52, W78–W82. doi: 10.1093/nar/gkae268
27
LetunicI.KhedkarS.BorkP. (2021). SMART: recent updates, new developments and status in 2020. Nucleic Acids Res.49, D458–d460. doi: 10.1093/nar/gkaa937
28
Lopes de CarvalhoL.Bligt-LindénE.RamaiahA.JohnsonM. S.SalminenT. A. (2019). Evolution and functional classification of mammalian copper amine oxidases. Mol. Phylogenet. Evol.139, 106571. doi: 10.1016/j.ympev.2019.106571
29
MaintzL.NovakN. (2007). Histamine and histamine intolerance. Am. J. Clin. Nutr.85, 1185–1196. doi: 10.1093/ajcn/85.5.1185
30
MaintzL.SchwarzerV.BieberT.van der VenK.NovakN. (2008). Effects of histamine and diamine oxidase activities on pregnancy: a critical review. Hum. Reprod. Update14, 485–495. doi: 10.1093/humupd/dmn014
31
Marchler-BauerA.DerbyshireM. K.GonzalesN. R.LuS.ChitsazF.GeerL. Y.et al. (2015). Hurwitz DI et al: CDD: NCBI’s conserved domain database. Nucleic Acids Res.43, D222–D226. doi: 10.1093/nar/gku1221
32
McGrathA. P.HilmerK. M.CollyerC. A.ShepardE. M.ElmoreB. O.BrownD. E.et al. (2009). Structure and inhibition of human diamine oxidase. Biochemistry48, 9810–9822. doi: 10.1021/bi9014192
33
MinhB. Q.SchmidtH. A.ChernomorO.SchrempfD.WoodhamsM. D.von HaeselerA.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
34
MurakawaT.BabaS.KawanoY.HayashiH.YanoT.KumasakaT.et al. (2019). In crystallo thermodynamic analysis of conformational change of the topaquinone cofactor in bacterial copper amine oxidase. Proc. Natl. Acad. Sci. United. States America116, 135–140. doi: 10.1073/pnas.1811837116
35
Murray-StewartT. R.WosterP. M.CaseroR. A.Jr. (2016). Targeting polyamine metabolism for cancer therapy and prevention. Biochem. J.473, 2937–2953. doi: 10.1042/BCJ20160383
36
PengW.LiW.SongN.TangZ.LiuJ.WangY.et al. (2021). Genome-wide characterization, evolution, and expression profile analysis of GATA transcription factors in brachypodium distachyon. Int. J. Mol. Sci.22, 1–14. doi: 10.3390/ijms22042026
37
PengY.ZhaoK.ZhengR.ChenJ.ZhuX.XieK.et al. (2024). A comprehensive analysis of auxin response factor gene family in melastoma dodecandrum genome. Int. J. Mol. Sci.25, 1–18. doi: 10.3390/ijms25020806
38
PietrangeliP.NoceraS.MondoviB.MorpurgoL. (2003). Is the catalytic mechanism of bacteria, plant, and mammal copper-TPQ amine oxidases identical? Biochim. Biophys. Acta1647, 152–156. doi: 10.1016/s1570-9639(03)00083-9
39
Planas-PortellJ.GallartM.TiburcioA. F.AltabellaT. (2013). Copper-containing amine oxidases contribute to terminal polyamine oxidation in peroxisomes and apoplast of Arabidopsis thaliana. BMC Plant Biol.13, 109. doi: 10.1186/1471-2229-13-109
40
PotterS. C.LucianiA.EddyS. R.ParkY.LopezR.FinnR. D. (2018). HMMER web server: 2018 update. Nucleic Acids Res.46, W200–w204. doi: 10.1093/nar/gky448
41
QuY.AnZ.ZhuangB.JingW.ZhangQ.ZhangW. (2014). Copper amine oxidase and phospholipase D act independently in abscisic acid (ABA)-induced stomatal closure in Vicia faba and Arabidopsis. J. Plant Res.127, 533–544. doi: 10.1007/s10265-014-0633-3
42
SalmiM.JalkanenS. (2019). Vascular adhesion protein-1: A cell surface amine oxidase in translation. Antioxid. Redox Signaling30, 314–332. doi: 10.1089/ars.2017.7418
43
SalmiM.KalimoK.JalkanenS. (1993). Induction and function of vascular adhesion protein-1 at sites of inflammation. J. Exp. Med.178, 2255–2260. doi: 10.1084/jem.178.6.2255
44
SalminenT. A.SmithD. J.JalkanenS.JohnsonM. S. (1998). Structural model of the catalytic domain of an enzyme with cell adhesion activity: human vascular adhesion protein-1 (HVAP-1) D4 domain is an amine oxidase. Protein Eng.11, 1195–1204. doi: 10.1093/protein/11.12.1195
45
SchwelbergerH. G. (2006). Origins of plasma amine oxidases in different mammalian species. Inflammation Res.55 Suppl 1, S57–S58. doi: 10.1007/s00011-005-0041-1
46
ShuklaV.FatimaT.GoyalR. K.HandaA. K.MattooA. K. (2020). Engineered ripening-specific accumulation of polyamines spermidine and spermine in tomato fruit upregulates clustered C/D box snoRNA gene transcripts in concert with ribosomal RNA biogenesis in the red ripe fruit. Plants (Basel Switzerland).9, 1–20. doi: 10.3390/plants9121710
47
SmithD. J.SalmiM.BonoP.HellmanJ.LeuT.JalkanenS. (1998). Cloning of vascular adhesion protein 1 reveals a novel multifunctional adhesion molecule. J. Exp. Med.188, 17–27. doi: 10.1084/jem.188.1.17
48
TavladorakiP.ConaA.AngeliniR. (2016). Copper-containing amine oxidases and FAD-dependent polyamine oxidases are key players in plant tissue differentiation and organ development. Front. Plant Sci.7, 824. doi: 10.3389/fpls.2016.00824
49
TavladorakiP.ConaA.FedericoR.TemperaG.ViceconteN.SaccoccioS.et al. (2012). Polyamine catabolism: target for antiproliferative therapies in animals and stress tolerance strategies in plants. Amino Acids42, 411–426. doi: 10.1007/s00726-011-1012-1
50
UpadhyayR. K.ShaoJ.MaulJ. E.SchombergH.HandaA. K.RobertsD. P.et al. (2024a). Unlocking the role of novel primary/di-amine oxidases in crop improvement: Tissue specificity leads to specific roles connected to abiotic stress, hormone responses and sensing nitrogen. J. Plant Physiol.303, 154374. doi: 10.1016/j.jplph.2024.154374
51
UpadhyayR. K.ShaoJ.RobertsG. E.MattooA. K. (2024b). Comparative genomics and evidence for an unusual polyamine oxidation pathway in aquatic duckweed (Spirodela polyrhiza L.). Curr. Plant Biol.39, 100359. doi: 10.1016/j.cpb.2024.100359
52
WangW.PaschalidisK.FengJ. C.SongJ.LiuJ. H. (2019). Polyamine catabolism in plants: A universal process with diverse functions. Front. Plant Sci.10, 561. doi: 10.3389/fpls.2019.00561
53
WangY.TangH.DebarryJ. D.TanX.LiJ.WangX.et al. (2012). MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res.40, e49. doi: 10.1093/nar/gkr1293
54
ZareiA.TrobacherC. P.CookeA. R.MeyersA. J.HallJ. C.ShelpB. J. (2015). Apple fruit copper amine oxidase isoforms: peroxisomal MdAO1 prefers diamines as substrates, whereas extracellular MdAO2 exclusively utilizes monoamines. Plant Cell Physiol.56, 137–147. doi: 10.1093/pcp/pcu155
55
ZhangZ. B.XiongT.WangX. J.ChenY. R.WangJ. L.GuoC. L.et al. (2024). Lineage-specific gene duplication and expansion of DUF1216 gene family in Brassicaceae. PloS One19, e0302292. doi: 10.1371/journal.pone.0302292
Summary
Keywords
CuAO, molecular evolution, origin, diversification, gene duplication
Citation
Zhang Z, Xiong T, Li K, Huang K, Wu S and Wu L (2025) The tree of life of copper-containing amine oxidases. Front. Plant Sci. 16:1544527. doi: 10.3389/fpls.2025.1544527
Received
13 December 2024
Accepted
24 March 2025
Published
24 April 2025
Volume
16 - 2025
Edited by
Laura Emma Maria Morello, National Research Council (CNR), Italy
Reviewed by
Rakesh K. Upadhyay, Bowie State University, United States
Ilaria Fraudentali, Roma Tre University, Italy
Updates
Copyright
© 2025 Zhang, Xiong, Li, Huang, Wu and Wu.
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: Zaibao Zhang, zaibaozhang79@163.com
Disclaimer
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.