Conquering new grounds: plant organellar C-to-U RNA editing factors can be functional in the plant cytosol.

Plant mitochondrial and chloroplast transcripts are subject to numerous events of specific cytidine-to-uridine (C-to-U) RNA editing to correct genetic information. Key protein factors for this process are specific RNA-binding pentatricopeptide repeat (PPR) proteins, which are encoded in the nucleus and post-translationally imported into the two endosymbiotic organelles. Despite hundreds of C-to-U editing sites in the plant organelles, no comparable editing has been found for nucleo-cytosolic mRNAs raising the question why plant RNA editing is restricted to chloroplasts and mitochondria. Here, we addressed this issue in the model moss Physcomitrium patens, where all PPR-type RNA editing factors comprise specific RNA-binding and cytidine deamination functionalities in single proteins. To explore whether organelle-type RNA editing can principally also take place in the plant cytosol, we expressed PPR56, PPR65 and PPR78, three editing factors recently shown to also function in a bacterial setup, together with cytosolic co-transcribed native targets in Physcomitrium. While we obtained unsatisfying results upon their constitutive expression, we found strong cytosolic RNA editing under hormone-inducible expression. Moreover, RNA-Seq analyses revealed varying numbers of up to more than 900 off-targets in other cytosolic transcripts. We conclude that PPR-mediated C-to-U RNA editing is not per se incompatible with the plant cytosol but that its limited target specificity has restricted its occurrence to the much less complex transcriptomes of mitochondria and chloroplast in the course of evolution.

Conservation of the moss RNA editing factor PPR78 despite the loss of its known cytidine-to-uridine editing sites is explained by a hidden extra target.

Cytidine (C)-to-uridine (U) RNA editing in plant organelles relies on specific RNA-binding pentatricopeptide repeat (PPR) proteins. In the moss Physcomitrium patens, all such RNA editing factors feature a C-terminal DYW domain that acts as the cytidine deaminase for C-to-U conversion. PPR78 of Physcomitrium targets 2 mitochondrial editing sites, cox1eU755SL and rps14eU137SL. Remarkably, the latter is edited to highly variable degrees in different mosses. Here, we aimed to unravel the coevolution of PPR78 and its 2 target sites in mosses. Heterologous complementation in a Physcomitrium knockout line revealed that the variable editing of rps14eU137SL depends on the PPR arrays of different PPR78 orthologues but not their C-terminal domains. Intriguingly, PPR78 has remained conserved despite the simultaneous loss of editing at both known targets among Hypnales (feather mosses), suggesting it serves an additional function. Using a recently established RNA editing assay in Escherichia coli, we confirmed site-specific RNA editing by PPR78 in the bacterium and identified 4 additional off-targets in the bacterial transcriptome. Based on conservation profiles, we predicted ccmFNeU1465RC as a candidate editing target of PPR78 in moss mitochondrial transcriptomes. We confirmed editing at this site in several mosses and verified that PPR78 targets ccmFNeU1465RC in the bacterial editing system, explaining the conservation and functional adaptation of PPR78 during moss evolution.

Beyond a PPR-RNA recognition code: Many aspects matter for the multi-targeting properties of RNA editing factor PPR56.

The mitochondrial C-to-U RNA editing factor PPR56 of the moss Physcomitrium patens is an RNA-binding pentatricopeptide repeat protein equipped with a terminal DYW-type cytidine deaminase domain. Transferred into Escherichia coli, PPR56 works faithfully on its two native RNA editing targets, nad3eU230SL and nad4eU272SL, and also converts cytidines into uridines at over 100 off-targets in the bacterial transcriptome. Accordingly, PPR56 is attractive for detailed mechanistic studies in the heterologous bacterial setup, allowing for scoring differential RNA editing activities of many target and protein variants in reasonable time. Here, we report (i) on the effects of numerous individual and combined PPR56 protein and target modifications, (ii) on the spectrum of off-target C-to-U editing in the bacterial background transcriptome for PPR56 and two variants engineered for target re-direction and (iii) on combinations of targets in tandem or separately at the 5'- and 3'-ends of large mRNAs. The latter experimentation finds enhancement of RNA editing at weak targets in many cases, including cox3eU290SF as a new candidate mitogenome target. We conclude that C-to-U RNA editing can be much enhanced by transcript features also outside the region ultimately targeted by PPRs of a plant editing factor, possibly facilitated by its enrichment or scanning along transcripts.

DYW cytidine deaminase domains have a long-range impact on RNA recognition by the PPR array of chimeric plant C-to-U RNA editing factors and strongly affect target selection.

The protein factors for the specific C-to-U RNA editing events in plant mitochondria and chloroplasts possess unique arrays of RNA-binding pentatricopeptide repeats (PPRs) linked to carboxy-terminal cytidine deaminase DYW domains via the extension motifs E1 and E2. The E1 and E2 motifs have distant similarities to tetratricopeptide repeats known to mediate protein-protein interactions but their precise function is unclear. Here, we investigate the tolerance of PPR56 and PPR65, two functionally characterized RNA editing factors of the moss Physcomitrium patens, for the creation of chimeras by variably replacing their C-terminal protein regions. Making use of a heterologous RNA editing assay system in Escherichia coli we find that heterologous DYW domains can strongly restrict or widen the spectrum of off-targets in the bacterial transcriptome for PPR56. Surprisingly, our data suggest that these changes are not only caused by the preference of a given heterologous DYW domain for the immediate sequence environment of the cytidine to be edited but also by a long-range impact on the nucleotide selectivity of the upstream PPRs.

Plant mitochondrial RNA editing factors can perform targeted C-to-U editing of nuclear transcripts in human cells.

RNA editing processes are strikingly different in animals and plants. Up to thousands of specific cytidines are converted into uridines in plant chloroplasts and mitochondria whereas up to millions of adenosines are converted into inosines in animal nucleo-cytosolic RNAs. It is unknown whether these two different RNA editing machineries are mutually incompatible. RNA-binding pentatricopeptide repeat (PPR) proteins are the key factors of plant organelle cytidine-to-uridine RNA editing. The complete absence of PPR mediated editing of cytosolic RNAs might be due to a yet unknown barrier that prevents its activity in the cytosol. Here, we transferred two plant mitochondrial PPR-type editing factors into human cell lines to explore whether they could operate in the nucleo-cytosolic environment. PPR56 and PPR65 not only faithfully edited their native, co-transcribed targets but also different sets of off-targets in the human background transcriptome. More than 900 of such off-targets with editing efficiencies up to 91%, largely explained by known PPR-RNA binding properties, were identified for PPR56. Engineering two crucial amino acid positions in its PPR array led to predictable shifts in target recognition. We conclude that plant PPR editing factors can operate in the entirely different genetic environment of the human nucleo-cytosol and can be intentionally re-engineered towards new targets.

Delayed differentiation of epidermal cells walls can underlie pedomorphosis in plants: the case of pedomorphic petals in the hummingbird-pollinated Caiophora hibiscifolia (Loasaceae, subfam. Loasoideae) species.

BACKGROUND: Understanding the relationship between macroevolutionary diversity and variation in organism development is an important goal of evolutionary biology. Variation in the morphology of several plant and animal lineages is attributed to pedomorphosis, a case of heterochrony, where an ancestral juvenile shape is retained in an adult descendant. Pedomorphosis facilitated morphological adaptation in different plant lineages, but its cellular and molecular basis needs further exploration. Plant development differs from animal development in that cells are enclosed by cell walls and do not migrate. Moreover, in many plant lineages, the differentiated epidermis of leaves, and leaf-derived structures, such as petals, limits organ growth. We, therefore, proposed that pedomorphosis in leaves, and in leaf-derived structures, results from delayed differentiation of epidermal cells with respect to reproductive maturity. This idea was explored for petal evolution, given the importance of corolla morphology for angiosperm reproductive success. RESULTS: By comparing cell morphology and transcriptional profiles between 5 mm flower buds and mature flowers of an entomophile and an ornitophile Loasoideae species (a lineage that experienced transitions from bee- to hummingbird-pollination), we show that evolution of pedomorphic petals of the ornithophile species likely involved delayed differentiation of epidermal cells with respect to flower maturity. We also found that developmental mechanisms other than pedomorphosis might have contributed to evolution of corolla morphology. CONCLUSIONS: Our results highlight a need for considering alternatives to the flower-centric perspective when studying the origin of variation in flower morphology, as this can be generated by developmental processes that are also shared with leaves.

Acetylation of conserved lysines fine-tunes mitochondrial malate dehydrogenase activity in land plants.

Plants need to rapidly and flexibly adjust their metabolism to changes of their immediate environment. Since this necessity results from the sessile lifestyle of land plants, key mechanisms for orchestrating central metabolic acclimation are likely to have evolved early. Here, we explore the role of lysine acetylation as a post-translational modification to directly modulate metabolic function. We generated a lysine acetylome of the moss Physcomitrium patens and identified 638 lysine acetylation sites, mostly found in mitochondrial and plastidial proteins. A comparison with available angiosperm data pinpointed lysine acetylation as a conserved regulatory strategy in land plants. Focusing on mitochondrial central metabolism, we functionally analyzed acetylation of mitochondrial malate dehydrogenase (mMDH), which acts as a hub of plant metabolic flexibility. In P. patens mMDH1, we detected a single acetylated lysine located next to one of the four acetylation sites detected in Arabidopsis thaliana mMDH1. We assessed the kinetic behavior of recombinant A. thaliana and P. patens mMDH1 with site-specifically incorporated acetyl-lysines. Acetylation of A. thaliana mMDH1 at K169, K170, and K334 decreases its oxaloacetate reduction activity, while acetylation of P. patens mMDH1 at K172 increases this activity. We found modulation of the malate oxidation activity only in A. thaliana mMDH1, where acetylation of K334 strongly activated it. Comparative homology modeling of MDH proteins revealed that evolutionarily conserved lysines serve as hotspots of acetylation. Our combined analyses indicate lysine acetylation as a common strategy to fine-tune the activity of central metabolic enzymes with likely impact on plant acclimation capacity.

DYW domain structures imply an unusual regulation principle in plant organellar RNA editing catalysis.

RNA editosomes selectively deaminate cytidines to uridines in plant organellar transcripts-mostly to restore protein functionality and consequently facilitate mitochondrial and chloroplast function. The RNA editosomal pentatricopeptide repeat proteins serve target RNA recognition, whereas the intensively studied DYW domain elicits catalysis. Here we present structures and functional data of a DYW domain in an inactive ground state and activated. DYW domains harbour a cytidine deaminase fold and a C-terminal DYW motif, with catalytic and structural zinc atoms, respectively. A conserved gating domain within the deaminase fold regulates the active site sterically and mechanistically in a process that we termed gated zinc shutter. Based on the structures, an autoinhibited ground state and its activation are cross-validated by RNA editing assays and differential scanning fluorimetry. We anticipate that, in vivo, the framework of an active plant RNA editosome triggers the release of DYW autoinhibition to ensure a controlled and coordinated cytidine deamination playing a key role in mitochondrial and chloroplast homeostasis.

One C-to-U RNA Editing Site and Two Independently Evolved Editing Factors: Testing Reciprocal Complementation with DYW-Type PPR Proteins from the Moss Physcomitrium (Physcomitrella) patens and the Flowering Plants Macadamia integrifolia and Arabidopsis.

Cytidine-to-uridine RNA editing is a posttranscriptional process in plant organelles, mediated by specific pentatricopeptide repeat (PPR) proteins. In angiosperms, hundreds of sites undergo RNA editing. By contrast, only 13 sites are edited in the moss Physcomitrium (Physcomitrella) patens Some are conserved between the two species, like the mitochondrial editing site nad5eU598RC. The PPR proteins assigned to this editing site are known in both species: the DYW-type PPR protein PPR79 in P. patens and the E+-type PPR protein CWM1 in Arabidopsis (Arabidopsis thaliana). CWM1 also edits sites ccmCeU463RC, ccmBeU428SL, and nad5eU609VV. Here, we reciprocally expressed the P. patens and Arabidopsis editing factors in the respective other genetic environment. Surprisingly, the P. patens editing factor edited all target sites when expressed in the Arabidopsis cwm1 mutant background, even when carboxy-terminally truncated. Conversely, neither Arabidopsis CWM1 nor CWM1-PPR79 chimeras restored editing in P. patens ppr79 knockout plants. A CWM1-like PPR protein from the early diverging angiosperm macadamia (Macadamia integrifolia) features a complete DYW domain and fully rescued editing of nad5eU598RC when expressed in P. patens. We conclude that (1) the independently evolved P. patens editing factor PPR79 faithfully operates in the more complex Arabidopsis editing system, (2) truncated PPR79 recruits catalytic DYW domains in trans when expressed in Arabidopsis, and (3) the macadamia CWM1-like protein retains the capacity to work in the less complex P. patens editing environment.

Single organelle function and organization as estimated from Arabidopsis mitochondrial proteomics.

Mitochondria host vital cellular functions, including oxidative phosphorylation and co-factor biosynthesis, which are reflected in their proteome. At the cellular level plant mitochondria are organized into hundreds of discrete functional entities, which undergo dynamic fission and fusion. It is the individual organelle that operates in the living cell, yet biochemical and physiological assessments have exclusively focused on the characteristics of large populations of mitochondria. Here, we explore the protein composition of an individual average plant mitochondrion to deduce principles of functional and structural organisation. We perform proteomics on purified mitochondria from cultured heterotrophic Arabidopsis cells with intensity-based absolute quantification and scale the dataset to the single organelle based on criteria that are justified by experimental evidence and theoretical considerations. We estimate that a total of 1.4 million protein molecules make up a single Arabidopsis mitochondrion on average. Copy numbers of the individual proteins span five orders of magnitude, ranging from >40 000 for Voltage-Dependent Anion Channel 1 to sub-stoichiometric copy numbers, i.e. less than a single copy per single mitochondrion, for several pentatricopeptide repeat proteins that modify mitochondrial transcripts. For our analysis, we consider the physical and chemical constraints of the single organelle and discuss prominent features of mitochondrial architecture, protein biogenesis, oxidative phosphorylation, metabolism, antioxidant defence, genome maintenance, gene expression, and dynamics. While assessing the limitations of our considerations, we exemplify how our understanding of biochemical function and structural organization of plant mitochondria can be connected in order to obtain global and specific insights into how organelles work.

Plant organellar RNA editing: what 30 years of research has revealed.

The central dogma in biology defines the flow of genetic information from DNA to RNA to protein. Accordingly, RNA molecules generally accurately follow the sequences of the genes from which they are transcribed. This rule is transgressed by RNA editing, which creates RNA products that differ from their DNA templates. Analyses of the RNA landscapes of terrestrial plants have indicated that RNA editing (in the form of C-U base transitions) is highly prevalent within organelles (that is, mitochondria and chloroplasts). Numerous C→U conversions (and in some plants also U→C) alter the coding sequences of many of the organellar transcripts and can also produce translatable mRNAs by creating AUG start sites or eliminating premature stop codons, or affect the RNA structure, influence splicing and alter the stability of RNAs. RNA-binding proteins are at the heart of post-transcriptional RNA expression. The C-to-U RNA editing process in plant mitochondria involves numerous nuclear-encoded factors, many of which have been identified as pentatricopeptide repeat (PPR) proteins that target editing sites in a sequence-specific manner. In this review we report on major discoveries on RNA editing in plant organelles, since it was first documented 30 years ago.

The Expansion and Diversification of Pentatricopeptide Repeat RNA-Editing Factors in Plants.

The RNA-binding pentatricopeptide repeat (PPR) family comprises hundreds to thousands of genes in most plants, but only a few dozen in algae, indicating massive gene expansions during land plant evolution. The nature and timing of these expansions has not been well defined due to the sparse sequence data available from early-diverging land plant lineages. In this study, we exploit the comprehensive OneKP datasets of over 1000 transcriptomes from diverse plants and algae toward establishing a clear picture of the evolution of this massive gene family, focusing on the proteins typically associated with RNA editing, which show the most spectacular variation in numbers and domain composition across the plant kingdom. We characterize over 2 250 000 PPR motifs in over 400 000 proteins. In lycophytes, polypod ferns, and hornworts, nearly 10% of expressed protein-coding genes encode putative PPR editing factors, whereas they are absent from algae and complex-thalloid liverworts. We show that rather than a single expansion, most land plant lineages with high numbers of editing factors have continued to generate novel sequence diversity. We identify sequence variations that imply functional differences between PPR proteins in seed plants versus non-seed plants and variations we propose to be linked to seed-plant-specific editing co-factors. Finally, using the sequence variations across the datasets, we develop a structural model of the catalytic DYW domain associated with C-to-U editing and identify a clade of unique DYW variants that are strong candidates as U-to-C RNA-editing factors, given their phylogenetic distribution and sequence characteristics.

Towards a plant model for enigmatic U-to-C RNA editing: the organelle genomes, transcriptomes, editomes and candidate RNA editing factors in the hornwort Anthoceros agrestis.

Hornworts are crucial to understand the phylogeny of early land plants. The emergence of 'reverse' U-to-C RNA editing accompanying the widespread C-to-U RNA editing in plant chloroplasts and mitochondria may be a molecular synapomorphy of a hornwort-tracheophyte clade. C-to-U RNA editing is well understood after identification of many editing factors in models like Arabidopsis thaliana and Physcomitrella patens, but there is no plant model yet to investigate U-to-C RNA editing. The hornwort Anthoceros agrestis is now emerging as such a model system. We report on the assembly and analyses of the A. agrestis chloroplast and mitochondrial genomes, their transcriptomes and editomes, and a large nuclear gene family encoding pentatricopeptide repeat (PPR) proteins likely acting as RNA editing factors. Both organelles in A. agrestis feature high amounts of RNA editing, with altogether > 1100 sites of C-to-U and 1300 sites of U-to-C editing. The nuclear genome reveals > 1400 genes for PPR proteins with variable carboxyterminal DYW domains. We observe significant variants of the 'classic' DYW domain, in the meantime confirmed as the cytidine deaminase for C-to-U editing, and discuss the first attractive candidates for reverse editing factors given their excellent matches to U-to-C editing targets according to the PPR-RNA binding code.

Plant-type pentatricopeptide repeat proteins with a DYW domain drive C-to-U RNA editing in Escherichia coli.

RNA editing converting cytidines into uridines is a hallmark of gene expression in land plant chloroplasts and mitochondria. Pentatricopeptide repeat (PPR) proteins have a key role in target recognition, but the functional editosome in the plant organelles has remained elusive. Here we show that individual Physcomitrella patens DYW-type PPR proteins alone can perform efficient C-to-U editing in Escherichia coli reproducing the moss mitochondrial editing. Single amino acid exchanges in the DYW domain abolish RNA editing, confirming it as the functional cytidine deaminase. The modification of RNA targets and the identification of numerous off-targets in the E. coli transcriptome reveal nucleotide identities critical for RNA recognition and cytidine conversion. The straightforward amenability of the new E. coli setup will accelerate future studies on RNA target recognition through PPRs, on the C-to-U editing deamination machinery and towards future establishment of transcript editing in other genetic systems.

Combination of the Endogenous lhcsr1 Promoter and Codon Usage Optimization Boosts Protein Expression in the Moss Physcomitrella patens.

The moss Physcomitrella patens is used both as an evo-devo model and biotechnological production system for metabolites and pharmaceuticals. Strong in vivo expression of genes of interest is important for production of recombinant proteins, e.g., selectable markers, fluorescent proteins, or enzymes. In this regard, the choice of the promoter sequence as well as codon usage optimization are two important inside factors to consider in order to obtain optimum protein accumulation level. To reliably quantify fluorescence, we transfected protoplasts with promoter:GFP fusion constructs and measured fluorescence intensity of living protoplasts in a plate reader system. We used the red fluorescent protein mCherry under 2x 35S promoter control as second reporter to normalize for different transfection efficiencies. We derived a novel endogenous promoter and compared deletion variants with exogenous promoters. We used different codon-adapted green fluorescent protein (GFP) genes to evaluate the influence of promoter choice and codon optimization on protein accumulation in P. patens, and show that the promoter of the gene of P. patens chlorophyll a/b binding protein lhcsr1 drives expression of GFP in protoplasts significantly (more than twofold) better than the commonly used 2x 35S promoter or the rice actin1 promoter. We identified a shortened 677 bp version of the lhcsr1 promoter that retains full activity in protoplasts. The codon optimized GFP yields significantly (more than twofold) stronger fluorescence signals and thus demonstrates that adjusting codon usage in P. patens can increase expression strength. In combination, new promotor and codon optimized GFP conveyed sixfold increased fluorescence signal.

A Single-Target Mitochondrial RNA Editing Factor of Funaria hygrometrica Can Fully Reconstitute RNA Editing at Two Sites in Physcomitrella patens.

Nuclear-encoded pentatricopeptide repeat (PPR) proteins are key factors for site-specific RNA editing, converting cytidines into uridines in plant mitochondria and chloroplasts. All editing factors in the model moss Physcomitrella patens have a C-terminal DYW domain with similarity to cytidine deaminase. However, numerous editing factors in flowering plants lack such a terminal DYW domain, questioning its immediate role in the pyrimidine base conversion process. Here we further investigate the Physcomitrella DYW-type PPR protein PPR_78, responsible for mitochondrial editing sites cox1eU755SL and rps14eU137SL. Complementation assays with truncated proteins demonstrate that the DYW domain is essential for full PPR_78 editing functionality. The DYW domain can be replaced, however, with its counterpart from another editing factor, PPR_79. The PPR_78 ortholog of the related moss Funaria hygrometrica fully complements the Physcomitrella mutant for editing at both sites, although the editing site in rps14 is lacking in Funaria. Editing factor orthologs in different taxa may thus retain editing capacity for multiple sites despite the absence of editing requirement.

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.

The Physarum polycephalum Genome Reveals Extensive Use of Prokaryotic Two-Component and Metazoan-Type Tyrosine Kinase Signaling.

Physarum polycephalum is a well-studied microbial eukaryote with unique experimental attributes relative to other experimental model organisms. It has a sophisticated life cycle with several distinct stages including amoebal, flagellated, and plasmodial cells. It is unusual in switching between open and closed mitosis according to specific life-cycle stages. Here we present the analysis of the genome of this enigmatic and important model organism and compare it with closely related species. The genome is littered with simple and complex repeats and the coding regions are frequently interrupted by introns with a mean size of 100 bases. Complemented with extensive transcriptome data, we define approximately 31,000 gene loci, providing unexpected insights into early eukaryote evolution. We describe extensive use of histidine kinase-based two-component systems and tyrosine kinase signaling, the presence of bacterial and plant type photoreceptors (phytochromes, cryptochrome, and phototropin) and of plant-type pentatricopeptide repeat proteins, as well as metabolic pathways, and a cell cycle control system typically found in more complex eukaryotes. Our analysis characterizes P. polycephalum as a prototypical eukaryote with features attributed to the last common ancestor of Amorphea, that is, the Amoebozoa and Opisthokonts. Specifically, the presence of tyrosine kinases in Acanthamoeba and Physarum as representatives of two distantly related subdivisions of Amoebozoa argues against the later emergence of tyrosine kinase signaling in the opisthokont lineage and also against the acquisition by horizontal gene transfer.

Molecular evidence for convergent evolution and allopolyploid speciation within the Physcomitrium-Physcomitrella species complex.

BACKGROUND: The moss Physcomitrella patens (Hedw.) Bruch & Schimp. is an important experimental model system for evolutionary-developmental studies. In order to shed light on the evolutionary history of Physcomitrella and related species within the Funariaceae, we analyzed the natural genetic diversity of the Physcomitrium-Physcomitrella species complex. RESULTS: Molecular analysis of the nuclear single copy gene BRK1 reveals that three Physcomitrium species feature larger genome sizes than Physcomitrella patens and encode two expressed BRK1 homeologs (polyploidization-derived paralogs), indicating that they may be allopolyploid hybrids. Phylogenetic analyses of BRK1 as well as microsatellite simple sequence repeat (SSR) data confirm a polyphyletic origin for three Physcomitrella lineages. Differences in the conservation of mitochondrial editing sites further support hybridization and cryptic speciation within the Physcomitrium-Physcomitrella species complex. CONCLUSIONS: We propose a revised classification of the previously described four subspecies of Physcomitrella patens into three distinct species, namely Physcomitrella patens, Physcomitrella readeri and Physcomitrella magdalenae. We argue that secondary reduction of sporophyte complexity in these species is due to the establishment of an ecological niche, namely spores resting in mud and possible spore dispersal by migratory birds. Besides the Physcomitrium-Physcomitrella species complex, the Funariaceae are host to their type species, Funaria hygrometrica, featuring a sporophyte morphology which is more complex. Their considerable developmental variation among closely related lineages and remarkable trait evolution render the Funariaceae an interesting group for evolutionary and genetic research.