Evolution of strigolactone receptors by gradual neo-functionalization of KAI2 paralogues
© Bennett et al. 2017
Received: 22 February 2017
Accepted: 13 June 2017
Published: 29 June 2017
Strigolactones (SLs) are a class of plant hormones that control many aspects of plant growth. The SL signalling mechanism is homologous to that of karrikins (KARs), smoke-derived compounds that stimulate seed germination. In angiosperms, the SL receptor is an α/β-hydrolase known as DWARF14 (D14); its close homologue, KARRIKIN INSENSITIVE2 (KAI2), functions as a KAR receptor and likely recognizes an uncharacterized, endogenous signal (‘KL’). Previous phylogenetic analyses have suggested that the KAI2 lineage is ancestral in land plants, and that canonical D14-type SL receptors only arose in seed plants; this is paradoxical, however, as non-vascular plants synthesize and respond to SLs.
We have used a combination of phylogenetic and structural approaches to re-assess the evolution of the D14/KAI2 family in land plants. We analysed 339 members of the D14/KAI2 family from land plants and charophyte algae. Our phylogenetic analyses show that the divergence between the eu-KAI2 lineage and the DDK (D14/DLK2/KAI2) lineage that includes D14 occurred very early in land plant evolution. We show that eu-KAI2 proteins are highly conserved, and have unique features not found in DDK proteins. Conversely, we show that DDK proteins show considerable sequence and structural variation to each other, and lack clearly definable characteristics. We use homology modelling to show that the earliest members of the DDK lineage structurally resemble KAI2 and that SL receptors in non-seed plants likely do not have D14-like structure. We also show that certain groups of DDK proteins lack the otherwise conserved MORE AXILLARY GROWTH2 (MAX2) interface, and may thus function independently of MAX2, which we show is highly conserved throughout land plant evolution.
Our results suggest that D14-like structure is not required for SL perception, and that SL perception has relatively relaxed structural requirements compared to KAI2-mediated signalling. We suggest that SL perception gradually evolved by neo-functionalization within the DDK lineage, and that the transition from KAI2-like to D14-like protein may have been driven by interactions with protein partners, rather than being required for SL perception per se.
KeywordsStrigolactone signalling Strigolactone evolution Phylogenetics Neo-functionalization
Plant hormones are a key link between environmental stimuli and development, allowing local information to be used systemically across the plant body. Strigolactones (SLs) are a recently identified class of terpenoid lactone hormones that neatly epitomize this concept. SLs are primarily synthesized by a core pathway involving a carotene isomerase (DWARF27), two carotenoid cleavage dioxygenases (CCD7 and CCD8)  and a cytochrome P450 enzyme (MAX1). SL synthesis is strongly upregulated by phosphate deficiency in the rhizosphere , increasing the pool of SL molecules in the root. In many flowering plants (angiosperms), SLs are exuded into the soil through the action of specific SL transporters and serve to attract mycorrhizal fungi ; the resulting symbioses provide the plants with phosphate in exchange for reduced carbon. SLs also act locally to regulate root system architecture; the precise effects seem to vary from species to species, but increased SL levels may promote increased nutrient foraging . Finally, a significant proportion of the SL pool produced in the root is transported into the shoot system via the xylem , where it has a well-defined set of effects on shoot growth and development [6, 7]. SL has an inhibitory effect on shoot branching, thereby coupling shoot growth to nutrient availability . SL responses thus form an integrated stimulus-response system acting over long distances both within the plant body and in its immediate environment.
Like several other plant hormonal signalling pathways, canonical SL signalling is mediated through ubiquitin-mediated degradation of target proteins (reviewed in ). The SL receptors for this signalling pathway are members of the DWARF14 (D14) class of α/β-hydrolase proteins, which are an unusual combination of enzyme and receptor [8, 9]. D14 proteins bind and then cleave SL molecules, producing a covalently linked intermediate molecule (CLIM) that is covalently bound to the receptor [8, 9]. SL signalling is mediated through the interaction of D14 with the MORE AXILLARY GROWTH2 (MAX2) class of F-box proteins, which forms part of an SCF (SKP1-CULLIN-F-BOX) E3 ubiquitin ligase [10–13]. Together, the covalent binding of CLIM and the interaction with SCFMAX2 allow D14 to undergo a stable conformational change that drives onward signalling [8, 9]. Although other targets have been proposed [14, 15], it is now clear that the principal proteolytic targets of SL signalling are proteins of the SMAX1-LIKE7/DWARF53 (SMXL7/D53) class [16–21]. The exact sequence of events is unclear, but it is probably after conformational change that D14 stably recruits SMXL7 to the complex; certainly, the D14-SMXL7 interaction is enhanced by SL [16, 17, 19, 20]. Events downstream of SMXL7 degradation are currently poorly defined; SMXL7 has been proposed to act both transcriptionally and non-transcriptionally [7, 22]. It may be that SMXL7 is a multi-functional protein that can regulate multiple cellular processes .
Intriguingly, a second pathway in angiosperms signals through SCFMAX2, forming a biochemical and evolutionary parallel to SL signalling. This pathway is defined by the KARRIKIN INSENSITIVE2 (KAI2) α/β-hydrolase protein, a close relative of D14. kai2 mutants have a range of developmental phenotypes [18, 21, 23] and are insensitive to the germination-promoting effects of smoke-derived ‘karrikins’ (KARs) . It has been hypothesized that karrikins promote germination by mimicking an as-yet-unidentified endogenous KAI2 ligand (‘KL’) [24, 25]. The KAI2 orthologue in rice (D14-LIKE) is also required for the establishment of mycorrhizal associations in the root system . It is currently unclear whether D14-LIKE perceives a fungal signal or endogenous KL in this context. As with D14, MAX2 (and orthologues) is required for both responses to karrikins and for other aspects of KAI2-dependent signalling [18, 21, 27]. Furthermore, the presumptive proteolytic targets of KAI2-SCFMAX2 signalling are close homologues of SMXL7; in Arabidopsis, these are SMXL2 and SMAX1 (SUPPRESSOR OF MAX2 1). Mutation of SMAX1 and SMXL2 suppresses the kai2-related phenotypes present in the max2 mutant, producing phenotypes that mimic constitutive karrikin responses [18, 28, 29]. In the Arabidopsis genome, there are further homologues of D14 and SMAX1, namely DWARF14-LIKE2 (DLK2) and SMXL3, SMXL4 and SMXL5, but the function of these proteins and their relationship to SL/KL signalling is currently unclear [23, 28].
The evolutionary history of SLs represents an intriguing and unresolved problem. SLs have been identified in most land plant groups, and in some related groups of charophyte algae . However, unambiguous CCD8 orthologues have not been identified in charophytes or liverworts  (a possible sister group to other land plants ). Moreover, ccd8 mutants in the moss Physcomitrella patens still produce some SLs , which suggests that there may be alternative pathways for SL synthesis [7, 33]. Even more uncertainty surrounds the origin of the canonical SL signalling pathway. Unambiguous D14 orthologues have only been identified in seed plants (gymnosperms and angiosperms), and they seem to be absent from mosses and liverworts [30, 34]. Conversely, it has been suggested that unambiguous KAI2 orthologues are present in charophytes, liverworts and mosses . This has led to the suggestion that KAI2 proteins could function as receptors for SLs in non-vascular plants, or that SL signalling occurs by non-canonical mechanisms in these lineages [7, 22]. Supporting the plausibility of the former hypothesis, it was recently shown that SL receptors evolved from KAI2 paralogues in parasitic plants within the Orobanchaceae [35–37]. In addition, MAX2 orthologues have so far only been identified in land plants , and while MAX2 is present in P. patens, Ppmax2 mutants do not resemble Ppccd8 mutants, suggesting that MAX2 may not be involved in SL signalling in mosses [22, 39]. Thus, even if KAI2 proteins can act as SL receptors in mosses, they may not act through SCFMAX2-mediated protein degradation. SMXL proteins are present in P. patens, but their function has not been investigated. Thus, while there is clear evidence for SL sensitivity in mosses, it is possible that this occurs through separate mechanisms to those in angiosperms. This would contrast strongly with the auxin signalling pathway for instance, which is completely conserved throughout land plants [40–42]. To resolve the evolutionary history of SL signalling, we have undertaken a major phylogenetic re-assessment of the D14/KAI2 family.
Preliminary analysis of the D14/KAI2 family
Major clades in the D14/KAI2 family
Multiple analyses support an early origin for the DDK super-clade
Irrespective of the underlying alignment and methodology, all analyses agreed on a basic topology for the family, with a deep duplication near the base of the land plants creating two super-clades. The first lineage contains KAI2 sequences from angiosperms and closely related sequences from gymnosperms, monilophytes, lycophytes, mosses and liverworts; we therefore named this clade eu-KAI2 (Table 1). The second super-clade contains sequences from mosses that have previously been described as KAI2-like [23, 34, 49], sequences from lycophytes and monilophytes that do not resemble known proteins, the previously characterized D14 and DLK2 genes from angiosperms and homologous genes from gymnosperms (Table 1). To reflect the mixed composition of this clade, we named it ‘DDK’ (for D14/DLK2/KAI2); we also used this name for the monilophyte and lycophyte sequences in the clade. The lycophyte DDK group contains the Selaginella moellendorffii gene previously described as ‘SmKAI2b’ , but we believe DDK designation better reflects the evolutionary context of these proteins. We observed some variation in the composition of the eu-KAI2 clade, partly as a result of the erratic behavior of the lycophyte KAI2 sequences. However, the moss KAI2E/F, lycophyte DDK, monilophyte DDK, gymnosperm D14, DLK4 and DLK23 and angiosperm D14 and DLK23 clades were associated into a single large clade in every analysis we performed, although the internal branching order did vary somewhat between analyses. This basic topology was evident even in very early analyses (Additional file 4).
Only two clades were inconsistently placed. The hornwort KAI2 clade is the most problematic in our analyses, mirroring the uncertainty about the position of the hornworts themselves in organismal phylogeny . In some analyses the hornwort KAI2 clade is placed in the eu-KAI2 lineage, between mosses and vascular plants (Additional file 3). Alternatively, it is also placed at the base of the eu-KAI2 lineage (Fig. 1) or as a sister clade to all other land plant D14/KAI2 sequences (Fig. 2, Additional file 2). None of these positions alter the interpretation of a deep duplication in the family, but they do affect its inferred timing. The liverwort KAI2B clade occurs either at the base of the DDK or eu-KAI2 lineages in different trees. In analyses performed without charophyte and lycophyte KAI2 sequences, it is always associated with the DDK lineage (Fig. 2, Additional file 3). This is also the case in some analyses including charophyte sequences (Additional file 2). The position at the base of the eu-KAI2 clade in some trees is likely to be erroneous and is probably caused by the slight misplacement of charophyte sequences. For instance, the liverwort-hornwort-liverwort branching order at the base of the eu-KAI2 clade in Fig. 1 is highly improbable. Rooting this tree with the hornwort KAI2 clade (to match Fig. 2) produces balanced eu-KAI2 and DDK clades, with realistic branching order, except for the inclusion of the charophyte sequences as an in-group within the DDK clade (Additional file 5). We believe the most parsimonious scenario is that KAI2B is part of the DDK clade.
Collectively, our phylogenetic analyses push the origin of the D14 lineage back much earlier than proposed in previous phylogenies that suggested an origin in the vascular plants  or within the seed plants . They resolve the enigmatic placement of SmKAI2b and divergent KAI2 sequences from P. patens in previous phylogenies [34, 49]. They also provide a convincing explanation for the presence of two distinct D14/KAI2 clades in most major plant groups. Key to this reconstruction topology is the placement of liverwort and moss clades with apparently KAI2-like primary protein structure (KAI2B and KAI2E/F respectively) in the DDK lineage. We wanted to test the robustness of this somewhat unexpected conclusion, and used a variety of methods to do so.
Non-parametric bootstrap analyses performed in GARLI did not provide very high levels of support for most of the nodes along the backbone of the tree (Fig. 1). However, bootstrap values were higher in reconstructions that excluded charophyte and lycophyte KAI2 sequences (Fig. 2). We next tested whether the recovered topology was stable to perturbations in the dataset. We re-ran our analysis multiple times, removing each DDK clade in turn (see Methods). Our analysis suggests that the placement of KAI2B is sensitive to the dataset used, but that the rest of the DDK clade is very stably associated (Additional file 6). Finally, we assessed whether our general topology is congruent with previous analyses. We observed that in , the Marchantia polymorpha KAI2A and KAI2B sequences do not group together, and neither do the P. patens KAI2C/D and KAI2E/F. This is consistent with our analyses. We repeated our analysis using a set of sequences pruned to match  and found essentially the same tree as in that study (Additional file 6). Furthermore, if we rooted the tree with a eu-KAI2 sequence, we observed essentially the same topology as in our study (Additional file 7). This shows that the difference in final topology between our study and  does not result from any particular methodological differences, but from our more densely populated sequence set.
Diverse evolutionary histories in the D14/KAI2 family
Conversely, the evolutionary history of the DDK clade is one of divergence and diversification. The liverwort and moss clades (KAI2B, KAI2E/F) are on relatively short branches (Fig. 3) and have been categorized previously as encoding KAI2-like proteins. The lycophyte and monilophyte ‘DDK’ proteins are neither obviously similar to the previously described KAI2, D14 or DLK2 protein types, nor indeed to each other. These clades also have long internal branch lengths, indicating a high degree of sequence divergence within the clades (Fig. 3). In the leptosporangiate fern core group there has been a duplication in the DDK lineage, and the resulting DDKA and DDKB protein types are strongly divergent both from each other and from other monilophyte DDK proteins. In seed plants, there are a number of major duplications and evidence for significant innovation in protein sequence (Fig. 3). In gymnosperms, we identified eu-D14 sequences that form a sister clade to the well-characterized angiosperm D14 clade. We also identified a second set of sequences in gymnosperms that are closely related to D14, which we named DWARF14-LIKE4 (DLK4). These form a sister clade to the gymnosperm/angiosperm eu-D14 clade, suggesting that the duplication that gave rise to DLK4 occurred before the separation of gymnosperms and angiosperms (Fig. 1). This in turn implies that the DLK4 clade has been lost from angiosperms (Fig. 3). Within the conifers there has been a major duplication in the DLK4 lineage giving rise to two sub-clades (DLK4A and DLK4B); since DLK4B is not found in Pinaceae, the separation of DLK4A and DLK4B seems to post-date the divergence of pines and other conifers (Fig. 3).
In angiosperms, we also discovered a third clade of proteins in addition to the expected D14 and DLK2 clades, which appeared as a sister clade to DLK2 in our analysis; we named these sequences DWARF14-LIKE3 (DLK3) (Fig. 1). Although our phylogenetic reconstruction suggests that the separation of DLK2 and DLK3 occurred before the radiation of extant angiosperms, the distribution of DLK3 sequences in our dataset suggests a slightly different history. We did not recover any DLK3-like sequences from the completed genome sequence of Amborella trichopoda (the sister group to all other angiosperms) or from the plants in the other early-diverging angiosperm orders (Nymphaeales, Austrobaileyales). We did identify unambiguous DLK3 sequences from the Chloroanthales and magnoliids, but not from any monocot species (including the fully sequenced genomes in Poaceae), despite extensive screening; we could however identify DLK2 sequences from across the monocot group. DLK3 sequences are present throughout the eudicots, although there have been sporadic losses, including some in the Brassicaceae. The exact inter-relationship of the major angiosperm lineages is currently uncertain, but one well-supported model is that monocots are sister to a clade containing magnoliids, Chloranthales and eudicots . Under this scenario, the distribution of genes suggests that the separation of the DLK2 and DLK3 lineages occurred after the divergence of monocots and other angiosperms (Fig. 3). Alternatively, DLK3 could have been lost from the monocot lineage. We also identified a group of gymnosperm proteins that form a sister group to the combined angiosperm DLK2-DLK3 clade, which we named DLK23. We also applied this name to the angiosperm proteins that pre-date the DLK2-DLK3 split, and to the wider seed plant clade containing all these proteins (Figs. 1 and 3).
Sequence conservation among D14/KAI2 proteins
To further understand the consequences of the evolutionary trajectories of the D14/KAI2 family members, we performed an in-depth analysis of their primary protein structure. Using our alignment, we identified a core set of 265 positions that occur in almost every D14/KAI2 protein (Fig. 4). The start and end positions of the polypeptide chain vary between individual sequences, but the majority of sequences are within the range –15 to 280. Extra amino acids are inserted within the core of the protein in some sequences; these are usually located outside secondary structural elements such as α-helices (Fig. 4). Most of these insertions are not conserved even between closely related sequences, although there are some exceptions. For instance, DDKB proteins from monilophytes have a conserved insertion of five amino acids after position 73.
Protein sequence conservation in D14/KAI2 proteins
Eu-KAI2 clade members have strong sequence conservation
Using this approach, we tested the hypothesis that evolution in the eu-KAI2 super-clade has generally been conservative. We analysed amino acid frequencies from 127 eu-KAI2 proteins and found that 22% of positions are invariant among eu-KAI2 proteins and 89% are conserved (Fig. 4, Table 2). By comparison, in the DDK super-clade only 5.6% of positions are invariant, with 63% conserved (Table 2). Indeed, the level of conservation across eu-KAI2 proteins as a whole is very comparable to conservation within taxon-level KAI2 clades. For instance, the angiosperm eu-KAI2 clade has 24% invariant positions and 94% conserved (Table 2). Together with the short branch lengths, the similarity in the level of between-clade and within-clade conservation in the eu-KAI2 super-clade supports the idea of a conservative evolutionary history.
Charophyte D14/KAI2 family members may encode proto-KAI2 proteins
Liverwort DDK clade members have conserved KAI2 structure
We next turned our attention to the DDK clade, which has lower overall amino acid conservation. We assessed whether the DDK proteins from liverworts (KAI2B), which have previously been characterized as KAI2-like, have conserved KAI2 features. We found that individual KAI2B proteins match the eu-KAI2 reference set at 29–33 out of 39 positions (Fig. 5a; Additional file 8). Although this is lower than eu-KAI2 proteins from liverworts (KAI2A), it suggests that these proteins could retain aspects of KAI2 primary protein structure. To test this idea, we generated homology models of liverwort KAI2B proteins using the crystal structure of karrikin-bound Arabidopsis thaliana KAI2 as a template . In each case, we found that the ligand binding pockets of KAI2B proteins are predicted to be essentially identical to those of Arabidopsis KAI2, and indeed liverwort KAI2A proteins (Fig. 6a–h; Additional files 9 and 10). Thus, while KAI2B proteins may be somewhat divergent relative to eu-KAI2 proteins, they probably still retain key features of eu-KAI2 structure.
Moss and lycophyte DDK clade members do not have KAI2- or D14-like sequences
Conversely, when we analysed the moss KAI2E and KAI2F proteins, we found that they only matched the KAI2 reference set at 22–24 positions (Fig. 5a, Additional file 8). This is a more considerable divergence from eu-KAI2 than liverwort KAI2B proteins and could imply a corresponding alteration in function. Indeed, structural modelling of the KAI2E/KAI2F proteins from P. patens has previously suggested that some of these proteins have altered ligand binding pockets relative to eu-KAI2 proteins in the same species . However, modelling of newly available KAI2E/F sequences from other mosses did not suggest major divergences from the KAI2 binding pocket (Fig. 6m–o).
Seed plant DLK23 and monilophyte DDK proteins may function independently of MAX2
Recent work has delineated the residues in D14 that are needed for interaction with MAX2-class F-box proteins . We confirm that these 18 residues are strongly conserved in D14 proteins, as suggested by . We also noted that 16 of those residues are very highly conserved in the eu-KAI2 super-clade, strongly suggesting that KAI2 proteins interact with MAX2 proteins through exactly the same interface as D14 (Additional file 11). However, the level of conservation across the D14/KAI2 family as a whole is considerably lower than in either the D14 or KAI2 groups, and we thus examined conservation of MAX2-interaction positions in other clades. Remarkably, we observed that of these 18 positions, 12 were not conserved in the highly divergent monilophyte DDKA and DDKB clades; 6 of these positions were not conserved in any monilophyte DDK protein (Additional file 11). Similarly, we found that 5, 7 and 8 of these positions are not conserved in DLK23, DLK2 and DLK3 proteins respectively (Additional file 11). Curiously, 5 of these positions are not conserved in the DLK4B, despite the MAX2 interface being otherwise conserved in the wider D14/DLK4 clade.
The monilophyte proteins occupy an intermediate position in the DDK clade, and we had therefore expected they would have protein sequences intermediate between the KAI2-like proteins in liverworts and eu-D14 proteins in seed plants. Individual monilophyte DDK proteins match the KAI2 reference set at 5–13 positions and the D14 reference set at 0–5 positions (Fig. 5), suggesting that, like lycophyte DDK proteins, they are not especially similar to characterized proteins, and have unique structural features. Indeed, homology modelling suggests that these proteins have quite variable ligand binding pockets that are generally larger than eu-KAI2 proteins but smaller than D14 proteins (Fig. 7e–h, Additional files 9 and 10). This is consistent with the general level of variation among monilophyte DDK proteins. Similarly, we observed that sequence conservation across the wider DLK23 clade is low; only 5% of positions are invariant, and only 60% conserved (Table 2). As would be expected, none of these proteins show affinity with KAI2 or D14 sequences (Fig. 5). It is therefore possible that loss of MAX2 interaction in these proteins has relaxed the structural requirements for protein function, resulting in divergent sequence characteristics.
The MAX2 family is highly conserved among land plants and charophyte algae
KAI2 signalling is highly conserved
Previous studies showed that proteins resembling KAI2 are found throughout land plants and in charophyte algae [23, 30, 34]. Consistent with this, we demonstrate that one of the two major clades in the land plant D14/KAI2 family contains only sequences that strongly resemble Arabidopsis KAI2. We demonstrate with very high resolution that these eu-KAI2 proteins are exceptionally conserved in protein sequence across the clade. Eu-KAI2 proteins have a clearly definable primary protein structure that is distinct from other members of the D14/KAI2 family, and their high levels of conservation arise from both shared-ancestral and shared-derived characteristics (Figs. 5 and 6). These data strongly suggest that there are very specific structural requirements for KAI2 function, and that these functional characteristics have been conserved throughout land plant evolution. Our results demonstrate that D14/KAI2 family proteins from charophytes do not quite meet the definition of eu-KAI2 proteins, but that they do have significant similarity with KAI2 proteins; we have thus categorized them as proto-KAI2. While the function and role of D14 in SL signalling are well understood, KAI2 proteins represent an enigma. In Arabidopsis, KAI2 is required for perception of karrikins, but has clearly defined developmental roles that are unrelated to karrikins; nor is Arabidopsis a naturally fire-following species [21, 23]. This has led to the hypothesis that KAI2 regulates development in response to an unknown endogenous ligand (KL), which is mimicked by karrikins [24, 25]. Consistent with an ancestral role of KL perception, expression of the eu-KAI2 protein from Selaginella moellendorffii (SmKAI2A) can partially rescue an Arabidopsis kai2 mutant but does not restore perception of karrikins . Identification of KL itself will be an important step in understanding the conserved function of KAI2 signalling across land plants .
An early origin for strigolactone signalling?
Previous analyses of the D14/KAI2 family have suggested that the origin of D14-type SL receptors is relatively recent, occurring within the vascular plant lineage, and perhaps restricted to seed plants [23, 30, 34]. Since SL sensitivity seems to be a widespread phenomenon in land plants and perhaps charophytes, this has led to significant speculation that non-canonical SL perception mechanisms exist in non-vascular plants [22, 33]. For instance, it has been suggested that KAI2 proteins could act as SL receptors in mosses and liverworts . Our analyses show that, as far as a distinct primary protein structure can be defined for eu-D14, such proteins do indeed only exist in seed plants. However, the separation of the DDK clade (of which eu-D14 proteins are members) from the eu-KAI2 clade occurred much earlier than previously suspected, at the base of the land plants. This raises the possibility that SL receptors might be a much earlier innovation in the D14/KAI2 family than previously suspected. The DDK protein from Selaginella moellendorffii (previously referred to as KAI2b) can hydrolyze SL-like stereoisomers of rac-GR24 , suggesting that it acts as an SL receptor. We show here that DDK proteins from lycophytes have little specific similarity to D14, which in turn suggests that other proteins in the clade could act as SL receptors despite their non-D14-like structure. However, understanding exactly when SL perception arose in the DDK lineage is contingent on understanding the evolution of land plants themselves. Although the phylogeny of vascular plants is well established, there is still considerable debate regarding the relationship of non-vascular plants, both to each other and to vascular plants. Depending on which scenario is correct, our understanding of the evolution of SL signalling may be considerably altered.
A more recent model of land plant evolution suggests that hornworts are the earliest-diverging group of land plants, and that liverworts and mosses form a clade that is sister to vascular plants (Fig. 9b) . The ‘hornworts-basal’ model is controversial, but consistent with it, we only identified a single clade of KAI2-like proteins from hornworts, which in some of our analyses place this clade as a sister clade to all other land plant D14/KAI2 sequences (Figs. 1 and 2). This would suggest that the duplication that created the eu-KAI2 and DDK lineages occurred after the separation of hornworts from all other land plants (Fig. 9b), although it should be noted that the recovery of a single hornwort clade could be due to the limitations of transcriptome databases. The close relationship of liverworts and mosses in this model (irrespective of their placement relative to hornworts) also has major implications for understanding the evolution of SL signalling. If this scenario is correct, then liverwort KAI2B and moss KAI2E/F are probable sister clades. Given the eu-KAI2 like structure of KAI2B protein, this would firmly imply that the ancestral state in the joint KAI2B-E/F clade would involve a KAI2-like binding pocket. If moss KAI2E/F proteins do indeed act as SL receptors, this would mean that SL-like binding pockets would have evolved twice independently in the DDK lineage, in mosses and vascular plants (Fig. 9b).
Our ability to precisely understand the origins of SL perception in the DDK lineage is thus currently limited by the lack of clarity regarding non-vascular plant phylogeny. It is nevertheless clear that the evolutionary trajectory of the DDK lineage is away from an initially KAI2-like structure, and that SL perception probably arose in the lineage at the latest in vascular plants. Given the high conservation between eu-KAI2 proteins, it is therefore very likely that that the majority of proteins in the DDK lineage are at least neo-functional with respect to KAI2. The primary question is thus whether they are neo-functional as SL receptors, or as something rather different. Our data suggest that the structural requirements for SL perception in vascular plants may be relatively relaxed, and even eu-D14 proteins only have limited shared-derived characteristics (Fig. 5a). We speculate that interactions with protein partners (such as SMXL proteins) may have driven the evolution of D14-like structure, rather than requirements for SL perception in itself.
MAX2-coupled signalling in the D14/KAI2 family
Alongside the origin of specific SL receptors, the evolution of SCFMAX2 coupling with D14/KAI2 signalling has also been a subject of debate. Two points have been emphasized; first, that proto-KAI2 proteins are present in charophyte algae, but that MAX2 homologues do not seem to be [30, 38]. Second, P. patens max2 mutants are reported to have a very different phenotype relative to P. patens SL synthesis mutants (no filamentous growth versus excessive filamentous growth), suggesting that they are not in the same pathway [32, 39]. On this basis, it has been suggested that SL signalling in non-vascular land plants might proceed by non-canonical mechanisms [7, 22]. Our data provide us with some insights in this respect. Firstly, the defined MAX2 interface found in D14 is highly conserved across most of the D14/KAI2 family, including in both eu-KAI2 and DDK proteins from liverworts, mosses and hornworts. It therefore seems likely that these proteins do indeed signal via MAX2 in non-vascular plants. We thus hypothesize that the reported max2 phenotype in P. patens arises from a lack of eu-KAI2 signalling, which in turn prevents expression of the SL-deficiency phenotype that would otherwise occur. Furthermore, our data show that the MAX2 interface is also conserved in charophyte D14/KAI2 proteins, tentatively suggesting the existence of MAX2-coupled signalling outside land plants. Consistent with this, we have identified an unambiguous MAX2-like protein in Coleochaete nitellarum. In contrast to the strong conservation of D14 and KAI2 proteins, we identified several clades of proteins (DLK2 and DLK3 from angiosperms, DLK23 and DLK4B from gymnosperms, DDKA/DDKB and probably all DDK proteins from monilophytes) that are strongly divergent at the positions that comprise the MAX2 interface. We find no evidence to suggest that these proteins might interact with specialized versions of MAX2. The MAX2 family has a very strongly conservative evolutionary history (Fig. 8), and there seems to have been a very strong selection pressure to maintain MAX2 as a single-copy gene. We only identified additional copies of MAX2 in a few angiosperm genomes known to have recent duplication events (e.g. Populus trichocarpa), and we did not identify additional or divergent copies of MAX2 in any non-angiosperm lineage. Thus, it seems highly likely these DDK proteins with non-conserved MAX2 interfaces signal independently of MAX2.
A diversity of small molecular receptors?
The DLK23 clade remains the most enigmatic set of proteins in the D14/KAI2 family. Not only do they probably lack the conserved MAX2 interface, but they are highly divergent from other D14/KAI2 proteins and have no known function. DLK2 in Arabidopsis does not seem to be a receptor for SL or KL, at least as far as can be defined genetically [21, 23]. One possibility is that the DLK23 proteins act as receptors for a novel ligand or perhaps multiple ligands. The DLK23 lineage in angiosperms has long internal branches, coupled with a lack of sequence conservation, but there is little evidence of gene loss. This suggests that the high degree of divergence does not simply represent drift in obsolete sequences. Rather, it may indicate continued innovation in the function of DLK23 proteins throughout angiosperm evolution, including the sub- or neo-functionalization process that led to independent DLK2 and DLK3 lineages. Since the DLK23 proteins from early-diverging angiosperms tend to group with eu-DLK2 species in phylogenetic analyses, this tentatively suggests that DLK2 maintained the original structure/function of DLK23, and that the DLK3 lineage is neo-functionalized. In addition to the DLK23 lineage, the fast-evolving DDK super-clade might contain further receptors for non-SL/KL ligands. For instance, since gymnosperms maintain conserved D14-type receptors, it is plausible that DLK4 proteins (and especially the more divergent DLK4B proteins) are not SL receptors. Our work broadens the structural biology platform for D14/KAI2 family members, and future work should provide very interesting insights into the ligand binding, structure and function of these diverse proteins, as well as their interactions with other SL signalling components.
We conclude that D14-like structure is not required for SL perception, and that SL perception has relatively relaxed structural requirements compared to KAI2-mediated signalling. We propose that SL perception gradually evolved by neo-functionalization within the DDK lineage, and that the transition from KAI2-like to D14-like protein may have been driven by interactions with protein partners, rather than being required for SL perception per se.
Bioinformatic retrieval of D14/KAI2 and MAX2 sequences
Members of the D14/KAI2 and MAX2 families were identified by BLAST searches against complete, annotated genomes from two major sources: Phytozome (https://phytozome.jgi.doe.gov/pz/portal.html) or the genome portals for individual species, for instance, the Amborella Genome Project (www.amborella.org). BLAST searches for D14/KAI2 sequences were performed using the full-length coding sequences of Arabidopsis thaliana D14, KAI2 and DLK2, using the BLASTN option. BLAST searches for MAX2 sequences were performed using the highly conserved C-terminal region of Arabidopsis thaliana MAX2, using the BLASTN option. Preliminary trees were assembled and used to guide the iterative interrogation of transcriptome databases, particularly those generated by the 1000 Plants (1KP) project (https://db.cngb.org/blast4onekp/). All sequences are listed in Additional file 13. For transcriptome datasets, we BLASTed each major taxonomic group separately. Where novel protein types were identified within a taxon (e.g. Angiosperm DLK3), we re-BLASTed the same taxonomic group with the novel sequence to increase the specificity of our searches. For non-annotated sequences from transcriptome datasets, we searched translations across all six reading frames to identify open reading frames (ORFs), and the longest ORFs were extracted for alignment.
Alignments were initially performed in BioEdit  using ClustalW . Full-length sequences from completed genomes were used for the initial alignment, which was manually refined as necessary. We then added sequences from transcriptome databases, many of which are incomplete, but the alignment of full-length sequences provided a scaffold to align these sequences correctly. For D14/KAI2 sequences the resultant alignment of 339 sequences is provided in Additional file 14. For MAX2 sequences, the resultant alignment of 57 sequences in provided in Additional file 15.
D14/KAI2 sequence analysis
For primary protein structure analyses, we focussed on positions in the alignment that are present in most sequences. We removed the non-conserved extensions at the N- and C-termini, producing an alignment with 265 core positions. We noted the positions of any non-conserved insertions within this core structure (Fig. 4) and then removed them prior to the final analyses. This 795-nucleotide alignment/265-amino acid alignment was used for analyses of primary protein structure (Figs. 4 and 5, Table 2, Additional file 11). Protein identity comparisons were performed in BioEdit using the ‘protein identity matrix’ function.
For the D14/KAI2 family we performed preliminary phylogenetic analyses to explore the topology of the tree and the effect of inclusion or exclusion of various groups of sequences. We removed 15 nucleotides (5 positions; 57–60, 252) from the 795-nucleotide alignment that were not well conserved across all sequences, leaving a ‘maximum’ phylogenetic alignment of 780 nucleotides. We implemented nucleotide-level maximum likelihood analyses in PhyML  and the Genetic Algorithm for Rapid Likelihood Inference (GARLI, v2.0) , using the GTR + G + I model of evolution. These analyses are generally congruent with subsequent analyses, but identified some problems with tree reconstruction, particularly with respect to the position of charophyte and lycophyte KAI2 sequences.
For final analyses, the alignment was manually modified in AliView v1.18-beta7 , and areas of ambiguous alignment were excluded from subsequent analyses. To determine the optimal model/partitioning scheme, we performed an exhaustive search in PartitionFinder v1.1.1 , with each of the three codon positions permitted its own parameters. All models were assessed, branch lengths were constrained to be proportional across partitions and the topology was fixed to that inferred by a preliminary GARLI v2.01 analysis with each codon position given its own GTR + I + G model and rates permitted to vary across partitions; the optimal scheme was selected by the Akaike information criterion (AIC) . Maximum likelihood tree searches were performed under this model (codon positions 1 and 2 with their own GTR + I + G sub-models and codon position 3 with a TVM + I + G sub-model; average rates permitted to vary across partitions) using GARLI v2.01, in the Cyberinfrastructure for Phylogenetic Research (CIPRES) Science Gateway . The GARLI tree searches were performed under the default settings with the exception that genthreshfortopoterm was increased to 40,000; these searches were performed from 48 different random addition sequence starting trees. Support was assessed with 528 bootstrap replicates in GARLI, under the same settings as the best-tree searches, but with each bootstrap search performed from 24 different random addition sequence starting trees. The resulting bootstrap support values were mapped onto our maximum likelihood phylogeny using the SumTrees v3.3.1 program in the DendroPy v3.12.0 package .
The MAX2 analyses followed the same workflow, with the exception that the optimal PartitionFinder scheme was for each codon position to have its own GTR + I + G model; searches for the best tree were performed from 10 different random addition starting trees, and 720 bootstrap searches were performed (each from two different random addition starting trees). Three sequences from 1KP accessions known to have contamination issues were pruned from the annotated tree using the Analyses of Phylogenetics and Evolution (APE) package in R [61, 62] prior to the preparation of Fig. 9.
Assessing tree robustness
We performed multiple analyses to test the robustness of our phylogenetic reconstructions, particularly the placement of KAI2B from liverworts and KAI2E/F from mosses within the DDK clade. Firstly, we removed each DDK clade from the alignment in turn, and re-ran the phylogenetic analysis in PhyML (Additional file 4). The 10 recovered trees have four commonalities: (1) KAI2B is always placed in the eu-KAI2 lineage (except in the ‘No KAI2’ tree), (2) the rest of the DDK clade is always stably grouped together (although there are some variations in the exact branching order within the clade), (3) the relative position of KAI2E/F is completely invariant (except in the ‘No KAI2E/F’ tree and (4) all of the trees place the eu-KAI2 lineage as a grade leading to the DDK clade. This latter point demonstrates that none of these trees are plausible in themselves, since the angiosperm eu-KAI2 clade is placed as a sister clade to the DDK clade containing moss, lycophyte, monilophyte, gymnosperms and angiosperm sequences. Secondly, we ran the analysis on an alignment cut down to match that of , using additional RbsQ (bacterial sigma factors with similarity to D14/KAI2 proteins) sequences identified in that study. If we rooted the resulting tree with RbsQ sequences, we observed the same basic topology as in . However, if we rooted with Selaginella moellendorffii KAI2, we obtained the same basic topology as in our main analyses, albeit with RbsQ as an in-group in the DDK lineage. Our analysis is thus congruent with the previous analysis in .
Protein homology modelling
KAI2 and DDK sequences were modelled using the SWISSMODEL server (http://swissmodel.expasy.org)  based on the ClustalW multiple sequence alignment of KAI2 and DDK sequences as described earlier in this manuscript. No further manipulation of the alignment was performed. Numerous KAI2 crystal structures were available for use as a model template [50, 64–66]; however, we chose the karrikin-bound A. thaliana structure (Protein Data Bank (PDB) code 4JYM)  as it was the most informative for probing the regions of the protein involved in ligand interaction. Modelled sequences share 37–71% sequence identity with A. thaliana KAI2 as computed by BioEdit  (Additional file 9). Protein structure and homology model figures were generated with PyMOL . Cavities within homology models were visualized using surface mode on the setting ’Cavities & Pockets (Culled)’ within PyMOL. Volume calculations were performed using the Computed Atlas of Surface Topography of proteins (CASTp) protein server  using a probe radius of 1.4 Å. Initial calculations of volume misleadingly included regions of the surface of the protein adjacent to the cavity opening. This problem was circumvented by artificially blocking the cavity opening with a free alanine residue which was not covalently attached to the protein molecule. This alanine was placed in the same xyz coordinates for all superposed homology models and crystal structures.
To independently confirm the robustness of the generated homology models, 10 representative sequences were also modelled using the I-TASSER server (http://zhanglab.ccmb.med.umich.edu/I-TASSER/) [69–71]. To confirm that both methods generated similar models, root-mean-square deviation (RMSD) values for the SWISSMODEL- and I-TASSER-generated models were then calculated for the 10 pairs of sequences using the SuperPose server based on the Cα coordinates (http://wishart.biology.ualberta.ca/superpose/) . The RMSD values confirmed that both model-generating servers converged on essentially the same result (Additional file 16).
We gratefully acknowledge the use of sequence data generated by members of the 1000 Plants (1KP) initiative, and in particular Michael Melkonian, Barbara Surek, Jim Leebens-Mack, Michael Deyholos, Douglas Soltis, Pamela Soltis, Anders Larsson, Lisa Pokory and Lisa DeGironimo.
TB is supported by RG160845 from the Royal Society. DS is supported by the Monell Foundation.
Availability of data and materials
All data generated or analysed during this study are included in this published article and its additional files.
TB and DN conceived and designed the study. CR, DS, SG and GK generated and provided sequence data. TB analysed sequence data and performed phylogenetic reconstructions with CR. RBD performed homology modelling. TB wrote the manuscript with input from all authors. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Al-Babili S, Bouwmeester HJ. Strigolactones, a novel carotenoid-derived plant hormone. Annu Rev Plant Biol. 2015;66:161–86.View ArticlePubMedGoogle Scholar
- López-Ráez JA, Charnikhova T, Gómez-Roldán V, Matusova R, Kohlen W, De Vos R, et al. Tomato strigolactones are derived from carotenoids and their biosynthesis is promoted by phosphate starvation. New Phytol. 2008;178:863–74.View ArticlePubMedGoogle Scholar
- Borghi L, Liu GW, Emonet A, Kretzschmar T, Martinoia E. The importance of strigolactone transport regulation for symbiotic signalling and shoot branching. Planta. 2016;243:1351–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Matthys C, Walton A, Struk S, Stes E, Boyer FD, Gevaert K, Goormachtig S. The whats, the wheres and the hows of strigolactone action in the roots. Planta. 2016;243:1327–37.View ArticlePubMedGoogle Scholar
- Kohlen W, Charnikhova T, Liu Q, Bours R, Domagalska MA, Beguerie S, et al. Strigolactones are transported through the xylem and play a key role in shoot architectural response to phosphate deficiency in nonarbuscular mycorrhizal host Arabidopsis. Plant Physiol. 2011;155:974–87.View ArticlePubMedGoogle Scholar
- Smith SM, Waters MT. Strigolactones: destruction-dependent perception? Curr Biol. 2012;22:R924–927.View ArticlePubMedGoogle Scholar
- Waters MT, Gutjahr C, Bennett T, Nelson D. Strigolactone signalling and evolution. Ann Rev Plant Biol. 2017;68:8.1–8.31.View ArticleGoogle Scholar
- de Saint GA, Clavé G, Badet-Denisot MA, Pillot JP, Cornu D, Le Caer JP, et al. An histidine covalent receptor and butenolide complex mediates strigolactone perception. Nat Chem Biol. 2016;12:787–94.View ArticleGoogle Scholar
- Yao R, Ming Z, Yan L, Li S, Wang F, Ma S, et al. DWARF14 is a non-canonical hormone receptor for strigolactone. Nature. 2016;536:469–73.View ArticlePubMedGoogle Scholar
- Stirnberg P, van De Sande K, Leyser HM. MAX1 and MAX2 control shoot lateral branching in Arabidopsis. Development. 2002;129:1131–41.PubMedGoogle Scholar
- Stirnberg P, Furner IJ, Leyser HMO. MAX2 participates in an SCF complex which acts locally at the node to suppress shoot branching. Plant J. 2007;50:80–94.View ArticlePubMedGoogle Scholar
- Hamiaux C, Drummond RS, Janssen BJ, Ledger SE, Cooney JM, Newcomb RD, Snowden KC. DAD2 is an alpha/beta hydrolase likely to be involved in the perception of the plant branching hormone, strigolactone. Curr Biol. 2012;22:2032–6.View ArticlePubMedGoogle Scholar
- Zhao LH, Zhou XE, Yi W, Wu Z, Liu Y, Kang Y, et al. Destabilization of strigolactone receptor DWARF14 by binding of ligand and E3-ligase signalling effector DWARF3. Cell Res. 2015;25:1219–36.View ArticlePubMedPubMed CentralGoogle Scholar
- Nakamura H, Xue YL, Miyakawa T, Hou F, Qin HM, Fukui K, et al. Molecular mechanism of strigolactone perception by DWARF14. Nat Commun. 2013;4:2613.PubMedGoogle Scholar
- Wang Y, Sun S, Zhu W, Jia K, Yang H, Wang X. Strigolactone/MAX2-induced degradation of brassinosteroid transcriptional effector BES1 regulates shoot branching. Dev Cell. 2013;27:681–8.View ArticlePubMedGoogle Scholar
- Jiang L, Liu X, Xiong G, Liu H, Chen F, Wang L, et al. DWARF 53 acts as a repressor of strigolactone signalling in rice. Nature. 2013;504:401–5.View ArticlePubMedGoogle Scholar
- Zhou F, Lin Q, Zhu L, Ren Y, Zhou K, Shabek N, et al. D14-SCF(D3)-dependent degradation of D53 regulates strigolactone signalling. Nature. 2013;504:406–10.View ArticlePubMedPubMed CentralGoogle Scholar
- Soundappan I, Bennett T, Morffy N, Liang Y, Stanga JP, Abbas A, et al. SMAX1-LIKE/D53 family members enable distinct MAX2-dependent responses to strigolactones and karrikins in Arabidopsis. Plant Cell. 2015;27:3143–59.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang L, Wang B, Jiang L, Liu X, Li X, Lu Z, et al. Strigolactone signalling in Arabidopsis regulates shoot development by targeting D53-like SMXL repressor proteins for ubiquitination and degradation. Plant Cell. 2015;27:3128–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Liang Y, Ward S, Li P, Bennett T, Leyser O. SMAX1-LIKE7 signals from the nucleus to regulate shoot development in Arabidopsis via partially EAR motif-independent mechanisms. Plant Cell. 2016;28:1581–601.PubMedPubMed CentralGoogle Scholar
- Bennett T, Liang Y, Seale M, Ward S, Müller D, Leyser O. Strigolactone regulates shoot development through a core signalling pathway. Biol Open. 2016;5:1806–20.PubMedPubMed CentralGoogle Scholar
- Bennett T, Leyser O. Strigolactone signalling: standing on the shoulders of DWARFs. Curr Opin Plant Biol. 2014;22:7–13.View ArticlePubMedGoogle Scholar
- Waters MT, Nelson DC, Scaffidi A, Flematti GR, Sun YK, Dixon KW, Smith SM. Specialisation within the DWARF14 protein family confers distinct responses to karrikins and strigolactones in Arabidopsis. Development. 2012;139:1285–95.View ArticlePubMedGoogle Scholar
- Flematti GR, Waters MT, Scaffidi A, Merritt DJ, Ghisalberti EL, Dixon KW, Smith SM. Karrikin and cyanohydrin smoke signals provide clues to new endogenous plant signalling compounds. Mol Plant. 2013;6:29–37.View ArticlePubMedGoogle Scholar
- Conn CE, Nelson DC. Evidence that KARRIKIN-INSENSITIVE2 (KAI2) receptors may perceive an unknown signal that is not karrikin or strigolactone. Front Plant Sci. 2016;6:1219.View ArticlePubMedPubMed CentralGoogle Scholar
- Gutjahr C, Gobbato E, Choi J, Riemann M, Johnston MG, Summers W, et al. Rice perception of symbiotic arbuscular mycorrhizal fungi requires the karrikin receptor complex. Science. 2015;350:1521–4.View ArticlePubMedGoogle Scholar
- Nelson DC, Scaffidi A, Dun EA, Waters MT, Flematti GR, Dixon KW, et al. F-box protein MAX2 has dual roles in karrikin and strigolactone signalling in Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2011;108:8897–902.View ArticlePubMedPubMed CentralGoogle Scholar
- Stanga JP, Smith SM, Briggs WR, Nelson DC. SUPPRESSOR OF MORE AXILLARY GROWTH2 1 controls seed germination and seedling development in Arabidopsis. Plant Physiol. 2013;163:318–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Stanga JP, Morffy N, Nelson DC. Functional redundancy in the control of seedling growth by the karrikin signalling pathway. Planta. 2016;243:1397–406.View ArticlePubMedGoogle Scholar
- Delaux PM, Xie X, Timme RE, Puech-Pages V, Dunand C, Lecompte E, et al. Origin of strigolactones in the green lineage. New Phytol. 2012;195:857–71.View ArticlePubMedGoogle Scholar
- Wickett NJ, Mirarab S, Nguyen N, Warnow T, Carpenter E, Matasci N, et al. Phylotranscriptomic analysis of the origin and early diversification of land plants. Proc Natl Acad Sci U S A. 2014;111:E4859–68.View ArticlePubMedPubMed CentralGoogle Scholar
- Proust H, Hoffmann B, Xie X, Yoneyama K, Schaefer DG, Yoneyama K, et al. Strigolactones regulate protonema branching and act as a quorum sensing-like signal in the moss Physcomitrella patens. Development. 2011;138:1531–9.View ArticlePubMedGoogle Scholar
- Waldie T, McCulloch H, Leyser O. Strigolactones and the control of plant development: lessons from shoot branching. Plant J. 2014;79:607–22.View ArticlePubMedGoogle Scholar
- Waters MT, Scaffidi A, Moulin SL, Sun YK, Flematti GR, Smith SM. A Selaginella moellendorffii ortholog of KARRIKIN INSENSITIVE2 functions in Arabidopsis development but cannot mediate responses to karrikins or strigolactones. Plant Cell. 2015;27:1925–44.View ArticlePubMedPubMed CentralGoogle Scholar
- Conn CE, Bythell-Douglas R, Neumann D, Yoshida S, Whittington B, Westwood JH, et al. Convergent evolution of strigolactone perception enabled host detection in parasitic plants. Science. 2015;349:540–3.View ArticlePubMedGoogle Scholar
- Toh S, Holbrook-Smith D, Stogios PJ, Onopriyenko O, Lumba S, Tsuchiya Y, et al. Structure-function analysis identifies highly sensitive strigolactone receptors in Striga. Science. 2015;350:203–7.View ArticlePubMedGoogle Scholar
- Tsuchiya Y, Yoshimura M, Sato Y, Kuwata K, Toh S, Holbrook-Smith D, et al. Probing strigolactone receptors in Striga hermonthica with fluorescence. Science. 2015;349:864–8.View ArticlePubMedGoogle Scholar
- Challis RJ, Hepworth J, Mouchel C, Waites R, Leyser O. A role for more axillary growth1 (MAX1) in evolutionary diversity in strigolactone signalling upstream of MAX2. Plant Physiol. 2013;161:1885–902.View ArticlePubMedPubMed CentralGoogle Scholar
- de Saint GA, Bonhomme S, Boyer FD, Rameau C. Novel insights into strigolactone distribution and signalling. Curr Opin Plant Biol. 2013;16:583–9.View ArticleGoogle Scholar
- Flores-Sandoval E, Eklund DM, Bowman JL. A simple auxin transcriptional response system regulates multiple morphogenetic processes in the liverwort Marchantia polymorpha. PLoS Genet. 2015;11:e1005207.View ArticlePubMedPubMed CentralGoogle Scholar
- Kato H, Ishizaki K, Kouno M, Shirakawa M, Bowman JL, Nishihama R, Kohchi T. Auxin-mediated transcriptional system with a minimal set of components is critical for morphogenesis through the life cycle in Marchantia polymorpha. PLoS Genet. 2015;11:e1005084.View ArticlePubMedPubMed CentralGoogle Scholar
- Lavy M, Prigge MJ, Tao S, Shain S, Kuo A, Kirchsteiger K, Estelle M. Constitutive auxin response in Physcomitrella reveals complex interactions between Aux/IAA and ARF proteins. Elife. 2016;5:e13325.View ArticlePubMedPubMed CentralGoogle Scholar
- Wodniok S, Brinkmann H, Glöckner G, Heidel AJ, Philippe H, Melkonian M, Becker B. Origin of land plants: do conjugating green algae hold the key? BMC Evol Biol. 2011;11:104.View ArticlePubMedPubMed CentralGoogle Scholar
- Timme RE, Bachvaroff TR, Delwiche CF. Broad phylogenomic sampling and the sister lineage of land plants. PLoS One. 2012;7:e29696.View ArticlePubMedPubMed CentralGoogle Scholar
- Doyle JA. Phylogenetic analyses and morphological innovations in land plants. In: Ambrose BA, Purugganan M, editors. The evolution of plant form. Hoboken: Blackwell; 2013. p. 1–50.Google Scholar
- Zwickl DJ. Genetic algorithm approaches for the phylogenetic analysis of large biological sequence datasets under the maximum likelihood criterion. Austin: The University of Texas Austin; 2006.Google Scholar
- Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Systematic Biol. 2010;59:307–21.View ArticleGoogle Scholar
- Cox CJ, Li B, Foster PG, Embley TM, Civán P. Conflicting phylogenies for early land plants are caused by composition biases among synonymous substitutions. Syst Biol. 2014;63:272–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Lopez-Obando M, Conn CE, Hoffmann B, Bythell-Douglas R, Nelson DC, Rameau C, Bonhomme S. Structural modelling and transcriptional responses highlight a clade of PpKAI2-LIKE genes as candidate receptors for strigolactones in Physcomitrella patens. Planta. 2016;243:1441–53.View ArticlePubMedGoogle Scholar
- Guo Y, Zheng Z, La Clair JJ, Chory J, Noel JP. Smoke-derived karrikin perception by the α/β-hydrolase KAI2 from Arabidopsis. Proc Natl Acad Sci U S A. 2013;110:8284–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang C, Liu Y, Li SS, Han GZ. Insights into the origin and evolution of the plant hormone signaling machinery. Plant Physiol. 2015;167:872–86.View ArticlePubMedPubMed CentralGoogle Scholar
- Sun YK, Flematti GR, Smith SM, Waters MT. Reporter gene-facilitated detection of compounds in Arabidopsis leaf extracts that activate the karrikin signalling pathway. Front Plant Sci. 2016;7:1799.PubMedPubMed CentralGoogle Scholar
- Qiu Y-L, Li L, Wang B, Chen Z, Knoop V, Groth-Malonek M, et al. The deepest divergences in land plants inferred from phylogenomic evidence. Proc Natl Acad Sci U S A. 2006;103:15511–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Hall T. BioEdit: biological sequence alignment editor. Carlsbad: Ibis Biosciences; 1997.Google Scholar
- Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22:4673–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Larsson A. AliView: a fast and lightweight alignment viewer and editor for large data sets. Bioinformatics. 2014;30:3276–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Lanfear R, Calcott B, Ho SYW, Guindon S. PartitionFinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol. 2012;29:1695–701.View ArticlePubMedGoogle Scholar
- Akaike H. A new look at the statistical model identification. IEEE Trans Autom Control. 1974;19:716–23.View ArticleGoogle Scholar
- Miller MA, Pfeiffer W, Schwartz T editors. Proceedings of the Gateway Computing Environments Workshop (GCE). New Orleans, LA. 2010.Google Scholar
- Sukumaran J, Holder MT. DendroPy: a Python library for phylogenetic computing. Bioinformatics. 2010;26:1569–71.View ArticlePubMedGoogle Scholar
- Paradis E, Claude J, Strimmer K. APE: Analyses of Phylogenetics and Evolution in R language. Bioinformatics. 2004;20:289–90. doi:https://doi.org/10.1093/bioinformatics/btg412.View ArticlePubMedGoogle Scholar
- Core TR. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2013.Google Scholar
- Bordoli L, Kiefer F, Arnold K, Benkert P, Battey J, Schwede T. Protein structure homology modeling using SWISS-MODEL workspace. Nat Protoc. 2009;4:1–13.View ArticlePubMedGoogle Scholar
- Bythell-Douglas R, Waters MT, Scaffidi A, Flematti GR, Smith SM, Bond CS. The structure of the karrikin-insensitive protein (KAI2) in Arabidopsis thaliana. PLoS One. 2013;8:e54758.View ArticlePubMedPubMed CentralGoogle Scholar
- Kagiyama M, Hirano Y, Mori T, Kim SY, Kyozuka J, Seto Y, Yamaguchi S, Hakoshima T. Structures of D14 and D14L in the strigolactone and karrikin signalling pathways. Genes Cells. 2013;18:147–60.View ArticlePubMedGoogle Scholar
- Zhao LH, Zhou XE, Wu ZS, Yi W, Xu Y, Li S, et al. Crystal structures of two phytohormone signal-transducing α/β hydrolases: karrikin-signalling KAI2 and strigolactone-signalling DWARF14. Cell Res. 2013;23:436–9.View ArticlePubMedPubMed CentralGoogle Scholar
- De Lano W. The Pymol molecular graphics system. San Carlos: Schrödinger, LLC; 2002.Google Scholar
- Dundas J, Ouyang Z, Tseng J, Binkowski A, Turpaz Y, Liang J. CASTp: Computed Atlas of Surface Topography of proteins with structural and topographical mapping of functionally annotated residues. Nucleic Acids Res. 2006;34:W116–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang J, Yan R, Roy A, Xu D, Poisson J, Zhang Y. The I-TASSER Suite: protein structure and function prediction. Nat Methods. 2015;12:7–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Roy A, Kucukural A, Zhang Y. I-TASSER: a unified platform for automated protein structure and function prediction. Nat Protoc. 2010;5:725–38.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang Y. I-TASSER server for protein 3D structure prediction. BMC Bioinf. 2008;9:40.View ArticleGoogle Scholar
- Maiti R, Domselaar GHV, Zhang H, Wishart DS. SuperPose: a simple server for sophisticated structural superposition. Nucleic Acids Res. 2004;32:W590–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Xu Y, Miyakawa T, Nakamura H, Nakamura A, Imamura Y, Asami T, et al. Structural basis of unique ligand specificity of KAI2-like protein from parasitic weed Striga hermonthica. Sci Rep. 2016;6:31386.Google Scholar