A membrane protein/signaling protein interaction network for Arabidopsis version AMPv2.

Interactions between membrane proteins and the soluble fraction are essential for signal transduction and for regulating nutrient transport. To gain insights into the membrane-based interactome, 3,852 open reading frames (ORFs) out of a target list of 8,383 representing membrane and signaling proteins from Arabidopsis thaliana were cloned into a Gateway-compatible vector. The mating-based split ubiquitin system was used to screen for potential protein-protein interactions (pPPIs) among 490 Arabidopsis ORFs. A binary robotic screen between 142 receptor-like kinases (RLKs), 72 transporters, 57 soluble protein kinases and phosphatases, 40 glycosyltransferases, 95 proteins of various functions, and 89 proteins with unknown function detected 387 out of 90,370 possible PPIs. A secondary screen confirmed 343 (of 386) pPPIs between 179 proteins, yielding a scale-free network (r(2) = 0.863). Eighty of 142 transmembrane RLKs tested positive, identifying 3 homomers, 63 heteromers, and 80 pPPIs with other proteins. Thirty-one out of 142 RLK interactors (including RLKs) had previously been found to be phosphorylated; thus interactors may be substrates for respective RLKs. None of the pPPIs described here had been reported in the major interactome databases, including potential interactors of G-protein-coupled receptors, phospholipase C, and AMT ammonium transporters. Two RLKs found as putative interactors of AMT1;1 were independently confirmed using a split luciferase assay in Arabidopsis protoplasts. These RLKs may be involved in ammonium-dependent phosphorylation of the C-terminus and regulation of ammonium uptake activity. The robotic screening method established here will enable a systematic analysis of membrane protein interactions in fungi, plants and metazoa.

Introduction

Membrane proteins play crucial roles in many biological processes. They control the cell's permeability to a myriad of compounds, mediate uptake and release of ions, metabolites and proteins across subcellular membranes, are involved in vesicle fusion, and are responsible for sensing chemical and physical stimuli from the environment (nutrients, hormones, light, pH, etc.).

Both regulation of transport activity and communication involve interactions of membrane proteins either with each other or with proteins in the adjacent soluble compartments. The development and use of a wide spectrum of new tools for testing protein–protein interactions has been key to much of the progress in understanding cellular functions. Despite the importance, we know only a fraction of the total interactome, and our knowledge is specifically limited in the case of membrane proteins. Large-scale protein interaction screens in yeast or animal systems, which in part made use of the classical yeast two hybrid (Y2H) system, have provided at least a partial overview over the potential interactome in these organisms (Rual et al., 2005; Yu et al., 2008; Simonis et al., 2009). The pioneering work of Vidal's group has provided most of the interactome data so far. At present, the human interactome comprises ∼130,000 binary interactions, with most of the interactome still remaining to be identified (Venkatesan et al., 2009). A bona fide plant interactome is still outstanding; therefore, systems approaches have to rely mainly on data for interactions of homologs from other organisms (Geisler-Lee et al., 2007). Particularly, knowledge of the membrane protein interactome is limited in all organisms, because conventional high throughput Y2H assays are not designed to detect potential protein–protein interactions (pPPIs) of membrane proteins.

To gain insights into the membrane-based interactome, the mating-based split ubiquitin system (mbSUS) was developed, which specifically detect interactions of membrane proteins and of membrane proteins with signaling proteins (Obrdlik et al., 2004; Miller et al., 2005). The split ubiquitin system concept relies on sequestration of a transcription factor to the membrane and its subsequent release when two proteins interact. It is based on peptide complementation and uses a split ubiquitin. The N-terminal domain of ubiquitin (Nub) when co-expressed with its C-terminal half (Cub) reconstitutes functional ubiquitin (Johnsson and Varshavsky, 1994). A Nub mutant, NubG (containing mutation Ile13Gly) with reduced affinity (compared to wild-type Nub) to the Cub moiety is unable to reconstitute functional ubiquitin unless brought into the vicinity of the Cub domain via interaction of two fusion partners (Johnsson and Varshavsky, 1994). An artificial transcription factor, PLV (protease A – LexA – VP16), is fused in frame to the C-terminal Cub moiety. When two proteins interact, the NubG and Cub moieties form a functional ubiquitin, and endogenous ubiquitin-specific proteases release the transcription factor into the cytosol. The transcription factor diffuses into the nucleus where it activates the transcription of the reporter genes (His3, Ade2, and LacZ). The split ubiquitin system was further improved for high throughput screens by introducing a mating approach (mbSUS) (Obrdlik et al., 2004). Subsequently, mbSUS was further modified to incorporate the Gateway recombination system for simplified cloning.

The split ubiquitin system has successfully been used to analyze pPPIs among 705 proteins annotated as integral to a cellular membrane in yeast (Miller et al., 2005). The screen identified 1,985 putative interactions among 536 proteins. Recently, a plant split ubiquitin system was developed and used to test interactions of translocon complex at the outer chloroplast membrane (Rahim et al., 2009). The split ubiquitin system identified the potential for oligomerization of plant transporters such as potassium channels, ammonium transporters, calcium/proton antiporters, H+/sucrose cotransporters, and a mammalian phosphate transporter (Schulze et al., 2003; Obrdlik et al., 2004; Gisler et al., 2008; Zhao et al., 2009). The identification of AMT oligomerization was a key step toward the characterization of the novel allosteric regulation of AMT activity by the C-terminus of neighboring subunits in a trimeric complex (Loqué et al., 2007). Extracellular ammonium was found to trigger phosphorylation of a specific threonine (T460) in the C-terminus of AMT1;1 (Lanquar et al., 2009). Additional phosphorylated sites were found in the C-terminus using phosphoproteomic studies (Lanquar et al., 2009). The split ubiquitin system may thus allow us to identify the respective protein kinases (Lanquar and Frommer, 2010). The split ubiquitin system was also successful in identifying a physical interaction between the dopamine transporter and the synaptic vesicle protein synaptogyrin-3 (Egana et al., 2009), as well as interactions between five different mammalian transporters and PDZ domain-containing partners (Gisler et al., 2008). Harter's group used mbSUS to show that ethylene receptors form homomeric and heteromeric protein complexes (Grefen et al., 2008). Similarly, an interactions between cytosolic glutamine synthetase of soybean root nodules and the aquaglyceroporin nodulin-26 was identified (Masalkar et al., 2010). In most of these examples, data obtained with the split ubiquitin system were confirmed independently by a variety of orthogonal interaction methods, supporting the reliability of this assay system.

The goal of the work presented here was to lay the foundation for a systematic screen of the Arabidopsis membrane protein interactome and its interface with key signaling proteins. After establishing a target list of genes, the first step was to clone the open reading frames (ORFs) from Arabidopsis thaliana. A Gateway-based strategy was used to create 3,852 Entry clones that were then mobilized into the mbSUS destination vectors. After yeast transformation, the Cub fusions were tested for “false positive” and “false negative” read-out. A pilot screen (AMPv1, Associome of membrane proteins version 1; Chen et al., personal communication) suggested that interaction tests in 96-well would not provide sufficient throughput for a systematic analysis of membrane signaling protein interactions from Arabidopsis. Here, each Cub clone was then mated against a Nub-fusion collection in a 384-well format and tested for interactions in a 1536-well format. In this screen (AMPv2), we tested 90,370 potential interactors and found 343 pPPIs between 179 proteins. Interactions between the ammonium transporter AMT1;1 and two receptor-like kinases (RLKs) identified in the screen were independently confirmed using a split luciferase assay.

Materials and Methods

Arabidopsis first-strand cDNA synthesis

RNA was extracted from Arabidopsis seedlings and flowers as described by Downing et al. (1992). DNA was removed by DNaseI treatment (Invitrogen) followed by ethanol precipitation. The polyA+ fraction was purified using DynaBeads Oligo(dT) (Invitrogen) and cDNA was synthesized using the SuperScript III kit (Invitrogen). mRNA for seedlings and flowers was isolated separately and mixed in a 1:1 ratio after the synthesis of the first-strand cDNA. The cloning scheme is illustrated in Figure 1.

Figure 1

Outline of the cloning and mbSUS screen. Selected ORFs were amplified from Arabidopsis seedling cDNA and cloned into pDONR221 (or pDONR221-f1) by Gateway recombination. ORFs that passed the quality control criteria were mobilized into Cub and NubG destination vectors and introduced into THY.AP4 and THY.AP5 yeast strains. Yeast expressing the Cub clones were pre-screened to remove “false positives” and “false negatives.” Remaining Cub clones were mated with NubG clones that had been arrayed into two 1536-well plates. Successful mating was selected on DS media. After 3 days, diploids were replica-plated to IS media selecting for activation of the His3 and Ade2 markers and growth was scored.

Cloning of arabidopsis ORFs

Primers for target genes were designed in bulk with a command-line version of Primer3 (Rozen and Skaletsky, 2000). ORFs (without stop codon) were amplified by PCR from Arabidopsis first-strand cDNA using KOD polymerase (Novagen). PCR products were isolated from crystal violet-stained (Rand, 1996) agarose gels using a NucleoSpin extraction kits (Macherey-Nagel). Purified PCR products were cloned into pDONR221 or pDONR221-f1 (f1 origin of replication was inserted counter-clockwise into the unique PciI restriction site after removing 5′-overhangs with T4 DNA polymerase). Colony PCR was used to identify plasmids containing an insert with the expected size (low resolution separation on 96 sample gels). Cloning products were sequenced from both ends with a median read length of 570 Phrep20 bases into the gene (Beckman Coulter Genomics). Plasmids with the correct sequences, or clones that contained multiple silent mutations but no more than a single amino acid change were isolated, verified by restriction digest analysis (PvuII) and transferred into the mbSUS destination vectors by in vitro LR recombination. Clones received from third parties (including ABRC; Underwood et al., 2006; Popescu et al., 2009) were reconfirmed by restriction analysis and/or re-sequencing.

Construction of gateway-compatible mbSUS destination vectors

The mbSUS interaction trap was originally created for cloning ORFs using in vivo recombination directly into the mbSUS vectors in yeast (Obrdlik et al., 2004). The vectors contained attB1 and attB2 site allowing for an easy transfer of the ORFs to other destination vectors. To increase the efficiency of cloning for this project, the KanMX containing vectors pXNgate22-HA and pMetYCgate (Obrdlik et al., 2004) were treated with BP clonase in the presence of pDONR221 to obtain mbSUS Gateway (GW) destination vectors. Bacteria containing plasmids (with Gateway cassette, ccdB and chloramphenicol resistance genes) were selected in ccdB survival cells. The new mbSUS Gateway destination vectors were named pXN22_GW and pMetYC_GW, respectively. Identity of the plasmids was confirmed by sequencing the borders of the region of interest. Plasmids are available through ABRC and maps are available through the project's website.

Yeast transformation

The yeast strains THY.AP4 [MATa ura3 leu2 lexA::lacZ::trp1 lexA::HIS3 lexA::ADE2] and THY.AP5 [MATα URA3 leu2 trp1 his3 loxP::ade2] (both derived from CEN.PK113-17D [MATα URA3 leu2 trp1 his3 ADE2]) were used in combination with pMetYC_GW (Cub Gateway destination vector), and pXN22_GW (NubG Gateway destination vector) respectively (Obrdlik et al., 2004). Yeast strains were transformed with respective vectors using the LiAc method (Gietz and Woods, 2002) optimized for high throughput. Briefly, strains were inoculated in 15 ml of YPD-adenine medium and diluted the following day to OD600 0.1 (final volume of 300 ml). Strains were incubated at 28°C until OD600 0.5–0.6 (∼6 h). Cells were harvested by centrifugation and resuspended in 100 mM LiAc in Tris–EDTA (10 mM Tris, 1 mM EDTA pH 7.5). The original 300 ml of culture resulting in 1.8 ml LiAc solution was used for one 96-well plate. In a PCR microplate (E & K Scientific), 1 μl plasmid DNA (50–100 ng/μl) was dispensed into each well of a 96-well microplate, 99 μl of a master mix [7 ml 50% (w/v) polyethylene glycol 4000 (Fluka #81240), 0.8 ml 1 M LiAc in TE, 1.8 ml cells, 0.3 ml freshly denatured ssDNA (8–10 mg/ml)] was added and incubated for 30 min at 28°C, followed by a heat shock at 42°C for 13 min. Cells were washed in water and resuspended in 75 μl water. Ten microliters of cells were plated on solid N (THY.AP5) or C (THY.AP4) media using a Hydra96 (Art Robbins Instruments, Sunnyvale). After growth for 3 days, cells were transferred to liquid N and C media in 96-well plates, respectively. After 3 days of growth in liquid culture, glycerol stocks were prepared and cells were transferred back to solid N and C media for further analysis. Four hundred and ninety unique ORFs (Table S1 in Supplementary Material) were mobilized into the mbSUS destination vectors by Gateway LR in vitro recombination. Membrane-localized proteins with cytosolic C-termini were fused to Cub-PLV: 443 ORFs were mobilized into the NubG- and 204 in the Cub-vector.

Yeast media

YPD media, synthetic complete (SC) media, and drop-out (minus histidine, adenine, uracil, tryptophan, leucine, and methionine) media were prepared according to standard protocols (Adams et al., 1998). Depending on the selection, the SC media was supplemented with 20 mg/L adenine hemisulfate, 20 mg/L histidine-HCl, 500 μM methionine, 20 mg/L tryptophan, 20 mg/L uracil, and/or 240 mg/L leucine. Cells harboring the Cub plasmids were grown on SC supplemented with tryptophan, histidine, adenine, and uracil (C media) and cells with the NubG plasmid on SC with leucine, adenine, and histidine (N media). Diploid cells were selected on SC media supplemented with adenine and histidine (DS media), while selection for pPPIs was performed on SC supplemented with 150 or 500 μM methionine (IS-150 and IS-500 media) or without methionine (IS-0 media). Cub fusions are under the control of the methionine repressible MET25 promoter; increasing the amount of methionine in the media decreases the expression level and hence, increases the stringency.

Identification of “false positive” and “false negative” candidates

Cub clones were pre-screened to identify “false positive” and “false negative” candidates (Figure 1). Cub clones in THY.AP4 were arrayed on C media in 384-well format using the RoToR HDA yeast pinning robot (Singer Instruments, UK). Recombinant THY.AP4 cells were mated with a yeast strain (THY.AP5) expressing either soluble NubWT (control for “false negatives”) or soluble NubG alone (control for “false positives”) on N media. Cells were allowed to mate on YPD media for 24 h at 28°C and then replica-plated onto DS media. Interactions were selected on IS-0, IS-150, and IS-500 media. Growth was documented by scanning plates on a flatbed scanner (CanoScan 8400F, Canon). Only Cub fusions showing no growth when mated with a strain containing soluble NubG (Cub-fusion that do not autoactivate the reporters) and growth for soluble NubWT (indicating protein expression) were used for the interaction screen.

Robotics

NubG clones were re-arrayed in 384-well format (RoToR HDA pinning robot, Singer Instruments, UK) and were mated with individual Cubs. For this purpose, individual Cub clones were grown in 15 ml C media, distributed into a 96-well plate and arrayed from liquid media onto solid C media in 384-well format. Plates containing cells with Nub and Cub clones were mated (see above). After 24 h at 28°C, cells were plated on DS media to select for diploid cells and grown for 3 days at 28°C. Cells were re-pinned in 1536-well format on IS media containing 0, 150, or 500 μM methionine for selection of pPPIs (IS-0, IS-150, and IS-500, respectively). Positive spots were scored. To independently verify the pPPIs, the original clones identified as part of the pPPI network were re-arrayed in 96-well format using a BioRobot 3000 (Qiagen) and tested for His3, Ade2, and the third marker LacZ.

β-galactosidase assay

The β-galactosidase assay, which tests for the third mbSUS marker lacZ, was performed with cells grown on a nitrocellulose-overlaid plate over an IS media plate (IS-0, IS-150, and IS-500). The β-galactosidase assay was performed as described (Obrdlik et al., 2004). The nitrocellulose filters were scanned on a flatbed scanner (CanoScan 8400F, Canon) and the intensity of the spots was quantified using GenePix v6.1 software (Molecular Devices, Sunnyvale) after conversion of the images to black and white. For each treatment (IS-0, IS-150, and IS-500 media), the average intensity and standard deviation were calculated; the cut-off value was set to the average minus one standard deviation. For a given pair, an intensity value higher than the cut-off for at least one treatment was considered a pPPI for the final network. The resulting network was named AMPv2 and can be queried through the project website.

Split luciferase interaction assay

Selected pPPIs were tested in Arabidopsis protoplasts using the orthogonal Renilla luciferase complementation assay (Fujikawa and Kato, 2007; Kato et al., 2010). SYP122 (AT3G52400)/VAMP721 (AT1G04750), and SYP122/PHT4 (AT2G38940) pairs were used as a positive and negative controls, respectively (Fujikawa and Kato, 2007). Arabidopsis protoplasts were isolated from leaves of 4-week-old plants according to Kato and Jonse (2010). Renilla luciferase activity (split luciferase assay) was measured in 96-well white plate using the ViviRen substrate according to Kato and Jonse (2010) using a microplate luminometer (Promega).

Results

The target gene list

To select proteins for the membrane protein/signaling protein interaction screen, 26,922 protein sequences from Arabidopsis (assembly TAIR 7) were grouped into protein families using the Markov Cluster (MCL) algorithm (Enright et al., 2002) based on pre-computed sequence similarity information. A total of 3,360 clusters containing 21,807 proteins, as well as 5,115 singletons were identified. The composition of each cluster was analyzed using multiple sequence alignments. Proteins that aligned only partially (less than ∼75%) with other proteins of a cluster were removed to ensure that each cluster contained proteins of one family. Clusters of interest were selected based on locus descriptions from TAIR, published data, and topology as predicted by the Aramemnon database (Schwacke et al., 2003) to select (1) plasma membrane proteins, and (2) proteins involved in signaling and in protein turnover. Certain classes of soluble proteins not involved in membrane signal transduction or ubiquitination were excluded, including transcription factors (families retrieved from TAIR), pentatricopeptide domain-containing proteins, plastidial and mitochondrial proteins. Eight hundred and forty-six clones were obtained from third parties involved in various processes (310 were included in the present screen) were included. In addition, 209 proteins that had been previously shown to localize to membranes by proteomic or cell biological approaches were added to the target list.

Protein topologies were assigned using published data and similarities to known proteins, or the number and orientation of transmembrane spans was predicted using Topcons (Bernsel et al., 2009), TMHMM (Krogh et al., 2001), and Phobius (Kall et al., 2007). Proteins predicted to have both N- and C-terminus outside the cytosol were eliminated leading to a final list of 8,383 candidates. This list contains 4,332 proteins with ≥1 transmembrane domains and are thus suitable as Nub and Cub fusions. Four hundred and thirty-two proteins were predicted to be myristoylated and can also be used as Nub and Cub fusions, while 4,051 soluble proteins were tested only as Nub fusions.

Cloning of arabidopsis ORFs for the mbSUS screen

The target list was used as the basis to clone a maximal number of ORFs. For cloning, RT-PCR was performed on Arabidopsis thaliana Col-0 cDNA from seedlings and flowers. A variety of Gateway™ cloning strategies was used with different success rates. Initially, ORFs were cloned using the TOPO®TA system in pCR8 vector (Invitrogen). However, the non-directional approach required additional analyses to identify clones suitable for screening. Cloning into pENTR-D (Invitrogen; directional cloning kit) did not yield a prevalence of correctly oriented clones; furthermore, unexpected 13 bp insertions at the 5′-end led to frame shifts. By comparison, in vitro recombination using BP cloning (Gateway, Invitrogen) was efficient. For Gateway cloning, primers with attB1 and attB2 sites were designed. Resulting clones were tested for correct insert size by restriction digests and sequenced from both ends. Clones selected for further processing were allowed to contain multiple silent mutations but no more than a single amino acid change. Based on these quality control (QC) criteria (size and sequence), the overall cloning success was ∼48%. For clones that passed QC, the error rate was 0.028%, or one error per 3,624 bases (292 errors in 1,924 clones). To date, 3,852 Entry clones (representing 3,571 unique AGIs) have been generated (846 genes as thirty party donations) and 2,106 ORFs have been made available through ABRC (cf. Table S2 in Supplementary Material). A set of Gateway-compatible mbSUS destination vectors was constructed, and the first 490 genes were mobilized into pMetYC_GW (204 Cub clones) and pNX22_GW destination vectors (443 Nub clones; Overlap: 157 ORFs present in both Nub and Cub vectors); this pilot study reports potential protein–protein interactions corresponding to a set of 490 proteins.

Identification of “false positives/negative” candidates

A pre-screen was performed to identify potential Cub fusions producing “false negative” and “false positive” reports. Correctly folded and efficiently expressed Cub fusions are expected to induce reporter activity when mated with cells expressing the soluble wild-type version of Nub (NubWT), which has a strong affinity for the Cub domain, (note that mbSUS uses a mutant NubG with reduced affinity for Cub; Johnsson and Varshavsky, 1994). Sources of potential “false negative” interactions include: gene products that expose the Cub-PLV domain to an intracellular compartment such as ER or Golgi and thus do not produce an interaction in the cytosol, or proteins in which the Cub-PLV domain is not folded properly and thus cannot be recognized by ubiquitin-specific proteases. “False positives” were identified through interaction tests with soluble NubG, specifically gene products that activated the reporters in the absence of an interaction partner. The pre-screen lead to elimination of 69 out of 204 ORFs cloned in the Cub-vector, resulting in 135 Cub clones for further analysis (Table S1 in Supplementary Material). The rate of clones eliminated in the pre-screen was slightly lower (66.2% retained vs. 51.7%) compared to the screen performed with endogenous yeast membrane proteins (Miller et al., 2005).

Robotic screen for interactions

A robotic mating procedure was established in 384-well plates. Subsequently, interaction tests were performed in 1536-well format using a Singer RoToR HDA replica-plating robot. To establish the conditions for the high throughput screen, 443 NubG fusions were first arrayed into 96-well format (seven plates), transferred into two 1536-well plates (in quadruplicate) and screened against the pre-screened 135 Cub fusions. Mating reactions were performed by robotic pinning of individual Cub-strains onto the arrayed Nub strains (1536-well format) on YPD plates. Twenty-four hours later, cells were transferred to DS media to select for diploid cells (a representative plate is shown in Figure 2A). After 3 days of growth at 28°C, cells were pinned robotically to IS media. Three methionine levels, which regulate the MET25 promoter driving the Cub fusions, were used to vary the stringency of the screening conditions (Figure 2B). Putative interactors were identified by visual inspection, based on complementation of the His3 marker (histidine auxotrophy) and the Ade2 marker (absence of red pigments). The interaction screen identified 386 pPPIs. One NubG clone, AT5G06100, a MYB transcription factor (consistent with being a “false positive” due to “auto-activation”; 86 interactors) was not included in further experiments. The remaining positive interaction pairs were re-arrayed and retested for complementation of his3, for red color of ade2 mutants and the third marker, LacZ. 88.9% of the initial positive pairs (343 out of 386) showed a positive read-out for the all three markers (representative plate, Figure 2C). Pairs that showed a positive read-out for all three markers (Table S3 in Supplementary Material) were entered into Cytoscape 2.6.3 (Shannon et al., 2003) to visualize the network (Figure 3). Network analysis was performed using NetworkAnalyzer 2.6.1 (Assenov et al., 2008). The final network, AMPv2, consists of 179 proteins linked by 343 edges; the node degree distribution follows the power-law with r2 = 0.863, indicating that the network is scale-free (Figure 4). From this network, a set of five subnets was selected covering RLKs (Figure 5), the seven-transmembrane domain receptors (7TM receptors; Figure 6), calcium and phosphoinositide related proteins (Figure 7), and an AMT ammonium transporter (Figure 8). Hundred and forty-one LRR-RLKs out of the 223 encoded by the genome and 8 MLOs (Mildew resistance locus o) out of 15 were represented in this screen. The pPPIs identified in this screen are discussed in more detail below.

Figure 2

pPPIs detected in the AMPv2 screen. (A) Selection of diploid cells on DS media supplemented with adenine and histidine. (B) Successful mating (His3, Ade2 positive) was selected on IS-0 media without met (right) and IS-500 with 500 μM met (left). (C) Representative image of the β-galactosidase assay performed on yeast grown on nitrocellulose-overlaid agar media. Positive pairs were picked by hand, re-arrayed and tested for the LacZ marker. Plate size for (A–C) was 125 mm × 85 mm.

Figure 3

AMPv2 pPPI network. Network of all pairs with a positive read-out for all three markers as visualized in Cytoscape. The network comprises 179 nodes and 343 edges. Symbol color corresponds to degree of interactions (color code on the left). Symbol shape indicates protein function (code at the top).

Figure 4

Representation of Nub and Cub clones in AMPv2 network and frequency distribution of the resulting network. (A) The mbSUS screen was performed with 490 genes of which 443 were mobilized into the Nub destination vector and 204 into the Cub destination vector (157 genes were common to both pools). The pPPI network between 179 nodes yielded 343 edges. (B) The power law frequency distribution of the network was obtained with the “NetworkAnalysis” tool of Cytoscape.

Figure 5

Sub-network of receptor kinase interactions (RLK). Labeling as in Figure 3. The blue borders surrounding 31 individual RLKs indicate evidence for phosphorylated sites according to PhosPhAt 3.0 (Durek et al., 2010).

Figure 6

Sub-network of 7TM receptors. (A) Members of the MLO receptor family are displayed as a sub-network. (B) GCR1 sub-network (AT1G48270). Labeling as in Figure 5.

Figure 7

Sub-network of calcium ATPases, phospholipase C2 and aquaporins. (A) ACA12 sub-network (AT3G63380), (B) PLC2 sub-network (AT3G08510), and (C) aquaporin sub-network. Labeling as in Figure 5.

Figure 8

Ammonium transporter sub-network. (A) AMT1;1 sub-network (AT4G13510). Labeling as in Figure 5. (B) Independent analysis of AMT1;1 interaction with two RLKs identified in AMPv2 using the split luciferase assay in Arabidopsis leaf protoplasts. Red asterisks in the diagram to the right indicate the position of the split luciferase fusions.

Receptor-like kinase potentially interacting with the ammonium transport AMT1

The screen identified three RLKs (AT2G28990, AT1G72180, and AT5G59650) that gave a positive read-out with the ammonium transporter AMT1;1 or a mutant form of AMT1;1 (T460A; Figure 8A). Two of the AMT/RLK interactions found in AMPv2 screen were tested independently in Arabidopsis leaf protoplasts using a split luciferase system (Fujiwara et al., 1993; Kato et al., 2010). An interaction test between an R-SNARE (VAMP721) with a Qa-SNARE (SYP122) was used as a positive control (Figure 8B). The luciferase activity obtained from coexpression of SYP122 with PHT1;4 (Arabidopsis phosphate transporter AT2G38940 that does not appear to interact with SYP122) served as a negative control (Figure 8B). Coexpression of SYP122 with VAMP721 produced 45-fold higher luciferase activity compared to coexpression of SYP122 with PHT1;4. Previous work had shown that AMT1 functions as an oligomer and that a specific interaction of the cytosolic C-terminus of each of the three subunits is required for transporter activity (Obrdlik et al., 2004; Loqué et al., 2007). The ability of AMT1;1 to oligomerize was confirmed here since coexpression of an AMT1;1-NLuc and a corresponding CLuc fusion reconstituted significant luciferase activity. The two RLKs (AT2G28990 and AT1G72180) were fused C-terminally with CLuc. An interaction test with AMT1;1-NLuc yielded significant reconstitution of luciferase activity relative to the negative control (Figure 8B). This result was confirmed in a reciprocal test of AMT1;1-CLuc with the RLKs as C-terminal NLuc fusions (Figure 8B). Three independent experiments yielded comparable results. Thus the split luciferase system was able to confirm the interaction of AMTs with each other as well as a potential interaction of AMT1 with two RLKs.

Discussion

Here we describe the cloning of 3,852 Arabidopsis ORFs and show the successful establishment of a high throughput membrane protein/signaling protein interaction screen using the mbSUS system in 1536-well format. To establish a high throughput screen for pPPIs, 490 ORFs were transferred into the yeast mbSUS vectors and used for an interaction screen (90,370 potential edges). A pre-screen aimed at putative “false positive” and “false negative” interactors eliminated 33.8% of the Cub fusions. Out of the remaining 59,805 potential PPIs between 443 Nub fusions and 135 Cub fusions, 343 pairs (0.57%) were positive in the interaction screen. Surprisingly, despite the small number of actual interactions, given the total theoretical number of possible interactions among 490 ORFs, a network (AMPv2) with small world properties was identified.

LRR-RLKs are expected to play important roles in the plasticity of plant development and in interactions with the environment. Due to the large representation of LRR-RLKs here, this work may provide new leads for characterizing the physiological role of RLKs in plants. Interestingly, the ammonium transporter AMT1;1 gave reproducible interactions with three RLKs. Two RLKs (AT1G72180 and AT2G28990) were further tested and independently confirmed by the split luciferase assay in Arabidopsis leaf protoplasts, supporting the reliability of the observed pPPIs. The cloning of a large fraction of the plant membrane/signaling proteins as provided here, together with the successful establishment of this screening platform will provide the foundation for a genomic scale interaction study of protein–protein interactions. The present platform enables a throughput of approx 75,000 pPPIs to be tested per week and provides a basis for a systematic binary analysis of interactions between and among membrane and signaling proteins in a multicellular organism.

Potential and limitations of the large-scale interactome screen using mbSUS

As outlined in the introduction, a wide spectrum of studies support the reliability of the split ubiquitin system as a discovery tool for protein interactions. Nevertheless, the system is based on overexpression and thus may produce “false positive” results. Since yeast is a heterologous system, the plant proteins may either be targeted incorrectly, or accumulate in the endoplasmic reticulum (ER) or in carmellae (Villalba et al., 1992). The heterologous system can potentially lead to accumulation of proteins either at unphysiologically high levels or colocalization of proteins that localize to different compartments in the plant. This can lead to an increased probability of random collisions producing “false positive” read-outs (Dünnwald et al., 1999). Moreover, some interactions may require post-translational modifications, e.g., phosphorylation as a prerequisite for 14-3-3 binding, a modification that may be lacking in yeast. Our protocol eliminates “autoactivators” and proteins with insufficient expression levels in a pre-screen (cf. 3.3), providing a systematic source of “false negatives.” A major advantage of large-scale screens will be that the data obtained are directly comparable, and that proteins with a propensity of giving “false positive” or “false negative” read-outs can be identified. The use of orthogonal interaction tests may be a means of testing the reliability of the mbSUS screens (Lalonde et al., 2008).

Interactions with receptor-like kinases

Membrane-localized RLKs recognize specific ligands via their N-terminal extracytosolic domain and transmit ligand-binding status across the transmembrane domain to activate the intracellular kinase. Typically, intermolecular autophosphorylation within dimers leads to subsequent transduction of the signal to downstream targets (Pawson and Nash, 2000). Plant genomes are unique in that they contain ∼500 genes encoding transmembrane RLKs (Arabidopsis 491, rice 807; Krupa et al., 2006). Plasma membrane-localized RLKs have been shown to function in development, as well as in hormone and defense signaling (De Smet et al., 2009; Tör et al., 2009; Zhao, 2009).

Although RLKs are thought to form heteromers, only a few protein-interactions of RLKs have been characterized in detail in plants. One of the best-characterized pairs is the brassinolide receptor BRI1 (AT4G39400) and its partner the BRI1-associated kinase BAK1 (AT4G33430). Both genes were present in our screen; however both were only tested as Nub fusions.

Here, the potential of 141 RLKs to interact with each other, or with other membrane and signaling proteins, was analyzed to gain insights into potential pathways in which these RLKs might be involved (e.g., RLKs potentially involved in regulating ammonium transporters; Figure 8). To analyze the pPPIs of RLKs in more detail, all nodes corresponding to RLKs were extracted from the network (Figure 5). The RLK sub-network comprises 80 RLKs forming 126 pPPIs. On average, RLKs share 3.0 edges with other proteins. The majority of RLKs (63) has four or fewer edges, whereas the LRR-RLK AT2G41820 had 21 edges (Figure 3). Three RLKs showed potential homo-oligomerization (AT5G48380, AT2G41820, and AT1G11130), 63 had heteromeric edges with other RLKs, and 13 RLKs potentially interacted exclusively with one other protein. Interestingly, 31 of the RLKs detected here had previously been found to be phosphorylated based on proteomics analyses (PhosPhAt 3.0; Durek et al., 2010), indicating that they might be substrates for the respective RLKs (blue dotted borders, Figure 5).

Interactions with 7TM receptors

The Arabidopsis genome encodes a single known Regulator of G-protein Signaling (RGS1, AT3G26090), single canonical Gα and Gβ subunits, and two Gγ subunits. The AMPv2 ORF set contained eighteen 7TM receptors: GCR1, which is a candidate G-protein-coupled receptor (GPCR) based on its sequence homology to Dictyostelium GPCRs, one “unknown,” two transferases, six members of the MtN3 family, and eight MLO receptors (Devoto et al., 2003). Of these 18 genes, 8 were present in the network, including GCR1 and five MLOs (for the AGI number, see Figure 6). All MLO interactions and their first neighbors were extracted (Figure 6A). Most striking is the finding that most MLOs are highly connected and are part of the same sub-network. MLO proteins have seven transmembrane domains (Devoto et al., 1999), but despite the topology reminiscent of GPCRs, previous reports had not shown an interaction of MLOs with heterotrimeric G-proteins (Kim et al., 2002; Chen et al., 2009). Here we show that MLO2 (AT1G11310) potentially interacts with GCR1 (AT1G48260, Figure 6A). Very few direct targets of G-protein signaling have been identified in plants, and our screen provides additional candidate targets, because we also observe potential interactions of GCR1 with the small G protein ROP5 (AT4G35950), the cell wall arabinogalactan protein AGP27 (AT3G06360), and two transporters (NRAMP2, AT1G47240, and ammonium transporter AMT1;1, AT4G13510; Figure 6B). It is interesting that three of these interactors (MLO2, NRAMP2, and AMT1;1) are also interaction partners of GCR1 (Figure 6B), suggesting the potential formation of a G-protein signaling complex.

Calcium and phosphoinositides

The second messenger calcium plays a central role in signal transduction. Calcium is involved in hormonal signaling, stress response signaling, plant cell development and host–pathogen interactions (McAinsh and Pittman, 2009). Changes in cytoplasmic calcium are generated by calcium channels and pumps, which mediate import and efflux of calcium across the plasma membrane, and release or storage from/in internal stores such as the ER or the vacuole (Mäser et al., 2001; Kim et al., 2010). In Arabidopsis, autoinhibited calcium ATPases (ACA) are encoded by 14 genes, and ACA proteins have been found at the plasma membrane and the ER (Boursiac and Harper, 2007). The plasma membrane ACA12 (AT3G63380, Boursiac and Harper, 2007) showed putative interactions with the receptor-like kinase RPK2|TOAD2 (AT3G02130), ROP5 (AT4G35950), and phospholipase PLC2 (AT3G08510, Figure 7A). All four genes in this sub-network are expressed in open flowers (Schmid et al., 2005). Thus the observed interactions can occur in planta supporting a role in flower-related processes. PLC2, a phosphoinositide-specific phospholipase (PI-PLC) catalyzes the hydrolysis of phosphoinositide [PI(4,5)P2] to inositol trisphosphate [Ins(1,4,5)P3] and diacylglycerol (DAG), which function as second messengers. Ins(1,4,5)P3 and DAG are known to stimulate release of Ca2+ from intracellular stores and Ca2+ entry in many eukaryotic cells (Taylor, 2002). Arabidopsis encodes seven PI-PLCs; here only PLC2 (AT3G08510) was tested. PLC2 has the potential to interact with the Ca2+-ATPase ACA12, 10 RLKs, three transporters, four unknown proteins, MLO2, a tetraspanin, and with CDP-cytidine synthase (Figure 7B).

The putative interactions between a small GTPase (ROP), phospholipase (PLC), a calcium pump (ACA), and a RLK (Figure 7B) may be linked functionally because both calcium and phosphoinositides play important roles in polarized growth. Moreover, all four proteins probably localize to the plasma membrane. PLC2 activity maintains polarity of RAC-ROP signaling, thereby affecting actin organization and membrane traffic. RAC-ROP signaling and Ins(1,4,5)P3 affect calcium, and calcium can affect PLC2. ACA12 appears to be the central connector of the four pPPIs. However it is important to note that not all possible binary interactions were tested. One may thus hypothesize that ACA12 can exist in a complex with the phosphoinositide metabolism enzymes and is involved in calcium flux across the plasma membrane. The interaction with the RLK RPK2|TOAD2 may function as a regulator of the participating enzyme activity, e.g., ACA12.

Aquaporin interactions

In addition to a crucial role in signaling, the plasma membrane is key to uptake and efflux of water, ions and metabolites. Many transporters have been shown to exist as oligomers and interact with regulators (Loqué et al., 2007). To identify clusters of interacting proteins that might form a complex or co-functioning module, we clustered the network with the MCL algorithm (Enright et al., 2002), resulting in 27 clusters. We performed GO enrichment analyses on the clusters containing at least three genes against all genes in the network as background using Ontologizer (Bauer et al., 2008) and found only one cluster that was enriched in relation to water deprivation (cluster 17). This group comprises three aquaporins, one cation/ion exchanger (CHX9, AT5G22910), and one protein of unknown function (Figure 7C). The unknown protein is predicted to contain four transmembrane spanning domains and has a DUF679 motif (domain of unknown function 679). Aquaporins of PIP1 and PIP2 families are known to homo- and hetero-oligomerize (Zelazny et al., 2007). We find interactions between PIP2;2 (AT2G37170) and PIP2;1 (AT3G53420), and between PIP2;1 and PIP2;5 (AT3G54820), as well as homomerization of PIP2;1. No other homodimer of PIPs could have been found because only PIP2;1 was present in the screen as both a Cub and Nub fusion.

Ammonium transporter interactions

Plant ammonium transporters function as trimers. The activity of the three subunits in the complex is allosterically regulated via trans-activation by the cytosolic C-terminus (Loqué et al., 2007). An increase of the extracellular ammonium concentrations was shown to trigger phosphorylation of a critical threonine (T460) in the trans-activation domain of the C-terminus (Loqué et al., 2007; Lanquar et al., 2009). Besides the regulatory site T460 in the C-terminus of AMT1;1 (AT4G13510), additional phosphorylation sites have been identified, suggesting that a suite of protein kinases regulates AMT activity (Lanquar et al., 2009). Because AMT1;1 phosphorylation is triggered by extracellular ammonium, RLKs might be prime candidate kinases involved in ammonium perception.

Here, AMT1;1 and a mutant form (T460A) were found to interact with three different RLKs (AT2G28990, AT1G72180, and AT5G59650) (Figure 8A). Two of the AMT/RLK interactions found in our mbSUS screen were tested independently in plant cells using a split luciferase system (Fujiwara et al., 1993; Kato et al., 2010). As may have been anticipated from the fact that AMTs oligomerize, coexpression of an AMT1;1-NLuc and a CLuc fusion reconstituted significant luciferase activity, at a level similar to a positive control interaction of an R-SNARE (VAMP721) with a Qa-SNARE (SYP122). Coexpression of both RLKs fused C-terminally with CLuc (AT2G28990 and AT1G72180) with AMT1;1-NLuc yielded significant luciferase reconstitution relative to the negative control. This result was confirmed in a reciprocal test of AMT1;1-CLuc with the RLKs as C-terminal NLuc fusions. Further analyses will be required to determine whether AMT1 is a substrate for the RLKs, and if so, which residues in AMT1 are phosphorylated.

In addition to the RLK interactions, the AMPv2 screen suggested that AMT may interact with three other transporters, four “unknowns,” and two signaling related proteins: GCR1 (AT1G48270) and CDP-cytidine synthase (Figure 8A).

Comparison of AMPv2 to soluble protein interaction data sets

We did not find any of the pPPIs in the networks of Popescu et al. (2009), Intact (1,692 proteins and 3,424 interactions), BioGrid (release 62; 1,645 proteins and 3,291 interactions) or TAIR (1,334 proteins and 2,178 interactions; as of June 2010). This may not be surprising because only 83 out of the 490 ORFs (as Nub and/or Cub) tested here were common to all the databases listed above. This difference may in part be due to the fact that AMPv2 focused on membrane and signaling proteins, whereas 89% of the ORFs in the above four databases are predicted to be soluble. The public databases contained eight putative interactions that could not be confirmed in this screen (Table 1). Furthermore, well-documented interactions such as between BRI1 and BAK1, BRI1 with the 14-3-3λ protein were not tested because all three were only present as Nub fusions (Gampala et al., 2007; Schulze et al., 2010). Further work will be required to determine whether these represent “false negatives” in AMPv2.

Table 1

Published interactions (Popescu, Intact, BioGrid, or TAIR) that were not found in our screen.

AT4G33430 (BAK1|SERK3) = AT1G17750 (RLK/Pelle)

AT4G33430 (BAK1|SERK3) = AT5G48380 (BIR1)

AT1G34210 (SERK2) = AT1G34210 (SERK2)

AT1G75820 (CLV1) = AT1G65380 (CLV2)

AT3G02520 (14-3-3 Nu) = AT1G71830 (SERK1)

AT4G39400 (Bri1|Bin1) = AT1G71830 (SERK1)

AT5G10450 (14-3-3 lambda) = AT1G71830 (SERK1)

AT1G75820 (CLV1) = AT1G75820 (CLV1)

Conclusion

Over 3,852 Gateway-compatible ORFs, encoding typically difficult to clone membrane protein genes, have been generated and made publicly available through ABRC. We have identified 343 potential protein–protein interactions among 179 proteins including 134 membrane proteins. Membrane protein and signaling protein interactions provide a basis for analyses into the physiological functions of plant membrane proteins. A larger scale screen, currently in progress, with a binary matrix using >2,000 proteins should reveal more information on the prevalence of “sticky” proteins and the overall topology of a membrane-based sub-interactome, and is expected to provide a wealth of leads for further biochemical, physiological, and bioinformatic analyses into plant membrane protein regulation and function.

Contributions

Sylvie Lalonde, Wolf B. Frommer, Sarah M. Assmann, June M. Kwak, Seung Y. Rhee, Julian I. Schroeder conceived the study; Sylvie Lalonde and Wolf B. Frommer supervised the project and wrote the manuscript; Guillaume Pilot, Rainer Schwacke, and Réjane Pratelli generated the target list of candidates; Antoinette Sero, Réjane Pratelli, Guillaume Pilot, Maria I. Sardi, Saman A. Parsa, Do-Young Kim, Biswa R. Acharya, Erica V. Stein, Heng-Chen Hu, Florent Villiers, Kouji Takeda, Yong S. Han, Yingzhen Yang generated the ORF/mbSUS clones, Antoinette Sero, Do-Young Kim, Maria I. Sardi, and Saman A. Parsa carried out the mbSUS screen; Sylvie Lalonde, Rainer Schwacke, Jin Chen, and Seung Y. Rhee performed the bioinformatics analyses, William Chiang and Naohiro Kato carried out the split luciferase assays; and Dominique Loqué established the 96-well yeast transformation.

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.

Supplementary Material

The Supplementary Material for this article can be found online at http://www.frontiersin.org/physiology/plantphysiology/paper/10.3389/fphys.2010.00024/