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

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.

Frequent allopolyploidy with distant progenitors in the moss genera Physcomitrium and Entosthodon (Funariaceae) identified via subgenome phasing of targeted nuclear genes

Allopolyploids represent a new frontier in species discovery among embryophytes. Within mosses, allopolyploid discovery is challenged by low morphological complexity. The rapid expansion of sequencing approaches in addition to computational developments to identifying genome merger and whole-genome duplication using variation among nuclear loci representing homeologs has allowed for increased allopolyploid discovery among mosses. Here, we test a novel approach to phasing homeologs within loci and phasing loci across subgenomes, or subgenome assignment, called Homologizer, in the family Funariaceae. We confirm the intergeneric hybrid nature of Entosthodon hungaricus, and the allopolyploid origin of Physcomitrium eurystomum and one population of Physcomitrium collenchymatum. We also reveal that hybridization gave rise to Physcomitrium immersum, as well as to yet unrecognized lineages sharing the phenotype of Physcomitrium pyriforme and Physcomitrium sphaericum. Our findings demonstrate the utility of our approach when working with polyploid genomes, and its value in identifying progenitor species using target capture data.

Phylogenomic structure and speciation in an emerging model: the Sphagnum magellanicum complex (Bryophyta)

Sphagnum magellanicum is one of two Sphagnum species for which a reference-quality genome exists to facilitate research in ecological genomics. Phylogenetic and comparative genomic analyses were conducted based on resequencing data from 48 samples and RADseq analyses based on 187 samples. We report herein that there are four clades/species within the S. magellanicum complex in eastern North America and that the reference genome belongs to Sphagnum divinum. The species exhibit tens of thousands (RADseq) to millions (resequencing) of fixed nucleotide differences. Two species, however, referred to informally as S. diabolicum and S. magni because they have not been formally described, are differentiated by only 100 (RADseq) to 1000 (resequencing) of differences. Introgression among species in the complex is demonstrated using D-statistics and f₄ ratios. One ecologically important functional trait, tissue decomposability, which underlies peat (carbon) accumulation, does not differ between segregates in the S. magellanicum complex, although previous research showed that many closely related Sphagnum species have evolved differences in decomposability/carbon sequestration. Phylogenetic resolution and more accurate species delimitation in the S. magellanicum complex substantially increase the value of this group for studying the early evolutionary stages of climate adaptation and ecological evolution more broadly.

Red Light and Glucose Enhance Cytokinin-Mediated Bud Initial Formation in Physcomitrium patens

Growth and development of Physcomitrium patens is endogenously regulated by phytohormones such as auxin and cytokinin. Auxin induces the transition of chloronema to caulonema. This transition is also regulated by additional factors such as quantity and quality of light, carbon supply, and other phytohormones such as strigolactones and precursors of gibberrelic acid. On the other hand, cytokinins induce the formation of bud initials following caulonema differentiation. However, the influence of external factors such as light or nutrient supply on cytokinin-mediated bud initial formation has not been demonstrated in Physcomitrium patens. This study deals with the effect of light quality and nutrient supply on cytokinin-mediated bud initial formation. Bud initial formation has been observed in wild type plants in different light conditions such as white, red, and blue light in response to exogenously supplied cytokinin as well as glucose. In addition, budding assay has been demonstrated in the cry1a mutant of Physcomitrium in different light conditions. The results indicate that carbon supply and red light enhance the cytokinin response, while blue light inhibits this process in Physcomitrium.

Autopolyploidization affects transcript patterns and gene targeting frequencies in Physcomitrella

In Physcomitrella, whole-genome duplications affected the expression of about 3.7% of the protein-encoding genes, some of them relevant for DNA repair, resulting in a massively reduced gene-targeting frequency. Qualitative changes in gene expression after an autopolyploidization event, a pure duplication of the whole genome (WGD), might be relevant for a different regulation of molecular mechanisms between angiosperms growing in a life cycle with a dominant diploid sporophytic stage and the haploid-dominant mosses. Whereas angiosperms repair DNA double-strand breaks (DSB) preferentially via non-homologous end joining (NHEJ), in the moss Physcomitrella homologous recombination (HR) is the main DNA-DSB repair pathway. HR facilitates the precise integration of foreign DNA into the genome via gene targeting (GT). Here, we studied the influence of ploidy on gene expression patterns and GT efficiency in Physcomitrella using haploid plants and autodiploid plants, generated via an artificial WGD. Single cells (protoplasts) were transfected with a GT construct and material from different time-points after transfection was analysed by microarrays and SuperSAGE sequencing. In the SuperSAGE data, we detected 3.7% of the Physcomitrella genes as differentially expressed in response to the WGD event. Among the differentially expressed genes involved in DNA-DSB repair was an upregulated gene encoding the X-ray repair cross-complementing protein 4 (XRCC4), a key player in NHEJ. Analysing the GT efficiency, we observed that autodiploid plants were significantly GT suppressed (p < 0.001) attaining only one third of the expected GT rates. Hence, an alteration of global transcript patterns, including genes related to DNA repair, in autodiploid Physcomitrella plants correlated with a drastic suppression of HR.

Axenic in vitro cultivation of 19 peat moss (Sphagnum L.) species as a resource for basic biology, biotechnology, and paludiculture

Sphagnum farming can substitute peat with renewable biomass and thus help mitigate climate change. Large volumes of the required founder material can only be supplied sustainably by axenic cultivation in bioreactors. We established axenic in vitro cultures from sporophytes of 19 Sphagnum species collected in Austria, Germany, Latvia, the Netherlands, Russia, and Sweden: S. angustifolium, S. balticum, S. capillifolium, S. centrale, S. compactum, S. cuspidatum, S. fallax, S. fimbriatum, S. fuscum, S. lindbergii, S. medium/divinum, S. palustre, S. papillosum, S. rubellum, S. russowii, S. squarrosum, S. subnitens, S. subfulvum and S. warnstorfii. These species cover five of the six European Sphagnum subgenera; namely, Acutifolia, Cuspidata, Rigida, Sphagnum and Squarrosa. Their growth was measured in suspension cultures, whereas their ploidy was determined by flow cytometry and compared with the genome size of Physcomitrella patens. We identified haploid and diploid Sphagnum species, found that their cells are predominantly arrested in the G1 phase of the cell cycle, and did not find a correlation between plant productivity and ploidy. DNA barcoding was achieved by sequencing introns of the BRK1 genes. With this collection, high-quality founder material for diverse large-scale applications, but also for basic Sphagnum research, is available from the International Moss Stock Center.

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.

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 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.

Characterisation of evolutionarily conserved key players affecting eukaryotic flagellar motility and fertility using a moss model

Defects in flagella/cilia are often associated with infertility and disease. Motile male gametes (sperm cells) are an ancestral eukaryotic trait that has been lost in several lineages like flowering plants. Here, we made use of a phenotypic male fertility difference between two moss (Physcomitrella patens) ecotypes to explore spermatozoid function. We compare genetic and epigenetic variation as well as expression profiles between the Gransden and Reute ecotype to identify a set of candidate genes associated with moss male infertility. We generated a loss-of-function mutant of a coiled-coil domain containing 39 (ccdc39) gene that is part of the flagellar hydin network. Defects in mammal and algal homologues of this gene coincide with a loss of fertility, demonstrating the evolutionary conservation of flagellar function related to male fertility across kingdoms. The Ppccdc39 mutant resembles the Gransden phenotype in terms of male fertility. Potentially, several somatic (epi-)mutations occurred during prolonged vegetative propagation of Gransden, causing regulatory differences of for example the homeodomain transcription factor BELL1. Probably these somatic changes are causative for the observed male fertility defect. We propose that moss spermatozoids might be employed as an easily accessible system to study male infertility of humans and animals in terms of flagellar structure and movement.

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.

RecQ Helicases Function in Development, DNA Repair, and Gene Targeting in Physcomitrella patens

RecQ DNA helicases are genome surveillance proteins found in all kingdoms of life. They are characterized best in humans, as mutations in RecQ genes lead to developmental abnormalities and diseases. To better understand RecQ functions in plants we concentrated on Arabidopsis thaliana and Physcomitrella patens, the model species predominantly used for studies on DNA repair and gene targeting. Phylogenetic analysis of the six P. patens RecQ genes revealed their orthologs in humans and plants. Because Arabidopsis and P. patens differ in their RecQ4 and RecQ6 genes, reporter and deletion moss mutants were generated and gene functions studied in reciprocal cross-species and cross-kingdom approaches. Both proteins can be found in meristematic moss tissues, although at low levels and with distinct expression patterns. PpRecQ4 is involved in embryogenesis and in subsequent development as demonstrated by sterility of ΔPpRecQ4 mutants and by morphological aberrations. Additionally, ΔPpRecQ4 displays an increased sensitivity to DNA damages and an increased rate of gene targeting. Therefore, we conclude that PpRecQ4 acts as a repressor of recombination. In contrast, PpRecQ6 is not obviously important for moss development or DNA repair but does function as a potent enhancer of gene targeting.

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.

Physcomitrella MADS-box genes regulate water supply and sperm movement for fertilization

MIKC classic (MIKC^C)-type MADS-box genes encode transcription factors that function in various developmental processes, including angiosperm floral organ identity. Phylogenetic analyses of the MIKC^C-type MADS-box family, including genes from non-flowering plants, suggest that the increased numbers of these genes in flowering plants is related to their functional divergence; however, their precise functions in non-flowering plants and their evolution throughout land plant diversification are unknown. Here, we show that MIKC^C-type MADS-box genes in the moss Physcomitrella patens function in two ways to enable fertilization. Analyses of protein localization, deletion mutants and overexpression lines of all six genes indicate that three MIKC^C-type MADS-box genes redundantly regulate cell division and growth in the stems for appropriate external water conduction, as well as the formation of sperm with motile flagella. The former function appears to be maintained in the flowering plant lineage, while the latter was lost in accordance with the loss of sperm.

Editing of Chloroplast rps14 by PPR Editing Factor EMB2261 Is Essential for Arabidopsis Development

RNA editing in plastids is known to be required for embryogenesis, but no single editing event had been shown to be essential. We show that the emb2261-2 mutation is lethal through a failure to express an editing factor that specifically recognizes the rps14-2 site. EMB2261 was predicted to bind the cis-element upstream of the rps14-2 site and genetic complementation with promoters of different strength followed by RNA-seq analysis was conducted to test the correlation between rps14-2 editing and EMB2261 expression. Rps14-2 is the only editing event in Arabidopsis chloroplasts that correlates with EMB2261 expression. Sequence divergence between the cis-element and the EMB2261 protein sequence in plants where rps14-2 editing is not required adds support to the association between them. We conclude that EMB2261 is the specificity factor for rps14-2 editing. This editing event converts P51 in Rps14 to L51, which is conserved among species lacking RNA editing, implying the importance of the editing event to Rps14 function. Rps14 is an essential ribosomal subunit for plastid translation, which, in turn, is essential for Arabidopsis embryogenesis.

Evolutionary dynamism in bryophytes: Phylogenomic inferences confirm rapid radiation in the moss family Funariaceae

Rapid diversifications of plants are primarily documented and studied in angiosperms, which are perceived as evolutionarily dynamic. Recent studies have, however, revealed that bryophytes have also undergone periods of rapid radiation. The speciose family Funariaceae, including the model taxon Physcomitrella patens, is one such lineage. Here, we infer relationships among major lineages within the Entosthodon-Physcomitrium complex from virtually complete organellar exomes (i.e., 123 genes) obtained through high throughput sequencing of genomic libraries enriched in these loci via targeted locus capture. Based on these extensive exonic data we (1) reconstructed a robust backbone topology of the Funariaceae, (2) confirmed the monophyly of Funaria and the polyphyly of Entosthodon, Physcomitrella, and Physcomitrium, and (3) argue for the occurrence of a rapid radiation within the Entosthodon-Physcomitrium complex that began 28 mya and gave rise more than half of the species diversity of the family. This diversification may have been triggered by a whole genome duplication and coincides with global Eocene cooling that continued through the Oligocene and Miocene. The Funariaceae join a growing list of bryophyte lineages whose history is marked by at least one burst of diversification, and our study thereby strengthens the view that bryophytes are evolutionarily dynamic lineages and that patterns and processes characterizing the evolution of angiosperms may be universal among land plants.

Polyploidy and interspecific hybridization: partners for adaptation, speciation and evolution in plants

BACKGROUND

Polyploidy or whole-genome duplication is now recognized as being present in almost all lineages of higher plants, with multiple rounds of polyploidy occurring in most extant species. The ancient evolutionary events have been identified through genome sequence analysis, while recent hybridization events are found in about half of the world's crops and wild species. Building from this new paradigm for understanding plant evolution, the papers in this Special Issue address questions about polyploidy in ecology, adaptation, reproduction and speciation of wild and cultivated plants from diverse ecosystems. Other papers, including this review, consider genomic aspects of polyploidy.

APPROACHES

Discovery of the evolutionary consequences of new, evolutionarily recent and ancient polyploidy requires a range of approaches. Large-scale studies of both single species and whole ecosystems, with hundreds to tens of thousands of individuals, sometimes involving 'garden' or transplant experiments, are important for studying adaptation. Molecular studies of genomes are needed to measure diversity in genotypes, showing ancestors, the nature and number of polyploidy and backcross events that have occurred, and allowing analysis of gene expression and transposable element activation. Speciation events and the impact of reticulate evolution require comprehensive phylogenetic analyses and can be assisted by resynthesis of hybrids. In this Special Issue, we include studies ranging in scope from experimental and genomic, through ecological to more theoretical.

CONCLUSIONS

The success of polyploidy, displacing the diploid ancestors of almost all plants, is well illustrated by the huge angiosperm diversity that is assumed to originate from recurrent polyploidization events. Strikingly, polyploidization often occurred prior to or simultaneously with major evolutionary transitions and adaptive radiation of species, supporting the concept that polyploidy plays a predominant role in bursts of adaptive speciation. Polyploidy results in immediate genetic redundancy and represents, with the emergence of new gene functions, an important source of novelty. Along with recombination, gene mutation, transposon activity and chromosomal rearrangement, polyploidy and whole-genome duplication act as drivers of evolution and divergence in plant behaviour and gene function, enabling diversification, speciation and hence plant evolution.

Selfing in Haploid Plants and Efficacy of Selection: Codon Usage Bias in the Model Moss Physcomitrella patens

A long-term reduction in effective population size will lead to major shift in genome evolution. In particular, when effective population size is small, genetic drift becomes dominant over natural selection. The onset of self-fertilization is one evolutionary event considerably reducing effective size of populations. Theory predicts that this reduction should be more dramatic in organisms capable for haploid than for diploid selfing. Although theoretically well-grounded, this assertion received mixed experimental support. Here, we test this hypothesis by analyzing synonymous codon usage bias of genes in the model moss Physcomitrella patens frequently undergoing haploid selfing. In line with population genetic theory, we found that the effect of natural selection on synonymous codon usage bias is very weak. Our conclusion is supported by four independent lines of evidence: 1) Very weak or nonsignificant correlation between gene expression and codon usage bias, 2) no increased codon usage bias in more broadly expressed genes, 3) no evidence that codon usage bias would constrain synonymous and nonsynonymous divergence, and 4) predominant role of genetic drift on synonymous codon usage predicted by a model-based analysis. These findings show striking similarity to those observed in AT-rich genomes with weak selection for optimal codon usage and GC content overall. Our finding is in contrast to a previous study reporting adaptive codon usage bias in the moss P. patens.

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.

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.

Comparative Cuticle Development Reveals Taller Sporophytes Are Covered by Thicker Calyptra Cuticles in Mosses

The calyptra is a maternal structure that protects the sporophyte offspring from dehydration, and positively impacts sporophyte survival and fitness in mosses. We explore the relationship between cuticle protection and sporophyte height as a proxy for dehydration stress in Funariaceae species with sporophytes across a range of sizes. Calyptrae and sporophytes from four species were collected from laboratory-grown populations at two developmental stages. Tissues were embedded, sectioned, and examined using transmission electron microscopy. Cuticle thickness was measured from three epidermal cells per organ for each individual and compared statistically. All four species have cuticles consisting of a cuticle proper and a cuticular layer on the calyptra and sporophyte at both developmental stages. Across species, shorter sporophytes are associated with smaller calyptra and thinner calyptra cuticles, whereas taller sporophytes are associated with larger calyptra and thicker calyptra cuticles. Independent of size, young sporophytes have a thin cuticle that thickens later during development, while calyptrae have a mature cuticle produced early during development that persists throughout development. This study adds to our knowledge of maternal effects influencing offspring survival in plants. Released from the pressures to invest in protection for their sporophyte offspring, maternal resources can be allocated to other processes that support sporophyte reproductive success. Using a comparative developmental framework enables us to broaden our understanding of cuticle development across species and provides structural evidence supporting the waterproofing role of the moss calyptra.

Endemism in the moss flora of North America

PREMISE OF THE STUDY

Identifying regions of high endemism is a critical step toward understanding the mechanisms underlying diversification and establishing conservation priorities. Here, we identified regions of high moss endemism across North America. We also identified lineages that contribute disproportionately to endemism and document the progress of efforts to inventory the endemic flora.

METHODS

To understand the documentation of endemic moss diversity in North America, we tabulated species publication dates to document the progress of species discovery across the continent. We analyzed herbarium specimen data and distribution data from the Flora of North America project to delineate major regions of moss endemism. Finally, we surveyed the literature to assess the importance of intercontinental vs. within-continent diversification for generating endemic species.

KEY RESULTS

Three primary regions of endemism were identified and two of these were further divided into a total of nine subregions. Overall endemic richness has two peaks, one in northern California and the Pacific Northwest, and the other in the southern Appalachians. Description of new endemic species has risen steeply over the last few decades, especially in western North America. Among the few studies documenting sister species relationships of endemics, recent diversification appears to have played a larger role in western North America, than in the east.

CONCLUSIONS

Our understanding of bryophyte endemism continues to grow rapidly. Large continent-wide data sets confirm early views on hotspots of endemic bryophyte richness and indicate a high rate of ongoing species discovery in North America.

Genomic scanning using AFLP to detect loci under selection in the moss Funaria hygrometrica along a climate gradient in the Sierra Nevada Mountains, Spain

The common cord moss Funaria hygrometrica has a worldwide distribution and thrives in a wide variety of environments. Here, we studied the genetic diversity in F. hygrometrica along an abiotic gradient in the Mediterranean high mountain of Sierra Nevada (Spain) using a genome scan method. Eighty-four samples from 17 locations from 24 to 2700 m were fingerprinted based on their amplified fragment length polymorphism (AFLP) banding pattern. Using PCA and Bayesian inference we found that the genetic diversity was structured in three or four clusters, respectively. Using a genome scan method we identified 13 outlier loci, which showed a signature of positive selection. Partial Mantel tests were performed between the Euclidean distance matrices of geographic and climatic variables, versus the pair-wise genetic distance of the AFLP dataset and AFLP-positive outliers dataset. AFLP-positive outlier data were significantly correlated with the gradient of the climatic variables, suggesting adaptive variation among populations of F. hygrometrica along the Sierra Nevada Mountains. We highlight the additional analyses necessary to identify the nature of these loci, and their biological role in the adaptation process.