Evolutionarily Conserved Principles Predict 3D Chromatin Organization.

Topologically associating domains (TADs), CTCF loop domains, and A/B compartments have been identified as important structural and functional components of 3D chromatin organization, yet the relationship between these features is not well understood. Using high-resolution Hi-C and HiChIP, we show that Drosophila chromatin is organized into domains we term compartmental domains that correspond precisely with A/B compartments at high resolution. We find that transcriptional state is a major predictor of Hi-C contact maps in several eukaryotes tested, including C. elegans and A. thaliana. Architectural proteins insulate compartmental domains by reducing interaction frequencies between neighboring regions in Drosophila, but CTCF loops do not play a distinct role in this organism. In mammals, compartmental domains exist alongside CTCF loop domains to form topological domains. The results suggest that compartmental domains are responsible for domain structure in all eukaryotes, with CTCF playing an important role in domain formation in mammals.

Different DNA repair pathways are involved in single-strand break-induced genomic changes in plants.

In nature, single-strand breaks (SSBs) in DNA occur more frequently (by orders of magnitude) than double-strand breaks (DSBs). SSBs induced by the CRISPR/Cas9 nickase at a distance of 50-100 bp on opposite strands are highly mutagenic, leading to insertions/deletions (InDels), with insertions mainly occurring as direct tandem duplications. As short tandem repeats are overrepresented in plant genomes, this mechanism seems to be important for genome evolution. We investigated the distance at which paired 5'-overhanging SSBs are mutagenic and which DNA repair pathways are essential for insertion formation in Arabidopsis thaliana. We were able to detect InDel formation up to a distance of 250 bp, although with much reduced efficiency. Surprisingly, the loss of the classical nonhomologous end joining (NHEJ) pathway factors KU70 or DNA ligase 4 completely abolished tandem repeat formation. The microhomology-mediated NHEJ factor POLQ was required only for patch-like insertions, which are well-known from DSB repair as templated insertions from ectopic sites. As SSBs can also be repaired using homology, we furthermore asked whether the classical homologous recombination (HR) pathway is involved in this process in plants. The fact that RAD54 is not required for homology-mediated SSB repair demonstrates that the mechanisms for DSB- and SSB-induced HR differ in plants.

Revealing atomic-scale molecular diffusion of a plant-transcription factor WRKY domain protein along DNA.

Transcription factor (TF) target search on genome is highly essential for gene expression and regulation. High-resolution determination of TF diffusion along DNA remains technically challenging. Here, we constructed a TF model system using the plant WRKY domain protein in complex with DNA from crystallography and demonstrated microsecond diffusion dynamics of WRKY on DNA by employing all-atom molecular-dynamics (MD) simulations. Notably, we found that WRKY preferentially binds to one strand of DNA with significant energetic bias compared with the other, or nonpreferred strand. The preferential DNA-strand binding becomes most prominent in the static process, from nonspecific to specific DNA binding, but less distinct during diffusive movements of the domain protein on the DNA. Remarkably, without employing acceleration forces or bias, we captured a complete one-base-pair stepping cycle of the protein tracking along major groove of DNA with a homogeneous poly-adenosine sequence, as individual hydrogen bonds break and reform at the protein-DNA binding interface. Further DNA-groove tracking motions of the protein forward or backward, with occasional sliding as well as strand crossing to minor groove of DNA, were also captured. The processive diffusion of WRKY along DNA has been further sampled via coarse-grained MD simulations. The study thus provides structural dynamics details on diffusion of a small TF domain protein, suggests how the protein approaches a specific recognition site on DNA, and supports further high-precision experimental detection. The stochastic movements revealed in the TF diffusion also provide general clues about how other protein walkers step and slide along DNA.

Pathway conversion enables a double-lock mechanism to maintain DNA methylation and genome stability.

The CMT2 and RNA-directed DNA methylation (RdDM) pathways have been proposed to separately maintain CHH methylation in specific regions of the Arabidopsis thaliana genome. Here, we show that dysfunction of the chromatin remodeler DDM1 causes hundreds of genomic regions to switch from CMT2 dependency to RdDM dependency in DNA methylation. These converted loci are enriched at the edge regions of long transposable elements (TEs). Furthermore, we found that dysfunction in both DDM1 and RdDM causes strong reactivation of TEs and a burst of TE transposition in the first generation of mutant plants, indicating that the DDM1 and RdDM pathways together are critical to maintaining TE repression and protecting genomic stability. Our findings reveal the existence of a pathway conversion-based backup mechanism to guarantee the maintenance of DNA methylation and genome integrity.

The intersection of DNA replication with antisense 3' RNA processing in Arabidopsis FLC chromatin silencing.

How noncoding transcription influences chromatin states is still unclear. The Arabidopsis floral repressor gene FLC is quantitatively regulated through an antisense-mediated chromatin silencing mechanism. The FLC antisense transcripts form a cotranscriptional R-loop that is dynamically resolved by RNA 3' processing factors (FCA and FY), and this is linked to chromatin silencing. Here, we investigate this silencing mechanism and show, using single-molecule DNA fiber analysis, that FCA and FY are required for unimpeded replication fork progression across the Arabidopsis genome. We then employ the chicken DT40 cell line system, developed to investigate sequence-dependent replication and chromatin inheritance, and find that FLC R-loop sequences have an orientation-dependent ability to stall replication forks. These data suggest a coordination between RNA 3' processing of antisense RNA and replication fork progression in the inheritance of chromatin silencing at FLC.

Back-spliced RNA from retrotransposon binds to centromere and regulates centromeric chromatin loops in maize.

In most plants, centromeric DNA contains highly repetitive sequences, including tandem repeats and retrotransposons; however, the roles of these sequences in the structure and function of the centromere are unclear. Here, we found that multiple RNA sequences from centromeric retrotransposons (CRMs) were enriched in maize (Zea mays) centromeres, and back-spliced RNAs were generated from CRM1. We identified 3 types of CRM1-derived circular RNAs with the same back-splicing site based on the back-spliced sequences. These circular RNAs bound to the centromere through R-loops. Two R-loop sites inside a single circular RNA promoted the formation of chromatin loops in CRM1 regions. When RNA interference (RNAi) was used to target the back-splicing site of the circular CRM1 RNAs, the levels of R-loops and chromatin loops formed by these circular RNAs decreased, while the levels of R-loops produced by linear RNAs with similar binding sites increased. Linear RNAs with only one R-loop site could not promote chromatin loop formation. Higher levels of R-loops and lower levels of chromatin loops in the CRM1 regions of RNAi plants led to a reduced localization of the centromeric H3 variant (CENH3). Our work reveals centromeric chromatin organization by circular CRM1 RNAs via R-loops and chromatin loops, which suggested that CRM1 elements might help build a suitable chromatin environment during centromere evolution. These results highlight that R-loops are integral components of centromeric chromatin and proper centromere structure is essential for CENH3 localization.

SDC mediates DNA methylation-controlled clock pace by interacting with ZTL in Arabidopsis.

Molecular bases of eukaryotic circadian clocks mainly rely on transcriptional-translational feedback loops (TTFLs), while epigenetic codes also play critical roles in fine-tuning circadian rhythms. However, unlike histone modification codes that play extensive and well-known roles in the regulation of circadian clocks, whether DNA methylation (5mC) can affect the circadian clock, and the associated underlying molecular mechanisms, remains largely unexplored in many organisms. Here we demonstrate that global genome DNA hypomethylation can significantly lengthen the circadian period of Arabidopsis. Transcriptomic and genetic evidence demonstrate that SUPPRESSOR OF drm1 drm2 cmt3 (SDC), encoding an F-box containing protein, is required for the DNA hypomethylation-tuned circadian clock. Moreover, SDC can physically interact with another F-box containing protein ZEITLUPE (ZTL) to diminish its accumulation. Genetic analysis further revealed that ZTL and its substrate TIMING OF CAB EXPRESSION 1 (TOC1) likely act downstream of DNA methyltransferases to control circadian rhythm. Together, our findings support the notion that DNA methylation is important to maintain proper circadian pace in Arabidopsis, and further established that SDC links DNA hypomethylation with a proteolytic cascade to assist in tuning the circadian clock.

INT-Hi-C reveals distinct chromatin architecture in endosperm and leaf tissues of Arabidopsis.

Higher-order chromatin structure undergoes striking changes in response to various developmental and environmental signals, causing distinct cell types to adopt specific chromatin organization. High throughput chromatin conformation capture (Hi-C) allows studying higher-order chromatin structure; however, this technique requires substantial amounts of starting material, which has limited the establishment of cell type-specific higher-order chromatin structure in plants. To overcome this limitation, we established a protocol that is applicable to a limited amount of nuclei by combining the INTACT (isolation of nuclei tagged in specific cell types) method and Hi-C (INT-Hi-C). Using this INT-Hi-C protocol, we generated Hi-C data from INTACT purified endosperm and leaf nuclei. Our INT-Hi-C data from leaf accurately reiterated chromatin interaction patterns derived from conventional leaf Hi-C data. We found that the higher-order chromatin organization of mixed leaf tissues and endosperm differs and that DNA methylation and repressive histone marks positively correlate with the chromatin compaction level. We furthermore found that self-looped interacting genes have increased expression in leaves and endosperm and that interacting intergenic regions negatively impact on gene expression in the endosperm. Last, we identified several imprinted genes involved in long-range and trans interactions exclusively in endosperm. Our study provides evidence that the endosperm adopts a distinct higher-order chromatin structure that differs from other cell types in plants and that chromatin interactions influence transcriptional activity.

Structural determinants for NF-Y subunit organization and NF-Y/DNA association in plants.

NF-Y transcription factor comprises three subunits: NF-YA, NF-YB and NF-YC. NF-YB and NF-YC dimerize through their histone fold domain (HFD), which can bind DNA in a non-sequence-specific fashion while serving as a scaffold for NF-YA trimerization. Upon trimerization, NF-YA specifically recognizes the CCAAT box sequence on promoters and enhancers. In plants, each NF-Y subunit is encoded by several genes giving rise to hundreds of potential heterotrimeric combinations. In addition, plant NF-YBs and NF-YCs interact with other protein partners to recognize a plethora of genomic motifs, as the CCT protein family that binds CORE sites. The NF-Y subunit organization and its DNA-binding properties, together with the NF-Y HFD capacity to adapt different protein modules, represent plant-specific features that play a key role in development, growth and reproduction. Despite their relevance, these features are still poorly understood at the molecular level. Here, we present the structures of Arabidopsis and rice NF-YB/NF-YC dimers, and of an Arabidopsis NF-Y trimer in complex with the FT CCAAT box, together with biochemical data on NF-Y mutants. The dimeric structures identify the key residues for NF-Y HFD stabilization. The NF-Y/DNA structure and the mutation experiments shed light on HFD trimerization interface properties and the NF-YA sequence appetite for the bases flanking the CCAAT motif. These data explain the logic of plant NF-Y gene expansion: the trimerization adaptability and the flexible DNA-binding rules serve the scopes of accommodating the large number of NF-YAs, CCTs and possibly other NF-Y HFD binding partners and a diverse audience of genomic motifs.

Hi-C analysis in Arabidopsis identifies the KNOT, a structure with similarities to the flamenco locus of Drosophila.

Chromosomes are folded, spatially organized, and regulated by epigenetic marks. How chromosomal architecture is connected to the epigenome is not well understood. We show that chromosomal architecture of Arabidopsis is tightly linked to the epigenetic state. Furthermore, we show how physical constraints, such as nuclear size, correlate with the folding principles of chromatin. We also describe a nuclear structure, termed KNOT, in which genomic regions of all five Arabidopsis chromosomes interact. These KNOT ENGAGED ELEMENT (KEE) regions represent heterochromatic islands within euchromatin. Similar to PIWI-interacting RNA clusters, such as flamenco in Drosophila, KEEs represent preferred landing sites for transposable elements, which may be part of a transposon defense mechanism in the Arabidopsis nucleus.

Genome-wide Hi-C analyses in wild-type and mutants reveal high-resolution chromatin interactions in Arabidopsis.

Chromosomes form 3D structures that are critical to the regulation of cellular and genetic processes. Here, we present a study of global chromatin interaction patterns in Arabidopsis thaliana. Our genome-wide approach confirmed interactions that were previously observed by other methods as well as uncovered long-range interactions such as those among small heterochromatic regions embedded in euchromatic arms. We also found that interactions are correlated with various epigenetic marks that are localized in active or silenced chromatin. Arabidopsis chromosomes do not contain large local interactive domains that resemble the topological domains described in animals but, instead, contain relatively small interactive regions scattered around the genome that contain H3K27me3 or H3K9me2. We generated interaction maps in mutants that are defective in specific epigenetic pathways and found altered interaction patterns that correlate with changes in the epigenome. These analyses provide further insights into molecular mechanisms of epigenetic regulation of the genome.

Mechanism of DNA methylation-directed histone methylation by KRYPTONITE.

In Arabidopsis, CHG DNA methylation is controlled by the H3K9 methylation mark through a self-reinforcing loop between DNA methyltransferase CHROMOMETHYLASE3 (CMT3) and H3K9 histone methyltransferase KRYPTONITE/SUVH4 (KYP). We report on the structure of KYP in complex with methylated DNA, substrate H3 peptide, and cofactor SAH, thereby defining the spatial positioning of the SRA domain relative to the SET domain. The methylated DNA is bound by the SRA domain with the 5mC flipped out of the DNA, while the H3(1-15) peptide substrate binds between the SET and post-SET domains, with the ε-ammonium of K9 positioned adjacent to bound SAH. These structural insights, complemented by functional data on key mutants of residues lining the 5mC and H3K9-binding pockets within KYP, establish how methylated DNA recruits KYP to the histone substrate. Together, the structures of KYP and previously reported CMT3 complexes provide insights into molecular mechanisms linking DNA and histone methylation.

Structural insights into target DNA recognition by R2R3-MYB transcription factors.

As the largest group of MYB family transcription factors, R2R3-MYB proteins play essential roles during plant growth and development. However, the structural basis underlying how R2R3-MYBs recognize the target DNA remains elusive. Here, we report the crystal structure of Arabidopsis WEREWOLF (WER), an R2R3-MYB protein, in complex with its target DNA. Structural analysis showed that the third α-helices in both the R2 and R3 repeats of WER fit in the major groove of the DNA, specifically recognizing the DNA motif 5'-AACNGC-3'. In combination with mutagenesis, in vitro binding and in vivo luciferase assays, we showed that K55, N106, K109 and N110 are critical for the function of WER. Although L59 of WER is not involved in DNA binding in the structure, ITC analysis suggested that L59 plays an important role in sensing DNA methylation at the fifth position of cytosine (5mC). Like 5mC, methylation at the sixth position of adenine (6mA) in the AAC element also inhibits the interaction between WER and its target DNA. Our study not only unravels the molecular basis of how WER recognizes its target DNA, but also suggests that 5mC and 6mA modifications may block the interaction between R2R3-MYB transcription factors and their target genes.

DNA-binding properties of the MADS-domain transcription factor SEPALLATA3 and mutant variants characterized by SELEX-seq.

KEY MESSAGE: We studied the DNA-binding profile of the MADS-domain transcription factor SEPALLATA3 and mutant variants by SELEX-seq. DNA-binding characteristics of SEPALLATA3 mutant proteins lead us to propose a novel DNA-binding mode. MIKC-type MADS-domain proteins, which function as essential transcription factors in plant development, bind as dimers to a 10-base-pair AT-rich motif termed CArG-box. However, this consensus motif cannot fully explain how the abundant family members in flowering plants can bind different target genes in specific ways. The aim of this study was to better understand the DNA-binding specificity of MADS-domain transcription factors. Also, we wanted to understand the role of a highly conserved arginine residue for binding specificity of the MADS-domain transcription factor family. Here, we studied the DNA-binding profile of the floral homeotic MADS-domain protein SEPALLATA3 by performing SELEX followed by high-throughput sequencing (SELEX-seq). We found a diverse set of bound sequences and could estimate the in vitro binding affinities of SEPALLATA3 to a huge number of different sequences. We found evidence for the preference of AT-rich motifs as flanking sequences. Whereas different CArG-boxes can act as SEPALLATA3 binding sites, our findings suggest that the preferred flanking motifs are almost always the same and thus mostly independent of the identity of the central CArG-box motif. Analysis of SEPALLATA3 proteins with a single amino acid substitution at position 3 of the DNA-binding MADS-domain further revealed that the conserved arginine residue, which has been shown to be involved in a shape readout mechanism, is especially important for the recognition of nucleotides at positions 3 and 8 of the CArG-box motif. This leads us to propose a novel DNA-binding mode for SEPALLATA3, which is different from that of other MADS-domain proteins known.

DNA sequence properties that predict susceptibility to epiallelic switching.

Transgenerationally heritable epialleles are defined by the stable propagation of alternative transcriptional states through mitotic and meiotic cell cycles. Given that the propagation of DNA methylation at CpG sites, mediated in Arabidopsis by MET1, plays a central role in epigenetic inheritance, we examined genomewide DNA methylation in partial and complete loss-of-function met1 mutants. We interpreted the data in relation to transgenerational epiallelic stability, which allowed us to classify chromosomal targets of epigenetic regulation into (i) single copy and methylated exclusively at CpGs, readily forming epialleles, and (ii) transposon-derived, methylated at all cytosines, which may or may not form epialleles. We provide evidence that DNA sequence features such as density of CpGs and genomic repetitiveness of the loci predispose their susceptibility to epiallelic switching. The importance and predictive power of these genetic features were confirmed by analyses of common epialleles in natural Arabidopsis accessions, epigenetic recombinant inbred lines (epiRILs) and also verified in rice.

Historical biogeography of Pomaderris (Rhamnaceae): Continental vicariance in Australia and repeated independent dispersals to New Zealand.

AIM: Gondwanan biogeographic patterns include a combination of old vicariance events following the breakup of the supercontinent, and more recent long-distance dispersals across the southern landmasses. Floristic relationships between Australia and New Zealand have mostly been attributed to recent dispersal events rather than vicariance. We assessed the biogeographic history of Pomaderris (Rhamnaceae), which occurs in both Australia and New Zealand, by constructing a time-calibrated molecular phylogeny to infer (1) phylogenetic relationships and (2) the relative contributions of vicariance and dispersal events in the biogeographic history of the genus. LOCATION: Australia and New Zealand. METHODS: Using hybrid capture and high throughput sequencing, we generated nuclear and plastid data sets to estimate phylogenetic relationships and fossil calibrated divergence time estimates for Pomaderris. BioGeoBEARS and biogeographical stochastic mapping (BSM) were used to assess the ancestral area of the genus and the relative contributions of vicariance vs dispersal, and the directionality of dispersal events. RESULTS: Our analyses indicate that Pomaderris originated in the Oligocene and had a widespread Australian distribution. Vicariance of western and eastern Australian clades coincides with the uplift of the Nullarbor Plain c. 14 Ma, followed by subsequent in-situ and within-biome diversification with little exchange across regions. A rapid radiation of southeastern Australian taxa beginning c. 10 Ma was the source for at least six independent long-distance dispersal events to New Zealand during the Pliocene-Pleistocene. MAIN CONCLUSIONS: Our study demonstrates the importance of dispersal in explaining not only the current cross-Tasman distributions of Pomaderris, but for the New Zealand flora more broadly. The pattern of multiple independent long-distance dispersal events for Pomaderris, without significant radiation within New Zealand, is congruent with other lowland plant groups, suggesting that this biome has a different evolutionary history compared with the younger alpine flora of New Zealand, which exhibits extensive radiations often following single long distance dispersal events.

Phylogenomic analyses based on genome-skimming data reveal cyto-nuclear discordance in the evolutionary history of Cotoneaster (Rosaceae).

As a consequence of hybridization, polyploidization, and apomixis, the genus Cotoneaster (Rosaceae) represents one of the most complicated and controversial lineages in Rosaceae, with ca. 370 species which have been classified into two subgenera and several sections, and is notorious for its taxonomic difficulty. The infrageneric relationships and taxonomy of Cotoneaster have remained poorly understood. Previous studies have focused mainly on natural hybridization involving only several species, and phylogeny based on very limited markers. In the present study, the sequences of complete chloroplast genomes and 204 low-copy nuclear genes of 72 accessions, representing 69 species as ingroups, were used to conduct the most comprehensive phylogenetic analysis so far for Cotoneaster. Based on the sequences of complete chloroplast genomes and many nuclear genes, our analyses yield two robust phylogenetic trees respectively. Chloroplast genome and nuclear data confidently resolved relationships of this genus into two major clades which largely supported current classification based on morphological evidence. However, conflicts between the chloroplast genome and low-copy nuclear phylogenies were observed in both the species level and clade level. Cyto-nuclear discordance in the phylogeny could be caused by frequent hybridization events and incomplete sorting lineage (ILS). In addition, our divergence-time analysis revealed an evolutionary radiation of the genus from late Miocene to date.

Combination of Sanger and target-enrichment markers supports revised generic delimitation in the problematic 'Urera clade' of the nettle family (Urticaceae).

Urera Gaudich, s.l. is a pantropical genus comprising c. 35 species of trees, shrubs, and vines. It has a long history of taxonomic uncertainty, and is repeatedly recovered as polyphyletic within a poorly resolved complex of genera in the Urticeae tribe of the nettle family (Urticaceae). To provide generic delimitations concordant with evolutionary history, we use increased taxonomic and genomic sampling to investigate phylogenetic relationships among Urera and associated genera. A cost-effective two-tier genome-sampling approach provides good phylogenetic resolution by using (i) a taxon-dense sample of Sanger sequence data from two barcoding regions to recover clades of putative generic rank, and (ii) a genome-dense sample of target-enrichment data for a subset of representative species from each well-supported clade to resolve relationships among them. The results confirm the polyphyly of Urera s.l. with respect to the morphologically distinct genera Obetia, Poikilospermum and Touchardia. Afrotropic members of Urera s.l. are recovered in a clade sister to the xerophytic African shrubs Obetia; and Hawaiian ones with Touchardia, also from Hawaii. Combined with distinctive morphological differences between Neotropical and African members of Urera s.l., these results lead us to resurrect the previously synonymised name Scepocarpus Wedd. for the latter. The new species epiphet Touchardia oahuensis T.Wells & A.K. Monro is offered as a replacement name for Touchardia glabra non H.St.John, and subgenera are created within Urera s.s. to account for the two morphologically distinct Neotropical clades. This new classification minimises taxonomic and nomenclatural disruption, while more accurately reflecting evolutionary relationships within the group.

Maize Endochitinase Expression in Response to Fall Armyworm Herbivory.

A large percentage of crop loss is due to insect damage, especially caterpillar damage. Plant chitinases are considered excellent candidates to combat these insects since they can degrade chitin in peritrophic matrix (PM), an important protective structure in caterpillar midgut. Compared to chemical insecticides, chitinases could improve host plant resistance and be both economically and environmentally advantageous. The focus of this research was to find chitinase candidates that could improve plant resistance by effectively limiting caterpillar damage. Five classes of endochitinase (I-V) genes were characterized in the maize genome, and we isolated and cloned four chitinase genes (chitinase A, chitinase B, chitinase I, and PRm3) present in two maize (Zea mays L.) inbred lines Mp708 and Tx601, with different levels of resistance to caterpillar pests. We also investigated the expression of these maize chitinases in response to fall armyworm (Spodoptera frugiperda, FAW) attack. The results indicated that both chitinase transcript abundance and enzymatic activity increased in response to FAW feeding and mechanical wounding. Furthermore, chitinases retained activity inside the caterpillar midgut and enzymatic activity was detected in the food bolus and frass. When examined under scanning electron microscopy, PMs from Tx601-fed caterpillars showed structural damage when compared to diet controls. Analysis of chitinase transcript abundance after caterpillar feeding and proteomic analysis of maize leaf trichomes in the two inbreds implicated chitinase PRm3 found in Tx601 as a potential insecticidal protein.

Structural basis of dimerization and dual W-box DNA recognition by rice WRKY domain.

In rice, the critical regulator of the salicylic acid signalling pathway is OsWRKY45, a transcription factor (TF) of the WRKY TF family that functions by binding to the W-box of gene promoters, but the structural basis of OsWRKY45/W-box DNA recognition is unknown. Here, we show the crystal structure of the DNA binding domain of OsWRKY45 (OsWRKY45-DBD, i.e. the WRKY and zinc finger domain) in complex with a W-box DNA. Surprisingly, two OsWRKY45-DBD molecules exchange ß4-ß5 strands to form a dimer. The domain swapping occurs at the hinge region between the ß3 and ß4 strands, and is bridged and stabilized by zinc ion via coordinating residues from different chains. The dimer contains two identical DNA binding domains that interact with the major groove of W-box DNA. In addition to hydrophobic and direct hydrogen bonds, water mediated hydrogen bonds are also involved in base-specific interaction between protein and DNA. Finally, we discussed the cause and consequence of domain swapping of OsWRKY45-DBD, and based on our work and that of previous studies present a detailed mechanism of W-box recognition by WRKY TFs. This work reveals a novel dimerization and DNA-binding mode of WRKY TFs, and an intricate picture of the WRKY/W-box DNA recognition.