Skip to main content
  • Research article
  • Open access
  • Published:

Genome-wide investigation and expression analysis suggest diverse roles and genetic redundancy of Pht1 family genes in response to Pi deficiency in tomato

Abstract

Background

Phosphorus (P) deficiency is one of the major nutrient stresses limiting plant growth. The uptake of P by plants is well considered to be mediated by a number of high-affinity phosphate (Pi) transporters belonging to the Pht1 family. Although the Pht1 genes have been extensively identified in several plant species, there is a lack of systematic analysis of the Pht1 gene family in any solanaceous species thus far.

Results

Here, we report the genome-wide analysis, phylogenetic evolution and expression patterns of the Pht1 genes in tomato (Solanum lycopersicum). A total of eight putative Pht1 genes (LePT1 to 8), distributed on three chromosomes (3, 6 and 9), were identified through extensive searches of the released tomato genome sequence database. Chromosomal organization and phylogenetic tree analysis suggested that the six Pht1 paralogues, LePT1/3, LePT2/6 and LePT4/5, which were assigned into three pairs with very close physical distance, were produced from recent tandem duplication events that occurred after Solanaceae splitting with other dicot families. Expression analysis of these Pht1 members revealed that except LePT8, of which the transcript was undetectable in all tissues, the other seven paralogues showed differential but partial-overlapping expression patterns. LePT1 and LePT7 were ubiquitously expressed in all tissues examined, and their transcripts were induced abundantly in response to Pi starvation; LePT2 and LePT6, the two paralogues harboring identical coding sequence, were predominantly expressed in Pi-deficient roots; LePT3, LePT4 and LePT5 were strongly activated in the roots colonized by arbuscular mycorrhizal fungi under low-P, but not high-P condition. Histochemical analysis revealed that a 1250-bp LePT3 promoter fragment and a 471-bp LePT5 promoter fragment containing the two elements, MYCS and P1BS, were sufficient to direct the GUS reporter expression in mycorrhizal roots and were limited to distinct cells harboring AM fungal structures. Additionally, the four paralogues, LePT1, LePT2, LePT6 and LePT7, were very significantly down-regulated in the mycorrhizal roots under low Pi supply condition.

Conclusions

The results obtained from this study provide new insights into the evolutionary expansion, functional divergence and genetic redundancy of the Pht1 genes in response to Pi deficiency and mycorrhizal symbiosis in tomato.

Background

Phosphorus (P) is one of the three most essential macronutrients required by plants. It is well recognized as serving a wide range of structural and biological roles, such as energy metabolism, signal transduction, biosynthesis of macromolecules, modulation of respiration, photosynthesis and other metabolic processes [1]. The primary source for P uptake by plants is orthophosphate (Pi) in soil. Due to the slow diffusion rate and chemical fixation, P is widely considered to be one of the most difficult nutrients for plants to forage, and often a major limiting factor to crop yields [2, 3].

The Pi concentration in soil solution is commonly no more than 10 μM, whereas plant cells need to maintain their cytoplasmic Pi concentrations at a millimolar range [4, 5], which determines the requirement of metabolic energy and specific transport systems for plants to acquire Pi from soils [6, 7]. In the past decades, aided by systematic studies of molecular biology and functional genomics of model plants, a great deal of knowledge about the mechanisms of Pi transport by plants has been accumulated, revealing that the uptake and subsequently redistribution of Pi within plants are mediated by a number of phosphate (Pi) transporters with different affinities that located in the plasma or organelle membranes [8].

The first gene encoding plant Pi transporter (AtPT1) was isolated from Arabidopsis[9] and showed high sequence identity to the genes encoding high-affinity Pi transporters in Saccharomyces cerevisiae (Pho84) [10] and in Glomus versiforme (GvPT) [11]. The later studies further led to the isolation of other eight homologues exhibiting substantial identities to the AtPT1 in the Arabidopsis genome [12], suggesting the expansion of Pi transporter genes in higher plants during evolution. By now, with the completion of whole genome analysis of model plants, such as Arabidopsis and rice, dozens of homologous genes encoding different affinities and groups of Pi transporters have been identified in various plant species by comparative genomic approaches [13]. Studies on the protein sequences and phylogenetic relatedness revealed that most of the Pi transporters identified so far are typical of H+/Pi symporters, and could be grouped into the high-affinity Pht1 family included in the super facilitator superfamily (MFS) [14–16].

Earlier studies on the regulation and tissue/cellular distribution indicated that members of the Pht1 family in many species are divergent in function and differentially expressed during plant development or in response to different P status [17, 18]. The relatively high levels of transcripts or proteins of some Pht1 genes in roots, especially in root epidermis and root hairs, in response to Pi deficiency well support a role of these genes in Pi capture and uptake [19, 20]. For example, in Arabidopsis, eight of the nine Pht1 genes were expressed in roots and two members, AtPT1 and AtPT4, had the highest expression levels in response to Pi deficiency. Knock out of either of the two genes showed significant defects in P uptake under a low Pi supply condition [21, 22]. In some cases, the transcripts of some Pht1 members are more widely distributed throughout plant tissues and showed less responses to Pi deficiency, providing strong evidence to support that some of the Pht1 members may be implicated in the internal mobilization of Pi, such as loading or unloading from the xylem or phloem and deposition into seeds or other storage organs [19, 23–25]. In addition to the Pi-responsive Pht1 genes, an increasing number of arbuscular mycorrhiza-induced Pi transporters belonging to the Pht1 family have been identified from several plant families, and their functions have been repeatedly documented to be associated with Pi uptake at the intraradical symbiotic interface [26–33].

Tomato, a member of the Solanaceae, is not only a world-wide major vegetable crop plant, but also a model plant for biological and genetic researches based on its relatively low-copy DNA sequence and the nearly complete genome sequencing [34]. Although previous studies have characterized the potential roles of a few individual Pht1 members in tomato [30, 35], there is a lack of genome-wide analysis of the Pht1 gene family in tomato and also in any other solanaceous species thus far. Moreover, compared to the other model species, such as Arabidopsis from Brassicaceae and rice from Gramineae, the evolutionary mechanisms, transcriptional regulation and possible functions of solanaceous Pht1 genes in Pi acquisition and mobilization still needs to be well explored [36–39].

In the current work, we reported the genome-wide identification and comparative characterization of Pht1 family genes in tomato and potato, and further investigated the expression patterns of tomato Pht1 genes in response to AM fungi inoculation under low- and high-P supply condition. The analysis in this study mainly focused on the chromosomal organization, phylogenetic evolution, tissue-specific expression and regulation of each member of the tomato Pht1 family. The results obtained from this study would not only strengthen our understanding on the molecular mechanisms underlying the evolutionary expansion, conservation and functional divergence of the Pht1 genes in tomato, but also provide valuable clues for the further comparative genomic studies across the whole Solanaceae family.

Results

Identification of Pht1 family genes in tomato

Previously, five Pht1 genes (three with full-length and two with partial mRNA sequence) encoding for putative high-affinity Phosphate (Pi) transporters (PT) in tomato have been reported [30, 40]. In order to determine whether there are any further members, as yet unidentified, comprising the tomato Pht1 family, the mRNA and amino acid sequences of Arabidopsis and rice Pht1 genes were employed for BLASTN and TBLASTN searches against the recently released tomato genomic sequence database (http://solgenomics.net/), which resulted in the identification of a total of eight non-allelic sequences as the putative tomato Pht1 genes (Additional file 1). BLAST searches of these sequences against the NCBI database demonstrated that five of the eight sequences were identical to the accessioned tomato Pht1 genes, LePT1 to 5. The rest three putative genes (named as LePT6 to 8), which were newly identified in this study, showed high levels of sequence identity to the known Pht1 genes from tomato and other plant species (Table 1). Moreover, LePT6, which represents a distinct locus, harbors its coding sequence identical to that of the known LePT2, but with much difference in un-translated regions between the two homologues (Additional file 2).

Table 1 Identity matrix for the eight putative tomato Pht1 genes

Comparative analysis of the full-length deduced polypeptides revealed that the eight Pht1 proteins contain 528-538 amino acids with 12 predicted transmembrane-spinning domains, similar to the molecular feature of Pht1 transporters from other plant species. Additionally, all the tomato Pht1 amino acid sequences share the consensus sites for phosphorylation by protein kinase C and casein kinase II and conserved residue for N-glycosylation (Figure 1). Using the DNAMAN multiple sequence alignment program, the conserved domain, GGDYPLSATIxSE, which have been suggested to be a typical signature of Pht1 proteins, was also identified in all of these proteins (Figure 1). These findings led to the suggestion that all the identified genes could be considered as tomato Pht1 genes.

Figure 1
figure 1

Predicted amino acid sequences of the eight tomato Pht1 genes LePT1 to 8 . Sequence alignment analysis was carried out using multiple alignment algorithm wrapped within the DNAMAN 7.0 program (http://www.lynnon.com/). Identical amino acids are shaded and gaps are indicated by dots. The consensus sites for phosphorylation by protein kinase C and casein kinase II are shown by two arrowheads with red colour and the conserved N-glycosylation residue is shown by a green arrowhead. The characteristic Pht1 signature was underlined. The transmembrane domains (broken underline) were predicted by the Toppred algorithm (http://bioweb.pasteur.fr/seqanal/interfaces/toppred.html).

A further blast searches against the tomato EST database at NCBI, SGN and TIGR database revealed that except LePT8, the other two newly identified Pht1 genes, LePT6 and LePT7, could matched perfectly to at least one significant EST sequences, indicating that the two genes, like their previously reported five paralogues, are transcriptionally active in a certain tissues. It should also be emphasized that except the eight Pht1 genes mentioned above, another sequence (named as LePTx in this study) identified in the tomato scaffold database searches also showed very high identity to the tomato LePT7 and LePT8 genes, but may be inactive due to the inclusion of some nonsense mutations and indels (insertions and deletions) within its putative coding region (Additional file 3), as well as to the absence of any EST sequence exactly matching.

Identification of tomato Pht1 homologues in potato and comparative analysis of these genes between the two solanaceous genomes

As a near complete set of potato genome sequences were also recently available at the SGN database [41], for further investigating the evolutionary conservation and divergence of Pht1 gene family between the two solanaceous species, the potato genome sequence database were extensively searched using the tomato Pht1 genes as queries, leading to the identification of a total of 10 distinct genes as putative potato Pht1 genes (Additional file 4). Sequence comparison of the potato Pht1 genes revealed similar amino acid sizes and high sequence identities to their corresponding orthologues from tomato (Additional files 4 and 5). It should be noted that there also exist two other sequences in the potato genome that showed substantial homology to the plant Pht1 genes, but may be pseudogenes (named as StPTx1 and StPTx2, respectively), due to the presence of some nonsense mutations and the inclusion of some indels (insertions and deletions) within their putative coding regions.

Similar to the high sequence identity between the tomato and potato Pht1 members, high conservation of chromosomal organization of the Pht1 homologues from the two solanaceous species could also be observed. Figure 2 shows the localizations of Pht1 genes on the tomato and potato chromosomes. It was revealed that the distribution of the tomato and potato Pht1 genes were obviously uneven, and concentrated on only three (3, 6, and 9) chromosomes of the two plants. In addition, the supposed three pseudogenes (LePTx, StPTx1 and StPTx2) were restrictedly assigned on the chromosome 9 of the two plants. Interestingly, except some individual members, such as LePT8 on tomato chromosome 6 and StPT9 on potato chromosome 9, most of the other Pht1 genes/pseudogenes on the corresponding chromosomes were distributed in clusters with very short physical distance (Additional file 1, Figure 2), suggesting that these clustered genes may be produced from independent tandem duplications during the evolution of Solanaceae Pht1 gene family.

Figure 2
figure 2

Distribution of Pht1 genes ( LePT1 to 8 and StPT1 to 10) on the tomato (T) and potato (P) chromosomes. Chromosome numbers are shown at the top of each bar. The arrows next to the gene names indicate the direction of transcription. LePTx, StPTx1 and StPTx2 are putative Pht1 pseudogenes residual in the tomato and potato genomes. The genes from other two families, Pht3 and Pht4, encoding putative Pi transporter or carrier with no homology to the Pht1 proteins, were also labeled on the corresponding chromosomes of the two plants.

In addition to the Pht1 gene themselves, the potential genes surrounding each of the Pht1 members were also carefully surveyed, resulting in the identification of several putative genes exhibiting substantial homology to plant Pht3 and Pht4 family genes on the corresponding chromosomes (Figure 2). By comparing the locations of these genes, such as LePht3;4 and StPht3;4, we confirmed the existence of two segmental inversions associated with the long arms of tomato/potato chromosomes 6 and 9, which resulted in the inverted linear order of the orthologous pairs of PT4/PT5 and PT1/PT3 on the corresponding chromosomal regions between the two species (Figure 2).

Phylogenetic analysis of Pht1 gene family in tomato and other plant species

In order to perform a comprehensive analysis of evolutionary relationships among Pht1 genes between tomato and other plant species, including the eight Pht1 proteins of tomato, a total of 90 plant Pht1 protein sequences, representing 11 species from four plant families, Gramineae, Brassicaceae, Leguminosae and Solanaceae, were aligned and used to construct an unrooted phylogenetic tree. As shown in Figure 3, except AtPT6 from Arabidopsis and HvPT8, OsPT13 and ZmPT5 from each of the three graminaceous species, barley, rice and maize, the other plant Pht1 proteins in the Neighbour-Joining tree were well clustered into four distinct groups, consisting of one dicot-specific group (I), one monocot-specific group (II) and two mixed groups (III and IV, respectively) with members from both dicots and monocots.

Figure 3
figure 3

Phylogenetic analysis of tomato Pht1 genes and other plant Pht1 homologs. An unrooted phylogenetic tree of the plant Pht1 proteins was constructed using the neighbor-joining method with MEGA 5.0 program. Transporters and corresponding plant species are: tomato, LePT1 to 8 [30, 42], this study; potato, StPT1 to 10 [30, 43, 44], this study; tobacco, NtPT1 to 5 [45, 46]; eggplant, SmPT1 to 5 [45, 46]; Arabidopsis thaliana, AtPT1 to 9 [47]; Medicago truncatula, MtPT1 to 6 [8, 48, 49]; Lotus japonicus, LjPT1 to 4 [31, 50]; Soybean, GmPT1 to 14 [13]; Rice, OsPT1 to 13 [28]; Barley, HvPT1 to 12 [14, 51]; Maize, ZmPT1 to 6 [18, 29].

The Group I harbors the proteins exclusively from the dicotyledonous species, in which they are subgrouped by phylogeny, and could be further classified into three subgroups (named as subgroup a, b and c, respectively). In addition, except the subgroup c, which only includes one member from each of the three species, Arabidopsis, Medicago and Lotus japonicus, both of the other two subgroups, a and b, contain multiple Pht1 members from the three plant families, Leguminosae, Brassicaceae and Solanaceae. For tomato, six of the eight Pht1 transporters (LePT1 to 3 and LePT6 to 8) fall into two of the three subgroups. Within subgroup a, the two tomato members, LePT1 and LePT3, group together with their orthologous pairs from potato, eggplant and tobacco, to the exclusion of other two paralogues, LePT7 and LePT8, which group together with other three potato homologues and forms a cluster with two soybean homologues, GmPT6 and GmPT14. Within the dicot subgroup b, the two tomato paralogues, LePT2 and LePT6, group closely, and cluster together with three potato homologues, StPT2, StPT6 and StPT7. The grouping of LePT2 and LePT6 was expected as the two paralogues contain identical coding sequences. The rest two paralogues, LePT4 and LePT5, were found to be assigned only into Group III, although the Group IV also contains the members from both dicots and monocots. Within Group III, the two genes, like their two paralogues LePT1 and LePT3 in Group I, also group together with other solanaceous orthologues by forming an independent Solanaceae clade consisting of two subclasses. Both of the subclasses contain the orthologous pairs of PT4 or PT5 from tomato, potato eggplant and tobacco, suggesting that the duplication events associated with the arising of PT4 and PT5, as well as PT1 and PT3 in tomato and other solanaceous species, occurred before the speciation of Solanaceae lineages. Additionally, most of the Pht1 proteins in the Group III, including the PT4 and PT5 orthologous pairs, have been experimentally evidenced to be strongly induced in the roots colonized by arbuscular mycorrhizal (AM) fungi [8]. Moreover, Pht1 members from Arabidopsis, of which the roots are unable to form AM symbiosis, are all absent from Group III. Interestingly, although the Group IV contains much fewer members as compared with the other three Groups, there exist two members from each of three species, Arabidopsis, rice and soybean, but no homologues from any of the solanaceous species in the Group IV, suggesting that the corresponding orthologues from solanaceous lineages have lost after Solanaceae separation with Brassicaceae and Leguminosae.

Expression analysis of the tomato Pht1 genes in different tissues under low-P condition

In this study, for gaining better understanding of the possible functions of specific Pht1 gene in tomato, the tissue-specific expression patterns of each tomato Pht1 gene were examined in various tissues, including roots, stems, young leaves, flowers, as well as fruits at young and ripe stages using Real-time RT-PCR. The quantitative data showed that except LePT4 and LePT8, of which the transcripts were not detectable in all tissues examined, the transcripts of other Pht1 paralogues were all detectable in a certain tissues and showed distinct but partially overlapping expression profiles (Figure 4).

Figure 4
figure 4

Tissue-specific expression analysis of tomato Pht1 genes. The RNA were prepared from different tissues, including roots (R), stems (S), young leaves (L), flowers (FL), as well as fruits at green (GF) and ripe (RF) stages. The relative expression levels of each of the tomato Pht1 genes were indicated as percentage of the constitutive Actin expression activity. Each bar was the mean of three biological replications with standard error.

LePT1 was expressed in all tissues examined, and its transcripts were detected abundantly in roots and leaves, and to a lesser extent in stems and flowers, as well as in fruits. The transcripts of LePT1 in green fruits were four times more than those in ripe fruits. In contrast to the ubiquitous expression profiles of LePT1, the expression of LePT2 showed relatively distinct tissue-specific profiles, with its transcripts intensively in roots and extremely faintly in some of other tissues, such as in green and ripe fruits. The expression patterns of the LePT3 and LePT5 were a little similar, as both of the two genes were expressed very weakly in all tissues. Even so, the highest transcript level for LePT5 was detected in ripe fruits, and was about ten times more than that in green fruits. LePT6, the closest fellow of LePT2 in phylogeny, was also dominantly expressed in roots, but with only one-third of the expression level of LePT2 in the root tissues. Additionally, very weak transcript levels of this gene were also detectable in stems and leaves. LePT7 was also ubiquitously expressed in all tissues, and had a very similar expression tendency, but significant lower expression levels in all tissues as compared to its paralogue, LePT1 (Figure 4). The differential but overlapping expression of the Pht1 genes well mirrors the evolutionary conservation and functional divergence of Pht1 transporters in tomato plants.

Expression analysis of tomato Pht1 genes in response to AMF colonization under low and high Pi supply conditions

Since the expression of some Pht1 genes in tomato and also in other plant species, have been characterized to be AM-inducible and Pi-responsive [21, 52], the relative expression levels of each tomato Pht1 member were thus further determined in roots and leaves in response to AM Fungi (Glomus intraradices) colonization under low (50 μM) and high (1 mM) Pi supply condition. As shown in Figure 5, colonization of AM fungi increased not only the biomass, but also the P concentration of the tomato plants under the low Pi supply condition; however, no significant difference of both the biomass and P concentration could be observed between the colonized and the noncolonized plants under the high Pi supply condition.

Figure 5
figure 5

Effects of mycorrhizal fungal colonization on tomato biomass (fresh weight, FW) and P concentrations under low Pi (0.05 mM) and high Pi (1 mM) supply condition. –P and + P represent supply of 0.05 mM and 1 mM Pi, respectively; –M and + M represent the inoculation with autoclaved and active inoculum containing arbuscular mycorrhizal fungi Glomus intraradices (Gi), respectively. Error bars indicate SE (n = 3).

qRT-PCR analysis revealed that except the three paralogues, LePT3, LePT4 and LePT5, of which the transcripts were strongly enhanced or specifically activated only in the inoculated roots under the low Pi supply condition, and LePT8, of which the transcripts were not detectable in both tissues under any treatments, the expression of the other four paralogues, LePT1, LePT2, LePT6 and LePT7, were significantly repressed under the high Pi supply condition (Figure 6). Such down-regulated cases occurred more conspicuously upon the two paralogues, LePT2 and LePT7, as their transcripts in both the root and leaf tissues were drastically decreased (LePT2) or even completely absent (LePT7) under high Pi conditions regardless of with or without AM colonization. In addition, very significant decrease of the transcript abundance of the four paralogues was also detected in both the roots and leaves of the colonized tomato plants as compared to those non-colonized controls under low Pi supply condition (Figure 6). The remarkable down regulation of these four members in response to high-P supply and AMF-colonization might be partially caused by the significant increase of P concentration in such treated tomato plants (Figure 5). Interestingly, although LePT2 and LePT6 were considered to be the closest related genes in tomato Pht1 family due to their identical coding sequence, the down regulation of LePT6 in roots in response to AM symbiosis under low Pi condition was much moderate than that of LePT2. Such discrepancy in expression levels strongly suggests that the regulatory components controlling the activation or suppression of LePT2 and LePT6 have divergent after the two paralogues produced from a relatively recent duplication event.

Figure 6
figure 6

Real-time RT-PCR analysis of the Pht1 genes in tomato roots and leaves in response to mycorrhizal fungi colonization under high and low Pi supply conditions. The plants were incubated for two months with a mycorrhizal inoculum containing Glomus intraradices. LP and HP represent supply with 0.05 mM and 1 mM Pi, respectively; -M and + M represent the inoculation with autoclaved and active inoculum, respectively. The relative expression levels of each tomato Pht1 gene was also shown as percentage of the constitutive Actin expression activity. Each bar was the mean of three biological replications with standard error.

The specialized expression profiles of the tomato Pht1 genes in response to AM symbiosis or different Pi status prompted us to investigate their promoter regions. As shown in Figure 7A, the number and localization of the two Pi-regulated (P1BS and W-box) [53, 54] and one AM-responsive elements (MYCS) [45, 55] differ widely in the promoter regions of these eight Pht1 genes, even though the coding sequence and expression profiles of some paralogues are highly conserved. However, similar as the AM-induced Pht1 genes in other dicot species, the MYCS motif was found to be present exclusively in the putative promoter regions of the three AM-activated Pht1 paralogues, LePT3, LePT4 and LePT5, and were located very closely to the Pi-regulated P1BS element [45]. Histochemical staining analysis further revealed that the LePT3 and LePT5 promoter regions (pLePT3 -1250 and pLePT5 -471 ) containing the two elements, MYCS and P1BS, were sufficient to direct β-glucuronidase (GUS) expression specifically in the mycorrhizal roots and were limited to distinct cells harboring AM fungal structures (arbuscules or intracellular hyphae) (Figure 7B), similar to the cellular distributions of their paralogue LePT4 and other AM-inducible Pht1 homologues from various other plant species reported previously [30, 55–57].

Figure 7
figure 7

Analysis of the tomato Pht1 gene promoters. (A) Comparative analysis of putative cis-regulatory elements responsible for the Pi- and AM-regulated expression between the eight tomato Pht1 promoters. Two previously reported Pi-responsive motifs (P1BS and W-box) and one AM-activated motif (MYCS) were searched using the DNA-pattern matching arithmetic (http://rsat.ulb.ac.be/rsat/). P1BS, GNATATNC; MYCS, TTCTTGTTC; W-box, TTGACY. (B) Histochemical analysis for the promoter activity of the two AM-induced Pht1 members, LePT3 and LePT5. (a-d) Localization of β-glucuronidase (GUS) activity (a and b, Magenta GUS; c and d, blue GUS) in mycorrhizal roots driven by the promoters of LePT3 (a, c) and LePT5 (b, d), respectively. (e, f) Co-localization of GUS activity (indicated by the purple color, from the overlay of the Magenta-GUS and Trypan Blue stains) showed that the LePT3and LePT5 promoter fragments (pLePT3 -1250 and pLePT5 -471 ) were sufficient to direct GUS expression in mycorrhizal roots and were confined to distinct cortical cells containing AM fungal structures (arbuscules or intracellular hyphae). Green arrows indicate arbuscule or arbusculate hyphae, yellow arrows indicate intracellular hyphae and red arrows indicate noncolonized cells.

Discussion

In recent studies, benefiting from the availability of whole genome sequence of model plants, dozens of genes belonging to the Pht1 family that encode putative high-affinity Pi transporters have been identified from various plant species using the comparative genome approaches. Tomato, a model species from the Solanaceae family, has been historically characterized to have at least five Pht1 genes [58]. Our present study, through extensive searches of available databases, led to the identification of a total of eight putative Pht1 genes in the tomato genome. As this is the first genome-wide analysis of the Pht1 gene family in any solanaceous species, the investigation of chromosomal organization, evolutionary relationships, as well as expression patterns of the tomato Pht1 genes in this study is of great significance and would offer a basis for better understanding the evolutionary mechanisms underlying the expansion, conservation and functional divergence of the Pht1 genes in the whole Solanaceae family.

Evolutionary expansion of the tomato Pht1 genes

Multigene families, in a general way, could arise through tandem duplications, resulting in a clustered occurrence, or through genome/segmental duplications, resulting in a discrete distribution of family members [59]. As most of the tomato Pht1 genes were assigned in clusters (such as PT1/PT3, PT2/PT6 and PT4/PT5), with not only very close physical localization (Figure 2), but also very high levels of sequence identity (Table 1, Additional file 5), it is strongly suggestive of that tandem duplications might be the major contributors to the expansion of the tomato Pht1 family. Additionally, since most of the tomato Pht1 members group together with their orthologues from other solanaceous species, such as potato, eggplant and tobacco by forming independent solanaceous clades to the exclusion of other dicot homologues (Figure 3), indicates that the duplications associated with the arising of the coupled paralogues such as PT1/PT3 and PT4/PT5 in solanaceous species, occurred before the speciation of solanaceous lineages from a common ancestor. Intriguingly, in viewing of the localization of tomato LePT2/6 and potato StPT2/6/7 on their corresponding chromosomes, it is tempting to make a tendentious conclusion that the duplication giving rise of the tomato LePT2 and LePT6 probably occurred before the split of tomato and potato. However, the distribution of the LePT2 and LePT6 in the terminal subclades of the phylogenetic tree, and the identical coding sequence shared by them well reflected that the two paralogues were produced from the more recent duplication events that occurred within the tomato lineage postdating it split from a common ancestor shared by potato.

It has been recently documented that the Solanum lineage genome has undergone two rounds of consecutive whole-genome triplication events, one that was ancient and shared with most dicot plant families, and one that was more recent and occurred before the divergence of tomato and potato lineages [34], which led to the hypothesis that segmental duplications produced by genome polyploidy may also exert important impact on the expansion of the Pht1 family. The loci of LePT2/6 on chromosome 3 and LePT4/5 on chromosome 6 flanked respectively by two paralogues, LePht3;3 and LePht3;4 (with non-homology to the Pht1 genes) (Figure 2) strongly suggests that the arising of the two pairs, PT2/PT6 and PT4/PT5 might originate from a segmental duplication, followed by two independent tandom duplications, which eventually resulted in the fixation of the two couples of Pht1 members on the chromosomes 3 and 6. As the PT4 and PT5 paralogues cluster together with other members from both dicots and monocots in the phylogenetic tree (Figure 3), indicating that the segmental duplication yielding the precursors of the two couples, PT2/PT6 and PT4/PT5, occurred before the divergence of monocots and dicots. With regard to the other four paralogues, the clustered LePT1 and LePT3, and the individual LePT7 and LePT8, they may be the result of several relatively recent segmental or single-gene duplication events that occurred before the speciation of tomato and potato lineages, and followed by at least one independent tandom duplication event (producing the two paralogues, LePT1 and LePT3). It has been well documented that genome polyploidizations are commonly accompanied by massive chromosomal rearrangements [60]. In our study, by comparing the phylogenetic tree and the chromosomal distribution of the tomato and potato Pht1 genes, two segmental inversions leading to the inconsistent linear orders of PT1/PT3 and PT4/PT5 were identified between the two solanaceous genomes, well supporting the very recent findings that at least nine large and several smaller inversions exist between the tomato and potato lineages [41].

Functional conservation and divergence of the tomato Pht1 gene family

It has been generally accepted that gene duplication followed by functional differentiation has performed a pivotal role in driving evolutionary novelty that allow plants to increase fitness to new environments [47]. To data, as the lack of genome-wide survey of Pht1 genes in any solanaceous species, there is no systematic analysis of tissue-specific expression patterns for the tomato Pht1 family so far. In our present work, we revealed that differential but partial overlapping expression of the Pht1 genes occurred in tomato, as did the members of this family in several other plant species, such as Arabidopsis, rice and soybean [13, 28, 43]. The specialized expression of these genes well mirrors the evolutionary divergence of regulatory elements that are required for controlling Pi uptake and mobilization within/across particular tissues or cells during tomato plant growth.

Earlier results, based on the study of tissue-specific expression and cellular distribution of Pht1 genes in several different plant families revealed that many of the Pht1 genes are expressed dominantly in roots, especially in root epidermis and root hairs, in response to P deprivation, suggesting a potential role of these genes in Pi capture and uptake [12, 19]. The transcription data obtained in this study indeed provide direct evidence for strong expression of most of the tomato Pht1 genes in the roots under low Pi supply condition (Figure 6). It was shown that although the transcripts of LePT1 in the roots and leaves significantly decreased in response to Pi sufficiency, a low level of constitutive expression could be detected throughout the plant, consistent with the expression patterns of its orthologue, StPT1, in potato [46], suggesting that LePT1 and its orthologues may be involved in not only uptake of Pi from soil solution but also redistribution of Pi within plants. LePT2 has been previously documented to be expressed exclusively in P-depleted roots. However, in the present work, a relatively weak but still observable transcription level could be detected in the roots irrigated with high Pi (1 mM) solution, similar results could also be observed from our previous studies on the LePT2 orthologues in other three solanaceous species, eggplant, pepper and tobacco [61]. In addition, very slight levels of the LePT2 transcripts were also detectable in stems, flowers and fruits at green and ripe stages under low Pi supply condition. Although LePT2 shares its coding sequence identical to its paralogue, LePT6, the transcript abundance of the two members were observably different whether under low Pi supply condition or in response to AM symbiosis. Such discrepancy between the two close paralogues may be caused by the inconsistent distributions of Pi-responsive elements, such as P1BS and W-box, in their promoters (Figure 7A) [53, 54, 62–64]. Even so, the identical protein activity and high degree of overlapping expression strongly implies the presence of functional redundancy between the two members.

With regard to the three AM-activated Pht1 paralogues, LePT3, LePT4 and LePT5, as their transcripts could be induced abundantly only in the inoculated roots under low Pi, but not high Pi condition, similar as the expression of their orthologues from tobacco, pepper, and eggplant, suggesting that the AM-activated expression of LePT3-5, like their solanaceous orthologues, might also be regulated by at least two conserved cis-elements, P1BS and MYCS [45]. Comparative screening of the putative promoter regions of the tomato Pht1 genes indeed led to the identification of the two motifs that were closely and exclusively localized in the promoters of the three AM-activated Pht1 genes (Figure 7A). The relatively high proportion of AM-induced Pht1 genes present/retain in solanaceous species compared to that in many of other mycorrhizal plants, such as rice, to a certain extent, reflects the importance of these genes in regulation of symbiosis during the Solanaceae evolution. In addition, the strong expression levels of these symbiosis-activated Pht1 genes in contrast to the remarkable down-regulation of the Pi transporters responsible for direct Pi uptake in mycorrhizal roots also provided strong evidence to support the earlier findings that symbiotic uptake pathway would contribute the majority of the accumulated Pi received by mycorrhizal tomato and other plants [65–67]. It is worth noting that although LePT4 was observed to be the gene that had the highest expression level than the other paralogues in tomato mycorrhizal roots under low Pi condition (Figure 6), mutation of the LePT4 expression in tomato virtually unaffected the establishment of AM symbiosis [30, 40], which seems contradictory to the recent findings that knock down/out of the AM-specific or up-regulated Pht1 genes in Lotus japonicus, Medicago and Astragalus sinicus significantly impaired both the development of AM interaction and symbiotic Pi uptake [31, 55, 68]. The absence of AM-associated phenotypes in the tomato lept4 knock out mutant was thus suggested to be the genetic redundancy within the tomato Pht1 gene family [30, 69]. This explanation might be reasonable as a result of that LePT4 and its AM-induced paralogue, LePT5, were considered to be diverged from a common precursor through tandom duplications (Figure 2), and thus it is unsurprising that LePT5 would share a similar physiological role with LePT4 in tomato. Interestingly, although there exist two AM-specific Pht1 genes, the strongly activated OsPT11 and poorly induced OsPT13, in the rice Pht1 family, no functional redundancy was observed between the two paralogues, as silencing whichever of the two paralogues caused a significant repression of AM symbiosis [70]. Since the rice OsPT11 and OsPT13 are distributed relatively distantly in phylogenetic tree, similar as the phylogenetic relationships between solanaceous PT3 and PT4/PT5 paralogues, it is tempting to raise a speculation that no functional redundancy might be existent between the two groups (PT3 and PT4/PT5) of Pht1 transporters in regulation of the development of AM symbiosis.

Interestingly, within some other plant species, such as in soybean, there exists specialized Pht1 member(s) that is (are) mainly expressed in sink tissues, such as in flowers, implying possible functions of these genes in Pi import from source to sink [13, 71]. However, in our present study, we did not observe any member in the tomato Pht1 family that is dominantly expressed in the flower or fruit tissue. The relatively high transcription levels of LePT1 and LePT7 in these sink tissues compared to the other six paralogues led to the suggestion that the two members, especially the LePT1, might have been evolved to meet the requirement of transporting Pi from the source to the sink organs or cells. One of the potential evolutionary fates of gene duplications is considered to silence (nonfunctionalization) one of the duplicate copies [72]. The inactivation of LePT8 in all the tissues led to the suggestion that LePT8 might be on the way to become a pseudogene. The presence of several fragments that showed substantial homology to plant Pht1 genes but with incomplete coding regions in the tomato and potato genomes also well supports the theory that genome polyploidizations as well as the following gene loss (diploidization) are common characters of plant genomes [73, 74].

Conclusions

Taken together, this study provided the first comprehensive analysis of the chromosomal organization, phylogenetic evolution and tissue-specific expression patterns for each member of the Pht1 family in tomato. The results presented here could offer a useful basis for future research work on better understanding the mechanisms underlying the evolutionary regulation of Pht1 genes in response to Pi deficiency and AM symbiosis during tomato growth. The high conservation not only in the coding sequence, but also in the chromosomal distribution between the tomato and potato Pht1 orthologues could also lend strong evidence to support the further comparative genomics analysis across the whole Solanaceae family. However, we also realize that although Pht1 family have been commonly considered to be a high-affinity Pi transporter family, an increasing number of members of this family in other plant species have been demonstrated exhibiting a low or even dual affinity for Pi uptake in heterologous yeast or oocyte expression system [75, 76]. Therefore, more information, especially the transport kinetics, cellular distributions and physiological phenotypes of knockout/knockdown mutants is needed in the near future to determine more precise functional roles for each of the tomato Pht1 genes.

Methods

Plant material and growth conditions

Tomato (Solanum lycopersicum cv. Micro-Tom), was used in this study. The seeds were surface-sterilized, germinated and maintained in tissue culture for three weeks using 1/2 MS medium supplemented with 1.5% sucrose. The aseptic plantlets were then transferred to pot culture for either with high-/low-P treatment or inoculation with AM fungi.

In pot culture, two plantlets were transplanted to a 3 dm3 plastic pot filled with sterilized sand. A sand-based inoculum containing Glomus intraradices (Gi) was used for inoculation. Each plants was inoculated with 5 g inoculum or autoclaved inoculum around the roots. The plants were grown in a growth room with a 14-h light period (28-30°C) and a 10-h dark period (18-20°C). The irrigating solution contained the following nutrients: 1 mM NH4NO3, 2 mM KNO3, 0.5 mM Ca(NO3)2, 0.25 mM CaCl2, 0.5 mM MgSO4, 20 μM Fe-EDTA, 9 μM MnCl2, 46 μM H3BO3, 8 μM ZnSO4, 3 μM CuSO4, 0.03 μM (NH4)2MoO4, and 1 mM (high-P treatment) or 0.05 mM Pi (low-P treatment) NaH2PO4. The experiment comprised 4 replicates for each treatment. After inoculation for six weeks, the plants were either harvested for collecting the root, stem and young leaf samples or for continuing growing for the later collection of flower and fruit samples at young and ripe stages. The collected samples were immediately frozen in liquid nitrogen and stored at -80°C for subsequent RNA isolation.

Identification of Pht1 genes in tomato and potato genomes

Members of Pht1 gene family in the tomato genome were identified using the BLASTN and TBLASTN algorithm wrapped in the BLAST 2.2.27+ applications. To identify the potential Pht1 genes in tomato genome, coding sequence (cds) and deduced protein sequences of the Arabidopsis and rice Pht1 genes were queried respectively in the tomato genomic sequence database downloaded from the Solanaceae Genomics Network (http://www.sgn.cornell.edu). Sequences with a query over 50% and e-value less than -10 were taken as the Pht1 candidates. All the obtained sequences were submitted to NCBI (http://www.ncbi.nlm.nih.gov/) and Pfam database (http://www.sanger.ac.uk/Software/Pfam/search.shtml) for further confirmative analysis. For chromosomal localization analysis, the tomato Pht1 candidates were further used as queries for BLASTN searches against the SGN Tomato Whole Genome Scaffolds data (2.40) (http://solgenomics.net/organism/Solanum_lycopersicum/genome).

For identification of the potential homologues of tomato Pht1 genes from potato genome, the potato genome sequence data downloaded from the SGN database (http://solgenomics.net/organism/Solanum_tuberosum/genome) were also extensively searched using the tomato Pht1 genes as queries. The naming of the potato Pht1 genes were partially based on their phylogenetic relationships with the tomato homologues.

Phylogenetic analysis

The sequence data used in this study were collected using a query search in the NCBI database using the known Pht1 family gene sequences from Arabidopsis and rice. Multiple sequence alignments were performed using the program ClustalX (version 1.8) with default gap penalties. An un-rooted phylogenetic tree was generated using the deduced amino acid sequences of Pht1 genes by neighbor-joining algorithms wrapped in the MEGA 5.1 phylogeny program (http://www.megasoftware.net). Bootstrap analysis was carried out with 1,000 replicates.

RNA extraction and real-time RT-PCR analysis

Total RNA was isolated from 100 mg of various tissue samples, including roots, stems, young leaves and flowers, using the guanidine thiocyanate extraction method with Trizol reagent (Invitrogen) and from fruit samples using CTAB-sour phenol extraction method as described by Chang [77]. After extraction, the RNA samples were treated with DNase I (TaKaRa) to eliminate the trace contaminants of genomic DNA. For conducting reverse transcription (RT) PCR analysis, approximately two micrograms of total RNA from each sample was used to synthesize first-strand cDNA using a reverse transcription kit (TaKaRa), and the synthesized cDNAs were used as templates in the following PCR reactions.

Real-time RT-PCR analysis was performed to relatively quantify the expression levels of tomato Pht1 genes in different tissues or in roots and leaves in response to mycorrhizal colonization at high and low P status. The reaction was conducted on the Applied Biosystems (ABI) Plus Real-Time PCR System using the SYBER premix ExTaq kit (TaKaRa). Relative quantification of the transcripts for each tomato Pht1 gene was standardized to the expression level of the tomato constitutive Actin gene, calculated by the formula Y = 10-(ΔCt/3) × 100% (ΔCt is the differences of cycle threshold value between the target Pht1 gene and the control Actin products) [47, 78]. The specificity of primer sets designed for the qRT-PCR (Additional file 6) was confirmed by sequencing after the PCR reaction.

Histochemical GUS staining and detection of mycorrhizal fungal colonization

A 1250-bp LePT3 promoter fragment and a 471-bp LePT5 promoter fragment immediately upstream of the translation initiator ATG were amplified and cloned into binary vector pBI12, respectively, to replace the CaMV35S promoter in front of the β-glucuronidase (GUS) gene. The resulting two constructs were designated as pLePT3 -1250 and pLePT5 -471 , respectively, and introduced into Agrobacterium tumefaciens strain EHA105 for genetic transformation.

Histochemical GUS staining of the fresh transgenic roots was performed as described previously [44]. For visualization of fungal structures, the Magenta-GUS stained root segments were treated with 10% KOH solution heated to 90°C for 1 h, and then neutralized with 1% HCl (v/v) solution for 5 min. The root materials were then counter-stained with 0.3% trypan blue solution for 2 h at 90°C. The co-localization of Magenta-GUS and the trypan blue staining were indicated by purple color. The stained materials were rinsed in 50% glycerol and photographed by a stereomicroscope with a color CCD camera (Olympus).

Abbreviations

PT:

Phosphate transporter

EST:

Expressed sequence tag

AM:

Arbuscular mycorrhizal

AMF:

Arbuscular mycorrhizal fungi

Indels:

Insertions and deletions.

References

  1. Mengel K, Kirkby EA: Principles of plant nutrition. Dordrecht, the Netherlands: Kluwer Academic Publishers 2001.

    Book  Google Scholar 

  2. Abel S, Ticconi CA, Delatorre CA: Phosphate sensing in higher plants. Physiol Plant. 2002, 115: 1-8. 10.1034/j.1399-3054.2002.1150101.x.

    Article  CAS  PubMed  Google Scholar 

  3. Gu M, Xu K, Chen AQ, Zhu YY, Tang GL, Xu GH: Expression analysis suggests potential roles of microRNAs for phosphate and arbuscular mycorrhizal signaling in Solanum lycopersicum. Physiol Plant. 2010, 138: 226-237. 10.1111/j.1399-3054.2009.01320.x.

    Article  CAS  PubMed  Google Scholar 

  4. Schachtman DP, Reid RJ, Ayling SM: Phosphorus uptake by plants: from soil to cell. Plant Physiol. 1998, 116: 447-453. 10.1104/pp.116.2.447.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Smith SE, Jakobsen I, Grønlund M, Smith FA: Roles of arbuscular mycorrhizas in plant phosphorus nutrition: interactions between pathways of phosphorus uptake in arbuscular mycorrhizal roots have important implications for understanding and manipulating plant phosphorus acquisition. Plant Physiol. 2011, 156: 1050-1057. 10.1104/pp.111.174581.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Harrison MJ: Signaling in the arbuscular mycorrhizal symbiosis. Annu Rev Microbiol. 2005, 50: 361-389.

    Google Scholar 

  7. Raghothama KG, Karthikeyan AS: Phosphate acquisition. Plant Soil. 2005, 274: 37-49. 10.1007/s11104-004-2005-6.

    Article  CAS  Google Scholar 

  8. Javot H, Pumplin N, Harrison MJ: Phosphate in the arbuscular mycorrhizal symbiosis: transport properties and regulatory roles. Plant Cell Environ. 2007, 30: 310-322. 10.1111/j.1365-3040.2006.01617.x.

    Article  CAS  PubMed  Google Scholar 

  9. Muchhal US, Pardo JM, Raghothama KG: Phosphate transporters from the higher plant Arabidopsis thaliana. Ment Health Nurs. 1996, 93: 10519-10523.

    CAS  Google Scholar 

  10. Bun-Ya M, Nishimura M, Harashima M, Harashima S, Oshima Y: The PHO84 gene of Saccharomyces cerevisiaenencodes an inorganic phosphate transporter. Mol Cell Biol. 1991, 11: 3229-3238.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  11. Harrison MJ, van Buuren ML: A phosphate transporter from the mycorrhizal fungus Glomus versiforme. Nature. 1995, 378: 626-629. 10.1038/378626a0.

    Article  CAS  PubMed  Google Scholar 

  12. Mudge SR, Rae AL, Diatloff E, Smith FW: Expression analysis suggests novel roles for members of the Pht1 family of phosphate transporters in Arabidopsis. Plant J. 2002, 31: 341-353. 10.1046/j.1365-313X.2002.01356.x.

    Article  CAS  PubMed  Google Scholar 

  13. Qin L, Guo YX, Chen LY, Liang RK, Gu M, Xu GH, Zhao J, Walk T, Liao H: Functional characterization of 14 Pht1 family genes in yeast and their expressions in response to nutrient starvation in soybean. PLoS One. 2012, 7: e47726-10.1371/journal.pone.0047726.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Rae AL, Cybinski DH, Jarmey JM, Smith FW: Characterization of two phosphate transporters from barley; evidence for diverse function and kinetic properties among members of the Pht1 family. Plant Mol Biol. 2003, 53: 27-36.

    Article  CAS  PubMed  Google Scholar 

  15. Preuss CP, Huang CY, Tyerman SD: Proton-coupled high-affinity phosphate transport revealed from heterologous characterization in Xenopus of barley-root plasma membrane transporter, HvPHT1;1. Plant Cell Environ. 2011, 34: 681-689. 10.1111/j.1365-3040.2010.02272.x.

    Article  CAS  PubMed  Google Scholar 

  16. Liu F, Chang XJ, Ye Y, Xie WB, Wu P, Lian XM: Comprehensive sequence and whole-life-cycle expression profile analysis of the phosphate transporter gene family in rice. Mol Plant. 2011, 4: 1105-1122. 10.1093/mp/ssr058.

    Article  CAS  PubMed  Google Scholar 

  17. Chiou TJ, Liu H, Harrison MJ: The spatial expression patterns of a phosphate transporter (MtPT1) from Medicago truncatula indicate a role in phosphate transport at the root/soil interface. Plant J. 2001, 25: 281-293. 10.1046/j.1365-313x.2001.00963.x.

    Article  CAS  PubMed  Google Scholar 

  18. Nagy R, Vasconcelos MJ, Zhao S, McElver J, Bruce W, Amrhein N, Raghothama KG, Bucher M: Differential regulation of five Pht1 phosphate transporters from maize (Zea mays L.). Plant Biol (Stuttg). 2006, 8: 186-197. 10.1055/s-2005-873052.

    Article  CAS  Google Scholar 

  19. Gordon-Weeks R, Tong YP, Davies TG, Leggewie G: Restricted spatial expression of a high-affinity phosphate transporter in potato roots. J Cell Sci. 2003, 116: 3135-3144. 10.1242/jcs.00615.

    Article  CAS  PubMed  Google Scholar 

  20. Xiao K, Liu J, Dewbre G, Harrison M, Wang ZY: Isolation and characterization of root-specific phosphate transporter promoters from Medicago truncatula. Plant Biol (Stuttg). 2007, 8: 439-449.

    Article  Google Scholar 

  21. Misson J, Thibaud MC, Bechtold N, Raghothama KG, Nussaume L: Transcriptional regulation and functional properties of Arabidopsis Pht1;4, a high affinity transporter contributing greatly to phosphate uptake in phosphate deprived plants. Plant Mol Biol. 2004, 55: 727-741. 10.1007/s11103-004-1965-5.

    Article  CAS  PubMed  Google Scholar 

  22. Shin H, Shin HS, Dewbre GR, Harrison MJ: Phosphate transport in Arabidopsis: Pht1;1 and Pht1;4 play a major role in phosphate acquisition from both low- and high-phosphate environments. Plant J. 2004, 39: 629-642. 10.1111/j.1365-313X.2004.02161.x.

    Article  CAS  PubMed  Google Scholar 

  23. Jia HF, Ren HY, Gu M, Zhao JN, Sun SB, Zhang X, Chen JY, Wu P, Xu GH: The phosphate transporter gene OsPht1;8 is involved in phosphate homeostasis in rice. Plant Physiol. 2011, 156: 1164-1175. 10.1104/pp.111.175240.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Nagarajan VK, Jain A, Poling MD, Lewis AJ, Raghothama KG, Smith AP: Arabidopsis Pht1;5 mobilizes phosphate between source and sink organs and influences the interaction between phosphate homeostasis and ethylene signaling. Plant Physiol. 2011, 156: 1149-1163. 10.1104/pp.111.174805.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Sun SB, Gu M, Cao Y, Huang XP, Zhang X, Ai PH, Zhao JN, Fan XR, Xu GH: A constitutive expressed phosphate transporter, OsPht1;1, modulates phosphate uptake and translocation in phosphate-replete rice. Plant Physiol. 2012, 159: 1571-1581. 10.1104/pp.112.196345.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Rausch C, Daram P, Brunner S, Jansa J, Laloi M, Leggewie G, Amrhein N, Bucher M: A phosphate transporter expressed in arbuscule-containing cells in potato. Nature. 2001, 414: 462-466. 10.1038/35106601.

    Article  CAS  PubMed  Google Scholar 

  27. Harrison MJ, Dewbre GR, Liu J: A phosphate transporter from Medicago truncatula involved in the acquisition of phosphate released by arbuscular mycorrhizal fungi. Plant Cell. 2002, 14: 2413-2429. 10.1105/tpc.004861.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Paszkowski U, Kroken S, Roux C, Briggs SP: Rice phosphate transporters include an evolutionarily divergent gene specifically activated in arbuscular mycorrhizal symbiosis. Proc Natl Acad Sci U S A. 2002, 99: 3324-13329. 10.1073/pnas.052161299.

    Article  Google Scholar 

  29. Glassop D, Smith S, Smith F: Cereal phosphate transporters associated with the mycorrhizal pathway of phosphate uptake into roots. Planta. 2005, 222: 688-698. 10.1007/s00425-005-0015-0.

    Article  CAS  PubMed  Google Scholar 

  30. Nagy R, Karandashov V, Chague V, Kalinkevich K, Tamasloukht MB, Xu G, Jakobsen I, Levy AA, Amrhein N, Bucher M: The characterization of novel mycorrhiza-specific phosphate transporters from Lycopersicon esculentum and Solanum tuberosum uncovers functional redundancy in symbiotic phosphate transport in solanaceous species. Plant J. 2005, 42: 236-250. 10.1111/j.1365-313X.2005.02364.x.

    Article  CAS  PubMed  Google Scholar 

  31. Maeda D, Ashida K, Iguchi K, Chechetka SA, Hijikata A, Okusako Y, Deguchi Y, Izui K, Hata S: Knockdown of an arbuscular mycorrhiza-inducible phosphate transporter gene of Lotus japonicus suppresses mutualistic symbiosis. Plant Cell Physiol. 2006, 47: 807-817. 10.1093/pcp/pcj069.

    Article  CAS  PubMed  Google Scholar 

  32. Wegmüller S, Svistoonoff S, Reinhardt D, Stuurman J, Amrhein N, Bucher M: A transgenic dTph1 insertional mutagenesis system for forward genetics in mycorrhizal phosphate transport of Petunia. Plant J. 2008, 54: 1115-1127. 10.1111/j.1365-313X.2008.03474.x.

    Article  PubMed  Google Scholar 

  33. Guether M, Balestrini R, Hannah M, He J, Udvardi MK, Bonfante P: Genome-wide reprogramming of regulatory networks, transport, cell wall and membrane biogenesis during arbuscular mycorrhizal symbiosis in Lotus japonicus. New Phytol. 2009, 182: 200-212. 10.1111/j.1469-8137.2008.02725.x.

    Article  CAS  PubMed  Google Scholar 

  34. The Tomato Genome Consortium: The tomato genome sequence provides insights into fleshy fruit evolution. Nature. 2012, 485: 635-641. 10.1038/nature11119.

    Article  Google Scholar 

  35. Daram P, Brunner S, Persson BL, Amrhein N, Bucher M: Functional analysis and cell-specific expression of a phosphate transporter from tomato. Planta. 1998, 206: 225-233. 10.1007/s004250050394.

    Article  CAS  PubMed  Google Scholar 

  36. Catarecha P, Segura MD, Franco-Zorrilla JM, Garcı’a-Ponce B, Lanza M, Solano R, Paz-Ares J, Leyva A: A mutant of the Arabidopsis phosphate transporter Pht1;1 displays enhanced arsenic accumulation. Plant Cell. 2007, 19: 1123-1133. 10.1105/tpc.106.041871.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  37. Wu P, Shou HX, Xu GH, Lian XM: Improvement of phosphorus efficiency in rice on the basis of understanding phosphate signaling and homeostasis. Curr Opin Plant Biol. 2013, 16: 205-212. 10.1016/j.pbi.2013.03.002.

    Article  CAS  PubMed  Google Scholar 

  38. Remy E, Cabrito TR, Batista RA, Teixeira MC, Sá-Correia I, Duque P: The Pht1;9 and Pht1;8 transporters mediate inorganic phosphate acquisition by the Arabidopsis thaliana root during phosphorus starvation. New Phytol. 2012, 195: 356-371. 10.1111/j.1469-8137.2012.04167.x.

    Article  CAS  PubMed  Google Scholar 

  39. Liu F, Wang ZY, Ren HY, Shen CJ, Li Y, Ling HQ, Wu CY, Lian XM, Wu P: OsSPX1 suppresses the function of OsPHR2 in the regulation of expression of OsPT2 and phosphate homeostasis in shoots of rice. Plant J. 2010, 62: 508-517. 10.1111/j.1365-313X.2010.04170.x.

    Article  CAS  PubMed  Google Scholar 

  40. Xu GH, Chague V, Melamed BC, Kapulnik Y, Jain A, Raghothama KG, Levy AA, Silbere A: Functional characterization of LePT4: a phosphate transporter in tomato with mycorrhiza-enhanced expression. J Exp Bot. 2007, 58: 2450-2491.

    Google Scholar 

  41. The Potato Genome Sequencing Consortium: Genome sequence and analysis of the tuber crop potato. Nature. 2011, 475: 189-195. 10.1038/nature10158.

    Article  Google Scholar 

  42. Liu C, Muchhal US, Uthappa M, Kononowicz AK, Raghothama KG: Tomato phosphate transporter genes are differentially regulated in plant tissues by phosphorus. Plant Physiol. 1998, 116: 91-99. 10.1104/pp.116.1.91.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  43. Poirier Y, Bucher M: Phosphate transport and homeostasis in Arabidopsis. The Arabidopsis Book. Edited by: Somerville CR, Meyerowitz EM. 2002, Rockville: American Society of Plant Biologists

    Google Scholar 

  44. Karandashov V, Nagy R, Wegmuller S, Amrhein N, Bucher M: Evolutionary conservation of a phosphate transporter in the arbuscular mycorrhizal symbiosis. Proc Natl Acad Sci U S A. 2004, 101: 6285-6290. 10.1073/pnas.0306074101.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. Chen AQ, Gu M, Sun SB, Zhu LL, Hong S, Xu GH: Identification of two conserved cis-acting elements, MYCS and P1BS, involved in the regulation of mycorrhiza-activated phosphate transporters in eudicot species. New Phytol. 2011, 189: 1157-1169. 10.1111/j.1469-8137.2010.03556.x.

    Article  CAS  PubMed  Google Scholar 

  46. Leggewie G, Wilmitzer L, Riesmeier JW: Two cDNAs from potato are able to complement a phosphate uptake-deficient yeast mutant: identification of phosphate transporters from higher plants. Plant Cell. 1997, 9: 381-392.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  47. Chen AQ, He SE, Li FF, Li Z, Ding MQ, Liu QP, Rong JK: Analyses of the sucrose synthase gene family in cotton: structure, phylogeny and expression patterns. BMC Plant Biol. 2012, 12: 85-10.1186/1471-2229-12-85.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  48. Liu H, Trieu AT, Blaylock LA, Harrison MJ: Cloning and characterization of two phosphate transporters from Medicago truncatula roots: regulation in response to phosphate and to colonization by arbuscular mycorrhizal (AM) fungi. Mol Plant Microbe Interact. 1998, 116: 14-22.

    Article  Google Scholar 

  49. Liu J, Versaw WK, Pumplin N, Gomez SK, Blaylock LA, Harrison MJ: Closely related members of the Medicago truncatula pht1 phosphate transporter gene family encode phosphate transporters with distinct biochemical activities. J Biol Chem. 2008, 283: 24673-24681. 10.1074/jbc.M802695200.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Volpe V, Dell’aglio E, Giovannetti M, Ruberti C, Costa A, Genre A, Guether M, Bonfante P: An AM-induced, MYB-family gene of Lotus japonicus (LjMAMI) affects root growth in an AM-independent manner. Plant J. 2013, 73: 442-455. 10.1111/tpj.12045.

    Article  CAS  PubMed  Google Scholar 

  51. Huang CY, Shirley N, Genc Y, Shi B, Langridge P: Phosphate utilization efficiency correlates with expression of low-affinity phosphate transporters and noncoding RNA, IPS1, in barley. Plant Physiol. 2011, 156: 1217-1229. 10.1104/pp.111.178459.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  52. Nagy R, Drissner D, Amrhein N, Jakobsen I, Bucher M: Mycorrhizal phosphate uptake pathway in tomato is phosphorus-repressible and transcriptional regulated. New Phytol. 2009, 181: 950-959. 10.1111/j.1469-8137.2008.02721.x.

    Article  CAS  PubMed  Google Scholar 

  53. Rubio V, Linhares F, Solano R, Martín AC, Iglesias J, Leyva A, Paz-Ares J: A conserved MYB transcription factor involved in phosphate starvation signaling both in vascular plant and in unicellular algae. Genes Dev. 2001, 15: 2122-2133. 10.1101/gad.204401.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  54. Devaiah BN, Karthikeyan AS, Raghothama KG: WRKY75 transcription factor is a modulator of phosphate acquisition and root development in Arabidopsis. Plant Physiol. 2007, 143: 1789-1801. 10.1104/pp.106.093971.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  55. Xie XN, Huang W, Liu FC, Tang NW, Liu Y, Lin H, Zhao B: Functional analysis of the novel mycorrhiza-specific phosphate transporter AsPT1 and PHT1 family from Astragalus sinicus during the arbuscular mycorrhizal symbiosis. New Phytol. 2013, 198: 836-852. 10.1111/nph.12188.

    Article  CAS  PubMed  Google Scholar 

  56. Karandashov V, Bucher M: Symbiotic phosphate transport in arbuscular mycorrhizas. Trends Plant Sci. 2005, 10: 22-29.

    Article  CAS  PubMed  Google Scholar 

  57. Loth-Pereda V, Orsini E, Courty PE, Lota F, Kohler A, Diss L, Blaudez D, Chalot M, Nehls U, Bucher M, Martin F: Structure and expression profile of the phosphate Pht1 transporter gene family in mycorrhizal Populus trichocarpa. Plant Physiol. 2011, 156: 2141-2154. 10.1104/pp.111.180646.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  58. Balestrini R, Gómez-Ariza J, Lanfranco L, Bonfante P: Laser microdissection reveals that transcripts for five plant and one fungal phosphate transporter genes are contemporaneously present in arbusculated cells. Mol Plant Microbe Interact. 2007, 20: 1055-1062. 10.1094/MPMI-20-9-1055.

    Article  CAS  PubMed  Google Scholar 

  59. Liu QP: Identification of rice TUBBY-like genes and their evolution. FEBS J. 2008, 275: 163-171. 10.1111/j.1742-4658.2007.06186.x.

    Article  CAS  PubMed  Google Scholar 

  60. Flagel LE, Wendel JF: Gene duplication and evolutionary novelty in plants. New Phytol. 2009, 183: 557-564. 10.1111/j.1469-8137.2009.02923.x.

    Article  PubMed  Google Scholar 

  61. Chen AQ, Hu J, Sun SB, Xu GH: Conservation and divergence of both phosphate- and mycorrhiza-regulated physiological responses and expression patterns of phosphate transporters in solanaceous species. New Phytol. 2007, 173: 817-831. 10.1111/j.1469-8137.2006.01962.x.

    Article  CAS  PubMed  Google Scholar 

  62. Schünmann PHD, Richardson AE, Vickers CE, Delhaize E: Promoter analysis of the barley Pht1;1 phosphate transporter gene identifies regions controlling root expression and responsiveness to phosphate deprivation. Plant Physiol. 2004, 136: 4205-4214. 10.1104/pp.104.045823.

    Article  PubMed Central  PubMed  Google Scholar 

  63. Sobkowiak L, Bielewicz D, Malecka EM, Jakobsen I, Albrechtsen M, Szweykowska-Kulinska Z, Pacak A: The role of the P1BS element containing promoter-driven genes in Pi transport and homeostasis in plants. Front Plant Sci. 2012, 3: 58-

    PubMed Central  CAS  PubMed  Google Scholar 

  64. Jain A, Nagarajan VK, Raghothama KG: Transcriptional regulation of phosphate acquisition by higher plants. Cell Mol Life Sci. 2012, 69: 3207-3224. 10.1007/s00018-012-1090-6.

    Article  CAS  PubMed  Google Scholar 

  65. Smith SE, Smith FA, Jakobsen I: Mycorrhizal fungi can dominate phosphate supply to plants irrespective of growth responses. Plant Physiol. 2003, 133: 16-20. 10.1104/pp.103.024380.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  66. Smith SE, Smith FA, Jakobsen I: Functional diversity in arbuscular mycorrhizal (AM) symbioses: the contribution of the mycorrhizal P uptake pathway is not correlated with mycorrhizal responses in growth or total P uptake. New Phytol. 2004, 162: 511-524. 10.1111/j.1469-8137.2004.01039.x.

    Article  Google Scholar 

  67. Smith SE, Smith FA: Roles of arbuscular mycorrhizas in plant nutrition and growth: new paradigms from cellular to ecosystem scales. Annu Rev Plant Biol. 2011, 62: 227-250. 10.1146/annurev-arplant-042110-103846.

    Article  CAS  PubMed  Google Scholar 

  68. Javot H, Penmetsa RV, Terzaghi N, Cook DR, Harrison MJ: A phosphate transporter indispensable for the arbuscular mycorrhizal Medicago truncatula symbiosis. Proc Natl Acad Sci U S A. 2007, 104: 1720-1725. 10.1073/pnas.0608136104.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  69. Yang SY, Paszkowski U: Phosphate import at the arbuscule: just a nutrient?. Mol Plant Microbe Interact. 2011, 24: 1296-1299. 10.1094/MPMI-06-11-0151.

    Article  CAS  PubMed  Google Scholar 

  70. Yang SY, Grønlund M, Jakobsen I, Grotemeyer MS, Rentsch D, Miyao A, Hirochika H, Kumar CS, Sundaresan V, Salamin N, Catausan S, Mattes N, Heuer S, Paszkowski U: Nonredundant regulation of rice arbuscular mycorrhizal symbiosis by two members of the phosphate transporter1 gene family. Plant Cell. 2012, 24: 4236-4251. 10.1105/tpc.112.104901.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  71. Fan C, Wang X, Hu RB, Wang YH, Xiao CW, Jiang Y, Zhang XM, Zheng CY, Fu YF: The pattern of Phosphate transporter 1 genes evolutionary divergence in Glycine max L. BMC Plant Biol. 2013, 13: 48-10.1186/1471-2229-13-48.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  72. Lynch M, Conery JS: The evutionary fate and consequences of duplicate genes. Science. 2000, 290: 1151-1155. 10.1126/science.290.5494.1151.

    Article  CAS  PubMed  Google Scholar 

  73. Guyot R, Keller B: Ancestral genome duplication in rice. Genome. 2004, 47: 610-614. 10.1139/g04-016.

    Article  CAS  PubMed  Google Scholar 

  74. Xie Z, Li X, Glover BJ, Bai S, Rao GY, Luo J, Yang J: Duplication and functional diversification of HAP3 genes leading to the origin of the seed-developmental regulatory gene, LEAFY COTYLEDON1 (LEC1), in nonseed plant genomes. Mol Biol Evol. 2008, 25: 1581-1592. 10.1093/molbev/msn105.

    Article  CAS  PubMed  Google Scholar 

  75. Preuss CP, Huang CY, Gilliham M, Tyerman SD: Channel-Like characteristics of the low-affinity barley phosphate transporter Pht1;6 when expressed in xenopus oocytes. Plant Physiol. 2010, 152: 1431-1441. 10.1104/pp.109.152009.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  76. Ai PH, Sun SB, Zhao JN: Two rice phosphate transporters, OsPht1;2 and OsPht1;6, have different functions and kinetic properties in uptakeand translocation. Plant J. 2009, 57: 798-809. 10.1111/j.1365-313X.2008.03726.x.

    Article  CAS  PubMed  Google Scholar 

  77. Chang SJ, Puryear J, Cairney J: A simple and efficient method for isolating RNA from pine trees. Plant Mol Biol. 1993, 11: 113-116. 10.1007/BF02670468.

    Article  CAS  Google Scholar 

  78. Li XB, Fan XP, Wang XL, Cai L, Yang WC: The cotton ACTIN1 gene is functionally expressed in fibers and participates in fiber elongation. Plant Cell. 2005, 17: 859-875. 10.1105/tpc.104.029629.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported by the National Natural Science Foundation of China (31372121, 31272225), the Basic Research Program of Jiangsu province in China (BK2012765), and A Foundation for the Author of National Excellent Doctoral Dissertation of PR China (FANEDD, 201264).

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Aiqun Chen or Xiao Chen.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

AQC and GHX contributed to the experimental design and manuscript drafting. MG contributed to the manuscript editing. XC, HMW and DHL performed the RNA extraction, primer design, RT-PCR validation. AQC and HYQ performed the bioinformatics analysis. All authors have read and approved the final manuscript.

Aiqun Chen, Xiao Chen contributed equally to this work.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

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 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.

Reprints and permissions

About this article

Cite this article

Chen, A., Chen, X., Wang, H. et al. Genome-wide investigation and expression analysis suggest diverse roles and genetic redundancy of Pht1 family genes in response to Pi deficiency in tomato. BMC Plant Biol 14, 61 (2014). https://doi.org/10.1186/1471-2229-14-61

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2229-14-61

Keywords