Insights into Land Plant Evolution Garnered from the Marchantia polymorpha Genome.

The evolution of land flora transformed the terrestrial environment. Land plants evolved from an ancestral charophycean alga from which they inherited developmental, biochemical, and cell biological attributes. Additional biochemical and physiological adaptations to land, and a life cycle with an alternation between multicellular haploid and diploid generations that facilitated efficient dispersal of desiccation tolerant spores, evolved in the ancestral land plant. We analyzed the genome of the liverwort Marchantia polymorpha, a member of a basal land plant lineage. Relative to charophycean algae, land plant genomes are characterized by genes encoding novel biochemical pathways, new phytohormone signaling pathways (notably auxin), expanded repertoires of signaling pathways, and increased diversity in some transcription factor families. Compared with other sequenced land plants, M. polymorpha exhibits low genetic redundancy in most regulatory pathways, with this portion of its genome resembling that predicted for the ancestral land plant. PAPERCLIP.

The landscape of transcription factor promoter activity during vegetative development in Marchantia.

Transcription factors (TFs) are essential for the regulation of gene expression and cell fate determination. Characterizing the transcriptional activity of TF genes in space and time is a critical step toward understanding complex biological systems. The vegetative gametophyte meristems of bryophytes share some characteristics with the shoot apical meristems of flowering plants. However, the identity and expression profiles of TFs associated with gametophyte organization are largely unknown. With only ∼450 putative TF genes, Marchantia (Marchantia polymorpha) is an outstanding model system for plant systems biology. We have generated a near-complete collection of promoter elements derived from Marchantia TF genes. We experimentally tested reporter fusions for all the TF promoters in the collection and systematically analyzed expression patterns in Marchantia gemmae. This allowed us to build a map of expression domains in early vegetative development and identify a set of TF-derived promoters that are active in the stem-cell zone. The cell markers provide additional tools and insight into the dynamic regulation of the gametophytic meristem and its evolution. In addition, we provide an online database of expression patterns for all promoters in the collection. We expect that these promoter elements will be useful for cell-type-specific expression, synthetic biology applications, and functional genomics.

The renaissance and enlightenment of Marchantia as a model system.

The liverwort Marchantia polymorpha has been utilized as a model for biological studies since the 18th century. In the past few decades, there has been a Renaissance in its utilization in genomic and genetic approaches to investigating physiological, developmental, and evolutionary aspects of land plant biology. The reasons for its adoption are similar to those of other genetic models, e.g. simple cultivation, ready access via its worldwide distribution, ease of crossing, facile genetics, and more recently, efficient transformation, genome editing, and genomic resources. The haploid gametophyte dominant life cycle of M. polymorpha is conducive to forward genetic approaches. The lack of ancient whole-genome duplications within liverworts facilitates reverse genetic approaches, and possibly related to this genomic stability, liverworts possess sex chromosomes that evolved in the ancestral liverwort. As a representative of one of the three bryophyte lineages, its phylogenetic position allows comparative approaches to provide insights into ancestral land plants. Given the karyotype and genome stability within liverworts, the resources developed for M. polymorpha have facilitated the development of related species as models for biological processes lacking in M. polymorpha.

Evolution in the Cycles of Life.

The life cycles of eukaryotes alternate between haploid and diploid phases, which are initiated by meiosis and gamete fusion, respectively. In both ascomycete and basidiomycete fungi and chlorophyte algae, the haploid-to-diploid transition is regulated by a pair of paralogous homeodomain protein encoding genes. That a common genetic program controls the haploid-to-diploid transition in phylogenetically disparate eukaryotic lineages suggests this may be the ancestral function for homeodomain proteins. Multicellularity has evolved independently in many eukaryotic lineages in either one or both phases of the life cycle. Organisms, such as land plants, exhibiting a life cycle whereby multicellular bodies develop in both the haploid and diploid phases are often referred to as possessing an alternation of generations. We review recent progress on understanding the genetic basis for the land plant alternation of generations and highlight the roles that homeodomain-encoding genes may have played in the evolution of complex multicellularity in this lineage.

KANADI promotes thallus differentiation and FR-induced gametangiophore formation in the liverwort Marchantia.

In angiosperms, KANADI transcription factors have roles in the sporophyte generation regulating tissue polarity, organogenesis and shade avoidance responses, but are not required during the gametophyte generation. Whether these roles are conserved in the gametophyte-dominant bryophyte lineages is unknown, which we examined by characterising the sole KANADI ortholog, MpKAN, in the liverwort Marchantia polymorpha. In contrast to angiosperm orthologs, MpKAN functions in the gametophyte generation in Marchantia, where it regulates apical branching and tissue differentiation, but does not influence tissue polarity in either generation. MpKAN can partially rescue the sporophyte polarity defects of kanadi mutants in Arabidopsis, indicating that MpKAN has conserved biochemical activity to its angiosperm counterparts. Mpkan loss-of-function plants display defects in far-red (FR) light responses. Mpkan plants have reduced FR-induced growth tropisms, have a delayed transition to sexual reproduction and fail to correctly form gametangiophores. Our results indicate that MpKAN is a modulator of FR responses, which may reflect a conserved role for KANADI across land plants. Under FR, MpKAN negatively regulates MpDELLA expression, suggesting that MpKAN and MpDELLA act in a pathway regulating FR responses, placing MpKAN in a gene regulatory network exhibiting similarities with those of angiosperms.

A transporter of 1-aminocyclopropane-1-carboxylic acid affects thallus growth and fertility in Marchantia polymorpha.

In seed plants, 1-aminocyclopropane-1-carboxylic acid (ACC) is the precursor of the plant hormone ethylene but also has ethylene-independent signaling roles. Nonseed plants produce ACC but do not efficiently convert it to ethylene. In Arabidopsis thaliana, ACC is transported by amino acid transporters, LYSINE HISTIDINE TRANSPORTER 1 (LHT1) and LHT2. In nonseed plants, LHT homologs have been uncharacterized. Here, we isolated an ACC-insensitive mutant (Mpain) that is defective in ACC uptake in the liverwort Marchantia polymorpha. Mpain contained a frameshift mutation (1 bp deletion) in the MpLHT1 coding sequence, and was complemented by expression of a wild-type MpLHT1 transgene. Additionally, ACC insensitivity was re-created in CRISPR/Cas9-Mplht1 knockout mutants. We found that MpLHT1 can also transport l-hydroxyproline and l-histidine. We examined the physiological functions of MpLHT1 in vegetative growth and reproduction based on mutant phenotypes. Mpain and Mplht1 plants were smaller and developed fewer gemmae cups compared to wild-type plants. Mplht1 mutants also had reduced fertility, and archegoniophores displayed early senescence. These findings reveal that MpLHT1 serves as an ACC and amino acid transporter in M. polymorpha and has diverse physiological functions. We propose that MpLHT1 contributes to homeostasis of ACC and other amino acids in M. polymorpha growth and reproduction.

Class C ARFs evolved before the origin of land plants and antagonize differentiation and developmental transitions in Marchantia polymorpha.

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.

Efficient and Inducible Use of Artificial MicroRNAs in Marchantia polymorpha.

We describe the efficient use of artificial microRNAs (amiRs) in Marchantia polymorpha using both endogenous and heterologous primary microRNA (pri-miR) hairpin backbones. Targeting of two transcription factor genes, MpARF1 and MpRR-B, mediating different hormonal responses, demonstrated that amiRs can create specific and reproducible physiological and morphological defects, facilitating interpretation of gene function. A third amiR was designed to target a gene encoding a component of the Polycomb recessive complex 2, MpE(z), and constitutive expression of this amiR results in sporeling lethality. Adaptation of an estrogen-inducible system allowed analysis of the phenotypic effects of induction of this amiR during other stages of the life cycle. We discuss the advantages and challenges of the use of amiRs as a tool for reverse genetic analysis in M. polymorpha.

A subclass II bHLH transcription factor in Marchantia polymorpha gives insight into the ancestral land plant trait of spore formation.

Sporopollenin is often said to be one of the toughest biopolymers known to man. The shift in dormancy cell wall deposition from around the diploid zygotes of charophycean algae to sporopollenin around the haploid spores of land plants essentially imparted onto land plants the gift of passive motility, a key acquisition that contributed to their vast and successful colonization across terrestrial habitats.1,2 A putative transcription factor controlling the land plant mode of sporopollenin deposition is the subclass II bHLHs, which are conserved and novel to land plants, with mutants of genes in angiosperms and mosses divulging roles relating to tapetum degeneration and spore development.3,4,5,6,7 We demonstrate that a subclass II bHLH gene, MpbHLH37, regulates sporopollenin biosynthesis and deposition in the model liverwort Marchantia polymorpha. Mpbhlh37 sporophytes show a striking loss of secondary wall deposits of the capsule wall, the elaters, and the spore exine, all while maintaining spore viability, identifying MpbHLH37 as a master regulator of secondary wall deposits of the sporophyte. Localization of MpbHLH37 to the capsule wall and elaters of the sporophyte directly designates these tissue types as a bona fide tapetum in liverworts, giving support to the notion that the presence of a tapetum is an ancestral land plant trait. Finally, as early land plant spore walls exhibit evidence of tapetal deposition,8,9,10,11,12 a tapetal capsule wall could have provided these plants with a developmental mechanism for sporopollenin deposition.

The Polycomb repressive complex 2 deposits H3K27me3 and represses transposable elements in a broad range of eukaryotes.

The mobility of transposable elements (TEs) contributes to evolution of genomes. Their uncontrolled activity causes genomic instability; therefore, expression of TEs is silenced by host genomes. TEs are marked with DNA and H3K9 methylation, which are associated with silencing in flowering plants, animals, and fungi. However, in distantly related groups of eukaryotes, TEs are marked by H3K27me3 deposited by the Polycomb repressive complex 2 (PRC2), an epigenetic mark associated with gene silencing in flowering plants and animals. The direct silencing of TEs by PRC2 has so far only been shown in one species of ciliates. To test if PRC2 silences TEs in a broader range of eukaryotes, we generated mutants with reduced PRC2 activity and analyzed the role of PRC2 in extant species along the lineage of Archaeplastida and in the diatom P. tricornutum. In this diatom and the red alga C. merolae, a greater proportion of TEs than genes were repressed by PRC2, whereas a greater proportion of genes than TEs were repressed by PRC2 in bryophytes. In flowering plants, TEs contained potential cis-elements recognized by transcription factors and associated with neighbor genes as transcriptional units repressed by PRC2. Thus, silencing of TEs by PRC2 is observed not only in Archaeplastida but also in diatoms and ciliates, suggesting that PRC2 deposited H3K27me3 to silence TEs in the last common ancestor of eukaryotes. We hypothesize that during the evolution of Archaeplastida, TE fragments marked with H3K27me3 were selected to shape transcriptional regulation, controlling networks of genes regulated by PRC2.

Gamete expression of TALE class HD genes activates the diploid sporophyte program in Marchantia polymorpha.

Eukaryotic life cycles alternate between haploid and diploid phases and in phylogenetically diverse unicellular eukaryotes, expression of paralogous homeodomain genes in gametes primes the haploid-to-diploid transition. In the unicellular chlorophyte alga Chlamydomonas, KNOX and BELL TALE-homeodomain genes mediate this transition. We demonstrate that in the liverwort Marchantia polymorpha, paternal (sperm) expression of three of five phylogenetically diverse BELL genes, MpBELL234, and maternal (egg) expression of both MpKNOX1 and MpBELL34 mediate the haploid-to-diploid transition. Loss-of-function alleles of MpKNOX1 result in zygotic arrest, whereas a loss of either maternal or paternal MpBELL234 results in variable zygotic and early embryonic arrest. Expression of MpKNOX1 and MpBELL34 during diploid sporophyte development is consistent with a later role for these genes in patterning the sporophyte. These results indicate that the ancestral mechanism to activate diploid gene expression was retained in early diverging land plants and subsequently co-opted during evolution of the diploid sporophyte body.

Oil Body Formation in Marchantia polymorpha Is Controlled by MpC1HDZ and Serves as a Defense against Arthropod Herbivores.

The origin of a terrestrial flora in the Ordovician required adaptation to novel biotic and abiotic stressors. Oil bodies, a synapomorphy of liverworts, accumulate secondary metabolites, but their function and development are poorly understood. Oil bodies of Marchantia polymorpha develop within specialized cells as one single large organelle. Here, we show that a class I homeodomain leucine-zipper (C1HDZ) transcription factor controls the differentiation of oil body cells in two different ecotypes of the liverwort M. polymorpha, a model genetic system for early divergent land plants. In flowering plants, these transcription factors primarily modulate responses to abiotic stress, including drought. However, loss-of-function alleles of the single ortholog gene, MpC1HDZ, in M. polymorpha did not exhibit phenotypes associated with abiotic stress. Rather, Mpc1hdz mutant plants were more susceptible to herbivory, and total plant extracts of the mutant exhibited reduced antibacterial activity. Transcriptomic analysis of the mutant revealed a reduction in expression of genes related to secondary metabolism that was accompanied by a specific depletion of oil body terpenoid compounds. Through time-lapse imaging, we observed that MpC1HDZ expression maxima precede oil body formation, indicating that MpC1HDZ mediates differentiation of oil body cells. Our results indicate that M. polymorpha oil bodies, and MpC1HDZ, are critical for defense against herbivory, but not for abiotic stress tolerance. Thus, C1HDZ genes were co-opted to regulate separate responses to biotic and abiotic stressors in two distinct land plant lineages.

A bicontinental origin of polyploid Australian/New Zealand Lepidium species (Brassicaceae)? Evidence from genomic in situ hybridization.

BACKGROUND AND AIMS: Incongruence between chloroplast and nuclear DNA phylogenies, and single additive nucleotide positions in internal transcribed spacer (ITS) sequences of polyploid Australian/New Zealand (NZ) Lepidium species have been used to suggest a bicontinental hybrid origin. This pattern was explained by two trans-oceanic dispersals of Lepidium species from California and Africa and subsequent hybridization followed by homogenization of the ribosomal DNA sequence either to the Californian (C-clade) or to the African ITS-type (A-clade) in two different ITS-lineages of Australian/NZ Lepidium polyploids. METHODS: Genomic in situ hybridization (GISH) was used to unravel the genomic origin of polyploid Australian/NZ Lepidium species. Fluorescence in situ hybridization (FISH) with ribosomal DNA (rDNA) probes was applied to test the purported ITS evolution, and to facilitate chromosome counting in high-numbered polyploids. KEY RESULTS: In Australian/NZ A-clade Lepidium polyploids, GISH identified African and Australian/NZ C-clade species as putative ancestral genomes. Neither the African nor the Californian genome were detected in Australian/NZ C-clade species and the Californian genome was not detected in Australian/NZ A-clade species. Five of the eight polyploid species (from 7x to 11x) displayed a diploid-like set of rDNA loci. Even the undecaploid species Lepidium muelleriferdinandi (2n = 11x = 88) showed only one pair of each rDNA repeat. In A-clade allopolyploids, in situ rDNA localization combined with GISH corroborated the presence of the African ITS-type. CONCLUSIONS: The nuclear genomes of African and Australian/NZ C-clade species were detected by GISH in allopolyploid Australian/NZ Lepidium species of the A-clade, supporting their hybrid origin. The presumed hybrid origin of Australian/NZ C-clade taxa could not be confirmed. Hence, it is assumed that Californian ancestral taxa experienced rapid radiation in Australia/NZ into extant C-clade polyploid taxa followed by hybridization with African species. As a result, A-clade allopolyploid Lepidium species share the Californian chloroplast type and the African ITS-type with the C-clade Australian/NZ polyploid and African diploid species, respectively.