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.

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.

What can hornworts teach us?

The hornworts are a small group of land plants, consisting of only 11 families and approximately 220 species. Despite their small size as a group, their phylogenetic position and unique biology are of great importance. Hornworts, together with mosses and liverworts, form the monophyletic group of bryophytes that is sister to all other land plants (Tracheophytes). It is only recently that hornworts became amenable to experimental investigation with the establishment of Anthoceros agrestis as a model system. In this perspective, we summarize the recent advances in the development of A. agrestis as an experimental system and compare it with other plant model systems. We also discuss how A. agrestis can help to further research in comparative developmental studies across land plants and to solve key questions of plant biology associated with the colonization of the terrestrial environment. Finally, we explore the significance of A. agrestis in crop improvement and synthetic biology applications in general.

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.

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.

Comparative transcriptomics reveals commonalities and differences in the genetic underpinnings of a floral dimorphism.

Distyly, a floral dimorphism associated with heteromorphic self-incompatibility and controlled by the S-locus supergene, evolved independently multiple times. Comparative analyses of the first transcriptome atlas for the main distyly model, Primula veris, with other distylous species produced the following findings. A set of 53 constitutively expressed genes in P. veris did not include any of the housekeeping genes commonly used to normalize gene expression in qPCR experiments. The S-locus gene CYPT acquired its role in controlling style elongation via a change in expression profile. Comparison of genes differentially expressed between floral morphs revealed that brassinosteroids and auxin are the main hormones controlling style elongation in P. veris and Fagopyrum esculentum, respectively. Furthermore, shared biochemical pathways might underlie the expression of distyly in the distantly related P. veris, F. esculentum and Turnera subulata, suggesting a degree of correspondence between evolutionary convergence at phenotypic and molecular levels. Finally, we provide the first evidence supporting the previously proposed hypothesis that distyly supergenes of distantly related species evolved via the recruitment of genes related to the phytochrome-interacting factor (PIF) signaling network. To conclude, this is the first study that discovered homologous genes involved in the control of distyly in distantly related taxa.

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.

Step-by-step protocol for the isolation and transient transformation of hornwort protoplasts.

Premise: A detailed protocol for the protoplast transformation of hornwort tissue is not yet available, limiting molecular biological investigations of these plants and comparative analyses with other bryophytes, which display a gametophyte-dominant life cycle and are critical to understanding the evolution of key land plant traits. Methods and Results: We describe a detailed protocol to isolate and transiently transform protoplasts of the model hornwort Anthoceros agrestis. The digestion of liquid cultures with Driselase yields a high number of viable protoplasts suitable for polyethylene glycol (PEG)-mediated transformation. We also report early signs of protoplast regeneration, such as chloroplast division and cell wall reconstitution. Conclusions: This protocol represents a straightforward method for isolating and transforming A. agrestis protoplasts that is less laborious than previously described approaches. In combination with the recently developed stable genome transformation technique, this work further expands the prospects of functional studies in this model hornwort.

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?

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.

Lipid exchanges drove the evolution of mutualism during plant terrestrialization.

Symbiosis with arbuscular mycorrhizal fungi (AMF) improves plant nutrition in most land plants, and its contribution to the colonization of land by plants has been hypothesized. Here, we identify a conserved transcriptomic response to AMF among land plants, including the activation of lipid metabolism. Using gain of function, we show the transfer of lipids from the liverwort Marchantia paleacea to AMF and its direct regulation by the transcription factor WRINKLED (WRI). Arbuscules, the nutrient-exchange structures, were not formed in loss-of-function wri mutants in M. paleacea, leading to aborted mutualism. Our results show the orthology of the symbiotic transfer of lipids across land plants and demonstrate that mutualism with arbuscular mycorrhizal fungi was present in the most recent ancestor of land plants 450 million years ago.

Charting the genomic landscape of seed-free plants.

During the past few years several high-quality genomes has been published from Charophyte algae, bryophytes, lycophytes and ferns. These genomes have not only elucidated the origin and evolution of early land plants, but have also provided important insights into the biology of the seed-free lineages. However, critical gaps across the phylogeny remain and many new questions have been raised through comparing seed-free and seed plant genomes. Here, we review the reference genomes available and identify those that are missing in the seed-free lineages. We compare patterns of various levels of genome and epigenomic organization found in seed-free plants to those of seed plants. Some genomic features appear to be fundamentally different. For instance, hornworts, Selaginella and most liverworts are devoid of whole-genome duplication, in stark contrast to other land plants. In addition, the distribution of genes and repeats appear to be less structured in seed-free genomes than in other plants, and the levels of gene body methylation appear to be much lower. Finally, we highlight the currently available (or needed) model systems, which are crucial to further our understanding about how changes in genes translate into evolutionary novelties.

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

Evolutionary History of the Marchantia polymorpha Complex.

The potential role of introgression in evolution has gained increased interest in recent years. Although some fascinating examples have been reported, more information is needed to generalize the importance of hybridization and introgression for adaptive divergence. As limited data exist on haploid dominant species, we analyzed genomes of three subspecies of the liverwort Marchantia polymorpha. We used available genomic data for subsp. ruderalis and carried out whole-genome (PacBio) sequencing for one individual each of subsp. montivagans and subsp. polymorpha as well as Illumina resequencing of additional genomes for all three subspecies. The three subspecies were compared against M. paleacea as outgroup. Our analyses revealed separation of the three taxa, but all three possible topologies were richly represented across the genomes, and the underlying divergence order less obvious. This uncertainty could be the result of the divergence of the three subspecies close in time, or that introgression has been frequent since divergence. In particular, we found that pseudo-chromosome 2 in subsp. montivagans was much more diverged than other parts of the genomes. This could either be explained by specific capture of chromosome 2 from an unknown related species through hybridization or by conservation of chromosome 2 despite intermittent or ongoing introgression affecting more permeable parts of the genomes. A higher degree of chromosomal rearrangements on pseudo-chromosome 2 support the second hypothesis. Species tree analyses recovered an overall topology where subsp. montivagans diverged first and subsp. ruderalis and subsp. polymorpha appeared as sister lineages. Each subspecies was associated with its own chloroplast and mitochondrial haplotype group. Our data suggest introgression but refute a previous hypothesis that subsp. ruderalis is a new stabilized hybrid between the other two subspecies.

Transcriptional Landscapes of Divergent Sporophyte Development in Two Mosses, Physcomitrium (Physcomitrella) patens and Funaria hygrometrica.

Understanding the molecular basis of morphological shifts is a fundamental question of evolutionary biology. New morphologies may arise through the birth/death of genes (gene gain/loss) or by reutilizing existing gene sets. Yet, the relative contribution of these two processes to radical morphological shifts is still poorly understood. Here, we use the model system of two mosses, Funaria hygrometrica and Physcomitrium (Physcomitrella) patens, to investigate the molecular mechanisms underlying contrasting sporophyte architectures. We used comparative analysis of time-series expression data for four stages of sporophyte development in both species to address this question in detail. We found that large-scale differences in sporophytic architecture are mainly governed by orthologous (i.e., shared) genes frequently experiencing temporal gene expression shifts between the two species. While the absolute number of species-specific genes expressed during sporophyte development is somewhat smaller, we observed a significant increase of their proportion in preferentially sporophyte expressed genes, suggesting a fundamental role in the sporophyte phase. However, further functional studies are necessary to determine their contribution to diverging sporophyte morphologies. Our results add to the growing set of studies suggesting that radical changes in morphology may rely on the heterochronic expression of conserved regulators.

Anthoceros genomes illuminate the origin of land plants and the unique biology of hornworts.

Hornworts comprise a bryophyte lineage that diverged from other extant land plants >400 million years ago and bears unique biological features, including a distinct sporophyte architecture, cyanobacterial symbiosis and a pyrenoid-based carbon-concentrating mechanism (CCM). Here, we provide three high-quality genomes of Anthoceros hornworts. Phylogenomic analyses place hornworts as a sister clade to liverworts plus mosses with high support. The Anthoceros genomes lack repeat-dense centromeres as well as whole-genome duplication, and contain a limited transcription factor repertoire. Several genes involved in angiosperm meristem and stomatal function are conserved in Anthoceros and upregulated during sporophyte development, suggesting possible homologies at the genetic level. We identified candidate genes involved in cyanobacterial symbiosis and found that LCIB, a Chlamydomonas CCM gene, is present in hornworts but absent in other plant lineages, implying a possible conserved role in CCM function. We anticipate that these hornwort genomes will serve as essential references for future hornwort research and comparative studies across land plants.

Development and characterization of novel SSR markers in the endangered endemic species Ferula sadleriana.

PREMISE: Ferula sadleriana (Apiaceae) is a polycarpic, perennial herb with a very limited range and small populations. It is listed as "endangered" on the IUCN Red List of Threatened Species. Microsatellite markers can contribute to conservation efforts by allowing the study of the genetic structure of its shrinking populations. METHODS AND RESULTS: We used a microsatellite-enriched library combined with an Illumina sequencing approach to develop simple sequence repeat markers in our target species. Out of 44 tested primer pairs, 22 provided specific products, and 13 showed heterologous amplification in the target species. Cross-species amplification was achieved at 20 and 19 loci in two congeneric species, F. soongarica and F. tatarica, respectively. CONCLUSIONS: The primers described here are the first tools that enable the population genetic characterization of F. sadleriana. Our results suggest a wider applicability in the genus Ferula.