Highlights d A database combining genomic information and chromatin profiles for Marchantia d Correlations between chromatin marks and transcription are conserved in land plants d A significant portion of constitutive heterochromatin is marked by H3K27me3 d Insights into the evolution of TAD organization in plants

Bryophytes (liverworts, hornworts and mosses) comprise the three earliest diverging lineages of land plants (embryophytes). Marchantia polymorpha, a complex thalloid Marchantiopsida liverwort that has been developed into a model genetic system, occupies a key phylogenetic position. Therefore, M. polymorpha is useful in studies aiming to elucidate the evolution of gene regulation mechanisms in plants. In this study, we used computational, transcriptomic, small RNA and degradome analyses to characterize microRNA (miRNA)-mediated pathways of gene regulation in M. polymorpha. The data have been integrated into the open access ContigViews-miRNA platform for further reference. In addition to core components of the miRNA pathway, 129 unique miRNA sequences, 11 of which could be classified into seven miRNA families that are conserved in embryophytes (miR166a, miR390, miR529c, miR171-3p, miR408a, miR160 and miR319a), were identified. A combination of computational and degradome analyses allowed us to identify and experimentally validate 249 targets. In some cases, the target genes are orthologous to those of other embryophytes, but in other cases, the conserved miRNAs target either paralogs or members of different gene families. In addition, the newly discovered Mpo-miR11707.1 and Mpo-miR11707.2 are generated from a common precursor and target MpARGONAUTE1 (LW1759). Two other newly discovered miRNAs, Mpo-miR11687.1 and Mpo-miR11681.1, target the MADS-box transcription factors MpMADS1 and MpMADS2, respectively. Interestingly, one of the pentatricopeptide repeat (PPR) gene family members, MpPPR_66 (LW9825), the protein products of which are generally involved in various steps of RNA metabolism, has a long stem–loop transcript that can generate Mpo-miR11692.1 to autoregulate MpPPR_66 (LW9825) mRNA. This study provides a foundation for further investigations of the RNA-mediated silencing mechanism in M. polymorpha as well as of the evolution of this gene silencing pathway in embryophytes.

A plethora of developmental and physiological processes in land plants is influenced by auxin, to a large extent via alterations in gene expression by AUXIN RESPONSE FACTORs (ARFs). The canonical auxin transcriptional response system is a land plant innovation, however, charophycean algae possess orthologues of at least some classes of ARF and AUXIN/INDOLE-3-ACETIC ACID (AUX/IAA) genes, suggesting that elements of the canonical land plant system existed in an ancestral alga. We reconstructed the phylogenetic relationships between streptophyte ARF and AUX/IAA genes and functionally characterized the solitary class C ARF, MpARF3, in Marchantia polymorpha. Phylogenetic analyses indicate that multiple ARF classes, including class C ARFs, existed in an ancestral alga. Loss- and gain-of-function MpARF3 alleles result in pleiotropic effects in the gametophyte, with MpARF3 inhibiting differentiation and developmental transitions in multiple stages of the life cycle. Although loss-of-function Mparf3 and Mpmir160 alleles respond to exogenous auxin treatments, strong miR-resistant MpARF3 alleles are auxin-insensitive, suggesting that class C ARFs act in a context-dependent fashion. We conclude that two modules independently evolved to regulate a pre-existing ARF transcriptional network. Whereas the auxin-TIR1-AUX/IAA pathway evolved to repress class A/B ARF activity, miR160 evolved to repress class C ARFs in a dynamic fashion.

Summary Phytoplasmas are bacterial plant pathogens which can induce severe symptoms including dwarfism, phyllody and virescence in an infected plant. Because phytoplasmas infect many important crops such as peanut and papaya they have caused serious agricultural losses. The phytoplasmal effector causing phyllody 1 (PHYL1) is an important phytoplasmal pathogenic factor which affects the biological function of MADS transcription factors by interacting with their K (keratin‐like) domain, thus resulting in abnormal plant developments such as phyllody. Until now, lack of information on the structure of PHYL1 has prevented a detailed understanding of the binding mechanism between PHYL1 and the MADS transcription factors. Here, we present the crystal structure of PHYL1 from peanut witches’‐broom phytoplasma (PHYL1PnWB). This protein was found to fold into a unique α‐helical hairpin with exposed hydrophobic residues on its surface that may play an important role in its biological function. Using proteomics approaches, we propose a binding mode of PHYL1PnWB with the K domain of the MADS transcription factor SEPALLATA3 (SEP3_K) and identify the residues of PHYL1PnWB that are important for this interaction. Furthermore, using surface plasmon resonance we measure the binding strength of PHYL1PnWB proteins to SEP3_K. Lastly, based on confocal images, we found that α‐helix 2 of PHYL1PnWB plays an important role in PHYL1‐mediated degradation of SEP3. Taken together, these results provide a structural understanding of the specific binding mechanism between PHYL1PnWB and SEP3_K.