Original Research ARTICLE
Using the yeast three-hybrid system to identify proteins that interact with a phloem-mobile mRNA
- 1 Plant Biology Major, Iowa State University, Ames, IA, USA
- 2 DuPont Company, Crop Protection, Stine-Haskell Research Center, Newark, DE, USA
- 3 Department of Biology and Physical Sciences, Chowan University, Murfreesboro, NC, USA
Heterografting and RNA transport experiments have demonstrated the long-distance mobility of StBEL5 RNA, its role in controlling tuber formation, and the function of the 503-nt 3′ untranslated region (UTR) of the RNA in mediating transport. Because the 3′ UTR of StBEL5 is a key element in regulating several aspects of RNA metabolism, a potato leaf cDNA library was screened using the 3′ UTR of StBEL5 as bait in the yeast three-hybrid (Y3H) system to identify putative partner RNA-binding proteins (RBPs). From this screen, 116 positive cDNA clones were isolated based on nutrient selection, HIS3 activation, and lacZ induction and were sequenced and classified. Thirty-five proteins that were predicted to function in either RNA- or DNA-binding were selected from this pool. Seven were monitored for their expression profiles and further evaluated for their capacity to bind to the 3′ UTR of StBEL5 using β-galactosidase assays in the Y3H system and RNA gel-shift assays. Among the final selections were two RBPs, a zinc finger protein, and one protein, StLSH10, from a family involved in light signaling. In this study, the Y3H system is presented as a valuable tool to screen and verify interactions between target RNAs and putative RBPs. These results can shed light on the dynamics and composition of plant RNA-protein complexes that function to regulate RNA metabolism.
For plant development, phloem plays important roles in not only transporting nutrients, but also as a conduit for moving signal RNAs and proteins. Full-length, phloem-mobile mRNAs function to integrate environmental cues for plant development via this long-distance signaling pathway (Haywood et al., 2005; Banerjee et al., 2006). There are three types of mobile RNAs in plants; (1) pathogenic viral and viroid RNAs, (2) small RNAs including siRNAs and microRNAs, and (3) full-length cellular RNA transcripts (Kehr and Buhtz, 2008). Although many full-length transcripts have been identified in the phloem, only a few of these transcripts have been confirmed to be mobile through the phloem translocation stream. The best examples of mobile RNAs are StBEL5 (Banerjee et al., 2006), CmGAI (Haywood et al., 2005), and the Arabidopsis FLOWERING LOCUS T (Li et al., 2011; Lu et al., 2012). StBEL5 is a transcription factor that works in tandem with Knotted1-types to regulate plant growth (Chen et al., 2003, 2004). RNA detection methods and heterografting experiments demonstrated that StBEL5 transcripts are present in phloem cells and move across a graft union to localize in stolon tips, the site of tuber induction (Banerjee et al., 2006). This movement of RNA originates in leaf veins and petioles and is induced by a short-day photoperiod, regulated by the untranslated regions, and correlated with enhanced tuber production (Banerjee et al., 2006, 2009). Long-distance movement of the RNA of GA INSENSITIVE (GAI) has also been clearly established in both cucumber and pumpkin (Haywood et al., 2005; Ham et al., 2009). Recent results suggest that in addition to FT protein, FT RNA may also be moving to shoot apices to contribute to systemic floral signaling (Li et al., 2011; Lu et al., 2012).
In general, RNA molecules are associated with RNA-binding proteins (RBPs) in the cell, and a number of RNA-protein interactions have been established. RBPs function in splicing, nuclear export, RNA transport and localization, translation, and stability (Dreyfuss et al., 2002; Fedoroff, 2002). RBPs are involved in coordinating gene expression and also influence the localization of protein synthesis (Lunde et al., 2007). For example, a polypyrimidine-tract binding protein (PTB), designated as CmRBP50, was reported as the core protein of a phloem-mobile ribonucleoprotein complex consisting of six RNAs, including CmGAI RNA, and 16 proteins in pumpkin phloem sap (Ham et al., 2009). Commonly, it is the UTRs that function via protein interactions in facilitating the cellular localization of a transcript (Jansen, 2001), in mediating its stability (Lee and Jeong, 2006), or in regulating the efficiency of translation (Barreau et al., 2006). Binding motifs have been identified in the RNAs of animals that function in recognizing RBPs (for review, see Jansen, 2001). These motifs are most predominant in the 3′ UTR (Saunders and Cohen, 1999; Corral-Debrinski et al., 2000; Thio et al., 2000). There are numerous examples demonstrating the importance of the 3′ UTRs in recognizing RBPs that regulate metabolism and movement (Ferrandon et al., 1994; Padmanabhan and Richter, 2006; Irion and St. Johnston, 2007). As a prime example in plants, the 3′ UTR of StBEL5 plays a significant role in mediating its long-distance transport, controlling translation, and regulating stability (Banerjee et al., 2006, 2009), suggesting the presence of cis-elements in this UTR that are recognized by RNA-binding partners.
Although there are several useful biochemical approaches to analyze RNA-protein interactions, the yeast three-hybrid (Y3H) system (Sengupta et al., 1996; Hook et al., 2005) represents a simple but powerful tool for searching a large collection of cDNAs to identify proteins that bind a specific RNA of interest (Cassiday and Maher III, 2003; Gonsalvez et al., 2003; Maniataki et al., 2003; Moore et al., 2003; Campalans et al., 2004; Hwang et al., 2005). Not only does it allow the identification of RNA-protein binding partners but also the dissection of higher-order RNA-protein complexes (Bernstein et al., 2002). Using the 503-nt 3′ UTR of StBEL5 as bait, the Y3H system was used with a potato leaf cDNA library for screening binding partners that may be involved in the metabolism of the full-length, mobile RNA, StBEL5. Initially, more than 100 cDNA clones were isolated from the screening based on nutrient selection and HIS3 and β-galactosidase activation. Seven proteins were selected based on their putative RNA-binding properties for further analyses and RNA gel-shift assays. These results clearly demonstrate the utility of the Y3H system in identifying candidate RBPs.
Materials and Methods
Constructs for the Y3H System
The DNA fragments encoding full-length and truncated forms (D1, T2, and UA baits) of the 3′ UTR of StBEL5 were amplified with gene-specific primer sets (Table A5 in Appendix) and cloned into pIIIA/MS2-1. The plasmids containing the hybrid RNA fused to the full-length 3′ UTR and the truncated forms were transformed into the YBZ-1 yeast strain. For screening, an amplified leaf cDNA library from potato (Solanum tuberosum cv Désirée) similar in design to the stolon library described by Chen et al. (2003) was used. This library was directionally cloned into pAD-GAL4-2.1 (Stratagene, La Jolla, CA, USA) and was a generous gift from Salomé Prat, Madrid, Spain. The YBZ-1 strain and the pIIIA/MS2-1 plasmid were graciously provided Dr. Marvin Wickens, University of Wisconsin, Madison.
Screening of RNA-binding Proteins
The YBZ-1 yeast strain containing the StBEL5 full-length 3′ UTR hybrid RNA was transformed with 60 μg of the potato cDNA library. The entire transformation mixture (about 10 mL) was spread onto plates containing SD/-his/-leu/-ura and 1.0 mM 3-aminotrizole (3-AT), a competitive inhibitor of the HIS3 gene product. To calculate the transformation efficiency, the transformation mixture was serially diluted (10−1, 10−2, 10−3) and grown on small SD/-leu/-ura plates. After 1st and 2nd rounds of screening on different concentration (0, 1, 5, 10, or 50 mM) of 3-AT-containing plates, positive colonies were applied to β-galactosidase assays in order to measure the induction of the lacZ reporter gene using a Yeast β-galactosidase Assay Kit (Pierce Biotechnology), according to the manufacturer’s protocol. From the assays, 116 colonies were selected. The 3′ UTR of StBEL5 with StPTB6-pAD (Mahajan et al., 2012) and an empty pAD were used as positive and negative controls, respectively.
Plasmid Rescue, Identification and Categorization of the Screened Clones
In order to identify the positive clones from the screening, yeast plasmid rescue was performed with E.Z.N.A.® Yeast Plasmid Kit (OMEGA bio-tek) with slight modifications. Positive yeast colonies were picked from the plate and inoculated in 3.0 ml of SD/-leu. Overnight grown cells were pelleted and incubated at 30°C for at least 30 min after resuspension in 480 μl Buffer SE/β-mercaptoethanol and 40 μl lyticase solutions. After incubation, yeast plasmid DNA was isolated by following the kit protocols. The rescued yeast plasmids were transformed into E. coli HB101 competent cells, and the plasmids were isolated and sequenced using pGAD-specific primers (Table A5 in Appendix) at the DNA Facility, Iowa State University. For putative identities of the clones, sequences of the cDNAs were analyzed at the Dana–Farber Cancer Institute (DFCI) Gene Index1 potato database. Translated protein sequences were obtained from Translate2 and ORF Finder3 and analyzed using BLAST on the TAIR database4. Functional categorization of selected proteins was performed using the MIPS Arabidopsis thaliana database5. The domains of the B5RBPs (Figure 3) were analyzed using SMART6 and NCBI’s Conserved Domain Search. For Figure 4A, amino acid sequences of Arabidopsis and potato LSH proteins were organized into a phylogenetic tree with the MEGA 4.0.2 package and the neighbor-joining program. The numbers listed at the branching points are boot-strapping values that indicate the level of significance (percentage) for the separation of two branches.
RNA Gel-Shift Assays
The PCR-amplified fragments with gene-specific primer sets (Table A5 in Appendix) were cloned into the pET-28a (+) plasmids after proper enzyme digestion to produce histidine tag (His)-fusion recombinant proteins (Figure A1 in Appendix). The constructs were transformed into E. coli BL21-Codon (DE3) cells (Stratagene). The recombinant proteins were induced with 0.4 mM IPTG, and purified using HisPur Cobalt Purification Kit (Pierce Biotechnology). For in vitro transcription to generate RNA probes, T3 promoter-containing sense primers were created by adding T3 sequences on the 5′ end of the sense primers of the target sequences (Table A5 in Appendix), and used for PCR amplification using Platinum Taq DNA Polymerase High Fidelity (Invitrogen). The gel-purified PCR product was transcribed using MEGAscript T3 (Ambion) incorporating biotin (biotin-11-UTP, Perkin Elmer) as described by the manufacturer’s manual. The biotin-labeled probe RNA was purified by gel purification using ZymocleanTM Gel RNA Recovery kit (ZymoResearch). Five femtomoles of biotin-labeled RNA probes were incubated with indicated amounts of purified recombinant proteins in the binding buffer provided by the Light Shift Chemiluminescent RNA EMSA kit (Pierce Biotechnology) on ice for 45 min. The RNA-protein complexes were separated in 2.5% agarose (for the full-length probe) or 5% polyacrylamide gel (for the IRE probe) and transferred onto BrightStar-Plus (Ambion) nylon membranes. The signal was detected using the EMSA kit according to the manufacturer’s manual.
Screening for RNA-binding Proteins Using the Yeast Three-Hybrid System
As eloquently explained by Bernstein et al. (2002), the Y3H system is based on two expression vectors, one for the RNA bait and the other for the protein target, and three-hybrid components. When bait RNA interacts with the target protein, the reporter genes, HIS3 and lacZ, are activated and can be readily detected by simple biochemical assays (Hook et al., 2005; Seay et al., 2006).
To isolate interacting proteins, the DNA fragment encoding 3′ UTR of StBEL5 as a bait was cloned into the MS2 portion of pIIIA/MS2-1 vector. The resulting plasmids carrying hybrid MS2-3′ UTR of StBEL5 were transformed into a yeast strain, YBZ-1, and the potato leaf cDNA library was sequentially transformed using conventional protocols with slight modifications (Bernstein et al., 2002; Seay et al., 2006). We analyzed approximately 6.5 × 105 yeast colonies, and the resulting transformed colonies were screened on SD/-his/-leu/-ura plates containing 1.0 mM 3-AT (Figure 1). From the first round, 448 colonies were selected as primary positives. Those selected positive colonies were replicated on SD/-his/-leu/-ura plates containing 1.0 mM 3-AT again to remove potential false positives. From these screenings, 281 colonies were chosen for further screening by using the two reporter genes, HIS3 and lacZ. SD/-his/-leu/-ura plates containing a series of 3-AT concentration (0, 1, 3, 5, 10, and 50 mM) were used for testing HIS3 expression. The 281 colonies were streaked on these plates, and 194 colonies were grown on SD/-his/-leu/-ura plates containing at least 5.0 mM 3-AT. Finally, 116 colonies were selected based on β-galactosidase and HIS3 activation and were sequenced and analyzed (Table A1 in Appendix). The overall strategy of the screening is summarized in Figure 1. As a comparison, using the yeast strain YBZ-1, 49 HIS3+ colonies were selected out of approximately 8.0 × 105 transformants according to Hook et al. (2005). The interactions of the 3′ UTR of StBEL5 with StPTB6 and the 3′ UTR of StBEL5 with empty pGAD vector were used as positive and negative controls for RNA-protein interaction, respectively. PTB proteins are multifunctional proteins that bind numerous mRNAs and are involved in a wide range of RNA metabolism, such as RNA stability, splicing, translational repression, and long-distance transport. There are six PTBs in the potato genome, and one, designated StPTB6, binds to untranslated regions of phloem-mobile mRNAs of potato (Mahajan et al., 2012).
Figure 1. Schematic diagram of the Y3H system for screening interacting partner proteins from a leaf cDNA library using the 503-nt 3′ UTR of StBEL5 as bait.
Characterization of Selected cDNA Clones
For identification of the screened colonies, their sequences were analyzed and BLAST was performed using the Arabidopsis and potato databases, and a total of 89 clones (76.7%) exhibited significant matches to previously characterized or known genes (Table A2 in Appendix). Thirteen clones were identified as redundant (30.2% redundancy), and therefore, a total of 94 unique singletons were isolated from the Y3H screening (listed in Tables A1– A3 in Appendix). These 13 included LSH3 (Light-dependent Short Hypocotyls3), LSH10 (Light-dependent Short Hypocotyls10), C3H zinc finger transcription factor, sucrose synthase4, a Transducin/WD40 repeat-like superfamily protein, LTP12 (Lipid Transfer Protein 12), ELI3-1 (elicitor-activated gene 3-1), X-ray induced transcript 1, glutamate-1-semialdehyde-2, 1-aminomutase, and some ribosomal proteins and unknown proteins (Table A3 in Appendix). Interestingly, LSH10, a close sequence match to AtLSH10 (AT2G42610), was identified from 10 clones along with another LSH member, LSH3, which was isolated twice. Twenty-seven clones (23.3%) were categorized as undefined clones, i.e., “unknown” or “no hit” from the database search (Table A2 in Appendix). Functional classification revealed a total of eight cDNAs that encoded proteins with DNA/RNA-binding properties. Nineteen clones encoded proteins that are components of the machinery for protein synthesis (16.4%) at the initiation and/or elongation of translation, such as eIF5A, transducin/WD40 repeat-like superfamily, heat shock protein 70, and an alpha-tubulin protein (Doroshenk et al., 2009; Lin et al., 2009; Tables A1 and A4 in Appendix). Overall, 35 cDNAs (30.1% of total screened clones) encoding proteins with RNA-binding properties were identified (Table A4 in Appendix).
Candidate Protein Partners for StBEL5 RNA
Based on activities of marker genes (HIS3 and lacZ), and their putative RNA-binding function, seven cDNA clones were selected for verification as protein partners of StBEL5, and designated as StBEL5 RNA-binding protein (B5RBP) one to seven (Figure 2A). Related proteins of three of the B5RBPs, B5RBP1, -4, and -6, were previously identified as RBPs in other species (Xu et al., 2004; Doroshenk et al., 2009; Ling et al., 2011), and three others, B5RBP2, -5, and -7, contain conserved RNA recognition motifs (RRMs; Figure 3). The most frequently identified clone from the Y3H screening, StLSH10 (B5RBP3), was also included (Table A3 in Appendix). Each of these proteins induced β-galactosidase activity in an interaction with the bait RNA to levels much higher than the negative controls and, in some cases, even higher than the positive control (Figures 2B,C).
Figure 2. List of seven B5RBPs (A) and a quantitative analyses of their interactions with the 503-nt 3 ′ UTR of StBEL5 using a β-galactosidase assay (B). The 3′ UTR of StBEL5 with StPTB6-pAD was used as a positive control (B), and an empty pAD (B) or the pIIIA/MS2-1 bait vectors (C) were used as negative controls. The numbers above each bar (B) represent the means for β-galactosidase activity in triplicate, and standard errors are shown for each mean (B,C). The RNA produced by the pIIIA/MS2-1 vector in yeast, including the 60-nt MS2 sites, is approximately 272 nt in length (Bernstein et al., 2002).
B5RBP1 encodes a sucrose synthase (SUS4), containing a sucrose synthase motif for the sucrose metabolic pathway and a glycosyltransferase motif for biosynthetic processes (Figure 3). There were three reasons to include SUS4. First, it was identified as a cytoskeleton-associated RBP from developing rice seeds (Doroshenk et al., 2009). Second, in Arabidopsis, SUS4 was detected in the companion cells of the phloem (Fallahi et al., 2008) and third, it plays an important role in starch metabolism in potato tubers (Fu et al., 1995; Zrenner et al., 1995; Bieniawska et al., 2007). With regard to the role of SUS4 in potato tuber development, high-level increases in SUS4 promoter activity were observed during early tuber formation (Fu et al., 1995) and transgenic lines overexpressing SUS4 produced enhanced tuber yields (Baroja-Fernandez et al., 2009). SUS4 could be an example of a multifunctional RBP directly involved in potato tuber development.
Figure 3. Protein structure and prediction of conserved domains of the seven B5RBPs. The deduced amino acid sequences of the B5RBPs were analyzed using SMART (http://smart.embl-heidelberg.de/) and NCBI’s Conserved Domain Search. RRM, RNA recognition motif; DUF, domain of unknown function; GRM, glycine-rich motif; KOW, from Kyrpides et al. (1996); CC, coiled-coil motif; ZF-CCCH, zinc finger, CCCH-type.
B5RBP2 encodes a glycine-rich RNA-binding protein that contains a RRM for RNA-binding and a glycine-rich motif (GRM, Figure 3). Plant proteins that contain a GRM are grouped into five classes based on structure. B5RBP2 is considered a class IV member because it contains a RRM (Mangeon et al., 2010). Class IV GRM-proteins are subdivided based on their links to osmotic stress, cold stress, flower timing, development, and responsiveness to abscisic acid. Interestingly, B5RBP2 is orthologous to AtGRP7 (AT2G21660) a protein related to a RBP found in pumpkin phloem sap (Lin et al., 2009). B5RBP2 may function in potato phloem sap by interacting with StBEL5 RNA to facilitate long-distance movement. AtGRP7 is also involved in the regulation of alternative splicing, ribosome function, and RNA metabolism (Wachter et al., 2012).
B5RBP3 is alight-dependent short hypocotyl (LSH10, AT2G42610) protein with a Domain of Unknown Function (DUF640, Figure 3). In Arabidopsis, there are 10 LSH genes and in potato, 15 (Figure 4A). The function of most these proteins, ranging in size from 164 to 219 aa, are unknown except for AtLSH1, -3, and -4. AtLSH1 is a nuclear protein in Arabidopsis with a nuclear localization signal (NLS) in the C-terminal region. It is involved in light regulation of seedling development (Zhao et al., 2004). Both AtLSH3 and -4 appear to have a role in Arabidopsis shoot and floral organ differentiation since constitutive expression of these genes resulted in abnormal development (Takeda et al., 2011). As described earlier, LSH10 was the most frequently selected cDNA from the screen (Table A3 in Appendix). β-galactosidase activity of the B5RBP3/StBEL5 interaction was several-fold greater than the other B5RBPs (Figure 2B). All of the potato LSH proteins exhibit a highly conserved internal region representative of a domain from the DUF640 superfamily flanked by sequence of considerable variance at both the amino- and carboxy-termini (Figure 4B).
Figure 4. Phylogenetic analysis (A) and alignment of the aa sequence (B) of the family of LSH proteins of potato (St). Included in the dendrogram (A) are the 10 LSH proteins of Arabidopsis (At). Arrows indicate the two StLSHs, LSH3 and -10, identified from the current yeast three-hybrid screening. The nuclear localization signal (NLS) is boxed (B). Asterisks below the aligned amino acids (B) indicate conserved residue identity.
B5RBP4 encodes a eukaryotic initiation factor 5A (eIF5A, AT1G13950) containing both KOW (acronym of the authors surname, Kyrpides et al., 1996; Figure 3) and eIF5A motifs for ribosome binding, RNA-binding, and translation activity (Figure 3). Studies in bacteria suggest that the KOW motif mediates RNA and protein interactions (Steiner et al., 2002). For simplicity, the eIF5A motif is referred to as a RNA-binding site (Figure 3) since the motif is characterized as a S1-like RNA-binding domain (Peat et al., 1998). eIF5A is a multifunctional protein involved in RNA-binding, processing, turnover, and transport from the nucleus to cytoplasm and in transcription and translation (Burd and Dreyfuss, 1994; Cusack, 1999; Xu and Chen, 2001). Recently, eIF5A in yeast was shown to stimulate protein synthesis but was not required for the process (Henderson and Hershey, 2011). In pumpkin phloem sap, CmeIF5A was detected as a component of the RBP50-based ribonucleoprotein complex (Ham et al., 2009). Further characterization of CmelF5A revealed that hypusination (lysine residue modification) was necessary for RNA-binding and protein interaction and that both hypusinated and non-hypusinated CmelF5A existed in the phloem (Ma et al., 2010). B5RBP4 may function through a similar mechanism in potato phloem sap to mediate the formation of ribonucleoprotein complexes.
B5RBP5 encodes a RNA-binding (RRM/RNA-Binding Domain/Ribonucleoprotein → RRM/RBD/RNP) protein family member containing a coiled-coil motif (CC) and a RRM (Figure 3) that is orthologous to At5G65260, an Arabidopsis RNA-binding protein. At5G65260 is designated as a polyadenylation factor that can bind to the poly (A) tail and control its length (Hunt et al., 2008). The Arabidopsis transcription factor Long Hypocotyl5 (HY5) that is involved in photomorphogenesis was shown to mediate the expression of At5G65260 (Lee et al., 2007). At5G65260 expression was down-regulated in a loss-of-function HY5 mutant. It is conceivable that B5RBP5 levels in potato may also be regulated by light.
B5RBP6 encodes a zinc finger (CCCH-type) family protein containing two ankyrin repeats for protein-protein interactions, and two zinc finger-C3H1 domains (ZF-CCCH) for zinc ion binding, and nucleic acid binding (Figure 3). Unlike other zinc finger proteins that generally function as DNA-binding proteins (Laity et al., 2001), CCCH zinc finger proteins bind to AU-rich elements of RNAs (Brown, 2005). A study in mouse revealed a role for a CCCH zinc finger protein (tristetraprolin) in mRNA decay (Lai et al., 2006). In Trypanosoma brucei, the causative agent of sleeping sickness, the CCCH zinc finger protein (Tb2C3H20) functions in mRNA stability (Ling et al., 2011). In Arabidopsis, two CCCH-type zinc finger genes, designated AtSZF1, and AtSZF2, were salt-inducible and mediated responsiveness to salt (Sun et al., 2007). From the current screen, two CCCH-type zinc finger proteins were identified (Table A4 in Appendix) that shared sequence similarity with AT2G40140 (AtSZF2) and AT5G58620 from Arabidopsis.
B5RBP7 encodes another RNA-binding (RRM/RBD/RNP motifs) family protein that shares sequence similarity with Arabidopsis AtRNP1 (AT4G14300) containing two RRM domains (Figure 3). AtRNP1 is a target of Arabidopsis transportin 1 (AtTRN1) that is an ortholog of the human nuclear import receptor transportin1 protein (Ziemienowicz et al., 2003). AtRNP1 may function as a shuttle protein moving RNAs between the nucleus and cytoplasm. Interestingly, AtGRP7 also interacted with AtTRN1. As discussed above, B5RBP2 is orthologous to AtGRP7 suggesting that B5RBP7 and B5RBP2 may associate with StBEL5 RNA as a tandem complex.
To assess transcript levels for select RBPs, expression values were obtained from the publicly available RNA-seq database from the RH genotype of the Potato Genome Sequencing Project (Xu et al., 2011). Abundance levels of StBEL5 and StHSP70 have been included as references. The potato HSP70 protein was selected during the screen (Table A4 in Appendix) and a HSP70-type was previously identified as a member of a phloem-mobile RNP complex in pumpkin (Ham et al., 2009). Relatively high and consistent levels of transcripts across all organs were observed for B5RBP2, -4, HSP70, and both RBPs, B5RBP5, and -7 (Table 1). The very high levels of eIF5A (B5RBP4) likely reflect its general, multifunctional role in several aspects of RNA metabolism (Zanelli and Valentini, 2007). B5RBP1 (sucrose synthase) scored abundant RNA levels in both stolons and young tubers indicative of its role in starch metabolism during tuber formation. The zinc finger CCCH protein (B5RBP6) was most abundant in petioles with a value of 334 FPKMs (fragments per kb per million mapped reads). A mobile RNA like StBEL5 is very abundant in petioles, an observation that is consistent with both its transcriptional source and its capacity to move long-distances through the phloem (Banerjee et al., 2006). Petioles serve two main functions: to provide support for the leaf lamina and to act as a protective sleeve for phloem cells that move sugar and signaling molecules (like RNA) from source leaves to sinks. RBPs are commonly detected in companion cells and sieve elements of leaf veins in position to chaperone mobile RNAs (Ham et al., 2009).
Table 1. Expression profile of select B5RBPs mined using the RNA-seq data from the publically available Potato Genome Database (Xu et al., 2011).
Despite the observation that B5RBP3 (LSH10) appeared 10 times from the Y3H screening (Table A3 in Appendix), its RNA levels were remarkably low in leaf RNA of the RH genotype (Table 1). Similar results for LSH10 RNA abundance levels were observed in RNA-seq data from the DM genotype (Xu et al., 2011). Among the proteins selected, this putative RBP exhibited the strongest induction of β-galactosidase activity (Figure 2). Whereas the second LSH protein selected in this screen, StLSH3, exhibited very low transcript values across all organs, transcript abundance values for B5RBP3 (Table 1) were extremely high in petioles (77 FPKMs), stolons (149 FPKMs), and young tubers (590 FPKMs). This latter transcript value was the second most abundant of any of the RNAs scored in this experiment.
Verification of RNA-protein Interaction of the B5RBPs
To validate the direct interaction of select proteins with the 3′ UTR of StBEL5, RNA gel-shift assays were performed with select B5RBPs. These assays were performed using biotin-incorporated RNA probes using the full-length 3′ UTR of StBEL5 and the purified recombinant B5RBP2, -3, and -5 proteins (Figure 5A). For showing specificity of the interaction, IRE RNA, which binds specifically with the iron responsive protein in the cell under iron-starved conditions, was used as a negative control (Figure 5B). Shifted bands were observed for all three interactions in a range of 30–250 nM of protein. B5RBP3 affected a shift with protein amounts as low as 30 nM, whereas shifted bands were clearly observed with the other two B5RBPs in the reactions containing 90 nM of protein. With comparable amounts of protein, no gel-shift was observed for the negative control, the iron responsive element (IRE, Figure 5B).
Figure 5. Mobility shift assays for in the vitro interaction of the 3′ UTR of StBEL5 with select B5RBPs. (A) Full-length 3′ UTR of StBEL5 with B5RBP2, -3, and -5. The iron response element (IRE) is included as a negative control (B). Approximately 5 fmole of biotin-labeled bait RNA and protein concentrations ranging from 30 to 250 nM were used in each reaction.
The 3′ UTR of StBEL5 is involved in several aspects of RNA metabolism and is replete with potential binding motifs (Banerjee et al., 2009). To identify shorter binding regions within the 3′ UTR that may be involved in protein/RNA interaction, truncated bait sequences were utilized in the β-galactosidase assay of the Y3H system (Figure 6). Three truncated sequences were used based on their conserved sequence motifs and their coverage of the UTR (Figure 7). The 5′ D1 sequence is enriched for CU motifs (underlined sequence, Figure 7). T2 contains several UAGU motifs (Figure 7, boxed), and the UA-bait sequence contains a number of uracil/adenine runs (underlined sequence, Figure 7). Overall, the greatest β-galactosidase activity was observed for the full-length 3′ UTR (Figure 6). Based on β-galactosidase activity, B5RBP3, -5, -6, and -7 exhibited the strongest interaction with sequence located toward the 5′ end of the UTR (D1 and T2 baits). B5RBP1 and -2 exhibited the strongest interaction with sequence located toward the 3′ end of the UTR (T2 and UA baits). B5RBP4, the potato ortholog of eIF5A, exhibited equivalent strength of activity with all three truncated baits suggesting a degree of non-specific binding. These results with the potato eIF5A are consistent with previous work showing that a pumpkin form of eIF5A exhibited RNA-binding that was non-sequence-specific in nature (Ma et al., 2010).
Figure 6. Analysis of β-galactosidase activity for the interaction of select StBEL5 RNA-binding proteins (B5RBPs) with sequences within the 3′ UTR of StBEL5. The values in each graph are normalized to activity relative to the full-length UTR. 3′ UTR, full-length UTR; D1; a 178-nt sequence within the UTR starting from the stop codon; T2, a UAGU-rich region within the UTR; UA, UA-rich region toward the 3′ end of the UTR; MS2-1, RNA sequence from the bait vector serving as a negative control. See Figure 7 for details on these truncated StBEL5 bait sequences.
Figure 7. Schematic diagram of the bait RNA sequences from the 3′ UTR of StBEL5 used in Figure 6. Included are three truncated sequences within the full-length 503-nt 3′ UTR: the 178-nt D1 sequence, the 106-nt T2 region, and the 123-nt UA-rich region (UA). Both UA-and CU-rich motifs are underlined. The UAGU motifs of the T2 sequence are boxed.
The Y3H system has been established as an efficient method for selecting protein partners of RNA from among thousands of putative partners, and for assaying binding affinity of specific RNA/protein interactions. With modifications, this system has been adapted for screening RNA/RNA interactions (Piganeau and Schroeder, 2006), for identifying protein/small signaling molecule complexes (Cottier et al., 2011), and for testing multi-component interactions (Bernstein et al., 2002). Although numerous false positives may arise during the screening process, there are several levels of selection that may be utilized to eliminate these. These include nutrient selection, HIS3 activation, addition of 3-aminotriazole, and 5-fluororotic acid to the media, and the induction of lacZ. As shown previously (Hook et al., 2005), both HIS3 and lacZ expression levels are directly related to binding affinity and may be used to assess the robustness of specific RNA/protein interactions (Mahajan et al., 2012) or to map specific motifs present in either bait or target sequences (Edwards et al., 2001; Mori et al., 2008; Stumpf et al., 2008).
Factors known to affect in vivo interactions include the intrinsic affinity of bait RNA and target protein, the length of the bait RNA, and accessibility of the insert to the target (Wurster and Maher III, 2010). RNA sequences less than 150 nt generally produce the most substantial and specific reporter activation. The inclusion of additional sequence can lead to a reduction in signal. In this study, the entire 503-nt 3′ UTR of StBEL5 was utilized to include a wide range of interactions. With this approach, subsequent analyses of shorter RNA sequences in RNA/protein interactions using either RNA gel-shift or Y3H assays may be used to identify the location of specific binding motifs. The source of the protein expression library will also play an important role in identification of protein partners. Here the use of a leaf cDNA library would provide wide coverage but could preclude the identification of important interactions with rare phloem-mobile proteins. Three of the RBPs compiled in the final list here (B5RBP1, -5, and -7) appear to function in the intracellular transport of RNAs. Two RBPs identified in the current study, eIF5A and HSP70, were isolated from a phloem-mobile RNP complex of pumpkin (Ham et al., 2009). A protein very close in sequence match to the glycine-rich RBP (B5RBP2) was also identified in pumpkin phloem sap (Lin et al., 2009). Both binding affinity and abundance of the target protein would impact the results of the Y3H screen. As demonstrated by StLSH10, despite its relatively low transcript level in leaves (Table 1), its affinity for the BEL5 UTR bait was demonstrated to be relatively high. The very high transcript levels observed for StLSH10 (B5RBP3) in stolons and young tubers (149 and 590 FPKMs, respectively) also suggest a possible role in tuber development. Screening a library of proteins from a sink organ like the tuberizing stolon would further expand our understanding of the RNP complex responsible for delivery of the mobile StBEL5 transcript to its functional site. Overall, these results demonstrate the utility of the Y3H system in screening and verifying interactions between target RNAs and putative RBPs.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Thanks to Kate Lueders for her valuable technical assistance and to Marv Wickens for graciously providing the Y3H vectors and yeast strains. Thanks also to Tian Lin and Pooja Sharma for their help in preparing the LSH phylogenetic tree and alignment. This research was supported by the NSF Plant Genome Research Program award no. DBI-0820659.
Banerjee, A. K., Yu, Y., Chatterjee, M., Suh, S. G., Miller, W. A., and Hannapel, D. J. (2006). Dynamics of a mobile RNA of potato involved in a long-distance signaling pathway. Plant Cell 18, 3443–3457.
Baroja-Fernandez, E., Munoz, F. J., Montero, M., Etxeberria, E., Sesma, M. T., Ovecka, M., Bahaji, A., Ezquer, I., Li, J., Prat, S., and Pozueta-Romero, J. (2009). Enhancing sucrose synthase activity in transgenic potato (Solanum tuberosum L.) tubers results in increased levels of starch, ADPglucose and UDPglucose and total yield. Plant Cell Physiol. 50, 1651–1662.
Bieniawska, Z., Paul Barratt, D. H., Garlick, A. P., Thole, V., Kruger, N. J., Martin, C., Zrenner, R., and Smith, A. M. (2007). Analysis of the sucrose synthase gene family in Arabidopsis. Plant J. 49, 810–828.
Campalans, A., Kondorosi, A., and Crespi, M. (2004). Enod40, a short open reading frame-containing mRNA, induces cytoplasmic localization of a nuclear RNA binding protein in Medicago truncatula. Plant Cell 16, 1047–1059.
Corral-Debrinski, M., Blugeon, C., and Jacq, C. (2000). In yeast, the 3′ untranslated region or the presequence of ATM1 is required for the exclusive localization of its mRNA to the vicinity of mitochondria. Mol. Cell. Biol. 20, 7881–7892.
Cottier, S., Mönig, T., Wang, Z., Svoboda, J., Boland, W., Kaiser, M., and Kombrink, E. (2011). The yeast three-hybrid system as an experimental platform to identify proteins interacting with small signaling molecules in plant cells: potential and limitations. Front. Plant Sci. 2:101. doi:10.3389/fpls.2011.00101
Doroshenk, K. A., Crofts, A. J., Morris, R. T., Wyrick, J. J., and Okita, T. W. (2009). Proteomic analysis of cytoskeleton-associated RNA binding proteins in developing rice seed. J. Proteome Res. 8, 4641–4653.
Fallahi, H., Scofield, G. N., Badger, M. R., Chow, W. S., Furbank, R. T., and Ruan, Y. L. (2008). Localization of sucrose synthase in developing seed and siliques of Arabidopsis thaliana reveals diverse roles for SUS during development. J. Exp. Bot. 59, 3283–3295.
Ferrandon, D., Elphick, L., Nusslein-Volhard, C., and St Johnston, D. (1994). Staufen protein associates with the 3′ UTR of bicoid mRNA to form particles that move in a microtubule-dependent manner. Cell 79, 1221–1223.
Fu, H., Kim, S. Y., and Park, W. D. (1995). High-level tuber expression and sucrose inducibility of a potato Sus4 sucrose synthase gene require 5′ and 3′ flanking sequences and the leader intron. Plant Cell 7, 1387–1394.
Gonsalvez, G. B., Lehmann, K. A., Ho, D. K., Stanitsa, E. S., Williamson, J. R., and Long, R. M. (2003). RNA-protein interactions promote asymmetric sorting of the ASH1 mRNA ribonucleoprotein complex. RNA 9, 1383–1399.
Ham, B. K., Brandom, J. L., Xoconostle-Cázares, B., Ringgold, V., Lough, T. J., and Lucas, W. J. (2009). A polypyrimidine tract binding protein, pumpkin RBP50, forms the basis of a phloem-mobile ribonucleoprotein complex. Plant Cell 21, 197–215.
Hunt, A. G., Xu, R., Addepalli, B., Rao, S., Forbes, K. P., Meeks, L. R., Xing, D., Mo, M., Zhao, H., Bandyopadhyay, A., Dampanaboina, L., Marion, A., Von Lanken, C., and Li, Q. Q. (2008). Arabidopsis mRNA polyadenylation machinery: comprehensive analysis of protein-protein interactions and gene expression profiling. BMC Genomics 9, 220. doi:10.1186/1471-2164-9-220
Hwang, M. S., Kim, S. H., Lee, J. H., Bae, J. M., Paek, K. H., and Park, Y. I. (2005). Evidence for interaction between the 2a polymerase protein and the 3a movement protein of cucumber mosaic virus. J. Gen. Virol. 86, 3171–3177.
Lai, W. S., Parker, J. S., Grissom, S. F., Stumpo, D. J., and Blackshear, P. J. (2006). Novel mRNA targets for tristetraprolin (TTP) identified by global analysis of stabilized transcripts in TTP-deficient fibroblast. Mol. Cell Biol. 26, 9196–9208.
Lee, J., He, K., Stolc, V., Lee, H., Figueroa, P., Gao, Y., Tongprasit, W., Zhao, H., Lee, I., and Deng, X. (2007). Analysis of transcription factor HY5 genomic binding sites revealed its hierarchical role in light regulation of development. Plant Cell 19, 731–749.
Li, C., Gu, M., Shi, N., Zhang, H., Yang, X., Osman, T., Liu, Y., Wang, H., Vatish, M., Jackson, S., and Hong, Y. (2011). Mobile FTmRNA contributes to the systemic florigen signalling in floral induction. Sci. Rep. 1, 73.
Lin, M. K., Lee, Y. J., Lough, T. J., Phinney, B. S., and Lucas, W. J. (2009). Analysis of the pumpkin phloem proteome provides insights into angiosperm sieve tube function. Mol. Cell Proteomics 8, 343–356.
Ma, Y., Miura, E., Ham, B. K., Cheng, H. W., Lee, Y. J., and Lucas, W. J. (2010). Pumpkin eIF5A isoforms interact with components of the translational machinery in the cucurbit sieve tube system. Plant J. 64, 536–550.
Maniataki, E., Martinez de Alba, E., Sägesser, R., Tabler, M., and Tsagris, M. (2003). Viroid RNA systemic spread may depend on the interaction of a 71-nucleotide bulged hair-pin with the host protein VirP1. RNA 9, 346–354.
Moore, F. L., Jaruzelska, J., Fox, M. S., Urano, J., Firpo, M. T., Turek, P. J., Dorfman, D. M., and Pera, R. A. (2003). HumanPumilio-2 is expressed in embryonic stemcells and germcells and interacts with DAZ (Deleted in AZoospermia) and DAZ-like proteins. Proc. Natl. Acad. Sci. U.S.A. 100, 538–543.
Peat, T. S., Newman, J., Waldo, G. S., Berendzen, J., and Terwilliger, T. C. (1998). Structure of translation initiation factor 5A from Pyrobaculum aerophilum at 1.75 A resolution. Structure 6, 1207–1214.
Xu, X., Pan, S., Cheng, S., Zhang, B., Mu, D., Ni, P., Zhang, G., Yang, S., Li, R., Wang, J., Orjeda, G., Guzman, F., Torres, M., Lozano, R., Ponce, O., Martinez, D., De la Cruz, G., Chakrabarti, S. K., Patil, V. U., Skryabin, K. G., Kuznetsov, B. B., Ravin, N. V., Kolganova, T. V., Beletsky, A. V., Mardanov, A. V., Di Genova, A., Bolser, D. M., Martin, D. M., Li, G., Yang, Y., Kuang, H., Hu, Q., Xiong, X., Bishop, G. J., Sagredo, B., Mejía, N., Zagorski, W., Gromadka, R., Gawor, J., Szczesny, P., Huang, S., Zhang, Z., Liang, C., He, J., Li, Y., He, Y., Xu, J., Zhang, Y., Xie, B., Du, Y., Qu, D., Bonierbale, M., Ghislain, M., Herrera Mdel, R., Giuliano, G., Pietrella, M., Perrotta, G., Facella, P., O’Brien, K., Feingold, S. E., Barreiro, L. E., Massa, G. A., Diambra, L., Whitty, B. R., Vaillancourt, B., Lin, H., Massa, A. N., Geoffroy, M., Lundback, S., DellaPenna, D., Buell, C. R., Sharma, S. K., Marshall, D. F., Waugh, R., Bryan, G. J., Destefanis, M., Nagy, I., Milbourne, D., Thomson, S. J., Fiers, M., Jacobs, J. M., Nielsen, K. L., Sønderkær, M., Iovene, M., Torres, G. A., Jiang, J., Veilleux, R. E., Bachem, C. W., de Boer, J., Borm, T., Kloosterman, B., van Eck, H., Datema, E., Hekkert, B. L., Goverse, A., van Ham, R. C., Visser, R. G., and Potato Genome Sequencing Consortium. (2011). Genome sequence and analysis of the tuber crop potato. Nature 475, 189–195.
Steiner, T., Kaiser, J. T., Marinkovic, S., Huber, R., and Wahl, M. C. (2002). Crystal structures of transcription factor NusG in light of its nucleic acid- and protein-binding activities. EMBO J. 21, 4641–4653.
Sun, J., Jiang, H., Xu, Y., Li, H., Wu, X., Xie, Q., and Li, C. (2007). The CCCH-type zinc finger proteins AtSZF1 and AtSZF2 regulate salt stress responses in Arabidopsis. Plant Cell Physiol. 48, 1148–1158.
Takeda, S., Hanano, K., Kariya, A., Shimizu, S., Zhao, L., Matsui, M., Tasaka, M., and Aida, M. (2011). Cup-shaped cotyledon1 transcription factor activates the expression of LSH4 and LSH3, two members of the ALOG gene family, in shoot organ boundary cells. Plant J. 66, 1066–1077.
Thio, G. L., Ray, R. P., Barcelo, G., and Schupbach, T. (2000). Localization of gurken RNA in Drosophila oogenesis requires elements in the 5′ and 3′ regions of the transcript. Dev. Biol. 15, 435–446.
Wachter, A., Rühl, C., and Stauffer, E. (2012). The role of polypyrimidine tract-binding proteins and other hnRNP proteins in plant splicing regulation. Front. Plant Sci. 3:81. doi:10.3389/fpls.2012.00081
Xu, A., and Chen, K. Y. (2001). Hypusine is required for a sequence-specific interaction of eukaryotic initiation factor 5A with post systematic evolution of ligands by exponential enrichment RNA. J. Biol. Chem. 276, 2555–2561.
Zhao, L., Nakazawa, M., Takase, T., Manabe, K., Kobayashi, M., Seki, M., Shinozaki, K., and Matsui, M. (2004). Overexpression of LSH1, a member of an uncharacterised gene family, causes enhanced light regulation of seedling development. Plant J. 37, 694–706.
Ziemienowicz, A., Haasen, D., Staiger, D., and Merkle, T. (2003). Arabidopsis transportin1 is the nuclear import receptor for the circadian clock-regulated RNA-binding protein AtGRP7. Plant Mol. Biol. 53, 201–212.
Zrenner, R., Salanoubat, M., Willmitzer, L., and Sonnewald, U. (1995). Evidence of the crucial role of sucrose synthase for sink strength using transgenic potato plants (Solanum tuberosum L.). Plant J. 7, 97–107.
Figure A1. Coomassie-stained gels of protein purification of three B5RBPs. Each protein was induced in E. coli, the soluble fraction was extracted, and the protein purified once or twice (B5RBP3) using the HisPur Cobalt purification kit (Pierce). The farthest right-hand lanes (arrows) represent from 15 to 100 ng of protein loaded and these samples were used for the RNA EMSAs. (+) indicates the addition of IPTG; (−) no IPTG; S, soluble fraction; IS, insoluble fraction.
Table A1. Full list of the screened genes with their putative identities and concentration of 3-AT for screening.
Keywords: mobile RNA, potato, Solanum tuberosum, yeast three-hybrid, StBEL5, BEL1 family
Citation: Cho SK, Kang I-H, Carr T and Hannapel DJ (2012) Using the yeast three-hybrid system to identify proteins that interact with a phloem-mobile mRNA. Front. Plant Sci. 3:189. doi: 10.3389/fpls.2012.00189
Received: 20 June 2012; Paper pending published: 16 July 2012;
Accepted: 02 August 2012; Published online: 27 August 2012.
Edited by:Yiguo Hong, Hangzhou Normal University, China
Reviewed by:Stephen Jackson, Warwick University, UK
Vipaporn Phuntumart, Bowling Green State University, USA
Copyright: © 2012 Cho, Kang, Carr and Hannapel. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in other forums, provided the original authors and source are credited and subject to any copyright notices concerning any third-party graphics etc.
*Correspondence: David J. Hannapel, Plant Biology Major, Iowa State University, Ames, IA 50011, USA. e-mail: email@example.com