Origin and evolution of the nuclear auxin response system

  1. Sumanth K Mutte
  2. Hirotaka Kato
  3. Carl Rothfels
  4. Michael Melkonian
  5. Gane Ka-Shu Wong
  6. Dolf Weijers  Is a corresponding author
  1. Wageningen University, Netherlands
  2. University of California, United States
  3. University of Cologne, Germany
  4. University of Alberta, Canada
  5. BGI-Shenzhen, China
8 figures, 1 table and 6 additional files

Figures

Figure 1 with 3 supplements
Proteins in nuclear auxin pathway; mechanism and origin of the domains.

(A, B) Scheme of NAP in land plants. In the absence of auxin, Aux/IAA inhibit ARF via their PB1 domains, and by recruiting the TPL co-repressor. Auxin stabilizes the interaction between Aux/IAA and TIR1/AFB, followed by proteasome-mediated degradation of Aux/IAA. (C) Domain structure of NAP components in land plants and presence of each domain in algae, as recovered in transcriptomes. (Sub-)domains are indicated by colors, that match those in (B). ^: basal to all B3-type transcription factors in land plants, *: difficult to assign to ARF or Aux/IAA family; #: forming basal clade to both TIR1/AFB and COI1 in land plants.

https://doi.org/10.7554/eLife.33399.003
Figure 1—figure supplement 1
The work flow of phylogenetic tree construction.
https://doi.org/10.7554/eLife.33399.004
Figure 1—figure supplement 2
Phylogenetic tree of ARF and RAV proteins.

Label color shows the taxonomic group of each protein as indicated in the box above. Numbers along with the branches indicate branch length. Orange circles indicate the bootstraps higher than 75. Colored boxes connected with gray bar shows the domain structure of each protein. Red: B3, green: DD2 + AD, blue: PB1, gray: AP2. The complete tree can be found at http://itol.embl.de/shared/dolfweijers (interactive Tree of Life; iTOL).

https://doi.org/10.7554/eLife.33399.005
Figure 1—figure supplement 3
DNA-binding domain of RAV proteins.

(A) Domain structure of RAV proteins in land plants. (B) Homology models for B3 domain of A. thaliana ARF1 (gray) and K. nitens proto-RAV are merged. Four serine residues which are critical for DNA binding of ARF is indicated as stick model. (C) Multiple alignment for B3 domain of ARFs and (proto-)RAVs. Numbering is based on the ARF1 protein of A. thaliana. Red and orange triangles indicate the residues which are important for DNA interaction in ARF proteins.

https://doi.org/10.7554/eLife.33399.006
Figure 2 with 1 supplement
Distribution of auxin signaling proteins precursors in algal lineages.

(A) Occurrence of NAP components in red algae, chlorophytes, and charophytes. Empty circles and filled circles indicate the absence and presence of that particular component, respectively. *: Land plants have defined three classes of ARFs, RAV without PB1, and separate TIR1/AFB and COI1 receptors. (B) Schematic illustration of the phylogenetic arrangement of RAV1, Aux/IAA and ARFs based on the DBD tree and PB1 tree. Note that only branches with strong bootstrap support are shown.

https://doi.org/10.7554/eLife.33399.007
Figure 2—figure supplement 1
Phylogenetic tree based on PB1 domain.

Colored branches indicate protein families. Orange: Chlorophytes, green: proto-RAV, blue: Aux/IAA, black: (proto-)ARF. Numbers along with the branches indicate bootstrap values. The complete tree can be found at http://itol.embl.de/shared/dolfweijers (interactive Tree of Life; iTOL).

https://doi.org/10.7554/eLife.33399.008
Homology models of ancestral ARF, Aux/IAA and TIR1/AFB proteins.

(A) Homology models for ARF DBDs. The crystal structure of Arabidopsis thaliana ARF1-DBD is shown on the left with important residues for DNA binding (top) and dimerization (bottom). Homology models for (proto-)ARFs are overlaid on A. thaliana ARF1 in right panels (brown). (B) Alignment of PB1 domain of (proto-)ARF, Aux/IAA and proto-RAV proteins. Numbering is based on the ARF5 protein of A. thaliana. Arrows and helices indicate β-sheets and α-helices in ARF5 and IAA17 of A.thaliana, respectively. Blue and red triangles indicate positive (+) and negative (-) faces, respectively. Golden asterisks represent the residues of polar interactions. (C) Homology models for TIR1/AFB and COI1 proteins. Left panel shows crystal structure of A. thaliana TIR1 from top view. Auxin-binding pocket of TIR1/AFB and jasmonate-binding pocket of COI1 are shown in right panels. Hormone-binding residues are indicated as stick model in TIR1 and COI1 of land plants. Blue, red or orange residues in the model for the Coleochaete irregularis protein indicate the residues aligned with hormone binding residues of TIR1, COI1 or both, respectively. Ci: Coloechaete irregularis, Da: Desmidium aptogonum, Kn: Klebsormidium nitens, Mc: Mesotaenium caldariorum, Mp: Marchantia polymorpha, Pk: Parachlorella kessleri, Ps: Pisum sativum, Sp: Spirogyra pratensis.

https://doi.org/10.7554/eLife.33399.009
Figure 4 with 3 supplements
Reconstruction of ancestral state of NAP components in plant evolution.

Phylogeny of taxonomic classes are shown in left. Time point of the lineage diversification was calculated using TimeTree database (Kumar et al., 2017). Black stars indicate whole genome duplication events (Jiao et al., 2011). Right: phylogenetic trees show the copy number and phylogenetic relationship of each protein family in the common ancestors. Each circle is colored according to protein type as indicated in the box. In the top row, numbers indicate which genes of Arabidopsis thaliana belong to each subfamily and red circles indicates missing subfamilies in A. thaliana. Note that only branches with strong bootstrap support are shown.

https://doi.org/10.7554/eLife.33399.010
Figure 4—figure supplement 1
Phylogenetic tree of Aux/IAA.

Label color shows the taxonomic group of each protein as indicated in top. Colored boxes connected with gray bar shows the domain structure of each protein. Magenta: domain I, yellow domain II, blue: PB1. Numbers along with the branches indicate branch length. Orange circles indicate bootstrap values higher than 75. The complete tree can be found at http://itol.embl.de/shared/dolfweijers (interactive Tree of Life; iTOL).

https://doi.org/10.7554/eLife.33399.011
Figure 4—figure supplement 2
Phylogenetic tree of the proteins containing F-box and LRR.

Colored branches indicate protein families. Green: TIR1/COI1 precursor of Charophytes, red: COI1, orange: TIR1/AFB, black: the others. Numbers along with the branches indicate bootstrap values. The complete tree can be found at http://itol.embl.de/shared/dolfweijers (interactive Tree of Life; iTOL).

https://doi.org/10.7554/eLife.33399.012
Figure 4—figure supplement 3
Phylogenetic tree of TIR1/AFB.

Label color shows the taxonomic group of each protein as indicated in left. Numbers along with the branches indicate branch length. Orange circles indicate bootstrap values higher than 75.

https://doi.org/10.7554/eLife.33399.013
Auxin-dependent gene regulation across basal plant species.

(A) Histograms represent the distribution of log2 fold change among differentially expressed genes on X-axis (Padj <0.01). Y-axis indicates the number of genes in each log2 fold-change bin. Pie charts indicate the total number of up- and down-regulated genes in each species. Circles in the top left of each graph indicate the number of NAP components. (B) Violin plots of log2 normalized expression values (DEseq2-based; y-axis) of 20 least auxin activated (Low) and 20 top-most auxin upregulated (High) genes in each six species. White dot indicated the median expression value.

https://doi.org/10.7554/eLife.33399.014
Figure 5—source data 1

Number of DEG in de novo assembly- or genome-based transcriptome analysis.

https://doi.org/10.7554/eLife.33399.015
Figure 5—source data 2

Summary statistics of comparative RNA-seq analysis.

https://doi.org/10.7554/eLife.33399.016
Figure 6 with 2 supplements
Identification of deeply conserved auxin-responsive genes.

(A) Auxin-dependence of six well-known angiosperm auxin-responsive gene families (top) surveyed from de novo assembly-based transcriptomes in six species. Each circle indicates a gene copy of each gene family. Red, blue and grey circle indicate up-, down- and non-regulated genes in response to auxin. #: no homologues were identified in our transcriptome possibly due to low expression, or they might be lost during evolution. (B) qPCR analysis of conserved auxin-responsive genes. Auxin treatment was performed in the same condition with RNA-seq experiment (10 μM 2,4-D for 1 hr). Relative expression values are normalized by the expression of EF1α in Marchantia polymorpha or the amount of total RNA in Anthoceros agrestis and Ceratopteris richardii. Each bar indicates average of expression with SD (biological replicates ≥3). *: p<0.01 (t-test).

https://doi.org/10.7554/eLife.33399.017
Figure 6—figure supplement 1
Network of up-regulated genes shared between different species upon auxin treatment.

Nodes represent the genes and edges represent the presence of BLAST similarity. Colors indicate the species in the legend above. Note that two edges connect nodes if the genes are bi-directional BLAST hits. See also Supplementary file 4.

https://doi.org/10.7554/eLife.33399.018
Figure 6—figure supplement 2
Network of down-regulated genes shared between different species upon auxin treatment.

Nodes represent the genes and edges represent the presence of BLAST similarity. Note that two edges connect nodes if the genes are bi-directional BLAST hits. Colors indicate the species in the legend above. See also Supplementary file 5.

https://doi.org/10.7554/eLife.33399.019
Figure 7 with 1 supplement
Genetic analysis of ancient components in Marchantia polymorpha.

(A) Diagrams of gene structure and CRISPR/Cas9-mediated mutation in ncIAA, ncARF and ARF3 loci. Arrowheads indicate sgRNAs target sites. Gray and black boxes indicate UTR and CDS, respectively. Red and blue bars indicate the region coding PB1 and DBD. (B) 10-day-old gemmalings grown without or with 3 μM 2,4-D. Scale bars: 5 mm. (C) Expression analysis of auxin-responsive genes in WT, nciaa, ncarf, and arf3 mutants by qPCR. 10-day-old gemmalings (nciaa and ncarf) or regenerating thalli (arf1 and arf3) were treated with 10 μM 2,4-D for 1 hr. Each bar indicates average ±SD (biological replicates = 3). Asterisks indicate significant differences. *: p<0.01 (Tukey test), **: p<0.05 (t-test). (D, E) Thallus tips grown for 2 weeks (D) and gemma cups (E) of WT and arf3ge1-1 mutant. arf3ge1-1 showed growth retardation and no mature gemmae, similar to the other alleles. (F) Expression analysis of proARF3:ARF3-Citrine in arf3ge2-1 background. Left and right panel show developing and mature gemmae, respectively. Scale bars: 5 mm in (B and D), 0.5 mm in (E), 50 μm in (F).

https://doi.org/10.7554/eLife.33399.020
Figure 7—figure supplement 1
CRISPR/Cas9-mediated mutagenesis in M.

polymorpha. (A–D) Mutations detected by sequencing analysis. The amino acid (AA) sequences encoded in WT are shown at the top. WT sequence is shown with the PAM sequence highlighted in red and the target sequence of sgRNA in blue. Purple bases indicate mutation. nciaa-6: 6 bp deletion and 20 bp insertion, nciaa-10: 776 bp deletion and 75 bp insertion, ncarf-2: 1 bp insertion, ncarf-10: 486 bp deletion and 6 bp insertion, arf3ge1-1: 5 bp deletion, arf3ge1-2: 11 bp deletion and 72 bp insertion, arf3ge2-1: 9 bp deletion. (E) Three-week-old gemmalings. Arrows indicate the thalli formed with up-side-down. (F) qPCR analysis on 10-day-old gemmalings with or without 10 μM 2,4-D treatment for 1 hr. Relative expression values are normalized by the expression of EF1α. Each bar indicates average with SD (biological replicate = 3). Each asterisk indicates significant difference between WT and mutants in the same condition (p<0.01, Tukey test). (G) Thallus tips of WT and arf3 mutants grown for 2 weeks. Scale bars = 5 mm.

https://doi.org/10.7554/eLife.33399.021

Tables

Key resources table
Reagent type (species)
or resource
DesignationSource or referenceIdentifiersAdditional information
Strain, strain background (Anthoceros agresitis)OxfordPMID: 26146510
Strain, strain background (Ceratopteris richardii)HnnPMID: 25886741
Strain, strain background (Klebsormidium nitens)NIES-2285National Institute of Environmental Studies (Japan)
Strain, strain background (Physcomitrella patens)GrandsdenPMID: 18079367
Strain, strain background (Spirogyra pratensis)UTEX928The University of Texas at Austin
Strain, strain background (Marchantia polymorpha)Tak-1PMID: 26020919
Genetic reagent (M. polymorpha)ncarf-2this papermutant of ncARF of M. polymorpha, Tak-1 background
Genetic reagent (M. polymorpha)ncarf-10this papermutant of ncARF of M. polymorpha, Tak-1 background
Genetic reagent (M. polymorpha)nciaa-6this papermutant of ncIAA ofM. polymorpha, Tak-1 background
Genetic reagent (M. polymorpha)nciaa-10this papermutant of ncIAA ofM. polymorpha, Tak-1 background
Genetic reagent (M. polymorpha)arf3ge1-1this papermutant of ARF3 ofM. polymorpha, Tak-1 background
Genetic reagent (M. polymorpha)arf3ge1-2this papermutant of ARF3 ofM. polymorpha, Tak-1 background
Genetic reagent (M. polymorpha)arf3ge2-1this papermutant of ARF3 ofM. polymorpha, Tak-1 background
Genetic reagent (M. polymorpha)arf1-4PMID: 29016901
Recombinant DNA reagentpMpGE_En03Addgene71535
Recombinant DNA reagentpMpGE_010Addgene71536
Recombinant DNA reagentpHKDW081this paperentry vector containing sgRNA for nciaa
Recombinant DNA reagentpHKDW084this paperentry vector containing sgRNA for ncarf
Recombinant DNA reagentpHKDW004this paperentry vector containing sgRNA for arf3ge1
Recombinant DNA reagentpHKDW005this paperentry vector containing sgRNA for arf3ge2
Commercial assay or kitTRIzol reagentThermo Fisher15596018
Commercial assay or kitRNeasy Plant Mini kitQIAGEN74904
Commercial assay or kitRNase-free DNase I setQIAGEN79254
Commercial assay or kitiScript cDNA Synthesis KitBio-Rad1708891
Commercial assay or kitiQ SYBR Green SupermixBio-Rad1708886
Chemical compound, drugGamborg B5 mediumDuchefa BiochemieG0209
Software, algorithmBLAST + v2.2.28PMID: 20003500
Software, algorithmTransDecoder (ver2.0.1)
Software, algorithmInterProScan database (ver5.19–58.0)PMID: 24451626
Software, algorithmMAFFTPMID: 23329690
Software, algorithmPhyutility (ver2.2.6)PMID: 18227120
Software, algorithmPartitionFinder (ver1.1.1)PMID: 22319168
Software, algorithmRAxML (ver8.1.20)PMID: 24451623
Software, algorithmiTOL (ver3)PMID: 27095192
Software, algorithmTrinityPMID: 23845962
Software, algorithmBowtie2PMID: 22388286
Software, algorithmCorsetPMID: 25063469
Software, algorithmDEseq2PMID: 25516281
Software, algorithmClustalOmegaPMID: 21988835
Software, algorithmEspriptPMID: 24753421
Software, algorithmModeller v9.17PMID: 27322406
Software, algorithmPyMOLSchrödinger
Software, algorithmCytoscapePMID: 14597658

Additional files

Supplementary file 1

Species used in phylogenomic analysis.

https://doi.org/10.7554/eLife.33399.022
Supplementary file 2

Multiple sequence alignments used in the study.

https://doi.org/10.7554/eLife.33399.023
Supplementary file 3

Primers used in this study.

https://doi.org/10.7554/eLife.33399.024
Supplementary file 4

Detail for network of up-regulated genes.

https://doi.org/10.7554/eLife.33399.025
Supplementary file 5

Detail for network of down-regulated genes.

https://doi.org/10.7554/eLife.33399.026
Transparent reporting form
https://doi.org/10.7554/eLife.33399.027

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Sumanth K Mutte
  2. Hirotaka Kato
  3. Carl Rothfels
  4. Michael Melkonian
  5. Gane Ka-Shu Wong
  6. Dolf Weijers
(2018)
Origin and evolution of the nuclear auxin response system
eLife 7:e33399.
https://doi.org/10.7554/eLife.33399