Crossroads of assembling a moss genome: navigating contaminants and horizontal gene transfer in the moss Physcomitrellopsis africana.

The first chromosome-scale reference genome of the rare narrow-endemic African moss Physcomitrellopsis africana is presented here. Assembled from 73x nanopore long reads and 163x BGI-seq short reads, the 414 Mb reference comprises 26 chromosomes and 22,925 protein-coding genes (BUSCO: C:94.8%[D:13.9%]). This genome holds two genes that withstood rigorous filtration of microbial contaminants, have no homolog in other land plants and are thus interpreted as resulting from two unique horizontal gene transfers from microbes. Further, Physcomitrellopsis africana shares 176 of the 273 published HGT candidates identified in Physcomitrium patens, but lacks 98 of these, highlighting that perhaps as many as 91 genes were acquired in P. patens in the last 40 million years following its divergence from its common ancestor with P. africana. These observations suggest rather continuous gene gains via HGT followed by potential losses, during the diversification of the Funariaceae. Our findings showcase both dynamic flux in plant HGTs over evolutionarily "short" timescales, alongside enduring impacts of successful integrations, like those still functionally maintained in extant Physcomitrellopsis africana. Furthermore, this study describes the informatic processes employed to distinguish contaminants from candidate HGT events.

Comprehensive phylogenomic time tree of bryophytes reveals deep relationships and uncovers gene incongruences in the last 500 million years of diversification.

PREMISE: Bryophytes form a major component of terrestrial plant biomass, structuring ecological communities in all biomes. Our understanding of the evolutionary history of hornworts, liverworts, and mosses has been significantly reshaped by inferences from molecular data, which have highlighted extensive homoplasy in various traits and repeated bursts of diversification. However, the timing of key events in the phylogeny, patterns, and processes of diversification across bryophytes remain unclear. METHODS: Using the GoFlag probe set, we sequenced 405 exons representing 228 nuclear genes for 531 species from 52 of the 54 orders of bryophytes. We inferred the species phylogeny from gene tree analyses using concatenated and coalescence approaches, assessed gene conflict, and estimated the timing of divergences based on 29 fossil calibrations. RESULTS: The phylogeny resolves many relationships across the bryophytes, enabling us to resurrect five liverwort orders and recognize three more and propose 10 new orders of mosses. Most orders originated in the Jurassic and diversified in the Cretaceous or later. The phylogenomic data also highlight topological conflict in parts of the tree, suggesting complex processes of diversification that cannot be adequately captured in a single gene-tree topology. CONCLUSIONS: We sampled hundreds of loci across a broad phylogenetic spectrum spanning at least 450 Ma of evolution; these data resolved many of the critical nodes of the diversification of bryophytes. The data also highlight the need to explore the mechanisms underlying the phylogenetic ambiguity at specific nodes. The phylogenomic data provide an expandable framework toward reconstructing a comprehensive phylogeny of this important group of plants.

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

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

Phylotranscriptomics of liverworts: revisiting the backbone phylogeny and ancestral gene duplications.

BACKGROUND AND AIMS: With some 7300 extant species, liverworts (Marchantiophyta) represent one of the major land plant lineages. The backbone relationships, such as the phylogenetic position of Ptilidiales, and the occurrence and timing of whole-genome duplications, are still contentious. METHODS: Based on analyses of the newly generated transcriptome data for 38 liverworts and complemented with those publicly available, we reconstructed the evolutionary history of liverworts and inferred gene duplication events along the 55 taxon liverwort species tree. KEY RESULTS: Our phylogenomic study provided an ordinal-level liverwort nuclear phylogeny and identified extensive gene tree conflicts and cyto-nuclear incongruences. Gene duplication analyses based on integrated phylogenomics and Ks distributions indicated no evidence of whole-genome duplication events along the backbone phylogeny of liverworts. CONCLUSIONS: With a broadened sampling of liverwort transcriptomes, we re-evaluated the backbone phylogeny of liverworts, and provided evidence for ancient hybridizations followed by incomplete lineage sorting that shaped the deep evolutionary history of liverworts. The lack of whole-genome duplication during the deep evolution of liverworts indicates that liverworts might represent one of the few major embryophyte lineages whose evolution was not driven by whole-genome duplications.

The Subantarctic Rayadito (Aphrastura subantarctica), a new bird species on the southernmost islands of the Americas.

We describe a new taxon of terrestrial bird of the genus Aphrastura (rayaditos) inhabiting the Diego Ramírez Archipelago, the southernmost point of the American continent. This archipelago is geographically isolated and lacks terrestrial mammalian predators as well as woody plants, providing a contrasted habitat to the forests inhabited by the other two Aphrastura spp. Individuals of Diego Ramírez differ morphologically from Aphrastura spinicauda, the taxonomic group they were originally attributed to, by their larger beaks, longer tarsi, shorter tails, and larger body mass. These birds move at shorter distances from ground level, and instead of nesting in cavities in trees, they breed in cavities in the ground, reflecting different life-histories. Both taxa are genetically differentiated based on mitochondrial and autosomal markers, with no evidence of current gene flow. Although further research is required to define how far divergence has proceeded along the speciation continuum, we propose A. subantarctica as a new taxonomic unit, given its unique morphological, genetic, and behavioral attributes in a non-forested habitat. The discovery of this endemic passerine highlights the need to monitor and conserve this still-pristine archipelago devoid of exotic species, which is now protected by the recently created Diego Ramírez Islands-Drake Passage Marine Park.

The Cycas genome and the early evolution of seed plants.

Cycads represent one of the most ancient lineages of living seed plants. Identifying genomic features uniquely shared by cycads and other extant seed plants, but not non-seed-producing plants, may shed light on the origin of key innovations, as well as the early diversification of seed plants. Here, we report the 10.5-Gb reference genome of Cycas panzhihuaensis, complemented by the transcriptomes of 339 cycad species. Nuclear and plastid phylogenomic analyses strongly suggest that cycads and Ginkgo form a clade sister to all other living gymnosperms, in contrast to mitochondrial data, which place cycads alone in this position. We found evidence for an ancient whole-genome duplication in the common ancestor of extant gymnosperms. The Cycas genome contains four homologues of the fitD gene family that were likely acquired via horizontal gene transfer from fungi, and these genes confer herbivore resistance in cycads. The male-specific region of the Y chromosome of C. panzhihuaensis contains a MADS-box transcription factor expressed exclusively in male cones that is similar to a system reported in Ginkgo, suggesting that a sex determination mechanism controlled by MADS-box genes may have originated in the common ancestor of cycads and Ginkgo. The C. panzhihuaensis genome provides an important new resource of broad utility for biologists.

Chromosome-Level Genome Assemblies of Two Hypnales (Mosses) Reveal High Intergeneric Synteny.

Mosses compose one of the three lineages of bryophytes. Today, about 13,000 species of mosses are recognized from across the globe, and at least one-third of this diversity composes the Hypnales, a lineage characterized by an early rapid radiation. We sequenced and de novo assembled the genomes of two hypnalean mosses, namely Entodon seductrix and Hypnum curvifolium, based on the 10x genomics and Hi-C data. The genome assemblies of E. seductrix and H. curvifolium comprise 348.4 and 262.0 Mb, respectively, estimated by k-mer analyses to represent 93.3% and 97.2% of their total genome size. Both genomes were assembled at the chromosome level, with scaffold N50 of 30.0 and 20.7 Mb, respectively. The annotated genome of E. seductrix comprises 25,801 protein-coding genes and that of H. curvifolium 29,077, estimated to represent 96.8% and 97.2%, respectively, of the total gene spaces based on BUSCO (Benchmarking Universal Single-Copy Ortholog) assessment. For both genomes, most contigs were anchored to the largest 11 pseudomolecules, corresponding to the 11 chromosomes of the two species, and each with a putative sex-related chromosome characterized by low gene density. The chromosomes of E. seductrix and H. curvifolium are highly syntenic, suggests limited architectural shifts occurred following the rapid radiation of the Hypnales. We compared their genomic features to the model moss Physcomitrium patens. The hypnalean moss genomes lack signatures of recent whole-genome duplication. The presented high-quality moss genomes provide new resources for comparative genomics to potentially unveil the genomic evolution of derived moss lineages.

Construction of DNA Tools for Hyperexpression in Marchantia Chloroplasts.

Chloroplasts are attractive platforms for synthetic biology applications since they are capable of driving very high levels of transgene expression, if mRNA production and stability are properly regulated. However, plastid transformation is a slow process and currently limited to a few plant species. The liverwort Marchantia polymorpha is a simple model plant that allows rapid transformation studies; however, its potential for protein hyperexpression has not been fully exploited. This is partially due to the fact that chloroplast post-transcriptional regulation is poorly characterized in this plant. We have mapped patterns of transcription in Marchantia chloroplasts. Furthermore, we have obtained and compared sequences from 51 bryophyte species and identified putative sites for pentatricopeptide repeat protein binding that are thought to play important roles in mRNA stabilization. Candidate binding sites were tested for their ability to confer high levels of reporter gene expression in Marchantia chloroplasts, and levels of protein production and effects on growth were measured in homoplastic transformed plants. We have produced novel DNA tools for protein hyperexpression in this facile plant system that is a test-bed for chloroplast engineering.

Plastid genomes and phylogenomics of liverworts (Marchantiophyta): Conserved genome structure but highest relative plastid substitution rate in land plants.

With some 7300 species of small nonvascular spore-producing plants, liverworts represent one of the major lineages of land plants. Although multi-locus molecular phylogenetic studies have elucidated relationships of liverworts at different taxonomic categories, the backbone phylogeny of liverworts is still to be fully resolved, especially for the placement of Ptilidiales and the relationships within Jungermanniales and Marchantiales. Here, we provided phylogenomic inferences of liverworts based on 42 newly sequenced and 24 published liverwort plastid genomes representing all but two orders of liverworts, and characterized the evolution of the plastome in liverworts. The structure of the plastid genome is overall conserved across the phylogeny of liverworts, with only two structural variants detected from simple thalloids, besides 18 out of 43 liverwort genera showing intron variations in their plastomes. Complex thalloid liverworts maintain the most plastid genes, and seem to undergo fewer gene deletions and pseudogenization events than other liverworts. Plastid phylogenetic inferences yielded mostly robustly supported relationships, and consistently resolved Ptilidiales as the sister to Porellales. The relative ratio of silent substitutions across the three genetic compartments (i.e., 1:15:10, for mitochondrial:plastid:nuclear) suggests that liverwort plastid genes have the potential to evolve faster than their nuclear counterparts, unlike in any other major land plant lineages where the mutation rate of nuclear genes overwhelm those of their plastid and mitochondrial counterparts.

High gene space divergence contrasts with frozen vegetative architecture in the moss family Funariaceae.

A new paradigm has slowly emerged regarding the diversification of bryophytes, with inferences from molecular data highlighting a dynamic evolution of their genome. However, comparative studies of expressed genes among closely related taxa is so far missing. Here we contrast the dimensions of the vegetative transcriptome of Funaria hygrometrica and Physcomitrium pyriforme against the genome of their relative, Physcomitrium (Physcomitrella) patens. These three species of Funariaceae share highly conserved vegetative bodies, and are partially sympatric, growing on mineral soil in mostly temperate regions. We analyzed the vegetative gametophytic transcriptome of F. hygrometrica and P. pyriforme and mapped short reads, transcripts, and proteins to the genome and gene space of P. patens. Only about half of the transcripts of F. hygrometrica map to their ortholog in P. patens, whereas at least 90% of those of P. pyriforme align to loci in P. patens. Such divergence is unexpected given the high morphological similarity of the gametophyte but reflects the estimated times of divergence of F. hygrometrica and P. pyriforme from P. patens, namely 55 and 20 mya, respectively. The newly sampled transcriptomes bear signatures of at least one, rather ancient, whole genome duplication (WGD), which may be shared with one reported for P. patens. The transcriptomes of F. hygrometrica and P. pyriforme reveal significant contractions or expansions of different gene families. While transcriptomes offer only an incomplete estimate of the gene space, the high number of transcripts obtained suggest a significant divergence in gene sequences, and gene number among the three species, indicative of a rather strong, dynamic genome evolution, shaped in part by whole, partial or localized genome duplication. The gene ontology of their specific and rapidly-evolving protein families, suggests that the evolution of the Funariaceae may have been driven by the diversification of metabolic genes that may optimize the adaptations to environmental conditions, a hypothesis well in line with ecological patterns in the genetic diversity and structure in seed plants.

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.

Population Genomics and Phylogeography of a Clonal Bryophyte With Spatially Separated Sexes and Extreme Sex Ratios.

The southern Appalachian (SA) is one of the most biodiversity-rich areas in North America and has been considered a refugium for many disjunct plant species, from the last glacial period to the present. Our study focuses on the SA clonal hornwort, Nothoceros aenigmaticus J. C. Villarreal & K. D. McFarland. This hornwort was described from North Carolina and is widespread in the SA, growing on rocks near or submerged in streams in six and one watersheds of the Tennessee (TR) and Alabama (AR) Rivers, respectively. Males and female populations occur in different watersheds, except in the Little Tennessee (TN) River where an isolated male population exists ca. 48 km upstream from the female populations. The sex ratio of 1:0 seems extreme in each population. In this study, we use nuclear and organellar microsatellites from 250 individuals from six watersheds (seven populations) in the SA region and two populations from Mexico (23 individuals). We, then, selected 86 individuals from seven populations and used genotyping by sequencing to sample over 600 bi-allelic markers. Our results suggest that the SA N. aenigmaticus and Mexican plants are a nested within a clade of sexual tropical populations. In the US populations, we confirm an extreme sex ratio and only contiguous US watersheds share genotypes. The phylogenetic analysis of SNP data resolves four clusters: Mexican populations, male plants (Little Pigeon and Pigeon river watersheds) and two clusters of female plants; one from the Little Tennessee and Hiwassee Rivers (TR) and the other from the Ocoee (TR) and Coosa (AR) Rivers. All clusters are highly differentiated (Fst values over 0.9). In addition, our individual assignment analyses and PCAs reflect the phylogenetic results grouping the SA samples in three clades and recovering males and female plants with high genetic differentiation (Fst values between 0.5 and 0.9 using microsatellites and bi-allelic markers). Our results point to Pleistocene events shaping the biogeographical pattern seen in US populations. The extreme sex ratio reflects isolation and highlights the high vulnerability of the populations in the SA.

Cophylogenetic patterns in algal symbionts correlate with repeated symbiont switches during diversification and geographic expansion of lichen-forming fungi in the genus Sticta (Ascomycota, Peltigeraceae).

Species in the fungal genus Sticta form symbiotic associations primarily with either green algae or cyanobacteria, but tripartite associations or photosymbiodemes involving both types of photobionts occur in some species. Sticta is known to associate with green algae in the genus Symbiochloris. However, previous studies have shown that algae from other genera, such as Heveochlorella, may also be suitable partners for Sticta. We examined the diversity of green algal partners in the genus Sticta and assessed the patterns of association between the host fungus and its algal symbiont. We used multi-locus sequence data from multiple individuals collected in Australia, Cuba, Madagascar, Mauritius, New Zealand, Reunion and South America to infer phylogenies for fungal and algal partners and performed tests of congruence to assess coevolution between the partners. In addition, event-based methods were implemented to examine which cophylogenetic processes have led to the observed association patterns in Sticta and its green algal symbionts. Our results show that in addition to Symbiochloris, Sticta associates with green algae from the genera Chloroidium, Coccomyxa, Elliptochloris and Heveochlorella, the latter being the most common algal symbiont associated with Sticta in this study. Geography plays a strong role in shaping fungal-algal association patterns in Sticta as mycobionts associate with different algal lineages in different geographic locations. While fungal and algal phylogenies were mostly congruent, event-based methods did not find any evidence for cospeciation between the partners. Instead, the association patterns observed in Sticta and associated algae, were largely explained by other cophylogenetic events such as host-switches, losses of symbiont and failure of the symbiont to diverge with its host. Our results also show that tripartite associations with green algae evolved multiple times in Sticta.

The Moss Physcomitrium (Physcomitrella) patens: A Model Organism for Non-Seed Plants.

Since the discovery two decades ago that transgenes are efficiently integrated into the genome of Physcomitrella patens by homologous recombination, this moss has been a premier model system to study evolutionary developmental biology questions, stem cell reprogramming, and the biology of nonvascular plants. P patens was the first non-seed plant to have its genome sequenced. With this level of genomic information, together with increasing molecular genetic tools, a large number of reverse genetic studies have propelled the use of this model system. A number of technological advances have recently opened the door to forward genetics as well as extremely efficient and precise genome editing in P patens Additionally, careful phylogenetic studies with increased resolution have suggested that P patens emerged from within Physcomitrium Thus, rather than Physcomitrella patens, the species should be named Physcomitrium patens Here we review these advances and describe the areas where P patens has had the most impact on plant biology.

The hornwort genome and early land plant evolution.

Hornworts, liverworts and mosses are three early diverging clades of land plants, and together comprise the bryophytes. Here, we report the draft genome sequence of the hornwort Anthoceros angustus. Phylogenomic inferences confirm the monophyly of bryophytes, with hornworts sister to liverworts and mosses. The simple morphology of hornworts correlates with low genetic redundancy in plant body plan, while the basic transcriptional regulation toolkit for plant development has already been established in this early land plant lineage. Although the Anthoceros genome is small and characterized by minimal redundancy, expansions are observed in gene families related to RNA editing, UV protection and desiccation tolerance. The genome of A. angustus bears the signatures of horizontally transferred genes from bacteria and fungi, in particular of genes operating in stress-response and metabolic pathways. Our study provides insight into the unique features of hornworts and their molecular adaptations to live on land.

Draft genome of the aquatic moss Fontinalis antipyretica (Fontinalaceae, Bryophyta).

Mosses comprise one of three lineages forming a sister group to extant vascular plants. Having emerged from an early split in the diversification of embryophytes, mosses may offer complementary insights into the evolution of traits following the transition to, and colonization of, land. Here, we report the draft nuclear genome of Fontinalis antipyretica (Fontinalaceae, Hypnales), a charismatic aquatic moss that is widespread in temperate regions of the Northern Hemisphere. We sequenced and de novo-assembled its genome using the 10X Genomics method. The genome comprises 385.2 Mbp, with a scaffold N50 of 45.8 Kbp. The assembly captured 87.2% of the 430 genes in the BUSCO Viridiplantae odb10 dataset. The newly generated F. antipyretica genome is the third moss genome, and the second seedless aquatic plant genome, to be sequenced and assembled to date.

Mitochondrial genomes of the early land plant lineage liverworts (Marchantiophyta): conserved genome structure, and ongoing low frequency recombination.

BACKGROUND: In contrast to the highly labile mitochondrial (mt) genomes of vascular plants, the architecture and composition of mt genomes within the main lineages of bryophytes appear stable and invariant. The available mt genomes of 18 liverwort accessions representing nine genera and five orders are syntenous except for Gymnomitrion concinnatum whose genome is characterized by two rearrangements. Here, we expanded the number of assembled liverwort mt genomes to 47, broadening the sampling to 31 genera and 10 orders spanning much of the phylogenetic breadth of liverworts to further test whether the evolution of the liverwort mitogenome is overall static. RESULTS: Liverwort mt genomes range in size from 147 Kb in Jungermanniales (clade B) to 185 Kb in Marchantiopsida, mainly due to the size variation of intergenic spacers and number of introns. All newly assembled liverwort mt genomes hold a conserved set of genes, but vary considerably in their intron content. The loss of introns in liverwort mt genomes might be explained by localized retroprocessing events. Liverwort mt genomes are strictly syntenous in genome structure with no structural variant detected in our newly assembled mt genomes. However, by screening the paired-end reads, we do find rare cases of recombination, which means multiple concurrent genome structures may exist in the vegetative tissues of liverworts. Our phylogenetic analyses of the nuclear encoded double stand break repair protein families revealed liverwort-specific subfamilies expansions. CONCLUSIONS: The low repeat recombination level, selection, along with the intensified nuclear surveillance, might together shape the structural evolution of liverwort mt genomes.

The Amount of RNA Editing Sites in Liverwort Organellar Genes Is Correlated with GC Content and Nuclear PPR Protein Diversity.

RNA editing occurs in the organellar mRNAs of all land plants but the marchantioid liverworts, making liverworts a perfect group for studying the evolution of RNA editing. Here, we profiled the RNA editing of 42 exemplars spanning the ordinal phylogenetic diversity of liverworts, and screened for the nuclear-encoded pentatricopeptide repeat (PPR) proteins in the transcriptome assemblies of these taxa. We identified 7,428 RNA editing sites in 128 organellar genes from 31 non-marchantioid liverwort species, and characterized 25,059 PPR protein sequences. The abundance of organellar RNA editing sites varies greatly among liverwort lineages, genes, and codon positions, and shows strong positive correlations with the GC content of protein-coding genes, and the diversity of the PLS class of nuclear PPR proteins.