- Research article
- Open Access
Genomic survey, characterization and expression profile analysis of the peptide transporter family in rice (Oryza sativa L.)
https://doi.org/10.1186/1471-2229-10-92
© Zhao et al; licensee BioMed Central Ltd. 2010
- Received: 22 January 2010
- Accepted: 20 May 2010
- Published: 20 May 2010
Abstract
Background
Peptide transporter (PTR) family whose member can transport di-/tripeptides and nitrate is important for plant growth and development. Although the rice (Oryza sativa L.) genome has been sequenced for a few years, a genomic survey, characterization and expression profile analysis of the PTR family in this species has not been reported.
Results
In this study, we report a comprehensive identification, characterization, phylogenetic and evolutionary analysis of 84 PTR family members in rice (OsPTR) as well as their whole-life expression patterns. Chromosomal distribution and sequence analysis indicate that nearly 70% of OsPTR members are involved in the tandem and segmental duplication events. It suggests that genome duplication might be a major mechanism for expansion of this family. Highly conserved motifs were identified in most of the OsPTR members. Meanwhile, expression profile of OsPTR genes has been analyzed by using Affymetrix rice microarray and real-time PCR in two elite hybrid rice parents, Minghui 63 and Zhenshan 97. Seven genes are found to exhibit either preferential or tissue-specific expression during different development stages of rice. Under phytohormone (NAA, GA3 and KT) and light/dark treatments, 14 and 17 OsPTR genes are differentially expressed respectively. Ka/Ks analysis of the paralogous OsPTR genes indicates that purifying selection plays an important role in function maintenance of this family.
Conclusion
These investigations add to our understanding of the importance of OsPTR family members and provide useful reference for selecting candidate genes for functional validation studies of this family in rice.
Keywords
- Tandem Duplication
- Segmental Duplication
- Peptide Transporter
- Nitrate Transporter
- Phytohormone Treatment
Background
Nutrient transport is essential to life and occurs in both prokaryotes and eukaryotes. It is well established that peptide transporters play an important role in the nutrition of bacteria, yeasts and animals [1]. However, the role of transporters of small peptides (2-6 amino acids) in plants is less well defined [2]. Peptide transporter (PTR) family (TC 2.A.17), also called the proton-dependent oligopeptide transport (POT) family [3, 4] is a family consists of di-/tripeptides transporters [5]. Besides, some plant PTRs acted as nitrate transporters are termed as the members of Nitrate Transporter 1 family (NRT1, [6]). PTR proteins in plants show significant sequence homology and contain several transmembrane (TM) regions, with a large hydrophilic loop among TMs. Generally, the identified PTR members in plants fall into three types based on the nature of their substrates: di-/tripeptides transporter, nitrate transporter and other substrates transporter [6].
AtPTR2 (At2g02040) was the first identified di-/tripeptides transporter in Arabidopsis [7–10] and had high mRNA expression levels in 3-d-germinating seed, root and young leaf [8]. The antisense plants of AtPTR2 exhibited delayed flowering time and arrested seed development [9, 11]. AtPTR1 (At3g54140) transported di-/tripeptides with low selectivity as well as substrates lacking a peptide bond. AtPTR1 was expressed in the vascular tissue throughout the plant, indicative of a role in long-distance transport [10]. AtPTR3 (At5g46050) was a salt stress and mechanical wounding inducible gene, JA and SA were both involved in regulation of it [12, 13]. AtPTR5 (At5g01180) which mediated high-affinity transport of dipeptides was most likely supplying peptides to maturating pollen, developing ovules and seeds. Overexpression of it resulted in enhanced shoot growth and increased N content [14]. The barley peptide transporter gene, HvPTR1, had seed-specific expression. Transport activity of HvPTR1 was regulated by phosphorylation in response to rising levels of amino acids in the germinating grain [15–17]. Functional di-/tripeptides transporters also have been reported in Nepenthes (NaNTR1, [18]), faba bean (VfPTR1, [19]) and Hakea actites (HaPTR4, [20]).
Although PTR is the abbreviation of peptide transporter, many cloned PTR members in plants are nitrate transporters, especially in Arabidopsis. AtNRT1.1 (At1g12110, CHL1), a dual-affinity nitrate transporter, was regulated by phosphorylation and involved in nascent organ development [21–26]. AtNRT1.2 (At1g69850, NTL1) encoded a low-affinity nitrate transporter. RNAi plants of it exhibited reduced nitrate-induced membrane depolarization and nitrate uptake activities [27]. The characterization of AtNRT1.4 (At2g26690) revealed the special role of petiole in nitrate homeostasis [28]. AtNRT1.5 (At1g32450) and AtNRT1.6 (At1g27080) were both plasma membrane proteins. AtNRT1.5 was expressed in root pericycle cells close to the xylem and participated in root xylem loading of nitrate [29] while AtNRT1.6 was only expressed in reproductive tissue and involved in delivering nitrate from maternal tissue to the developing embryo [30]. CsNitr1, a nitrite transporter in Cucumis sativus, had two isoforms: CsNitr1-S and CsNitr1-L. In contrast to that CsNitr1-S enhanced the efflux of nitrite from the cell, CsNitr1-L might load cytosolic nitrite into chloroplast stroma [31].
Amino acids and other substrates can also be transported by some PTRs. BnNRT1.2, isolated from Brassica napus could transport both nitrate and histidine [32]. AgDCAT1 localized at the symbiotic interface was a dicarboxylates transporter in Alnus glutinosa and the mRNA of AgDCAT1 was only detected in the nodules [33].
In rice, only two PTR genes have been functionally verified. OsNRT1 (LOC_Os03g13274), which displayed low-affinity nitrate transport activity was constitutively expressed in the most external layer of the root, epidermis and root hair [34]. SP1 (LOC_Os11g12740) which determined the panicle size had high expression level in the phloem of the branches of young panicles. Phylogenetic analysis implied that SP1 might be a nitrate transporter, however, neither nitrate nor other compounds transport activity could be obtained from it [35].
As the main staple food for a large segment of the world population, rice also serves as a model plant for monocotyledon species [36, 37]. Nutrient transport is critical for plant growth and development, in which PTR family members may play important roles. This work focused on a comprehensive identification, characterization, phylogenetic and evolutionary analysis of all the PTR family members in rice (OsPTR) as well as their whole-life expression profiling. Moreover, it has been demonstrated that the expression of some OsPTR genes is regulated phytohormone (NAA, KT and GA3) and light/dark treatments.
Results
Collection and identification of PTR genes in rice
Domain search (PF00854) of MSU Rice Genome Annotation Project (MSU RGAP) database released 125 sequences, out of which 91 were unique. Keywords "PTR", "peptide transporter" and "proton-dependent oligopeptide transport" searching identified 58, 85 and 10 genes respectively. BLASTP and TBLASTN searches against rice genome in three databases: MSU RGAP, NCBI and KOME released 95, 85 and 116 sequences respectively. By removal of the same sequences and different transcripts of the same gene, we initially identified 93 putative PTR genes in rice. Protein sequences of these genes were subjected to Pfam and SMART searches for the presence of PTR domains. As a result, seven genes whose protein contained only partial PTR domain and two genes annotated as retrotransposons were excluded from further analysis. Taken together, a total of 84 genes were predicted to encode PTR proteins in rice. Eighty-three genes had corresponding locus IDs in MSU RGAP database and one gene (Os01g0960900) was only found in NCBI. For convenience, all the "LOC_" prefix of MSU RGAP locus IDs were omitted representing OsPTR genes or proteins for further analysis. Meanwhile, similar methods were used and 53 PTR family members in Arabidopsis (AtPTR) were identified. The detailed information about the full-length cDNA, transcript, chromosomal position, BAC accession, gene structure for each OsPTR gene and characteristics of corresponding proteins can be found in Additional file 1. Protein structures predicted by SMART are shown in Additional file 2. The putative transmembrane (TM) regions in each OsPTR protein are listed in Additional file 3.
Chromosomal localization and gene duplication
Chromosomal localization and gene duplication events of OsPTR genes. White ellipses on the chromosomes (vertical bars) indicate the position of centromeres. The scale on the left is in megabases (Mb). Respective chromosome numbers are indicated at the top of each bar. The chromosome order has been arranged to bring duplicated regions in the vicinity. The segmental duplication genes have been connected by straight line. In case of tandemly duplicated genes present in clusters, the genes are marked with rectangles of different colors.
Analysis of the MSU RGAP segmental duplication database revealed that 21 OsPTR genes (11 pairs) could be assigned to rice segmental duplication blocks (Figure 1). Os02g47090 participated in two segmental duplications with Os04g50940 and Os10g42870 respectively. The overall identities of the protein sequences of these genes ranged from 37.6% to 76.1%. A total of 44 genes (14 clusters, nine genes also involved in segmental duplication) were considered as tandem duplications according to the criterion adopted in our analysis, including the largest cluster of nine genes on chromosome 1, two four-gene clusters on chromosome 10 and 5, five three-gene clusters on chromosome 3, 4, 6 and 11, six two-gene clusters on chromosome 1, 5, 6, 10 and 12 (Figure 1). The protein homology of these genes varied from 22.9% to 89.5%. Interestingly, all the tandemly duplicated genes in the same cluster had the same direction of transcription except Os06g13200 and Os06g13210. This might suggest the conserved behavior of tandem duplications in this family. Of the 84 OsPTR genes, 66.7% (56 of 84) involved in the duplication events. Moreover, most of the duplicated genes had relative close phylogenetic relationship (see below). Therefore, segmental and large-scale tandem duplication events appeared to have exclusively contributed to the expansion of the OsPTR gene family. Similar manner of member duplication had also been described in some large gene families in rice, for instance, the ankyrin repeat gene family [38] and the bHLH family [39].
Phylogenetic analysis and multiple sequence alignment
Phylogenetic relationship and subfamily formation of OsPTR proteins. Bootstrap values from 1000 replicates are indicated at each node. Scale bar represents 0.1 amino acid substitution per site. The proteins on the tree can be divided into five distinct subfamilies (I-V). Subfamily I, II and III are further divided into seven groups (I-1, I-2, I-3, II-1, II-2 and III-1, III-2). The paralogs whose corresponding genes involved in the duplication are marked with triangles (blue for tandem duplication and red for segmental duplication). The branches of different subfamilies are marked by different colors.
In order to investigate the relationship between the function and phylogenetic subfamily formation of PTR members as well as identify some orthologous genes, a combined phylogenetic tree with OsPTR and AtPTR proteins was also established (Additional file 4). As a result, similar subfamilies were formed compared to the tree of OsPTR. Each subfamily contained both Arabidopsis and rice PTR members. However, most of the members were clustered in species-specific distinct clades, and only four pairs of orthologs (Os01g01360 and At5g13400, Os01g37590 and At2g26690, Os06g15370 and At1g68570, Os05g27010 and At2g40460) could be figured out. This result indicated that the main characteristics of PTR family in rice and Arabidopsis were formed before the split of monocotyledonous and dicotyledonous plants and then evolved separately in a species-specific manner. Moreover, the difference in the total number of OsPTR and AtPTR was mainly due to the variation in the number of subfamily I genes, which were 30 in rice and eight in Arabidopsis. Of the 30 subfamily I PTRs in rice, 21 were tandemly or segmentally duplicated genes.
Phylogenetic relationship and subfamily formation of OsPTR proteins. Identical, conservative and block of similar amino acid residues are shaded in black, charcoal gray and gray, respectively. The TM regions are marked by blue rectangles. Green lines represent the hydrophilic loops in the PTR domain. The red arrows indicate the start and end position of PTR domain. The three conserved motifs are marked by red rectangles.
HMM logos of conserved motifs in OsPTR and AtPTR. X-axis indicates the conserved sequences of the motif. Y-axis is scales of relative entropy which reflect the conservative rate of each amino acid. The names of the motifs are marked under the graphics. (a) Logos of three motifs in OsPTR. (b) Logos of three motifs in AtPTR.
Expression profiling of OsPTRgenes in whole-life cycle of rice
Expression profile analysis of OsPTR genes. (a) Hierarchical cluster display of expression profile for 81 OsPTR genes with corresponding probes in MH63. Color key represents log2 expression values. Samples are mentioned at the bottom: S, germinating seed at 72-h of imbibitions; Se1, Seedling of 3-day after sowing; Se2, Seedling of trefoil stage; Sh, R, shoot and root of seedling with two tillers, respectively; L1, leaf at secondary branch primordium differentiation stage; L2, leaf at meiosis stage; She1, sheath at secondary branch primordium differentiation stage; She2, sheath at meiosis stage; Ste1, stem at 5-day before heading; Ste2, stem at heading stage; FL1, flag leaf at 5-day before heading; FL2, flag leaf at 14-day after flowering; P1, P2, panicle at secondary branch and pistil/stamen primordium differentiation stage, respectively; P3, panicle at pollen-mother cell formation stage; P4, panicle at meiosis stage; P5, panicle at heading stage; H, hull of one day before flowering; Sta, stamen of one day before flowering; Spi, spikelet of 3-day after flowering; E1, E2, E3, endosperm of 7-day, 14-day and 21-day after pollination, respectively. On the left, cluster dendrogram is shown. On the right, five groups are marked with different rectangles. (b) Real-time PCR results of three OsPTR genes which have no probe in MH63. X-axis represents the stages mentioned above. Y-axis represents relative expression level compared with control. Reference samples (marked with asterisk) are used to illustrate the expression pattern more clearly.
Based on hierarchical clustering, expression pattern of OsPTR genes could be divided into five major groups. Group 1 consisted of 25 genes which showed high expression levels in all the tissues expect six genes that showed relatively low expression in certain tissues (Os03g60850 in Sta, Os03g04570 in Ste1, P1, Sta and E3, Os06g15370 in FL1, FL2 and Ste1, Os10g02240 in P1-P3, Os10g33170 in Se, Sta, E1 and E3 and Os12g44100 in P1-P3). Conversely, group 2 contained 10 genes that showed low expression level in the whole life cycle. Group 3 was comprised of 16 genes that showed relatively high expression level in both vegetative tissues and hull. All these genes had low expression level in early panicle and endosperm stages. Group 4 comprised of 16 genes with medium expression level in all the tissues except that Os04g56560 had high expression in root, Os01g65120 in E2, Os11g23890 in E1-E2 and Os01g55600 in E1-E3. The rest 14 genes belonged to group 5 although these genes had different expression patterns. For example, Os04g41400 had high expression signal in root and stamen, Os04g59480 showed high expression in FL2, Os11g12740 had high expression in panicle and stem. Interestingly, Os01g65110, Os01g55610, Os11g17970, Os10g42900 and Os06g21900 all showed relatively high expression in endosperms. The expression patterns of most OsPTR genes in Zhenshan 97 (ZS97) were similar to that in MH63 and the corresponding hierarchical cluster was shown in Additional file 6.
Expression pattern of preferentially and tissue-specifically expressed OsPTR genes in MH63 and ZS97. X-axis indicates representative samples and Y-axis is scales of expression level. Error bars showing the standard errors are based on biological replicates. (a) Preferentially expressed genes. (b) Tissue-specifically expressed genes. Detailed information of samples is listed in legend of Figure 5.
The hierarchical cluster of the AtPTR genes in 24 development stages released the expression profile of 51 genes except At1g72125, At3g45660 and At3g54450 (Additional file 7). Obviously, the expression patterns of AtPTR genes fell into the similar groups as OsPTR genes. Some tissue-specific expression AtPTR genes could also be identified, for example, At1g27040, At4g21680 and At3g25260 in seed and siliques, At3g45710 and At3g45700 in root. Most of the AtPTR genes had relatively low expression level in mature pollen except At1g62200, At5g28470 and At5g01180. These genes were most likely supplying nutrition to maturating pollen. Some of the orthologous OsPTR and AtPTR genes showed the similar expression pattern such as Os06g15370 and At1g68570, Os05g27010 and At2g40460. They both had high expression level in all the tissues. However, no more apparent relationship between the expression pattern and sequence similarities of OsPTR and AtPTR members could be pointed out.
Expression regulation of OsPTRgenes under phytohormone and light/dark treatments
Differential expressions shown by OsPTR genes in response to three phytohormone (NAA, GA3 and KT) treatments. X-axis represents the differentially expressed genes. Y-axis represents average expression values obtained from microarrays. Each treatment is mentioned at the top of each diagram.
Histogram diagrams of differentially expressed OsPTR genes regulated by light/dark treatments in MH63 and ZS97. X-axis represents the differentially expressed genes. Y-axis represents average expression values gained from microarrays. Each treatment is mentioned at the top of each diagram. (a) Differential expression genes in plumule of 48 h after emergence with light/dark treatments. (b) Differential expression genes in radicle of 48 h after emergence with light/dark treatments.
Ka/Ks analysis and express pattern divergence of paralogous OsPTRgenes involved in duplication
Ka/Ks and divergence analysis of paralogous OsPTR genes involved in duplication
Paralog pair | Ka | Ks | Ka/Ks | Purifying selection | Diver- gence model | Duplicate type | Sub- family | Protein identity (%) | Interval genes |
---|---|---|---|---|---|---|---|---|---|
Os02g37040 | 0.119 | 1.333 | 0.0892 | yes | NF | Segmental | V | 76.1 | / |
Os04g39030 | |||||||||
Os02g48570 | 0.247 | 10.084 | 0.0245 | yes | NF | Segmental | I | 55.0 | / |
Os06g21900 | |||||||||
Os03g04570 | 0.137 | 2.312 | 0.0593 | yes | SF | Segmental | II | 70.6 | / |
Os10g33170 | |||||||||
Os01g55600 | 0.412 | 57.191 | 0.0072 | yes | D | Tandem | IV | 44.0 | 0 |
Os01g55610 | |||||||||
Os01g65100 | 0.178 | 1.721 | 0.1035 | yes | NF | Tandem | II | 64.3 | 0 |
Os01g65110 | |||||||||
Os01g65150 | 0.090 | 0.370 | 0.2429 | yes | D | Tandem | II | 82.5 | 0 |
Os01g65169 | |||||||||
Os01g65200 | 0.304 | 1.141 | 0.2662 | yes | NOF | Tandem | II | 30.1 | 0 |
Os01g65210 | |||||||||
Os01g67630 | 0.087 | 0.525 | 0.1653 | yes | NF | Tandem | I | 61.6 | 0 |
Os01g67640 | |||||||||
Os03g13240 | 0.233 | 1.981 | 0.1175 | yes | SF | Tandem | I | 62.6 | 0 |
Os03g13250 | |||||||||
Os04g41400 | 0.153 | 0.975 | 0.1565 | yes | SF | Tandem | III | 75.2 | 3 |
Os04g41450 | |||||||||
Os04g50940 | 0.085 | 0.611 | 0.1397 | yes | SF | Tandem | I | 82.1 | 0 |
Os04g50950 | |||||||||
Os05g33960 | 0.183 | 1.262 | 0.1447 | yes | NOF | Tandem | IV | 62.6 | 3 |
Os05g34000 | |||||||||
Os05g34010 | 0.183 | 3.869 | 0.0473 | yes | NOF | Tandem | IV | 65.4 | 1 |
Os05g34030 | |||||||||
Os05g35594 | 0.247 | 1.884 | 0.1313 | yes | / | Tandem | II | 55.4 | 3 |
Os05g35650 | |||||||||
Os06g13200 | 0.038 | 0.274 | 0.1371 | yes | / | Tandem | II | 89.5 | 0 |
Os06g13210 | |||||||||
Os06g49220 | 0.082 | 0.387 | 0.2124 | yes | NOF | Tandem | I | 79.9 | 2 |
Os06g49250 | |||||||||
Os10g02240 | 0.290 | 1.856 | 0.1565 | yes | NF | Tandem | I | 56.5 | 1 |
Os10g02260 | |||||||||
Os11g18044 | 0.149 | 0.675 | 0.2201 | yes | SF | Tandem | III | 65.4 | 3 |
Os11g18110 | |||||||||
Os12g44100 | 0.177 | 0.592 | 0.2988 | yes | NF | Tandem | IV | 70.9 | 0 |
Os12g44110 |
Expression divergence of paralogous OsPTR genes involved in duplication. X-axis indicates representative samples and Y-axis is scales of expression level. (a) Just formed due to duplication. (b) Nonfunctionalization. (c) Subfunctionalization. (d) Neofunctionalization.
Discussion
Di-/tripeptides and nitrate transport is important for various aspects of plant growth and development, in which peptide transporter family members may play important roles. Considering the potential functional diversity of PTR members associated with nutrient transport and only a few PTR family members have been described in higher plants, it was timely and quite relevant for us to characterize the rice PTR family. Combination of the whole genome sequence data, sequences alignment and high-throughput expression analysis undertaken in this investigation provides an insight into all the 84 members of PTR family in rice.
Duplication played an important role in the expansion of the OsPTRfamily
It is well known that gene duplication is important for gene family evolution and diversity, which can occur via three major mechanisms: segmental duplication, tandem duplication and retroposition [42]. Due to repetitive plant genome duplications, certain plant genomes are likely to harbor relatively large members of a gene family, as compared to animal genomes [43]. The PTR families in rice and Arabidopsis are very much in line with the statement and showed increased family members in order to adapt for specific functions. The number of 84 OsPTR and 53 AtPTR members is far greater as compared to the low number of six in human, four in C. elegans, three in Drosophila and two in yeast [44]. It may indicate the importance of these transporters in plant growth and metabolism. Nearly 70% of the 84 OsPTR genes are duplicated genes, 21 genes (25%) are found to be involved in the segmental duplication while 44 (53.4%) participated in tandem duplication. Our analysis also indicated that in subfamily I, duplication events were largely responsible for the number difference of the member between OsPTR and AtPTR. It demonstrates that the OsPTR family may undergo extensive and rapid genomic changes which can maintain relative large number of member and lead to expression divergences between duplicate genes, further contributing to the establishment of novel gene functions [45].
Phylogenetic relationship and functional diversity among the OsPTR members
The phylogenetic analysis reveals that the OsPTR family can be further divided into five subfamilies. The combined phylogenetic tree of OsPTR and AtPTR showed similar subfamily formation whereas most PTR members were clustered in species-specific clades and only four pairs of orthologs could be identified. It suggests that the members of the PTR family have diverged from each other primarily due to speciation and then undergone significant changes since the monocot/dicot split.
However, highly conserved motifs have been identified in most members of AtPTR and OsPTR. We further adopted the sequence alignment including the PTR members in plant, animal, human and prokaryotic organisms. It was surprising to find that the three motifs could be found in most eukaryotic PTR members but changed a lot or lost in the prokaryotic PTR members, especially motif 2 and motif 3. Moreover, some members such as Os01g65120, Os04g50930, Os12g13790 and At2g38100 which had far phylogenetic relationship from other PTR members, had great amino acid variation in motif 2 and motif 3. We inferred that motif 2 and motif 3 might be important characteristic logos of the eukaryotic PTR members and played considerable roles in the function formation of the PTR family. Further research may employ the point mutations introduced into these motifs, which should cause a loss of function of the corresponding protein. This experimentation could indicate the essentiality of these motifs for the formation of the transporter pore and the interaction with the substrates.
The sequence similarity among the OsPTR members seemed to show little relationship with their putative transporter type classification. In addition, we found that gene structures, intron phases of OsPTR genes and distribution of TM regions in the proteins had no relationship with their phylogenetic subfamily and transporter type classifications. As most OsPTR proteins had PTR as their only recognizable domain (See Additional file 2), we adopted the MEME method [46], in order to identify some new motifs or domains. However, the conserved motifs investigated by MEME were all TM regions (data not shown). We infer that the substrate selection of OsPTR proteins may be decided by their three dimensional structures. The determination of the di-/tripeptide, nitrate and other type transporters in this family by mutagenesis or crystal structure studies should solve the puzzle that di-/tripeptides, nitrate and other substrates of PTR proteins are very different in structure but they share the same transporter family - PTR [6].
Expression pattern divergence of the OsPTRgenes
Microarray provides a rich resource for investigating the evolution of gene expression. The temporal and spatial expression patterns of OsPTR genes are useful information for establishing their putative functions. Our microarray analysis indicated that the expression patterns of the 81 OsPTR genes could be classified into five groups. We also identified some preferential or tissue-specific expression OsPTR genes. Of the 81 genes analyzed in this study, only seven were found to exhibit either preferential or tissue-specific expression during different development stages of MH63 and ZS97. Interestingly, these genes were most preferentially or specifically expressed in root and stamen. Os02g47090, Os04g50940 and Os04g50950 were preferentially expressed in root. Moreover, these three genes were in the same phylogenetic clade that belonged to the group three of subfamily I. In Arabidopsis, the expression of AtNRT1.1 was found concentrated in the tips of primary and lateral roots [26]. AtNRT1.2 was primarily expressed in root hairs and the epidermis in both young and mature regions of roots [27]. OsNRT1 was constitutively expressed in the most external layer of the root, epidermis and root hair [34]. Our data showed that OsNRT1 had high expression levels in all the tissues especially in root. Os03g13240 and Os03g13250 were the paralogs of OsNRT1. Moreover, these three genes formed a tandem duplication cluster. Os03g13240 and Os03g13250 showed similar expression patterns with OsNRT1 and Os03g13250 had the highest expression level in root. It can be inferred that these root-preferentially or -specifically expressed OsPTR genes might be nitrate transporters and play important roles in root formation and development. AtNRT1.4 was expressed in the leaf petiole and participated in regulating leaf nitrate homeostasis [28]. The ortholog of AtNRT1.4 in rice was Os01g37590 which continuously expressed in vegetative tissues with high level in leaves. Expression of Os04g50940 was specific in flag leaf of 14-day after flowering both in MH63 and ZS97. Therefore, we infer that Os01g37590 and Os04g50940 might be participated in the development of leaf as well as the nutrient transporting in and/or out of leaf.
In addition, we also identified three genes (Os01g65110, Os04g41410 and Os12g44100) that were preferentially expressed in stamen and one gene (Os01g55600) in endosperm. One interesting phenomenon in our dendrogram was that three-quarter of the OsPTR genes had significant low expression level in panicles of different development stages. SP1 whose transcription accumulated in panicle was one of the exceptions. SP1 was highly expressed in the phloem of the branches of young panicles, which was consistent with the predicted function of SP1 and the phenotype of sp1 mutant [35]. It can be demonstrated that some PTR members might also play roles in nutrient transport in reproductive organ development.
Potential involvement of OsPTR members in phytohormone and light/dark regulations
Little attention has been paid on the functions of PTR members in phytohormone and light/dark regulations. Until now, AtNRT1.1 was the only reported PTR gene that was involved in these regulations. AtNRT1.1 was regulated by auxin in both shoots and roots and involved in light-induced stomatal opening [47, 48]. However, our microarray data showed that the expression patterns of more than one-third OsPTR genes changed significantly under phytohormone and light/dark treatments. Fourteen OsPTR genes are regulated by at least one phytohormone and 17 OsPTR genes are involved in the regulation of light/dark treatment in MH63 and ZS97. From these results, it can be speculated that these differentially expressed OsPTR genes may play an important role in the transport of regulatory molecules (e.g. hormone-peptide conjugates) and other metabolites induced by regulation of phytohormone and light/dark treatments. For example, IAA is an important plant growth phytohormone and 95% of the IAA pool is sequestered in an inactivated conjugated form, e.g. IAA can be conjugated to single amino acids, peptides or proteins via amide bonds. These conjugates may be involved in regulating phytohormone activity in response to plant growth and development, as well as facilitating hormone transport [49]. It seems that the proteins of those OsPTR genes regulated by NAA might have IAA conjugates transport activity.
Fate of the paralogous OsPTRgenes after duplication
From the Ka/Ks analysis of the paralogous OsPTR genes, it is obvious that purifying selection plays an important role in function maintenance of this family. However, comparison of the expression data implied that after duplication, the paralogous pairs may have been evolved in different hypotheses to gain functional divergence. Expression data shows that some of these duplicated genes have undergone subfunctionalization and neofunctionalization, which contribute to the possibility of maintenance and divergence of a paralogous set of genes. Some paralogs have adopted nonfunctionalization which contributes to the likelihood of loss of one member in a pair of paralogous genes. However, sufficient experimentation is required to prove these hypotheses.
Conclusion
In conclusion, the results of this study display the genomic framework, characterization, conserved motifs, duplication manner and evolutionary analysis of the 84 OsPTR members, along with their expression profiles during whole-life stages as well as under phytohormone and light/dark treatments. These data will provide the basis for understanding the evolutionary history of PTR members and their roles in rice growth and development. The findings in our work would be useful in selecting candidate genes for functional validation studies of PTR members in rice. However, future research by adopting RNAi/overexpression strategies or insertion mutagenesis is required to elucidate the precise role of individual OsPTR gene.
Methods
Database screening and mining
Three methods were used to identify all putative PTR members in rice. First, domain search was executed on website of MSU RGAP http://rice.plantbiology.msu.edu/domain_search.shtml. Second, key words "PTR", "peptide transporter" and "proton-dependent oligopeptide transport" were used as queries to search against the annotation database of rice genome (Release 6.1 of MSU RGAP). Third, BLASTP and TBLASTN searches with the conserved sequence of PTR download from Pfam http://pfam.sanger.ac.uk/ as query were carried out in three databases: MSU RGAP http://rice.plantbiology.msu.edu/, National Center for Biotechnology Information (NCBI, http://www.ncbi.nlm.nih.gov/) and the Knowledge-Based Oryza Molecular Biological Encyclopedia (KOME, http://cdna01.dna.affrc.go.jp/cDNA). The parameters of blast searches were set as follows: max target sequences-500 and expect value less than 10. The Pfam search (E value < 1.0) and SMART http://smart.embl-heidelberg.de/ search (normal mode, including outlier homologues and homologues of known structure, Pfam domains, signal peptides, internal repeats and intrinsic protein disorder) were used to confirm the presence of PTR domain in each predicted protein. Information about the transcript, chromosomal position, full-length cDNA, BAC accession, gene structure for each gene and characteristics of proteins were obtained from MSU RGAP. The structure of each OsPTR protein was predicted in SMART. To predict the putative TM regions in each OsPTR protein, the TMHMM Server v2.0 http://www.cbs.dtu.dk/services/TMHMM/ was applied with default settings.
Chromosomal localization and gene duplication
The physical positions of OsPTR genes were used to map these genes onto the corresponding rice chromosomes. The distribution of OsPTR genes on chromosomes was drawn by using the MapChart software [50] and modified manually with annotation. The presence of OsPTR genes on duplicated chromosomal segments was investigated by searching the segmental genome duplication of rice at MSU RGAP http://rice.plantbiology.msu.edu/segmental_dup/index.shtml with the maximum length distance permitted between collinear gene pairs of 500 kb. Gene distributed nearby and separated by five or fewer genes were considered to be tandem duplicates.
Phylogenetic analysis and sequence alignment
ClustalX version 1.83 [51] was used to perform the multiple sequence alignment and generate the unrooted phylogenetic tree based on the protein sequences of OsPTR members with neighbor-joining method [52]. Bootstrap analysis was performed using 1000 replicates. The phylogenetic tree was displayed and annotated using MEGA software version 4 [53]. The combined tree with OsPTR and AtPTR proteins was generated using the same method. The multiple sequence alignment generated by ClustalX was modified manually and the conserved sequences were annotated according to the results of SMART and Pfam searches. The sequence alignment of the conserved motifs was submitted to the hmmbuild program of Mobyle [54] to build a profile HMM. Then the profile HMM was submitted into LogoMat-M [55] to draw the HMM logo.
Rice materials and treatments
RNA samples used for microarray hybridization were obtained from two elite hybrid rice parents MH63 and ZS97 growing in normal field conditions at different developmental stages. For phytohormone treatments, seedlings at trefoil stage of MH63 and ZS97 were both treated with solutions of 0.1 mM NAA, GA3 and KT respectively. Samples were collected at the time points of 5, 15, 30 and 60 minutes after treatment and the samples under the same phytohormone treatment of different time points were mixed together. For light/dark treatments, seedlings at the plumule and radicle stages of MH63 and ZS97 were placed under 48 hours continuous light or darkness (when plumule was about 1-2 cm and radicle was about 2-3 cm in length) after emergence and harvested at the assigned time respectively.
Microarray expression data analysis
Expression profile data of OsPTR genes was extracted from CREP database http://crep.ncpgr.cn of our lab that composed of hybridization of RNA samples from 39 tissues/organs covering the whole life cycle of rice with the Affymetrix rice microarray. The microarray data have been submitted under the accession number GSE19024 [56] into the NCBI Gene Expression Omnibus (GEO). Thirteen vegetative and 11 reproductive samples of different developmental stages that could cover the whole life cycle of rice were used for OsPTR expression profile analysis in this study. The detailed information of samples used in analysis is listed in Additional file 9.
Probe sets with the entire set of 11 probes (8-10 in some cases) present on the array with 100% identity over the corresponding genes were considered to be significant. Wherever more than one probe set was available for one gene, the probe set located on the 3' end was given preference. After normalization and variance stabilization, the average signal value of replicates for each probe set was used for analysis. Expression values of each gene were logarithmized in hierarchical cluster analysis and the heatmap were performed using R version 2.9.0 http://www.r-project.org/ with the method of complete linkage. The expression profile data of AtPTR genes was extracted from AtGenExpress (http://www.weigelworld.org/resources/microarray/AtGenExpress, [57]). Hierarchical cluster analysis and the heatmap drawing were used the same method as OsPTR.
To identify preferentially expressed genes, a student-t test was adopted. A gene was considered as preferential expression gene only if the signal value of the gene in a designated sample was more than two-fold and had a P value less than 0.05 as compared with the other samples. Also, a check of our whole genome expression profile data of rice for different tissues was used to identify tissue-specific expression genes [56]. Under phytohormone and light/dark treatments, genes that were up- or down-regulated more than two-fold with P value less than 0.05 as compared to control were considered as differential expression genes.
Quantitative real-time PCR analysis
Real-time PCR reactions were carried out using the same RNA samples that were used for microarray. The sequences of primers used for real-time PCR are: Os01g65140 (F: AACTGGCTTCTGGCGAAAGAA, R: CGGAAGCAGCTTGACGAATAA), Os05g35650 (F: GGCAGCTGGTTCTCCAACAA, R: CGGCGATCAGCCAGTAGAAG), Os06g13200 (F: GCTCATGTTCGCGGTCATCT, R: TCCATCAACATCCCCTGCTT) and Actin1 (F: TGGCATCTCTCAGCACATTCC, R: TGCACAATGGATGGGTCAGA). First-stand cDNA was synthesized from DNaseI-treated total RNA using Superscript III reverse transcriptase (Invitrogen). Real-time PCR were performed in 25 μl reaction mixture with 1 μl of first-strand cDNA, 12.5 μl 2× SYBR Premix® Ex Taq™ (TaKaRa), 0.5 μl 50× ROX reference dye II and 0.5 μl of each primer. The reactions were carried out using ABI PRISM 7500 system (Applied Biosystems). The following program was used: 95°C for 10 s; 45 cycles of 95°C for 5 s, 60°C for 34 s. Rice Actin1 was used as internal control. The relative expression levels were analyzed as described previously [58].
Ka and Kscalculations
The paralogous PTR genes in rice were inferred from the phylogenetic tree analysis. Pairwise aligned paralogous protein sequences generated by the ClustalW2 http://www.ebi.ac.uk/Tools/clustalw2/ were used as the alignment guides. PAL2NAL (http://www.bork.embl.de/pal2nal/, [59]) was used to convert the multiple sequence alignment of proteins and the corresponding mRNA sequences into a codon alignment. Then the calculation of non-synonymous (Ka) and synonymous (Ks) substitution rates were carried out on the PAL2NAL website automatically by using the codeml program in PAML [60]. By performing the EMBOSS Pairwise Alignment Algorithms on EBI http://www.ebi.ac.uk/Tools/emboss/align/, we analyzed the protein identities between paralogs of OsPTR.
Declarations
Acknowledgements
We thank Professor Qifa Zhang for his encouragement and direction. We also are grateful to the reviewers and editors for their helpful comments and criticisms. This research was supported by grants from the National Special Key Project of China on Functional Genomics of Major Plants and Animals, and the National Natural Science Foundation of China.
Authors’ Affiliations
References
- Saier MH: Families of transmembrane transporters selective for amino acids and their derivatives. Microbiology. 2000, 146: 1775-1795.PubMedView ArticleGoogle Scholar
- Waterworth Wm, Bray CM: Enigma Variations for Peptides and Their Transporters in Higher Plants. Ann Bot. 2006, 98: 1-8. 10.1093/aob/mcl099.PubMedPubMed CentralView ArticleGoogle Scholar
- Steiner HY, Naider F, Becker JM: The PTR family: a new group of peptide transporters. Mol Microbiol. 1995, 16: 825-834. 10.1111/j.1365-2958.1995.tb02310.x.PubMedView ArticleGoogle Scholar
- Paulsen IT, Skurray RA: The POT family of transport proteins. Trends Biochem Sci. 1994, 19: 404-10.1016/0968-0004(94)90087-6.PubMedView ArticleGoogle Scholar
- Saier MH, Tran CV, Barabote RD: TCDB: the Transporter Classification Database for membrane transport protein analyses and information. Nucleic Acids Res. 2006, 34: D181-186. 10.1093/nar/gkj001.PubMedPubMed CentralView ArticleGoogle Scholar
- Tsay YF, Chiu CC, Tsai CB, Ho CH, Hsu PK: Nitrate transporters and peptide transporters. FEBS Letters. 2007, 581: 2290-2300. 10.1016/j.febslet.2007.04.047.PubMedView ArticleGoogle Scholar
- Frommer WB, Hummel S, Rentsch D: Cloning of an Arabidopsis histidine transporting protein related to nitrate and peptide transporters. FEBS Letters. 1994, 347: 185-189. 10.1016/0014-5793(94)00533-8.PubMedView ArticleGoogle Scholar
- Rentsch D, Laloi M, Rouhara I, Schmelzer E, Delrot S, Frommer WB: NTR1 encodes a high affinity oligopeptide transporter in Arabidopsis. FEBS Letters. 1995, 370: 264-268. 10.1016/0014-5793(95)00853-2.PubMedView ArticleGoogle Scholar
- Song W, Steiner HY, Zhang L, Naider F, Stacey G, Becker JM: Cloning of a second Arabidopsis peptide transport gene. Plant Physiol. 1996, 110: 171-178. 10.1104/pp.110.1.171.PubMedPubMed CentralView ArticleGoogle Scholar
- Dietrich D, Hammes U, Thor K, Suter-Grotemeyer M, Fluckiger R, Slusarenko AJ, Ward JM, Rentsch D: AtPTR1, a plasma membrane peptide transporter expressed during seed germination and in vascular tissue of Arabidopsis. Plant J. 2004, 40: 488-499. 10.1111/j.1365-313X.2004.02224.x.PubMedView ArticleGoogle Scholar
- Song W, Koh S, Czako M, Marton L, Drenkard E, Becker JM, Stacey G: Antisense expression of the peptide transport gene AtPTR2-B delays flowering and arrests seed development in transgenic Arabidopsis plants. Plant Physiol. 1997, 114: 927-935. 10.1104/pp.114.3.927.PubMedPubMed CentralView ArticleGoogle Scholar
- Karim S, Lundh D, Holmstrom KO, Mandal A, Pirhonen M: Structural and functional characterization of AtPTR3, a stress-induced peptide transporter of Arabidopsis. J Mol Model. 2005, 11: 226-236. 10.1007/s00894-005-0257-6.PubMedView ArticleGoogle Scholar
- Karim S, Holmstrom KO, Mandal A, Dahl P, Hohmann S, Brader G, Palva ET, Pirhonen M: AtPTR3, a wound-induced peptide transporter needed for defence against virulent bacterial pathogens in Arabidopsis. Planta. 2007, 225: 1431-1445. 10.1007/s00425-006-0451-5.PubMedView ArticleGoogle Scholar
- Komarova NY, Thor K, Gubler A, Meier S, Dietrich D, Weichert A, Suter Grotemeyer M, Tegeder M, Rentsch D: AtPTR1 and AtPTR5 transport dipeptides in planta. Plant Physiol. 2008, 148: 856-869. 10.1104/pp.108.123844.PubMedPubMed CentralView ArticleGoogle Scholar
- West CE, Waterworth WM, Stephens SM, Smith CP, Bray CM: Cloning and functional characterisation of a peptide transporter expressed in the scutellum of barley grain during the early stages of germination. Plant J. 1998, 15: 221-229. 10.1046/j.1365-313X.1998.00199.x.PubMedView ArticleGoogle Scholar
- Waterworth WM, West CE, Bray CM: The barley scutellar peptide transporter: biochemical characterization and localization to the plasma membrane. J Exp Bot. 2000, 51: 1201-1209. 10.1093/jexbot/51.348.1201.PubMedView ArticleGoogle Scholar
- Waterworth WM, Ashley MK, West CE, Sunderland PA, Bray CM: A role for phosphorylation in the regulation of the barley scutellar peptide transporter HvPTR1 by amino acids. J Exp Bot. 2005, 56: 1545-1552. 10.1093/jxb/eri149.PubMedView ArticleGoogle Scholar
- Schulze W, Frommer WB, Ward JM: Transporters for ammonium, amino acids and peptides are expressed in pitchers of the carnivorous plant Nepenthes. Plant J. 1999, 17: 637-646. 10.1046/j.1365-313X.1999.00414.x.PubMedView ArticleGoogle Scholar
- Miranda M, Borisjuk L, Tewes A, Dietrich D, Rentsch D, Weber H, Wobus U: Peptide and amino acid transporters are differentially regulated during seed development and germination in faba bean. Plant Physiol. 2003, 132: 1950-1960. 10.1104/pp.103.024422.PubMedPubMed CentralView ArticleGoogle Scholar
- Paungfoo-Lonhienne C, Schenk PM, Lonhienne TG, Brackin R, Meier S, Rentsch D, Schmidt S: Nitrogen affects cluster root formation and expression of putative peptide transporters. J Exp Bot. 2009, 60: 2665-2776. 10.1093/jxb/erp111.PubMedPubMed CentralView ArticleGoogle Scholar
- Tsay YF, Schroeder JI, Feldmann KA, Crawford NM: The herbicide sensitivity gene CHL1 of Arabidopsis encodes a nitrate-inducible nitrate transporter. Cell. 1993, 72: 705-713. 10.1016/0092-8674(93)90399-B.PubMedView ArticleGoogle Scholar
- Huang NC, Chiang CS, Crawford NM, Tsay YF: CHL1 encodes a component of the low-affinity nitrate uptake system in Arabidopsis and shows cell type-specific expression in roots. Plant Cell. 1996, 8: 2183-2191. 10.1105/tpc.8.12.2183.PubMedPubMed CentralView ArticleGoogle Scholar
- Liu KH, Huang CY, Tsay YF: CHL1 is a dual-affinity nitrate transporter of Arabidopsis involved in multiple phases of nitrate uptake. Plant Cell. 1999, 11: 865-874. 10.1105/tpc.11.5.865.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang R, Liu D, Crawford NM: The Arabidopsis CHL1 protein plays a major role in high-affinity nitrate uptake. Proc Natl Acad Sci USA. 1998, 95: 15134-15139. 10.1073/pnas.95.25.15134.PubMedPubMed CentralView ArticleGoogle Scholar
- Liu KH, Tsay YF: Switching between the two action modes of the dual-affinity nitrate transporter CHL1 by phosphorylation. EMBO J. 2003, 22: 1005-1013. 10.1093/emboj/cdg118.PubMedPubMed CentralView ArticleGoogle Scholar
- Guo FQ, Wang R, Chen M, Crawford NM: The Arabidopsis dual-affinity nitrate transporter gene AtNRT1.1 (CHL1) is activated and functions in nascent organ development during vegetative and reproductive growth. Plant Cell. 2001, 13: 1761-1777. 10.1105/tpc.13.8.1761.PubMedPubMed CentralView ArticleGoogle Scholar
- Huang NC, Liu KH, Lo HJ, Tsay YF: Cloning and functional characterization of an Arabidopsis nitrate transporter gene that encodes a constitutive component of low-affinity uptake. Plant Cell. 1999, 11: 1381-1392. 10.1105/tpc.11.8.1381.PubMedPubMed CentralView ArticleGoogle Scholar
- Chiu CC, Lin CS, Hsia AP, Su RC, Lin HL, Tsay YF: Mutation of a nitrate transporter, AtNRT1:4, results in a reduced petiole nitrate content and altered leaf development. Plant Cell Physiol. 2004, 45: 1139-1148. 10.1093/pcp/pch143.PubMedView ArticleGoogle Scholar
- Lin SH, Kuo HF, Canivenc G, Lin CS, Lepetit M, Hsu PK, Tillard P, Lin HL, Wang YY, Tsai CB, Gojon A, Tsay YF: Mutation of the Arabidopsis NRT1.5 nitrate transporter causes defective root-to-shoot nitrate transport. Plant Cell. 2008, 20: 2514-2528. 10.1105/tpc.108.060244.PubMedPubMed CentralView ArticleGoogle Scholar
- Almagro A, Lin SH, Tsay YF: Characterization of the Arabidopsis nitrate transporter NRT1.6 reveals a role of nitrate in early embryo development. Plant Cell. 2008, 20: 3289-3299. 10.1105/tpc.107.056788.PubMedPubMed CentralView ArticleGoogle Scholar
- Sugiura M, Georgescu MN, Takahashi M: A nitrite transporter associated with nitrite uptake by higher plant chloroplasts. Plant Cell Physiol. 2007, 48: 1022-1035. 10.1093/pcp/pcm073.PubMedView ArticleGoogle Scholar
- Zhou JJ, Theodoulou FL, Muldin I, Ingemarsson B, Miller AJ: Cloning and functional characterization of a Brassica napus transporter that is able to transport nitrate and histidine. J Biol Chem. 1998, 273: 12017-12023. 10.1074/jbc.273.20.12017.PubMedView ArticleGoogle Scholar
- Jeong J, Suh S, Guan C, Tsay YF, Moran N, Oh CJ, An CS, Demchenko KN, Pawlowski K, Lee Y: A nodule-specific dicarboxylate transporter from alder is a member of the peptide transporter family. Plant Physiol. 2004, 134: 969-978. 10.1104/pp.103.032102.PubMedPubMed CentralView ArticleGoogle Scholar
- Lin CM, Koh S, Stacey G, Yu SM, Lin TY, Tsay YF: Cloning and functional characterization of a constitutively expressed nitrate transporter gene, OsNRT1, from rice. Plant Physiol. 2000, 122: 379-388. 10.1104/pp.122.2.379.PubMedPubMed CentralView ArticleGoogle Scholar
- Li S, Qian Q, Fu Z, Zeng D, Meng X, Kyozuka J, Maekawa M, Zhu X, Zhang J, Li J, Wang Y: Short panicle1 encodes a putative PTR family transporter and determines rice panicle size. Plant J. 2009, 58: 592-605. 10.1111/j.1365-313X.2009.03799.x.PubMedView ArticleGoogle Scholar
- Zhang Q: Strategies for developing Green Super Rice. Proc Natl Acad Sci USA. 2007, 104: 16402-16409. 10.1073/pnas.0708013104.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang Q, Li J, Xue Y, Han B, Deng XW: Rice 2020: A Call for an international coordinated effort in rice functional genomics. Mol Plant. 2008, 1: 715-719. 10.1093/mp/ssn043.PubMedView ArticleGoogle Scholar
- Huang J, Zhao X, Yu H, Ouyang Y, Wang L, Zhang Q: The ankyrin repeat gene family in rice: genome-wide identification, classification and expression profiling. Plant Mol Biol. 2009, 71: 207-226. 10.1007/s11103-009-9518-6.PubMedView ArticleGoogle Scholar
- Li X, Duan X, Jiang H, Sun Y, Tang Y, Yuan Z, Guo J, Liang W, Chen L, Yin J, Ma H, Wang J, Zhang D: Genome-wide analysis of basic/helix-loop-helix transcription factor family in rice and Arabidopsis. Plant Physiol. 2006, 141: 1167-1184. 10.1104/pp.106.080580.PubMedPubMed CentralView ArticleGoogle Scholar
- Hurst LD: The Ka/Ks ratio: diagnosing the form of sequence evolution. Trends in Genetics. 2002, 18: 486-487. 10.1016/S0168-9525(02)02722-1.PubMedView ArticleGoogle Scholar
- Duarte JM, Cui L, Wall PK, Zhang Q, Zhang X, Leebens-Mack J, Ma H, Altman N, dePamphilis CW: Expression pattern shifts following duplication indicative of subfunctionalization and neofunctionalization in regulatory genes of Arabidopsis. Mol Biol Evol. 2006, 23: 469-478. 10.1093/molbev/msj051.PubMedView ArticleGoogle Scholar
- Kong H, Landherr LL, Frohlich MW, Leebens-Mack J, Ma H, dePamphilis CW: Patterns of gene duplication in the plant SKP1 gene family in angiosperms: evidence for multiple mechanisms of rapid gene birth. Plant J. 2007, 50: 873-885. 10.1111/j.1365-313X.2007.03097.x.PubMedView ArticleGoogle Scholar
- Lockton S, Gaut BS: Plant conserved non-coding sequences and paralogue evolution. Trends Genet. 2005, 21: 60-65. 10.1016/j.tig.2004.11.013.PubMedView ArticleGoogle Scholar
- Chiang CS, Stacey G, Tsay YF: Mechanisms and functional properties of two peptide transporters, AtPTR2 and fPTR2. J Biol Chem. 2004, 279: 30150-30157. 10.1074/jbc.M405192200.PubMedView ArticleGoogle Scholar
- Moore RC, Purugganan MD: The evolutionary dynamics of plant duplicate genes. Curr Opin Plant Biol. 2005, 8: 122-128. 10.1016/j.pbi.2004.12.001.PubMedView ArticleGoogle Scholar
- Bailey TL, Williams N, Misleh C, Li WW: MEME: discovering and analyzing DNA and protein sequence motifs. Nucleic Acids Res. 2006, 34: W369-373. 10.1093/nar/gkl198.PubMedPubMed CentralView ArticleGoogle Scholar
- Guo FQ, Wang R, Crawford NM: The Arabidopsis dual-affinity nitrate transporter gene AtNRT1.1 (CHL1) is regulated by auxin in both shoots and roots. J Exp Bot. 2002, 53: 835-844. 10.1093/jexbot/53.370.835.PubMedView ArticleGoogle Scholar
- Guo FQ, Young J, Crawford NM: The nitrate transporter AtNRT1.1 (CHL1) functions in stomatal opening and contributes to drought susceptibility in Arabidopsis. Plant Cell. 2003, 15: 107-117. 10.1105/tpc.006312.PubMedPubMed CentralView ArticleGoogle Scholar
- Bajguz A, Piotrowska A: Conjugates of auxin and cytokinin. Phytochemistry. 2009, 70: 957-969. 10.1016/j.phytochem.2009.05.006.PubMedView ArticleGoogle Scholar
- Voorrips RE: MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered. 2002, 93: 77-78. 10.1093/jhered/93.1.77.PubMedView ArticleGoogle Scholar
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25: 4876-4882. 10.1093/nar/25.24.4876.PubMedPubMed CentralView ArticleGoogle Scholar
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.PubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.PubMedView ArticleGoogle Scholar
- Neron B, Menager H, Maufrais C, Joly N, Maupetit J, Letort S, Carrere S, Tuffery P, Letondal C: Mobyle: a new full web bioinformatics framework. Bioinformatics. 2009, 25: 3005-3011. 10.1093/bioinformatics/btp493.PubMedPubMed CentralView ArticleGoogle Scholar
- Schuster-Bockler B, Schultz J, Rahmann S: HMM Logos for visualization of protein families. BMC Bioinformatics. 2004, 5: 7-10.1186/1471-2105-5-7.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang L, Xie W, Chen Y, Tang W, Yang J, Ye R, Liu L, Lin Y, Xu C, Xiao J, Zhang Q: A dynamic gene expression atlas covering the entire life cycle of rice. The Plant Journal. 2010, 61: 752-766. 10.1111/j.1365-313X.2009.04100.x.PubMedView ArticleGoogle Scholar
- Schmid M, Davison TS, Henz SR, Pape UJ, Demar M, Vingron M, Scholkopf B, Weigel D, Lohmann JU: A gene expression map of Arabidopsis thaliana development. Nat Genet. 2005, 37: 501-506. 10.1038/ng1543.PubMedView ArticleGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.PubMedView ArticleGoogle Scholar
- Suyama M, Torrents D, Bork P: PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic Acids Res. 2006, 34: W609-612. 10.1093/nar/gkl315.PubMedPubMed CentralView ArticleGoogle Scholar
- Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24: 1586-1591. 10.1093/molbev/msm088.PubMedView ArticleGoogle Scholar
Copyright
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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.