Molecular evolution and functional characterisation of an ancient phenylalanine ammonia-lyase gene (NnPAL1) from Nelumbo nucifera: novel insight into the evolution of the PAL family in angiosperms
© Wu et al.; licensee BioMed Central Ltd. 2014
Received: 24 March 2014
Accepted: 28 April 2014
Published: 9 May 2014
Phenylalanine ammonia-lyase (PAL; E.C.126.96.36.199) is a key enzyme of the phenylpropanoid pathway in plant development, and it catalyses the deamination of phenylalanine to trans-cinnamic acid, leading to the production of secondary metabolites. This enzyme has been identified in many organisms, ranging from prokaryotes to higher plants. Because Nelumbo nucifera is a basal dicot rich in many secondary metabolites, it is a suitable candidate for research on the phenylpropanoid pathway.
Three PAL members, NnPAL1, NnPAL2 and NnPAL3, have been identified in N. nucifera using genome-wide analysis. NnPAL1 contains two introns; however, both NnPAL2 and NnPAL3 have only one intron. Molecular and evolutionary analysis of NnPAL1 confirms that it is an ancient PAL member of the angiosperms and may have a different origin. However, PAL clusters, except NnPAL1, are monophyletic after the split between dicots and monocots. These observations suggest that duplication events remain an important occurrence in the evolution of the PAL gene family. Molecular assays demonstrate that the mRNA of the NnPAL1 gene is 2343 bp in size and encodes a 717 amino acid polypeptide. The optimal pH and temperature of the recombinant NnPAL1 protein are 9.0 and 55°C, respectively. The NnPAL1 protein retains both PAL and weak TAL catalytic activities with Km values of 1.07 mM for L-phenylalanine and 3.43 mM for L-tyrosine, respectively. Cis-elements response to environmental stress are identified and confirmed using real-time PCR for treatments with abscisic acid (ABA), indoleacetic acid (IAA), ultraviolet light, Neurospora crassa (fungi) and drought.
We conclude that the angiosperm PAL genes are not derived from a single gene in an ancestral angiosperm genome; therefore, there may be another ancestral duplication and vertical inheritance from the gymnosperms. The different evolutionary histories for PAL genes in angiosperms suggest different mechanisms of functional regulation. The expression patterns of NnPAL1 in response to stress may be necessary for the survival of N. nucifera since the Cretaceous Period. The discovery and characterisation of the ancient NnPAL1 help to elucidate PAL evolution in angiosperms.
KeywordsNelumbo nucifera Phenylalanine ammonia-lyase gene Evolution Expression Bioinformatics analysis
The phenylpropanoid pathway is an important branch of the plant secondary metabolism pathways that produces many essential secondary metabolites. In this pathway, secondary metabolic products, such as lignin, flavonoids and coumarins, play important roles in plant growth, development, mechanical support, and disease resistance [1, 2]. Phenylalanine ammonia-lyase (PAL; E.C.188.8.131.52) is the first and key enzyme between primary and secondary metabolism, and it catalyses the biotransformation of L-phenylalanine to trans-cinnamic acid. The synthesis of many secondary metabolites, such as flavonoids, flavonols, anthocyanins, condensed tannins, lignins, coumarins, and ubiquinone occur downstream of the phenylpropanoid pathway, [3–6] and is controlled by PAL. Koukol and Conn reported the first plant PAL in 1961. Currently, it is known that the PAL is widely found in all higher plants, a few fungi, and a single prokaryote, Streptomyces, but not animals . Furthermore, PAL shows potential to treat human phenylketonuria, an inborn error of phenylalanine metabolism . Several studies [9–12] have shown that the PALs from Rhodotorula photosynthetic bacteria and monocot plants also utilise tyrosine in addition to phenylalanine; however, the dicot PALs only utilise Phe efficiently. During the past four decades, many PAL genes have been cloned and studied from various plants, such as Ginkgo biloba, Ephedra sinica, Oryza sativa, Isatis indigotica, Arabidopsis thaliana, Jatropha curcas, and Lycoris radiate, and the first crystal structure of a plant PAL was determined from parsley (Petroselinum crispum) . PAL exists as a small multigene family, consisting of 2–6 members; however, some species contain additional member, such as potato (~40 copies)  and tomato (~26 copies) . During the evolution of higher plants, the plant PAL genes diversified into various functions in each species, such as Arabidopsis thaliana. Another important ammonia lyase, histidine ammonia-lyase (HAL), is found in prokaryotes and animals and plays roles in the general histidine degradation pathway. The crystal structure of HAL from Pseudomonas putida revealed its catalytic mechanism of novel polypeptide modification . Despite large differences in the primary sequence of proteins, PAL functions as a tetramer, similar to HAL in vivo. Presumably, PAL developed from HAL when fungi and plants diverged from the other kingdoms [7, 25].
Nelumbo nucifera (Nelumbo, Nelumbonaceae) (2n = 16) is a perennial aquatic plant with ornamental flowers of medicinal and phylogenetic importance. N. nucifera produces a series of important secondary metabolites, including alkaloids, flavonoids, steroids, triterpenoids, glycosides and polyphenols . The N. nucifera secondary metabolites have a wide range of medical functions and also play important roles in the response to environmental stress, such as pathogen attack and ultraviolet damage. For example, it has been reported that benzylisoquinoline alkaloids and flavonoids from the leaves of N. nucifera are a potential candidate for HIV therapy . Nelumbo has survived since the Late Cretaceous, along with a number of other relicts, including Ginkgo, Sequoia, Metasequoia, and Liriodendron. It remains to be determined the mechanism by which PAL evolution has allowed N. nucifera to adapt to harsh environmental stress. Along with the N. nucifera genome project [29, 30], high-throughput sequencing data will provide a foundation for identifying the key genes in metabolic pathways. However, related research for N. nucifera is very limited.
In this study, three intact PAL genes in N. nucifera, NnPAL1, NnPAL2 and NnPAL3 are identified by genome-wide analysis. NnPAL1 is an ancient PAL member in angiosperms. The objective of this study is to determine the evolutionary origin, gene structure, function, and expression patterns of this gene under stress conditions.
Genomic identification and exon/intron structure analysis of the PAL gene family in N. nucifera
Using BLASTP to search the protein database in NCBI, we found that NnPAL1 is more similar to the PAL genes of gymnosperms (73% identity to GbPAL, ABU49842.1; 72% identity to PmPAL, ACS28225.2; and 69% identity to EsPAL, BAG74771.1) than dicots (63% BnPAL, ABC69916.1; 64% AtPAL NP_181241.1; and 63% DcPAL, BAC56977.1) (Additional file 2: Figure S2). This is contrary to the phylogeny of N. nucifera in the plant kingdom . However, the deduced NnPAL1 protein has the same nine strictly conserved residues, Y112,L140,S204,N260,Q348,Y351,R354,F400, and Q488, that are found in PcPAL of Petroselinum crispum. A typical phenylalanine and histidine ammonia-lyase signature (GTITASGDLVPLSYIA) also exists at position 199–214 (Additional file 3: Figure S3).
Evolutionary analysis of NnPAL1 in N. nucifera
Identification of the PAL gene family from the Phytozome database in this study
PAL seq ID
Physcomitrella patens (10)
Pp1s22_3V6.1 Pp1s32_311V6.1 Pp1s36_253V6.1 Pp1s5_197V6.1 Pp1s20_305V6.1 Pp1s494_3V6.1 Pp1s500_4V6.1 Pp1s43_88V6.1 Pp1s43_67V6.1 Pp1s52_44V6.1
Pinus taeda (5)
Pteda1143311 Pteda17307 Pteda9006 Pteda28316 Pteda34319
Bradi5g15830.1 Bradi3g48840.1 Bradi3g47120.1 Bradi3g49260.1 Bradi3g49270.1 Bradi3g49250.2
Oryza sativa (9)
Os02g41630.2 Os02g41650.1 Os02g41670.1 Os02g41680.1 Os04g43760.1 Os04g43800.1 Os05g35290.1 Os11g48110.1 Os12g33610.1
Sb04g026520.1 Sb04g026560.1 Sb04g026510.1 Sb06g022740.1 Sb06g022750.1 Sb01g014020.1
Setaria italica (6)
Si016478m Si016504m Si016467m Si009345m Si009509m Si012256m
Zea mays (11)
GRMZM2G441347_T01 GRMZM2G118345_T01 GRMZM2G447436_T01 GRMZM2G063917_T01 GRMZM2G160541_T01 GRMZM2G081582_T01 GRMZM2G326335_T01 GRMZM2G334660_T01 GRMZM2G170692_T01 GRMZM2G074604_T01 GRMZM2G029048_T01
Nelumbo nucifera (3)
NnPAL1 NnPAL2 NnPAL3
Aquilegia coerulea (2)
Arabidopsis thaliana (4)
AT2G37040.1(AtPAL1) AT3G53260.1 (AtPAL2) AT5G04230.2 (AtPAL3) AT3G10340.1 (AtPAL4)
Cucumis sativus (8)
Cucsa.124460.1 Cucsa.124480.1 Cucsa.124470.1 Cucsa.124500.1 Cucsa.124490.1 Cucsa.385970.1 Cucsa.124510.1 Cucsa.137590.1
Glycine max (5)
Glyma13g20800.1 Glyma03g33880.1 Glyma02g47940.1 Glyma20g32135.1 Glyma10g35381.1
Mimulus guttatus (3)
mgv1a001582m mgv1a019971m mgv1a002104m
Populus trichocarpa (4)
Potri.010G224200.1 Potri.010G224100.1 Potri.006G126800.1 Potri.008G038200.1
Vitis vinifera (10)
GSVIVT01024306001 GSVIVT01016257001 GSVIVT01024292001 GSVIVT01025214001 GSVIVT01024294001 GSVIVT01024305001 GSVIVT01024315001 GSVIVT01025703001 GSVIVT01024295001 GSVIVT01024293001
Isolation and bioinformatics characterisation of the full-length NnPAL1 cDNA in N. nucifera
NnPAL1 has a unique gene structure and phylogenetic position. To determine whether NnPAL1 became a pseudogene during evolution, isolation of the full-length NnPAL1 cDNA is performed from the transcripts of tender leaves. The partial cDNA is obtained by DOP-PCR with degenerate primers. A full-length cDNA containing an open reading frame of 2151 bp is then produced using 5′-RACE and 3′-RACE. Using BLASTN to search the whole genome sequence of N. nucifera, we confirm that the newly cloned ancient PAL gene exactly matched the NnPAL1 sequence.
Based on the crystal structure of PcPAL (1 W27), the SWISS-MODEL software is used to predict the three-dimensional structure of the NnPAL1 protein. The result indicate that NnPAL1 comprises an MIO domain, core domain and shielding domain . Moreover, a highly conserved Ala-Ser-Gly triad  that can be converted autocatalytically is also identified within NnPAL1 (Figure 3B). The results of the bioinformatics prediction and structural analysis of NnPAL1 indicate that NnPAL1 has similar structural features to the reported angiosperm PAL proteins.
Purification and functional characterisation of recombinant NnPAL1
Expression profile of NnPAL1 under stress conditions
Identification of the PAL family in N. nucifera from whole genomic sequences
The genomic DNA used for de novo sequencing is extracted from the clean shoots of N. nucifera. We use sixteen assembled virtual chromosomes of high quality as the resource for the PAL search. In previous reports of PAL from higher plants, all functionally identified PAL genes [13–23] encode approximately 700 amino acids and contain the characteristic conserved GTITASGDLVPLSYIA motif. Therefore, we set the sizes of the PAL family to larger than 500 amino acids with the GTITASGDLVPLSYIA signature. Three PAL genes, NnPAL1, NnPAL2 and NnPAL3, are located in the well-defined regions of the assembled sequences. Consistent with the phylogeny of angiosperms, NnPAL2 and NnPAL3 are similar to the PAL from dicots. However, NnPAL1 is similar to the PAL from gymnosperms. The full-length cDNA of NnPAL1 is cloned from the RNA transcripts of tender leaves using RACE method. NnPAL1 is transcribed with an intact open reading frame, suggesting that it does not become a pseudogene during evolution.
NnPAL1 from the genuine PAL family of N. nucifera
N. nucifera is a perennial aquatic plant. Therefore, obtaining pure tissues is a prerequisite for molecular biology experiments. An endophyte is a bacterial or fungal microorganism, which colonises inter- and/or intracellularly inside the healthy tissues of the host plant . We are careful to remove the residues from both shoots and leaves. To confirm that NnPAL1 is a member of the PAL family of N. nucifera and not endophytes, we performed several experiments.
First, we determine the location of NnPAL1 in Vchr3 and extract the upstream sequence (31,942 bps) and downstream sequence (26,288 bps) flanking NnPAL1 (Additional file 6). Then, we performed a discontiguous megablast search against the nucleotide collection database in NCBI to search for homologous regions. In the upstream 1–5000 bps of NnPAL1, we find out a highly homologous region to dicots. The sequences with the first three highest scores are uncharacterised mRNA from Vitis vinifera, the mRNA of a tetratricopeptide repeat-containing family protein in Populus trichocarpa, and an mRNA of a conserved hypothetical protein in Ricinus communis. In the downstream 10,000-15,000 bps of NnPAL1, we find out a homologous partial coding sequence of the GWD gene for alpha-glucan water dikinase from N. nucifera. A homologous region to dicots is also identified in the downstream 20,000-25,000 bps. The sequences with the first three highest scores are the mRNA of a zinc finger CCCH domain-containing protein 17-like in Citrus sinensis and the mRNA of a zinc finger family protein in Populus trichocarpa (Additional file 7). As a basal dicot, the flanking sequences of N. nucifera NnPAL1 show homology to sequences of all the other dicots.
Evolution of NnPAL1 in N. nucifera during the evolution of plants
In this study, three PALs, NnPAL1, NnPAL2 and NnPAL3, are identified using the database of whole genomic sequences as a resource. In previous reports, the angiosperm PAL had phase 2 introns at an Arg codon of [15, 31, 37, 38], but the gymnosperm PAL had no intron . Similarly, both NnPAL2 and NnPAL3 have only one intron of phase 2, whereas NnPAL1 has two introns of phase 0. This result demonstrates that NnPAL1 has unique gene structure that is different from NnPAL2, NnPAL3 and other PAL genes from angiosperms. This difference between NnPAL1 and other angiosperm PALs suggests that it is an ancient gene with a different evolutionary origin.
PAL and HAL are members of the lyase class I_like superfamily of enzymes, which catalyse similar beta-elimination reactions and are active as homotetramers. PAL and HAL diverged from each other when fungi and plants diverged from the other kingdoms (Figure 6). Because of their similar structures, PAL is derived from the His ammonia-lyase. HAL is a basic enzyme, participating in a central metabolic pathway, and PAL is derived from HAL to fulfil specific functions.
PAL is a ubiquitous higher-plant enzyme that catalyses the nonoxidative deamination of phenylalanine to trans-cinnamic acid. However, the origin and evolution of the PAL gene family in seed plants (Spermatophyta) have not been determined [31, 39]. Currently, two major mechanisms are responsible for the evolution and functional divergence of genes. One evolutionary mechanism is called HGT (horizontal gene transfer) and refers to the movement of genes between different species . HGT events occur only in plant mitochondrial genes [41–43], and rarely in nuclear genes . The other evolutionary mechanism is gene duplication, which is the main mechanism for evolutionary innovations and functional divergence . Based on morphological characteristics and molecular data, gymnosperms are considered ancestral to the angiosperms . At least three ancestral duplication events of PAL occurred, leading to three clades of gymnosperm PAL genes, gymnosperm-I, gymnosperm-II and gymnosperm-III. It appears that angiosperms diverged from gymnosperm III when only one paralogue PAL gene is retained within the angiosperms . In this study, we construct PAL phylogenetic trees that include the PAL gene families from Pinus taeda (gymnosperm I, II, III), monocots and dicots according to the gene sequences of the sequenced species (Figure 2 and Additional file 5: Figure S5). The phylogenetic trees show that NnPAL1 is clustered together with Pteda1143311 and Pteda17307 (gymnosperm II); however, NnPAL2 and NnPAL3 are clustered with dicots with high bootstrap and posterior probability values. Perhaps, NnPAL1 has a different evolutionary origin from NnPAL2 and NnPAL3. Except for NnPAL1, the other PAL clusters are monophyletic after the split between dicots and monocots (Figure 2 and Additional file 5: Figure S5). However, the PAL from one species is clustered together with the other species rather than with a single species. This result indicates that duplication events are important in the evolution of PAL genes after the split between dicots and monocots [46, 47].
During evolution, NnPAL1 is found to be an ancient member of the PAL family that has been retained in angiosperms. A different evolutionary history for PAL genes in angiosperms suggests different mechanisms of functional regulation. In the phylogenetic trees of PAL (Figure 2), NnPAL1 is not found where expected. Interestingly, NnPAL1 shows high homology to Pteda1143311 and Pteda17307 from Pinus taeda, and Pinus taeda is also rich in various secondary metabolites. There may be a shared secondary metabolite produced by NnPAL1 or Pteda1143311 and Pteda17307. Moreover, this specific product may protect N. nucifera and Pinus taeda from similarly extreme environments. NnPAL1 may have been essential for N .nucifera to survive in harsh environments during the Cretaceous period.
We speculate that the angiosperm PAL is not of monophyletic origin. Ancestral gene duplication and vertical inheritance from gymnosperms may occur during evolution from parent to offspring. In gymnosperms another paralogue of the ancient PAL exists that is retained prior to the formation of angiosperms. NnPAL1 may be derived from the product of gymnosperm-II PAL. Discovering a functional NnPAL1 indicates that angiosperm PAL genes are not derived from a single gene in the ancestral angiosperm genome. However, similar modification sites and structure to other angiosperm PALs suggest that NnPAL1 can catalyse the deamination of phenylalanine to trans-cinnamic acid and is involved in the phenylpropanoid pathway.
Functional characterisation and expression patterns of NnPAL1
During evolution, NnPAL1 remained functional with both PAL and TAL activities. Compared with other PALs cloned from other plants [14, 18], both the PAL and TAL activities of NnPAL1 show higher Km values, which can be explained as follows: (I) during the evolution of angiosperms, the function of the most ancient PAL (NnPAL1) is gradually replaced by a new PAL; (II) NnPAL1 has many posttranslational modification sites, which may be involved in the subunit turnover of NnPAL1 in vivo, and prokaryote expression systems lack multiple protein modifications, which affect enzyme protein stability; and (III) the ancient NnPAL1 has evolved a novel function required for other metabolic pathways . The optimum pH is 9.0 and the optimum temperature is 55°C, which is similar to other PALs from higher plants . The expression patterns are validated by real-time PCR. In response to environmental stress during the Cretaceous period, N. nucifera is eventually trapped in aquatic areas of Asia . This type of environment makes secondary metabolites important for N. nucifera because they protect the species from various stimuli. PAL expression is regulated by various factors at the transcriptional level. Bioinformatics analysis of the upstream cis-elements in NnPAL1 identified several related regulatory elements, such as the MYB binding site involved in drought-inducibility, the auxin-responsive element (IAA), the fungal elicitor responsive element, the cis-acting elements involved in abscisic acid responsiveness (ABA), and the light responsive element in NnPAL1 gene. All the treatments used in this study cause increases of the NnPAL1 transcripts, which suggests that NnPAL1 is regulated by these elements. The ancient NnPAL1 of N. nucifera is involved in the response to stressful environments, which makes N. nucifera the representative of plants that survived from the Cretaceous period .
Using comparative genomics and phylogenetic analyses, three PAL members, NnPAL1, NnPAL2 and NnPAL3, are identified. The distinction between NnPAL1 and other angiosperm PALs suggests that NnPAL1 is not derived from a PAL paralogue of a gymnosperm leading to angiosperms. We postulate that there may be another ancestral duplication event and vertical inheritance from the gymnosperms. The ancient PAL NnPAL1 from N. nucifera is characterised at both the RNA and protein levels in vitro. The unique biochemical characteristics of N. nucifera may allow it to overcome the harsh environment. Additionally, as a basal dicot, N. nucifera is a perennial aquatic plant with agricultural, evolutionary and medicinal importance [26, 27]. Polyphenolic compounds in N. nucifera have important pharmacological and physiological activities. The discovery and characterisation of an ancient NnPAL1 provides new insight into PAL evolution in angiosperms and may also lead to improved function through the genetic engineering of N. nucifera.
Identification of the PAL gene family in N. nucifera
High purity DNA is extracted from clean and tender shoots of N. nucifera, and is used for de novo sequencing. For the de novo assembly, 16.4 Gb of filtered data with 15-fold depth is used. Sixteen virtual chromosomes (2n = 16) are assembled with high quality.
A local database with the sequences of sixteen virtual chromosomes are constructed on the Bio-Linux platform;
tBLASTN is conducted with a cut-off E value of 1e-20 in the local database with each member of the AtPAL family;
the aligned frames containing a highly conserved phenylalanine and histidine ammonia-lyase signature (GTITASGDLVPLSYIA) are selected for further analysis;
the related genome sequences with intact open reading frames are located in the well-defined region of assembled sequences and are extracted;
the extracted codes larger than 500 amino acids are selected and annotated.
Identification of gene families in other plants and construction of the PAL phylogenetic tree
Based on a conserved phenylalanine and histidine ammonia-lyase signature, the PAL families from other plants except for Pinus taeda are identified and downloaded from Phytozome (http://www.phytozome.net) with a cut-off E value of 1e−20. The PAL proteins from Pinus taeda are deduced from their RNA transcripts . The analysed species (Table 1) are as follows: one Bryophyta (Physcomitrella patens), one gymnosperm (Pinus taeda), five monocotyledons (Brachypodium distachyon, Oryza sativa, Sorghum bicolor, Setaria italica, Zea mays), and eight dicots (Aquilegia caerulea, Arabidopsis thaliana, Cucumis sativus, Glycine max, Mimulus guttatus, Populus trichocarpa, Vitis vinifera, and Nelumbo nucifera). The protein sequences are aligned with the CLUSTALW program  with manual adjustments. The phylogenetic trees are simultaneously inferred from the protein alignment using three methods as follows: the NJ (Neighbour-joining) tree with the JTT model, the ML (maximum likelihood) tree with LG model, and BI (Bayesian inference) tree with the GTR model, are generated with Mega 5 , PhyML 3.0  and Mrbayes 3.2 , respectively. The bootstrap values are set 1000 for the neighbour-joining and maximum likelihood tree. For Bayesian inference, we sample every 10 generations for 300,000 total generations on two independent parallel runs of the Monte Carlo Markov Chains. Then, the average standard deviation of the split frequencies is calculated to check the convergence of the two runs.
Plant material, cloning and expression vectors
N. nucifera mature seeds are harvested from East Lake of Wuhan, China. The tender leaves are collected when the seedling germinated from seeds in the greenhouse. E.coli top10 (TaKaRa, Dalian, China) is used as the host for plasmid pMD18-T vector (TaKaRa, Dalian, China) amplification. E.coli BL21(DE3) is selected as the host for pET-28a(+) expression vector.
Isolation of the full-length NnPAL1 cDNA
The total RNA is isolated from the leaves using a modified CTAB method . The first strand cDNA is produced by RT-PCR using reverse transcriptase (MBI Fermentas). Two degenerate primers, (Nf-F) 5′-GCNTCNGGHGAYYTDGTBCC-3′ and (Nf-R) 5′-ARNCCBARDGART TWACATC-3′, are designed according to the highly conserved regions of the plant PAL for both the amino acid and nucleotide sequences. The partial cDNA of PAL is amplified using the following conditions: initial denaturing at 94°C for 4 min, followed by 35 cycles of 94°C for 40 s, 59°C for 40 s, and 72°C for 80 s, with a final extension at 72°C for 10 min. The target fragment is checked on 1% agarose gel and purified with a Gel Extraction Kit (BioDev-Tech, Beijing, China). The purified products are then ligated into the pMD18-T Easy vector (TaKaRa, Dalian, China), transformed into competent E. coli Top10 and sequenced on an ABI 3730.
The full-length cDNA of the PAL gene is isolated by 3′- and 5′-RACE using the RACE Kit (TaKaRa, Dalian, China). Based on the sequenced DNA fragment, four gene-specific primers, 3′GSP1 (5′-CTGGACTACGGATTCAAGGGTG-3′), 3′GSP2 (5′-TCAGTATTTGGCAAACC CAGTCA-3′), 5′GSP1(5′-AGCATCACTTCGCAGAACATCG-3′) and 5′GSP2 (5′-GTACGGAC CTTGGAGTTGGGAC-3′), are designed for the 3′- and 5′-RACE experiments, respectively. A 860-bp fragment and 791-bp fragment are then obtained by 3′-RACE and 5′-RACE, respectively. The full-length coding cDNA of the 2154-bps is amplified and sequenced using two gene-specific primers, 5′-GAATTCATGGTTGCAGGGGCCGAGATAG-3′ and 5′-CCCTCGAGCACAAGAAGGCAACACCAAAGT-3′.
Bioinformatics analysis of NnPAL1
The amino acid sequence and protein analysis of NnPAL1 are performed with the ExPASy tools (http://us.expasy.org/tools) and NCBI server (http://www.ncbi.nlm.nih.gov/). The possible posttranslational modification sites are predicted by PROSITE (http://prosite.expasy.org/). The prediction of secondary structure and trans-membrane helices in the PAL protein are performed with SOPMA (http://pbil.ibcp.fr/htm/index.php) and the TMHMM Server v. 2.0 (http://www.cbs.dtu.dk/services/TMHMM-2.0/), respectively. Homology modelling is performed with Swiss-Model (http://swissmodel.expasy.org/) and is based on the PAL crystal structure from Petroselinum crispum.
Expression of NnPAL1 in E. coli
Primers NnP1 (5′-GAATTCATGGTTGCAGGGGCCGAGATAG-3′, the italics is EcoRI restriction site) and NnP2 (5′-CCCTCGAGCACAAGAAGGCAACACCAAAGT-3′, the italics is XhoI restriction site) are used to amplify the NnPAL1 gene coding sequence. The PCR products are digested with EcoRI and XhoI and then inserted into pET28a(+) expression vector. The recombinant plasmid NnPAL1-pET28a(+) is transformed into the BL21 strain and sequenced to confirm the correct ORF of NnPAL1.
The transformant with the correct NnPAL1-pET28a(+) is selected and cultured in Luria–Bertani (LB) medium containing 50 μg/ml kanamycin at 37°C until the OD600 reached 0.6. Protein expression is induced with 0.5 mM isopropyl β-D-1-thiogalactopyranoside (IPTG) at 16°C for 12 h. The recombinant proteins are purified on a Ni-NTA agarose column and eluted with a step gradient of imidazole buffers (10 mM, 50 mM, 100 mM, and 200 mM). The purity of the recombinant protein is verified by SDS-PAGE. The 200 mM fractions are dialysed with Spectra/Por Membranes (MWCO: 8,000-14,000) in dialysis buffer.
Enzyme activity assay for the recombinant NnPAL1 protein
The protein concentrations are determined with the G250 dye-binding method  using bovine serum albumin as the protein standard. The enzyme activity of the recombinant NnPAL1 is assayed by measuring the trans-cinnamic acid formation at 290 nm  and p-coumaric acid formation at 310 nm . The PAL activity and TAL activity is expressed in nkat (nanomole of trans-cinnamic acid/p-coumaric acid formed per second).
To determine the optimum temperature and optimum pH for enzyme activity, several assays are performed at pH 8.5 for 30 min at varying temperatures (4, 23, 30, 37, 45, 50, 55, 60, 70, 80 and 90°C), and at 37°C for 30 min with buffer of various of pH (5, 6, 7, 7.5, 8, 8.5, 9, 10, 11), respectively. The reactions are performed in 150 μl reaction mixtures with 6 μg recombinant NnPAL1, 15 mM L-phenylalanine and 50 mM Tris–HCl (pH 8.5), and are terminated with the addition of concentrated HCl .
To determine the kinetic parameters and substrate specificity, 150 μl reaction mixtures containing 6 μg recombinant NnPAL1 proteins, 50 mM Tris–HCl (pH 8.5) and a range of L-phenylalanine (0.15-15 mM) or L-tyrosine (0.3-2 mM) concentrations are used. Hyperbolic plots and double reciprocal plots (Lineweaver–Burk plot) are used to calculate the Km (Michaelis-Menten constant) using the Michaelis-Menten equation .
Cis-regulatory element analysis and expression of NnPAL1 by quantitative real-time PCR
The 5′ upstream region of NnPAL1 is characterised using BLASTN against the whole genomic sequence of N. nucifera with NnPAL1 gene. We predicted the cis-elements by submitting 5′ fragment to PlantCARE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/).
Two-week-old leaves are treated with 250 μM abscisic acid (ABA), 100 ng/ml IAA, ultraviolet light, Neurospora crassa (fungi) and drought according to the cis-elements. The treated leaves are harvested and immediately frozen after 0, 4, 8, and 12 h. The total RNA is isolated from the leaves and treated with RNase-free DNase I to avoid DNA contamination. Real time RT-PCR analysis of NnPAL1 is performed on an Applied Biosystems StepOne Plus, using β-actin gene (β-actin-F: 5′-CCTGATGGGCAAGTGATT-3′, β-actin-R: 5′-GCTCATACGGTCAG CAATA-3′) as an internal control for all the samples.
One ancient member from PAL gene family of Nelumbo nucifera
Degenerate oligonucleotide primer PCR
Rapid amplification of cDNA ends
The authors thank Dr Huabin Zhao in Wuhan University for helps in evolutionary analyses. This research is financially supported by National Natural Science Foundation of China (31271310).
- Dixon RA, Paiva NL: Stress-induced phenylpropanoid metabolism. Plant Cell. 1995, 7: 1085-1097. 10.1105/tpc.7.7.1085.PubMedPubMed CentralView ArticleGoogle Scholar
- Pellegrini L, Rohfritsch O, Fritig B, Legrand M: Phenylalanine ammonia-lyase in tobacco: molecular cloning and gene expression during the hypersensitive reaction to tobacco mosaic virus and the response to a fungal elicitor. Plant Physiol. 1994, 106: 877-886. 10.1104/pp.106.3.877.PubMedPubMed CentralView ArticleGoogle Scholar
- Hamberger B, Ellis M, Friedmann M, de Azevedo SC, Barbazuk B, Douglas CJ: Genome-wide analyses of phenylpropanoid-related genes in Populus trichocarpa, Arabidopsis thaliana, and Oryza sativa: the Populus lignin toolbox and conservation and diversification of angiosperm gene families. Can J Bot. 2007, 85 (12): 1182-1201. 10.1139/B07-098.View ArticleGoogle Scholar
- Naoumkina MA, Zhao Q, Gallego-Giraldo L, Dai X, Zhao PX, Dixon RA: Genome-wide analysis of phenylpropanoid defence pathways. Mol Plant Pathol. 2010, 11 (6): 829-846.PubMedGoogle Scholar
- Raes J, Rohde A, Christensen JH, Van de Peer Y, Boerjan W: Genome-wide characterization of the lignification toolbox in Arabidopsis. Plant Physiol. 2003, 133 (3): 1051-1071. 10.1104/pp.103.026484.PubMedPubMed CentralView ArticleGoogle Scholar
- Tsai CJ, Harding SA, Tschaplinski TJ, Lindroth RL, Yuan Y: Genome-wide analysis of the structural genes regulating defense phenylpropanoid metabolism in Populus. New Phytol. 2006, 172 (1): 47-62. 10.1111/j.1469-8137.2006.01798.x.PubMedView ArticleGoogle Scholar
- MacDonald MJ, D’Cunha GB: A modern view of phenylalanine ammonia lyase. Biochem Cell Biol. 2007, 85: 273-282. 10.1139/O07-018.PubMedView ArticleGoogle Scholar
- Kim W, Erlandsen H, Surendran S, Stevens RC, Gamez A, Michols-Matalon K, Tyring SK, Matalon R: Trends in enzyme therapy forphenylketonuria. Mol Ther. 2004, 10: 220-224. 10.1016/j.ymthe.2004.05.001.PubMedView ArticleGoogle Scholar
- Kyndt JA, Meyer TE, Cusanovich MA, Van Beeumen JJ: Characterization of a bacterial tyrosine ammonia lyase, a biosynthetic enzyme for the photoactive yellow protein. FEBS Lett. 2002, 512: 240-244. 10.1016/S0014-5793(02)02272-X.PubMedView ArticleGoogle Scholar
- Moffitt MC, Louie GV, Bowman ME, Pence J, Noel JP, Moore BS: Discovery of two cyanobacterial phenylalanine ammonia lyases: kinetic and structural characterization. Biochemistry. 2007, 46: 1004-1012. 10.1021/bi061774g.PubMedPubMed CentralView ArticleGoogle Scholar
- Ogata K, Uchiyama K, Yamada H: Metabolism of aromatic amino acid in microorganisms. Part I: formation of cinnamic acid from phenylalanine. Agric Biol Chem. 1967, 31: 200-206. 10.1271/bbb1961.31.200.View ArticleGoogle Scholar
- Xiang L, Moore BS: Biochemical characterization of a prokaryotic phenylalanine ammonia lyase. J Bacteriol. 2005, 187: 4286-4289. 10.1128/JB.187.12.4286-4289.2005.PubMedPubMed CentralView ArticleGoogle Scholar
- Xu F, Cai R, Cheng S, Du H, Wang Y, Cheng S: Molecular cloning, characterization and expression of phenylalanine ammonia-lyase gene from Ginkgo biloba. Afr J Biotechnol. 2008, 7: 721-729.Google Scholar
- Okada T, Mikage M, Sekita S: Molecular characterization of the phenylalanine Ammonia-Lyase from Ephedra sinica. Biol Pharm Bull. 2008, 31: 2194-2199. 10.1248/bpb.31.2194.PubMedView ArticleGoogle Scholar
- Minami E, Ozeki Y, Matsuoka M, Koizuka N, Tanaka Y: Structure and some characterization of the gene for phenylalanine ammonialyase from rice plants. Eur J Biochem. 1989, 185: 19-25. 10.1111/j.1432-1033.1989.tb15075.x.PubMedView ArticleGoogle Scholar
- Lu BB, Du Z, Ding RX, Zhang L, Yu XJ, Liu CH, Chen WS: Cloning and characterization of a differentially expressed phenylalanine ammonia-lyase gene (IiPAL) after genome duplication from tetraploid Isatis indigotica Fort. J Integr Plant Biol. 2006, 48 (12): 1439-1449. 10.1111/j.1744-7909.2006.00363.x.View ArticleGoogle Scholar
- Wanner LA, Li G, Ware D, Somssich IE, Davis KR: The phenylalanine ammonia lyase gene family in Arabidopsis thaliana. Plant Mol Biol. 1995, 27: 327-338. 10.1007/BF00020187.PubMedView ArticleGoogle Scholar
- Gao JH, Zhang SW, Cai F, Zheng XJ, Lin N, Qin XB, Ou YC, Gu XP, Zhu XH, Xu Y, Chen F: Characterization, and expression profile of a phenylalanine ammonia lyase gene from Jatropha curcas L. Mol Biol Rep. 2012, 39: 3443-3452. 10.1007/s11033-011-1116-4.PubMedView ArticleGoogle Scholar
- Jiang YM, Xia N, Li XD, Shen WB, Liang LJ, Wang CY, Wang R, Peng F, Xia B: Molecular cloning and characterization of a phenylalanine ammonia-lyase gene (LrPAL) from Lycoris radiate. Mol Biol Rep. 2011, 38: 1935-1940. 10.1007/s11033-010-0314-9.PubMedView ArticleGoogle Scholar
- Ritter H, Schulz GE: Structural basis for the entrance into the phenylpropanoid metabolism catalyzed by phenylalanine ammonia-lyase. Plant Cell. 2004, 16: 3426-3436. 10.1105/tpc.104.025288.PubMedPubMed CentralView ArticleGoogle Scholar
- Joos HJ, Hahlbrock K: Phenylalanine ammonia-lyase in potato (Solanum tuberosumL.). Genomic complexity, structural comparison of two selected genes and modes of expression. Eur J Biochem. 1992, 204: 621-629. 10.1111/j.1432-1033.1992.tb16675.x.PubMedView ArticleGoogle Scholar
- Chang A, Lim MH, Lee SW, Robb EJ, Nazar RN: Tomato phenylalanine ammonia-lyase gene family, highly redundant but strongly underutilized. J Biol Chem. 2008, 283: 33591-33601. 10.1074/jbc.M804428200.PubMedPubMed CentralView ArticleGoogle Scholar
- Huang JL, Gu M, Lai ZB, Fan BF, Shi K, Zhou YH, Yu JQ, Chen ZX: Functional analysis of the arabidopsis PAL Gene family in plant growth, development, and response to environmental stress. Plant Physiol. 2010, 153: 1526-1538. 10.1104/pp.110.157370.PubMedPubMed CentralView ArticleGoogle Scholar
- Schwede TF, Re’tey J, Schulz GE: Crystal structure of histidine ammonia-lyase revealing a novel polypeptide modification as the catalytic electrophile. Biochemistry. 1999, 38: 5355-5361. 10.1021/bi982929q.PubMedView ArticleGoogle Scholar
- Watanabe SK, Hernandez-Velazco G, Iturbe-Chinas F, Lopez-Munguia A: Phenylalanine ammonia lyase from Sporidiobolus pararoseus and Rhodosporidium toruloides: application for phenylalanine and tyrosine deamination. World J Microbiol Biotechnol. 1992, 8: 406-410. 10.1007/BF01198755.View ArticleGoogle Scholar
- Mukherjee PK, Mukherjee D, Maji AK, Rai S, Heinrich M: The sacred lotus (Nelumbo nucifera)-phytochemical and therapeutic profile. J Pharm Pharmacol. 2009, 61 (4): 407-422.PubMedGoogle Scholar
- Kashiwada Y, Aoshima A, Ikeshiro Y, Chen YP, Furukawa H, Itoigawa M, Fujioka T, Mihashi K, Cosentino LM, Morris-Natschke SL, Lee KH: Anti-HIV benzylisoquinoline alkaloids and flavonoids from the leaves of Nelumbo nucifera, and structure-activity correlations with related alkaloids. Bioorg Med Chem. 2005, 13: 443-448. 10.1016/j.bmc.2004.10.020.PubMedView ArticleGoogle Scholar
- Hsu J: Late cretaceous and cenozoic vegetation in China, emphasizing their connections with north America. Ann Mo Bot Gard. 1983, 70: 490-508. 10.2307/2992084.View ArticleGoogle Scholar
- Ming R, Vanburen R, Liu Y, Yang M, Han Y, Li LT, Zhang Q, Kim MJ, Schatz MC, Campbell M, Li J, Bowers JE, Tang H, Lyons E, Ferguson AA, Narzisi G, Nelson DR, Blaby-Haas CE, Gschwend AR, Jiao Y, Der JP, Zeng F, Han J, Min XJ, Hudson KA, Singh R, Grennan AK, Karpowicz SJ, Watling JR, Ito K, et al: Genome of the long-living sacred lotus (Nelumbo nucifera Gaertn.). Genome Biol. 2013, 14: R41-10.1186/gb-2013-14-5-r41.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang Y, Fan G, Liu Y, Sun F, Shi C, Liu X, Peng J, Chen W, Huang X, Cheng S, Liu Y, Liang X, Zhu H, Bian C, Zhong L, Lv T, Dong H, Liu W, Zhong X, Chen J, Quan Z, Wang Z, Tan B, Lin C, Mu F, Xu X, Ding Y, Guo AY, Wang J, Ke W: The sacred lotus genome provides insights into the evolution of flowering plants. Plant J. 2013, 76: 557-567. 10.1111/tpj.12313.PubMedView ArticleGoogle Scholar
- Schmidt K, Heberle B, Kurrasch J, Nehls R, Stahl DJ: Suppression of phenylalanine ammonia lyase expression in sugar beet by the fungal pathogen Cercospora beticola is mediated at the core promoter of the gene. Plant Mol Biol. 2004, 55: 835-852. 10.1007/s11103-005-2141-2.PubMedView ArticleGoogle Scholar
- Angiosperm Phylogeny Group: An update of the Angiosperm Phylogeny Group classification for the orders and families of flowering plants: APG III. Bot J Linn Soc. 2003, 141: 399-436.View ArticleGoogle Scholar
- Rösler J, Krefel F, Amrhein N, Sohmid J: Maize phenylalanine ammonia-lyase activity. Plant Physiol. 1997, 113: 175-179. 10.1104/pp.113.1.175.PubMedPubMed CentralView ArticleGoogle Scholar
- Bagal UR, Leebens-Mack JH, Lorenz WW, Dean JFD: The phenylalanine ammonia lyase (PAL) gene family shows a gymnosperm-specific lineage. BMC Genomics. 2012, 13 (Suppl.3): S1-PubMedPubMed CentralGoogle Scholar
- Hsieh LS, Hsieh YL, Yeh CS, Cheng CY, Yang CC, Lee PD: Molecular characterization of a phenylalanine ammonia-lyase gene (BoPAL1) from B.oldhamii. Mol Biol Rep. 2011, 38: 283-290. 10.1007/s11033-010-0106-2.PubMedView ArticleGoogle Scholar
- Tan RX, Zou WX: Endophytes: a rich source of functional metabolites. Nat Prod Rep. 2001, 18: 448-459. 10.1039/b100918o.PubMedView ArticleGoogle Scholar
- Fiona CC, Laurence BD, Norman GL: The Arabidopsis phenylalanine ammonia lyase gene family: kinetic characterization of the four PAL isoforms. Phytochemistry. 2004, 65: 1557-1564. 10.1016/j.phytochem.2004.05.006.View ArticleGoogle Scholar
- Lee SW, Robb J, Nazar RN: Truncated phenylalanine ammonia-lyase expression in tomato (Lycopersicon esculentum). J Biol Chem. 1992, 267: 11824-11830.PubMedGoogle Scholar
- Chaw SM, Zhaekikh A, Sung HM, Lau TC, Li WH: Molecular phylogeny of extant gymnosperms and seed plant evolution: analysis of nuclear 18S rRNA sequences. Mol Biol Evol. 1997, 14: 56-68. 10.1093/oxfordjournals.molbev.a025702.PubMedView ArticleGoogle Scholar
- Keeling PJ, Palmer JD: Horizontal gene transfer in eukaryotic evolution. Nat Rev Genet. 2008, 9: 605-618. 10.1038/nrg2386.PubMedView ArticleGoogle Scholar
- Bergthorsson U, Adams KL, Thomason B, Palmer JD: Widespread horizontal transfer of mitochondrial genes in flowering plants. Nature. 2003, 424: 197-201. 10.1038/nature01743.PubMedView ArticleGoogle Scholar
- Bergthorsson U, Richardson AO, Young GJ, Goertzen LR, Palmer JD: Massive horizontal transfer of mitochondrial genes from diverse land plant donors to the basal angiosperm Amborella. Proc Natl Acad Sci USA. 2004, 101: 17747-17752. 10.1073/pnas.0408336102.PubMedPubMed CentralView ArticleGoogle Scholar
- Hao W, Richardson AO, Zheng Y, Palmer JD: Gorgeous mosaic of mitochondrial genes created by horizontal transfer and gene conversion. Proc Natl Acad Sci USA. 2010, 107: 21576-21581. 10.1073/pnas.1016295107.PubMedPubMed CentralView ArticleGoogle Scholar
- Rumpho ME, Worful JM, Lee J, Kannan K, Tyler MS, Bhattacharya D, Moustafa M, Manhart JR: Horizontal gene transfer of the algal nuclear gene psbO to the photosynthetic sea slug Elysia chlorotica. Proc Natl Acad Sci USA. 2008, 105: 17867-17871. 10.1073/pnas.0804968105.PubMedPubMed CentralView ArticleGoogle Scholar
- Yang J, Huang JX, Gu HY, Zhong Y, Yang ZH: Duplication and adaptive evolution of the chalcone synthase genes of Dendranthema (Asteraceae). Mol Biol Evol. 2002, 19: 1752-1759. 10.1093/oxfordjournals.molbev.a003997.PubMedView ArticleGoogle Scholar
- Kumar A, Ellis BE: The phenylalanine ammonia-lyase gene family in raspberry. Structure, expression, and evolution. Plant Physiol. 2001, 127: 230-239. 10.1104/pp.127.1.230.PubMedPubMed CentralView ArticleGoogle Scholar
- Röther D, Poppe L, Morlock G, Viergutz S, Rétey J: An active site homology model of phenylalanine ammonia-lyase from Petroselinum crispum. Eur J Biochem. 2002, 269: 3065-3075. 10.1046/j.1432-1033.2002.02984.x.PubMedView ArticleGoogle Scholar
- Allwood EG, Davies DR, Gerrish C, Ellis BE, Bolwell GP: Phosphorylation of phenylalanine ammonia-lyase: evidence for a novel protein kinase and identification of the phosphorylated residue. FEBS Lett. 1999, 457: 47-52. 10.1016/S0014-5793(99)00998-9.PubMedView ArticleGoogle Scholar
- Xu EY, Moore FL, Reijo Pera RA: A gene family required for human germ cell development evolved from an ancient meiotic gene conserved in all metazoans. Proc Natl Acad Sci USA. 2001, 98: 7414-7419. 10.1073/pnas.131090498.PubMedPubMed CentralView ArticleGoogle Scholar
- Li JK, Zhou EX, Li DX, Huang SQ: Multiple northern refugia for Asian sacred lotus, an aquatic plant with characteristics of ice-age endurance. Aust J Bot. 2010, 58: 463-472. 10.1071/BT10002.View ArticleGoogle Scholar
- Thompson JD, Gibson TJ, Higgins DG: Multiple sequence alignment using ClustalW and ClustalX. Curr Protoc Bioinformatics. 2002, 00: 2.3.1-2.3.22.Google Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using likelihood, distance, and parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.PubMedPubMed CentralView ArticleGoogle Scholar
- Guindon S, Delsuc F, Dufayard JF, Gascuel O: Estimating maximum likelihood phylogenies with PhyML. Methods Mol Biol. 2009, 537: 113-137. 10.1007/978-1-59745-251-9_6.PubMedView ArticleGoogle Scholar
- Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.PubMedView ArticleGoogle Scholar
- Gambino G, Perrone I, Gribaudo I: A Rapid and effective method for RNA extraction from different tissues of grapevine and other woody plants. Phytochem Anal. 2008, 19: 520-525. 10.1002/pca.1078.PubMedView ArticleGoogle Scholar
- Bradford MM: A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. Anal Biochem. 1976, 72: 248-254. 10.1016/0003-2697(76)90527-3.PubMedView ArticleGoogle Scholar
- D’Cunha GB, Satyanarayan S, Nair PM: Purification of phenylalanine ammonia lyase from Rhodotorula glutinis. Phytochemistry. 1996, 42: 17-20. 10.1016/0031-9422(95)00914-0.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.