# Evolution of the Auxin Response Factors from charophyte ancestors

Raquel Martin-Arevalillo, Emmanuel Thevenon, Fanny Jégu, Thomas

Vinos-Poyo, Teva Vernoux, François Parcy, Renaud Dumas

► To cite this version:

Raquel Martin-Arevalillo, Emmanuel Thevenon, Fanny Jégu, Thomas Vinos-Poyo, Teva Vernoux, et al.. Evolution of the Auxin Response Factors from charophyte ancestors. PLoS Genetics, 2019, 15 (9), pp.e1008400. 10.1371/journal.pgen.1008400 . hal-02329218

HAL Id: hal-02329218

https://hal.science/hal-02329218

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Evolution of the Auxin Response Factors from charophyte ancestors

Raquel Martin-Arevalillo¹, Emmanuel Thévenon², Fanny Jégu ID², Thomas Vinos-Poyo², Teva Vernoux ID¹, François Parcy², Renaud Dumas ID²*

1 Laboratoire de Reproduction et Développement des Plantes, Univ. Lyon, ENS de Lyon, UCB Lyon1, CNRS, INRA, Lyon, France, 2 Univ. Grenoble Alpes, CNRS, CEA, INRA, IRIG-DBSCI-LPCV, Grenoble, France

* renaud.dumas@cea.fr

<figure><img src="image_2.png" /></figure>

## Abstract

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.

<figure><img src="image_3.png" /></figure>

### OPEN ACCESS

**Citation:** Martin-Arevalillo R, Thévenon E, Jégu F, Vinos-Poyo T, Vernoux T, Parcy F, et al. (2019) Evolution of the Auxin Response Factors from charophyte ancestors. PLoS Genet 15(9): e1008400. https://doi.org/10.1371/journal.pgen.1008400

**Editor:** Jason Reed, University of North Carolina, UNITED STATES

**Received:** July 17, 2019

**Accepted:** September 4, 2019

**Published:** September 25, 2019

**Copyright:** © 2019 Martin-Arevalillo et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the manuscript and its Supporting Information files.

**Funding:** This work was supported by the Agence Nationale de la Recherche [ANR-12-BSV6-0005 Auxiflo to R.D., T.V., and F.P.][ANR-18-CE12-0014 ChromAuxi to R.M.A., R.D., T.V., and F.P.], a PhD fellowship from the University Grenoble Alpes [R. M.A.] and the Grenoble Alliance for Cell and Structural Biology [ANR-10-LABX-49-01 to F.P. and R.D.]. The funders had no role in study design,

## Author summary

Plants transition from water to land was determining for the history of our planet, since it led to atmospheric and soil condition changes that promoted the appearance of other life forms. This transition initiated around 1 billion years ago from a Charophyte algae lineage that acquired features allowing it to adapt to the very different terrestrial conditions. Land plants coordinate their development with external stimuli through signalling mechanisms triggered by plant hormones. Therefore, evolution of these molecules and their signalling pathways likely played an important role in the aquatic to terrestrial move. In this manuscript we study the origin of auxin signalling, a plant hormone implicated in all plant developmental steps. Our studies suggest that out of the three families of proteins originally proposed to trigger auxin signalling in land plants, only one existed in Charophyte
data collection and analysis, decision to publish, or
preparation of the manuscript.

**Competing interests:** The authors have declared
that no competing interests exist.

ancestors as a likely transcriptional repressor independent of auxin. We show that despite millions of years of evolution, this family of proteins has conserved its biochemical and structural properties that are found today in land plants. The results presented here provide an insight on how hormone signalling pathways could have evolved by co-opting a pre-existing hormone-independent transcriptional regulatory mechanism.

**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. [1–5]. “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) [6,7]. “Late divergent” clades include Charophyceae and Coleochaetophyceae that share features with land plants (S1 Fig), [8,9] but also Zygnematophyceae, that despite their simple structure are considered sisters to land plants according to recent phylogenetic studies [10,11].

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 [1]. 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) [12,13]. 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 [14]. 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) [1,17–19].

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 [1,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 [23–26] and recruit TOPLESS (TPL)/TOPLESS-RE-LATED 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 [30–32] 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 [34–36].

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
<figure><img src="image_4.png" /><figcaption>**Fig 1.** 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 [9,42]. Although some publications placed Mesostigmaphoticae and Chlorokybophycae in a unified clade large differences in sequences and morphological traits argues for two different clades [2,6,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.</figcaption></figure>

https://doi.org/10.1371/journal.pgen.1008400.g001

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 [1,20,37,39–41]. How these ancestral components evolved to form the land plants NAP remains an open question. Through the structural, biochemical and phylogenetic characterization 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) [11,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 [1]. 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 *Klebsormidium 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.

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
<figure><img src="image_5.png" /><figcaption>**Fig 2.** 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 7 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).</figcaption></figure>

https://doi.org/10.1371/journal.pgen.1008400.g002

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.
<figure><img src="image_6.png" /><figcaption>**Fig 3.** 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.</figcaption></figure>

https://doi.org/10.1371/journal.pgen.1008400.g003

## 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<sup>-F74A</sup>, mutated in the hydrophobic EAR peptide binding groove (Fig 3I) [49]. Moreover, mutations in CaARF KLFG (CaARF<sup>-L523S/F524S</sup>) or LPLLPS (CaARF<sup>-L492A/L493A</sup>) 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” char-
ophytes, 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 sce-
nario, 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 resi-
dues (S3 Fig), and implies repeated loss of class A/B ARFs from Chlorokybophyceae to Coleo-
chaetophyceae (S6 Fig). Further identification of class C ARFs in the “early divergent”
charophytes (Klebsormidiophyceae [1] 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 par-
simonious 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” charo-
phytes (S6 Fig). This hypothesis implies only a few class C ARF gene losses in some Klebsormi-
diophyceae, 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 [1,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 [1,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).

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 Multalin (http://multalin.toulouse.inra.fr/multalin/) and ESPrit (http://esprit.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
<figure><img src="image_7.png" /><figcaption>**Fig 4.** 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.</figcaption></figure>

https://doi.org/10.1371/journal.pgen.1008400.g004

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 OD₆₀₀nm 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 1 mM 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. Land plants evolutionary line from charophyte algae.** Acquisition of more complex structures and features similar to those found in land plants is observed along charophytes evolutionary line (Adapted from [9]).
(PPTX)

**S2 Fig. Double AuxREs possible configurations.** ARFs binding sites are double sites which can be Direct Repeats (DRs), with the binding sites located in the same DNA brand, Everted Repeats (ER) or Inverted Repeats (IR), with the binding sites in different DNA brands.
(PPTX)

**S3 Fig. ARF-DBDs alignment.** Predicted charophyte ARF-DBDs aligned to DBDs of classes A (A.tha-5, 6, 7, 8, 19 and M.pol-1), B (A.tha-1, 2–4, 9, 11–15, 18, 20–22 and M.pol-2) and C (A. tha-10, 16 and 17 and M.pol-3) ARFs from *A. thaliana* and *M. polymorpha*. Class-C ARFs present an insertion in the DBD that is located inside the second part of the dimerization domain (DDII) (underlined) described for the ARFs (in between residues 260 and 280 in the alignment, referenced to A.tha-5). In agreement, ancestral ARFs sharing this insertion (C.atm, Entr, N.mir, C.irr, C.scu, M.end, S.pra) were classed with class-C ARFs, whereas proto-ARFs lacking this insertion (C. orb and Moug) belonged to class A/B (See phylogeny in Supplemental S4 Fig).
Note that the M. vir sequence [1] does not contain the consensus B3 DNA binding sequence. Abbreviations used in the alignment: A.tha, *A. thaliana*; M.pol, *M. polymorpha*; M.vir, *M. viride*; C.atm, *C.atmophyticus*; Entr, *Entransia*; N.mir, *N. mirabilis*; C.irr, *C. irregularis*; C.scu, *C. scutata*; C.orb, *C. orbicularis*; Moug, *Mougeotia*; M.end, *M. endlicheranium*; S.pra, *S. pratensis*. Discontinuous underlines mark the regions involved in dimerization (DDI and DDII); violet underline marks B3 domain; black stars indicate the residues implicated in the interaction with AuxREs. Arrows point at C.atm and M.vir sequences. The incomplete sequence of the DBD of the class C ARF from *C. orbicularis* (GBSL01007362) was not added in the alignment.
(PPTX)

**S4 Fig. Residues conservation in KnRAV-DBD.** A, B3RAV domains alignment. Black stars indicate residues implicated in the interaction with DNA. WNSSQS, amino acids characteristic of B3RAV TFs [45], are conserved in KnRAV (residues numbering referred to *A.thaliana* RAV1, A.tha-RAV1). B, AP2 domains alignment. Black stars indicate residues implicated in the interaction with DNA [53] (residues numbering referred to *A. thaliana* ERF1). Abbreviations used in the alignment: A.tha, *A. thaliana*; M.pol, *M. polymorpha*; K.nit, *K. nitens*.
(PPTX)

**S5 Fig. Phylogenetic classification of proto-ARFs in charophyte organisms.** Maximum likelihood tree (built with MEGA software from DBD sequences) showing ARF evolutionary clades (A, B, C and A/B). Bootstrap values are shown next to branches. Abbreviations: A.tha, *A. thaliana*; M.pol, *M. polymorpha*; M.vir, *M. viride*; C.atm, *C. atmophyticus*; Entr, *Entransia*; N.mir, *N. mirabilis*; C.orb, *C. orbicularis*; C.irr, *C. irregularis*; C.scu, *C. scutata*; Moug, *Mougeotia*; M.end, *M. endlicheranium*; S.pra, *S. pratensis*. The incomplete sequence of the DBD of the class C ARF from *C. orbicularis* (GBSL01007362) was not added in the analysis.
(PPTX)

**S6 Fig. Different evolution hypotheses.** A. ARF C and A/B originated from a common ancestor that had already diverged in an early charophyte and evolved independently in later clades,
with subsequent losses in different clades/species. **B.** The presence of ARF C homologues from the first clades of charophytes evolutionary line suggests this subfamily or a closely-related one (C-like), as the common ancestor for current charophytean A/B and C ARFs. In both scenarios duplication of A/B into A and B happened in land plants. **C.** Phylogenetic tree generated by Maximum likelihood (phylogeny.fr [54,55]) that supports charophytes C clade as ancestor of charophyte and land plants ARF subfamilies.
(PPTX)

**S7 Fig.** *Modelled DBDs from charophyte ancestral class C ARFs.* Modelled structures superposed to ARF1 DBD structure (4LDX, in grey [46]). **A.** *Entransia* ARF-DBD, model in red. **B.** *N. mirabilis* ARF-C-DBD, model in green. **C.** *C. scutata* ARF-DBD, model in blue. **D.** *M. endlicheranium* ARF-DBD, model in yellow. **E.** *S. pratensis* ARF-DBD, model in purple. Yellow dots indicate the site of the insertion characteristic of ARF-C class members positioned either inside or at the end of the helix belonging to the DDII, depending on the model.
(PPTX)

**S8 Fig.** *Charophyte TPL N-ter alignment.* TPL homologues found in charophytes aligned to *A. thaliana* (A.tha) and *M. polymorpha* (M.pol) TPL N-ter. Residues involved in the interaction with EAR motifs indicated with a triangle. Residues involved in TPL tetramerization indicated with a star. In the alignment: C.atm, *C.atmophyticus*; Entr, *Entransia*; N.mir, *N. mirabilis*; C.irr, *C. irregularis*; C.scu, *C. scutata*; C.orb, *C. orbicularis*; Moug, *Mougeotia*; M.end, *M. endlicheranium*; S.pra, *S. pratensis*.
(PPTX)

**S1 Table.** *Different B3 subfamilies and their DNA specificities.* B3 domains were classed into B3ABI3 (not reported in this manuscript), B3RAV and B3ARF subfamilies according to the residues present in the predicted DNA-interacting loop. Indicated in the table the characteristic amino acidic sequence of the DNA-interacting loop and the DNA binding sequence for each B3 subfamily [45,56].
(DOCX)

**S2 Table.** *RAV and ARF sequences identified in the OneKP [11], Marchantia.info*<sup>+</sup> [1,40], or *Klebsormidium nitens* genome<sup>**</sup> [7] databases.
(DOCX)

**S3 Table.** *SEC-MALLS molecular weight determination of KnRAV-PB1, CaARF-PB1 and AtARF5-PB1 domains.*
(DOCX)

**S4 Table.** *TPL homologues found in charophyte organisms.* Accession numbers for transcripts or proteins and the databases used for each search are indicated. Amino acidic sequences were obtained by transcripts translation, except for *K.nitens*-RAV protein, obtained from PlantTFDB. Predicted domains are indicated with a tick.
(DOCX)

**S5 Table.** *ARF Charophytes potential EAR motifs.* Accession numbers for transcripts or proteins and the databases used for each search are indicated. Potential EAR motifs in the Middle Regions (MR) were searched for each protein, with the MR corresponding to the sequence in between the DBD domains and the PB1 domains. Possible EAR motifs were identified as potential TPL-recruitment sites based on the EAR/EAR-like motifs described in TPL interac-tome publication [35,36].
(DOCX)
**S6 Table. Primers for Domains amplification.** Note that CaARF DBD and PB1 domains were isolated from CaARF through restriction sites introduced in CaARF synthetic cDNA (Thermofisher).
(DOCX)

**S7 Table. Buffers used for purification.**
(DOCX)

**S8 Table. DNA sequence probes used for EMSA assays.** Sequences in bold indicating the binding sites and in italic the mutated sites.
(DOCX)

**Acknowledgments**

We thank Mathilde Armand and Elena Suber for their technical assistance, Roger Miras (Laboratoire Chimie et Biologie des Métaux Grenoble) for help with SEC-MALLS experiments and David Cobessi and Charlie Scutt for helpful discussion.

**Author Contributions**

**Conceptualization:** Raquel Martin-Arevalillo, Renaud Dumas.

**Data curation:** Raquel Martin-Arevalillo, Emmanuel Thévenon, Fanny Jégu, Thomas Vinos-Poyo, Renaud Dumas.

**Formal analysis:** Raquel Martin-Arevalillo, Emmanuel Thévenon, Fanny Jégu, Thomas Vinos-Poyo, Renaud Dumas.

**Funding acquisition:** Teva Vernoux, François Parcy, Renaud Dumas.

**Investigation:** Raquel Martin-Arevalillo, Emmanuel Thévenon, Thomas Vinos-Poyo, Renaud Dumas.

**Methodology:** Raquel Martin-Arevalillo, Emmanuel Thévenon, Fanny Jégu, Thomas Vinos-Poyo, Renaud Dumas.

**Project administration:** Renaud Dumas.

**Resources:** Renaud Dumas.

**Supervision:** Renaud Dumas.

**Validation:** Raquel Martin-Arevalillo, Emmanuel Thévenon, Thomas Vinos-Poyo, Teva Vernoux, François Parcy, Renaud Dumas.

**Visualization:** Raquel Martin-Arevalillo, Renaud Dumas.

**Writing – original draft:** Raquel Martin-Arevalillo, Renaud Dumas.

**Writing – review & editing:** Raquel Martin-Arevalillo, Teva Vernoux, François Parcy, Renaud Dumas.

**References**

1. 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. https://doi.org/10.1016/j.cell.2017.09.030 PMID: 28985561

2. Delwiche CF, Cooper ED. The evolutionary origin of a terrestrial flora. Current Biology. 2015; 25(19): R899–910. https://doi.org/10.1016/j.cub.2015.08.029 PMID: 26439353
3. Domozych DS, Popper ZA, Sørensen I. Charophytes: Evolutionary Giants and Emerging Model Organisms. Frontiers in plant science. 2016; 7(October):1470. https://doi.org/10.3389/fpls.2016.01470 PMID: 27777578

4. 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. https://doi.org/10.1016/j.cub.2010.11.035 PMID: 21145743

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

6. 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

7. 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

8. Domozych DS, Domozych CE. Multicellularity in green algae: upsizing in a walled complex. Frontiers in Plant Science. 2014; 5(649):1-8. https://doi.org/10.3389/fpls.2014.00649 PMID: 25477895

9. Umen JG. Green Algae and the Origins of Multicellularity in the Plant Kingdom. Cold Spring Harbor Perspectives in Biology. 2014; 6(11):1-27. https://doi.org/10.1101/cshperspect.a016170 PMID: 25324214

10. 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. https://doi.org/10.1371/journal.pone.0029696 PMID: 22253761

11. 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. https://doi.org/10.1073/pnas.1323926111 PMID: 25355905

12. Lavy M, Estelle M. Mechanisms of auxin signaling. Development. 2016; 143:3226-9. https://doi.org/10.1242/dev.131870 PMID: 27624827

13. Leyser O. Auxin Signaling. Plant Physiology. 2018; 176(1):465-79. https://doi.org/10.1104/pp.17.00765 PMID: 28818861

14. 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. https://doi.org/10.1093/molbev/mss220 PMID: 22977118

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. https://doi.org/10.1105/tpc.4.10.1251 PMID: 1359917

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. https://doi.org/10.1016/j.tplants.2012.09.001 PMID: 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. https://doi.org/10.1038/s41467-018-06977-6 PMID: 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. https://doi.org/10.1016/j.cell.2016.04.038 PMID: 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. https://doi.org/10.1016/j.molp.2018.09.010 PMID: 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. https://doi.org/10.7554/eLife.33399.001 PMID: 29580381

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. https://doi.org/10.1104/pp.114.247403 PMID: 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. https://doi.org/10.1105/tpc.114.132753 PMID: 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. https://doi.org/10.1038/ncomms4617 PMID: 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. https://doi.org/10.1073/pnas.1400074111 PMID: 24706860

25. Parcy F, Vernoux T, Dumas R. A Glimpse beyond Structures in Auxin-Dependent Transcription. Trends in Plant Science. 2016; 21:574-83. https://doi.org/10.1016/j.tplants.2016.02.002 PMID: 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. https://doi.org/10.1016/j.plantsci.2012.04.003 PMID: 22608522

27. Szemenyei H, Hannon M, Long J A. TOPLESS Mediates Auxin-Dependent Transcriptional Repression During Arabidopsis Embryogenesis. Science. 2008; 319:1384. https://doi.org/10.1126/science. 1151461 PMID: 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. https://doi.org/10.7554/ eLife.09269.001 PMID: 26460543

29. Dinesh DC, Villalobos LI a C, Abel S. Structural Biology of Nuclear Auxin Action. Trends in Plant Science. 2016; 21(4):302-16. https://doi.org/10.1016/j.tplants.2015.10.019 PMID: 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. https://doi.org/10.1038/msb.2011.39 PMID: 21734647

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. https://doi.org/10.1038/nmeth.4343 PMID: 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). https://doi.org/10.7554/eLife.13325.001 PMID: 27247276

34. Chandler JW. Auxin response factors. Plant, cell & environment. 2016; 39:1014-28. https://doi.org/10.1111/pce.12662 PMID: 26487015

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. https://doi.org/10.1104/pp.111. 186999 PMID: 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. https://doi.org/10.4161/psb. 19283 PMID: 22476455

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. https://doi.org/10.1371/journal.pgen.1005207 PMID: 26020649

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. https://doi.org/10.1371/journal.pgen. 1005365 PMID: 26020919

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. https://doi.org/10.1016/j.pbi.2018.09.009 PMID: 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. https://doi.org/10.1111/nph. 15090 PMID: 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. https://doi.org/ 10.1104/pp.17.00274 PMID: 28533212

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). https://doi.org/10.3389/fpls.2016.00697 PMID: 27252715

43. Turmel M, Otis C, Lemieux C. Tracing the Evolution of Streptophyte Algae and Their. Genome Biology and Evolution. 2013; 5(10):1817-35. https://doi.org/10.1093/gbe/evt135 PMID: 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. https://doi.org/10.1186/2047-217X-3-17 PMID: 25625010

45. Swaminathan K, Peterson K, Jack T. The plant B3 superfamily. Trends in Plant Science. 2008; 13 (12):647–55. https://doi.org/10.1016/j.tplants.2008.09.006 PMID: 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. https://doi.org/10.1016/j.cell.2013.12.027 PMID: 24485461

47. Kelley LA, Mezulis S, Yates CM, Wass MN, Sternberg MJE. The Phyre2 web portal for protein modeling, prediction and analysis. Nature Protocols. 2017; 10(6):845–58. https://doi.org/10.1038/nprot.2015.053 PMID: 25950237

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. https://doi.org/10.1093/nar/27.2.470 PMID: 9862967

49. Martin-Arevalillo R, Nanao MH, Larrieu A, Vinos-poyo T, Mast D, et al. Structure of the Arabidopsis TOPLESS corepressor provides insight into the evolution of transcriptional repression. PNAS. 2017; 114(30):8107–12. https://doi.org/10.1073/pnas.1703054114 PMID: 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. https://doi.org/10.3389/fpls.2018.00372 PMID: 29619039

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. https://doi.org/10.1126/sciadv.1500107 PMID: 26601214

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. https://doi.org/10.1007/s11103-019-00842-w PMID: 30783952

53. Allen MD, Yamasaki K, Ohme-Takagi M, Tateno M, Suzuki M. A novel mode of DNA recognition by a β-sheet revealed by the solution structure of the GCC-box binding domain in complex with DNA. The EMBO Journal. 1998; 17(18):5484–96. https://doi.org/10.1093/emboj/17.18.5484 PMID: 9736626

54. Dereeper A, Guignon V, Blanc G, Audic S, Buffet S, Chevenet F, et al. Phylogeny. fr: robust phylogenetic analysis for the. Nucleic Acids Research. 2008; 36(April):465–9. https://doi.org/10.1093/nar/gkn180 PMID: 18424797

55. Dereeper A, Audic S, Claverie J, Blanc G. BLAST-EXPLORER helps you building datasets for phylogenetic analysis. BMC Evolutionary Biology. 2010; 10(8):8–13. https://doi.org/10.1186/1471-2148-10-8 PMID: 20067610

56. Golovenko D, Manakova E, Zakrys L, Zaremba M, Sasnauskas G, Gražulis S, et al. Structural insight into the specificity of the B3 DNA-binding domains provided by the co-crystal structure of the C-terminal fragment of Bfil restriction enzyme. Nucleic Acids Research. 2014; 42:4113–22. https://doi.org/10.1093/nar/gkt1368 PMID: 24423868