Tetramer formation in Arabidopsis MADS domain proteins: analysis of a protein-protein interaction network
© Espinosa-Soto et al.; licensee BioMed Central Ltd. 2014
Received: 21 September 2013
Accepted: 2 January 2014
Published: 27 January 2014
MADS domain proteins are transcription factors that coordinate several important developmental processes in plants. These proteins interact with other MADS domain proteins to form dimers, and it has been proposed that they are able to associate as tetrameric complexes that regulate transcription of target genes. Whether the formation of functional tetramers is a widespread property of plant MADS domain proteins, or it is specific to few of these transcriptional regulators remains unclear.
We analyzed the structure of the network of physical interactions among MADS domain proteins in Arabidopsis thaliana. We determined the abundance of subgraphs that represent the connection pattern expected for a MADS domain protein heterotetramer. These subgraphs were significantly more abundant in the MADS domain protein interaction network than in randomized analogous networks. Importantly, these subgraphs are not significantly frequent in a protein interaction network of TCP plant transcription factors, when compared to expectation by chance. In addition, we found that MADS domain proteins in tetramer-like subgraphs are more likely to be expressed jointly than proteins in other subgraphs. This effect is mainly due to proteins in the monophyletic MIKC clade, as there is no association between tetramer-like subgraphs and co-expression for proteins outside this clade.
Our results support that the tendency to form functional tetramers is widespread in the MADS domain protein-protein interaction network. Our observations also suggest that this trend is prevalent, or perhaps exclusive, for proteins in the MIKC clade. Because it is possible to retrodict several experimental results from our analyses, our work can be an important aid to make new predictions and facilitates experimental research on plant MADS domain proteins.
KeywordsMADS domain proteins Protein-protein interaction network Transcription factors Tetramers Subgraph abundance Arabidopsis thaliana
MADS domain proteins comprise a family of eukaryotic transcription factors that are very important for plant and animal development [1–4]. In plants, MADS domain transcription factors have diversified extensively after multiple independent duplication events, and they now play important roles in many aspects of development [2, 5, 6]. They are involved in many processes, such as floral timing, floral organ identity specification, gametophyte, fruit and root development, among others [2, 7–10]. The importance of plant MADS domain proteins is reflected in their evolutionary history: some events in the history of this protein family correlate with major evolutionary modifications in plant morphogenesis. Moreover, there is evidence of positive selection having acted at particular aminoacids of these plant proteins at the times of these evolutionary events [5, 6, 11].
MADS domain proteins are grouped into two distinct monophyletic clades, type I and type II, that split before the separation of animal and plant lineages . Each of these two clades includes smaller sub-clades . Within type II proteins, the plant-specific MIKC monophyletic clade stands out. MIKC proteins have a distinctive structure in which four different domains have been recognized. The first is the MADS (M) domain, that is widely conserved, common to all MADS domain proteins, and participates in binding to DNA. MIKC proteins also have a poorly conserved Intervening (I) domain, which is important for the specificity of dimer formation. Next we find a Keratin-like (K) domain that is involved in protein-protein interactions. Finally, there is a C-terminal (C) domain that may carry transcriptional activation domains. Most MADS domain proteins that are functionally characterized belong to the MIKC clade [6, 7, 9].
Importantly, MADS domain proteins do not exert their functions as monomers, but rather they form multimeric protein complexes with other MADS domain proteins [13–16]. A comprehensive understanding of how MADS domain proteins interact with multiple molecules in a concerted manner would provide a significant advance in understanding the molecular basis of the wide range of developmental mechanisms in which these proteins participate.
At least some MADS domain proteins build tetramers with other MADS domain proteins to regulate transcription of their target genes, as proposed originally by the ‘quartet model’ [17, 18]. The tetramer binds two different sites in the regulatory sequence of target genes and forces looping of DNA stretches, consequently regulating the transcription of the associated genes [19–21]. For instance, the floral homeotic MADS domain proteins APETALA3 (AP3) and PISTILLATA (PI) bind two other MADS domain proteins, APETALA1 (AP1) and SEPALLATA3 (SEP3), to determine petal cell identity in the flower of the thale cress Arabidopsis thaliana[22–24]. It is noteworthy that, so far, all experimentally validated tetramers comprise proteins from the MIKC clade exclusively. It is still unknown if tetramer formation is a widespread property of many plant MADS domain proteins, or if it is specific to a few of them.
Here we address if tetramer formation is widespread among MADS domain proteins in the plant model system Arabidopsis thaliana. We do so by analyzing the structure of a network of physical interactions among Arabidopsis MADS domain proteins. We built the network with data from high-throughput yeast two and three-hybrid studies that have uncovered a near-complete map of interactions among MADS domain proteins in Arabidopsis[15, 25]. We supplemented these data with additional information from various sources and analyzed the resulting protein-protein interaction network.
First, we studied the abundance of different subgraphs in the network. Each subgraph represents a pattern of connections among proteins in a set, in a manner independent of the identity of each protein. The abundance of particular subgraphs in different kinds of biological networks may indicate that natural selection has favored such a connection pattern across evolution. This would be the case, for example, when elements connected in that manner jointly perform a function that confers an advantage to an organism. Consider the case of biological networks where transmission of information is important, such as neuron nets or signal transduction networks. In this kind of networks a particular subgraph, known as feedforward loop, appears with a frequency higher than expected by chance [26–28]. Such a connection pattern is seemingly useful to perform distinct signal-processing tasks [29, 30].
If the ability to build functional tetramers is a widespread property of plant MADS domain transcription factors, we expect that subgraphs compatible with tetramer formation will appear more often than expected by chance in the MADS-domain protein-protein interaction network. Experimental research has provided valuable information regarding how some MADS domain proteins bind other such proteins to form tetramers. These studies thus hint on the properties of the subgraphs that could represent this kind of protein complexes. Apparently, each MADS domain protein binds two other proteins in the tetramer through different protein domains. A MADS domain protein binds, through one of its domains, a partner to form a dimer. It also establishes another interaction, through a different domain, to another MADS domain protein in a second dimer. For example, consider the experimental analyses of the roles of different parts of AGAMOUS (AG), another floral homeotic MADS domain protein. Such studies have shown that two adjacent regions of AG, the MADS domain and I region, are sufficient for DNA binding and for the formation of either homodimers  or heterodimers with the MADS domain protein SEPALLATA1 (SEP1) . However, these two protein regions are not sufficient for a functional AG protein. The C region is also required, most likely for binding other proteins to form higher-order protein complexes . Indeed, Fan and collaborators found that the K and C domains of AG can bind other MADS domain proteins, such as AGAMOUSLIKE-6 (AGL6), SEP1, SEPALLATA2 or SEPALLATA3 (SEP1-3) . Another example concerns the three predicted α-helices (K1, K2, and K3) in the K domain of the PI protein. While K1 is more important for the PI-AP3 interaction, K3 is more important for PI’s interaction with SEP3 . Finally, Melzer et al. have shown that a part of the K domain in the MADS domain protein SEPALLATA3 (SEP3) is not essential for the formation of SEP3 homodimers . However, the same SEP3 K domain is required to mediate interactions between a SEP3 homodimer and the APETALA3-PISTILLATA heterodimer , or another SEP3 homodimer .
The experimental evidence mentioned in the preceding paragraph suggests that functional MADS domain tetramers require that each protein binds, through different protein regions, two other proteins in the tetramer. Thus, we paid special attention in our analysis to those four-node subgraphs where each node interacts with two other nodes in the subgraph. We found that the network of protein-protein interactions among MADS domain proteins in Arabidopsis thaliana contains significantly more such tetramer-compatible subgraphs than expected by chance. In contrast, tetramer-compatible subgraphs are not significantly abundant in an analogous network for a different family of plant transcription factors, the TCP protein family [35, 36]. Moreover, MADS domain proteins in a tetramer-compatible subgraph are more likely to be co-expressed than proteins in other subgraphs. In addition, we show that this association between tetramer-like subgraphs and joint expression is evident only for proteins in the MIKC clade. We also analyzed several aspects of the organization of the whole network of interactions among MADS domain proteins. Taken together, our results suggest that the formation of tetramers is a widespread property of Arabidopsis MIKC MADS domain proteins. We finally discuss how our work may be an important aid towards the design of new experiments to uncover functional tetramers in planta.
A network of interactions among MADS domain proteins
The MADS network comprises 79 nodes, including four proteins that result from alternative splicing (ABS-2, SEP4-2, AGL74-2 and AGL74-N). There are 312 protein-protein interactions in the network. The most highly connected node (AGL74-N) has 34 interactions, and the median number of interactions per protein is 6.
Some connection patterns appear with a frequency higher than expected by chance in the MADS network
Abundance of frequent subgraphs in the MADS network as compared to that of randomized networks
Abundance of connection patterns in the MADS network
Throughout this work we have assumed that most positive results in yeast-three hybrid assays reflect pairwise protein-protein interactions. Consider an assay where proteins X and Z function as bait and prey, respectively, and that produces a positive result only in the presence of a third protein Y. In this case, we assume pairwise interactions between X and Y on the one hand, and between Y and Z, on the other (see Methods). Alternatively, such a positive result could indicate that these proteins bind only in the presence of all three proteins. We have dismissed the latter alternative, as a simplifying assumption. Nonetheless, it is noteworthy that more than two thirds (58/84) of the pairwise interactions that we derived from a large scale three-hybrid study  were also supported by two-hybrid data in the study by de Folter and collaborators . Moreover, if the formation of higher-order complexes, as those inferred from yeast three-hybrid studies, depended on non-pairwise interactions between a dimer and a third protein there would be specific implications that are not fulfilled. Namely, if the aforementioned hypothetical proteins X and Y form the dimer that binds in a non-pairwise manner to protein Z, we would also recover the interaction between X and Y from yeast two-hybrid studies, but not the interaction between Y and Z. In other words, a scenario where non-pairwise interactions are prevalent in higher-order protein complexes implies that, for most positive yeast three-hybrid results, only one pairwise interaction is also supported by yeast two-hybrid data. In contrast, our assumption predicts that for most positive yeast three-hybrid results, there are two interactions supported by yeast two-hybrid assays. We thus assessed how many of the sets of three proteins that produce a positive yeast three-hybrid in a large-scale study  reflect these patterns in the large-scale two-hybrid study . Of the 116 sets of three proteins that yield a three-hybrid positive result, 86 (74.1%) have two interactions also supported by two-hybrid data, consistently with our hypothesis. In contrast, 26 (22.4%) and four (3.4%) have only one or none interactions, respectively, supported by two-hybrid data.
Our observations suggest that most interactions inferred from yeast three-hybrid data are pairwise, as we assumed. Notwithstanding, we addressed whether our observations in Tables 1 and 2 depend on this assumption. We assembled a new network that includes only interactions from the aforementioned two-hybrid study . This network thus contains only interactions confirmed to be pairwise. We compared subgraph counts of this network to that of 104 of randomized networks that preserved the number of nodes, connections, and connectivity distribution. The results are qualitatively the same to those that we obtained from the analysis of the MADS network (compare Tables S2 and S3 in Additional file 1 to Tables 1 and 2, respectively). Thus, the statistically significant abundance of tetramer-like connection patterns in the MADS network does not depend on the assumption that interactions derived from three-hybrid data are pairwise. In conclusion, our results on the analysis of the MADS network are robust to the removal of all the interactions derived from three-hybrid data.
Tetramer-like subgraphs are not abundant in interaction networks of the TCP plant protein family
It may be that the abundance of tetramer-like connection patterns is not specific of MADS domain proteins, but rather a general trend in plant protein interaction networks. If this is the case, the abundance of tetramer-like subgraphs would be more easily explained by factors common to the evolution of plant protein interaction networks rather than by functional constraints favoring MADS domain proteins to form tetrameric complexes. Secondly, the observed abundance of tetramer-like connection patterns could be due to the nature of the yeast-based protein-protein interaction assay. In this case, similar trends and distribution of subgraphs are expected in a large-scale yeast n-hybrid data set for another protein family.
Abundance of tetramer-like connection patterns in the TCP network as compared to that of randomized networks
The numbers of proteins and interactions in the TCP and MADS networks are so contrasting that one may think that network size could underlie the differences in subgraph abundance that we found between these two networks. This is unlikely, since we compare each network with randomized versions with the same number of nodes, connections, and connectivity distribution. Notwithstanding, we addressed this possibility by sampling randomly 20-node sub-networks from the MADS network, keeping only those sub-networks that had the same connectivity distribution as the TCP network. In that manner, from a sample of more than 106 sets of MADS domain proteins, we obtained 50 sub-networks. These sub-networks contained experimentally validated interactions among MADS domain proteins, but also the same number of nodes and connectivity distribution as the TCP network. The 50 sub-networks taken from the MADS network have an average of 216.22 tetramer-like connection patterns (median = 216.5; std. dev. = 19.37). In contrast, the 104 randomized networks have a mean number of only 161.1 tetramer-like connection patterns (median = 160; std. dev. = 22.9). The difference is statistically significant, according to a Mann-Whitney U test (U = 11.39; p < 2.6 × 10-308). In sum, tetramer-like connection patterns are significantly abundant in partitions of the MADS network with the same connectivity distribution as the TCP network, but not in the TCP network.
Co-expression is frequent for MADS domain proteins in tetramer-like connection patterns
In order to form functional tetrameric complexes, MADS domain proteins in a tetramer-like connection patterns must be expressed jointly in vivo. To address whether this is usually the case, and also whether there is a difference in the degree of co-expression with proteins in other four-node subgraphs, we considered gene expression data reported in the developmental data set of the AtGenExpress project  (see Methods).
Abundance of tetramer-like subgraphs in the MADS exp network as compared to that of randomized networks
All tetramer-like subgraphs
Co-expressed tetramer-like subgraphs
We repeated the analyses described in the previous paragraph, but using the At-TAX data base , instead of the AtGenExpress data set. The At-TAX database resulted from a whole genome tiling array for Arabidopsis. It contains expression data for more genes, but for fewer tissues, than the AtGenExpress data set. We call the network that resulted after elimination of proteins lacking expression data in this set the MADS At-TAX network. It contains 74 nodes and 254 edges. This network also has a significant abundance of tetramer-like connection patterns, when compared to randomized networks with the same connectivity distribution. Again, this result holds after dismissing those subgraphs including proteins whose encoding mRNA never coincides in at least one tissue. Despite the At-TAX database contains expression data for more genes encoding MADS domain proteins, it is less suitable than the AtGenExpress data set for the analysis of co-expression of different proteins. Given that the At-TAX database includes information for only ten different wild-type plant tissues, there are many pairs of co-expression false negatives. Consider, for instance, those sets of four proteins connected according to a tetramer-like connection pattern. The MADS exp network contains 1351 such sets and the MADS At-TAX contains 1849 of them. However, after dismissing those sets where evidence of co-expression is lacking, the MADS exp network contains 1134 but the MADS At-TAX includes only 574 of them. Thus, we performed the remaining analyses only using the AtGenExpress data set and the MADS exp network.
Recent experimental research suggests that SEP1-3 proteins have an important role facilitating the formation of MADS domain protein tetramers [20, 21, 25]. The composition of the list of prospective tetramers supports the prominence of these proteins for the assembly of tetrameric complexes: SEP1-3 proteins appear 840 times in the list of prospective tetramers. While SEP3 and SEP1 appear each 385 times, SEP2 appears only 70 times. It is noteworthy that other MADS domain proteins also occur many times in the list of prospective tetramers. AGL21 appears 376 times, SOC1 457 times, AGL6 290 times and AGL24 278 times. SEP1-3 proteins perform their functions in the floral meristem . It may well be that other proteins that appear recurrently in the list of prospective tetramers play a SEP-like role in the formation of tetramers outside the flower.
Proteins in prospective tetramers are mostly MIKC MADS domain proteins
Until now, all experimentally validated MADS domain protein tetramers consist exclusively of proteins in the MIKC clade [20–24]. The K domain appears exclusively in MIKC proteins [4, 12] and, importantly, it has been implicated in several protein-protein interactions [14, 20, 21, 33, 34, 43]. For these reasons, it may well be that the formation of MADS domain protein tetramers occurs mainly for MIKC proteins. We thus addressed whether this could be the case.
Tetramer-like subgraphs composed of MIKC proteins exclusively are more likely to be jointly expressed
Not co-expressed tetramer-like subgraphs (T)
Co-expressed tetramer-like subgraphs (C)
Although our results support the formation of functional tetramers for MIKC proteins but not for non-MIKC proteins, two confounding factors must be first taken into account. The first is that the AtGenExpress data set includes expression data for a larger fraction of MIKC proteins when compared to non-MIKC proteins. The MADS network comprises 35 MIKC proteins and 44 non-MIKC proteins. In contrast, the MADS exp network, in which we discarded proteins encoded by genes lacking expression data, includes 32 MIKC but only 23 non-MIKC proteins. We thus performed an additional analysis to control for this factor. In this analysis, we first sampled randomly 15 non-MIKC and 15 MIKC proteins from the MADS exp network. The AtGenExpress data set contained expression information for all 30 genes encoding for such proteins. We next counted the number of sets of four of these 30 proteins where: i) the four proteins interacted according to a tetramer-like connection pattern, and ii) the genes encoding for such four proteins were co-expressed in at least one tissue. We performed 104 independent samplings. Even though the number of MIKC and non-MIKC proteins was the same in each sample, we found immense differences in the number of MIKC and non-MIKC prospective tetramers. On average, in each sample we found 0.6 prospective tetramers composed exclusively of non-MIKC proteins (median = 0, std. dev. = 0.94). In contrast, we found a mean number of 49.1 prospective tetramers comprising only MIKC proteins (median = 37, std. dev. = 44.14). We also counted prospective tetramers comprising both MIKC and non-MIKC. In this case we found only an average of 0.36 (median = 0, std. dev. = 0.64). The difference between the abundance of MIKC-only and non-MIKC-only prospective tetramers is highly significant, according to a Mann-Whitney U test (U = 117.5; p < 2.6 × 10-308).
Tetramer-like subgraphs comprising only non-MIKC proteins are not more likely to be co-expressed than other subgraphs
Not co-expressed (A)
Taken together, the results presented in this section do not support that non-MIKC proteins have an important role in the formation of functional tetramers. Rather, they support that MIKC proteins are predominantly, perhaps exclusively, involved in the formation of MADS domain protein tetramers.
Retrodiction of experimental results
Experimental studies have found several examples of MADS domain proteins that assemble into functional tetramers (see Table S4 in Additional file 1). We checked whether those experimentally supported tetramers appear in our list of prospective tetramers (Additional file 2).
Based on genetic data Gregis et al. infer that the proteins AP1, AGAMOUS-LIKE 24 (AGL24), SHORT VEGETATIVE PHASE (SVP) and CAULIFLOWER (CAL) participate in a functional tetramer . Indeed, this set of proteins appears in our list of prospective tetramers, thus supporting that this set of proteins may build a functional tetramer. The story is not as straightforward for the remaining experimentally-supported tetramers for the following reasons: First, several of these tetramers include the protein APETALA3 (AP3). In the data sets from which we built our network, AP3 only interacts with PISTILLATA (PI). Thus, in our study, AP3 can not take part in a tetramer-like subgraph, defined by each node having two interactions with different proteins. Second, one tetramer included a protein, AGAMOUS-LIKE 30 (AGL30), that has no interactions in our data set. Third, a protein appears twice in a validated tetramer. Thus, a three-node subgraph would represent one such tetramer. This occurs for AG, AP1 and SEP3. The latter is the protein that appears in most of the validated tetramers.
In order to partially overcome the limitations mentioned in the previous paragraph, we considered that SEP1 is a very close paralogue of SEP3 [4, 5]. Moreover, SEP1 and SEP3 are highly redundant , and they share many of the proteins with which they interact [15, 25]. Therefore, we assumed that SEP1 could replace one SEP3 instance in a tetramer. It is noteworthy that for all cases where the only limitation was that SEP3 appeared twice in a validated tetramer, we found a match in the list of prospective tetramers where SEP1 replaced one of the SEP3 instances.
Loss of function mutations of several different genes that encode for MADS domain proteins do not display obvious phenotypic alterations. The reason is that they are redundant with other MADS domain proteins so that phenotypic effects appear only in multiple mutant combinations . Identification of redundant genes may not be easy. Thus, we attempted to retrodict experimentally validated combinations of redundant genes as a means to evaluate the potential usefulness of our results as an aid for the design of new experiments.
First we searched in our list of prospective tetramers for all those that included a protein of our interest, X. For one of the X-including prospective tetramers, e.g. X-A-B-C, we scanned again the whole list of prospective tetramers to identify all those protein sets that share all proteins with the exception of X. In other words, we identified sets Y-A-B-C, Z-A-B-C, and so on. We preserved each such set only if the gene encoding for the protein that replaced X was co-expressed, in at least one tissue, with those genes of all four proteins in the original X-including prospective tetramer. Gene expression was determined from microarray expression data  (see Methods). We repeated the above for each of the remaining X-including prospective tetramers. We expected that proteins that replaced X most often in the prospective tetramers were more likely to be redundant with X. Indeed, we found that this was the case, as we describe below.
Pelaz and collaborators found that the SEP1-3 proteins are highly redundant . After following the algorithm described in the preceding paragraph, we found that SEP3 most often replaces SEP1 in prospective tetramers, both SEP1 and SEP3 substitute SEP2 more frequently, and SEP1 takes the place of SEP3 with greatest probability. Another relevant example concerns the protein CAL. Loss of function cal mutants do not display an obvious phenotype, unless combined with the ap1 mutation . In our analysis, there are three proteins that replace CAL most often in the same number of prospective tetramers: AP1, AGL6, and SEP1. Thus, also in this case we could identify the protein, AP1, that substitutes for the function of CAL. To our knowledge, there is no experimental evidence that loss of function agl6 or sep1 mutations could enhance the cal phenotype, but this remains as an open possibility and a prediction derived from our analyses. Finally, the SHATTERPROOF1/2 (SHP1/2) proteins are also redundant . Again, SHP2 is the protein that substitutes SHP1 most times in our list of prospective tetramers. These results support that our list of prospective tetramers and simple algorithms as that described above may be of value to infer relationships among MADS domain proteins and design experiments accordingly.
Homodimers in the formation of tetramers
So far, we have addressed how the patterns of interactions within sets of four different proteins may be indicative of MADS-domain protein heterotetramers. Hence, our analyses have left out of the picture those tetramers that incorporate two copies of the same protein. Our method, based on the analysis of pairwise interaction patterns, is unable to identify some of such tetramers. Consider for instance a hypothetical protein A that forms two kinds of dimers, one with protein B and other with protein C. The A molecule in the A-B dimer binds C in the A-C dimer, and the A molecule in the A-C dimer binds B in the other dimer to form a tetramer. In our network we would only observe interactions between A and B and between A and C (B-A-C), which is a connection pattern that appears in all non-disjoint 3 node sets. This argument is even stronger when we consider a tetramer in which two A-B dimers form a tetramer through interactions between A and B molecules in different dimers. Here, we would be unable to distinguish the interaction pattern associated to this tetramer (A-B) from any pairwise interaction in our network. Thus, our method is clearly not useful to associate a connection pattern to some of the possible tetramers that include two or more copies of a single protein. However, there are some cases where this association can be accomplished. These cases involve tetramers that incorporate homodimers.
As assumed for heterotetramers, we expect that each protein molecule in a tetramer that incorporates homodimers also binds two other protein molecules. The only difference to heterotetramers is that one such interaction must occur between the identical protein molecules that form the homodimer. Consider first a tetramer that contains a single homodimer. The corresponding 3-node subgraph is a triangle in which each node binds the other two nodes, and at least one node has a self interaction. For brevity, we will refer to 3-node connection patterns that embed this subgraph as three-node tetramer patterns (3NTP). Now consider a tetramer that incorporates two different homodimers. In this case, the subgraph corresponds to a pair of interacting nodes, each with self-interactions. We will call these subgraphs two-node tetramer patterns (2NTP). Subgraphs that include either of these interaction patterns are not significantly abundant in the MADS network when compared to randomized networks with the same connectivity distribution. However, it must be taken into account that few MADS domain proteins self-interact. Indeed, self-interactions occur significantly less often in the MADS network (6 self-interactions) than in randomized networks (a median of 10 self-interactions) (p = 0.045). It is thus relevant comparing the counts of these subgraphs in the MADS network to those in randomized networks with the same number of self-interactions as the MADS network. In this case, the number of occurrences of 2NTP’s and 3NTP’s is significantly higher in the MADS network. Out of 104 randomized networks with the same connectivity distribution and the same number of nodes, interactions and self-interactions as the MADS network, there are only two that have as many 2NTP’s as the MADS network (p = 0.0002) and none with as many 3NTP’s as the MADS network (p < 10-4). While the MADS network has 14 and 235 2NTP’s and 3NTP’s, respectively, randomized networks have a median of 8 2NTP’s and 160 3NTP’s. Hence, when we control for the occurrence of self-interactions, subgraphs associated to homodimer-including tetramers are overrepresented in the MADS network.
We also evaluated whether the mRNA coding for proteins in homodimer-including tetramers is co-expressed more often than for proteins in other subgraphs. The MADS exp network contains 14 2NTP’s, all of which have proteins with jointly expressed mRNA in at least one tissue. In contrast, there are 147 two-node subgraphs that are not 2NTP’s in the MADS network, but only for 131 of them the pair of proteins are jointly expressed, as inferred from RNA expression patterns. In other words, while proteins in 100% (14/14) 2NTP’s are jointly expressed, proteins in only 89% (131/147) of non-2NTP’s are co-expressed. The difference is not statistically significant according to a Pearson’s χ 2 test (χ 2 = 1.69; d f = 1; p = 0.193). However, this lack of significance may not be conclusive because of the low number of 2NTP’s in the MADS network.
As for 3NTP’s, there are 229 in the MADS network, and 213 of them are composed of proteins with mRNA co-expressed in one or more tissues. At the same time, the MADS network has 922 non-3NTP 3-node subgraphs, but only 752 of them are jointly expressed. Hence, the proteins in 93% of 3NTP’s are co-expressed in at least one tissue, while proteins in only 81% (752/922) of non-3NTP’s are expressed jointly. Indeed, there is a significant statistical association between joint expression and 3NTP’s (χ 2 = 17.75; d f = 1; p = 2.5 × 10-5).
Large and complex networks of distinct kinds of molecules orchestrate developmental mechanisms. High throughput tools, such as microarray experiments [40, 47] or large- scale yeast two-hybrid studies [15, 48, 49], are enabling the empirical identification of such networks. In order to study those systems, computational and mathematical tools are being used to organize and test the sufficiency of available experimental information to reproduce observed phenomena [50, 51].
Many networks share global features, such as a power-law connectivity distribution, or the property of having short paths between any two elements [52–55]. But several types of networks are different at an intermediate scale [27, 56]. For example, some networks contain certain subgraphs that are overrepresented with respect to randomly generated networks with the same number of nodes, edges and degree distribution [26–28]. Hence, a particular subgraph composition may characterize a certain kind of network . Such deviations from randomness may arise as a side-effect of a network’s generative rules . For example, Artzy-Randrup and collaborators considered a set of nodes, each fixed in a spatial location. Connecting such nodes in a random fashion, but with a bias towards nearer nodes, produces a statistically significant abundance of the same subgraphs that are overrepresented in the network of spatially-aggregated neurons of the nematod Caenorhabditis elegans. Another example concerns gene regulatory networks. Kuo et al. analyzed a gene regulatory network model in which interactions are defined by complementarity of binary ‘protein’ and ‘cis-regulatory’ sequences. They found that neutral processes, like duplication and mutation, suffice to produce an abundance of subgraphs similar to that of regulatory networks in yeast and Escherichia coli. However, in other cases, the abundance of some subgraphs may be attributed to adaptive optimization of the biological process controlled by the network, if the subgraph’s connection pattern confers an advantage to the organism [27, 60, 61]. This is very likely the case where alternative neutral scenarios, like duplication of entire subgraphs, have been discarded , or for subgraphs with dynamic properties that may increase the efficiency of tasks like information transfer [29, 30]. In sum, because neutral and adaptive explanations may be valid under different circumstances, the possible role of natural selection as an explanation for network structure should be tested for each case.
The formation of transcription factor tetramers is not exclusive of MADS domain proteins. Several transcriptional regulators do so, both in prokaryotes and eukaryotes. Pertinent examples include the lactose repressor LacI , Stat5 proteins , the IclR protein TtgV , bHLH proteins as in the Max-Myc heterotetramer , or proteins in the p53 family of transcription factors [66, 67]. In the present work we show that subgraphs compatible with MADS domain protein heterotetramers appear significantly more often than as expected by chance in the Arabidopsis network of physical interactions between these proteins. Formation of transcription factor oligomers may increase specificity of protein-DNA interaction  or allow synergistic binding to distinct DNA motifs . For these reasons, tetramer formation in MADS domain proteins could be beneficial . In fact, several of our observations suggest that the abundance of tetramer-like subgraphs is due to the advantages that the formation of tetramers of MADS domain proteins may confer to plants. First, the same kind of subgraphs is not significantly abundant in a network of physical interactions for a different family of plant transcription factors. These results are consistent with the hypothesis that the abundance of tetramer-like subgraphs can not be attributed to a common mechanism in the evolution of plant transcription factors. Moreover, we found that the mRNA of proteins in tetramer-like subgraphs is co-expressed more frequently than that of proteins in other subgraphs. Such an enhanced coordination of expression is expected if these subgraphs are associated to proteins that jointly exert their functions. One could think that sometimes gene duplication can increase the counts of what we call tetramer-like subgraphs, even without forming new tetramers. For instance, consider a protein X that can bind either protein Y or Z. If the gene coding for X is duplicated, then the new gene X’ would also have interactions with Y and Z. The resulting connection pattern would be a tetramer-like subgraph that does not correspond to a real functional tetramer. However, it must also be acknowledged that other equally probable duplications increase the frequencies of other subgraphs. Take as an example a set of four proteins in which A binds B, B binds C and C binds D. This interaction pattern is the linear subgraph A-B-C-D. Say D duplicates to produce D’ and D”. Then we would get twice the linear subgraph: A-B-C-D’ and A-B-C-D”. This argument is equally valid for duplication of the remaining three nodes. Although many subgraphs can increase their frequency because of duplication, not all of them are significantly abundant in the MADS network, when compared to what is expected by chance. Moreover, it is also noteworthy that joint-expression, as inferred from mRNA patterns, occurs significantly more often for tetramer-like subgraphs but not other subgraphs. In sum, our results suggest that formation of functional tetramers, as proposed by the ‘quartet model’ [17, 18], is a property of a large group of Arabidopsis MADS domain proteins, and is not restricted to the few experimentally validated cases [20–24].
In addition, we found that proteins in the monophyletic MIKC clade that take part in a tetramer-like subgraph will very likely be co-expressed in at least one Arabidopsis tissue. In contrast, we found no association between tetramer-like subgraphs and joint expression for sets of MADS domain proteins outside of the MIKC clade. These observations further support that functional tetramers comprise MIKC proteins exclusively. Our results, albeit forcefully based on a small number of non-MIKC proteins with mRNA expression patterns, are statistically sound. The prediction that MIKC proteins are the ones mainly involved in the formation of functional tetramers is consistent with the fact that all MADS domain tetramers that have been experimentally verified up to now comprise only MIKC proteins [20–24]. It well may be that the presence of the K domain, a specific trait of MIKC proteins that evolved along the plant lineage , facilitates tetramer formation. In fact, many of the interactions between proteins in tetramers occur through their K domains, as described in the Background section [14, 20, 21, 33, 34].
We also used the AtGenExpress data  in combination with the protein interaction information to extract protein sets in tetramer-like subgraphs for which the encoding genes are jointly expressed in at least one Arabidopsis tissue. In that manner we built a list of prospective tetramers (Additional file 2). Despite the overlap in expression, we cannot exclude that some of such protein sets do not represent functional tetramers occurring in planta. One likely cause would be that some interactions in a prospective tetramer may be mutually exclusive. This would occur if one protein binds two other proteins in the set through the same protein domain. Furthermore, the list of prospective tetramers is still rather long, though it is not unreasonable to think that a large number of functional tetramers occur in plants. However, it should be acknowledged that this is a major reduction of the search space, from 794 ∼ 39 × 106, considering only the 79 proteins in the MADS network, to just over a thousand. Moreover, analyses of this list of prospective tetramers allow the retrodiction of several experimental results, thus demonstrating its predictive value.
Most of the interactions with which we assembled the MADS network are pairwise, as they were mostly derived from yeast two-hybrid studies. One could argue that pairwise interactions do not provide information regarding higher-order protein complexes. Indeed, this would be an issue if each of the interactions between proteins in a tetramer only occurred in the presence of three or more of the tetramer’s constituent proteins. For the sake of simplicity, we assumed this was not the case. Notwithstanding, there is substantial evidence that backs up our assumption. This evidence concerns yeast two-hybrid data that supports the majority of the interactions that we derived from yeast three-hybrid assays. Moreover, we found the same significantly abundant subgraphs in our MADS network and in a network that excludes interactions supported only by yeast three-hybrid data (Tables S2 and S3 in Additional file 1). These observations substantiate our claim that pairwise interactions have a prevalent role in the formation of higher-order complexes of MADS domain proteins.
In sum, our work supports that the formation of functional tetramers is a widespread property of plant MADS domain transcription factors, specifically of MIKC proteins. Experimental research will be required to settle the issue, which is critical in the understanding of several mechanisms in plant development at the molecular level. However, the computational analysis that we put forward facilitates this endeavor by pointing those protein sets with higher chances of participating in functional tetrameric complexes.
Assemblage of the MADS network
In the MADS network, nodes represent the different MADS domain proteins in the plant Arabidopsis thaliana and links joining two nodes represent the potential of two proteins to interact physically. The main source of data to assemble the network is the study by de Folter et al. , that provides a semi-exhaustive map of protein-protein interactions between MADS domain proteins based on yeast two-hybrid assays. Another important source of information is a series of yeast three-hybrid assays for Arabidopsis MADS domain proteins . We took additional experimental evidence from the sources listed in Table S1 in Additional file 1. We did not consider those Arabidopsis MADS domain proteins for which we found no evidence of protein-protein interactions.
For every positive yeast two-hybrid result, an interaction between the two proteins involved was added. For every positive yeast three-hybrid experiment, two new interactions were added. For example, consider the yeast three-hybrid experiment involving the hypothetical proteins A, B, and C, where A was fused to a DNA-binding domain and C to a transactivation domain. If the result was positive and there was no yeast two-hybrid experiment indicating an interaction between A and C, then the interactions between A and B, and between B and C were included.
In order to obtain the MADS exp network, we eliminated from the MADS network all those proteins that lack expression data in the AtGenExpress data set. We also discarded those proteins that only interacted with proteins that lack expression data, as they would be devoid of interactions in the MADS exp network. We followed an analogous procedure to obtain the MADS At-TAX network.
We built a catalogue with all possible manners of connecting up to four nodes in a single component (Figure 2). In other words, the connection patterns that we follow are not disjoint. We accomplished this by exhaustive enumeration and dismissing those connection patterns that do not form a single component and those that were isomorphic to connection patterns already in the catalogue. Two subgraphs are isomorphic if there is a one-to-one correspondence between their node sets which preserves adjacency .
In order to count the number of n-node connection patterns in a given network, we determined for each possible n-node set in the network whether it formed a single component. If that was the case, we searched for its isomorphic subgraph in our catalogue to increase its count.
In order to obtain an ensemble of randomized networks comparable to the MADS network we performed the ‘switching’ algorithm described by Milo et al. . For each randomized network in the ensemble we started with the original MADS network as seed, and iteratively we randomly picked two connections and exchanged two of their ends. For example, assume we pick connections joining proteins A and B, on one hand, and proteins C and D, on the other. We then replaced those connections for one joining A and D, and another joining C and B. In this manner, each node preserves its number of connections, but now they are linked to different nodes. We switched pairs of connections s times, with s taken with an uniform probability distribution from the interval [100E, 200E]. Here, E stands for the number of connections in the MADS network. The large value of s guarantees that all traces of the original structure of the MADS network is lost in its randomized counterparts. At the same time, this procedure preserves the number of connections that each node has, and thus the network’s connectivity distribution.
For building randomized networks that, in addition to maintaining the connectivity distribution, preserve the number of self-interactions we modified slightly the procedure described above. The modification merely consisted on prohibiting both the formation and loss of self-interactions.
We obtained expression data for Arabidopsis MADS-box genes from the developmental data set of the AtGenExpress expression atlas (http://weigelworld.org/resources/microarray/AtGenExpress) . We excluded from our analyses those samples that came from mutant backgrounds. We used gcRMA to obtain expression data. We considered a certain gene in a given sample as expressed only if the log2 of its average expression was greater than a pre-specified threshold. We followed Immink et al.  and set this threshold to four for all of our analyses.
We thank three anonymous reviewers, who provided useful comments to improve our manuscript. CE-S acknowledges financial support from a postdoctoral fellowship from Mexican National Council of Science and Technology (CONACyT), from UASLP’s FAI (Fondo de Apoyo para la Investigación) [C13-FAI-03-11.11], and from Programa para el Mejoramiento del Profesorado (PROMEP/103.5/13/6575). SdF acknowledges The Company of Biologists for a Development Traveling Fellowship. The work in the SdF lab was financed by the CONACyT grant 82826 and 177739. Research in ERA-B laboratory was supported by: CONACyT (81433; 180098; 180380; 167705; 152649; 105678; Red Temática de Investigación: “Complejidad, Ciencia y Sociedad”, 124909, 81542), PAPIIT, UNAM (IN204011-3; IN226510-3; IN203113-3; IB201212-2, IN229009-3) and UC-MEXUS ECO-IE415 grants. ERA-B was supported by the Miller Institute for Basic Research in Science, University of California, Berkeley, California, United States, while spending a sabbatical leave in the laboratory of Chelsea Spetch at UC-B. RGHI and GCA are supported by a grant from the Netherlands Consortium for Systems Biology (NCSB).
- Gramzow L, Ritz MS, Theißen G: On the origin of MADS-domain transcription factors. Trends Genet. 2010, 26 (4): 149-153. 10.1016/j.tig.2010.01.004. doi:10.1016/j.tig.2010.01.004View ArticlePubMedGoogle Scholar
- Gramzow L, Theißen G: A hitchhiker’s guide to the MADS world of plants. Genome Biol. 2010, 11 (6): 214-10.1186/gb-2010-11-6-214.PubMed CentralView ArticlePubMedGoogle Scholar
- Lawton-Rauh AL, Alvarez-Buylla ER, Purugganan MD: Molecular evolution of flower development. Trends Ecol Evol. 2000, 15 (4): 144-149. 10.1016/S0169-5347(99)01816-9.View ArticlePubMedGoogle Scholar
- Parenicová L, de Folter S, Kieffer M, Horner DS, Favalli C, Busscher J, Cook HE, Ingram RM, Kater MM, Davies B, Angenent GC, Colombo L: Molecular and phylogenetic analyses of the complete MADS-box transcription factor family in Arabidopsis: New openings to the MADS world. Plant Cell. 2003, 15 (7): 1538-1551. 10.1105/tpc.011544. doi:10.1105/tpc.011544.thesePubMed CentralView ArticlePubMedGoogle Scholar
- Martínez-Castilla LP, Alvarez-Buylla ER: Adaptive evolution in the Arabidopsis MADS-box gene family inferred from its complete resolved phylogeny. Proc Natl Acad Sci USA. 2003, 100 (23): 13407-13412. 10.1073/pnas.1835864100. doi:10.1073/pnas.1835864100PubMed CentralView ArticlePubMedGoogle Scholar
- Theißen G, Becker A, Di Rosa A, Kanno A, Kim JT, Münster T, Winter K-U, Saedler H: A short history of MADS-box genes in plants. Plant Mol Biol. 2000, 42 (1): 115-149. 10.1023/A:1006332105728.View ArticlePubMedGoogle Scholar
- Masiero S, Colombo L, Grini PE, Schnittger A, Kater MM: The emerging importance of type I MADS box transcription factors for plant reproduction. Plant Cell. 2011, 23 (3): 865-872. 10.1105/tpc.110.081737. doi:10.1105/tpc.110.081737PubMed CentralView ArticlePubMedGoogle Scholar
- Airoldi CA, Davies B: Gene duplication and the evolution of plant MADS-box transcription factors. J Genet Genomics. 2012, 39 (4): 157-165. 10.1016/j.jgg.2012.02.008. doi:10.1016/j.jgg.2012.02.008View ArticlePubMedGoogle Scholar
- Smaczniak C, Immink RGH, Angenent GC, Kaufmann K: Developmental and evolutionary diversity of plant MADS-domain factors: insights from recent studies. Development. 2012, 139 (17): 3081-3098. 10.1242/dev.074674. doi:10.1242/dev.074674View ArticlePubMedGoogle Scholar
- Herrera-Ubaldo H, Zanchetti E, Colombo L, de Folter S: Protein interactions guiding carpel and fruit development in Arabidopsis. Plant Biosyst. in press. doi:10.1080/11263504.2013.870255
- Hernández-Hernández T, Martínez-Castilla LP, Alvarez-Buylla ER: Functional diversification of B MADS-box homeotic regulators of flower development: Adaptive evolution in protein-protein interaction domains after major gene duplication events. Mol Biol Evol. 2007, 24 (2): 465-481. doi:10.1093/molbev/msl182View ArticlePubMedGoogle Scholar
- Alvarez-Buylla ER, Pelaz S, Liljegren SJ, Gold SE, Burgeff C, Ditta GS, Ribas de Pouplana L, Martínez-Castilla LP, Yanofsky MF: An ancestral MADS-box gene duplication occurred before the divergence of plants and animals. Proc Natl Acad Sci USA. 2000, 97 (10): 5328-5333. 10.1073/pnas.97.10.5328.PubMed CentralView ArticlePubMedGoogle Scholar
- Riechmann JL, Krizek BA, Meyerowitz EM: Dimerization specificity of Arabidopsis MADS domain homeotic proteins APETALA1, APETALA3, PISTILLATA, and AGAMOUS. Proc Natl Acad Sci USA. 1996, 93 (10): 4793-4798. 10.1073/pnas.93.10.4793.PubMed CentralView ArticlePubMedGoogle Scholar
- Egea-Cortines M, Saedler H, Sommer H: Ternary complex formation between the MADS-box proteins SQUAMOSA, DEFICIENS and GLOBOSA is involved in the control of floral architecture in Antirrhinum majus. EMBO J. 1999, 18 (19): 5370-5379. 10.1093/emboj/18.19.5370. doi:10.1093/emboj/18.19.5370PubMed CentralView ArticlePubMedGoogle Scholar
- de Folter S, Immink RGH, Kieffer M, Parenicová L, Henz SR, Weigel D, Busscher M, Kooiker M, Colombo L, Kater MM, Davies B, Angenent GC: Comprehensive interaction map of the Arabidopsis MADS box transcription factors. Plant Cell. 2005, 17 (5): 1424-1433. 10.1105/tpc.105.031831. doi:10.1105/tpc.105.031831.1PubMed CentralView ArticlePubMedGoogle Scholar
- Immink RGH, Kaufmann K, Angenent GC: The ‘ABC’ of MADS domain protein behaviour and interactions. Semin Cell Dev Biol. 2010, 21 (1): 87-93. 10.1016/j.semcdb.2009.10.004. doi:10.1016/j.semcdb.2009.10.004View ArticlePubMedGoogle Scholar
- Theißen G: Development of floral organ identity: stories from the MADS house. Curr Opin Plant Biol. 2001, 4 (1): 75-85. 10.1016/S1369-5266(00)00139-4.View ArticlePubMedGoogle Scholar
- Theißen G, Saedler H: Floral quartets. Nature. 2001, 409: 469-471. 10.1038/35054172.View ArticlePubMedGoogle Scholar
- de Folter S, Angenent GC: trans meets cis in MADS science. Trends Plant Sci. 2006, 11 (5): 224-231. 10.1016/j.tplants.2006.03.008. doi:10.1016/j.tplants.2006.03.008View ArticlePubMedGoogle Scholar
- Melzer R, Theißen G: Reconstitution of ‘floral quartets’ in vitro involving class B and class E floral homeotic proteins. Nucleic Acids Res. 2009, 37 (8): 2723-2736. 10.1093/nar/gkp129. doi:10.1093/nar/gkp129PubMed CentralView ArticlePubMedGoogle Scholar
- Melzer R, Verelst W, Theißen G: The class E floral homeotic protein SEPALLATA3 is sufficient to loop DNA in ‘floral quartet’-like complexes in vitro. Nucleic Acids Res. 2009, 37 (1): 144-157. 10.1093/nar/gkn900. doi:10.1093/nar/gkn900PubMed CentralView ArticlePubMedGoogle Scholar
- Honma T, Goto K: Complexes of MADS-box proteins are sufficient to convert leaves into floral organs. Nature. 2001, 409: 525-529. 10.1038/35054083.View ArticlePubMedGoogle Scholar
- Pelaz S, Tapia-López R, Alvarez-Buylla ER, Yanofsky MF: Conversion of leaves into petals in Arabidopsis. Curr Biol. 2001, 11 (3): 182-184. 10.1016/S0960-9822(01)00024-0.View ArticlePubMedGoogle Scholar
- Smaczniak C, Immink RGH, Muiño JM, Blanvillain R, Busscher M, Busscher-Lange J, Dinh QP, Liu S, Westphal AH, Boeren S, Parcy F, Xu L, Carles CC, Angenent GC, Kaufmann K: Characterization of MADS-domain transcription factor complexes in Arabidopsis flower development. Proc Natl Acad Sci USA. 2012, 109 (5): 1560-1565. 10.1073/pnas.1112871109. doi:10.1073/pnas.1112871109/-/DCSupplemental.http://www.pnas.org/cgi/doi/10.1073/pnas.1112871109 PubMed CentralView ArticlePubMedGoogle Scholar
- Immink RGH, Tonaco IAN, de Folter S, Shchennikova A, van Dijk ADJ, Busscher-Lange J, Borst JW, Angenent GC: SEPALLATA3: the ‘glue’ for MADS box transcription factor complex formation. Genome Biol. 2009, 10 (2): 24-10.1186/gb-2009-10-2-r24. doi:10.1186/gb-2009-10-2-r24View ArticleGoogle Scholar
- Shen-Orr SS, Milo R, Mangan S, Alon U: Network motifs in the transcriptional regulation network of Escherichia coli. Nat Genet. 2002, 31 (1): 64-68. 10.1038/ng881. doi:10.1038/ng881View ArticlePubMedGoogle Scholar
- Milo R, Shen-Orr SS, Itzkovitz S, Kashtan N, Chklovskii D, Alon U: Network motifs: simple building blocks of complex networks. Science. 2002, 298 (5594): 824-827. 10.1126/science.298.5594.824. doi:10.1126/science.298.5594.824View ArticlePubMedGoogle Scholar
- Milo R, Itzkovitz S, Kashtan N, Levitt R: Superfamilies of evolved and designed networks. Science. 2004, 303 (March): 8-10.Google Scholar
- Mangan S, Alon U: Structure and function of the feed-forward loop network motif. Proc Natl Acad Sci USA. 2003, 100 (21): 11980-11985. 10.1073/pnas.2133841100. doi:10.1073/pnas.2133841100PubMed CentralView ArticlePubMedGoogle Scholar
- Alon U: Network motifs: theory and experimental approaches. Nat Rev Genet. 2007, 8 (6): 450-461. 10.1038/nrg2102. doi:10.1038/nrg2102View ArticlePubMedGoogle Scholar
- Mizukami Y, Huang H, Tudor M, Hu Y, Ma H: Functional domains of the floral regulator AGAMOUS: characterization of the DNA binding domain and analysis of dominant negative mutations. Plant Cell. 1996, 8 (5): 831-845. doi:10.1105/tpc.8.5.831PubMed CentralView ArticlePubMedGoogle Scholar
- Huang H, Tudor M, Su T, Zhang Y, Hu Y, Ma H: DNA binding properties of two Arabidopsis MADS domain proteins binding consensus and dimer formation. Plant Cell. 1996, 8 (1): 81-94. doi:10.1105/tpc.8.1.81PubMed CentralView ArticlePubMedGoogle Scholar
- Fan H-Y, Hu Y, Tudor M, Ma H: Specific interactions between the K domains of AG and AGLs, members of the MADS domain family of DNA binding proteins. Plant J. 1997, 12 (5): 999-1010. 10.1046/j.1365-313X.1997.12050999.x.View ArticlePubMedGoogle Scholar
- Yang Y, Jack T: Defining subdomains of the K domain important for protein-protein interactions of plant MADS proteins. Plant Mol Biol. 2004, 55 (1): 45-59. 10.1007/s11103-004-0416-7. doi:10.1007/s11103-004-0416-7View ArticlePubMedGoogle Scholar
- Martín-Trillo M, Cubas P: TCP genes: a family snapshot ten years later. Trends Plant Sci. 2010, 15 (1): 31-39. 10.1016/j.tplants.2009.11.003. doi:10.1016/j.tplants.2009.11.003View ArticlePubMedGoogle Scholar
- Danisman S, van derWal F, Dhondt S, Waites R, de Folter S, Bimbo A, van Dijk ADJ, Muino JM, Cutri L, Dornelas MC, Angenent GC, Immink RGH: Arabidopsis class I and class II TCP transcription factors regulate jasmonic acid metabolism and leaf development antagonistically. Plant Physiol. 2012, 159 (4): 1511-1523. 10.1104/pp.112.200303. doi:10.1104/pp.112.200303PubMed CentralView ArticlePubMedGoogle Scholar
- Favaro R, Pinyopich A, Battaglia R, Kooiker M, Borghi L, Ditta G, Yanofsky MF, Kater MM, Colombo L: MADS-box protein complexes control carpel and ovule development in Arabidopsis. Plant Cell. 2003, 15 (11): 2603-2611. 10.1105/tpc.015123. doi:10.1105/tpc.015123.2PubMed CentralView ArticlePubMedGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc B. 1995, 57 (1): 289-300.Google Scholar
- Danisman S, van Dijk ADJ, Bimbo A, van der Wal F, Hennig L, de Folter S, Angenent GC, Immink RGH: Analysis of functional redundancies within the Arabidopsis TCP transcription factor family. J Exp Bot. 2013, 64 (18): 5673-5685. 10.1093/jxb/ert337. doi:10.1093/jxb/ert337PubMed CentralView ArticlePubMedGoogle Scholar
- Schmid M, Davison TS, Henz SR, Pape UJ, Demar M, Vingron M, Schölkopf B, Weigel D, Lohmann JU: A gene expression map of Arabidopsis thaliana development. Nat Genet. 2005, 37 (5): 501-506. 10.1038/ng1543. doi:10.1038/ng1543View ArticlePubMedGoogle Scholar
- Laubinger S, Zeller G, Henz SR, Sachsenberg T, Widmer CK, Naouar N, Vuylsteke M, Schölkopf B, Rätsch G, Weigel D: At-TAX: a whole genome tiling array resource for developmental expression analysis and transcript identification in Arabidopsis thaliana. Genome Biol. 2008, 9 (7): 112-10.1186/gb-2008-9-7-r112. doi:10.1186/gb-2008-9-7-r112View ArticleGoogle Scholar
- Pelaz S, Ditta GS, Baumann E, Wisman E, Yanofsky MF: B and C floral organ identity functions require SEPALLATA MADS-box genes. Nature. 2000, 405: 200-203. 10.1038/35012103.View ArticlePubMedGoogle Scholar
- Davies B, Egea-Cortines M, de Andrade Silva E, Saedler H, Sommer H: Multiple interactions amongst floral homeotic MADS box proteins. EMBO J. 1996, 15 (16): 4330-4343.PubMed CentralPubMedGoogle Scholar
- Gregis V, Sessa A, Colombo L, Kater MM: AGAMOUS-LIKE24 and SHORT VEGETATIVE PHASE determine floral meristem identity in Arabidopsis. Plant J. 2008, 56 (6): 891-902. 10.1111/j.1365-313X.2008.03648.x. doi:10.1111/j.1365-313X.2008.03648.xView ArticlePubMedGoogle Scholar
- Kempin S, Savidge B, Yanofsky MF: Molecular basis of the cauliflower phenotype in Arabidopsis. Science. 1995, 267 (5197): 522-525. 10.1126/science.7824951. doi:10.1126/science.7824951View ArticlePubMedGoogle Scholar
- Liljegren SJ, Ditta GS, Eshed Y, Savidge B, Bowman JL, Yanofsky MF: SHATTERPROOF MADS-box genes control seed dispersal in Arabidopsis. Nature. 2000, 404 (6779): 766-770. 10.1038/35008089. doi:10.1038/35008089View ArticlePubMedGoogle Scholar
- Gómez-Mena C, de Folter S, Costa MMR, Angenent GC, Sablowski R: Transcriptional program controlled by the floral homeotic gene AGAMOUS during early organogenesis. Development. 2005, 132 (3): 429-438. 10.1242/dev.01600. doi:10.1242/dev.01600View ArticlePubMedGoogle Scholar
- Ito T, Chiba T, Ozawa R, Yoshida M, Hattori M, Sakaki Y: A comprehensive two-hybrid analysis to explore the yeast protein interactome. Proc Natl Acad Sci USA. 2001, 98 (8): 4569-4574. 10.1073/pnas.061034498. doi:10.1073/pnas.061034498PubMed CentralView ArticlePubMedGoogle Scholar
- Uetz P, Giot L, Cagney G, Mansfield TA, Judson RS, Knight JR, Lockshon D, Narayan V, Srinivasan M, Pochart P, Qureshi-Emili A, Li Y, Godwin B, Conover D, Kalbfleisch T, Vijayadamodar G, Yang M, Johnston M, Fields S, Rothberg JM: A comprehensive analysis of protein-protein interactions in Saccharomyces cerevisiae. Nature. 2000, 403 (6770): 623-627. 10.1038/35001009. doi:10.1038/35001009View ArticlePubMedGoogle Scholar
- Alvarez-Buylla ER, Benítez M, Dávila EB, Chaos A, Espinosa-Soto C, Padilla-Longoria P: Gene regulatory network models for plant development. Curr Opin Plant Biol. 2007, 10 (1): 83-91. 10.1016/j.pbi.2006.11.008. doi:10.1016/j.pbi.2006.11.008View ArticlePubMedGoogle Scholar
- van Dijk ADJ, Morabito G, Fiers M, van Ham RCHJ, Angenent GC, Immink RGH: Sequence motifs in MADS transcription factors responsible for specificity and diversification of protein-protein interaction. PLoS Comput Biol. 2010, 6 (11): 1001017-10.1371/journal.pcbi.1001017. doi:10.1371/journal.pcbi.1001017View ArticleGoogle Scholar
- Albert R: Scale-free networks in cell biology. J Cell Sci. 2005, 118 (Pt 21): 4947-4957. doi:10.1242/jcs.02714View ArticlePubMedGoogle Scholar
- Albert R: Network inference, analysis, and modeling in systems biology. Plant Cell. 2007, 19 (11): 3327-3338. 10.1105/tpc.107.054700. doi:10.1105/tpc.107.054700PubMed CentralView ArticlePubMedGoogle Scholar
- Barabási A-L: Linked. The New Science of Networks. 2002: 280, Cambridge: Perseus PublishingGoogle Scholar
- Barabási A-L, Oltvai ZN: Network biology: understanding the cell’s functional organization. Nat Rev Genet. 2004, 5 (2): 101-113. 10.1038/nrg1272. doi:10.1038/nrg1272View ArticlePubMedGoogle Scholar
- Maslov S, Sneppen K: Specificity and stability in topology of protein networks. Science. 2002, 296 (5569): 910-913. 10.1126/science.1065103. doi:10.1126/science.1065103View ArticlePubMedGoogle Scholar
- Solé RV, Valverde S: Are network motifs the spandrels of cellular complexity?. Trends Ecol Evol. 2006, 21 (8): 419-422. 10.1016/j.tree.2006.05.013. doi:10.1016/j.tree.2006.06.001View ArticlePubMedGoogle Scholar
- Artzy-Randrup Y, Fleishman SJ, Ben-Tal N, Stone L: Comment on “Network motifs: simple building blocks of complex networks” and “Superfamilies of evolved and designed networks”. Science. 2004, 305: 1107-View ArticlePubMedGoogle Scholar
- Kuo PD, Banzhaf W, Leier A: Network topology and the evolution of dynamics in an artificial genetic regulatory network model created by whole genome duplication and divergence. Biosystems. 2006, 85 (3): 177-200. 10.1016/j.biosystems.2006.01.004. doi:10.1016/j.biosystems.2006.01.004View ArticleGoogle Scholar
- Alon U: Biological networks: the tinkerer as an engineer. Science. 2003, 301 (5641): 1866-1867. 10.1126/science.1089072. doi:10.1126/science.1089072View ArticlePubMedGoogle Scholar
- Conant GC, Wagner A: Convergent evolution of gene circuits. Nat Genet. 2003, 34 (3): 264-266. 10.1038/ng1181. doi:10.1038/ng1181View ArticlePubMedGoogle Scholar
- Ramot R, Kishore Inampudi K, Wilson CJ: Lactose repressor experimental folding landscape: fundamental functional unit and tetramer folding mechanisms. Biochemistry. 2012, 51 (38): 7569-7579. 10.1021/bi300545f. doi:10.1021/bi300545fView ArticlePubMedGoogle Scholar
- John S, Vinkemeier U, Soldaini E, Darnell JE, Leonard WJ: The significance of tetramerization in promoter recruitment by Stat5. Mol Cell Biol. 1999, 19 (3): 1910-1918.PubMed CentralView ArticlePubMedGoogle Scholar
- Fillet S, Krell T, Morel B, Lu D, Zhang X, Ramos JL: Intramolecular signal transmission in a tetrameric repressor of the IclR family. Proc Natl Acad Sci USA. 2011, 108 (37): 15372-15377. 10.1073/pnas.1018894108. doi:10.1073/pnas.1018894108/-/ DCSupplemental.http://www.pnas.org/cgi/doi/10.1073/pnas.1018894108 PubMed CentralView ArticlePubMedGoogle Scholar
- Nair SK, Burley SK: X-ray structures of Myc-Max and Mad-Max recognizing DNA. Molecular bases of regulation by proto-oncogenic transcription factors. Cell. 2003, 112 (2): 193-205. 10.1016/S0092-8674(02)01284-9.View ArticlePubMedGoogle Scholar
- Joerger AC, Rajagopalan S, Natan E, Veprintsev DB, Robinson CV, Fersht AR: Structural evolution of p53, p63, and p73: implication for heterotetramer formation. Proc Natl Acad Sci USA. 2009, 106 (42): 17705-17710. 10.1073/pnas.0905867106. doi:10.1073/pnas.0905867106PubMed CentralView ArticlePubMedGoogle Scholar
- Natan E, Joerger AC: Structure and kinetic stability of the p63 tetramerization domain. J Mol Biol. 2012, 415 (3): 503-513. 10.1016/j.jmb.2011.11.007. doi:10.1016/j.jmb.2011.11.007PubMed CentralView ArticlePubMedGoogle Scholar
- Harary F: Graph Theory. 1969: 280, Addison-Wesley: ReadingGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.