Identification of reference genes for real-time quantitative PCR experiments in the liverwort Marchantia polymorpha

Analysis of auxin responses in the fern Ceratopteris richardii identifies the developmental phase as a major determinant for response properties

The auxin signaling molecule regulates a range of plant growth and developmental processes. The core transcriptional machinery responsible for auxin-mediated responses is conserved across all land plants. Genetic, physiological and molecular exploration in bryophyte and angiosperm model species have shown both qualitative and quantitative differences in auxin responses. Given the highly divergent ontogeny of the dominant gametophyte (bryophytes) and sporophyte (angiosperms) generations, however, it is unclear whether such differences derive from distinct phylogeny or ontogeny. Here, we address this question by comparing a range of physiological, developmental and molecular responses to auxin in both generations of the model fern Ceratopteris richardii. We find that auxin response in Ceratopteris gametophytes closely resembles that of a thalloid bryophyte, whereas the sporophyte mimics auxin response in flowering plants. This resemblance manifests both at the phenotypic and transcriptional levels. Furthermore, we show that disrupting auxin transport can lead to ectopic sporophyte induction on the gametophyte, suggesting a role for auxin in the alternation of generations. Our study thus identifies developmental phase, rather than phylogeny, as a major determinant of auxin response properties in land plants.

Streamlined regulation of chloroplast development in the liverwort Marchantia polymorpha

Chloroplasts develop from undifferentiated plastids in response to light. In angiosperms, after the perception of light, the Elongated Hypocotyl 5 (HY5) transcription factor initiates photomorphogenesis, and two families of transcription factors known as GOLDEN2-LIKE (GLK) and GATA are considered master regulators of chloroplast development. In addition, the MIR171-targeted SCARECROW-LIKE GRAS transcription factors also impact chlorophyll biosynthesis. The extent to which these proteins carry out conserved roles in non-seed plants is not known. Using the model liverwort Marchantia polymorpha, we show that GLK controls chloroplast biogenesis, and HY5 shows a small conditional effect on chlorophyll content. Chromatin immunoprecipitation sequencing (ChIP-seq) revealed that MpGLK has a broader set of targets than has been reported in angiosperms. We also identified a functional GLK homolog in green algae. In summary, our data support the hypothesis that GLK carries out a conserved role relating to chloroplast biogenesis in land plants and green algae.

MYB-related transcription factors control chloroplast biogenesis

Chloroplast biogenesis is dependent on master regulators from the GOLDEN2-LIKE (GLK) family of transcription factors. However, glk mutants contain residual chlorophyll, indicating that other proteins must be involved. Here, we identify MYB-related transcription factors as regulators of chloroplast biogenesis in the liverwort Marchantia polymorpha and angiosperm Arabidopsis thaliana. In both species, double-mutant alleles in MYB-related genes show very limited chloroplast development, and photosynthesis gene expression is perturbed to a greater extent than in GLK mutants. Genes encoding enzymes of chlorophyll biosynthesis are controlled by MYB-related and GLK proteins, whereas those allowing CO2 fixation, photorespiration, and photosystem assembly and repair require MYB-related proteins. Regulation between the MYB-related and GLK transcription factors appears more extensive in A. thaliana than in M. polymorpha. Thus, MYB-related and GLK genes have overlapping as well as distinct targets. We conclude that MYB-related and GLK transcription factors orchestrate chloroplast development in land plants.

Elucidation of arsenic detoxification mechanism in Marchantia polymorpha: The role of ACR3

The arsenic-specific ACR3 transporter plays pivotal roles in As detoxification in yeast and a group of ancient tracheophytes, the ferns. Despite putative ACR3 genes being present in the genomes of bryophytes, whether they have the same relevance also in this lineage is currently unknown. In this study, we characterized the MpACR3 gene from the bryophyte Marchantia polymorpha L. through a multiplicity of functional approaches ranging from phylogenetic reconstruction, expression analysis, loss- and gain-of-function as well as genetic complementation with an MpACR3 gene tagged with a fluorescent protein. Genetic complementation demonstrates that MpACR3 plays a pivotal role in As tolerance in M. polymorpha, with loss-of-function Mpacr3 mutants being hypersensitive and MpACR3 overexpressors more tolerant to As. Additionally, MpACR3 activity regulates intracellular As concentration, affects its speciation and controls the levels of intracellular oxidative stress. The MpACR3::3xCitrine appears to localize at the plasma membrane and possibly in other endomembrane systems. Taken together, these results demonstrate the pivotal function of ACR3 detoxification in both sister lineages of land plants, indicating that it was present in the common ancestor to all embryophytes. We propose that Mpacr3 mutants could be used in developing countries as low-cost and low-technology visual bioindicators to detect As pollution in water.

Thermospermine Is an Evolutionarily Ancestral Phytohormone Required for Organ Development and Stress Responses in Marchantia Polymorpha

Thermospermine suppresses auxin-inducible xylem differentiation, whereas its structural isomer, spermine, is involved in stress responses in angiosperms. The thermospermine synthase, ACAULIS5 (ACL5), is conserved from algae to land plants, but its physiological functions remain elusive in non-vascular plants. Here, we focused on MpACL5, a gene in the liverwort Marchantia polymorpha, that rescued the dwarf phenotype of the acl5 mutant in Arabidopsis. In the Mpacl5 mutants generated by genome editing, severe growth retardation was observed in the vegetative organ, thallus, and the sexual reproductive organ, gametangiophore. The mutant gametangiophores exhibited remarkable morphological defects such as short stalks, fasciation and indeterminate growth. Two gametangiophores fused together, and new gametangiophores were often initiated from the old ones. Furthermore, Mpacl5 showed altered responses to heat and salt stresses. Given the absence of spermine in bryophytes, these results suggest that thermospermine has a dual primordial function in organ development and stress responses in M. polymorpha. The stress response function may have eventually been assigned to spermine during land plant evolution.

Divergent evolution of the alcohol-forming pathway of wax biosynthesis among bryophytes

The plant cuticle is a hydrophobic barrier, which seals the epidermal surface of most aboveground organs. While the cuticle biosynthesis of angiosperms has been intensively studied, knowledge about its existence and composition in nonvascular plants is scarce. Here, we identified and characterized homologs of Arabidopsis thaliana fatty acyl-CoA reductase (FAR) ECERIFERUM 4 (AtCER4) and bifunctional wax ester synthase/acyl-CoA:diacylglycerol acyltransferase 1 (AtWSD1) in the liverwort Marchantia polymorpha (MpFAR2 and MpWSD1) and the moss Physcomitrium patens (PpFAR2A, PpFAR2B, and PpWSD1). Although bryophyte harbor similar compound classes as described for angiosperm cuticles, their biosynthesis may not be fully conserved between the bryophytes M. polymorpha and P. patens or between these bryophytes and angiosperms. While PpFAR2A and PpFAR2B contribute to the production of primary alcohols in P. patens, loss of MpFAR2 function does not affect the wax profile of M. polymorpha. By contrast, MpWSD1 acts as the major wax ester-producing enzyme in M. polymorpha, whereas mutations of PpWSD1 do not affect the wax ester levels of P. patens. Our results suggest that the biosynthetic enzymes involved in primary alcohol and wax ester formation in land plants have either evolved multiple times independently or undergone pronounced radiation followed by the formation of lineage-specific toolkits.

MpTGA, together with MpNPR, regulates sexual reproduction and independently affects oil body formation in Marchantia polymorpha

In angiosperms, basic leucine-zipper (bZIP) TGACG-motif-binding (TGA) transcription factors (TFs) regulate developmental and stress-related processes, the latter often involving NON EXPRESSOR OF PATHOGENESIS-RELATED GENES (NPR) coregulator interactions. To gain insight into their functions in an early diverging land-plant lineage, the single MpTGA and sole MpNPR genes were investigated in the liverwort Marchantia polymorpha. We generated Marchantia MpTGA and MpNPR knockout and overexpression mutants and conducted morphological, transcriptomic and expression studies. Furthermore, we investigated MpTGA interactions with wild-type and mutagenized MpNPR and expanded our analyses including TGA TFs from two streptophyte algae. Mptga mutants fail to induce the switch from vegetative to reproductive development and lack gametangiophore formation. MpTGA and MpNPR proteins interact and Mpnpr mutant analysis reveals a novel coregulatory NPR role in sexual reproduction. Additionally, MpTGA acts independently of MpNPR as a repressor of oil body (OB) formation and can thereby affect herbivory. The single MpTGA TF exerts a dual role in sexual reproduction and OB formation in Marchantia. Common activities of MpTGA/MpNPR in sexual development suggest that coregulatory interactions were established after emergence of land-plant-specific NPR genes and contributed to the diversification of TGA TF functions during land-plant evolution.

MiRNAs differentially expressed in vegetative and reproductive organs of Marchantia polymorpha - insights into their expression pattern, gene structures and function

MicroRNAs regulate gene expression affecting a variety of plant developmental processes. The evolutionary position of Marchantia polymorpha makes it a significant model to understand miRNA-mediated gene regulatory pathways in plants. Previous studies focused on conserved miRNA-target mRNA modules showed their critical role in Marchantia development. Here, we demonstrate that the differential expression of conserved miRNAs among land plants and their targets in selected organs of Marchantia additionally underlines their role in regulating fundamental developmental processes. The main aim of this study was to characterize selected liverwort-specific miRNAs, as there is a limited knowledge on their biogenesis, accumulation, targets, and function in Marchantia. We demonstrate their differential accumulation in vegetative and generative organs. We reveal that all liverwort-specific miRNAs examined are encoded by independent transcriptional units. MpmiR11737a, MpmiR11887 and MpmiR11796, annotated as being encoded within protein-encoding genes, have their own independent transcription start sites. The analysis of selected liverwort-specific miRNAs and their pri-miRNAs often reveal correlation in their levels, suggesting transcriptional regulation. However, MpmiR11796 shows a reverse correlation to its pri-miRNA level, suggesting post-transcriptional regulation. Moreover, we identify novel targets for selected liverwort-specific miRNAs and demonstrate an inverse correlation between their expression and miRNA accumulation. In the case of one miRNA precursor, we provide evidence that it encodes two functional miRNAs with two independent targets. Overall, our research sheds light on liverwort-specific miRNA gene structure, provides new data on their biogenesis and expression regulation. Furthermore, identifying their targets, we hypothesize the potential role of these miRNAs in early land plant development and functioning.

Control of stem cell behavior by CLE-JINGASA signaling in the shoot apical meristem in Marchantia polymorpha

Land plants undergo indeterminate growth by the activity of meristems in both gametophyte (haploid) and sporophyte (diploid) generations. In the sporophyte of the flowering plant Arabidopsis thaliana, the apical meristems are located at the shoot and root tips in which a number of regulatory gene homologs are shared for their development, implying deep evolutionary origins. However, little is known about their functional conservation with gametophytic meristems in distantly related land plants such as bryophytes, even though genomic studies have revealed that the subfamily-level diversity of regulatory genes is mostly conserved throughout land plants. Here, we show that a NAM/ATAF/CUC (NAC) domain transcription factor, JINGASA (MpJIN), acts downstream of CLAVATA3 (CLV3)/ESR-related (CLE) peptide signaling and controls stem cell behavior in the gametophytic shoot apical meristem of the liverwort Marchantia polymorpha. In the meristem, strong MpJIN expression was associated with the periclinal cell division at the periphery of the stem cell zone (SCZ), whereas faint MpJIN expression was found at the center of the SCZ. Time course observation indicates that the MpJIN-negative cells are lost from the SCZ and respecified de novo at two separate positions during the dichotomous branching event. Consistently, the induction of MpJIN results in ectopic periclinal cell division in the SCZ and meristem termination. Based on the comparative expression data, we speculate that the function of JIN/FEZ subfamily genes was shared among the shoot apical meristems in the gametophyte and sporophyte generations in early land plants but was lost in certain lineages, including the flowering plant A. thaliana.

Phosphatidic acid phosphohydrolase modulates glycerolipid synthesis in Marchantia polymorpha and is crucial for growth under both nutrient-replete and -deficient conditions

MAIN CONCLUSION

The phosphatidic acid phosphohydrolase of Marchantia polymorpha modulates plastid glycolipid synthesis through the ER pathway and is essential for normal plant development regardless of nutrient availability. Membrane lipid remodeling is one of the strategies plant cells use to secure inorganic phosphate (Pi) for plant growth, but many aspects of the molecular mechanism and its regulation remain unclear. Here we analyzed membrane lipid remodeling using a non-vascular plant, Marchantia polymorpha. The lipid composition and fatty acid profile during Pi starvation in M. polymorpha revealed a decrease in phospholipids and an increase in both galactolipids and betaine lipids. In Arabidopsis thaliana, phosphatidic acid phosphohydrolase (PAH) is involved in phospholipid degradation and is crucial for tolerance to both Pi and nitrogen starvation. We produced two M. polymorpha PAH (MpPAH) knockout mutants (Mppah-1 and Mppah-2) and found that, unlike Arabidopsis mutants, Mppah impaired plant growth with shorter rhizoids compared with wild-type plants even under nutrient-replete conditions. Mutation of MpPAH did not significantly affect the mole percent of each glycerolipid among total membrane glycerolipids from whole plants under both Pi-replete and Pi-deficient conditions. However, the fatty acid composition of monogalactosyldiacylglycerol indicated that the amount of plastid glycolipids produced through the endoplasmic reticulum pathway was suppressed in Mppah mutants. Phospholipids accumulated in the mutants under N starvation. These results reveal that MpPAH modulates plastid glycolipid synthesis through the endoplasmic reticulum pathway more so than what has been observed for Arabidopsis PAH; moreover, unlike Arabidopsis, MpPAH is crucial for M. polymorpha growth regardless of nutrient availability.

Insights into cryptochrome modulation of ABA signaling to mediate dormancy regulation in Marchantia polymorpha

The acquisition of dormancy capabilities has enabled plants to survive in adverse terrestrial environmental conditions. Dormancy accumulation and release is coupled with light signaling, which is well studied in Arabidopsis, but it is unclear in the distant nonvascular relative. We study the characteristics and function on dormancy regulation of a blue light receptor cryptochrome in Marchantia polymorpha (MpCRY). Here, we identified MpCRY via bioinformatics and mutant complement analysis. The biochemical characteristics were assessed by multiple protein-binding assays. The function of MpCRY in gemma dormancy was clarified by overexpression and mutation of MpCRY, and its mechanism was analyzed via RNA sequencing and quantitative PCR analyses associated with hormone treatment. We found that the unique MpCRY protein in M. polymorpha undergoes both blue light-promoted interaction with itself (self-interaction) and blue light-dependent phosphorylation. MpCRY has the specific characteristics of blue light-induced nuclear localization and degradation. We further demonstrated that MpCRY transcriptionally represses abscisic acid (ABA) signaling-related gene expression to suppress gemma dormancy, which is dependent on blue light signaling. Our findings indicate that MpCRY possesses specific biochemical and molecular characteristics, and modulates ABA signaling under blue light conditions to regulate gemma dormancy in M. polymorpha.

Transcriptome changes in chlorsulfuron-treated plants are caused by acetolactate synthase inhibition and not induction of a herbicide detoxification system in Marchantia polymorpha

A sensing mechanism in mammals perceives xenobiotics and induces the transcription of genes encoding proteins that detoxify these molecules. However, it is unclear if plants sense xenobiotics, and activate an analogous signalling system leading to their detoxification. Using the liverwort Marchantia polymorpha, we tested the hypothesis that there is a sensing system in plants that perceives herbicides resulting in the increased transcription of genes encoding proteins that detoxify these herbicides. Consistent with the hypothesis, we show that chlorsulfuron-treatment induces changes in the M. polymorpha transcriptome. However, these transcriptome changes do not occur in chlorsulfuron (CS)-treated target site resistant mutants, where the gene encoding the target carries a mutation that confers resistance to chlorsulfuron. Instead, we show that inactivation of the chlorsulfuron target, acetolactate synthase (ALS) (also known as acetohydroxyacid synthase (AHAS)), is required for the transcriptome response. These data demonstrate that the transcriptome changes in chlorsulfuron-treated plants are caused by disrupted amino acid synthesis and metabolism resulting from acetolactate synthase inhibition, and indicate that the transcriptome changes are not caused by a herbicide sensing mechanism.

Plant lineage-specific PIKMIN1 drives APC/CCCS52A2 E3-ligase activity-dependent cell division

The anaphase-promoting complex/cyclosome (APC/C) marks key cell cycle proteins for proteasomal breakdown, thereby ensuring unidirectional progression through the cell cycle. Its target recognition is temporally regulated by activating subunits, one of which is called CELL CYCLE SWITCH 52 A2 (CCS52A2). We sought to expand the knowledge on the APC/C by using the severe growth phenotypes of CCS52A2-deficient Arabidopsis (Arabidopsis thaliana) plants as a readout in a suppressor mutagenesis screen, resulting in the identification of the previously undescribed gene called PIKMIN1 (PKN1). PKN1 deficiency rescues the disorganized root stem cell phenotype of the ccs52a2-1 mutant, whereas an excess of PKN1 inhibits the growth of ccs52a2-1 plants, indicating the need for control of PKN1 abundance for proper development. Accordingly, the lack of PKN1 in a wild-type background negatively impacts cell division, while its systemic overexpression promotes proliferation. PKN1 shows a cell cycle phase-dependent accumulation pattern, localizing to microtubular structures, including the preprophase band, the mitotic spindle, and the phragmoplast. PKN1 is conserved throughout the plant kingdom, with its function in cell division being evolutionarily conserved in the liverwort Marchantia polymorpha. Our data thus demonstrate that PKN1 represents a novel, plant-specific protein with a role in cell division that is likely proteolytically controlled by the CCS52A2-activated APC/C.

Control of vegetative reproduction in Marchantiapolymorpha by the KAI2-ligand signaling pathway

In vegetative reproduction of Marchantia polymorpha (M. polymorpha), propagules, called gemmae, are formed in gemma cups. Despite its significance for survival, control of gemma and gemma cup formation by environmental cues is not well understood. We show here that the number of gemmae formed in a gemma cup is a genetic trait. Gemma formation starts from the central region of the floor of the gemma cup, proceeds to the periphery, and terminates when the appropriate number of gemmae is initiated. The MpKARRIKIN INSENSITIVE2 (MpKAI2)-dependent signaling pathway promotes gemma cup formation and gemma initiation. The number of gemmae in a cup is controlled by modulating the ON/OFF switch of the KAI2-dependent signaling. Termination of the signaling results in the accumulation of MpSMXL, a suppressor protein. In the Mpsmxl mutants, gemma initiation continues, leading to the formation of a highly increased number of gemmae in a cup. Consistent with its function, the MpKAI2-dependent signaling pathway is active in gemma cups where gemmae initiate, as well as in the notch region of the mature gemma and midrib of the ventral side of the thallus. In this work, we also show that GEMMA CUP-ASSOCIATED MYB1 works downstream of this signaling pathway to promote gemma cup formation and gemma initiation. We also found that the availability of potassium affects gemma cup formation independently from the KAI2-dependent signaling pathway in M. polymorpha. We propose that the KAI2-dependent signaling pathway functions to optimize vegetative reproduction by adapting to the environment in M. polymorpha.

A Pipeline to Monitor Proteasome Homeostasis in Plants

The proteasome is a key component for regulation of protein turnover across kingdoms. The proteasome has been shown to be involved in or affected by various stress conditions in multiple model organisms in plants. As such, studying proteasome homeostasis is crucial to understand its participation in different cellular conditions. However, the involvement of the proteasome in many cellular processes and its interplay with other degradation pathways hamper the interpretation of experiments based on a single approach. Thus, it is crucial to formulate a framework to investigate proteasome dynamics in different model organisms including plants. Here, we describe a pipeline to monitor proteasome homeostasis using four different methods including (i) luminescent-based proteasome activity measurement, (ii) immunoblot analysis of ubiquitinated proteins, (iii) evaluation of proteasome subunit protein levels, and (iv) monitoring of the proteasome stress regulon on mRNA levels using quantitative real-time PCR (polymerase chain reaction).

Functional analysis of the BEige and Chediak-Higashi domain gene MpSPIRRIG in Marchantia polymorpha

BEige and Chediak-Higashi domain containing proteins (BDCPs) have been described to function in membrane-dependent processes in eukaryotes. This role was also observed for the BDCP SPIRRIG (SPI) in the model plant Arabidopsis thaliana in the context of cell morphogenesis. Additionally, AtSPI was found to control salt stress resistance by mediating mRNA stability and salt stress-dependent processing body formation. In this work, we utilize an evolutionarily comparative approach to unravel conserved, basal BDCP functions in the liverwort Marchantia polymorpha. Our phenotypic and physiological analyses show that MpSPI is involved in cell morphogenesis and salt resistance regulation, indicating that both functions are evolutionarily conserved between the two species. Co-localization was found with endosomal and P-body markers, suggesting links to membrane-dependent processes and mRNA metabolism. Finally, we present transcriptomics data showing that AtSPI and MpSPI regulate orthologous genes in A. thaliana and M. polymorpha.

The wound-activated ERF15 transcription factor drives Marchantia polymorpha regeneration by activating an oxylipin biosynthesis feedback loop

The regenerative potential in response to wounding varies widely among species. Within the plant lineage, the liverwort Marchantia polymorpha displays an extraordinary regeneration capacity. However, its molecular pathways controlling the initial regeneration response are unknown. Here, we demonstrate that the MpERF15 transcription factor gene is instantly activated after wounding and is essential for gemmaling regeneration following tissue incision. MpERF15 operates both upstream and downstream of the MpCOI1 oxylipin receptor by controlling the expression of oxylipin biosynthesis genes. The resulting rise in the oxylipin dinor-12-oxo-phytodienoic acid (dn-OPDA) levels results in an increase in gemma cell number and apical notch organogenesis, generating highly disorganized and compact thalli. Our data pinpoint MpERF15 as a key factor activating an oxylipin biosynthesis amplification loop after wounding, which eventually results in reactivation of cell division and regeneration. We suggest that the genetic networks controlling oxylipin biosynthesis in response to wounding might have been reshuffled over evolution.

PIF-independent regulation of growth by an evening complex in the liverwort Marchantia polymorpha

Previous studies in the liverwort Marchantia polymorpha have shown that the putative evening complex (EC) genes LUX ARRHYTHMO (LUX) and ELF4-LIKE (EFL) have a function in the liverwort circadian clock. Here, we studied the growth phenotypes of MpLUX and MpEFL loss-of-function mutants, to establish if PHYTOCHROME-INTERACTING FACTOR (PIF) and auxin act downstream of the M. polymorpha EC in a growth-related pathway similar to the one described for the flowering plant Arabidopsis. We examined growth rates and cell properties of loss-of-function mutants, analyzed protein-protein interactions and performed gene expression studies using reporter genes. Obtained data indicate that an EC can form in M. polymorpha and that this EC regulates growth of the thallus. Altered auxin levels in Mplux mutants could explain some of the phenotypes related to an increased thallus surface area. However, because MpPIF is not regulated by the EC, and because Mppif mutants do not show reduced growth, the growth phenotype of EC-mutants is likely not mediated via MpPIF. In Arabidopsis, the circadian clock regulates elongation growth via PIF and auxin, but this is likely not an evolutionarily conserved growth mechanism in land plants. Previous inventories of orthologs to Arabidopsis clock genes in various plant lineages showed that there is high levels of structural differences between clocks of different plant lineages. Here, we conclude that there is also variation in the output pathways used by the different plant clocks to control growth and development.

Distinct Functions of the Atypical Terminal Hydrophilic Domain of the HKT Transporter in the Liverwort Marchantia polymorpha

K+/Na+ homeostasis is important for land plants, particularly under salt stress. In this study, the structure and ion transport properties of the high-affinity K+ transporter (HKT) of the liverwort Marchantia polymorpha were investigated. Only one HKT gene, MpHKT1, was identified in the genome of M. polymorpha. Phylogenetic analysis of HKT proteins revealed that non-seed plants possess HKTs grouped into a clade independent of the other two clades including HKTs of angiosperms. A distinct long hydrophilic domain was found in the C-terminus of MpHKT1. Complementary DNA (cDNA) of truncated MpHKT1 (t-MpHKT1) encoding the MpHKT_Δ596-812 protein was used to examine the functions of the C-terminal domain. Both MpHKT1 transporters fused with enhanced green fluorescent protein at the N-terminus were localized to the plasma membrane when expressed in rice protoplasts. Two-electrode voltage clamp experiments using Xenopus laevis oocytes indicated that MpHKT1 mediated the transport of monovalent alkali cations with higher selectivity for Na+ and K+, but truncation of the C-terminal domain significantly reduced the transport activity with a decrease in the Na+ permeability. Overexpression of MpHKT1 or t-MpHKT1 in M. polymorpha conferred accumulation of higher Na+ levels and showed higher Na+ uptake rates, compared to those of wild-type plants; however, phenotypes with t-MpHKT1 were consistently weaker than those with MpHKT1. Together, these findings suggest that the hydrophilic C-terminal domain plays a unique role in the regulation of transport activity and ion selectivity of MpHKT1.

Normal oil body formation in Marchantia polymorpha requires functional coat protein complex I proteins

Eukaryotic cells possess endomembrane organelles equipped with specific sets of proteins, lipids, and polysaccharides that are fundamental for realizing each organelle's specific function and shape. A tightly regulated membrane trafficking system mediates the transportation and localization of these substances. Generally, the secretory/exocytic pathway is responsible for transporting cargo to the plasma membrane and/or the extracellular space. However, in the case of oil body cells in the liverwort Marchantia polymorpha, the oil body, a liverwort-unique organelle, is thought to be formed by secretory vesicle fusion through redirection of the secretory pathway inside the cell. Although their formation mechanism remains largely unclear, oil bodies exhibit a complex and bumpy surface structure. In this study, we isolated a mutant with spherical oil bodies through visual screening of mutants with abnormally shaped oil bodies. This mutant harbored a mutation in a coat protein complex I (COPI) subunit MpSEC28, and a similar effect on oil body morphology was also detected in knockdown mutants of other COPI subunits. Fluorescently tagged MpSEC28 was localized to the periphery of the Golgi apparatus together with other subunits, suggesting that it is involved in retrograde transport from and/or in the Golgi apparatus as a component of the COPI coat. The Mpsec28 mutants also exhibited weakened stiffness of the thalli, suggesting impaired cell-cell adhesion and cell wall integrity. These findings suggest that the mechanism of cell wall biosynthesis is also involved in shaping the oil body in M. polymorpha, supporting the redirection of the secretory pathway inward the cell during oil body formation.

Roles of a Cryptochrome in Carbon Fixation and Sucrose Metabolism in the Liverwort Marchantia polymorpha

In vascular plants, cryptochromes acting as blue-light photoreceptors have various functions to adapt plants to the fluctuating light conditions on land, while the roles of cryptochromes in bryophytes have been rarely reported. In this study, we investigated functions of a single-copy ortholog of cryptochrome (MpCRY) in the liverwort Marchantia polymorpha. Knock-out of MpCRY showed that a large number of the mutant plants exhibited asymmetric growth of thalli under blue light. Transcriptome analyses indicated that MpCRY is mainly involved in photosynthesis and sugar metabolism. Further physiological analysis showed that Mpcry mutant exhibited a reduction in CO₂ uptake and sucrose metabolism. In addition, exogenous application of sucrose or glucose partially restored the symmetrical growth of the Mpcry mutant thalli. Together, these results suggest that MpCRY is involved in the symmetrical growth of thallus and the regulation of carbon fixation and sucrose metabolism in M. polymorpha.

The BES1/BZR1-family transcription factor MpBES1 regulates cell division and differentiation in Marchantia polymorpha

Brassinosteroids (BRs) play essential roles in growth and development in seed plants;¹ disturbances in BR homeostasis lead to altered mitotic activity in meristems^2,3 and organ boundaries^4,5 and to changes in meristem determinacy.⁶ An intricate signaling cascade linking the perception of BRs at the plasma membrane to the regulation of master transcriptional regulators belonging to the BEH, for BES1 homologues, family⁷ has been described in great detail in model angiosperms. Homologs of these transcription factors are present in streptophyte algae and in land plant lineages where BR signaling or function is absent or has not yet been characterized. The genome of the bryophyte Marchantia polymorpha does not encode for BR receptors but includes one close ortholog of Arabidopsis thaliana BRI1-EMS-SUPPRESSOR 1 (AtBES1)⁸ and Arabidopsis thaliana BRASSINAZOLE-RESISTANT 1 (AtBZR1),⁹ MpBES1. Altered levels of MpBES1 severely compromised cell division and differentiation, resulting in stunted thalli that failed to differentiate adult tissues and reproductive organs. The transcriptome of Mpbes1 knockout plants revealed a significant overlap with homologous functions controlled by AtBES1 and AtBZR1, suggesting that members of this gene family share a subset of common targets. Indeed, MpBES1 behaved as a gain-of-function substitute of AtBES1/AtBZR1 when expressed in Arabidopsis, probably because it mediates conserved functions but evades the regulatory mechanisms that native counterparts are subject to. Our results show that this family of transcription factors plays an ancestral role in the control of cell division and differentiation in plants and that BR signaling likely co-opted this function and imposed additional regulatory checkpoints upon it.

DE-ETIOLATED1 has a role in the circadian clock of the liverwort Marchantia polymorpha

Previous studies of plant circadian clock evolution have often relied on clock models and genes defined in Arabidopsis. These studies identified homologues with seemingly conserved function, as well as frequent gene loss. In the present study, we aimed to identify candidate clock genes in the liverwort Marchantia polymorpha using a more unbiased approach. To identify genes with circadian rhythm we sequenced the transcriptomes of gemmalings in a time series in constant light conditions. Subsequently, we performed functional studies using loss-of-function mutants and gene expression reporters. Among the genes displaying circadian rhythm, a homologue to the transcriptional co-repressor Arabidopsis DE-ETIOLATED1 showed high amplitude and morning phase. Because AtDET1 is arrhythmic and associated with the morning gene function of AtCCA1/LHY, that lack a homologue in liverworts, we functionally studied DET1 in M. polymorpha. We found that the circadian rhythm of MpDET1 expression is disrupted in loss-of-function mutants of core clock genes and putative evening-complex genes. MpDET1 knock-down in turn results in altered circadian rhythm of nyctinastic thallus movement and clock gene expression. We could not detect any effect of MpDET1 knock-down on circadian response to light, suggesting that MpDET1 has a yet unknown function in the M. polymorpha circadian clock.

Distinct Functions of Ethylene and ACC in the Basal Land Plant Marchantia polymorpha

Ethylene is a gaseous phytohormone involved in various physiological processes, including fruit ripening, senescence, root hair development and stress responses. Recent genomics studies have suggested that most homologous genes of ethylene biosynthesis and signaling are conserved from algae to angiosperms, whereas the function and biosynthesis of ethylene remain unknown in basal plants. Here, we examined the physiological effects of ethylene, an ethylene precursor, 1-aminocyclopropane-1-carboxylic acid (ACC) and an inhibitor of ethylene perception, silver thiosulfate (STS), in a basal land plant, Marchantia polymorpha. M. polymorpha plants biosynthesized ethylene, and treatment with high concentrations of ACC slightly promoted ethylene production. ACC remarkably suppressed the growth of thalli (vegetative organs) and rhizoids (root-hair-like cells), whereas exogenous ethylene slightly promoted thallus growth. STS suppressed thallus growth and induced ectopic rhizoid formation on the dorsal surface of thalli. Thus, ACC and ethylene have different effects on the vegetative growth of M. polymorpha. We generated single and double mutants of ACC synthase-like (ACSL) genes, MpACSL1 and MpACSL2. The mutants did not show obvious defects in thallus growth, ACC content and ethylene production, indicating that MpACSL genes are not essential for the vegetative growth and biosynthesis of ACC and ethylene. Gene expression analysis suggested the involvement of MpACSL1 and MpACSL2 in stress responses. Collectively, our results imply ethylene-independent function of ACC and the absence of ACC-mediated ethylene biosynthesis in M. polymorpha.

The liverwort oil body is formed by redirection of the secretory pathway

Eukaryotic cells acquired novel organelles during evolution through mechanisms that remain largely obscure. The existence of the unique oil body compartment is a synapomorphy of liverworts that represents lineage-specific acquisition of this organelle during evolution, although its origin, biogenesis, and physiological function are yet unknown. We find that two paralogous syntaxin-1 homologs in the liverwort Marchantia polymorpha are distinctly targeted to forming cell plates and the oil body, suggesting that these structures share some developmental similarity. Oil body formation is regulated by an ERF/AP2-type transcription factor and loss of the oil body increases M. polymorpha herbivory. These findings highlight a common strategy for the acquisition of organelles with distinct functions in plants, via periodical redirection of the secretory pathway depending on cellular phase transition.

Transcriptional and Morpho-Physiological Responses of Marchantia polymorpha upon Phosphate Starvation

Phosphate (Pi) is a pivotal nutrient that constraints plant development and productivity in natural ecosystems. Land colonization by plants, more than 470 million years ago, evolved adaptive mechanisms to conquer Pi-scarce environments. However, little is known about the molecular basis underlying such adaptations at early branches of plant phylogeny. To shed light on how early divergent plants respond to Pi limitation, we analyzed the morpho-physiological and transcriptional dynamics of Marchantia polymorpha upon Pi starvation. Our phylogenomic analysis highlights some gene networks present since the Chlorophytes and others established in the Streptophytes (e.g., PHR1–SPX1 and STOP1–ALMT1, respectively). At the morpho-physiological level, the response is characterized by the induction of phosphatase activity, media acidification, accumulation of auronidins, reduction of internal Pi concentration, and developmental modifications of rhizoids. The transcriptional response involves the induction of MpPHR1, Pi transporters, lipid turnover enzymes, and MpMYB14, which is an essential transcription factor for auronidins biosynthesis. MpSTOP2 up-regulation correlates with expression changes in genes related to organic acid biosynthesis and transport, suggesting a preference for citrate exudation. An analysis of MpPHR1 binding sequences (P1BS) shows an enrichment of this cis regulatory element in differentially expressed genes. Our study unravels the strategies, at diverse levels of organization, exerted by M. polymorpha to cope with low Pi availability.

Ancestral function of the phytochelatin synthase C-terminal domain in inhibition of heavy metal-mediated enzyme overactivation

Phytochelatin synthases (PCSs) play essential roles in detoxification of a broad range of heavy metals in plants and other organisms. Until now, however, no PCS gene from liverworts, the earliest branch of land plants and possibly the first one to acquire a PCS with a C-terminal domain, has been characterized. In this study, we isolated and functionally characterized the first PCS gene from a liverwort, Marchantia polymorpha (MpPCS). MpPCS is constitutively expressed in all organs examined, with stronger expression in thallus midrib. The gene expression is repressed by Cd2+ and Zn2+. The ability of MpPCS to increase heavy metal resistance in yeast and to complement cad1-3 (the null mutant of the Arabidopsis ortholog AtPCS1) proves its function as the only PCS from M. polymorpha. Site-directed mutagenesis of the most conserved cysteines of the C-terminus of the enzyme further uncovered that two twin-cysteine motifs repress, to different extents, enzyme activation by heavy metal exposure. These results highlight an ancestral function of the PCS elusive C-terminus as a regulatory domain inhibiting enzyme overactivation by essential and non-essential heavy metals. The latter finding may be relevant for obtaining crops with decreased root to shoot mobility of cadmium, thus preventing its accumulation in the food chain.

Nyctinastic thallus movement in the liverwort Marchantia polymorpha is regulated by a circadian clock

The circadian clock coordinates an organism's growth, development and physiology with environmental factors. One illuminating example is the rhythmic growth of hypocotyls and cotyledons in Arabidopsis thaliana. Such daily oscillations in leaf position are often referred to as sleep movements or nyctinasty. Here, we report that plantlets of the liverwort Marchantia polymorpha show analogous rhythmic movements of thallus lobes, and that the circadian clock controls this rhythm, with auxin a likely output pathway affecting these movements. The mechanisms of this circadian clock are partly conserved as compared to angiosperms, with homologs to the core clock genes PRR, RVE and TOC1 forming a core transcriptional feedback loop also in M. polymorpha.

Systematic Tools for Reprogramming Plant Gene Expression in a Simple Model, Marchantia polymorpha

We present the OpenPlant toolkit, a set of interlinked resources and techniques to develop Marchantia as testbed for bioengineering in plants. Marchantia is a liverwort, a simple plant with an open form of development that allows direct visualization of gene expression and dynamics of cellular growth in living tissues. We describe new techniques for simple and efficient axenic propagation and maintenance of Marchantia lines with no requirement for glasshouse facilities. Marchantia plants spontaneously produce clonal propagules within a few weeks of regeneration, and lines can be amplified million-fold in a single generation by induction of the sexual phase of growth, crossing, and harvesting of progeny spores. The plant has a simple morphology and genome with reduced gene redundancy, and the dominant phase of its life cycle is haploid, making genetic analysis easier. We have built robust Loop assembly vector systems for nuclear and chloroplast transformation and genome editing. These have provided the basis for building and testing a modular library of standardized DNA elements with highly desirable properties. We have screened transcriptomic data to identify a range of candidate genes, extracted putative promoter sequences, and tested them in vivo to identify new constitutive promoter elements. The resources have been combined into a toolkit for plant bioengineering that is accessible for laboratories without access to traditional facilities for plant biology research. The toolkit is being made available under the terms of the OpenMTA and will facilitate the establishment of common standards and the use of this simple plant as testbed for synthetic biology.

MpFEW RHIZOIDS1 miRNA-Mediated Lateral Inhibition Controls Rhizoid Cell Patterning in Marchantia polymorpha

Lateral inhibition patterns differentiated cell types among equivalent cells during development in bacteria, metazoans, and plants. Tip-growing rhizoid cells develop among flat epidermal cells in the epidermis of the early-diverging land plant Marchantia polymorpha. We show that the majority of rhizoid cells develop individually, but some develop in linear, one-dimensional groups (chains) of between 2 and 7 rhizoid cells in wild-type plants. The distribution of rhizoid cells can be accounted for within a simple cellular automata model of lateral inhibition. The model predicted that in the absence of lateral inhibition, two-dimensional rhizoid cell groups (clusters) form. These can be larger than those formed with lateral inhibition. M. polymorpha rhizoid differentiation is positively regulated by the ROOT HAIR DEFECTIVE SIX-LIKE1 (MpRSL1) basic-helix-loop-helix (bHLH) transcription factor, which is directly repressed by the FEW RHIZOIDS1 (MpFRH1) microRNA (miRNA). To test if MpFRH1 miRNA acts during lateral inhibition, we generated loss-of-function (lof) mutants without the MpFRH1 miRNA. Two-dimensional clusters of rhizoids develop in Mpfrh1^lof mutants as predicted by the model for plants that lack lateral inhibition. Furthermore, two-dimensional clusters of up to 9 rhizoid cells developed in the Mpfrh1^lof mutants compared to a maximum number of 7 observed in wild-type groups. The higher steady-state levels of MpRSL1 mRNA in Mpfrh1^lof mutants indicate that MpFRH1-mediated lateral inhibition involves the repression of MpRSL1 activity. Together, the modeling and genetic data indicate that MpFRH1 miRNA mediates lateral inhibition by repressing MpRSL1 during pattern formation in the M. polymorpha epidermis.

The RopGEF KARAPPO Is Essential for the Initiation of Vegetative Reproduction in Marchantia polymorpha

Many plants can reproduce vegetatively, producing clonal progeny from vegetative cells; however, little is known about the molecular mechanisms underlying this process. Liverwort (Marchantia polymorpha), a basal land plant, propagates asexually via gemmae, which are clonal plantlets formed in gemma cups on the dorsal side of the vegetative thallus [1]. The initial stage of gemma development involves elongation and asymmetric divisions of a specific type of epidermal cell, called a gemma initial, which forms on the floor of the gemma cup [2, 3]. To investigate the regulatory mechanism underlying gemma development, we focused on two allelic mutants in which no gemma initial formed; these mutants were named karappo, meaning "empty." We used whole-genome sequencing of both mutants and molecular genetic analysis to identify the causal gene, KARAPPO (KAR), which encodes a ROP guanine nucleotide exchange factor (RopGEF) carrying a plant-specific ROP nucleotide exchanger (PRONE) catalytic domain. In vitro GEF assays showed that the full-length KAR protein and the PRONE domain have significant GEF activity toward MpROP, the only ROP GTPase in M. polymorpha. Moreover, genetic complementation experiments showed a significant role for the N- and C-terminal variable regions in gemma development. Our investigation demonstrates an essential role for KAR/RopGEF in the initiation of plantlet development from a differentiated cell, which may involve cell-polarity formation and subsequent asymmetric cell division via activation of ROP signaling, implying a similar developmental mechanism in vegetative reproduction of various land plants.

Conserved Biochemical Defenses Underpin Host Responses to Oomycete Infection in an Early-Divergent Land Plant Lineage

The expansion of plants onto land necessitated the evolution of robust defense strategies to protect against a wide array of microbial invaders. Whereas host responses to microbial colonization are extensively explored in evolutionarily young land plant lineages such as angiosperms, we know relatively little about plant-pathogen interactions in early-diverging land plants thought to better represent the ancestral state. Here, we define the transcriptional and proteomic response of the early-divergent liverwort Marchantia polymorpha to infection with the oomycete pathogen Phytophthora palmivora. We uncover a robust molecular response to oomycete colonization in Marchantia that consists of conserved land plant gene families. Direct macroevolutionary comparisons of host infection responses in Marchantia and the model angiosperm Nicotiana benthamiana further reveal a shared set of orthologous microbe-responsive genes that include members of the phenylpropanoid metabolic pathway. In addition, we identify a role for the Marchantia R2R3-MYB transcription factor MpMyb14 in activating phenylpropanoid (flavonoid) biosynthesis during oomycete infection. Mpmyb14 mutants infected with P. palmivora fail to activate phenylpropanoid biosynthesis gene expression and display enhanced disease susceptibility compared to wild-type plants. Conversely, the ectopic induction of MpMyb14 led to the accumulation of anthocyanin-like pigments and dramatically enhanced liverwort resistance to P. palmivora infection. Collectively, our results demonstrate that the Marchantia response to oomycete infection displays evolutionarily conserved features indicative of an ancestral pathogen deterrence strategy centered on phenylpropanoid-mediated biochemical defenses.

Phytochrome and Light Signaling in Marchantia

Bryophytes, which comprise liverworts, mosses, and hornworts, are one of the earliest diverging lineages of extant land plants and a key plant group for understanding evolutionary aspects of land plant adaptation. Marchantia polymorpha, a liverwort, has recently been established as a model plant species having molecular genetic tractability. In M. polymorpha, phytochrome is encoded by a single-copy gene, MpPHY, with Mpphy regulating various physiological responses through PHYTOCHROME INTERACTING FACTOR (PIF)-mediated transcriptional regulation. The phytochrome signaling system of M. polymorpha, with its single Mpphy and single PIF (MpPIF), is relatively simple compared with other model plants carrying multiple phytochromes and PIFs. Consequently, investigation of phytochrome signaling using M. polymorpha may provide novel insights into fundamental mechanisms and roles of phytochrome during the course of land plant evolution. This chapter provides a number of basic procedures, along with some tips, for designing and performing experiments with M. polymorpha to study phytochrome signaling.

Neofunctionalisation of basic helix-loop-helix proteins occurred when embryophytes colonised the land

ROOT HAIR DEFECTIVE SIX-LIKE (RSL) genes control the development of structures from single cells at the surface of embryophytes (land plants) such as rhizoids and root hairs. RSL proteins constitute a subclass (VIIIc) of the basic helix-loop-helix (bHLH) class VIII transcription factor family. The Charophyceae form the only class of streptophyte algae with tissue-like structures and rhizoids. To determine if the function of RSL genes in the control of cell differentiation in embryophytes was inherited from a streptophyte algal ancestor, we identified the single class VIII bHLH gene from the charophyceaen alga Chara braunii (CbbHLHVIII). CbbHLHVIII is sister to the RSL proteins; they constitute a monophyletic group. Expression of CbbHLHVIII does not compensate for loss of RSL functions in Marchantia polymorpha or Arabidopsis thaliana. In C. braunii CbbHLHVIII is expressed at sites of morphogenesis but not in rhizoids. This finding indicates that C. braunii class VIII protein is functionally different from land plant RSL proteins. This result suggests that the function of RSL proteins in cell differentiation at the plant surface evolved by neofunctionalisation in the land plants lineage after its divergence from its last common ancestor with C. braunii, at or before the colonisation of the land by embryophytes.

Physiological function of photoreceptor UVR8 in UV-B tolerance in the liverwort Marchantia polymorpha

The physiological importance of MpUVR8 in UV-B resistance and translocation in a UV-B-dependent manner from the cytosol into the nucleus is characterized in Marchantia polymorpha. UV RESISTANCE LOCUS 8 (UVR8) is an ultraviolet-B (UV-B) light receptor functioning for UV-B sensing and tolerance in Arabidopsis thaliana and other species. It is unclear whether UVR8 physiologically functions in UV-B-induced defense responses in Marchantia polymorpha, which belongs to the earliest diverging group of embryophyte lineages. Here, we demonstrate that UVR8 has a physiological function in UV-B tolerance and that there is a UVR8-dependent pathway involved. In addition, a UVR8-independent pathway is revealed. We examine the tissue-specific expression pattern of M. polymorpha UVR8 (MpUVR8), showing that it is highly expressed in the apical notch in thalli and gametangiophores, as well as in antheridial and archegonial heads. Furthermore, Mpuvr8^KO plant transformants, in which the MpUVR8 locus was disrupted, were produced and analyzed to understand the physiological and molecular function of MpUVR8. Analysis using these plants indicates the important roles of MpUVR8 and MpUVR8-regulated genes, and of MpUVR8-independent pathways in UV-B tolerance. Subcellular localization of Citrine-fused MpUVR8 in M. polymorpha cells was also investigated. It was found to translocate from the cytosol into the nucleus in response to UV-B irradiation. Our findings indicate strong conservation of the physiological function of UVR8 and the molecular mechanisms for UVR8-dependent signal transduction through regulation of gene expression in embryophytes.

Reproductive Induction is a Far-Red High Irradiance Response that is Mediated by Phytochrome and PHYTOCHROME INTERACTING FACTOR in Marchantia polymorpha

Land plants have evolved a series of photoreceptors to precisely perceive environmental information. Among these, phytochromes are the sole photoreceptors for red light (R) and far-red light (FR), and play pivotal roles in modulating various developmental processes. Most extant land plants possess multiple phytochromes that probably evolved from a single phytochrome in the common ancestor of land plants. However, the ancestral phytochrome signaling mechanism remains unknown due to a paucity of knowledge regarding phytochrome functions in basal land plants. It has recently been reported that Mpphy, a single phytochrome in the liverwort Marchantia polymorpha, regulates typical photoreversible responses collectively classified as low fluence response (LFR). Here, we show that Mpphy also regulates the gametangiophore formation analogous to the mode of action of the far-red high irradiance response (FR-HIR) in angiosperms. Our phenotypic analyses using mutant plants obtained by CRISPR/Cas9-based genome editing revealed that MpFHY1, an ortholog of FAR-RED ELONGATED HYPOCOTYL1, as well as Mpphy is critical for the FR-HIR signaling in M. polymorpha. In addition, knockout of MpPIF, a single PHYTOCHROME INTERACTING FACTOR gene in M. polymorpha, completely abolished the FR-HIR-dependent gametangiophore formation, while overexpression of MpPIF accelerated the response. FR-HIR-dependent transcriptional regulation was also disrupted in the Mppif mutant. Our findings suggest that plants had already acquired the FR-HIR signaling mediated by phytochrome and PIF at a very early stage during the course of land plant evolution, and that a single phytochrome in the common ancestor of land plants could mediate both LFR and FR-HIR.

An Evolutionarily Conserved Abscisic Acid Signaling Pathway Regulates Dormancy in the Liverwort Marchantia polymorpha

Dormancy is a key process allowing land plants to adapt to changing conditions in the terrestrial habitat, allowing the cessation of growth in response to environmental or physiological cues, entrance into a temporary quiescent state, and subsequent reactivation of growth in more favorable environmental conditions [1-3]. Dormancy may be induced seasonally, sporadically (e.g., in response to drought), or developmentally (e.g., seeds and apical dominance). Asexual propagules, known as gemmae, derived via clonal reproduction in bryophytes, are often dormant until displaced from the parent plant. In the liverwort Marchantia polymorpha, gemmae are produced within specialized receptacles, gemma cups, located on the dorsal side of the vegetative thallus [4]. Mature gemmae are detached from the parent plant but may remain in the cup, with gemma growth suppressed as long as the gemmae remain in the gemma cup and the parental plant is alive [5]. Following dispersal of gemmae from gemma cups by rain, the gemmae germinate in the presence of light and moisture, producing clonal offspring [6]. In land plants, the plant hormone abscisic acid (ABA) regulates many aspects of dormancy and water balance [7]. Here, we demonstrate that ABA plays a central role in the control of gemma dormancy as transgenic M. polymorpha gemmae with reduced sensitivity to ABA fail to establish and/or maintain dormancy. Thus, the common ancestor of land plants used the ABA signaling module to regulate germination of progeny in response to environmental cues, with both gemmae and seeds being derived structures co-opting an ancestral response system.

Negative regulation of conserved RSL class I bHLH transcription factors evolved independently among land plants

Basic helix-loop-helix transcription factors encoded by RSL class I genes control a gene regulatory network that positively regulates the development of filamentous rooting cells – root hairs and rhizoids – in land plants. The GLABRA2 transcription factor negatively regulates these genes in the angiosperm Arabidopsis thaliana. To find negative regulators of RSL class I genes in early diverging land plants we conducted a mutant screen in the liverwort Marchantia polymorpha. This identified FEW RHIZOIDS1 (MpFRH1) microRNA (miRNA) that negatively regulates the RSL class I gene MpRSL1. The miRNA and its mRNA target constitute a feedback mechanism that controls epidermal cell differentiation. MpFRH1 miRNA target sites are conserved among liverwort RSL class I mRNAs but are not present in RSL class I mRNAs of other land plants. These findings indicate that while RSL class I genes are ancient and conserved, independent negative regulatory mechanisms evolved in different lineages during land plant evolution.

Plants colonised the land sometime more than 500 million years ago. The ancestors of the first land plants were algae that were most likely simple with a few different types of cell. Yet, when faced with the challenges of life on land, plants evolved new cell types and specialised structures with roles such as anchorage, nutrient uptake and gas exchange.

Many of these specialised structures, including the root hairs and rhizoids that allow plants to collect water and minerals from the soil, first develop as outgrowths from cells in the outer layer of the plant. An ancient and conserved mechanism activates the development of these outgrowths via genes belonging to a group known as RSL class I.

In the flowering plant Arabidopsis thaliana, a protein switches off RSL class I genes in a subset of these outer cells, to stop too many root hairs forming. To see whether this kind of negative regulation is also conserved among land plants, Honkanen et al. looked for regulators of RSL class I genes in liverworts. Small and without flowers, liverworts are a group of plants that first appeared during the earliest stages of land plant evolution.

Honkanen et al. discovered that RSL class I genes in liverworts are negatively regulated by a molecule named FEW RHIZOIDS1 (or FRH1). However, rather than being a protein, FRH1 is a microRNA – a short strand of genetic code that reduces how much protein is produced from a given gene. The FRH1 microRNA is conserved among liverworts and most likely evolved very early in the history of these plants. The findings indicate that different groups of land plants have evolved different negative regulators to control the conserved genes behind some of the specialised structures crucial to life on land.

Identification of reference genes for RT-qPCR in the Antarctic moss Sanionia uncinata under abiotic stress conditions

Sanionia uncinata is a dominant moss species in the maritime Antarctic. Due to its high adaptability to harsh environments, this extremophile plant has been considered a good target for studying the molecular adaptation mechanisms of plants to a variety of environmental stresses. Despite the importance of S. uncinata as a representative Antarctic plant species for the identification and characterization of genes associated with abiotic stress tolerance, suitable reference genes, which are critical for RT-qPCR analyses, have not yet been identified. In this report, 11 traditionally used and 6 novel candidate reference genes were selected from transcriptome data of S. uncinata and the expression stability of these genes was evaluated under various abiotic stress conditions using three statistical algorithms; geNorm, NormFinder, and BestKeeper. The stability ranking analysis selected the best reference genes depending on the stress conditions. Among the 17 candidates, the most stable references were POB1 and UFD2 for cold stress, POB1 and AKB for drought treatment, and UFD2 and AKB for the field samples from a different water contents in Antarctica. Overall, novel genes POB1 and AKB were the most reliable references across all samples, irrespective of experimental conditions. In addition, 6 novel candidate genes including AKB, POB1 and UFD2, were more stable than the housekeeping genes traditionally used for internal controls, indicating that transcriptome data can be useful for identifying novel robust normalizers. The reference genes validated in this study will be useful for improving the accuracy of RT-qPCR analysis for gene expression studies of S. uncinata in Antarctica and for further functional genomic analysis of bryophytes.

MpAMT1;2 from Marchantia polymorpha is a High-Affinity, Plasma Membrane Ammonium Transporter

Plant ammonium transporters in the AMT/MEP/Rh (ammonium transporter/methylammonium and ammonium permease/Rhesus factor) superfamily have only been previously characterized in flowering plants (angiosperms). Plant AMT1s are electrogenic, while plant AMT2s are electroneutral, and MEP and Rh transporters in other organisms are electroneutral. We analyzed the transport function of MpAMT1;2 from the basal land plant Marchantia polymorpha, a liverwort. MpAMT1;2 was shown to localize to the plasma membrane in Marchantia gametophyte thallus by stable transformation using a C-terminal citrine fusion. MpAMT1;2 expression was studied using quantitative real-time PCR and shown to be higher when plants were N deficient and lower when plants were grown on media containing ammonium, nitrate or the amino acid glutamine. Expression in Xenopus oocytes and analysis by electrophysiology revealed that MpAMT1;2 is an electrogenic ammonium transporter with a very high affinity for ammonium (7 µM at pH 5.6 and a membrane potential of -137 mV). A conserved inhibitory phosphorylation site identified in angiosperm AMT1s is also present in all AMT1s in Marchantia. Here we show that a phosphomimetic mutation T475D in MpAMT1;2 completely inhibits ammonium transport activity. The results indicate that MpAMT1;2 may be important for ammonium uptake into cells in the Marchantia thallus.

Origin and evolution of the nuclear auxin response system

The small signaling molecule auxin controls numerous developmental processes in land plants, acting mostly by regulating gene expression. Auxin response proteins are represented by large families of diverse functions, but neither their origin nor their evolution is understood. Here, we use a deep phylogenomics approach to reconstruct both the origin and the evolutionary trajectory of all nuclear auxin response protein families. We found that, while all subdomains are ancient, a complete auxin response mechanism is limited to land plants. Functional phylogenomics predicts defined steps in the evolution of response system properties, and comparative transcriptomics across six ancient lineages revealed how these innovations shaped a sophisticated response mechanism. Genetic analysis in a basal land plant revealed unexpected contributions of ancient non-canonical proteins in auxin response as well as auxin-unrelated function of core transcription factors. Our study provides a functional evolutionary framework for understanding diverse functions of the auxin signal.

Across all kingdoms of life, signaling molecules like hormones, for example, control many aspects of the lives of organisms, including how they grow and develop. Cells have dedicated proteins that can recognize the signaling molecules, relay the information, and respond to the signal, for example by switching genes on or off. Such response systems usually consist of multiple components, and, throughout evolution, these response components have regularly been copied such that many species have multiple different versions of each one.

Auxin is a plant hormone that controls virtually all growth and developmental processes in plants, including many yield traits in crops. However, no one knows why it is involved in so many processes. This is partly because it is not clear how the response system for this central signaling molecule was first born, or how it has increased in its complexity.

To address this, Mutte, Kato et al. explored the genetic information of more than a thousand plant species, including algae, which span more than 700 million years of evolution. Their analysis showed that all auxin response components were assembled from pieces of much older genes, but that they first came together when plants conquered land. Indeed, the auxin response appears to have developed on top of a pre-existing genetic regulator that is still present in modern-day algae. Mutte, Kato et al. then used experiments to show how stepwise increases in the number and types of auxin response components have shaped sophisticated, complex responses in land plants, and to demonstrate how ancient components control auxin response.

Together these findings provide a framework for understanding the many functions of auxin in plants, and how this came to be. They also show how complexity can be accomplished in a signal response pathway, and how diversity evolves in gene families. Similar studies on other response systems in plants and beyond are likely to help reveal common principles of hormone response evolution and diversification of gene regulation systems.

Origin and evolution of the nuclear auxin response system

The small signaling molecule auxin controls numerous developmental processes in land plants, acting mostly by regulating gene expression. Auxin response proteins are represented by large families of diverse functions, but neither their origin nor their evolution is understood. Here, we use a deep phylogenomics approach to reconstruct both the origin and the evolutionary trajectory of all nuclear auxin response protein families. We found that, while all subdomains are ancient, a complete auxin response mechanism is limited to land plants. Functional phylogenomics predicts defined steps in the evolution of response system properties, and comparative transcriptomics across six ancient lineages revealed how these innovations shaped a sophisticated response mechanism. Genetic analysis in a basal land plant revealed unexpected contributions of ancient non-canonical proteins in auxin response as well as auxin-unrelated function of core transcription factors. Our study provides a functional evolutionary framework for understanding diverse functions of the auxin signal.

An evolutionarily conserved NIMA-related kinase directs rhizoid tip growth in the basal land plant Marchantia polymorpha

Tip growth is driven by turgor pressure and mediated by the polarized accumulation of cellular materials. How a single polarized growth site is established and maintained is unclear. Here, we analyzed the function of NIMA-related protein kinase 1 (MpNEK1) in the liverwort Marchantia polymorpha In the wild type, rhizoid cells differentiate from the ventral epidermis and elongate through tip growth to form hair-like protrusions. In Mpnek1 knockout mutants, rhizoids underwent frequent changes in growth direction, resulting in a twisted and/or spiral morphology. The functional MpNEK1-Citrine protein fusion localized to microtubule foci in the apical growing region of rhizoids. Mpnek1 knockouts exhibited increases in both microtubule density and bundling in the apical dome of rhizoids. Treatment with the microtubule-stabilizing drug taxol phenocopied the Mpnek1 knockout. These results suggest that MpNEK1 directs tip growth in rhizoids through microtubule organization. Furthermore, MpNEK1 expression rescued ectopic outgrowth of epidermal cells in the Arabidopsis thaliana nek6 mutant, strongly supporting an evolutionarily conserved NEK-dependent mechanism of directional growth. It is possible that such a mechanism contributed to the evolution of the early rooting system in land plants.

Extensive epigenetic reprogramming during the life cycle of Marchantia polymorpha

Background

In plants, the existence and possible role of epigenetic reprogramming has been questioned because of the occurrence of stably inherited epialleles. Evidence suggests that epigenetic reprogramming does occur during land plant reproduction, but there is little consensus on the generality and extent of epigenetic reprogramming in plants. We studied DNA methylation dynamics during the life cycle of the liverwort Marchantia polymorpha. We isolated thalli and meristems from male and female gametophytes, archegonia, antherozoids, as well as sporophytes at early and late developmental stages, and compared their DNA methylation profiles.

Results

Of all cytosines tested for differential DNA methylation, 42% vary significantly in their methylation pattern throughout the life cycle. However, the differences are limited to few comparisons between specific stages of the life cycle and suggest four major epigenetic states specific to sporophytes, vegetative gametophytes, antherozoids, and archegonia. Further analyses indicated clear differences in the mechanisms underlying reprogramming in the gametophytic and sporophytic generations, which are paralleled by differences in the expression of genes involved in DNA methylation. Differentially methylated cytosines with a gain in methylation in antherozoids and archegonia are enriched in the CG and CHG contexts, as well as in gene bodies and gene flanking regions. In contrast, gain of DNA methylation during sporophyte development is mostly limited to the CHH context, LTR retrotransposons, DNA transposons, and repeats.

Conclusion

We conclude that epigenetic reprogramming occurs at least twice during the life cycle of M. polymorpha and that the underlying mechanisms are likely different between the two events.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1383-z) contains supplementary material, which is available to authorized users.

Evolutionary conservation of structure and function of the UVR8 photoreceptor from the liverwort Marchantia polymorpha and the moss Physcomitrella patens

The ultraviolet-B (UV-B) photoreceptor UV RESISTANCE LOCUS 8 (UVR8) mediates photomorphogenic responses to UV-B in Arabidopsis through differential gene expression, but little is known about UVR8 in other species. Bryophyte lineages were the earliest diverging embryophytes, thus being the first plants facing the UV-B regime typical of land. We therefore examined whether liverwort and moss species have functional UVR8 proteins and whether they are regulated similarly to Arabidopsis UVR8. We examined the expression, dimer/monomer status, cellular localisation and function of Marchantia polymorpha and Physcomitrella patens UVR8 in experiments with bryophyte tissue and expression of green fluorescent protein (GFP)-UVR8 fusions in Nicotiana leaves and transgenic Arabidopsis. P. patens expresses two UVR8 genes that encode functional proteins, whereas the single M. polymorpha UVR8 gene expresses two transcripts by alternative splicing that encode functional UVR8 variants. P. patens UVR8 proteins form dimers that monomerise and accumulate in the nucleus following UV-B exposure, similar to Arabidopsis UVR8, but M. polymorpha UVR8 has weaker dimers and the proteins appear more constitutively nuclear. We conclude that liverwort and moss species produce functional UVR8 proteins. Although there are differences in expression and regulation of P. patens and M. polymorpha UVR8, the mechanism of UVR8 action is strongly conserved in evolution.

Early evolution of the land plant circadian clock

While angiosperm clocks can be described as an intricate network of interlocked transcriptional feedback loops, clocks of green algae have been modelled as a loop of only two genes. To investigate the transition from a simple clock in algae to a complex one in angiosperms, we performed an inventory of circadian clock genes in bryophytes and charophytes. Additionally, we performed functional characterization of putative core clock genes in the liverwort Marchantia polymorpha and the hornwort Anthoceros agrestis. Phylogenetic construction was combined with studies of spatiotemporal expression patterns and analysis of M. polymorpha clock gene mutants. Homologues to core clock genes identified in Arabidopsis were found not only in bryophytes but also in charophytes, albeit in fewer copies. Circadian rhythms were detected for most identified genes in M. polymorpha and A. agrestis, and mutant analysis supports a role for putative clock genes in M. polymorpha. Our data are in line with a recent hypothesis that adaptation to terrestrial life occurred earlier than previously expected in the evolutionary history of charophyte algae. Both gene duplication and acquisition of new genes was important in the evolution of the plant circadian clock, but gene loss has also contributed to shaping the clock of bryophytes.

Functional PTB phosphate transporters are present in streptophyte algae and early diverging land plants

Two inorganic phosphate (Pi) uptake mechanisms operate in streptophytes and chlorophytes, the two lineages of green plants. PHOSPHATE TRANSPORTER B (PTB) proteins are hypothesized to be the Na⁺ /Pi symporters catalysing Pi uptake in chlorophytes, whereas PHOSPHATE TRANSPORTER 1 (PHT1) proteins are the H⁺ /Pi symporters that carry out Pi uptake in angiosperms. PHT1 proteins are present in all streptophyte lineages. However, Pi uptake in streptophyte algae and marine angiosperms requires Na⁺ influx, suggesting that Na⁺ /Pi symporters also function in some streptophytes. We tested the hypothesis that Na⁺ /Pi symporters exist in streptophytes. We identified PTB sequences in streptophyte genomes. Core PTB proteins are present at the plasma membrane of the liverwort Marchantia polymorpha. The expression of M. polymorpha core PTB proteins in the Saccharomyces cerevisiae pho2 mutant defective in high-affinity Pi transport rescues growth in low-Pi environments. Moreover, levels of core PTB mRNAs of M. polymorpha and the streptophyte alga Coleochaete nitellarum are higher in low-Pi than in Pi-replete conditions, consistent with a role in Pi uptake from the environment. We conclude that land plants inherited two Pi uptake mechanisms - mediated by the PTB and PHT1 proteins, respectively - from their streptophyte algal ancestor. Both systems operate in parallel in extant early diverging land plants.

MpWIP regulates air pore complex development in the liverwort Marchantia polymorpha

The colonisation of the land by plants was accompanied by the evolution of complex tissues and multicellular structures comprising different cell types as morphological adaptations to the terrestrial environment. Here, we show that the single WIP protein in the early-diverging land plant Marchantia polymorpha L. is required for the development of the multicellular gas exchange structure: the air pore complex. This 16-cell barrel-shaped structure surrounds an opening between epidermal cells that facilitates the exchange of gases between the chamber containing the photosynthetic cells inside the plant and the air outside. MpWIP is expressed in cells of the developing air pore complex and the morphogenesis of the complex is defective in plants with reduced MpWIP function. The role of WIP proteins in the control of different multicellular structures in M. polymorpha and the flowering plant Arabidopsis thaliana suggests that these proteins controlled the development of multicellular structures in the common ancestor of land plants. We hypothesise that WIP genes were subsequently co-opted in the control of morphogenesis of novel multicellular structures that evolved during the diversification of land plants.

Highlighted Article: The MpWIP gene controls the development of the air pore complex – a multicellular structure that increases CO₂ uptake – in the early-diverging land plant Marchantia polymorpha.

Selection and Validation of Reference Genes for Gene Expression Analysis in Vigna angularis Using Quantitative Real-Time RT-PCR

Adzuki bean (Vigna angularis) is one of the most important legume crops in Asian countries like China, Japan and Korea due to its nutritious protein and starch contents. In spite of its economic importance, gene expression analysis system for gene function verification of adzuki bean is still absent. Therefore, reference genes for gene expression analysis based on the quantitative real time PCR (qRT-PCR) were screened in current study. A total of nine general housekeeping genes, including ACT, Fbox, ZMPP, GAPDH, EF, PP2A, UBC, UBN and PTB were evaluated for their expression stability by qRT-PCR in four adzuki bean cultivars, three different tissues, four abiotic stress and one biotic stress. The best group of candidates as reference genes were as follows: PTB and ACT for different cultivars; EF and UBN for different tissues; ACT and ZMPP for biotic stress and waterlogging stress; Fbox and UBC for salinity-alkalinity stress; Fbox and PTB for drought stress. Our results will provide a more accurate and reliable normalization of qRT-PCR data in adzuki bean.

Phytochrome Signaling Is Mediated by PHYTOCHROME INTERACTING FACTOR in the Liverwort Marchantia polymorpha

Phytochromes are red light (R) and far-red light (FR) receptors that play important roles in many aspects of plant growth and development. Phytochromes mainly function in the nucleus and regulate sets of genes by inhibiting negatively acting basic helix-loop-helix transcription factors named PHYTOCHROME INTERACTING FACTORs (PIFs) in Arabidopsis thaliana Although R/FR photoreversible responses and phytochrome genes are well documented in diverse lineages of plants, the extent to which phytochrome signaling is mediated by gene regulation beyond angiosperms remains largely unclear. Here, we show that the liverwort Marchantia polymorpha, an emerging model basal land plant, has only one phytochrome gene, Mp-PHY, and only one PIF gene, Mp-PIF These genes mediate typical low fluence responses, which are reversibly elicited by R and FR, and regulate gene expression. Mp-phy is light-stable and translocates into the nucleus upon irradiation with either R or FR, demonstrating that the single phytochrome Mp-phy exhibits combined biochemical and cell-biological characteristics of type I and type II phytochromes. Mp-phy photoreversibly regulates gemma germination and downstream gene expression by interacting with Mp-PIF and targeting it for degradation in an R-dependent manner. Our findings suggest that the molecular mechanisms for light-dependent transcriptional regulation mediated by PIF transcription factors were established early in land plant evolution.