Biological implications of the occurrence of 32 members of the XTH (xyloglucan endotransglucosylase/hydrolase) family of proteins in the bryophyte Physcomitrella patens

This comprehensive overview of the xyloglucan endotransglucosylase/hydrolase (XTH) family of genes and proteins in bryophytes, based on research using genomic resources that are newly available for the moss Physcomitrella patens, provides new insights into plant evolution. In angiosperms, the XTH genes are found in large multi‐gene families, probably reflecting the diverse roles of individual XTHs in various cell types. As there are fewer cell types in P. patens than in angiosperms such as Arabidopsis and rice, it is tempting to deduce that there are fewer XTH family genes in bryophytes. However, the present study unexpectedly identified as many as 32 genes that potentially encode XTH family proteins in the genome of P. patens, constituting a fairly large multi‐gene family that is comparable in size with those of Arabidopsis and rice. In situ localization of xyloglucan endotransglucosylase activity in this moss indicates that some P. patens XTH proteins exhibit biochemical functions similar to those found in angiosperms, and that their expression profiles are tissue‐dependent. However, comparison of structural features of families of XTH genes between P. patens and angiosperms demonstrated the existence of several bryophyte‐specific XTH genes with distinct structural and functional features that are not found in angiosperms. These bryophyte‐specific XTH genes might have evolved to meet morphological and functional needs specific to the bryophyte. These findings raise interesting questions about the biological implications of the XTH family of proteins in non‐seed plants.     

A pseudomolecule‐scale genome assembly of the liverwort Marchantia polymorpha

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

The cellulose synthase (<Emphasis Type="Italic">CESA</Emphasis>) gene superfamily of the moss <Emphasis Type="Italic">Physcomitrella patens</Emphasis>

The CESA gene superfamily of Arabidopsis and other seed plants comprises the CESA family, which encodes the catalytic subunits of cellulose synthase, and eight families of CESA-like (CSL) genes whose functions are largely unknown. The CSL genes have been proposed to encode processive β-glycosyl transferases that synthesize noncellulosic cell wall polysaccharides. BLAST searches of EST and shotgun genomic sequences from the moss Physcomitrella patens (Hedw.) B.S.G. were used to identify genes with high similarity to vascular plant CESAs, CSLAs, CSLCs, and CSLDs. However, searches using Arabidopsis CSLBs, CSLEs, and CSLGs or rice CSLFs or CSLHs as queries identified no additional CESA superfamily members in P. patens, indicating that this moss lacks representatives of these families. Intron insertion sites are highly conserved between Arabidopsis and P. patens in all four shared gene families. However, phylogenetic analysis strongly supports independent diversification of the shared families in mosses and vascular plants. The lack of orthologs of vascular plant CESAs in the P. patens genome indicates that the divergence of mosses and vascular plants predated divergence and specialization of CESAs for primary and secondary cell wall syntheses and for distinct roles within the rosette terminal complexes. In contrast to Arabidopsis, the CSLD family is highly represented among P. patens ESTs. This is consistent with the proposed function of CSLDs in tip growth and the central role of tip growth in the development of the moss protonema. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users. Accession numbers: DQ417756, DQ417757, DQ898284–6, DQ898147–54, DQ902545–51.     

CYP710A genes encoding sterol C22-desaturase in <Emphasis Type="Italic">Physcomitrella patens</Emphasis> as molecular evidence for the evolutionary conservation of a sterol biosynthetic pathway in plants

We have characterized cytochromes P450, CYP710A13, and CYP710A14, as the sterol C22-desaturase in the moss Physcomitrella patens. GC–MS analyses demonstrated that P. patens accumulated stigmasterol as the major sterol (56–60% of total sterol) and sitosterol to a lesser extent (8–12%); this sterol profile contrasts with those in higher plants accumulating stigmasterol as a minor component. Recombinant CYP710A13 and CYP710A14 proteins prepared using a baculovirus/insect cell system exhibited the C22-desaturase activity with β-sitosterol to produce stigmasterol, while campesterol and 24-epi-campesterol were not accepted as the substrates. The K _m values for β-sitosterol of CYP710A13 (1.0 ± 0.043 μM) and CYP710A14 (2.1 ± 0.17 μM) were at comparable levels of those reported with higher plant CYP710A proteins. In Arabidopsis T87 cells over-expressing CYP710A14, stigmasterol contents reached a level 20- to 72-fold higher than those in the basal level of T87 cells, confirming the C22-desaturase activity of this P450 enzyme. The occurrence of the end-products together with the enzymes involved in the last step of the pathway substantiated the presence of an entire sterol biosynthetic pathway in P. patens, providing evidence for the conservation of the sterol biosynthetic pathway through the evolutionary process of land plants. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Does the core circadian clock in the moss <Emphasis Type="Italic">Physcomitrella patens</Emphasis> (Bryophyta) comprise a single loop?

Background  

The endogenous circadian clock allows the organism to synchronize processes both to daily and seasonal changes. In plants, many metabolic processes such as photosynthesis, as well as photoperiodic responses, are under the control of a circadian clock. Comparative studies with the moss Physcomitrella patens provide the opportunity to study many aspects of land plant evolution. Here we present a comparative overview of clock-associated components and the circadian network in the moss P. patens.     

Role of ENA ATPase in Na<Superscript>+</Superscript> efflux at high pH in bryophytes

Potassium or Na⁺ efflux ATPases, ENA ATPases, are present in all fungi and play a central role in Na⁺ efflux and Na⁺ tolerance. Flowering plants lack ENA ATPases but two ENA ATPases have been identified in the moss Physcomitrella patens, PpENA1 and PpENA2. PpENA1 mediates Na⁺ efflux in Saccharomyces cerevisiae. To propose a general function of ENA ATPases in bryophytes it was necessary to demonstrate that these ATPases mediate Na⁺ efflux in planta and that they exist in more bryophytes than P. patens. For these demonstrations (1) we cloned a third ATPase from P. patens, PpENA3, and studied the expression pattern of the three PpENA genes; (2) we constructed and studied the single and double Δppena1 and Δppena2 mutants; and (3) we cloned two ENA ATPases from the liverwort Marchantia polymorpha, MpENA1 and MpENA2, and expressed them in S. cerevisiae. The results from the first two approaches revealed that the expression of ENA ATPases was greatly enhanced at high pH and that Na⁺ efflux at high pH depended on PpENA1. The ENA1 ATPase of M. polymorpha suppressed the defective growth of a S. cerevisiae mutant at high K⁺ or Na⁺ concentrations, especially at high K⁺.     

Functional analyses of the ABI1-related protein phosphatase type 2C reveal evolutionarily conserved regulation of abscisic acid signaling between Arabidopsis and the moss <Emphasis Type="Italic">Physcomitrella patens</Emphasis>

We employed a comparative genomic approach to understand protein phosphatase 2C (PP2C)-mediated abscisic acid (ABA) signaling in the moss Physcomitrella patens. Ectopic expression of Arabidopsis (Arabidopsis thaliana) abi1-1, a dominant mutant allele of ABI1 encoding a PP2C involved in the negative regulation of ABA signaling, caused ABA insensitivity of P. patens both in gene expression of late embryogenesis abundant (LEA) genes and in ABA-induced protonemal growth inhibition. The transgenic abi1-1 plants showed decreased ABA-induced freezing tolerance, and decreased tolerance to osmotic stress. Analyses of the P. patens genome revealed that only two (PpABI1A and PpABI1B) PP2C genes were related to ABI1. In the ppabi1a null mutants, ABA-induced expression of LEA genes was elevated, and protonemal growth was inhibited with lower ABA concentration compared to the wild type. Moreover, ABA-induced freezing tolerance of the ppabi1a mutants was markedly enhanced. We provide the genetic evidence that PP2C-mediated ABA signaling is evolutionarily conserved between Arabidopsis and P. patens. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Accession Numbers: PpABI1A-AB369256, PpABI1B-AB369255, pphn39k21-AB369257.     

Coding sequence divergence between two closely related plant species: Arabidopsis thaliana and Brassica rapa ssp. pekinensis

To characterize the coding-sequence divergence of closely related genomes, we compared DNA sequence divergence between sequences from a Brassica rapa ssp. pekinensis EST library isolated from flower buds and genomic sequences from Arabidopsis thaliana. The specific objectives were (i) to determine the distribution of and relationship between K _a and K _s, (ii) to identify genes with the lowest and highest K _a:K _s values, and (iii) to evaluate how codon usage has diverged between two closely related species. We found that the distribution of K _a:K _s was unimodal, and that substitution rates were more variable at nonsynonymous than synonymous sites, and detected no evidence that K _a and K _s were positively correlated. Several genes had K _a:K _s values equal to or near zero, as expected for genes that have evolved under strong selective constraint. In contrast, there were no genes with K _a:K _s >1 and thus we found no strong evidence that any of the 218 sequences we analyzed have evolved in response to positive selection. We detected a stronger codon bias but a lower frequency of GC at synonymous sites in A. thaliana than B. rapa. Moreover, there has been a shift in the profile of most commonly used synonymous codons since these two species diverged from one another. This shift in codon usage may have been caused by stronger selection acting on codon usage or by a shift in the direction of mutational bias in the B. rapa phylogenetic lineage.     

Turnover of recently assimilated carbon in arctic bryophytes

Carbon (C) allocation and turnover in arctic bryophytes is largely unknown, but their response to climatic change has potentially significant impacts on arctic ecosystem C budgets. Using a combination of pulse-chase experiments and a newly developed model of C turnover in bryophytes, we show significant differences in C turnover between two contrasting arctic moss species (Polytrichum piliferum and Sphagnum fuscum). ¹³C abundance in moss tissues (measured up to 1 year) and respired CO₂ (traced over 5 days) were used to parameterise the bryophyte C model with four pools representing labile and structural C in photosynthetic and stem tissue. The model was optimised using an Ensemble Kalman Filter to ensure a focus on estimating the confidence intervals (CI) on model parameters and outputs. The ratio of aboveground NPP:GPP in Polytrichum piliferum was 23% (CI 9–35%), with an average turnover time of 1.7 days (CI 1.1–2.5 days). The aboveground NPP:GPP ratio in Sphagnum fuscum was 43% (CI 19–65%) with an average turnover time of 3.1 days (CI 1.6–6.1 days). These results are the first to show differences in C partitioning between arctic bryophyte species in situ and highlight the importance of modelling C dynamics of this group separately from vascular plants for a realistic representation of vegetation in arctic C models.     

The bryophytes Physcomitrium patens and Marchantia polymorpha as model systems for studying evolutionary cell and developmental biology in plants

Bryophytes are nonvascular spore-forming plants. Unlike in flowering plants, the gametophyte (haploid) generation of bryophytes dominates the sporophyte (diploid) generation. A comparison of bryophytes with flowering plants allows us to answer some fundamental questions raised in evolutionary cell and developmental biology. The moss Physcomitrium patens was the first bryophyte with a sequenced genome. Many cell and developmental studies have been conducted in this species using gene targeting by homologous recombination. The liverwort Marchantia polymorpha has recently emerged as an excellent model system with low genomic redundancy in most of its regulatory pathways. With the development of molecular genetic tools such as efficient genome editing, both P. patens and M. polymorpha have provided many valuable insights. Here, we review these advances with a special focus on polarity formation at the cell and tissue levels. We examine current knowledge regarding the cellular mechanisms of polarized cell elongation and cell division, including symmetric and asymmetric cell division. We also examine the role of polar auxin transport in mosses and liverworts. Finally, we discuss the future of evolutionary cell and developmental biological studies in plants.

A review of the cell biological and developmental mechanisms of bryophytes, including Physcomitrium patens and Marchantia polymorpha.     

Prenylation is required for polar cell elongation,cell adhesion,and differentiation in Physcomitrella patens

Protein prenylation is required for a variety of growth and developmental processes in flowering plants. Here we report the consequences of loss of function of all known prenylation subunits in the moss Physcomitrella patens. As in Arabidopsis, protein farnesyltransferase and protein geranylgeranyltransferase type I are not required for viability. However, protein geranylgeranyltransferase type I activity is required for cell adhesion, polar cell elongation, and cell differentiation. Loss of protein geranylgeranyltransferase activity results in colonies of round, single‐celled organisms that resemble unicellular algae. The loss of protein farnesylation is not as severe but also results in polar cell elongation and differentiation defects. The complete loss of Rab geranylgeranyltransferase activity appears to be lethal in P. patens. Labeling with antibodies to cell wall components support the lack of polarity establishment and the undifferentiated state of geranylgeranyltransferase type I mutant plants. Our results show that prenylated proteins play key roles in P. patens development and differentiation processes.     

Heterologous expression of cDNAs encoding monodehydroascorbate reductases from the moss, <Emphasis Type="Italic">Physcomitrella patens</Emphasis> and characterization of the expressed enzymes

Monodehydroascorbate reductase (MDHAR; EC 1.6.5.4) catalyses the reduction of the monodehydroascorbate (MDHA) radical to ascorbate, using NADH or NADPH as an electron donor, and is believed to be involved in maintaining the reactive oxygen scavenging capability of plant cells. This key enzyme in the ascorbate-glutathione cycle has been studied here in the moss Physcomitrella patens, which is tolerant to a range of abiotic stresses and is increasingly used as a model plant. In the present study, three cDNAs encoding different MDHAR isoforms of 47 kDa were identified in P. patens, and found to exhibit enzymic characteristics similar to MDHARs in vascular plants despite low-sequence identity and a distant evolutionary relationship between the species. The three cDNAs for the P. patens MDHAR enzymes were expressed in Escherichia coli and the active enzymes were purified and characterized. Each recombinant protein displayed an absorbance spectrum typical of flavoenzymes and contained a single non-covalently bound FAD coenzyme molecule. The K _m and k _cat values for the heterologously expressed PpMDHAR enzymes ranged from 8 to 18 μM and 120–130 s⁻¹, respectively, using NADH as the electron donor. The K _m values were at least an order of magnitude higher for NADPH. The K _m values for the MDHA radical were ∼0.5–1.0 μM for each of the purified enzymes, and further kinetic analyses indicated that PpMDHARs follow a ‘ping–pong’ kinetic mechanism. In contrast to previously published data, site-directed mutagenesis indicated that the conserved cysteine residue is not directly involved in the reduction of MDHA.     

Diversification of ftsZ During Early Land Plant Evolution

The plastid division proteins FtsZ are encoded by a small nuclear gene family in land plants. Although it has been shown for some of the gene products that they are imported into plastids and function in plastid division, the evolution and function of this gene family and their products remain to be unraveled. Here we present two new ftsZ genes from the moss Physcomitrella patens and compare the genomic structure of members of the two plant ftsZ gene families. Comparison of sequence features and phylogenetic analyses confirm the presence of two clusters of paralogues in land plants and demonstrate that these genes were duplicated before the divergence of mosses, ferns and seed plants.     

Analysis of synonymous codon usage patterns in different plant mitochondrial genomes

Codon usage in mitochondrial genome of the six different plants was analyzed to find general patterns of codon usage in plant mitochondrial genomes. The neutrality analysis indicated that the codon usage patterns of mitochondrial genes were more conserved in GC content and no correlation between GC12 and GC3. T and A ending codons were detected as the preferred codons in plant mitochondrial genomes. The Parity Rule 2 plot analysis showed that T was used more frequently than A. The EN_C-plot showed that although a majority of the points with low EN_C values were lying below the expected curve, a few genes lied on the expected curve. Correspondence analysis of relative synonymous codon usage yielded a first axis that explained only a partial amount of variation of codon usage. These findings suggest that natural selection is likely to be playing a large role in codon usage bias in plant mitochondrial genomes, but not only natural selection but also other several factors are likely to be involved in determining the selective constraints on codon bias in plant mitochondrial genomes. Meantime, 1 codon (P. patens), 6 codons (Z. mays), 9 codons (T. aestivum), 15 codons (A. thaliana), 15 codons (M. polymorpha) and 15 codons (N. tabacum) were defined as the preferred codons of the six plant mitochondrial genomes.     

Adaptation to metal-contaminated soils in populations of the moss,Ceratodon purpureus: vegetative growth and reproductive expression

Many observations suggest that morphological evolution occurs slowly in bryophytes, and this has been suggested to reflect low genetic diversity within species. Isozyme studies, however, stand in apparent contrast and have shown that bryophytes can contain high levels of genetic variability within and among populations. In light of this conflict, we tested the potential of the moss, Ceratodon purpureus, to undergo adaptive change (i.e., ecotypic differentiation) in response to soils that have been contaminated with high levels of metals for 90 years by measuring gametophytic growth and reproductive expression under experimental conditions. Variation in protonemal growth in sterile culture indicates that plants from one population growing on contaminated soil near a smelter are significantly more tolerant of zinc, cadmium, and lead than plants from uncontaminated sites. Results from a common garden experiment, in which plants were grown on soil from the smelter site, indicate that plants from near the smelter are significantly more tolerant of contaminated soils than plants from uncontaminated sites for vegetative growth. The same experiment suggests that plants from the smelter site are also more tolerant in terms of gametangial production (although we could not test this statistically). Our results demonstrate that C. purpureus has been able to undergo relatively rapid evolution in response to strong selective pressures.     

<Emphasis Type="Italic">Germin-like protein</Emphasis> gene family of a moss, <Emphasis Type="Italic">Physcomitrella patens</Emphasis>, phylogenetically falls into two characteristic new clades

We identified 77 EST clones encoding germin-like proteins (GLPs) from a moss, Physcomitrella patens in a database search. These Physcomitrella GLPs (PpGLPs) were separated into seven groups based on DNA sequence homology. Phylogenetic analysis showed that these groups were divided into two novel clades clearly distinguishable from higher plant germins and GLPs, named bryophyte subfamilies 1 and 2. PpGLPs belonging to bryophyte subfamilies 1 lacked two cysteines at the conserved positions observed in higher plant germins or GLPs. PpGLPs belonging to bryophyte subfamily 2 contained two cysteines as observed in higher plant germins and GLPs. In bryophyte subfamily 1, 12 amino acids, in which one of two cysteines is included, were deleted between boxes A and B. Further, we determined the genomic structure of all of seven PpGLP genes. The sequences of PpGLPs of bryophyte subfamily 1 contained one or two introns, whereas those of bryophyte subfamily 2 contained no introns. Other GLPs from bryophytes, a liverwort GLP from Marchantia polymorpha, and two moss GLPs from Barbula unguiculata and Ceratodon purpureus also fell into bryophyte subfamily 1 and bryophyte subfamily 2, respectively. No higher plant germins and GLPs were grouped into the bryophyte subfamilies 1 and 2 by our analysis. Moreover, we revealed that PpGLP6 had manganese-containing extracellular superoxide dismutase activity. These results indicated that bryophyte possess characteristic GLPs, which phylogenetically are clearly distinguishable from higher plant GLPs.