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.

Convergent evolution of fern nectaries facilitated independent recruitment of ant-bodyguards from flowering plants.

Plant-herbivore interactions reciprocally influence species' evolutionary trajectories. These interactions have led to many physical and chemical defenses across the plant kingdom. Some plants have even evolved indirect defense strategies to outsource their protection to ant bodyguards by bribing them with a sugary reward (nectar). Identifying the evolutionary processes underpinning these indirect defenses provide insight into the evolution of plant-animal interactions. Using a cross-kingdom, phylogenetic approach, we examined the convergent evolution of ant-guarding nectaries across ferns and flowering plants. Here, we discover that nectaries originated in ferns and flowering plants concurrently during the Cretaceous, coinciding with the rise of plant associations in ants. While nectaries in flowering plants evolved steadily through time, ferns showed a pronounced lag of nearly 100 My between their origin and subsequent diversification in the Cenozoic. Importantly, we find that as ferns transitioned from the forest floor into the canopy, they secondarily recruited ant bodyguards from existing ant-angiosperm relationships.

Genomes of multicellular algal sisters to land plants illuminate signaling network evolution.

Zygnematophyceae are the algal sisters of land plants. Here we sequenced four genomes of filamentous Zygnematophyceae, including chromosome-scale assemblies for three strains of Zygnema circumcarinatum. We inferred traits in the ancestor of Zygnematophyceae and land plants that might have ushered in the conquest of land by plants: expanded genes for signaling cascades, environmental response, and multicellular growth. Zygnematophyceae and land plants share all the major enzymes for cell wall synthesis and remodifications, and gene gains shaped this toolkit. Co-expression network analyses uncover gene cohorts that unite environmental signaling with multicellular developmental programs. Our data shed light on a molecular chassis that balances environmental response and growth modulation across more than 600 million years of streptophyte evolution.

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.

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.

Twelve new metagenome-assembled genomes from non-axenic culture of Griffithsia monilis (Rhodophyta).

We report 12 metagenome-assembled genomes from a non-axenic culture of the red alga Griffithsia monilis Harvey, some of which are distantly related to publicly available genomes.

A structure of the relict phycobilisome from a thylakoid-free cyanobacterium.

Phycobilisomes (PBS) are antenna megacomplexes that transfer energy to photosystems II and I in thylakoids. PBS likely evolved from a basic, inefficient form into the predominant hemidiscoidal shape with radiating peripheral rods. However, it has been challenging to test this hypothesis because ancestral species are generally inaccessible. Here we use spectroscopy and cryo-electron microscopy to reveal a structure of a "paddle-shaped" PBS from a thylakoid-free cyanobacterium that likely retains ancestral traits. This PBS lacks rods and specialized ApcD and ApcF subunits, indicating relict characteristics. Other features include linkers connecting two chains of five phycocyanin hexamers (CpcN) and two core subdomains (ApcH), resulting in a paddle-shaped configuration. Energy transfer calculations demonstrate that chains are less efficient than rods. These features may nevertheless have increased light absorption by elongating PBS before multilayered thylakoids with hemidiscoidal PBS evolved. Our results provide insights into the evolution and diversification of light-harvesting strategies before the origin of thylakoids.

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.

Genome sequence of a symbiotic cyanobacterium from the flowering plant Gunnera tinctoria.

Metagenomic analysis of the symbiotic cyanobacteria colonies within Gunnera tinctoria stems revealed a new strain of Nostoc. Here, we report its genome sequence.

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.

Balancing read length and sequencing depth: Optimizing Nanopore long-read sequencing for monocots with an emphasis on the Liliales.

Premise: We present approaches used to generate long-read Nanopore sequencing reads for the Liliales and demonstrate how modifications to standard protocols directly impact read length and total output. The goal is to help those interested in generating long-read sequencing data determine which steps may be necessary for optimizing output and results. Methods: Four species of Calochortus (Liliaceae) were sequenced. Modifications made to sodium dodecyl sulfate (SDS) extractions and cleanup protocols included grinding with a mortar and pestle, using cut or wide-bore tips, chloroform cleaning, bead cleaning, eliminating short fragments, and using highly purified DNA. Results: Steps taken to maximize read length can decrease overall output. Notably, the number of pores in a flow cell is correlated with the overall output, yet we did not see an association between the pore number and the read length or the number of reads produced. Discussion: Many factors contribute to the overall success of a Nanopore sequencing run. We showed the direct impact that several modifications to the DNA extraction and cleaning steps have on the total sequencing output, read size, and number of reads generated. We show a tradeoff between read length and the number of reads and, to a lesser extent, the total sequencing output, all of which are important factors for successful de novo genome assembly.

Structure and dynamics of the Arabidopsis O-fucosyltransferase SPINDLY.

SPINDLY (SPY) in Arabidopsis thaliana is a novel nucleocytoplasmic protein O-fucosyltransferase (POFUT), which regulates diverse developmental processes. Sequence analysis indicates that SPY is distinct from ER-localized POFUTs and contains N-terminal tetratricopeptide repeats (TPRs) and a C-terminal catalytic domain resembling the O-linked-N-acetylglucosamine (GlcNAc) transferases (OGTs). However, the structural feature that determines the distinct enzymatic selectivity of SPY remains unknown. Here we report the cryo-electron microscopy (cryo-EM) structure of SPY and its complex with GDP-fucose, revealing distinct active-site features enabling GDP-fucose instead of UDP-GlcNAc binding. SPY forms an antiparallel dimer instead of the X-shaped dimer in human OGT, and its catalytic domain interconverts among multiple conformations. Analysis of mass spectrometry, co-IP, fucosylation activity, and cryo-EM data further demonstrates that the N-terminal disordered peptide in SPY contains trans auto-fucosylation sites and inhibits the POFUT activity, whereas TPRs 1-5 dynamically regulate SPY activity by interfering with protein substrate binding.

Chromosome-level genomes of multicellular algal sisters to land plants illuminate signaling network evolution.

The filamentous and unicellular algae of the class Zygnematophyceae are the closest algal relatives of land plants. Inferring the properties of the last common ancestor shared by these algae and land plants allows us to identify decisive traits that enabled the conquest of land by plants. We sequenced four genomes of filamentous Zygnematophyceae (three strains of Zygnema circumcarinatum and one strain of Z. cylindricum) and generated chromosome-scale assemblies for all strains of the emerging model system Z. circumcarinatum. Comparative genomic analyses reveal expanded genes for signaling cascades, environmental response, and intracellular trafficking that we associate with multicellularity. Gene family analyses suggest that Zygnematophyceae share all the major enzymes with land plants for cell wall polysaccharide synthesis, degradation, and modifications; most of the enzymes for cell wall innovations, especially for polysaccharide backbone synthesis, were gained more than 700 million years ago. In Zygnematophyceae, these enzyme families expanded, forming co-expressed modules. Transcriptomic profiling of over 19 growth conditions combined with co-expression network analyses uncover cohorts of genes that unite environmental signaling with multicellular developmental programs. Our data shed light on a molecular chassis that balances environmental response and growth modulation across more than 600 million years of streptophyte evolution.

PURC Provides Improved Sequence Inference for Polyploid Phylogenetics and Other Manifestations of the Multiple-Copy Problem.

Inferring the true biological sequences from amplicon mixtures remains a difficult bioinformatics problem. The traditional approach is to cluster sequencing reads by similarity thresholds and treat the consensus sequence of each cluster as an "operational taxonomic unit" (OTU). Recently, this approach has been improved by model-based methods that correct PCR and sequencing errors in order to infer "amplicon sequence variants" (ASVs). To date, ASV approaches have been used primarily in metagenomics, but they are also useful for determining homeologs in polyploid organisms. To facilitate the usage of ASV methods among polyploidy researchers, we incorporated ASV inference alongside OTU clustering in PURC v2.0, a major update to PURC (Pipeline for Untangling Reticulate Complexes). In addition, PURC v2.0 features faster demultiplexing than the original version and has been updated to be compatible with Python 3. In this chapter we present results indicating that using the ASV approach is more likely to infer the correct biological sequences in comparison to the earlier OTU-based PURC and describe how to prepare sequencing data, run PURC v2.0 under several different modes, and interpret the output.

Dynamic genome evolution in a model fern.

The large size and complexity of most fern genomes have hampered efforts to elucidate fundamental aspects of fern biology and land plant evolution through genome-enabled research. Here we present a chromosomal genome assembly and associated methylome, transcriptome and metabolome analyses for the model fern species Ceratopteris richardii. The assembly reveals a history of remarkably dynamic genome evolution including rapid changes in genome content and structure following the most recent whole-genome duplication approximately 60 million years ago. These changes include massive gene loss, rampant tandem duplications and multiple horizontal gene transfers from bacteria, contributing to the diversification of defence-related gene families. The insertion of transposable elements into introns has led to the large size of the Ceratopteris genome and to exceptionally long genes relative to other plants. Gene family analyses indicate that genes directing seed development were co-opted from those controlling the development of fern sporangia, providing insights into seed plant evolution. Our findings and annotated genome assembly extend the utility of Ceratopteris as a model for investigating and teaching plant biology.

Systematics and Plastome Evolution in Schizaeaceae.

While the family Schizaeaceae (Schizaeales) represents only about 0.4% of the extant fern species diversity, it differs from other ferns greatly in gross morphologies, niche preferences, and life histories. One of the most notable features in this family is its mycoheterotrophic life style in the gametophytic stage, which appears to be associated with extensive losses of plastid genes. However, the limited number of sequenced plastomes, and the lack of a well-resolved phylogenetic framework of Schizaeaceae, makes it difficult to gain any further insight. Here, with a comprehensive sampling of ~77% of the species diversity of this family, we first inferred a plastid phylogeny of Schizaeaceae using three DNA regions. To resolve the deep relationships within this family, we then reconstructed a plastome-based phylogeny focusing on a selection of representatives that covered all the major clades. From this phylogenomic backbone, we traced the evolutionary histories of plastid genes and examined whether gene losses were associated with the evolution of gametophytic mycoheterotrophy. Our results reveal that extant Schizaeaceae is comprised of four major clades-Microschizaea, Actinostachys, Schizaea, and Schizaea pusilla. The loss of all plastid NADH-like dehydrogenase (ndh) genes was confirmed to have occurred in the ancestor of extant Schizaeaceae, which coincides with the evolution of mycoheterotrophy in this family. For chlorophyll biosynthesis genes (chl), the losses were interpreted as convergent in Schizaeaceae, and found not only in Actinostachys, a clade producing achlorophyllous gametophytes, but also in S. pusilla with chlorophyllous gametophytes. In addition, we discovered a previously undescribed but phylogenetically distinct species hidden in the Schizaea dichotoma complex and provided a taxonomic treatment and morphological diagnostics for this new species-Schizaea medusa. Finally, our phylogenetic results suggest that the current PPG I circumscription of Schizaea is non-monophyletic, and we therefore proposed a three-genus classification moving a subset of Schizaea species sensu PPG I to a third genus-Microschizaea.