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

Interactions between U and V sex chromosomes during the life cycle of Ectocarpus

In many animals and flowering plants, sex determination occurs in the diploid phase of the life cycle with XX/XY or ZW/ZZ sex chromosomes. However, in early diverging plants and most macroalgae, sex is determined by female (U) or male (V) sex chromosomes in a haploid phase called the gametophyte. Once the U and V chromosomes unite at fertilization to produce a diploid sporophyte, sex determination no longer occurs, raising key questions about the fate of the U and V sex chromosomes in the sporophyte phase. Here, we investigate genetic and molecular interactions of the UV sex chromosomes in both the haploid and diploid phases of the brown alga Ectocarpus. We reveal extensive developmental regulation of sex chromosome genes across its life cycle and implicate the TALE-HD transcription factor OUROBOROS in suppressing sex determination in the diploid phase. Small RNAs may also play a role in the repression of a female sex-linked gene, and transition to the diploid sporophyte coincides with major reconfiguration of histone H3K79me2, suggesting a more intricate role for this histone mark in Ectocarpus development than previously appreciated.

The SpRY Cas9 variant release the PAM sequence constraint for genome editing in the model plant Physcomitrium patens

Genome editing via CRISPR/Cas has enabled targeted genetic modifications in various species, including plants. The requirement for specific protospacer-adjacent motifs (PAMs) near the target gene, as seen with Cas nucleases like SpCas9, limits its application. PAMless SpCas9 variants, designed with a relaxed PAM requirement, have widened targeting options. However, these so-call PAMless SpCas9 still show variation of editing efficiency depending on the PAM and their efficiency lags behind the native SpCas9. Here we assess the potential of a PAMless SpCas9 variant for genome editing in the model plant Physcomitrium patens. For this purpose, we developed a SpRYCas9i variant, where expression was optimized, and tested its editing efficiency using the APT as a reporter gene. We show that the near PAMless SpRYCas9i effectively recognizes specific PAMs in P. patens that are not or poorly recognized by the native SpCas9. Pattern of mutations found using the SpRYCas9i are similar to the ones found with the SpCas9 and we could not detect off-target activity for the sgRNAs tested in this study. These findings contribute to advancing versatile genome editing techniques in plants.

Repeated co-option of HMG-box genes for sex determination in brown algae and animals

In many eukaryotes, genetic sex determination is not governed by XX/XY or ZW/ZZ systems but by a specialized region on the poorly studied U (female) or V (male) sex chromosomes. Previous studies have hinted at the existence of a dominant male-sex factor on the V chromosome in brown algae, a group of multicellular eukaryotes distantly related to animals and plants. The nature of this factor has remained elusive. Here, we demonstrate that an HMG-box gene acts as the male-determining factor in brown algae, mirroring the role HMG-box genes play in sex determination in animals. Over a billion-year evolutionary timeline, these lineages have independently co-opted the HMG box for male determination, representing a paradigm for evolution's ability to recurrently use the same genetic "toolkit" to accomplish similar tasks.

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

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

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.

Telomere Length Variation in Model Bryophytes

The ends of linear chromosomes of most eukaryotes consist of protein-bound DNA arrays called telomeres, which play essential roles in protecting genome integrity. Despite general evolutionary conservation in function, telomeric DNA is known to drastically vary in length and sequence between different eukaryotic lineages. Bryophytes are a group of early diverging land plants that include mosses, liverworts, and hornworts. This group of ancient land plants recently emerged as a new model for important discoveries in genomics and evolutionary biology, as well as for understanding plant adaptations to a terrestrial lifestyle. We measured telomere length in different ecotypes of model bryophyte species, including Physcomitrium patens, Marchantia polymorpha, Ceratodon purpureus, and in Sphagnum isolates. Our data indicate that all analyzed moss and liverwort genotypes have relatively short telomeres. Furthermore, all analyzed ecotypes and isolates of model mosses and liverworts display evidence of substantial natural variation in telomere length. Interestingly, telomere length also differs between male and female strains of the dioecious liverwort M. polymorpha and dioecious moss C. purpureus. Given that bryophytes are extraordinarily well adapted to different ecological niches from polar to tropical environments, our data will contribute to understanding the impact of natural telomere length variation on evolutionary adaptations in this ancient land plant lineage.

Improved prime editing allows for routine predictable gene editing in Physcomitrium patens

Efficient and precise gene editing is the gold standard of any reverse genetic study. The recently developed prime editing approach, a modified CRISPR/Cas9 [clustered regularly interspaced palindromic repeats (CRISPR)/CRISPR-associated protein] editing method, has reached the precision goal but its editing rate can be improved. We present an improved methodology that allows for routine prime editing in the model plant Physcomitrium patens, whilst exploring potential new prime editing improvements. Using a standardized protoplast transfection procedure, multiple prime editing guide RNA (pegRNA) structural and prime editor variants were evaluated targeting the APT reporter gene through direct plant selection. Together, enhancements of expression of the prime editor, modifications of the 3' extension of the pegRNA, and the addition of synonymous mutation in the reverse transcriptase template sequence of the pegRNA dramatically improve the editing rate without affecting the quality of the edits. Furthermore, we show that prime editing is amenable to edit a gene of interest through indirect selection, as demonstrated by the generation of a Ppdek10 mutant. Additionally, we determine that a plant retrotransposon reverse transcriptase enables prime editing. Finally, we show for the first time the possibility of performing prime editing with two independently coded peptides.

Protocol: an improved method for inducing sporophyte generation in the model moss Physcomitrium patens under nitrogen starvation

BACKGROUND

Land plants exhibit a haplodiplontic life cycle, whereby multicellular bodies develop in both the haploid and diploid generations. The early-diverging land plants, known as bryophytes, have a haploid-dominant life cycle, in which a short-lived multicellular body in the diploid generation, known as the sporophyte, develops on the maternal haploid gametophyte tissues. The moss Physcomitrium (Physcomitrella) patens has become one of the most powerful model systems in evolutionary plant developmental studies. To induce diploid sporophytes of P. patens, several protocols are implemented. One of the conventional approaches is to grow approximately one-month-old gametophores for another month on Jiffy-7 pellets made from the peat moss that is difficult to fully sterilize. A more efficient method to obtain all tissues throughout the life cycle should accelerate studies of P. patens.

RESULTS

Here, we investigated the effect of nitrogen conditions on the growth and development of P. patens. We provide an improved protocol for the sporophyte induction of P. patens using a BCD-based solid culture medium without Jiffy-7 pellets, based on the finding that the formation of gametangia and subsequent sporophytes is promoted by nitrogen-free growth conditions. The protocol consists of two steps; first, culture the protonemata and gametophores on nitrogen-rich medium under continuous light at 25 °C, and then transfer the gametophores onto nitrogen-free medium under short-day and at 15 °C for sporophyte induction. The protocol enables to shorten the induction period and reduce the culture space.

CONCLUSIONS

Our more efficient and shortened protocol for inducing the formation of sporophytes will contribute to future studies into the fertilization or the diploid sporophyte generation of P. patens.

DELLA proteins regulate spore germination and reproductive development in Physcomitrium patens

Proteins of the DELLA family integrate environmental signals to regulate growth and development throughout the plant kingdom. Plants expressing non-degradable DELLA proteins underpinned the development of high-yielding 'Green Revolution' dwarf crop varieties in the 1960s. In vascular plants, DELLAs are regulated by gibberellins, diterpenoid plant hormones. How DELLA protein function has changed during land plant evolution is not fully understood. We have examined the function and interactions of DELLA proteins in the moss Physcomitrium (Physcomitrella) patens, in the sister group of vascular plants (Bryophytes). PpDELLAs do not undergo the same regulation as flowering plant DELLAs. PpDELLAs are not degraded by diterpenes, do not interact with GID1 gibberellin receptor proteins and do not participate in responses to abiotic stress. PpDELLAs do share a function with vascular plant DELLAs during reproductive development. PpDELLAs also regulate spore germination. PpDELLAs interact with moss-specific photoreceptors although a function for PpDELLAs in light responses was not detected. PpDELLAs likely act as 'hubs' for transcriptional regulation similarly to their homologues across the plant kingdom. Taken together, these data demonstrate that PpDELLA proteins share some biological functions with DELLAs in flowering plants, but other DELLA functions and regulation evolved independently in both plant lineages.

Co-action of COP1, SPA and cryptochrome in light signal transduction and photomorphogenesis of the moss Physcomitrium patens

The Arabidopsis COP1/SPA ubiquitin ligase suppresses photomorphogenesis in darkness. In the light, photoreceptors inactivate COP1/SPA to allow a light response. While SPA genes are specific to the green lineage, COP1 also exists in humans. This raises the question of when in evolution plant COP1 acquired the need for SPA accessory proteins. We addressed this question by generating Physcomitrium Ppcop1 mutants and comparing their visible and molecular phenotypes with those of Physcomitrium Ppspa mutants. The phenotype of Ppcop1 nonuple mutants resembles that of Ppspa mutants. Most importantly, both mutants produce green chloroplasts in complete darkness. They also exhibit dwarfed gametophores, disturbed branching of protonemata and absent gravitropism. RNA-sequencing analysis indicates that both mutants undergo weak constitutive light signaling in darkness. PpCOP1 and PpSPA proteins form a complex and they interact via their WD repeat domains with the VP motif of the cryptochrome CCE domain in a blue light-dependent manner. This resembles the interaction of Arabidopsis SPA proteins with Arabidopsis CRY1, and is different from that with Arabidopsis CRY2. Taken together, the data indicate that PpCOP1 and PpSPA act together to regulate growth and development of Physcomitrium. However, in contrast to their Arabidopsis orthologs, PpCOP1 and PpSPA proteins execute only partial suppression of light signaling in darkness. Hence, additional repressors may exist that contribute to the repression of a light response in dark-exposed Physcomitrium.

A Physcomitrella PIN protein acts in spermatogenesis and sporophyte retention

The auxin efflux PIN-FORMED (PIN) proteins are conserved in all land plants and important players in plant development. In the moss Physcomitrella (Physcomitrium patens), three canonical PINs (PpPINA-C) are expressed in the leafy shoot (gametophore). PpPINA and PpPINB show functional activity in vegetative growth and sporophyte development. Here, we examined the role of PpPINC in the life cycle of Physcomitrella. We established reporter and knockout lines for PpPINC and analysed vegetative and reproductive tissues using microscopy and transcriptomic sequencing of moss gametangia. PpPINC is expressed in immature leaves, mature gametangia and during sporophyte development. The sperm cells (spermatozoids) of pinC knockout mutants exhibit increased motility and an altered flagella phenotype. Furthermore, the pinC mutants have a higher portion of differentially expressed genes related to spermatogenesis, increased fertility and an increased abortion rate of premeiotic sporophytes. Here, we show that PpPINC is important for spermatogenesis and sporophyte retention. We propose an evolutionary conserved way of polar growth during early moss embryo development and sporophyte attachment to the gametophore while suggesting the mechanical function in sporophyte retention of a ring structure, the Lorch ring.

Establishment and optimization of a new model organism to study early land plant evolution: Germination, cultivation and oospore variation of Chara braunii Gmelin, 1826

For studying land plant evolution, the establishment and optimization of model organisms representing streptophytic algae, sister to land plants, is essential. Long-term cultivation experiments with Chara braunii S276 were performed over 8 years, since 4 years (Nov. 2018) under constant conditions. Additionally, short-term experiments for optimization of culture conditions were performed with three strains of C. braunii (S276, NIES-1604 and Lausiger Teiche, LaT-2708). Germination success after application of sterilization agents, addition of gibberellic acid and under different incubation conditions with respect to pre-treatment, irradiance regime and substrate was investigated in order to develop protocols for generative cultivation of at least unialgal cultures. The resulting cultivation protocols for C. braunii S276, allowing maintenance of vegetative as well as generative cultures are presented in detail, including protocols for germination induction and growth of sterilized and unsterilized oospores.

An overview of bioinformatics, genomics, and transcriptomics resources for bryophytes

Bryophytes are useful models for the study of plant evolution, development, plant-fungal symbiosis, stress responses, and gametogenesis. Additionally, their dominant haploid gametophytic phase makes them great models for functional genomics research, allowing straightforward genome editing and gene knockout via CRISPR or homologous recombination. Until 2016, however, the only bryophyte genome sequence published was that of Physcomitrium patens. Throughout recent years, several other bryophyte genomes and transcriptome datasets became available, enabling better comparative genomics in evolutionary studies. The increase in the number of bryophyte genome and transcriptome resources available has yielded a plethora of annotations, databases, and bioinformatics tools to access the new data, which covers the large diversity of this clade and whose biology comprises features such as association with arbuscular mycorrhiza fungi, sex chromosomes, low gene redundancy, or loss of RNA editing genes for organellar transcripts. Here we provide a guide to resources available for bryophytes with regards to genome and transcriptome databases and bioinformatics tools.

Two fingerprinting sets for Humulus lupulus based on KASP and microsatellite markers

Verification of clonal identity of hop (Humulus lupulus L.) cultivars within breeding programs and germplasm collections is vital to conserving genetic resources. Accurate and economic DNA-based tools are needed in dioecious hop to confirm identity and parentage, neither of which can be reliably determined from morphological observations. In this study, we developed two fingerprinting sets for hop: a 9-SSR fingerprinting set containing high-core repeats that can be run in a single PCR reaction and a kompetitive allele specific PCR (KASP) assay of 25 single nucleotide polymorphisms (SNPs). The SSR set contains a sex-linked primer pair, HI-AGA7, that was used to genotype 629 hop accessions from the US Department of Agriculture (USDA) National Clonal Germplasm Repository (NCGR), the USDA Forage Seed and Cereal Research (FSCR), and the University of Nebraska-Lincoln (UNL) collections. The SSR set identified unique genotypes except for 89 sets of synonymous samples. These synonyms included: cultivars with different designations, the same cultivars from different sources, heat-treated clones, and clonal variants. Population structure analysis clustered accessions into wild North American (WNA) and cultivated groups. Diversity was slightly higher in the cultivated samples due to larger sample size. Parentage and sib-ship analyses were used to identify true-to-type cultivars. The HI-AGA7 marker generated two male- and nine female-specific alleles among the cultivated and WNA samples. The SSR and KASP fingerprinting sets were compared in 190 samples consisting of cultivated and WNA accession for their ability to confirm identity and assess diversity and population structure. The SSR fingerprinting set distinguished cultivars, selections and WNA accessions while the KASP assays were unable to distinguish the WNA samples and had lower diversity estimates than the SSR set. Both fingerprinting sets are valuable tools for identity confirmation and parentage analysis in hop for different purposes. The 9-SSR assay is cost efficient when genotyping a small number of wild and cultivated hop samples (<96) while the KASP assay is easy to interpret and cost efficient for genotyping a large number of cultivated samples (multiples of 96).

Prime Editing in the model plant Physcomitrium patens and its potential in the tetraploid potato

Since its discovery and first applications for genome editing in plants, the clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 technology has revolutionized plant research and precision crop breeding. Although the classical CRISPR-Cas9 system is a highly efficient tool for disruptive targeted mutagenesis, this system is mostly inefficient for the introduction of precise and predictable nucleotide substitutions. Recently, Prime Editing technology has been developed, allowing the simultaneous generation of nucleotide transitions and transversions but also short defined indels. In this study, we report on the successful use of Prime Editing in two plants of interest: the plant model Physcomitrium patens and the tetraploid and highly heterozygous potato (Solanum tuberosum). In both cases editing rates were lower than with other CRISPR-Cas9 based techniques, but we were able to successfully introduce nucleotide transversions into targeted genes, a unique feature of Prime Editing. Additionally, the analysis of potential off-target mutation sites in P. patens suggested very high targeting fidelity in this organism. The present work paves the way for the use Prime Editing in Physcomitrium patens and potato, however highlighting the limitations that need to be overcome for more efficient precision plant breeding.

A vertically transmitted amalgavirus is present in certain accessions of the bryophyte Physcomitrium patens

In the last few years, next-generation sequencing techniques have started to be used to identify new viruses infecting plants. This has allowed to rapidly increase our knowledge on viruses other than those causing symptoms in economically important crops. Here we used this approach to identify a virus infecting Physcomitrium patens that has the typical structure of the double-stranded RNA endogenous viruses of the Amalgaviridae family, which we named Physcomitrium patens amalgavirus 1, or PHPAV1. PHPAV1 is present only in certain accessions of P. patens, where its RNA can be detected throughout the cell cycle of the plant. Our analysis demonstrates that PHPAV1 can be vertically transmitted through both paternal and maternal germlines, in crosses between accessions that contain the virus with accessions that do not contain it. This work suggests that PHPAV1 can replicate in genomic backgrounds different from those that actually contain the virus and opens the door for future studies on virus-host coevolution.

The nuclear GUCT domain-containing DEAD-box RNA helicases govern gametophytic and sporophytic development in Physcomitrium patens

KEY MESSAGE

In Physcomitrium patens, PpRH1/PpRH2 are GUCT-domain-containing DEAD-BOX RNA helicases localize to the nucleus. They are implicated in cell and tissue development in all stages of the moss life cycle.

ABSTRACT

The DEAD-box-containing RNA helicase family encompasses a large and functionally important group of enzymes involved in cellular processes committed to the metabolism of RNA, including its transcription, processing, transport, translation and decay. Studies indicate this protein family has implied roles in plant vegetative and reproductive developmental processes as well as response to environmental stresses such has cold and high salinity. We focus here on a small conserved sub-group of GUCT domain-containing RNA helicase in the moss Physcomitrium patens. Phylogenetic analysis shows that RNA helicases containing the GUCT domain form a distinct conserved clade across the green lineage. In this clade, the P. patens genome possesses two closely related paralogues RNA helicases predicted to be nuclear, PpRH1 and PpRH2. Using in-locus gene fluorescent tagging we show that PpRH1 is localized to the nucleus in protonema. Analysis of PpRH1 and PpRH2 deletions, individually and together, indicates their potential roles in protonema, gametophore and sporophyte cellular and tissue development in P. patens. Additionally, the ultrastructural analysis of phyllid chloroplasts in Δrh2 and Δrh1/2 shows distinct starch granule accumulation under standard growth conditions associated with changes in photosynthetic activity parameters. We could not detect effects of either temperature or stress on protonema growth or PpRH1 and PpRH2 expression. Together, these results suggest that nuclear GUCT-containing RNA helicases play a role primarily in developmental processes directly or indirectly linked to photosynthesis activity in the moss P. patens.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s11103-021-01152-w.