Hornworts comprise a bryophyte lineage that diverged from other extant land plants >400 million years ago and bears unique biological features, including a distinct sporophyte architecture, cyanobacterial symbiosis and a pyrenoid-based carbonconcentrating mechanism (CCM). Here, we provide three high-quality genomes of Anthoceros hornworts. Phylogenomic analyses place hornworts as a sister clade to liverworts plus mosses with high support. The Anthoceros genomes lack repeat-dense centromeres as well as whole-genome duplication, and contain a limited transcription factor repertoire. Several genes involved in angiosperm meristem and stomatal function are conserved in Anthoceros and upregulated during sporophyte development, suggesting possible homologies at the genetic level. We identified candidate genes involved in cyanobacterial symbiosis and found that LCIB, a Chlamydomonas CCM gene, is present in hornworts but absent in other plant lineages, implying a possible conserved role in CCM function. We anticipate that these hornwort genomes will serve as essential references for future hornwort research and comparative studies across land plants.

Summary Hornworts are crucial to understand the phylogeny of early land plants. The emergence of ‘reverse’ U‐to‐C RNA editing accompanying the widespread C‐to‐U RNA editing in plant chloroplasts and mitochondria may be a molecular synapomorphy of a hornwort–tracheophyte clade. C‐to‐U RNA editing is well understood after identification of many editing factors in models like Arabidopsis thaliana and Physcomitrella patens, but there is no plant model yet to investigate U‐to‐C RNA editing. The hornwort Anthoceros agrestis is now emerging as such a model system. We report on the assembly and analyses of the A. agrestis chloroplast and mitochondrial genomes, their transcriptomes and editomes, and a large nuclear gene family encoding pentatricopeptide repeat (PPR) proteins likely acting as RNA editing factors. Both organelles in A. agrestis feature high amounts of RNA editing, with altogether > 1100 sites of C‐to‐U and 1300 sites of U‐to‐C editing. The nuclear genome reveals > 1400 genes for PPR proteins with variable carboxyterminal DYW domains. We observe significant variants of the ‘classic’ DYW domain, in the meantime confirmed as the cytidine deaminase for C‐to‐U editing, and discuss the first attractive candidates for reverse editing factors given their excellent matches to U‐to‐C editing targets according to the PPR‐RNA binding code.

Summary Marchantia polymorpha has recently become a prime model for cellular, evo‐devo, synthetic biological, and evolutionary investigations. We present a pseudomolecule‐scale assembly of the M. polymorpha genome, making comparative genome structure analysis and classical genetic mapping approaches feasible. We anchored 88% of the M. polymorpha draft genome to a high‐density linkage map resulting in eight pseudomolecules. We found that the overall genome structure of M. polymorpha is in some respects different from that of the model moss Physcomitrella patens. Specifically, genome collinearity between the two bryophyte genomes and vascular plants is limited, suggesting extensive rearrangements since divergence. Furthermore, recombination rates are greatest in the middle of the chromosome arms in M. polymorpha like in most vascular plant genomes, which is in contrast with P. patens where recombination rates are evenly distributed along the chromosomes. Nevertheless, some other properties of the genome are shared with P. patens. As in P. patens, DNA methylation in M. polymorpha is spread evenly along the chromosomes, which is in stark contrast with the angiosperm model Arabidopsis thaliana, where DNA methylation is strongly enriched at the centromeres. Nevertheless, DNA methylation and recombination rate are anticorrelated in all three species. Finally, M. polymorpha and P. patens centromeres are of similar structure and marked by high abundance of retroelements unlike in vascular plants. Taken together, the highly contiguous genome assembly we present opens unexplored avenues for M. polymorpha research by linking the physical and genetic maps, making novel genomic and genetic analyses, including map‐based cloning, feasible.

Understanding the molecular basis of morphological shifts is a fundamental question of evolutionary biology. New morphologies may arise through the birth/death of genes (gene gain/loss) or by reutilizing existing gene sets. Yet, the relative contribution of these two processes to radical morphological shifts is still poorly understood. Here, we use the model system of two mosses, Funaria hygrometrica and Physcomitrium (Physcomitrella) patens, to investigate the molecular mechanisms underlying contrasting sporophyte architectures. We used comparative analysis of time-series expression data for four stages of sporophyte development in both species to address this question in detail. We found that large-scale differences in sporophytic architecture are mainly governed by orthologous (i.e., shared) genes frequently experiencing temporal gene expression shifts between the two species. While the absolute number of species-specific genes expressed during sporophyte development is somewhat smaller, we observed a significant increase of their proportion in preferentially sporophyte expressed genes, suggesting a fundamental role in the sporophyte phase. However, further functional studies are necessary to determine their contribution to diverging sporophyte morphologies. Our results add to the growing set of studies suggesting that radical changes in morphology may rely on the heterochronic expression of conserved regulators.

Premise A detailed protocol for the protoplast transformation of hornwort tissue is not yet available, limiting molecular biological investigations of these plants and comparative analyses with other bryophytes, which display a gametophyte‐dominant life cycle and are critical to understanding the evolution of key land plant traits. Methods and Results We describe a detailed protocol to isolate and transiently transform protoplasts of the model hornwort Anthoceros agrestis . The digestion of liquid cultures with Driselase yields a high number of viable protoplasts suitable for polyethylene glycol (PEG)‐mediated transformation. We also report early signs of protoplast regeneration, such as chloroplast division and cell wall reconstitution. Conclusions This protocol represents a straightforward method for isolating and transforming A. agrestis protoplasts that is less laborious than previously described approaches. In combination with the recently developed stable genome transformation technique, this work further expands the prospects of functional studies in this model hornwort.