An ancient genome duplication contributed to the abundance of metabolic genes in the moss Physcomitrella patens

Molecular Mechanisms of Plant Responses to Copper: From Deficiency to Excess

Copper (Cu) is an essential nutrient for plant growth and development. This metal serves as a constituent element or enzyme cofactor that participates in many biochemical pathways and plays a key role in photosynthesis, respiration, ethylene sensing, and antioxidant systems. The physiological significance of Cu uptake and compartmentalization in plants has been underestimated, despite the importance of Cu in cellular metabolic processes. As a micronutrient, Cu has low cellular requirements in plants. However, its bioavailability may be significantly reduced in alkaline or organic matter-rich soils. Cu deficiency is a severe and widespread nutritional disorder that affects plants. In contrast, excessive levels of available Cu in soil can inhibit plant photosynthesis and induce cellular oxidative stress. This can affect plant productivity and potentially pose serious health risks to humans via bioaccumulation in the food chain. Plants have evolved mechanisms to strictly regulate Cu uptake, transport, and cellular homeostasis during long-term environmental adaptation. This review provides a comprehensive overview of the diverse functions of Cu chelators, chaperones, and transporters involved in Cu homeostasis and their regulatory mechanisms in plant responses to varying Cu availability conditions. Finally, we identified that future research needs to enhance our understanding of the mechanisms regulating Cu deficiency or stress in plants. This will pave the way for improving the Cu utilization efficiency and/or Cu tolerance of crops grown in alkaline or Cu-contaminated soils.

Near telomere-to-telomere genome of the model plant Physcomitrium patens

The model plant Physcomitrium patens has played a pivotal role in enhancing our comprehension of plant evolution and development. However, the current genome harbours numerous regions that remain unfinished and erroneous. To address these issues, we generated an assembly using Oxford Nanopore reads and Hi-C mapping. The assembly incorporates telomeric and centromeric regions, thereby establishing it as a near telomere-to-telomere genome except a region in chromosome 1 that is not fully assembled due to its highly repetitive nature. This near telomere-to-telomere genome resolves the chromosome number at 26 and provides a gap-free genome assembly as well as updated gene models to aid future studies using this model organism.

Genome-wide identification of the GRF family in sweet orange (Citrus sinensis) and functional analysis of the CsGRF04 in response to multiple abiotic stresses

BACKGROUND

Citrus is one of the most valuable fruits worldwide and an economic pillar industry in southern China. Nevertheless, it frequently suffers from undesirable environmental stresses during the growth cycle, which severely restricts the growth, development and yield of citrus. In plants, the growth-regulating factor (GRF) family of transcription factors (TF) is extensively distributed and plays an vital part in plant growth and development, hormone response, as well as stress adaptation. However, the systematic identification and functional analysis of GRF TFs in citrus have not been reported.

RESULTS

Here, a genome-wide identification of GRF TFs was performed in Citrus sinensis, 9 members of CsGRFs were systematically identified and discovered to be scattered throughout 5 chromosomes. Subsequently, physical and chemical properties, phylogenetic relationships, structural characteristics, gene duplication events, collinearity and cis-elements of promoter were elaborately analyzed. In particular, the expression patterns of the CsGRF genes in response to multiple phytohormone and abiotic stress treatments were investigated. Predicated on this result, CsGRF04, which exhibited the most differential expression pattern under multiple phytohormone and abiotic stress treatments was screened out. Virus-induced gene silencing (VIGS) technology was utilized to obtain gene silenced plants for CsGRF04 successfully. After the three stress treatments of high salinity, low temperature and drought, the CsGRF04-VIGS lines showed significantly reduced resistance to high salinity and low temperature stresses, but extremely increased resistance to drought stress.

CONCLUSIONS

Taken together, our findings systematically analyzed the genomic characterization of GRF family in Citrus sinensis, and excavated a CsGRF04 with potential functions under multiple abiotic stresses. Our study lay a foundation for further study on the function of CsGRFs in abiotic stress and hormone signaling response.

An enhancer trap system to track developmental dynamics in Marchantia polymorpha

A combination of streamlined genetics, experimental tractability and relative morphological simplicity compared to vascular plants makes the liverwort Marchantia polymorpha an ideal model system for studying many aspects of plant biology. Here we describe a transformation vector combining a constitutive fluorescent membrane marker with a nuclear marker that is regulated by nearby enhancer elements and use this to produce a library of enhancer trap lines for Marchantia. Screening gemmae from these lines allowed the identification and characterization of novel marker lines, including markers for rhizoids and oil cells. The library allowed the identification of a margin tissue running around the thallus edge, highlighted during thallus development. The expression of this marker is correlated with auxin levels. We generated multiple markers for the meristematic apical notch region, which have different spatial expression patterns, reappear at different times during meristem regeneration following apical notch excision and have varying responses to auxin supplementation or inhibition. This reveals that there are proximodistal substructures within the apical notch that could not be observed otherwise. We employed our markers to study Marchantia sporeling development, observing meristem emergence as defining the protonema-to-prothallus stage transition, and subsequent production of margin tissue during the prothallus stage. Exogenous auxin treatment stalls meristem emergence at the protonema stage but does not inhibit cell division, resulting in callus-like sporelings with many rhizoids, whereas pharmacologically inhibiting auxin synthesis and transport does not prevent meristem emergence. This enhancer trap system presents a useful resource for the community and will contribute to future Marchantia research.

An mTRAN-mRNA interaction mediates mitochondrial translation initiation in plants

Plant mitochondria represent the largest group of respiring organelles on the planet. Plant mitochondrial messenger RNAs (mRNAs) lack Shine-Dalgarno-like ribosome-binding sites, so it is unknown how plant mitoribosomes recognize mRNA. We show that "mitochondrial translation factors" mTRAN1 and mTRAN2 are land plant-specific proteins, required for normal mitochondrial respiration chain biogenesis. Our studies suggest that mTRANs are noncanonical pentatricopeptide repeat (PPR)-like RNA binding proteins of the mitoribosomal "small" subunit. We identified conserved Adenosine (A)/Uridine (U)-rich motifs in the 5' regions of plant mitochondrial mRNAs. mTRAN1 binds this motif, suggesting that it is a mitoribosome homing factor to identify mRNAs. We demonstrate that mTRANs are likely required for translation of all plant mitochondrial mRNAs. Plant mitochondrial translation initiation thus appears to use a protein-mRNA interaction that is divergent from bacteria or mammalian mitochondria.

CIPK-B is essential for salt stress signalling in Marchantia polymorpha

Calcium signalling is central to many plant processes, with families of calcium decoder proteins having expanded across the green lineage and redundancy existing between decoders. The liverwort Marchantia polymorpha has fast become a new model plant, but the calcium decoders that exist in this species remain unclear. We performed phylogenetic analyses to identify the calcineurin B-like (CBL) and CBL-interacting protein kinase (CIPK) network of M. polymorpha. We analysed CBL-CIPK expression during salt stress, and determined protein-protein interactions using yeast two-hybrid and bimolecular fluorescence complementation. We also created genetic knockouts using CRISPR/Cas9. We confirm that M. polymorpha has two CIPKs and three CBLs. Both CIPKs and one CBL show pronounced salt-responsive transcriptional changes. All M. polymorpha CBL-CIPKs interact with each other in planta. Knocking out CIPK-B causes increased sensitivity to salt, suggesting that this CIPK is involved in salt signalling. We have identified CBL-CIPKs that form part of a salt tolerance pathway in M. polymorpha. Phylogeny and interaction studies imply that these CBL-CIPKs form an evolutionarily conserved salt overly sensitive pathway. Hence, salt responses may be some of the early functions of CBL-CIPK networks and increased abiotic stress tolerance required for land plant emergence.

Crossroads in the evolution of plant specialized metabolism

The monophyletic group of embryophytes (land plants) stands out among photosynthetic eukaryotes: they are the sole constituents of the macroscopic flora on land. In their entirety, embryophytes account for the majority of the biomass on land and constitute an astounding biodiversity. What allowed for the massive radiation of this particular lineage? One of the defining features of all land plants is the production of an array of specialized metabolites. The compounds that the specialized metabolic pathways of embryophytes produce have diverse functions, ranging from superabundant structural polymers and compounds that ward off abiotic and biotic challenges, to signaling molecules whose abundance is measured at the nanomolar scale. These specialized metabolites govern the growth, development, and physiology of land plants-including their response to the environment. Hence, specialized metabolites define the biology of land plants as we know it. And they were likely a foundation for their success. It is thus intriguing to find that the closest algal relatives of land plants, freshwater organisms from the grade of streptophyte algae, possess homologs for key enzymes of specialized metabolic pathways known from land plants. Indeed, some studies suggest that signature metabolites emerging from these pathways can be found in streptophyte algae. Here we synthesize the current understanding of which routes of the specialized metabolism of embryophytes can be traced to a time before plants had conquered land.

Genome-wide identification, phylogeny, and expression analysis of GRF transcription factors in pineapple (Ananas comosus)

Background

Pineapple is the only commercially grown fruit crop in the Bromeliaceae family and has significant agricultural, industrial, economic, and ornamental value. GRF (growth-regulating factor) proteins are important transcription factors that have evolved in seed plants (embryophytes). They contain two conserved domains, QLQ (Gln, Leu, Gln) and WRC (Trp, Arg, Cys), and regulate multiple aspects of plant growth and stress response, including floral organ development, leaf growth, and hormone responses. The GRF family has been characterized in a number of plant species, but little is known about this family in pineapple and other bromeliads.

Main discoveries

We identified eight GRF transcription factor genes in pineapple, and phylogenetic analysis placed them into five subfamilies (I, III, IV, V, VI). Segmental duplication appeared to be the major contributor to expansion of the AcGRF family, and the family has undergone strong purifying selection during evolution. Relative to that of other gene families, the gene structure of the GRF family showed less conservation. Analysis of promoter cis-elements suggested that AcGRF genes are widely involved in plant growth and development. Transcriptome data and qRT-PCR results showed that, with the exception of AcGRF5, the AcGRFs were preferentially expressed in the early stage of floral organ development and AcGRF2 was strongly expressed in ovules. Gibberellin treatment significantly induced AcGRF7/8 expression, suggesting that these two genes may be involved in the molecular regulatory pathway by which gibberellin promotes pineapple fruit expansion.

Conclusion

AcGRF proteins appear to play a role in the regulation of floral organ development and the response to gibberellin. The information reported here provides a foundation for further study of the functions of AcGRF genes and the traits they regulate.

Comprehensive evolutionary analysis of the TCP gene family: Further insights for its origin, expansion, and diversification

TCP proteins are plant-specific transcription factors, which are involved in a broad range of physiological processes of plant growth and development. However, the origin and evolutionary history of this gene family is not fully resolved. Here, we present a genome-wide survey of TCP genes in 59 species (including 42 genomes and 17 transcriptomes) covering all main lineages of green plants, and reconstruct the evolutionary history of this gene family. Our results suggested that the origin of TCP genes predated the emergence of land plants, possibly in the common ancestor of Phragmoplastophyta. The TCP gene family gradually experienced a continuous expansion and grew from a few members in algae, moss and lycophytes to dozens, and sometimes over 50 members in angiosperms. Phylogenetic analysis indicated that at least four subclades (Class I and three subclades of Class II) have been occurred in the ancestor of spermatophyte (seed plant). Both dispersed duplication and segmental duplication or whole-genome duplication (WGD) contributed significantly to the expansion of the TCP gene family over the course of evolution. Our findings provide a comprehensive evolutionary analysis of the TCP gene family and highlight the importance of gene duplications in the evolution of this plant-specific transcription factors.

An empirical analysis of mtSSRs: could microsatellite distribution patterns explain the evolution of mitogenomes in plants?

Microsatellites (SSRs) are tandem repeat sequences in eukaryote genomes, including plant cytoplasmic genomes. The mitochondrial genome (mtDNA) has been shown to vary in size, number, and distribution of SSRs among different plant groups. Thus, SSRs contribute with genomic diversity in mtDNAs. However, the abundance, distribution, and evolutionary significance of SSRs in mtDNA from a wide range of algae and plants have not been explored. In this study, the mtDNAs of 204 plant and algal species were investigated related to the presence of SSRs. The number of SSRs was positively correlated with genome size. Its distribution is dependent on plant and algal groups analyzed, although the cluster analysis indicates the conservation of some common motifs in algal and terrestrial plants that reflect common ancestry of groups. Many SSRs in coding and non-coding regions can be useful for molecular markers. Moreover, mitochondrial SSRs are highly abundant, representing an important source for natural or induced genetic variation, i.e., for biotechnological approaches that can modulate mtDNA gene regulation. Thus, this comparative study increases the understanding of the plant and algal SSR evolution and brings perspectives for further studies.

Super determinant1A, a RAWULdomain-containing protein, modulates axillary meristem formation and compound leaf development in tomato

Shoot branching and complex leaf development relies on the establishment of boundaries that precedes the formation of axillary meristems (AMs) and leaflets. The tomato (Solanum lycopersicum) super determinant mutant is compromised in both processes, due to a mutation in Sde1A. Sde1A encodes a protein with a RAWUL domain, which is also present in Polycomb Group Repressive Complex 1 (PRC1) RING finger proteins and WD Repeat Domain 48 proteins. Genetic analysis revealed that Sde1A and Bmi1A cooperate, whereas Bmi1C antagonizes both activities, indicating the existence of functionally opposing PRC1 complexes that interact with Sde1A. Sde1A is expressed at early stages of boundary development in a small group of cells in the center of the leaf-axil boundary, but its activity is required for meristem formation at later stages. This suggests that Sde1A and Bmi1A promote AM formation and complex leaf development by safeguarding a pool of cells in the developing boundary zones. Genetic and protein interaction analyses showed that Sde1A and Lateral suppressor (Ls) are components of the same genetic pathway. In contrast to ls, sde1a mutants are not compromised in inflorescence branching, suggesting that Sde1A is a potential target for breeding tomato cultivars with reduced side-shoot formation during vegetative development.

Evolution of an assembly factor-based subunit contributed to a novel NDH-PSI supercomplex formation in chloroplasts

Chloroplast NADH dehydrogenase-like (NDH) complex is structurally related to mitochondrial Complex I and forms a supercomplex with two copies of Photosystem I (the NDH-PSI supercomplex) via linker proteins Lhca5 and Lhca6. The latter was acquired relatively recently in a common ancestor of angiosperms. Here we show that NDH-dependent Cyclic Electron Flow 5 (NDF5) is an NDH assembly factor in Arabidopsis. NDF5 initiates the assembly of NDH subunits (PnsB2 and PnsB3) and Lhca6, suggesting that they form a contact site with Lhca6. Our analysis of the NDF5 ortholog in Physcomitrella and angiosperm genomes reveals the subunit PnsB2 to be newly acquired via tandem gene duplication of NDF5 at some point in the evolution of angiosperms. Another Lhca6 contact subunit, PnsB3, has evolved from a protein unrelated to NDH. The structure of the largest photosynthetic electron transport chain complex has become more complicated by acquiring novel subunits and supercomplex formation with PSI.

The chloroplast NDH complex interacts with Photosystem I to form the NDH-PSI supercomplex. Here the authors show that Arabidopsis NDF5 shares a common ancestor with the NDH subunit PnsB2 and acts as an NDH assembly factor initiating the assembly of PnsB2 and the evolutionarily distinct PnsB3.

Gene-rich UV sex chromosomes harbor conserved regulators of sexual development

Nonrecombining sex chromosomes, like the mammalian Y, often lose genes and accumulate transposable elements, a process termed degeneration. The correlation between suppressed recombination and degeneration is clear in animal XY systems, but the absence of recombination is confounded with other asymmetries between the X and Y. In contrast, UV sex chromosomes, like those found in bryophytes, experience symmetrical population genetic conditions. Here, we generate nearly gapless female and male chromosome-scale reference genomes of the moss Ceratodon purpureus to test for degeneration in the bryophyte UV sex chromosomes. We show that the moss sex chromosomes evolved over 300 million years ago and expanded via two chromosomal fusions. Although the sex chromosomes exhibit weaker purifying selection than autosomes, we find that suppressed recombination alone is insufficient to drive degeneration. Instead, the U and V sex chromosomes harbor thousands of broadly expressed genes, including numerous key regulators of sexual development across land plants.

Carotenoid Biosynthesis in Liverworts

In higher plants, there are many studies on carotenoid biosynthetic pathways and their relevant genes. On the other hand, few researches exist on carotenoid biosynthesis in early-land plants containing liverworts, mosses, and ferns. Thus, the evolutionary history of carotenoid biosynthesis genes in land plants has remained unclear. A liverwort Marchantia polymorpha is thought to be one of the first land plants, since this plant remains a primitive figure. Moreover, this liverwort is regarded as the model plant of bryophytes due to several reasons. In this chapter, we review carotenoid biosynthesis in liverworts and discuss the functional evolution and evolutionary history of carotenogenic genes in land plants.

Antenna arrangement and energy-transfer pathways of PSI-LHCI from the moss Physcomitrella patens

Plants harvest light energy utilized for photosynthesis by light-harvesting complex I and II (LHCI and LHCII) surrounding photosystem I and II (PSI and PSII), respectively. During the evolution of green plants, moss is at an evolutionarily intermediate position from aquatic photosynthetic organisms to land plants, being the first photosynthetic organisms that landed. Here, we report the structure of the PSI–LHCI supercomplex from the moss Physcomitrella patens (Pp) at 3.23 Å resolution solved by cryo-electron microscopy. Our structure revealed that four Lhca subunits are associated with the PSI core in an order of Lhca1–Lhca5–Lhca2–Lhca3. This number is much decreased from 8 to 10, the number of subunits in most green algal PSI–LHCI, but the same as those of land plants. Although Pp PSI–LHCI has a similar structure as PSI–LHCI of land plants, it has Lhca5, instead of Lhca4, in the second position of Lhca, and several differences were found in the arrangement of chlorophylls among green algal, moss, and land plant PSI–LHCI. One chlorophyll, PsaF–Chl 305, which is found in the moss PSI–LHCI, is located at the gap region between the two middle Lhca subunits and the PSI core, and therefore may make the excitation energy transfer from LHCI to the core more efficient than that of land plants. On the other hand, energy-transfer paths at the two side Lhca subunits are relatively conserved. These results provide a structural basis for unravelling the mechanisms of light-energy harvesting and transfer in the moss PSI–LHCI, as well as important clues on the changes of PSI–LHCI after landing.

Deciphering OPDA Signaling Components in the Momilactone-Producing Moss Calohypnum plumiforme

Jasmonic acid (JA) and its biologically active form jasmonoyl-L-isoleucine (JA-Ile) regulate defense responses to various environmental stresses and developmental processes in plants. JA and JA-Ile are synthesized from α-linolenic acids derived from membrane lipids via 12-oxo-phytodienoic acid (OPDA). In the presence of JA-Ile, the COI1 receptor physically interacts with JAZ repressors, leading to their degradation, resulting in the transcription of JA-responsive genes by MYC transcription factors. Although the biosynthesis of JA-Ile is conserved in vascular plants, it is not recognized by COI1 in bryophytes and is not biologically active. In the liverwort Marchantia polymorpha, dinor-OPDA (dn-OPDA), a homolog of OPDA with two fewer carbons, and its isomer dn-iso-OPDA accumulate after wounding and are recognized by COI1 to activate downstream signaling. The moss Calohypnum plumiforme produces the antimicrobial-specialized metabolites, momilactones. It has been reported that JA and JA-Ile are not detected in C. plumiforme and that OPDA, but not JA, can induce momilactone accumulation and the expression of these biosynthetic genes, suggesting that OPDA or its derivative is a biologically active molecule in C. plumiforme that induces chemical defense. In the present study, we investigated the biological functions of OPDA and its derivatives in C. plumiforme. Searching for the components potentially involving oxylipin signaling from transcriptomic and genomic data revealed that two COI1, three JAZ, and two MYC genes were present. Quantification analyses revealed that OPDA and its isomer iso-OPDA accumulated in larger amounts than dn-OPDA and dn-iso-OPDA after wounding. Moreover, exogenously applied OPDA, dn-OPDA, or dn-iso-OPDA induced the transcription of JAZ genes. These results imply that OPDA, dn-OPDA, and/or their isomers potentially act as biologically active molecules to induce the signaling downstream of COI1-JAZ. Furthermore, co-immunoprecipitation analysis showed the physical interaction between JAZs and MYCs, indicating the functional conservation of JAZs in C. plumiforme with other plants. These results suggest that COI1-JAZ-MYC mediated signaling is conserved and functional in C. plumiforme.

DNA methylation mutants in Physcomitrella patens elucidate individual roles of CG and non-CG methylation in genome regulation

Cytosine (DNA) methylation in plants regulates the expression of genes and transposons. While methylation in plant genomes occurs at CG, CHG, and CHH sequence contexts, the comparative roles of the individual methylation contexts remain elusive. Here, we present Physcomitrella patens as the second plant system, besides Arabidopsis thaliana, with viable mutants with an essentially complete loss of methylation in the CG and non-CG contexts. In contrast to A. thaliana, P. patens has more robust CHH methylation, similar CG and CHG methylation levels, and minimal cross-talk between CG and non-CG methylation, making it possible to study context-specific effects independently. Our data found CHH methylation to act in redundancy with symmetric methylation in silencing transposons and to regulate the expression of CG/CHG-depleted transposons. Specific elimination of CG methylation did not dysregulate transposons or genes. In contrast, exclusive removal of non-CG methylation massively up-regulated transposons and genes. In addition, comparing two exclusively but equally CG- or CHG-methylated genomes, we show that CHG methylation acts as a greater transcriptional regulator than CG methylation. These results disentangle the transcriptional roles of CG and non-CG, as well as symmetric and asymmetric methylation in a plant genome, and point to the crucial role of non-CG methylation in genome regulation.

High gene space divergence contrasts with frozen vegetative architecture in the moss family Funariaceae

A new paradigm has slowly emerged regarding the diversification of bryophytes, with inferences from molecular data highlighting a dynamic evolution of their genome. However, comparative studies of expressed genes among closely related taxa is so far missing. Here we contrast the dimensions of the vegetative transcriptome of Funaria hygrometrica and Physcomitrium pyriforme against the genome of their relative, Physcomitrium (Physcomitrella) patens. These three species of Funariaceae share highly conserved vegetative bodies, and are partially sympatric, growing on mineral soil in mostly temperate regions. We analyzed the vegetative gametophytic transcriptome of F. hygrometrica and P. pyriforme and mapped short reads, transcripts, and proteins to the genome and gene space of P. patens. Only about half of the transcripts of F. hygrometrica map to their ortholog in P. patens, whereas at least 90% of those of P. pyriforme align to loci in P. patens. Such divergence is unexpected given the high morphological similarity of the gametophyte but reflects the estimated times of divergence of F. hygrometrica and P. pyriforme from P. patens, namely 55 and 20 mya, respectively. The newly sampled transcriptomes bear signatures of at least one, rather ancient, whole genome duplication (WGD), which may be shared with one reported for P. patens. The transcriptomes of F. hygrometrica and P. pyriforme reveal significant contractions or expansions of different gene families. While transcriptomes offer only an incomplete estimate of the gene space, the high number of transcripts obtained suggest a significant divergence in gene sequences, and gene number among the three species, indicative of a rather strong, dynamic genome evolution, shaped in part by whole, partial or localized genome duplication. The gene ontology of their specific and rapidly-evolving protein families, suggests that the evolution of the Funariaceae may have been driven by the diversification of metabolic genes that may optimize the adaptations to environmental conditions, a hypothesis well in line with ecological patterns in the genetic diversity and structure in seed plants.

Single Nucleotide Polymorphism Charting of P. patens Reveals Accumulation of Somatic Mutations During in vitro Culture on the Scale of Natural Variation by Selfing

Introduction: Physcomitrium patens (Hedw.) Mitten (previously known as Physcomitrella patens) was collected by H.L.K. Whitehouse in Gransden Wood (Huntingdonshire, United Kingdom) in 1962 and distributed across the globe starting in 1974. Hence, the Gransden accession has been cultured in vitro in laboratories for half a century. Today, there are more than 13 different pedigrees derived from the original accession. Additionally, accessions from other sites worldwide were collected during the last decades. Methods and Results: In this study, 250 high throughput RNA sequencing (RNA-seq) samples and 25 gDNA samples were used to detect single nucleotide polymorphisms (SNPs). Analyses were performed using five different P. patens accessions and 13 different Gransden pedigrees. SNPs were overlaid with metadata and known phenotypic variations. Unique SNPs defining Gransden pedigrees and accessions were identified and experimentally confirmed. They can be successfully employed for PCR-based identification. Conclusion: We show independent mutations in different Gransden laboratory pedigrees, demonstrating that somatic mutations occur and accumulate during in vitro culture. The frequency of such mutations is similar to those observed in naturally occurring populations. We present evidence that vegetative propagation leads to accumulation of deleterious mutations, and that sexual reproduction purges those. Unique SNP sets for five different P. patens accessions were isolated and can be used to determine individual accessions as well as Gransden pedigrees. Based on that, laboratory methods to easily determine P. patens accessions and Gransden pedigrees are presented.

The Ceratodon purpureus genome uncovers structurally complex, gene rich sex chromosomes.. [PMID: 0 DOI: 10.1101/2020.07.03.163634]

Non-recombining sex chromosomes, like the mammalian Y, often lose genes and accumulate transposable elements, a process termed degeneration1,2. The correlation between suppressed recombination and degeneration is clear in animal XY systems1,2, but the absence of recombination is confounded with other asymmetries between the X and Y. In contrast, UV sex chromosomes, like those found in bryophytes, experience symmetrical population genetic conditions3,4. Here we test for degeneration in the bryophyte UV sex chromosome system through genomic comparisons with new female and male chromosome-scale reference genomes of the moss Ceratodon purpureus. We show that the moss sex chromosomes evolved over 300 million years ago and expanded via two chromosomal fusions. Although the sex chromosomes show signs of weaker purifying selection than autosomes, we find suppressed recombination alone is insufficient to drive gene loss on sex-specific chromosomes. Instead, the U and V sex chromosomes harbor thousands of broadly-expressed genes, including numerous key regulators of sexual development across land plants.

The Moss Physcomitrium (Physcomitrella) patens: A Model Organism for Non-Seed Plants

Since the discovery two decades ago that transgenes are efficiently integrated into the genome of Physcomitrella patens by homologous recombination, this moss has been a premier model system to study evolutionary developmental biology questions, stem cell reprogramming, and the biology of nonvascular plants. P patens was the first non-seed plant to have its genome sequenced. With this level of genomic information, together with increasing molecular genetic tools, a large number of reverse genetic studies have propelled the use of this model system. A number of technological advances have recently opened the door to forward genetics as well as extremely efficient and precise genome editing in P patens Additionally, careful phylogenetic studies with increased resolution have suggested that P patens emerged from within Physcomitrium Thus, rather than Physcomitrella patens, the species should be named Physcomitrium patens Here we review these advances and describe the areas where P patens has had the most impact on plant biology.

Recurrent genome duplication events likely contributed to both the ancient and recent rise of ferns

Ferns, the second largest group of vascular plants, originated ~400 million years ago (Mya). They became dominant in the ancient Earth landscape before the angiosperms and are still important in current ecosystems. Many ferns have exceptionally high chromosome numbers, possibly resulting from whole-genome duplications (WGDs). However, WGDs have not been investigated molecularly across fern diversity. Here we detected and dated fern WGDs using a phylogenomic approach and by calculating synonymous substitution rates (Ks). We also investigated a possible correlation between proposed WGDs and shifts in species diversification rates. We identified 19 WGDs: three ancient events along the fern phylogenetic backbone that are shared by 66%-97% of extant ferns, with additional lineage-specific WGDs for eight orders, providing strong evidence for recurring genome duplications across fern evolutionary history. We also observed similar Ks peak values for more than half of these WGDs, with multiple WGDs occurring close to the Cretaceous (~145-66 Mya). Despite the repeated WGD events, the biodiversity of ferns declined during the Cretaceous, implying that other factors probably contributed to the floristic turnover from ferns to angiosperms. This study provides molecular evidence for recurring WGDs in ferns and offers important clues to the genomic evolutionary history of ferns.

Current Challenges in Studying Alternative Splicing in Plants: The Case of Physcomitrella patens SR Proteins

To colonize different terrestrial habitats, early land plants had to overcome the challenge of coping with harsh new environments. Alternative splicing - an RNA processing mechanism through which splice sites are differentially recognized, originating multiple transcripts and potentially different proteins from the same gene - can be key for plant stress tolerance. Serine/arginine-rich (SR) proteins constitute an evolutionarily conserved family of major alternative splicing regulators that in plants subdivides into six subfamilies. Despite being well studied in animals and a few plant species, such as the model angiosperm Arabidopsis thaliana and the crop Oryza sativa, little is known of these splicing factors in early land plants. Establishing the whole complement of SR proteins in different species is essential to understand the functional and evolutionary significance of alternative splicing. An in silico search for SR proteins in the extant moss Physcomitrella patens revealed inconsistencies both in the published data and available databases, likely arising from automatic annotation lacking adequate manual curation. These misannotations interfere with the description not only of the number and subfamily classification of Physcomitrella SR proteins but also of their domain architecture, potentially hindering the elucidation of their molecular functions. We therefore advise caution when looking into P. patens genomic resources. Our systematic survey nonetheless confidently identified 16 P. patens SR proteins that fall into the six described subfamilies and represent counterparts of well-established members in Arabidopsis and rice. Intensified research efforts should disclose whether SR proteins were already determining alternative splicing modulation and stress tolerance in early land plants.

Different Families of Retrotransposons and DNA Transposons Are Actively Transcribed and May Have Transposed Recently in Physcomitrium (Physcomitrella) patens

Similarly to other plant genomes of similar size, more than half of the genome of P. patens is covered by Transposable Elements (TEs). However, the composition and distribution of P. patens TEs is quite peculiar, with Long Terminal Repeat (LTR)-retrotransposons, which form patches of TE-rich regions interleaved with gene-rich regions, accounting for the vast majority of the TE space. We have already shown that RLG1, the most abundant TE in P. patens, is expressed in non-stressed protonema tissue. Here we present a non-targeted analysis of the TE expression based on RNA-Seq data and confirmed by qRT-PCR analyses that shows that, at least four LTR-RTs (RLG1, RLG2, RLC4 and tRLC5) and one DNA transposon (PpTc2) are expressed in P. patens. These TEs are expressed during development or under stresses that P. patens frequently faces, such as dehydratation/rehydratation stresses, suggesting that TEs have ample possibilities to transpose during P. patens life cycle. Indeed, an analysis of the TE polymorphisms among four different P. patens accessions shows that different TE families have recently transposed in this species and have generated genetic variability that may have phenotypic consequences, as a fraction of the TE polymorphisms are within or close to genes. Among the transcribed and mobile TEs, tRLC5 is particularly interesting as it concentrates in a single position per chromosome that could coincide with the centromere, and its expression is specifically induced in young sporophyte, where meiosis takes place.

Mosses in biotechnology

Biotechnological exploitation of mosses has several aspects, for example, the use of moss extracts or the whole plant for diverse industrial applications as well as their employment as production platforms for valuable metabolites or pharmaceutical proteins, especially using the genetically and developmentally best-characterised model moss Physcomitrella patens. Whole moss plants, in particular peat mosses (Sphagnum spec.), are useful for environmental approaches, biomonitoring of environmental pollution and CO₂-neutral 'farming' on rewetted bogs to combat climate change. In addition, the lifestyle of mosses suggests the evolution of genes necessary to cope with biotic and abiotic stress situations, which could be applied to crop plants, and their structural features bear an inspiring potential for biomimetics approaches.

A CRISPR/LbCas12a-based method for highly efficient multiplex gene editing in Physcomitrella patens

Due to their high efficiency, specificity, and flexibility, programmable nucleases, such as those of the clustered regularly interspaced short palindromic repeats (CRISPR)/Cas12a (Cpf1) system, have greatly expanded the applicability of editing the genomes of various organisms. Genes from different gene families or genes with redundant functions in the same gene family can be examined by assembling multiple CRISPR RNAs (crRNAs) in a single vector. However, the activity and efficiency of CRISPR/Cas12a in the non-vascular plant Physcomitrella patens are largely unknown. Here, we demonstrate that LbCas12a together with its mature crRNA can target multiple loci simultaneously in P. patens with high efficiency via co-delivery of LbCas12a and a crRNA expression cassette in vivo. The mutation frequencies induced by CRISPR/LbCas12a at a single locus ranged from 26.5 to 100%, with diverse deletions being the most common type of mutation. Our method expands the repertoire of genome editing tools available for P. patens and facilitates the creation of loss-of-function mutants of multiple genes from different gene families.

Mutations in Glucan, Water Dikinase Affect Starch Degradation and Gametophore Development in the Moss Physcomitrella patens

The role of starch degradation in non-vascular plants is poorly understood. To expand our knowledge of this area, we have studied this process in Physcomitrella patens. This has been achieved through examination of the step known to initiate starch degradation in angiosperms, glucan phosphorylation, catalysed by glucan, water dikinase (GWD) enzymes. Phylogenetic analysis indicates that GWD isoforms can be divided into two clades, one of which contains GWD1/GWD2 and the other GWD3 isoforms. These clades split at a very early stage within plant evolution, as distinct sequences that cluster within each were identified in all major plant lineages. Of the five genes we identified within the Physcomitrella genome that encode GWD-like enzymes, two group within the GWD1/GWD2 clade and the others within the GWD3 clade. Proteins encoded by both loci in the GWD1/GWD2 clade, named PpGWDa and PpGWDb, are localised in plastids. Mutations of either PpGWDa or PpGWDb reduce starch phosphate abundance, however, a mutation at the PpGWDa locus had a much greater influence than one at PpGWDb. Only mutations affecting PpGWDa inhibited starch degradation. Mutants lacking this enzyme also failed to develop gametophores, a phenotype that could be chemically complemented using glucose supplementation within the growth medium.

Distinct types of short open reading frames are translated in plant cells

Genomes contain millions of short (<100 codons) open reading frames (sORFs), which are usually dismissed during gene annotation. Nevertheless, peptides encoded by such sORFs can play important biological roles, and their impact on cellular processes has long been underestimated. Here, we analyzed approximately 70,000 transcribed sORFs in the model plant Physcomitrella patens (moss). Several distinct classes of sORFs that differ in terms of their position on transcripts and the level of evolutionary conservation are present in the moss genome. Over 5000 sORFs were conserved in at least one of 10 plant species examined. Mass spectrometry analysis of proteomic and peptidomic data sets suggested that tens of sORFs located on distinct parts of mRNAs and long noncoding RNAs (lncRNAs) are translated, including conserved sORFs. Translational analysis of the sORFs and main ORFs at a single locus suggested the existence of genes that code for multiple proteins and peptides with tissue-specific expression. Functional analysis of four lncRNA-encoded peptides showed that sORFs-encoded peptides are involved in regulation of growth and differentiation in moss. Knocking out lncRNA-encoded peptides resulted in a decrease of moss growth. In contrast, the overexpression of these peptides resulted in a diverse range of phenotypic effects. Our results thus open new avenues for discovering novel, biologically active peptides in the plant kingdom.

Evolution of TOR-SnRK dynamics in green plants and its integration with phytohormone signaling networks

The target of rapamycin (TOR)-sucrose non-fermenting 1 (SNF1)-related protein kinase 1 (SnRK1) signaling is an ancient regulatory mechanism that originated in eukaryotes to regulate nutrient-dependent growth. Although the TOR-SnRK1 signaling cascade shows highly conserved functions among eukaryotes, studies in the past two decades have identified many important plant-specific innovations in this pathway. Plants also possess SnRK2 and SnRK3 kinases, which originated from the ancient SnRK1-related kinases and have specialized roles in controlling growth, stress responses and nutrient homeostasis in plants. Recently, an integrative picture has started to emerge in which different SnRKs and TOR kinase are highly interconnected to control nutrient and stress responses of plants. Further, these kinases are intimately involved with phytohormone signaling networks that originated at different stages of plant evolution. In this review, we highlight the evolution and divergence of TOR-SnRK signaling components in plants and their communication with each other as well as phytohormone signaling to fine-tune growth and stress responses in plants.

Evaluation of the diversity and phylogenetic implications of NAC transcription factor members of four reference species from the different embryophytic plant groups

NAC transcription factors (TFs) are one of the largest and important TF family that are involved in the regulation of plant growth and development. They are characterized by a highly conserved N-terminal domain and a variable C-terminal domain. In the present study, the amino acid sequences of NAC TFs from four embryophytic plant species viz. Arabidopsis thaliana (Angiosperm), Picea abies (Gymnosperm), Selaginella moellendorffii (Pteridophyte) and Physcomitrella patens (Bryophyte) as reference of the different plant groups were collected from the Plant Transcription Factor Database (PTFD) and the phylogenetic relationships were evaluated. The phylogenetic tree revealed that the majority of the NAC members were interspersed in the major subgroups that indicated the expansion of the NAC members predates the speciation events. Thirty one (31), five (05), one (1) and ten (10) paralog pairs were determined respectively for Arabidopsis, Picea, Selaginella and Physcomitrella. The structure-function relationship of paralog pairs were inferred from the phylogenetic tree of combined set of paralogous gene pairs by studying the prevalence of flanking regions and motif analysis of the NAC proteins. The motif analysis revealed the presence of an N-terminal conserved domain, a characteristic of the majority of NAC family proteins. Conserved motifs in the C-terminal region were absent in the majority of the protein sequences except few members in Arabidopsis and Physcomitrella. Also the time of gene duplication of the paralog pairs were calculated that revealed the duplication events occurred between 4.48 and 45.94 MYA Arabidopsis, 167.57-532.86 MYA in Picea, and 29.12-53.53 MYA in Physcomitrella.

Abscisic acid, cold and salt stimulate conserved metabolic regulation in the moss Physcomitrella patens

Salt and cold are major abiotic stresses that have adverse effects on plant growth and development. To cope with these stresses and their detrimental effects plants have evolved several metabolic, biochemical and physiological processes that are mainly triggered and mediated by the plant hormone abscisic acid (ABA). To elucidate the metabolic responses of the moss Physcomitrella patens, which serves as a model plant for abiotic stress adaptation, we performed GC-MS-based metabolic profiling of plants challenged for 5 and 28 h with either salt, cold or ABA. Our results indicate significant changes in the accumulation of several sugars including maltose, isomaltose and trehalose, amino acids including arginine, histidine, ornithine, tryptophan and tyrosine, and organic acids mainly citric acid and malonic acid. The metabolic responses provoked by ABA, cold and salt show considerable similarities. The accumulation of certain metabolites positively correlates with gene expression data whereas some metabolites do not show correlation with cognate transcript abundance. To place our results into an evolutionary context we compared the ABA- and stress-induced metabolic changes in moss to available metabolic profiles of the seed plant Arabidopsis thaliana. We detected considerable conservation between the species, indicating early evolution of stress-associated metabolic adaptations that probably occurred at the plant water-to-land transition.

Assessing the Performance of Ks Plots for Detecting Ancient Whole Genome Duplications

Genomic data have provided evidence of previously unknown ancient whole genome duplications (WGDs) and highlighted the role of WGDs in the evolution of many eukaryotic lineages. Ancient WGDs often are detected by examining distributions of synonymous substitutions per site (Ks) within a genome, or “Ks plots.” For example, WGDs can be detected from Ks plots by using univariate mixture models to identify peaks in Ks distributions. We performed gene family simulation experiments to evaluate the effects of different Ks estimation methods and mixture models on our ability to detect ancient WGDs from Ks plots. The simulation experiments, which accounted for variation in substitution rates and gene duplication and loss rates across gene families, tested the effects of WGD age and gene retention rates following WGD on inferring WGDs from Ks plots. Our simulations reveal limitations of Ks plot analyses. Strict interpretations of mixture model analyses often overestimate the number of WGD events, and Ks plot analyses typically fail to detect WGDs when ≤10% of the duplicated genes are retained following the WGD. However, WGDs can accurately be characterized over an intermediate range of Ks. The simulation results are supported by empirical analyses of transcriptomic data, which also suggest that biases in gene retention likely affect our ability to detect ancient WGDs. Although our results indicate mixture model results should be interpreted with great caution, using node-averaged Ks estimates and applying more appropriate mixture models can improve the accuracy of detecting WGDs.

Physcomitrella patens, a versatile synthetic biology chassis

During three decades the moss Physcomitrella patens has been developed to a superb green cell factory with the first commercial products on the market. In the past three decades the moss P. patens has been developed from an obscure bryophyte to a model organism in basic biology, biotechnology, and synthetic biology. Some of the key features of this system include a wide range of Omics technologies, precise genome-engineering via homologous recombination with yeast-like efficiency, a certified good-manufacturing-practice production in bioreactors, successful upscaling to 500 L wave reactors, excellent homogeneity of protein products, superb product stability from batch-to-batch, and a reliable procedure for cryopreservation of cell lines in a master cell bank. About a dozen human proteins are being produced in P. patens as potential biopharmaceuticals, some of them are not only similar to their animal-produced counterparts, but are real biobetters with superior performance. A moss-made pharmaceutical successfully passed phase 1 clinical trials, a fragrant moss, and a cosmetic moss-product is already on the market, highlighting the economic potential of this synthetic biology chassis. Here, we focus on the features of mosses as versatile cell factories for synthetic biology and their impact on metabolic engineering.

Identification and Analysis of OVATE Family Members from Genome of the Early Land Plants Provide Insights into Evolutionary History of OFP Family and Function

Mosses, liverworts, hornworts and lycophytes represent transition stages between the aquatic to terrestrial/land plants. Several morphological and adaptive novelties driven by genomic components including emergence and expansion of new or existing gene families have played a critical role during and after the transition, and contributed towards successful colonization of terrestrial ecosystems. It is crucial to decipher the evolutionary transitions and natural selection on the gene structure and function to understand the emergence of phenotypic and adaptive diversity. Plants at the "transition zone", between aquatic and terrestrial ecosystem, are also the most vulnerable because of climate change and may contain clues for successful mitigation of the challenges of climate change. Identification and comparative analyses of such genetic elements and gene families are few in mosses, liverworts, hornworts and lycophytes. Ovate family proteins (OFPs) are plant-specific transcriptional repressors and are acknowledged for their roles in important growth and developmental processes in land plants, and information about the functional aspects of OFPs in early land plants is fragmentary. As a first step towards addressing this gap, a comprehensive in silico analysis was carried out utilizing publicly available genome sequences of Marchantia polymorpha (Mp), Physcomitrella patens (Pp), Selaginella moellendorffii (Sm) and Sphagnum fallax (Sf). Our analysis led to the identification of 4 MpOFPs, 19 PpOFPs, 6 SmOFPs and 3 SfOFPs. Cross-genera analysis revealed a drastic change in the structure and physiochemical properties in OFPs suggesting functional diversification and genomic plasticity during the evolutionary course. Knowledge gained from this comparative analysis will form the framework towards deciphering and dissection of their developmental and adaptive role/s in early land plants and could provide insights into evolutionary strategies adapted by land plants.

The single berberine bridge enzyme homolog of Physcomitrella patens is a cellobiose oxidase

The berberine bridge enzyme from the California poppy Eschscholzia californica (EcBBE) catalyzes the oxidative cyclization of (S)‐reticuline to (S)‐scoulerine, that is, the formation of the berberine bridge in the biosynthesis of benzylisoquinoline alkaloids. Interestingly, a large number of BBE‐like genes have been identified in plants that lack alkaloid biosynthesis. This finding raised the question of the primordial role of BBE in the plant kingdom, which prompted us to investigate the closest relative of EcBBE in Physcomitrella patens (PpBBE1), the most basal plant harboring a BBE‐like gene. Here, we report the biochemical, structural, and in vivo characterization of PpBBE1. Our studies revealed that PpBBE1 is structurally and biochemically very similar to EcBBE. In contrast to EcBBE, we found that PpBBE1 catalyzes the oxidation of the disaccharide cellobiose to the corresponding lactone, that is, PpBBE1 is a cellobiose oxidase. The enzymatic reaction mechanism was characterized by a structure‐guided mutagenesis approach that enabled us to assign a catalytic role to amino acid residues in the active site of PpBBE1. In vivo experiments revealed the highest level of PpBBE1 expression in chloronema, the earliest stage of the plant's life cycle, where carbon metabolism is strongly upregulated. It was also shown that the enzyme is secreted to the extracellular space, where it may be involved in later steps of cellulose degradation, thereby allowing the moss to make use of cellulose for energy production. Overall, our results suggest that the primordial role of BBE‐like enzymes in plants revolved around primary metabolic reactions in carbohydrate utilization.

Database

Structural data are available in the PDB under the accession numbers 6EO4 and 6EO5.

Evaluation of synthetic promoters in Physcomitrella patens

Securing a molecular toolbox including diverse promoters is essential for genome engineering. However, native promoters have limitations such as the available number or the length of the promoter. In this work, three short synthetic promoters were characterized by using the yellow fluorescent protein Venus. All of the tested promoters were active and showed higher mRNA expression than housekeeping gene PpAct7, and similar protein expression level to the AtUBQ10 promoter. This study shows that few cis-regulatory elements are enough to establish a strong promoter for continuous expression of genes in plants. Along with this, the study enhance the number of available promotors to be used in P. patens. It also demonstrates the potential to construct multiple non-native promoters on demand, which would aid to resolve one bottleneck in multiple pathway expression in P. patens and other plants.

RAB GTPases in the Basal Land Plant Marchantia polymorpha

The RAB GTPase is an evolutionarily conserved machinery component of membrane trafficking, which is the fundamental system for cell viability and higher order biological functions. The composition of RAB GTPases in each organism is closely related to the complexity and organization of the membrane trafficking pathway, which has been developed uniquely to realize the organism-specific membrane trafficking system. Comparative genomics has suggested that terrestrialization and/or multicellularization were associated with the expansion of membrane trafficking pathways in green plants, which has yet to be validated in basal land plant lineages. To obtain insight into the diversification of membrane trafficking systems in green plants, we analyzed RAB GTPases encoded in the genome of the liverwort Marchantia polymorpha in a comprehensive manner. We isolated all genes for RAB GTPases in Marchantia and analyzed their expression patterns and subcellular localizations in thallus cells. While a majority of MpRAB GTPases exhibited a ubiquitous expression pattern, specific exceptions were also observed; MpRAB2b, which contains a sequence similar to an intraflagellar transport protein at the C-terminal region; and MpRAB23, which has been secondarily lost in angiosperms, were specifically expressed in the male reproductive organ. MpRAB21, which is another RAB GTPase whose homolog is absent in Arabidopsis, exhibited endosomal localization with RAB5 members in Marchantia. These results suggest that Marchantia possesses unique membrane trafficking pathways involving a unique repertoire of RAB GTPases.

The Physcomitrella patens chromosome-scale assembly reveals moss genome structure and evolution

The draft genome of the moss model, Physcomitrella patens, comprised approximately 2000 unordered scaffolds. In order to enable analyses of genome structure and evolution we generated a chromosome-scale genome assembly using genetic linkage as well as (end) sequencing of long DNA fragments. We find that 57% of the genome comprises transposable elements (TEs), some of which may be actively transposing during the life cycle. Unlike in flowering plant genomes, gene- and TE-rich regions show an overall even distribution along the chromosomes. However, the chromosomes are mono-centric with peaks of a class of Copia elements potentially coinciding with centromeres. Gene body methylation is evident in 5.7% of the protein-coding genes, typically coinciding with low GC and low expression. Some giant virus insertions are transcriptionally active and might protect gametes from viral infection via siRNA mediated silencing. Structure-based detection methods show that the genome evolved via two rounds of whole genome duplications (WGDs), apparently common in mosses but not in liverworts and hornworts. Several hundred genes are present in colinear regions conserved since the last common ancestor of plants. These syntenic regions are enriched for functions related to plant-specific cell growth and tissue organization. The P. patens genome lacks the TE-rich pericentromeric and gene-rich distal regions typical for most flowering plant genomes. More non-seed plant genomes are needed to unravel how plant genomes evolve, and to understand whether the P. patens genome structure is typical for mosses or bryophytes.

Evolution of nuclear auxin signaling: lessons from genetic studies with basal land plants

Auxin plays critical roles in growth and development through the regulation of cell differentiation, cell expansion, and pattern formation. The auxin signal is mainly conveyed through a so-called nuclear auxin pathway involving the receptor TIR1/AFB, the transcriptional co-repressor AUX/IAA, and the transcription factor ARF with direct DNA-binding ability. Recent progress in sequence information and molecular genetics in basal plants has provided many insights into the evolutionary origin of the nuclear auxin pathway and its pleiotropic roles in land plant development. In this review, we summarize the latest knowledge of the nuclear auxin pathway gained from studies using basal plants, including charophycean green algae and two major model bryophytes, Marchantia polymorpha and Physcomitrella patens. In addition, we discuss the functional implication of the increase in genetic complexity of the nuclear auxin pathway during land plant evolution.

Genome-wide analysis of carbohydrate-active enzymes in Pyramimonas parkeae (Prasinophyceae)

The wall-less green flagellate Pyramimonas parkeae is classified in clade I of the prasinophytes, a paraphyletic assemblage representing the last common ancestor of Viridiplantae, a monophyletic group composed of the green algae and land plants. Consequently, P. parkeae and other prasinophytes illuminate early-evolved Viridiplantae traits likely fundamental in the systems biology of green algae and land plants. Cellular structure and organellar genomes of P. parkeae are now well understood, and transcriptomic sequence data are also publically available for one strain of this species, but corresponding nuclear genomic sequence data are lacking. For this reason, we obtained shotgun genomic sequence and assembled a draft nuclear genome for P. parkeaeNIES254 to use along with existing transcriptomic sequence to focus on carbohydrate-active enzymes. We found that the P. parkeae nuclear genome encodes carbohydrate-active protein families similar to those previously observed for other prasinophytes, green algae, and early-diverging embryophytes for which full nuclear genomic sequence is publically available. Sequences homologous to genes related to biosynthesis of starch and cell wall carbohydrates were identified in the P. parkeae genome, indicating molecular traits common to Viridiplantae. For example, the P. parkeae genome includes sequences clustering with bacterial genes that encode cellulose synthases (Bcs), including regions coding for domains common to bacterial and plant cellulose synthases; these new sequences were incorporated into phylogenies aimed at illuminating the evolutionary history of cellulose production by Viridiplantae. Genomic sequences related to biosynthesis of xyloglucans, pectin, and starch likewise shed light on the origin of key Viridiplantae traits.

Evolution of intrinsic disorder in eukaryotic proteins

Conformational flexibility conferred though regions of intrinsic structural disorder allows proteins to behave as dynamic molecules. While it is well-known that intrinsically disordered regions can undergo disorder-to-order transitions in real-time as part of their function, we also are beginning to learn more about the dynamics of disorder-to-order transitions along evolutionary time-scales. Intrinsically disordered regions endow proteins with functional promiscuity, which is further enhanced by the ability of some of these regions to undergo real-time disorder-to-order transitions. Disorder content affects gene retention after whole genome duplication, but it is not necessarily conserved. Altered patterns of disorder resulting from evolutionary disorder-to-order transitions indicate that disorder evolves to modify function through refining stability, regulation, and interactions. Here, we review the evolution of intrinsically disordered regions in eukaryotic proteins. We discuss the interplay between secondary structure and disorder on evolutionary time-scales, the importance of disorder for eukaryotic proteome expansion and functional divergence, and the evolutionary dynamics of disorder.

Diverse genome organization following 13 independent mesopolyploid events in Brassicaceae contrasts with convergent patterns of gene retention

Hybridization and polyploidy followed by genome-wide diploidization had a significant impact on the diversification of land plants. The ancient At-α whole-genome duplication (WGD) preceded the diversification of crucifers (Brassicaceae). Some genera and tribes also experienced younger, mesopolyploid WGDs concealed by subsequent genome diploidization. Here we tested if multiple base chromosome numbers originated due to genome diploidization after independent mesopolyploid WGDs and how diploidization affected post-polyploid gene retention. Sixteen species representing 10 Brassicaceae tribes were analyzed by comparative chromosome painting and/or whole-transcriptome analysis of gene age distributions and phylogenetic analyses of gene duplications. Overall, we found evidence for at least 13 independent mesopolyploidies followed by different degrees of diploidization across the Brassicaceae. New mesotetraploid events were uncovered for the tribes Anastaticeae, Iberideae and Schizopetaleae, and mesohexaploid WGDs for Cochlearieae and Physarieae. In contrast, we found convergent patterns of gene retention and loss among these independent WGDs. Our combined analyses of genomic data for Brassicaceae indicate that extant chromosome number variation in many plant groups, and especially monophyletic taxa with multiple base chromosome numbers, can result from clade-specific genome duplications followed by diploidization. Our observation of parallel gene retention and loss across multiple independent WGDs provides one of the first multi-species tests of the predictability of patterns of post-polyploid genome evolution.

Light-harvesting antenna complexes in the moss Physcomitrella patens: implications for the evolutionary transition from green algae to land plants

Plants have successfully adapted to a vast range of terrestrial environments during their evolution. To elucidate the evolutionary transition of light-harvesting antenna proteins from green algae to land plants, the moss Physcomitrella patens is ideally placed basally among land plants. Compared to the genomes of green algae and land plants, the P. patens genome codes for more diverse and redundant light-harvesting antenna proteins. It also encodes Lhcb9, which has characteristics not found in other light-harvesting antenna proteins. The unique complement of light-harvesting antenna proteins in P. patens appears to facilitate protein interactions that include those lost in both green algae and land plants with regard to stromal electron transport pathways and photoprotection mechanisms. This review will highlight unique characteristics of the P. patens light-harvesting antenna system and the resulting implications about the evolutionary transition during plant terrestrialization.

Sexual reproduction, sporophyte development and molecular variation in the model moss Physcomitrella patens: introducing the ecotype Reute

Rich ecotype collections are used for several plant models to unravel the molecular causes of phenotypic differences, and to investigate the effects of environmental adaption and acclimation. For the model moss Physcomitrella patens collections of accessions are available, and have been used for phylogenetic and taxonomic studies, for example, but few have been investigated further for phenotypic differences. Here, we focus on the Reute accession and provide expression profiling and comparative developmental data for several stages of sporophyte development, as well as information on genetic variation via genomic sequencing. We analysed cross-technology and cross-laboratory data to define a confident set of 15 mature sporophyte-specific genes. We find that the standard laboratory strain Gransden produces fewer sporophytes than Reute or Villersexel, although gametangia develop with the same time course and do not show evident morphological differences. Reute exhibits less genetic variation relative to Gransden than Villersexel, yet we found variation between Gransden and Reute in the expression profiles of several genes, as well as variation hot spots and genes that appear to evolve under positive Darwinian selection. We analyzed expression differences between the ecotypes for selected candidate genes in the GRAS transcription factor family, the chalcone synthase family and in genes involved in cell wall modification that are potentially related to phenotypic differences. We confirm that Reute is a P. patens ecotype, and suggest its use for reverse-genetics studies that involve progression through the life cycle and multiple generations.

Evolution of the 3R-MYB Gene Family in Plants

Plant 3R-MYB transcription factors are an important subgroup of the MYB super family in plants; however, their evolutionary history and functions remain poorly understood. We identified 225 3R-MYB proteins from 65 plant species, including algae and all major lineages of land plants. Two segmental duplication events preceding the common ancestor of angiosperms have given rise to three subgroups of the 3R-MYB proteins. Five conserved introns in the domain region of the 3R-MYB genes were identified, which arose through a step-wise pattern of intron gain during plant evolution. Alternative splicing (AS) analysis of selected species revealed that transcripts from more than 60% of 3R-MYB genes undergo AS. AS could regulate transcriptional activity for some of the plant 3R-MYBs by generating different regulatory motifs. The 3R-MYB genes of all subgroups appear to be enriched for Mitosis-Specific Activator element core sequences within their upstream promoter region, which suggests a functional involvement in cell cycle. Notably, expression of 3R-MYB genes from different species exhibits differential regulation under various abiotic stresses. These data suggest that the plant 3R-MYBs function in both cell cycle regulation and abiotic stress response, which may contribute to the adaptation of plants to a sessile lifestyle.

Characterization of Phytochrome Interacting Factors from the Moss Physcomitrella patens Illustrates Conservation of Phytochrome Signaling Modules in Land Plants

Across the plant kingdom, phytochrome (PHY) photoreceptors play an important role during adaptive and developmental responses to light. In Arabidopsis thaliana, light-activated PHYs accumulate in the nucleus, where they regulate downstream signaling components, such as phytochrome interacting factors (PIFs). PIFs are transcription factors that act as repressors of photomorphogenesis; their inhibition by PHYs leads to substantial changes in gene expression. The nuclear function of PHYs, however, has so far been investigated in only a few non-seed plants. Here, we identified putative target genes of PHY signaling in the moss Physcomitrella patens and found light-regulated genes that are putative orthologs of PIF-controlled genes in Arabidopsis. Phylogenetic analyses revealed that an ancestral PIF-like gene was already present in streptophyte algae, i.e., before the water-to-land transition of plants. The PIF homologs in the genome of P. patens resemble Arabidopsis PIFs in their protein domain structure, molecular properties, and physiological effects, albeit with notable differences in the motif-dependent PHY interaction. Our results suggest that P. patens PIFs are involved in PHY signaling. The PHY-PIF signaling node that relays light signals to target genes has been largely conserved during land plant evolution, with evidence of lineage-specific diversification.

Adaptation Mechanisms in the Evolution of Moss Defenses to Microbes

Bryophytes, including mosses, liverworts and hornworts are early land plants that have evolved key adaptation mechanisms to cope with abiotic stresses and microorganisms. Microbial symbioses facilitated plant colonization of land by enhancing nutrient uptake leading to improved plant growth and fitness. In addition, early land plants acquired novel defense mechanisms to protect plant tissues from pre-existing microbial pathogens. Due to its evolutionary stage linking unicellular green algae to vascular plants, the non-vascular moss Physcomitrella patens is an interesting organism to explore the adaptation mechanisms developed in the evolution of plant defenses to microbes. Cellular and biochemical approaches, gene expression profiles, and functional analysis of genes by targeted gene disruption have revealed that several defense mechanisms against microbial pathogens are conserved between mosses and flowering plants. P. patens perceives pathogen associated molecular patterns by plasma membrane receptor(s) and transduces the signal through a MAP kinase (MAPK) cascade leading to the activation of cell wall associated defenses and expression of genes that encode proteins with different roles in plant resistance. After pathogen assault, P. patens also activates the production of ROS, induces a HR-like reaction and increases levels of some hormones. Furthermore, alternative metabolic pathways are present in P. patens leading to the production of a distinct metabolic scenario than flowering plants that could contribute to defense. P. patens has acquired genes by horizontal transfer from prokaryotes and fungi, and some of them could represent adaptive benefits for resistance to biotic stress. In this review, the current knowledge related to the evolution of plant defense responses against pathogens will be discussed, focusing on the latest advances made in the model plant P. patens.

Plant Evo-Devo: How Tip Growth Evolved

Apical elongation of polarized plant cells (tip growth) occurs in root hairs of flowering plants and in rhizoids of bryophytes. A new report shows that the formation of these cells relies on genes already present in the first land plants.