Fern cell walls and the evolution of arabinogalactan proteins in streptophytes.

Significant changes have occurred in plant cell wall composition during evolution and diversification of tracheophytes. As the sister lineage to seed plants, knowledge on the cell wall of ferns is key to track evolutionary changes across tracheophytes and to understand seed plant-specific evolutionary innovations. Fern cell wall composition is not fully understood, including limited knowledge of glycoproteins such as the fern arabinogalactan proteins (AGPs). Here, we characterize the AGPs from the leptosporangiate fern genera Azolla, Salvinia, and Ceratopteris. The carbohydrate moiety of seed plant AGPs consists of a galactan backbone including mainly 1,3- and 1,3,6-linked pyranosidic galactose, which is conserved across the investigated fern AGPs. Yet, unlike AGPs of angiosperms, those of ferns contained the unusual sugar 3-O-methylrhamnose. Besides terminal furanosidic arabinose, Ara (Araf), the main linkage type of Araf in the ferns was 1,2-linked Araf, whereas in seed plants 1,5-linked Araf is often dominating. Antibodies directed against carbohydrate epitopes of AGPs supported the structural differences between AGPs of ferns and seed plants. Comparison of AGP linkage types across the streptophyte lineage showed that angiosperms have rather conserved monosaccharide linkage types; by contrast bryophytes, ferns, and gymnosperms showed more variability. Phylogenetic analyses of glycosyltransferases involved in AGP biosynthesis and bioinformatic search for AGP protein backbones revealed a versatile genetic toolkit for AGP complexity in ferns. Our data reveal important differences across AGP diversity of which the functional significance is unknown. This diversity sheds light on the evolution of the hallmark feature of tracheophytes: their elaborate cell walls.

An optimized transformation protocol for Anthoceros agrestis and three more hornwort species.

Land plants comprise two large monophyletic lineages, the vascular plants and the bryophytes, which diverged from their most recent common ancestor approximately 480 million years ago. Of the three lineages of bryophytes, only the mosses and the liverworts are systematically investigated, while the hornworts are understudied. Despite their importance for understanding fundamental questions of land plant evolution, they only recently became amenable to experimental investigation, with Anthoceros agrestis being developed as a hornwort model system. Availability of a high-quality genome assembly and a recently developed genetic transformation technique makes A. agrestis an attractive model species for hornworts. Here we describe an updated and optimized transformation protocol for A. agrestis, which can be successfully used to genetically modify one more strain of A. agrestis and three more hornwort species, Anthoceros punctatus, Leiosporoceros dussii, and Phaeoceros carolinianus. The new transformation method is less laborious, faster, and results in the generation of greatly increased numbers of transformants compared with the previous method. We have also developed a new selection marker for transformation. Finally, we report the development of a set of different cellular localization signal peptides for hornworts providing new tools to better understand the hornwort cell biology.

Morphological, genetic and ecological divergence in near-cryptic bryophyte species widespread in the Holarctic: the Dicranum acutifolium complex (Dicranales) revisited in the Alps.

There is mounting evidence that reproductively isolated, but morphologically weakly differentiated species (so-called cryptic species) represent a substantial part of biological diversity, especially in bryophytes. We assessed the evolutionary history and ecological differentiation of a species pair, Dicranum brevifolium and D. septentrionale, which have overlapping ranges in the Holarctic. Despite their morphological similarity, we found similar genetic differentiation as between morphologically well-differentiated Dicranum species. Moreover, we detected gene tree discordance between plastid and nuclear markers, but neither of the two datasets resolved the two as sister species. The signal in trnL-trnF better reflects the morphological and ecological affinities and indicates a close relationship while ITS sequence data resolved the two taxa as phylogenetically distantly related. The discordance is probably unrelated to the ecological differentiation of D. septentrionale to colonise subneutral to alkaline substrates (vs. acidic in D. brevifolium), because this ability is rare in the genus and shared with D. acutifolium. This taxon is the closest relative of D. septentrionale according to the trnL-trnF data and does not share the discordance in ITS. We furthermore demonstrate that beside D. acutifolium, both D. septentrionale and D. brevifolium occur in the Alps but D. brevifolium is most likely rarer. Based on morphological analyses including factor analysis for mixed data of 45 traits we suggest treating the latter two as near-cryptic species and we recommend verifying morphological determinations molecularly.

Challenging the term symbiosis in plant-microbe associations to create an understanding across sciences.

Scientific progress relies on clear and consistent definitions for effective communication and collaboration. The term "symbiosis" in the context of plant-microbe associations suffers from diverse interpretations, leading to ambiguity in classification of these associations. This review elaborates on the issue, proposing an inclusive definition as well as a keyword.

Comparative Genomics Elucidates the Origin of a Supergene Controlling Floral Heteromorphism.

Supergenes are nonrecombining genomic regions ensuring the coinheritance of multiple, coadapted genes. Despite the importance of supergenes in adaptation, little is known on how they originate. A classic example of supergene is the S locus controlling heterostyly, a floral heteromorphism occurring in 28 angiosperm families. In Primula, heterostyly is characterized by the cooccurrence of two complementary, self-incompatible floral morphs and is controlled by five genes clustered in the hemizygous, ca. 300-kb S locus. Here, we present the first chromosome-scale genome assembly of any heterostylous species, that of Primula veris (cowslip). By leveraging the high contiguity of the P. veris assembly and comparative genomic analyses, we demonstrated that the S-locus evolved via multiple, asynchronous gene duplications and independent gene translocations. Furthermore, we discovered a new whole-genome duplication in Ericales that is specific to the Primula lineage. We also propose a mechanism for the origin of S-locus hemizygosity via nonhomologous recombination involving the newly discovered two pairs of CFB genes flanking the S locus. Finally, we detected only weak signatures of degeneration in the S locus, as predicted for hemizygous supergenes. The present study provides a useful resource for future research addressing key questions on the evolution of supergenes in general and the S locus in particular: How do supergenes arise? What is the role of genome architecture in the evolution of complex adaptations? Is the molecular architecture of heterostyly supergenes across angiosperms similar to that of Primula?

Different molecular changes underlie the same phenotypic transition: Origins and consequences of independent shifts to homostyly within species.

The repeated transition from outcrossing to selfing is a key topic in evolutionary biology. However, the molecular basis of such shifts has been rarely examined due to lack of knowledge of the genes controlling these transitions. A classic example of mating system transition is the repeated shift from heterostyly to homostyly. Occurring in 28 angiosperm families, heterostyly is characterized by the reciprocal position of male and female sexual organs in two (or three) distinct, usually self-incompatible floral morphs. Conversely, homostyly is characterized by a single, self-compatible floral morph with reduced separation of male and female organs, facilitating selfing. Here, we investigate the origins of homostyly in Primula vulgaris and its microevolutionary consequences by integrating surveys of the frequency of homostyles in natural populations, DNA sequence analyses of the gene controlling the position of female sexual organs (CYPᵀ), and microsatellite genotyping of both progeny arrays and natural populations characterized by varying frequencies of homostyles. As expected, we found that homostyles displace short-styled individuals, but long-style morphs are maintained at low frequencies within populations. We also demonstrated that homostyles repeatedly evolved from short-styled individuals in association with different types of loss-of-function mutations in CYPᵀ. Additionally, homostyly triggers a shift to selfing, promoting increased inbreeding within and genetic differentiation among populations. Our results elucidate the causes and consequences of repeated transitions to homostyly within species, and the putative mechanisms precluding its fixation in P. vulgaris. This study represents a benchmark for future analyses of losses of heterostyly in other angiosperms.

Extensive Genome-Wide Phylogenetic Discordance Is Due to Incomplete Lineage Sorting and Not Ongoing Introgression in a Rapidly Radiated Bryophyte Genus.

The relative importance of introgression for diversification has long been a highly disputed topic in speciation research and remains an open question despite the great attention it has received over the past decade. Gene flow leaves traces in the genome similar to those created by incomplete lineage sorting (ILS), and identification and quantification of gene flow in the presence of ILS is challenging and requires knowledge about the true phylogenetic relationship among the species. We use whole nuclear, plastid, and organellar genomes from 12 species in the rapidly radiated, ecologically diverse, actively hybridizing genus of peatmoss (Sphagnum) to reconstruct the species phylogeny and quantify introgression using a suite of phylogenomic methods. We found extensive phylogenetic discordance among nuclear and organellar phylogenies, as well as across the nuclear genome and the nodes in the species tree, best explained by extensive ILS following the rapid radiation of the genus rather than by postspeciation introgression. Our analyses support the idea of ancient introgression among the ancestral lineages followed by ILS, whereas recent gene flow among the species is highly restricted despite widespread interspecific hybridization known in the group. Our results contribute to phylogenomic understanding of how speciation proceeds in rapidly radiated, actively hybridizing species groups, and demonstrate that employing a combination of diverse phylogenomic methods can facilitate untangling complex phylogenetic patterns created by ILS and introgression.

The first step into phenolic metabolism in the hornwort Anthoceros agrestis: molecular and biochemical characterization of two phenylalanine ammonia-lyase isoforms.

MAIN CONCLUSION: Two isoforms of phenylalanine ammonia-lyase (PAL) have been isolated as cDNA sequences from the hornwort Anthoceros agrestis. The encoded enzymes convert L-phenylalanine and to lower extents L-tyrosine and L-histidine. Thus, the functional presence of the general phenylpropanoid pathway in one of the earliest land plant groups is established. The hornwort Anthoceros agrestis has an elaborated phenolic metabolism resulting in phenolic compounds, such as rosmarinic acid or megacerotonic acid. The general phenylpropanoid pathway is involved in the biosynthesis of these compounds. Two phenylalanine ammonia-lyase (PAL) genes, AaPAL1 and AaPAL2, have been identified in Anthoceros agrestis and the protein with an N-terminal 6xHis-tag heterologously synthesized in Escherichia coli for a full biochemical characterization. Both PAL proteins accept L-phenylalanine, L-tyrosine as well as L-histidine as substrates, although the activity is explicitly the highest with L-phenylalanine. Km values as well as catalytic efficiencies were determined for phenylalanine (Km AaPAL1 39 µM, AaPAL2 18 µM) and tyrosine (Km AaPAL1 3.3 mM, AaPAL2 3.5 mM). In suspension cultures of Anthoceros agrestis, PAL genes were transcribed in parallel to rosmarinic acid (RA) accumulation and both showed highest abundance in the early growth phase. In a phylogenetic tree, both AaPAL amino acid sequences grouped within a clade with PAL amino acid sequences of diverse origin ranging from non-vascular to vascular plants, while most PALs from eudicots and monocots were mainly found in two other clades. The similarity of the hornwort PAL amino acid sequences to PAL sequences from vascular plants is more than 80% showing a strong conservation within the land plants. With this characterization of PALs from Anthoceros agrestis together with former investigations concerning cinnamic acid 4-hydroxylase and 4-coumaric acid CoA-ligase, the functional presence of the general phenylpropanoid pathway in this hornwort is proven.

A pseudomolecule-scale genome assembly of the liverwort Marchantia polymorpha.

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

The hornworts: morphology, evolution and development.

Extant land plants consist of two deeply divergent groups, tracheophytes and bryophytes, which shared a common ancestor some 500 million years ago. While information about vascular plants and the two of the three lineages of bryophytes, the mosses and liverworts, is steadily accumulating, the biology of hornworts remains poorly explored. Yet, as the sister group to liverworts and mosses, hornworts are critical in understanding the evolution of key land plant traits. Until recently, there was no hornwort model species amenable to systematic experimental investigation, which hampered detailed insight into the molecular biology and genetics of this unique group of land plants. The emerging hornwort model species, Anthoceros agrestis, is instrumental in our efforts to better understand not only hornwort biology but also fundamental questions of land plant evolution. To this end, here we provide an overview of hornwort biology and current research on the model plant A. agrestis to highlight its potential in answering key questions of land plant biology and evolution.

An Agrobacterium-mediated stable transformation technique for the hornwort model Anthoceros agrestis.

Despite their key phylogenetic position and their unique biology, hornworts have been widely overlooked. Until recently there was no hornwort model species amenable to systematic experimental investigation. Anthoceros agrestis has been proposed as the model species to study hornwort biology. We have developed an Agrobacterium-mediated method for the stable transformation of A. agrestis, a hornwort model species for which a genetic manipulation technique was not yet available. High transformation efficiency was achieved by using thallus tissue grown under low light conditions. We generated a total of 274 transgenic A. agrestis lines expressing the ß-glucuronidase (GUS), cyan, green, and yellow fluorescent proteins under control of the CaMV 35S promoter and several endogenous promoters. Nuclear and plasma membrane localization with multiple color fluorescent proteins was also confirmed. The transformation technique described here should pave the way for detailed molecular and genetic studies of hornwort biology, providing much needed insight into the molecular mechanisms underlying symbiosis, carbon-concentrating mechanism, RNA editing and land plant evolution in general.

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.

Organellomic data sets confirm a cryptic consensus on (unrooted) land-plant relationships and provide new insights into bryophyte molecular evolution.

PREMISE: Phylogenetic trees of bryophytes provide important evolutionary context for land plants. However, published inferences of overall embryophyte relationships vary considerably. We performed phylogenomic analyses of bryophytes and relatives using both mitochondrial and plastid gene sets, and investigated bryophyte plastome evolution. METHODS: We employed diverse likelihood-based analyses to infer large-scale bryophyte phylogeny for mitochondrial and plastid data sets. We tested for changes in purifying selection in plastid genes of a mycoheterotrophic liverwort (Aneura mirabilis) and a putatively mycoheterotrophic moss (Buxbaumia), and compared 15 bryophyte plastomes for major structural rearrangements. RESULTS: Overall land-plant relationships conflict across analyses, generally weakly. However, an underlying (unrooted) four-taxon tree is consistent across most analyses and published studies. Despite gene coverage patchiness, relationships within mosses, liverworts, and hornworts are largely congruent with previous studies, with plastid results generally better supported. Exclusion of RNA edit sites restores cases of unexpected non-monophyly to monophyly for Takakia and two hornwort genera. Relaxed purifying selection affects multiple plastid genes in mycoheterotrophic Aneura but not Buxbaumia. Plastid genome structure is nearly invariant across bryophytes, but the tufA locus, presumed lost in embryophytes, is unexpectedly retained in several mosses. CONCLUSIONS: A common unrooted tree underlies embryophyte phylogeny, [(liverworts, mosses), (hornworts, vascular plants)]; rooting inconsistency across studies likely reflects substantial distance to algal outgroups. Analyses combining genomic and transcriptomic data may be misled locally for heavily RNA-edited taxa. The Buxbaumia plastome lacks hallmarks of relaxed selection found in mycoheterotrophic Aneura. Autotrophic bryophyte plastomes, including Buxbaumia, hardly vary in overall structure.

The Physcomitrella patens gene atlas project: large-scale RNA-seq based expression data.

High-throughput RNA sequencing (RNA-seq) has recently become the method of choice to define and analyze transcriptomes. For the model moss Physcomitrella patens, although this method has been used to help analyze specific perturbations, no overall reference dataset has yet been established. In the framework of the Gene Atlas project, the Joint Genome Institute selected P. patens as a flagship genome, opening the way to generate the first comprehensive transcriptome dataset for this moss. The first round of sequencing described here is composed of 99 independent libraries spanning 34 different developmental stages and conditions. Upon dataset quality control and processing through read mapping, 28 509 of the 34 361 v3.3 gene models (83%) were detected to be expressed across the samples. Differentially expressed genes (DEGs) were calculated across the dataset to permit perturbation comparisons between conditions. The analysis of the three most distinct and abundant P. patens growth stages - protonema, gametophore and sporophyte - allowed us to define both general transcriptional patterns and stage-specific transcripts. As an example of variation of physico-chemical growth conditions, we detail here the impact of ammonium supplementation under standard growth conditions on the protonemal transcriptome. Finally, the cooperative nature of this project allowed us to analyze inter-laboratory variation, as 13 different laboratories around the world provided samples. We compare differences in the replication of experiments in a single laboratory and between different laboratories.

The lichen symbiosis re-viewed through the genomes of Cladonia grayi and its algal partner Asterochloris glomerata.

BACKGROUND: Lichens, encompassing 20,000 known species, are symbioses between specialized fungi (mycobionts), mostly ascomycetes, and unicellular green algae or cyanobacteria (photobionts). Here we describe the first parallel genomic analysis of the mycobiont Cladonia grayi and of its green algal photobiont Asterochloris glomerata. We focus on genes/predicted proteins of potential symbiotic significance, sought by surveying proteins differentially activated during early stages of mycobiont and photobiont interaction in coculture, expanded or contracted protein families, and proteins with differential rates of evolution. RESULTS: A) In coculture, the fungus upregulated small secreted proteins, membrane transport proteins, signal transduction components, extracellular hydrolases and, notably, a ribitol transporter and an ammonium transporter, and the alga activated DNA metabolism, signal transduction, and expression of flagellar components. B) Expanded fungal protein families include heterokaryon incompatibility proteins, polyketide synthases, and a unique set of G-protein α subunit paralogs. Expanded algal protein families include carbohydrate active enzymes and a specific subclass of cytoplasmic carbonic anhydrases. The alga also appears to have acquired by horizontal gene transfer from prokaryotes novel archaeal ATPases and Desiccation-Related Proteins. Expanded in both symbionts are signal transduction components, ankyrin domain proteins and transcription factors involved in chromatin remodeling and stress responses. The fungal transportome is contracted, as are algal nitrate assimilation genes. C) In the mycobiont, slow-evolving proteins were enriched for components involved in protein translation, translocation and sorting. CONCLUSIONS: The surveyed genes affect stress resistance, signaling, genome reprogramming, nutritional and structural interactions. The alga carries many genes likely transferred horizontally through viruses, yet we found no evidence of inter-symbiont gene transfer. The presence in the photobiont of meiosis-specific genes supports the notion that sexual reproduction occurs in Asterochloris while they are free-living, a phenomenon with implications for the adaptability of lichens and the persistent autonomy of the symbionts. The diversity of the genes affecting the symbiosis suggests that lichens evolved by accretion of many scattered regulatory and structural changes rather than through introduction of a few key innovations. This predicts that paths to lichenization were variable in different phyla, which is consistent with the emerging consensus that ascolichens could have had a few independent origins.

A fern AINTEGUMENTA gene mirrors BABY BOOM in promoting apogamy in Ceratopteris richardii.

Asexual reproduction is widespread in land plants, including ferns where 10% of all species are obligate asexuals. In these ferns, apogamous sporophytes are generated directly from gametophytes, bypassing fertilization. In the model fern Ceratopteris richardii, a sexual species, apogamy can be induced by culture on high sugar media. BABY BOOM (BBM) genes in angiosperms are known to promote somatic embryogenesis, which like apogamy produce sporophytes without fertilization. Here, a Brassica napus BBM (BnBBM) was used to investigate genetic similarity between apogamy in ferns and somatic embryogenesis in angiosperms. A C. richardii transcriptome was constructed from which one AINTEGUMENTA-LIKE unigene, CrANT, was identified. Whole mount in situ hybridization showed that CrANT is expressed in sperm and fertilized eggs. Phylogenetic analysis grouped CrANT with other non-seed-plant ANT genes to the euANT clade but in a branch separate from BBM genes. Overexpression of CrANT or BnBBM promotes apogamy in C. richardii without sugar supplement. CrANT knockdown gametophytes responded weakly to sugar for apogamy promotion. Theses results suggest some genetic conservation between apogamy and somatic embryogenesis and that such asexual reproduction may be ancient.

The Sphagnome Project: enabling ecological and evolutionary insights through a genus-level sequencing project.

Considerable progress has been made in ecological and evolutionary genetics with studies demonstrating how genes underlying plant and microbial traits can influence adaptation and even 'extend' to influence community structure and ecosystem level processes. Progress in this area is limited to model systems with deep genetic and genomic resources that often have negligible ecological impact or interest. Thus, important linkages between genetic adaptations and their consequences at organismal and ecological scales are often lacking. Here we introduce the Sphagnome Project, which incorporates genomics into a long-running history of Sphagnum research that has documented unparalleled contributions to peatland ecology, carbon sequestration, biogeochemistry, microbiome research, niche construction, and ecosystem engineering. The Sphagnome Project encompasses a genus-level sequencing effort that represents a new type of model system driven not only by genetic tractability, but by ecologically relevant questions and hypotheses.

How Do Cold-Adapted Plants Respond to Climatic Cycles? Interglacial Expansion Explains Current Distribution and Genomic Diversity in Primula farinosa L.

Understanding the effects of past climatic fluctuations on the distribution and population-size dynamics of cold-adapted species is essential for predicting their responses to ongoing global climate change. In spite of the heterogeneity of cold-adapted species, two main contrasting hypotheses have been proposed to explain their responses to Late Quaternary glacial cycles, namely, the interglacial contraction versus the interglacial expansion hypotheses. Here, we use the cold-adapted plant Primula farinosa to test two demographic models under each of the two alternative hypotheses and a fifth, null model. We first approximate the time and extent of demographic contractions and expansions during the Late Quaternary by projecting species distribution models across the last 72 ka. We also generate genome-wide sequence data using a Reduced Representation Library approach to reconstruct the spatial structure, genetic diversity, and phylogenetic relationships of lineages within P. farinosa. Finally, by integrating the results of climatic and genomic analyses in an Approximate Bayesian Computation framework, we propose the most likely model for the extent and direction of population-size changes in $P$. farinosa through the Late Quaternary. Our results support the interglacial expansion of $P$. farinosa, differing from the prevailing paradigm that the observed distribution of cold-adapted species currently fragmented in high altitude and latitude regions reflects the consequences of postglacial contraction processes.

Divergent evolution and niche differentiation within the common peatmoss Sphagnum magellanicum.

PREMISE OF THE STUDY: Populations with phenotypic polymorphism in discrete characters may be good models for investigating genome evolution and speciation. Sphagnum magellanicum Brid. is found throughout the northern hemisphere, and despite considerable variation in morphological characters, it is considered one of the least taxonomically controversial peatmoss species. We have observed two main morphs of the species associated with different microhabitats. Here we investigated the genomic and environmental basis of this intraspecific morphological variation. METHODS: We conducted transplant and common garden experiments to test whether the two morphs are genetically differentiated. We then used RAD-sequencing to quantify the genomic divergence between the morphs and approximate Bayesian computation (ABC) to infer the most likely demographic scenario explaining the genome-wide differentiation of the two morphs. KEY RESULTS: We found that genomic differentiation between the two morphs is unexpectedly high and that several of the differentiated morphological characters have a genetic basis. Using simulation approaches, we found support for a scenario of ancient divergence followed by recent secondary contact. CONCLUSIONS: We show that the two morphs represent the two main genetic clusters previously found worldwide. Our results demonstrate that relatively minor morphological differentiation in a presumed phenotypically plastic peatmoss may be associated with massive divergence across the genome.

Analyses of transcriptome sequences reveal multiple ancient large-scale duplication events in the ancestor of Sphagnopsida (Bryophyta).

The goal of this research was to investigate whether there has been a whole-genome duplication (WGD) in the ancestry of Sphagnum (peatmoss) or the class Sphagnopsida, and to determine if the timing of any such duplication(s) and patterns of paralog retention could help explain the rapid radiation and current ecological dominance of peatmosses. RNA sequencing (RNA-seq) data were generated for nine taxa in Sphagnopsida (Bryophyta). Analyses of frequency plots for synonymous substitutions per synonymous site (Ks ) between paralogous gene pairs and reconciliation of 578 gene trees were conducted to assess evidence of large-scale or genome-wide duplication events in each transcriptome. Both Ks frequency plots and gene tree-based analyses indicate multiple duplication events in the history of the Sphagnopsida. The most recent WGD event predates divergence of Sphagnum from the two other genera of Sphagnopsida. Duplicate retention is highly variable across species, which might be best explained by local adaptation. Our analyses indicate that the last WGD could have been an important factor underlying the diversification of peatmosses and facilitated their rise to ecological dominance in peatlands. The timing of the duplication events and their significance in the evolutionary history of peat mosses are discussed.