Extraordinary preservation of gene collinearity over three hundred million years revealed in homosporous lycophytes.

Homosporous lycophytes (Lycopodiaceae) are a deeply diverged lineage in the plant tree of life, having split from heterosporous lycophytes (Selaginella and Isoetes) ~400 Mya. Compared to the heterosporous lineage, Lycopodiaceae has markedly larger genome sizes and remains the last major plant clade for which no chromosome-level assembly has been available. Here, we present chromosomal genome assemblies for two homosporous lycophyte species, the allotetraploid Huperzia asiatica and the diploid Diphasiastrum complanatum. Remarkably, despite that the two species diverged ~350 Mya, around 30% of the genes are still in syntenic blocks. Furthermore, both genomes had undergone independent whole genome duplications, and the resulting intragenomic syntenies have likewise been preserved relatively well. Such slow genome evolution over deep time is in stark contrast to heterosporous lycophytes and is correlated with a decelerated rate of nucleotide substitution. Together, the genomes of H. asiatica and D. complanatum not only fill a crucial gap in the plant genomic landscape but also highlight a potentially meaningful genomic contrast between homosporous and heterosporous species.

Biolistics-mediated transformation of hornworts and its application to study pyrenoid protein localization.

Hornworts are a deeply diverged lineage of bryophytes that are sister to mosses and liverworts. Hornworts have an array of unique features that can be leveraged to illuminate not only the early evolution of land plants, but also alternative paths for nitrogen and carbon assimilation via cyanobacterial symbiosis and a pyrenoid-based CO2-concentrating mechanism (CCM), respectively. Despite this, hornworts are one of the few plant lineages with limited available genetic tools. Here we report an efficient biolistics method for generating transient-expression and stable transgenic lines in the model hornwort, Anthoceros agrestis. An average of 569 (± 268) cells showed transient expression per bombardment, with green fluorescent protein expression observed within 48-72 hours. A total of 81 stably transformed lines were recovered across three separate experiments, averaging six lines per bombardment. We followed the same method to transiently transform nine additional hornwort species, and obtained stable transformants from one. This method was further used to verify the localization of Rubisco and Rubisco activase in pyrenoids, which are central proteins for CCM function. Together, our biolistics approach offers key advantages over existing methods as it enables rapid transient expression and can be applied to widely diverse hornwort species.

Population genomics of the Isoetes appalachiana (Isoetaceae) complex supports a 'diploids-first' approach to conservation.

BACKGROUND AND AIMS: Allopolyploidy is an important driver of diversification and a key contributor to genetic novelty across the tree of life. However, many studies have questioned the importance of extant polyploid lineages, suggesting that the vast majority may constitute evolutionary 'dead ends'. This has important implications for conservation efforts where polyploids and diploid progenitors often compete for wildlife management resources. Isoetes appalachiana is an allotetraploid that is broadly distributed throughout the eastern USA alongside its diploid progenitors, I. valida and I. engelmannii. As such, this species complex provides an excellent opportunity to investigate the processes that underpin the formation and survival of allopolyploid lineages. METHODS: Here we utilized RADseq and whole-chloroplast sequencing to unravel the demographic and evolutionary history of hybridization in this widespread species complex. We developed a modified protocol for phasing RADseq loci from an allopolyploid in order to examine each progenitor's genetic contribution independently in a phylogenetic context. Additionally, we conducted population-level analyses to examine genetic diversity and evidence of gene flow within species. KEY RESULTS: Isoetes appalachiana is the product of multiple phylogenetic origins, suggesting that formation and establishment of allopolyploids are common in this group. Hybridization appears to be unidirectional, with I. engelmannii consistently being the maternal progenitor. Additionally, we find that polyploid lineages are genetically isolated, rarely if ever experiencing gene flow between geographically distinct populations. CONCLUSIONS: Allopolyploid lineages of I. appalachiana appear to form frequently and experience a high degree of genetic isolation following formation. Thus, our results appear to corroborate the hypothesis that the vast majority of recently formed polyploids may represent evolutionary dead ends. However, this does not necessarily lessen the evolutionary importance or ecological impact of polyploidy per se. Accordingly, we propose a conservation strategy that prioritizes diploid taxa, thus preserving downstream processes that recurrently generate allopolyploid diversity.

The ancestral chromatin landscape of land plants.

Recent studies have shown that correlations between chromatin modifications and transcription vary among eukaryotes. This is the case for marked differences between the chromatin of the moss Physcomitrium patens and the liverwort Marchantia polymorpha. Mosses and liverworts diverged from hornworts, altogether forming the lineage of bryophytes that shared a common ancestor with land plants. We aimed to describe chromatin in hornworts to establish synapomorphies across bryophytes and approach a definition of the ancestral chromatin organization of land plants. We used genomic methods to define the 3D organization of chromatin and map the chromatin landscape of the model hornwort Anthoceros agrestis. We report that nearly half of the hornwort transposons were associated with facultative heterochromatin and euchromatin and formed the center of topologically associated domains delimited by protein coding genes. Transposons were scattered across autosomes, which contrasted with the dense compartments of constitutive heterochromatin surrounding the centromeres in flowering plants. Most of the features observed in hornworts are also present in liverworts or in mosses but are distinct from flowering plants. Hence, the ancestral genome of bryophytes was likely a patchwork of units of euchromatin interspersed within facultative and constitutive heterochromatin. We propose this genome organization was ancestral to land plants.

An ancient route towards salicylic acid and its implications for the perpetual Trichormus-Azolla symbiosis.

Despite its small size, the water fern Azolla is a giant among plant symbioses. Within each of its leaflets, a specialized leaf cavity is home to a population of nitrogen-fixing cyanobacteria (cyanobionts). Although a number of plant-cyanobiont symbioses exist, Azolla is unique in that its symbiosis is perpetual: the cyanobionts are inherited during sexual and vegetative propagation. What underpins the communication between the two partners? In angiosperms, the phytohormone salicylic acid (SA) is a well-known regulator of plant-microbe interactions. Using high-performance liquid chromatography-tandem mass spectrometry, we pinpoint the presence of SA in the fern. Comparative genomics and phylogenetics on SA biosynthesis genes across Chloroplastida reveal that the entire Phenylalanine ammonia-lyase-dependent pathway likely existed in the last common ancestor of land plants. Indeed, Azolla filiculoides secondarily lost its isochorismate synthase but has the genetic competence to derive SA from benzoic acid; the presence of SA in artificially cyanobiont-free Azolla supports the existence of this route. Global gene expression data and SA levels from cyanobiont-containing and -free A. filiculoides link SA synthesis with the symbioses: SA appears to induce cyanobacterial proliferation, whereas removal of the symbiont results in reduced SA levels in a nitrogen-dependent manner.

Nostoc Talks Back: Temporal Patterns of Differential Gene Expression During Establishment of Anthoceros-Nostoc Symbiosis.

Endosymbiotic associations between hornworts and nitrogen-fixing cyanobacteria form when the plant is limited for combined nitrogen (N). We generated RNA-seq data to examine temporal gene expression patterns during the culturing of N-starved Anthoceros punctatus in the absence and the presence of symbiotic cyanobacterium Nostoc punctiforme. In symbiont-free A. punctatus gametophytes, N starvation caused downregulation of chlorophyll content and chlorophyll fluorescence characteristics as well as transcription of photosynthesis-related genes. This downregulation was reversed in A. punctatus cocultured with N. punctiforme, corresponding to the provision by the symbiont of N2-derived NH4+, which commenced within 5 days of coculture and reached a maximum by 14 days. We also observed transient increases in transcription of ammonium and nitrate transporters in a N. punctiforme-dependent manner as well as that of a SWEET transporter that was initially independent of N2-derived NH4+. The temporal patterns of differential gene expression indicated that N. punctiforme transmits signals that impact gene expression to A. punctatus both prior to and after its provision of fixed N. This study is the first illustrating the temporal patterns of gene expression during establishment of an endosymbiotic nitrogen-fixing association in this monophyletic evolutionary lineage of land plants. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Dynamic plastid and mitochondrial genomes in Chaetopeltidales (Chlorophyceae) and characterization of a new chlorophyte taxon.

PREMISE: Chaetopeltidales is a poorly characterized order in the Chlorophyceae, with only two plastid and no mitochondrial genomes published. Here we describe a new taxon in Chaetopeltidales, Gormaniella terricola gen. et sp. nov. and characterize both of its organellar genomes. METHODS: Gormaniella terricola was inadvertently isolated from a surface-sterilized hornwort thallus. Light microscopy was used to characterize its vegetative morphology. Organellar genomes were assembled, annotated, and analyzed using a variety of software packages. RESULTS: The mitochondrial genome (66,927 bp) represents the first complete mitochondrial genome published for Chaetopeltidales. The chloroplast genome, measuring 428,981 bp, is one of the largest plastid genomes published to date and shares this large size and an incredible number of short, dispersed repeats with the other sequenced chloroplast genomes in Chaetopeltidales. Despite these shared features, the chloroplast genomes of Chaetopeltidales appear to be highly rearranged when compared to one another, with numerous inversions, translocations, and duplications, suggesting a particularly dynamic chloroplast genome. Both the chloroplast and mitochondrial genomes of G. terricola contain a number of mobile group I and group II introns, which appear to have invaded separately. Three of the introns within the mitochondrial genome encode homing endonucleases that are phylogenetically nested within those found in fungi, rather than algae, suggesting a possible case of horizontal gene transfer. CONCLUSIONS: These results help to shed light on a poorly understood group of algae and their unusual organellar genomes, raising additional questions about the unique patterns of genome evolution within Chaetopeltidales.

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.

The diversity and community structure of symbiotic cyanobacteria in hornworts inferred from long-read amplicon sequencing.

PREMISE: Nitrogen-fixing endosymbioses with cyanobacteria have evolved independently in five very different plant lineages. Expanding knowledge of these symbioses promises to improve the understanding of symbiosis evolution and broaden the toolkit for agricultural engineering to reduce artificial fertilizer use. Here we focused on hornworts, a bryophyte lineage in which all members host cyanobacteria, and investigated factors shaping the diversity of their cyanobiont communities. METHODS: We sampled hornworts and adjacent soils in upstate New York throughout the hornwort growing season. We included all three sympatric hornwort species in the area, allowing us to directly compare partner selectivity. To profile cyanobacteria communities, we established a metabarcoding protocol targeting rbcL-X with PacBio long reads. RESULTS: The hornwort cyanobionts detected were phylogenetically diverse, including clades that do not contain other known plant symbionts. We found significant overlap between hornwort cyanobionts and soil cyanobacteria, a pattern not previously reported in other plant-cyanobacteria symbioses. Cyanobiont communities differed between host plants only centimeters apart, but we did not detect an effect of sampling time or host species on the cyanobacterial community structure. CONCLUSIONS: This study expands the phylogenetic diversity of known symbiotic cyanobacteria. Our analyses suggest that hornwort cyanobionts have a tight connection to the soil background, and we found no evidence that time within growing season, host species, or distance at the scale of meters strongly govern cyanobacteria community assembly. This study provides a critical foundation for further study of the ecology, evolution, and interaction dynamics of plant-cyanobacteria symbiosis.

On the evolutionary significance of horizontal gene transfers in plants.

Horizontal gene transfer (HGT) has long been seen as a crucial process in the evolution of prokaryotic species, but until recently it was thought to have little, if any, effect on the evolution of eukaryotic life forms. Detecting and describing HGT events in eukaryotes is difficult, making this phenomenon at times controversial. However, modern advances in genomics and bioinformatics have radically altered our view of HGT in eukaryotes, especially in plants. It now appears that HGT to and from plant lineages is more common than previously suspected. Importantly, the transfer of functional nuclear genes with adaptive significance has been reported in numerous taxa. Here we review several recent studies that have found evidence of the horizontal transfer of nuclear genes, and argue that HGT has undoubtedly had profound impacts on plant evolution as a whole.

A global phylogeny of Stegnogramma ferns (Thelypteridaceae): generic and sectional revision, historical biogeography and evolution of leaf architecture.

The thelypteroid fern genus Stegnogramma s.l. contains around 18-35 species and has a global, cross-continental distribution ranging from tropical to temperate regions. Several genera and infrageneric sections have been recognized previously in Stegnogramma s.l., but their phylogenetic relationships are still unclear. In this study, we present a global phylogeny of Stegnogramma s.l. with the most comprehensive sampling to date and aim to pinpoint the phylogenetic positions of biogeographically and taxonomically important taxa. Based on the reconstructed historical biogeography and character evolution, we propose a new (infra)generic classification and discuss the diversification of Stegnogramma s.l. in a biogeographical context. New names or combinations are made for 12 (infra)species, including transferring the monotypic species of Craspedosorus to Leptogramma. Finally, we discuss a possible link between leaf architecture and ecological adaptation, and hypothesize that the increase in leaf dissection and free-vein proportion is an adaptive feature to cool climates in Stegnogramma s.l.

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

Genetic Analysis of Physcomitrella patens Identifies ABSCISIC ACID NON-RESPONSIVE, a Regulator of ABA Responses Unique to Basal Land Plants and Required for Desiccation Tolerance.

The anatomically simple plants that first colonized land must have acquired molecular and biochemical adaptations to drought stress. Abscisic acid (ABA) coordinates responses leading to desiccation tolerance in all land plants. We identified ABA nonresponsive mutants in the model bryophyte Physcomitrella patens and genotyped a segregating population to map and identify the ABA NON-RESPONSIVE (ANR) gene encoding a modular protein kinase comprising an N-terminal PAS domain, a central EDR domain, and a C-terminal MAPKKK-like domain. anr mutants fail to accumulate dehydration tolerance-associated gene products in response to drought, ABA, or osmotic stress and do not acquire ABA-dependent desiccation tolerance. The crystal structure of the PAS domain, determined to 1.7-Å resolution, shows a conserved PAS-fold that dimerizes through a weak dimerization interface. Targeted mutagenesis of a conserved tryptophan residue within the PAS domain generates plants with ABA nonresponsive growth and strongly attenuated ABA-responsive gene expression, whereas deleting this domain retains a fully ABA-responsive phenotype. ANR orthologs are found in early-diverging land plant lineages and aquatic algae but are absent from more recently diverged vascular plants. We propose that ANR genes represent an ancestral adaptation that enabled drought stress survival of the first terrestrial colonizers but were lost during land plant evolution.

Super-resolution ribosome profiling reveals unannotated translation events in Arabidopsis.

Deep sequencing of ribosome footprints (ribosome profiling) maps and quantifies mRNA translation. Because ribosomes decode mRNA every 3 nt, the periodic property of ribosome footprints could be used to identify novel translated ORFs. However, due to the limited resolution of existing methods, the 3-nt periodicity is observed mostly in a global analysis, but not in individual transcripts. Here, we report a protocol applied to Arabidopsis that maps over 90% of the footprints to the main reading frame and thus offers super-resolution profiles for individual transcripts to precisely define translated regions. The resulting data not only support many annotated and predicted noncanonical translation events but also uncover small ORFs in annotated noncoding RNAs and pseudogenes. A substantial number of these unannotated ORFs are evolutionarily conserved, and some produce stable proteins. Thus, our study provides a valuable resource for plant genomics and an efficient optimization strategy for ribosome profiling in other organisms.

Microbial-type terpene synthase genes occur widely in nonseed land plants, but not in seed plants.

The vast abundance of terpene natural products in nature is due to enzymes known as terpene synthases (TPSs) that convert acyclic prenyl diphosphate precursors into a multitude of cyclic and acyclic carbon skeletons. Yet the evolution of TPSs is not well understood at higher levels of classification. Microbial TPSs from bacteria and fungi are only distantly related to typical plant TPSs, whereas genes similar to microbial TPS genes have been recently identified in the lycophyte Selaginella moellendorffii The goal of this study was to investigate the distribution, evolution, and biochemical functions of microbial terpene synthase-like (MTPSL) genes in other plants. By analyzing the transcriptomes of 1,103 plant species ranging from green algae to flowering plants, putative MTPSL genes were identified predominantly from nonseed plants, including liverworts, mosses, hornworts, lycophytes, and monilophytes. Directed searching for MTPSL genes in the sequenced genomes of a wide range of seed plants confirmed their general absence in this group. Among themselves, MTPSL proteins from nonseed plants form four major groups, with two of these more closely related to bacterial TPSs and the other two to fungal TPSs. Two of the four groups contain a canonical aspartate-rich "DDxxD" motif. The third group has a "DDxxxD" motif, and the fourth group has only the first two "DD" conserved in this motif. Upon heterologous expression, representative members from each of the four groups displayed diverse catalytic functions as monoterpene and sesquiterpene synthases, suggesting these are important for terpene formation in nonseed plants.

Admixture, evolution, and variation in reproductive isolation in the Boechera puberula clade.

BACKGROUND: Hybridization is very common in plants, and the incorporation of new alleles into existing lineages (i.e. admixture) can blur species boundaries. However, admixture also has the potential to increase standing genetic variation. With new sequencing methods, we can now study admixture and reproductive isolation at a much finer scale than in the past. The genus Boechera is an extraordinary example of admixture, with over 400 hybrid derivates of varying ploidy levels. Yet, few studies have assessed admixture in this genus on a genomic scale. RESULTS: In this study, we used Genotyping-by-Sequencing (GBS) to clarify the evolution of the Boechera puberula clade, whose six members are scattered across the western United States. We further assessed patterns of admixture and reproductive isolation within the group, including two additional species (B. stricta and B. retrofracta) that are widespread across North America. Based on 14,815 common genetic variants, we found evidence for some cases of hybridization. We find evidence of both recent and more ancient admixture, and that levels of admixture vary across species. CONCLUSIONS: We present evidence for a monophyletic origin of the B. puberula group, and a split of B. puberula into two subspecies. Further, when inferring reproductive isolation on the basis of presence and absence of admixture, we found that the accumulation of reproductive isolation between species does not seem to occur linearly with time since divergence in this system. We discuss our results in the context of sexuality and asexuality in Boechera.