#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Evolution of the Auxin Response Factors from charophyte ancestors


Authors: Raquel Martin-Arevalillo aff001;  Emmanuel Thévenon aff002;  Fanny Jégu aff002;  Thomas Vinos-Poyo aff002;  Teva Vernoux aff001;  François Parcy aff002;  Renaud Dumas aff002
Authors place of work: Laboratoire de Reproduction et Développement des Plantes, Univ. Lyon, ENS de Lyon, UCB Lyon1, CNRS, INRA, Lyon, France aff001;  Univ. Grenoble Alpes, CNRS, CEA, INRA, IRIG-DBSCI-LPCV, Grenoble, France aff002
Published in the journal: Evolution of the Auxin Response Factors from charophyte ancestors. PLoS Genet 15(9): e32767. doi:10.1371/journal.pgen.1008400
Category: Research Article
doi: https://doi.org/10.1371/journal.pgen.1008400

Summary

Auxin is a major developmental regulator in plants and the acquisition of a transcriptional response to auxin likely contributed to developmental innovations at the time of water-to-land transition. Auxin Response Factors (ARFs) Transcription Factors (TFs) that mediate auxin-dependent transcriptional changes are divided into A, B and C evolutive classes in land plants. The origin and nature of the first ARF proteins in algae is still debated. Here, we identify the most ‘ancient’ ARF homologue to date in the early divergent charophyte algae Chlorokybus atmophyticus, CaARF. Structural modelling combined with biochemical studies showed that CaARF already shares many features with modern ARFs: it is capable of oligomerization, interacts with the TOPLESS co-repressor and specifically binds Auxin Response Elements as dimer. In addition, CaARF possesses a DNA-binding specificity that differs from class A and B ARFs and that was maintained in class C ARF along plants evolution. Phylogenetic evidence together with CaARF biochemical properties indicate that the different classes of ARFs likely arose from an ancestral proto-ARF protein with class C-like features. The foundation of auxin signalling would have thus happened from a pre-existing hormone-independent transcriptional regulation together with the emergence of a functional hormone perception complex.

Keywords:

DNA-binding proteins – Sequence alignment – Sequence motif analysis – Algae – Sequence databases – Arabidopsis thaliana – Auxins – Plant evolution

Introduction

Charophytes diverged from chlorophyte algae more than a billion years ago (y.a.) and led to land plants emergence around 450 million y.a. [26]. “Early divergent” clades display a range of body complexity going from unicellular algae in Mesostigmatophyceae and Chlorokybophyceae, to multicellular filaments in Klebsormidiophyceae (Fig 1; S1 Fig) [7,8]. “Late divergent” clades include Charophyceae and Coleochaetophyceae that share features with land plants (S1 Fig), [9,10] but also Zygnematophyceae, that despite their simple structure are considered sisters to land plants according to recent phylogenetic studies [11,12].

Fig. 1. Charophytes B3- and PB1-domain containing proteins along evolution of charophytes algae and land plants.
Charophytes B3- and PB1-domain containing proteins along evolution of charophytes algae and land plants.
On the left, an illustrated view of charophytes evolutionary line from chlorophytes to land plants adapted from [10,42]. Although some publications placed Mesostigmaphoticae and Chlorokybophycae in a unified clade large differences in sequences and morphological traits argues for two different clades [3,7,42,43]. The presence of class C, A/B, A or B ARFs is indicated by blue, orange, red or green circles whereas yellow circles correspond to RAV proteins. Identities of proteins reported in the figure are shown in Supplementary S2 Table.

Given the importance of the phytohormone auxin in plant development, the acquisition of its signalling pathway allowing for auxin-dependent changes in transcription is thought to have been a milestone in the water-to-land transition [2]. In land plants, this signalling pathway, called the Nuclear Auxin Pathway (NAP), relies on three main protein families: TIR1/AFB (Transport Inhibitor Response 1/Auxin Signalling F-box) co-receptors, Aux/IAA transcriptional repressors (Auxin/Indole-3-Acetic Acid Protein) and ARF (Auxin Response Factors) Transcription Factors (TFs) [13,14]. ARFs have been classified into three evolutive classes, A, B and C. Class A includes activator ARFs whereas classes B and C contain repressor ARFs [1]. ARF interaction with DNA is mediated by their B3 domain (B3ARF). Such domain is also present in ABI3 (Abscisic Acid insensitive 3, B3ABI3) and RAV (Related to ABI/VP1, B3RAV) plant TFs but with different DNA binding specificities [15,16]. In the ARF family, the B3 domain is embedded in a larger N-terminal DNA Binding Domain (DBD) that includes a Dimerization Domain (DD). As dimers, ARFs bind double AuxREs (Auxin Response Elements) sites arranged in three possible orientations: Direct Repeat (DR), Everted Repeat (ER) and Inverted Repeat (IR) (S2 Fig) [2,1719].

In charophyte algae and the bryophyte Marchantia polymorpha, the B3RAV and B3ARF domains are often associated with a C-terminal PB1 oligomerization domain, a landmark of most ARF TFs in higher plants but that was lost by RAV TFs from tracheophytes [2,20]. This shared B3 + PB1 domain composition led to the initial hypothesis that ARF could have arisen from RAV [21]. In the NAP, the PB1 domain mediates ARF homo-oligomerization and hetero-oligomerization with Aux/IAAs [22]. Under low auxin concentrations, Aux/IAAs bind activator ARFs through their PB1 domain [2326] and recruit TOPLESS (TPL)/TOPLESS-RELATED co-repressors, leading to the formation of a repressor complex on regulatory sequences of auxin-responsive genes [27]. When auxin levels increase, the hormone-mediated interaction between Aux/IAA and TIR1/AFB leads to Aux/IAA proteasomal degradation, unlocking activator ARFs and inducing transcription [28,29]. Most class B and C ARF members have limited interaction capacities with Aux/IAAs [3032] and are proposed to regulate auxin transcriptional responses in an auxin-independent manner, possibly by competitive binding with class A ARFs on DNA regulatory sequences [33,34]. Since some of class B and C ARFs can interact directly with TPL, formation of co-repressor complexes was proposed as another possible mechanism for transcriptional repression of auxin target genes [3436].

The NAP was established at the beginning of land plants history. In the bryophyte M. polymorpha for example, the 3 families of NAP proteins are present (with one member of each ARF class) as well as the TPL co-repressor [37,38]. Recent studies showed the existence of two ARF subfamilies in charophytes, class C and class A/B [20], but the absence of functional TIR1/AFB and Aux/IAAs suggested that a fully functional NAP did not exist before land plants [2,20,37,3941]. How these ancestral components evolved to form the land plants NAP remains an open question. Through the structural, biochemical and phylogenetic characterisation of a proto-ARF from an early divergent charophyte we set a scenario of how the co-option of ancestral mechanisms of transcriptional control possibly led to the evolution of hormone signalling pathways in plants.

Results

Identification of proto-ARF and proto-RAV in early divergent charophytes

To understand the evolution of ARFs, we first characterized the biochemical properties of proto-ARFs and closely related proto-RAVs from early divergent charophytes. We searched for B3 homologues in charophyte transcripts databases (OneKp and Marchantia.info) [12,44] and classified them as B3RAV or B3ARF, depending on the residues signature of their predicted DBDs (S1 Table) [45]. B3RAV domains were frequently associated with an APETALA2 (AP2) domain and/or PB1 domains in the basal charophyte Mesostigma viride and all later clades (Fig 1; S2 Table). M. viride also has an ARF homologue (GBSK01006108.1) devoid of a PB1 domain [2]. Its DBD was reliably modelled as an ARF (100% confidence with AtARF1 [46,47]), but it lacks most residues involved in the interaction with AuxREs (S3 Fig) and thus does not qualify as a functional ARF. The proto-ARF of the earliest diverging green Charophyte algae with predicted functional B3ARF and PB1 domains was found in C. atmophyticus. Other ARF homologues were also present in all later diverging clades (Fig 1; S3 Fig).

DNA binding specificities and oligomerization potential of proto-RAV and proto-ARF

We determined the properties of “ancestral” RAV and ARF proteins, focusing on K. nitens proto-RAV (containing predicted AP2, B3RAV and PB1) (KnRAV, kfl00094_0070) and C. atmophyticus proto-ARF (CaARF, AZZW-2021616). The predicted B3 domains of KnRAV and CaARF display the signature residues typical of B3RAV and B3ARF, respectively (S1 Table; S3 and S4 Figs) suggesting that their divergent DNA binding specificities were already established in charophytes. To test this hypothesis, we characterized the binding of their DBD against the canonical DNA binding sites identified in angiosperms for ABI3, RAV and ARF TFs. KnRAV specifically bound the AP2/B3RAV bipartite element described for Arabidopsis thaliana RAV TFs (Fig 2A) [48]. CaARF interacted strongly with double AuxRE sites (DR or ER, Fig 2B) but not with a single AuxRE site suggesting that the DBD of CaARF binds DR and ER motifs as a dimer without the help of the Middle Region (MR) and the PB1 domain. Altogether, these results confirm that RAV and ARF DNA binding preferences were established in basal charophytes and maintained along evolution.

Fig. 2. Proto-RAV and proto-ARF biochemical properties.
Proto-RAV and proto-ARF biochemical properties.
(A), KnRAV-DBD (AP2/B3 domains) or (B), Ca-ARF-DBD interactions in EMSA with specific single binding sites for different B3 families (ABI3, ARF and RAV) [45], AP2, a composite site for AP2-B3/RAV family, AP2(-)/RAV mutated on AP2, AP2/RAV(-) mutated on RAV, AP2(-)/RAV(-) mutated on both AP2 and RAV, DR5 (Direct Repeats AuxRE spaced by 5 nucleotides) and ER7 (Everted Repeats spaced by 8 nucleotides) ARF motifs. DNA sequences are indicated in S8 Table. Proteins added at 0 and 0.5 μM. (C-D), Structure models for KnRAV-PB1 (C, green) and CaARF-PB1 (D, blue) superposed to AtARF5-PB1 structure (cyan) (PDB code 4CHK [23]). Conserved residues indicated refer to AtARF5-PB1 structure. Positive and negative signs indicate potential interaction surfaces for oligomerization. (E), SEC-MALLS KnRAV-PB1 molecular weight determination for four protein concentrations (from dark to light green: 5 mg/mL, 2.5 mg/mL, 1.25 mg/mL and 0.625 mg/mL). (F), SEC-MALLs CaARF-PB1 molecular weight determination (from dark to light blue: 5 mg/mL, 2.5 mg/mL, 1.25 mg/mL and 0.625 mg/mL.

Next, we studied the oligomerization capacity of their PB1 domain. Based on AtARF5 PB1 structure [23,47], the PB1 domains of KnRAV and CaARF were modelled as type I/II PB1 with electrostatic oligomerization potential (Fig 2C and 2D). Molecular weight determination of KnRAV-PB1 and CaARF-PB1 by Size Exclusion Chromatography combined with Multi-Angle-Light Scattering (SEC-MALLS) experimentally validated that both domains form oligomeric complexes (Fig 2E and 2F) but with a lower oligomeric potential than AtARF5-PB1 (S3 Table). Charophycean algae therefore appear to possess proto-RAV and proto-ARF proteins with oligomerization potential and diverging DNA binding specificities (Fig 2A and 2B; S3 and S4 Figs).

Evolution of ARF DNA binding specificity from early divergent charophytes to land plants

To further characterize the biophysical properties of proto-ARFs, we determined the predicted structure of CaARF DBD and showed that it was reliably modeled (99% confidence; Phyre 2) with AtARF1 and AtARF5 DBDs [46,47] except for an additional disordered region in CaARF present within the DD (Fig 3A). Similar disordered regions were found as a characteristic feature of all class C ARFs (Fig 3B and 3C; S3 and S7 Figs). In agreement with this, our phylogenetic studies position CaARF within clade C (S5 Fig). Such insertions are expected to modify class C DNA binding compared to A and B ARFs. We tested this hypothesis using ER motifs with different spacing (ER4-9). Unlike Arabidopsis AtARF2 (class B) and AtARF5 (class A) that largely prefer ER7/8 motifs (Fig 3D and 3E), CaARF showed promiscuous binding to ER4-9 but did not interact with a single AuxRE motif (Fig 3D–3F) confirming its interaction with ER motifs as a dimer. Arabidopsis class C AtARF10 behaves similarly to CaARF (Fig 3G). This shows that CaARF has a relaxed DNA specificity allowing binding to ER binding sites with various distances between the monomeric sites and that this specificity was maintained in class C ARF along plants evolution. The presence of a specific disordered region (Fig 3A–3C; S3 and S7 Figs) in class C ARF DBDs suggests a possible role in their relaxed specificity, that remains to be tested.

Fig. 3. Ancestral class C ARFs exhibit different DNA specificities and TPL recruitment potential.
Ancestral class C ARFs exhibit different DNA specificities and TPL recruitment potential.
(A), CaARF-DBD modelled structure (dark blue) superposed to ARF1 dimer (cyan with DD in yellow and B3 in red indicated on one monomer (4LDX, [46])). (B-C), Class-C MpARF3 (orange, B) and AtARF10 (green, C) DBD structure models superposed to ARF1 (cyan). Yellow dots indicate the position of the additional disordered region, not modelled. (D-G), AtARF2 (D), AtARF5 (E), CaARF (F) and AtARF10 (G) binding to ER DNA binding sites with a spacing ranging from 4 to 9 nucleotides. Positive signs indicate lanes where protein was added (AtARF2, AtARF5 and AtARF10 at 0.25 μM; CaARF at 1 μM). (H), C. atmophyticus predicted TPL (AZZW-2021890) N-ter structure model (dark blue) superposed with AtTPL N-ter structure (cyan) (5NQV, [49]). (I), Co-purification interaction assays between MBP-tagged CaARF (and mutants: CaARFmL523S/F524S; CaARFmL492A/L493A) and AtTPL202 (and mutant AtTPL202mF74A). Complexes were bound to Dextrin-Sepharose columns. SDS-PAGEs show the proteins eluted from the column. AtTPL202* is a control indicating the size of the TPL N-ter protein.

Interaction with co-repressors in early divergent charophytes

As mentioned before, certain land plants ARF proteins have the capacity to interact directly or indirectly with the TPL co-repressor [35,36,50]. We wondered when in evolution this interaction was first established. Direct TPL-recruitment usually involves two different amino acid regions in the Middle Region (MR) of repressor ARFs: the EAR-motif (ERF-associated Amphiphilic Repression motif with LxLxL sequence or its LxLxPP variant) and the BRD domain (B3 Repression Domain with the K/RLFG sequence) [35,36], the BRD domain also being found in RAV proteins. CaARF-MR presents two potential repression regions with an EAR-like motif (LPLLPS, similar to LxLxPP) and a BRD domain (KLFG). Since TPL EAR-interacting-region (TPL N-terminal, TPL-N) is extremely conserved between charophytes and land plants [49,51] (Fig 3H; S8 Fig; S4 Table), we used A. thaliana TPL-N (AtTPL202) to assay the TPL/CaARF interaction. CaARF interacted with AtTPL202 in co-purification assays and this interaction was lost with AtTPL202-F74A, mutated in the hydrophobic EAR peptide binding groove (Fig 3I) [49]. Moreover, mutations in CaARF KLFG (CaARF-L523S/F524S) or LPLLPS (CaARF-L492A/L493A) weakened the interaction with AtTPL202, indicating that both sites might participate to TPL-N recruitment. The binding of the BRD domain of CaARF differs from that of the RAV1 of A. thaliana which interacts with the C-terminal part of TPL [52], suggesting different TPL recruitment mechanisms for these two protein families. The presence of similar TPL-recruitment sequences in proto-ARFs of different charophytes clades ARFs (S5 Table) suggests that they might also interact with TPL.

Discussion

The present biochemical characterization of CaARF, a proto-ARF from an “early divergent” charophyte, identifies this protein as class C ARF, in agreement with our phylogenetic analyses (S5 Fig). Mutte et al. (2018) proposed the existence of two ARF classes in “late divergent” charophytes, C and A/B, deriving from a common ancestor that diverged in an ancient charophyte clade [20]. Based on phylogenetic analyses showing that class C ARF is sister to classes A and B, and on the identification of a M. viride sequence classified as a class A/B, Flores-Sandoval et al. (2018) proposed a similar scenario where the divergence between classes A and B and class C occurred prior to the diversification of extant streptophytes [40]. This plausible scenario, built before the identification of class C ARFs in “early divergent” charophytes, is based on an unusual M. viride sequence that does not exhibit the conserved ARF DNA binding residues (S3 Fig), and implies repeated loss of class A/B ARFs from Chlorokybophyceae to Coleochaetophyceae (S6 Fig). Further identification of class C ARFs in the “early divergent” charophytes (Klebsormidiophyceae [2] and Chlorokybophyceae (this work)) and the presence of both classes C and A/B in the “late divergent” C. orbicularis suggest a second and more parsimonious scenario in which class A/B ARF members come from an ancestral proto-ARF, belonging to class C or class C-like that existed before the emergence of “late divergent” charophytes (S6 Fig). This hypothesis implies only a few class C ARF gene losses in some Klebsormidiophyceae, Coleochaetophyceae and Zygnematophyceae species. Still, all these scenarios need to be taken with caution as they are based on transcriptomic datasets and could be challenged when genomic sequences become available.

When comparing C and A/B clades we found a disordered region within the predicted DD of ancestral and land plants clade C ARFs that is not present in clade A/B neither in land plants clades A and B. We speculate that during the duplication event leading to A/B emergence from clade C, the loss of this disordered sequence occurred. The DNA interaction experiments presented in this manuscript suggest that this event might have contributed to the acquisition of a more restricted DNA specificity of class A and B ARFs for ER motifs.

Apart from the similar behaviour observed for CaARF and AtARF10 when binding to DNA, ancestral clade C ARFs already presented PB1 oligomerization potential and interaction with the co-repressor TPL. The conservation of these properties along evolution is consistent with experiments conducted on Marchantia showing partial complementation of the loss of function MpARF3 by class C AtARF10 [40]. Moreover, these biochemical facts are instructive on several aspects of the evolution of the NAP in plants. First, proto-ARFs being able to interact with AuxREs supports that the NAP could have co-opted sets of genes already regulated by ARFs in charophytes, as suggested in other studies [2,20,39,40]. In this context, the emergence of the A/B clade with a different DNA binding behaviour could have allowed to target a more specific set of genes. Second, proto-ARF interaction with TPL provides functional evidence for a role for class C ARFs as transcriptional repressors. Putative TPL interaction motifs are also present in proto-RAV and most proto-ARFs across charophytes, which includes class A/B ARFs. The capacity to recruit TPL co-repressors could thus be an ancestral property of RAV and ARF TFs.

From these observations, we propose ARFs recruitment of co-repressor complexes to AuxREs promoter elements as a primitive and conserved mechanism predating the NAP. The absence of a functional TIR/AFB-Aux/IAA co-receptor [2,20,41] indicates that this primitive system was auxin-independent. These observations are consistent with a series of experiments in Marchantia showing that auxin-responsive genes show similar transcriptional responses in WT and MpARF3 mutants [20,40]. Alongside the diversification of ARF DNA binding specificity, emergence of the auxin perception complex in the first land plants turned ARFs-regulated genes into auxin-responsive genes through ARF-Aux/IAA-TIR/AFB interactions evolution (Fig 4).

Fig. 4. Model for ARFs and NAP evolution.
Model for ARFs and NAP evolution.
First ARF homologues that appeared in plants evolutionary line are class C ARFs, found in early charophytes and in all subsequent clades till land plants. Charophycean class C ARF ancestor recognizes different AuxRE motifs as a dimer, can oligomerize via its PB1 domain and interact with the corepressor TPL. In charophyte clades that diverged later, the class A/B ARF subfamily emerged, likely from class C. This new subfamily evolved a different DNA specificity through an alternative dimerization mode. The division of A/B into classes A and B appeared only in land plants together with TIR1-AuxIAAs perception system that confers auxin sensibility to the charophycean ARF-regulatory mechanism leading to a functional NAP.

Our work thus allowed proposing a scenario where the evolution of the binding specificity of an ancestral TF together with the emergence of a functional hormone perception complex create a hormone signalling pathway. This scenario offers a better understanding of how hormone signalling pathways can evolve from pre-existing mechanisms of transcriptional regulation independent of any hormone signalling.

Materials & methods

Protein homologues search and classification

Potential homologs of the NAP components were searched by sequence homology to the corresponding NAP proteins from M. polymorpha. Blasts were done using different databases: OneKp, PlantTFDB and Marchantia.info. Due to the lack of proteomic data in charophyte organisms, we carried out tblastn. Each transcript was then translated using Expasy Translate tool. Sequences resulting from this search were classed using protein sequence alignments and phylogenetical studies. Protein sequences alignments were done with Multialin (http://multalin.toulouse.inra.fr/multalin/) and ESPrit (http://espript.ibcp.fr/ESPript/ESPript/) online tools. Phylogenetic analyses were conducted using predicted DBDs from charophyte proto-ARFs and DBDs belonging to A. thaliana and M. polymorpha ARFs. Phylogenies were done with MEGA and Phylogeny.fr software using Maximum likelihood algorithm.

3D structure modelling

Protein structure modelling was done with Phyre2 online tool [47]. Three-dimensional structures were visualized with PyMOL software (www.pymol.org).

Plasmids construction for expression in E. coli

cDNA sequences coding for potential ancestors and the corresponding mutants were constructed as synthetic DNA (Thermofisher). KnRAV and CaARF (full-length, fragments (CaARF-DBD (residues 1–421), CaARF-PB1 (residues 644–750), KnRAV-DBD (residues 256–523), KnRAV-PB1 (residues 724–798)) or mutants) coding sequences were cloned into pETM40 plasmid (EMBL) that contains a MBP-tag in the N-terminal region except for PB1 domains from both proteins that were cloned into pETM11 (EMBL) that confers a N-terminal His-tag.

KnRAV and CaARF specific domains were isolated by PCR from synthetic cDNA sequences (S6 Table). Full-length ARF2, ARF5 and ARF10 were cloned into pHMGWA vectors (Addgene) containing N-terminal His-MBP-His tags.

Protein expression and purification

All proteins were expressed in Escherichia coli BL21 strain. Bacteria cultures were grown with the appropriate antibiotics at 37°C until they achieved an OD600nm of 0.6–0.9. Protein expression was induced with isopropyl-β-D-1-thyogalactopiranoside (IPTG) at a final concentration of 400 μM at 18°C overnight. Bacteria cultures were centrifuged, and the pellets were resuspended and sonicated in the buffers indicated in S7 Table.

After centrifugation, soluble fractions of KnRAV, KnRAV-DBD, CaARF, CaARF-DBD and CaARF mutants were loaded on Dextrin-Sepharose (GE Healthcare) column previously equilibrated in buffer A (S7 Table). After column washing, proteins were eluted in buffer A with maltose 10 mM (S6 Table).

PB1 domains of KnRAV and CaARF as well as full-length proteins ARF2, ARF5 and ARF10 were purified on Nickel-Sepharose (GE Healthcare) columns previously equilibrated in the appropriate buffers (S6 Table). After protein binding, columns were washed with 30 mM imidazole to remove all proteins non-specifically bound to the column. Proteins were eluted in the corresponding buffer containing 300 mM imidazole (S6 Table). His-tags of PB1 domains were cleaved by TEV protease (5% w/w) overnight at 4°C followed by incubation at 20°C for 2 h for SEC-MALLS experiments.

AtTPL202 and mutants were purified as explained in Martin-Arevalillo et al., 2017 [49]. Following purification step, all proteins were dialyzed for 15 h at 4°C in their purification buffers, frozen in liquid nitrogen and conserved at -80°C until used.

EMSA DNA binding tests

DNA probes were artificially designed based on the DNA binding site for each TF (S8 Table) (Eurofins). Oligonucleotides for the sense strand were designed with an overhanging G in 5’ that allows the labelling of the DNA (S8 Table). Annealing of the oligonucleotides and Cy5-labelling of the probes were performed as described in Stigliani et al.,(2019) [19]. Electrophoretic Mobility Shift Assays (EMSA), were done on native 2% agarose gels prepared with TBE buffer 0.5X. Gels were pre-run in TBE buffer 0.5X at 90 V for 90 min at 4°C. Protein-DNA mixes contained Salmon and Herring Sperm competitor DNA (final concentration 0.07 mg/ml) and labelled DNA (final concentration 20 nM) in the interaction buffer (20 mM HEPES pH 7.8; 50 mM KCl; 100 mM Tris-HCl pH 8.0, 2.5% glycerol; 1 mM DTT). Mixes were incubated in darkness for 30 min at 4°C and next loaded in the gels. Gels were run for 1 hour at 90 V at 4°C in TBE 0.5X. DNA-protein interactions were visualized with Cy5-exposition filter (Biorad ChemiDoc MP Imaging System).

Co-purification protein-protein interaction tests

For protein-protein interaction analyses, complexes between potential interaction partners were first formed by mixing MBP-tagged CaARF (wt or mutants) (90 μg) with His-tagged AtTPL202 (and mutants) (70 μg) in CAPS 20 mM pH 9.6; Tris-HCl 100 mM pH 8; NaCl 50 mM; TCEP 1 mM buffer for 1 h at 4°C. Complexes formed were fixed through the MBP tag to Dextrin-Sepharose columns previously equilibrated with CAPS 20 mM pH 9.6; Tris-HCl 100 mM pH 8; NaCl 50 mM; TCEP 0.1 mM buffer. After incubation of the complexes with Dextrin-Sepharose for 30 min at 4°C, nonspecific interactions were removed by a washing step with the same buffer. Protein complexes were eluted with 200 μl of the same buffer containing 10 mM of maltose. MBP was used as control for unspecific interactions. The eluted fractions were analysed by SDS-page polyacrylamide gel electrophoresis 12%.

SEC-MALLS Native molecular mass determination

Molecular weights were determined by Size-Exclusion Chromatography-Multi Angle Light Scattering (SEC-MALLS) on an analytical Superdex-S200 increase (GE Healthcare) connected to an in-line MALLS spectrometer (DAWN HELEOS II, Wyatt Instruments). Analytical size exclusion chromatography was performed at 25°C at a rate of 0.5 mL/min for untagged PB1 domains resulting from TEV cleavage. Untagged KnRAV-PB1 MW determination was carried out in CAPS 100 mM pH 9.6; TCEP 1mM buffer, whereas Tris-HCl 20 mM pH 8; TCEP 1 mM was used for untagged CaARF-PB1 and AtARF5. The refractive index measured with in-line refractive index detector (Optirex, Wyatt Instruments) was used to follow the differential refractive index relative to the solvent. Molecular masses calculation was done with the Debye model using ASTRA version 5.3.4.20 (Wyatt Instruments) and a theoretical dn/dc value of 0.185 mL/g.

Supporting information

S1 Fig [pptx]
Land plants evolutionary line from charophyte algae.

S2 Fig [drs]
Double AuxREs possible configurations.

S3 Fig [ddii]
ARF-DBDs alignment.

S4 Fig [pptx]
Residues conservation in KnRAV-DBD.

S5 Fig [pptx]
Phylogenetic classification of proto-ARFs in charophyte organisms.

S6 Fig [pptx]
Different evolution hypotheses.

S7 Fig [pptx]
Modelled DBDs from charophyte ancestral class C ARFs.

S8 Fig [pptx]
Charophyte TPL N-ter aligment.

S1 Table [docx]
Different B3 subfamilies and their DNA specificities.

S2 Table [docx]
RAV and ARF sequences identified in the OneKP [], Marchantia.info* [,], or genome** [] databases.

S3 Table [docx]
SEC-MALLS molecular weight determination of KnRAV-PB1, CaARF-PB1 and AtARF5-PB1 domains.

S4 Table [docx]
TPL homologues found in charophyte organisms.

S5 Table [mr]
ARF Charophytes potential EAR motifs.

S6 Table [docx]
Primers for Domains amplification.

S7 Table [docx]
Buffers used for purification.

S8 Table [docx]
DNA sequence probes used for EMSA assays.

S1 References [docx]
Supplemental references.


Zdroje

1. Finet C, Berne-Dedieu A, Scutt CP, Marletaz F. Evolution of the ARF gene family in land plants: Old domains, new tricks. Molecular Biology and Evolution. 2013;30:45–56. doi: 10.1093/molbev/mss220 22977118

2. Bowman JL, Kohchi T, Yamato KT, Jenkins J, Shu S, Ishizaki K, et al. Insights into Land Plant Evolution Garnered from the Marchantia polymorpha Genome. Cell. 2017 Oct;171(2):287–304.e15. doi: 10.1016/j.cell.2017.09.030 28985561

3. Delwiche CF, Cooper ED. The evolutionary origin of a terrestrial flora. Current Biology. 2015;25(19):R899–910. doi: 10.1016/j.cub.2015.08.029 26439353

4. Domozych DS, Popper ZA, Sørensen I. Charophytes: Evolutionary Giants and Emerging Model Organisms. Frontiers in plant science. 2016;7(October):1470. doi: 10.3389/fpls.2016.01470 27777578

5. Finet C, Timme RE, Delwiche CF, Marlétaz F. Multigene Phylogeny of the Green Lineage Reveals the Origin and Diversification of Land Plants. Current Biology. 2010 Dec;20(24):2217–22. doi: 10.1016/j.cub.2010.11.035 21145743

6. Vries J De Archibald JM. Plant evolution: landmarks on the path to terrestrial life. New Phytologist. 2018;217(4):1428–34. doi: 10.1111/nph.14975 29318635

7. Lemieux C, Otis C, Turmel M. A clade uniting the green algae Mesostigma viride and Chlorokybus atmophyticus represents the deepest branch of the Streptophyta in chloroplast genome-based phylogenies. BMC Biology [Internet]. 2007 Dec [cited 2019 Jun 3];5(1). Available from: https://bmcbiol.biomedcentral.com/articles/10.1186/1741-7007-5-2

8. Hori K, Maruyama F, Fujisawa T, Togashi T, Yamamoto N, Seo M, et al. Klebsormidium flaccidum genome reveals primary factors for plant terrestrial adaptation. Nature Communications [Internet]. 2014 Dec [cited 2018 Nov 21];5(1). Available from: http://www.nature.com/articles/ncomms4978

9. Domozych DS, Domozych CE. Multicellularity in green algae: upsizing in a walled complex. Frontiers in Plant Science. 2014;5(649):1–8.

10. Umen JG. Green Algae and the Origins of Multicellularity in the Plant Kingdom. Cold Spring Harbor Perspectives in Biology. 2014;6(11):1–27.

11. Timme RE, Bachvaroff TR, Delwiche CF. Broad Phylogenomic Sampling and the Sister Lineage of Land Plants. Joly S, editor. PLoS ONE. 2012 Jan 13;7(1):e29696. doi: 10.1371/journal.pone.0029696 22253761

12. Wickett NJ, Mirarab S, Nguyen N, Warnow T, Carpenter E, Matasci N. Phylotranscriptomic analysis of the origin and early diversification of land plants. PNAS. 2014;111(45):4859–68.

13. Lavy M, Estelle M. Mechanisms of auxin signaling. Development. 2016;143:3226–9. doi: 10.1242/dev.131870 27624827

14. Leyser O. Auxin Signaling. Plant Physiology. 2018;176(1):465–79. doi: 10.1104/pp.17.00765 28818861

15. Giraudat J, Hauge BM, Valon C, Smalle J, Parcy F, Goodman HM. Isolation of the Arabidopsis ABI3 Gene by Positional Cloning. THE PLANT CELL ONLINE. 1992 Oct 1;4(10):1251–61.

16. Yamasaki K, Kigawa T, Seki M, Shinozaki K, Yokoyama S. DNA-binding domains of plant-specific transcription factors: structure, function, and evolution. Trends in Plant Science. 2013;18(5):267–76. doi: 10.1016/j.tplants.2012.09.001 23040085

17. Galli M, Khakhar A, Lu Z, Chen Z, Sen S, Joshi T, et al. The DNA binding landscape of the maize AUXIN RESPONSE FACTOR family. Nature Communications. 2018 Oct 30;9(1):4526. doi: 10.1038/s41467-018-06977-6 30375394

18. O’Malley RC, Huang SC, Song L, Lewsey MG, Bartlett A, Nery JR, et al. Cistrome and Epicistrome Features Shape the Regulatory DNA Landscape. Cell. 2016 May;165(5):1280–92. doi: 10.1016/j.cell.2016.04.038 27203113

19. Stigliani A, Martin-Arevalillo R, Lucas J, Bessy A, Vinos-Poyo T, Mironova V, et al. Capturing Auxin Response Factors Syntax Using DNA Binding Models. Molecular Plant. 2019 Jun;12(6):822–32. doi: 10.1016/j.molp.2018.09.010 30336329

20. Mutte SK, Kato H, Rothfels C, Melkonian M, Wong GK, Weijers D. Origin and evolution of the nuclear auxin response system. eLife. 2018;March 27(7):1–25.

21. Wang C, Liu Y, Li S-S, Han G-Z. Insights into the Origin and Evolution of the Plant Hormone Signaling Machinery. Plant Physiology. 2015;167(3):872–86. doi: 10.1104/pp.114.247403 25560880

22. Guilfoyle TJ. The PB1 Domain in Auxin Response Factor and Aux/IAA Proteins: A Versatile Protein Interaction Module in the Auxin Response. The Plant Cell. 2015;27:33–43. doi: 10.1105/tpc.114.132753 25604444

23. Nanao MH, Vinos-Poyo T, Brunoud G, Thévenon E, Mazzoleni M, Mast D, et al. Structural basis for oligomerization of auxin transcriptional regulators. Nature communications. 2014;5:3617. doi: 10.1038/ncomms4617 24710426

24. Korasick DA, Westfall CS, Goo S, Nanao MH, Dumas R, Hagen G. Molecular basis for AUXIN RESPONSE FACTOR protein interaction and the control of auxin response repression. PNAS. 2014;111(14):5427–32. doi: 10.1073/pnas.1400074111 24706860

25. Parcy F, Vernoux T, Dumas R. A Glimpse beyond Structures in Auxin-Dependent Transcription. Trends in Plant Science. 2016;21:574–83. doi: 10.1016/j.tplants.2016.02.002 26994657

26. Guilfoyle TJ, Hagen G. Getting a grasp on domain III/IV responsible for Auxin Response Factor–IAA protein interactions. Plant Science. 2012 Jul;190:82–8. doi: 10.1016/j.plantsci.2012.04.003 22608522

27. Szemenyei H, Hannon M, Long J a. TOPLESS Mediates Auxin-Dependent Transcriptional Repression During Arabidopsis Embryogenesis. Science. 2008;319:1384. doi: 10.1126/science.1151461 18258861

28. Wu MF, Yamaguchi N, Xiao J, Bargmann B, Estelle M, Sang Y, et al. Auxin-regulated chromatin switch directs acquisition of flower primordium founder fate. eLife. 2015;13(4):1–20.

29. Dinesh DC, Villalobos LI a C, Abel S. Structural Biology of Nuclear Auxin Action. Trends in Plant Science. 2016;21(4):302–16. doi: 10.1016/j.tplants.2015.10.019 26651917

30. Vernoux T, Brunoud G, Farcot E, Morin V, Van den Daele H, Legrand J, et al. The auxin signalling network translates dynamic input into robust patterning at the shoot apex. Molecular Systems Biology. 2014 Apr 16;7(1):508–508.

31. Piya S, Shrestha SK, Binder B, Stewart CN, Hewezi T. Protein-protein interaction and gene co-expression maps of ARFs and Aux/IAAs in Arabidopsis. Frontiers in Plant Science [Internet]. 2014 Dec 23 [cited 2018 Mar 14];5. Available from: http://journal.frontiersin.org/article/10.3389/fpls.2014.00744/abstract

32. Trigg SA, Garza RM, MacWilliams A, Nery JR, Bartlett A, Castanon R, et al. CrY2H-seq: a massively multiplexed assay for deep-coverage interactome mapping. Nature Methods. 2017 Aug;14(8):819–25. doi: 10.1038/nmeth.4343 28650476

33. Lavy M, Prigge MJ, Tao S, Shain S, Kuo A, Kirchsteiger K, et al. Constitutive auxin response in Physcomitrella reveals complex interactions between Aux/IAA and ARF proteins. eLife. 2016;1(5).

34. Chandler JW. Auxin response factors. Plant, cell & environment. 2016;39:1014–28.

35. Causier B, Ashworth M, Guo W, Davies B. The TOPLESS Interactome: A Framework for Gene Repression in Arabidopsis. Plant Physiology. 2012;158(January):423–38. doi: 10.1104/pp.111.186999 22065421

36. Causier B, Lloyd J, Stevens L, Davies B. TOPLESS co-repressor interactions and their evolutionary conservation in plants. Plant Signaling & Behavior. 2012;7(3):325–8.

37. Flores-Sandoval E, Eklund DM, Bowman JL. A Simple Auxin Transcriptional Response System Regulates Multiple Morphogenetic Processes in the Liverwort Marchantia polymorpha. PLoS Genetics. 2015;11(5):1–26.

38. Kato H, Ishizaki K, Kouno M, Shirakawa M, Bowman JL, Nishihama R, et al. Auxin-Mediated Transcriptional System with a Minimal Set of Components Is Critical for Morphogenesis through the Life Cycle in Marchantia polymorpha. PLoS Genetics. 2015;11(5):1–26.

39. Bowman JL, Briginshaw LN, Fisher TJ, Flores-Sandoval E. Something ancient and something neofunctionalized—evolution of land plant hormone signaling pathways. Current Opinion in Plant Biology. 2019 Feb;47:64–72. doi: 10.1016/j.pbi.2018.09.009 30339930

40. Flores-Sandoval E, Eklund DM, Hong S-F, Alvarez JP, Fisher TJ, Lampugnani ER, et al. Class C ARFs evolved before the origin of land plants and antagonize differentiation and developmental transitions in Marchantia polymorpha. New Phytologist. 2018 Jun;218(4):1612–30. doi: 10.1111/nph.15090 29574879

41. Ohtaka K, Hori K, Kanno Y, Seo M, Ohta H. Primitive Auxin Response without TIR1 and Aux / IAA in the Charophyte Alga Klebsormidium nitens 1. Plant Physiology. 2017;174(July):1621–32.

42. Lemieux C, Otis C, Turmel M. Comparative Chloroplast Genome Analyses of Streptophyte Green Algae Uncover Major Structural Alterations in the Coleochaetophyceae and Zygnematophyceae. Frontiers in plant science. 2016;7(697).

43. Turmel M, Otis C, Lemieux C. Tracing the Evolution of Streptophyte Algae and Their. Genome Biology and Evolution. 2013;5(10):1817–35. doi: 10.1093/gbe/evt135 24022472

44. Matasci N, Hung L-H, Yan Z, Carpenter EJ, Wickett NJ, Mirarab S, et al. Data access for the 1,000 Plants (1KP) project. GIGA SCIENCE. 2014;3(17):1–10.

45. Swaminathan K, Peterson K, Jack T. The plant B3 superfamily. Trends in Plant Science. 2008;13(12):647–55. doi: 10.1016/j.tplants.2008.09.006 18986826

46. Boer DR, Freire-Rios A, Van Den Berg WAM, Saaki T, Manfield IW, Kepinski S, et al. Structural basis for DNA binding specificity by the auxin-dependent ARF transcription factors. Cell. 2014;156(3):577–89. doi: 10.1016/j.cell.2013.12.027 24485461

47. Kelley LA, Mezulis S, Yates CM, Wass MN, Sternberg MJE. The Phyre2 web portal for protein modelling, prediction and analysis. Nature Protocols. 2017;10(6):845–58.

48. Kagaya Y, Ohmiya K, Hattori T. RAV1, a novel DNA-binding protein, binds to bipartite recognition sequence through two distinct DNA-binding domains uniquely found in higher plants. Nucleic Acids Research. 1999;27(2):470–8. doi: 10.1093/nar/27.2.470 9862967

49. Martin-Arevalillo R, Nanao MH, Larrieu A, Vinos-poyo T, Mast D. Structure of the Arabidopsis TOPLESS corepressor provides insight into the evolution of transcriptional repression. PNAS. 2017;114(30):8107–12. doi: 10.1073/pnas.1703054114 28698367

50. Choi H, Seo M, Cho H. Two TPL-Binding Motifs of ARF2 Are Involved in Repression of Auxin Responses. Frontiers in Plant Science. 2018;9(372):1–9.

51. Ke J, Ma H, Gu X, Thelen A, Brunzelle JS, Li J. Structural basis for recognition of diverse transcriptional repressors by the TOPLESS family of corepressors. Science Advances. 2015;1(6):1–12.

52. Collins J, O’Grady K, Chen S, Gurley W. The C-terminal WD40 repeats on the TOPLESS co-repressor function as a protein–protein interaction surface. Plant Molecular Biology. 2019 May;100(1–2):47–58. doi: 10.1007/s11103-019-00842-w 30783952

Štítky
Genetika Reprodukční medicína

Článek vyšel v časopise

PLOS Genetics


2019 Číslo 9
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

Důležitost adherence při depresivním onemocnění
nový kurz
Autoři: MUDr. Eliška Bartečková, Ph.D.

Koncepce osteologické péče pro gynekology a praktické lékaře
Autoři: MUDr. František Šenk

Sekvenční léčba schizofrenie
Autoři: MUDr. Jana Hořínková, Ph.D.

Hypertenze a hypercholesterolémie – synergický efekt léčby
Autoři: prof. MUDr. Hana Rosolová, DrSc.

Multidisciplinární zkušenosti u pacientů s diabetem
Autoři: Prof. MUDr. Martin Haluzík, DrSc., prof. MUDr. Vojtěch Melenovský, CSc., prof. MUDr. Vladimír Tesař, DrSc.

Všechny kurzy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#