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.

Divergent outcomes of genetic conflict on the UV sex chromosomes of Marchantia polymorpha and Ceratodon purpureus.

In species with separate sexes, the genome must produce two distinct developmental programs. Sexually dimorphic development may be controlled by either sex-limited loci or biased expression of loci transmitted through both sexes. Variation in the gene content of sex-limited chromosomes demonstrates that eukaryotic species differ markedly in the roles of these two mechanisms in governing sexual dimorphism. The bryophyte model systems Marchantia polymorpha and Ceratodon purpureus provide a particularly striking contrast. Although both species possess a haploid UV sex chromosome system, in which females carry a U chromosome and males carry a V, M. polymorpha relies on biased autosomal expression, while in C. purpureus, sex-linked genes drive dimorphism. Framing these genetic architectures as divergent outcomes of genetic conflict highlights comparative genomic analyses to better understand the evolution of sexual dimorphism.

Heavy metal tolerance in Scopelophila cataractae: Transcriptomic and epigenetic datasets.

Studying how different plant groups deal with heavy metal exposure is crucial to improve our understanding of the diversity of molecular mechanisms involved in plant stress response. Here, we used RNA sequencing (RNA-seq) and epigenotyping by sequencing (epiGBS) to assess gene expression and DNA methylation changes respectively in plants from four populations of the metallophyte moss Scopelophila cataractae treated with Cd or Cu in the laboratory. We built RNA-seq and epiGBS sequencing libraries from control and treated samples from each population and sequenced them using Illumina HiSeq 3000 (PE-150 bp) and Illumina HiSeq X-Ten System (PE-150 bp) respectively. For the RNA-seq data, we performed a read quality filter, mapped the reads to the de novo transcriptome created with Trinity, and estimated transcript abundance for each sample. For the epiGBS data, we used a custom pipeline (https://doi.org/10.5281/zenodo.7040291) to map the reads to a de novo reference genome and performed strand-specific nucleotide (single nucleotide polymorphisms, SNPs) and methylation (single cytosine methylation polymorphisms, SMPs) variant calling. We filtered out SNPs and SMPs with low coverage within (positions with <10 sequencing reads per sample) and across samples (positions with poor representation on the full set of samples). Finally, we performed pairwise comparisons between control and treated samples from each population and identified differentially expressed genes and differentially methylated cytosines associated to heavy metal exposure. We payed particular attention to the different responses of the more and the less tolerant populations of S. cataractae. These datasets could contribute to future comparative studies of abiotic stress response across plant groups.

Microevolutionary variation in molar morphology of Onychomys leucogaster decoupled from genetic structure.

In neutral models of quantitative trait evolution, both genetic and phenotypic divergence scale as random walks, producing a correlation between the two measures. However, complexity in the genotype-phenotype map may alter the correlation between genotypic and phenotypic divergence, even when both are evolving neutrally or nearly so. Understanding this correlation between phenotypic and genetic variation is critical for accurately interpreting the fossil record. This study compares the geographic structure and scaling of morphological variation of the shape of the first lower molar of 77 individuals of the northern grasshopper mouse Onychomys leucogaster to genome-wide SNP variation in the same sample. We found strong genetic structure but weak or absent morphological structure indicating that the scaling of each type of variation is decoupled from one another. Low PST values relative to FST values are consistent with a lack of morphological divergence in contrast to genetic divergence between groups. This lack of phenotypic structure and the presence of notable within-sample phenotypic variance are consistent with uniform selection or constraints on molar shape across a wide geographic and environmental range. Over time, this kind of decoupling may result in patterns of phenotypic stasis masking underlying genetic patterns.

An NGS-Based Phylogeny of Orthotricheae (Orthotrichaceae, Bryophyta) With the Proposal of the New Genus Rehubryum From Zealandia.

Phylogenomic data increase the possibilities of resolving the evolutionary and systematic relationships among taxa. This is especially valuable in groups with few and homoplasious morphological characters, in which systematic and taxonomical delimitations have been traditionally difficult. Such is the case of several lineages within Bryophyta, like Orthotrichaceae, the second most diverse family of mosses. Members of tribe Orthotricheae are common in temperate and cold regions, as well as in high tropical mountains. In extratropical areas, they represent one of the main components of epiphytic communities, both in dry and oceanic or hyperoceanic conditions. The epiphytic environment is considered a hostile one for plant development, mainly due to its low capacity of moisture retention. Thus, the diversification of the Orthotrichaceae in this environment could be seen as striking. Over the last two decades, great taxonomic and systematic progresses have led to a rearrangement at the generic level in this tribe, providing a new framework to link environment to patterns of diversification. Here, we use nuclear loci targeted with the GoFlag 408 enrichment probe set to generate a well-sampled phylogeny with well-supported suprageneric taxa and increasing the phylogenetic resolution within the two recognized subtribes. Specifically, we show that several genera with Ulota-like morphology jointly constitute an independent lineage. Within this lineage, the recently described Atlantichella from Macaronesia and Western Europe appears as the sister group of Ulota bellii from Zealandia. This latter species is here segregated in the new genus Rehubryum. Assessment of the ecological and biogeographical affinities of the species within the phylogenetic framework suggests that niche adaptation (including climate and substrate) may be a key evolutionary driver that shaped the high diversification of Orthotricheae.

Measuring volatile emissions from moss gametophytes: A review of methodologies and new applications.

Mosses inhabit nearly all terrestrial ecosystems and engage in important interactions with nitrogen-fixing microbes, sperm-dispersing arthropods, and other plants. It is hypothesized that these interactions could be mediated by biogenic volatile organic compounds (BVOCs). Moss BVOCs may play fundamental roles in influencing local ecologies, such as biosphere-atmosphere-hydrosphere communications, physiological and evolutionary dynamics, plant-microbe interactions, and gametophyte stress physiology. Further progress in quantifying the composition, magnitude, and variability of moss BVOC emissions, and their response to environmental drivers and metabolic requirements, is limited by methodological and analytical challenges. We review several sampling techniques with various analytical approaches and describe best practices in generating moss gametophyte BVOC measures. We emphasize the importance of characterizing the composition and magnitude of moss BVOC emissions across a variety of species to better inform and stimulate important cross-disciplinary studies. We conclude by highlighting how current methods could be employed, as well as best practices for choosing methodologies.

Gene-rich UV sex chromosomes harbor conserved regulators of sexual development.

Nonrecombining sex chromosomes, like the mammalian Y, often lose genes and accumulate transposable elements, a process termed degeneration. The correlation between suppressed recombination and degeneration is clear in animal XY systems, but the absence of recombination is confounded with other asymmetries between the X and Y. In contrast, UV sex chromosomes, like those found in bryophytes, experience symmetrical population genetic conditions. Here, we generate nearly gapless female and male chromosome-scale reference genomes of the moss Ceratodon purpureus to test for degeneration in the bryophyte UV sex chromosomes. We show that the moss sex chromosomes evolved over 300 million years ago and expanded via two chromosomal fusions. Although the sex chromosomes exhibit weaker purifying selection than autosomes, we find that suppressed recombination alone is insufficient to drive degeneration. Instead, the U and V sex chromosomes harbor thousands of broadly expressed genes, including numerous key regulators of sexual development across land plants.

The relationship of C and N stable isotopes to high-latitude moss-associated N2 fixation.

Moss-associated N2 fixation by epiphytic microbes is a key biogeochemical process in nutrient-limited high-latitude ecosystems. Abiotic drivers, such as temperature and moisture, and the identity of host mosses are critical sources of variation in N2 fixation rates. An understanding of the potential interaction between these factors is essential for predicting N inputs as moss communities change with the climate. To further understand the drivers and results of N2 fixation rate variation, we obtained natural abundance values of C and N isotopes and an associated rate of N2 fixation with 15N2 gas incubations in 34 moss species collected in three regions across Alaska, USA. We hypothesized that δ15N values would increase toward 0‰ with higher N2 fixation to reflect the increasing contribution of fixed N2 in moss biomass. Second, we hypothesized that δ13C and N2 fixation would be positively related, as enriched δ13C signatures reflect abiotic conditions favorable to N2 fixation. We expected that the magnitude of these relationships would vary among types of host mosses, reflecting differences in anatomy and habitat. We found little support for our first hypothesis, with only a modest positive relationship between N2 fixation rates and δ15N in a structural equation model. We found a significant positive relationship between δ13C and N2 fixation only in Hypnales, where the probability of N2 fixation activity reached 95% when δ13C values exceeded - 30.4‰. We conclude that moisture and temperature interact strongly with host moss identity in determining the extent to which abiotic conditions impact associated N2 fixation rates.

Area from image analyses accurately estimates dry-weight biomass of juvenile tissue from the moss Ceratodon purpureus.

PREMISE: Mosses have long served as models for studying many areas of plant biology. Investigators have used two-dimensional measurements of juvenile growth from photographs as a surrogate for dry-weight biomass. The relationship between area and biomass, however, has not been critically evaluated. METHODS: Here we grew axenic tissue cultures of 10 Ceratodon purpureus isolates to study the relationship between these parameters. We measured area and biomass on replicate cultures with two distinct starting inoculum sizes each week for three weeks. We then examined the correlation between area and biomass as well as the influence of variation in inoculum size on both parameters. RESULTS: We found a strong correlation between area and biomass after two weeks of growth. Furthermore, we found inoculum size affected biomass during the first week of growth but not in subsequent weeks and inoculum size had no detectable effect on area. DISCUSSION: These analyses provide experimental confirmation that area is a suitable proxy for biomass and provide clear guidelines for when inoculum size variation may affect downstream growth estimates.

The genetic architecture of sexual dimorphism in the moss Ceratodon purpureus.

A central problem in evolutionary biology is to identify the forces that maintain genetic variation for fitness in natural populations. Sexual antagonism, in which selection favours different variants in males and females, can slow the transit of a polymorphism through a population or can actively maintain fitness variation. The amount of sexually antagonistic variation to be expected depends in part on the genetic architecture of sexual dimorphism, about which we know relatively little. Here, we used a multivariate quantitative genetic approach to examine the genetic architecture of sexual dimorphism in a scent-based fertilization syndrome of the moss Ceratodon purpureus. We found sexual dimorphism in numerous traits, consistent with a history of sexually antagonistic selection. The cross-sex genetic correlations (rmf) were generally heterogeneous with many values indistinguishable from zero, which typically suggests that genetic constraints do not limit the response to sexually antagonistic selection. However, we detected no differentiation between the female- and male-specific trait (co)variance matrices (Gf and Gm, respectively), meaning the evolution of sexual dimorphism may be constrained. The cross-sex cross-trait covariance matrix B contained both symmetric and asymmetric elements, indicating that the response to sexually antagonistic or sexually concordant selection, and the constraint to sexual dimorphism, are highly dependent on the traits experiencing selection. The patterns of genetic variances and covariances among these fitness components is consistent with partly sex-specific genetic architectures having evolved in order to partially resolve multivariate genetic constraints (i.e. sexual conflict), enabling the sexes to evolve towards their sex-specific multivariate trait optima.

Microarthropod contributions to fitness variation in the common moss Ceratodon purpureus.

The evolution of sustained plant-animal interactions depends critically upon genetic variation in the fitness benefits from the interaction. Genetic analyses of such interactions are limited to a few model systems, in part because genetic variation may be absent or the interacting species may be experimentally intractable. Here, we examine the role of sperm-dispersing microarthropods in shaping reproduction and genetic variation in mosses. We established experimental mesocosms with known moss genotypes and inferred the parents of progeny from mesocosms with and without microarthropods, using a pooled sequencing approach. Moss reproductive rates increased fivefold in the presence of microarthropods, relative to control mesocosms. Furthermore, the presence of microarthropods increased the total number of reproducing moss genotypes, and changed the rank-order of fitness of male and female moss genotypes. Interestingly, the genotypes that reproduced most frequently did not produce sporophytes with the most spores, highlighting the challenge of defining fitness in mosses. These results demonstrate that microarthropods provide a fitness benefit for mosses, and highlight the potential for biotic dispersal agents to alter fitness among moss genotypes.

A target enrichment probe set for resolving the flagellate land plant tree of life.

PREMISE: New sequencing technologies facilitate the generation of large-scale molecular data sets for constructing the plant tree of life. We describe a new probe set for target enrichment sequencing to generate nuclear sequence data to build phylogenetic trees with any flagellate land plants, including hornworts, liverworts, mosses, lycophytes, ferns, and all gymnosperms. METHODS: We leveraged existing transcriptome and genome sequence data to design the GoFlag 451 probes, a set of 56,989 probes for target enrichment sequencing of 451 exons that are found in 248 single-copy or low-copy nuclear genes across flagellate plant lineages. RESULTS: Our results indicate that target enrichment using the GoFlag451 probe set can provide large nuclear data sets that can be used to resolve relationships among both distantly and closely related taxa across the flagellate land plants. We also describe the GoFlag 408 probes, an optimized probe set covering 408 of the 451 exons from the GoFlag 451 probe set that is commercialized by RAPiD Genomics. CONCLUSIONS: A target enrichment approach using the new probe set provides a relatively low-cost solution to obtain large-scale nuclear sequence data for inferring phylogenetic relationships across flagellate land plants.

The bacterial communities of Alaskan mosses and their contributions to N2-fixation.

BACKGROUND: Mosses in high-latitude ecosystems harbor diverse bacterial taxa, including N2-fixers which are key contributors to nitrogen dynamics in these systems. Yet the relative importance of moss host species, and environmental factors, in structuring these microbial communities and their N2-fixing potential remains unclear. We studied 26 boreal and tundra moss species across 24 sites in Alaska, USA, from 61 to 69° N. We used cultivation-independent approaches to characterize the variation in moss-associated bacterial communities as a function of host species identity and site characteristics. We also measured N2-fixation rates via 15N2 isotopic enrichment and identified potential N2-fixing bacteria using available literature and genomic information. RESULTS: Host species identity and host evolutionary history were both highly predictive of moss microbiome composition, highlighting strong phylogenetic coherence in these microbial communities. Although less important, light availability and temperature also influenced composition of the moss microbiome. Finally, we identified putative N2-fixing bacteria specific to some moss hosts, including potential N2-fixing bacteria outside well-studied cyanobacterial clades. CONCLUSIONS: The strong effect of host identity on moss-associated bacterial communities demonstrates mosses' utility for understanding plant-microbe interactions in non-leguminous systems. Our work also highlights the likely importance of novel bacterial taxa to N2-fixation in high-latitude ecosystems. Video Abstract.

Bryophytes are not early diverging land plants.

Phylogenetic trees have permeated biology. However, an understanding of how to interpret phylogenies has lagged behind, notably in publications outside of evolutionary biology. Here I argue that some language commonly used in plant systematics has contributed to the confusion by describing phylogenetic trees using intuitive but misleading terms reminiscent of Aristotle's Scala Naturae. These terms (perhaps inadvertently) misrepresent evolution, not as a process acting on all living species, but rather as a progression of successively diverging lineages leading to a group that represents a subjectively defined endpoint. My goal here is to show how thinking of the tree of life in terms of early-diverging lineages and higher groups can distort evolutionary literacy, confound interdisciplinary communication, and potentially bias research agendas. I focus on the relationship between bryophytes and angiosperms as a case study, but the theme applies to all branches of the tree of life. Fortunately, evolutionary biologists have developed an easily understood alternative framework - tree thinking - which I highlight as a means to promote a clear understanding of phylogenies across sub-disciplines of biology, and between practicing biologists and students, or members the public which funds much of our work.

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

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

Hitchhiking the high seas: Global genomics of rafting crabs.

Population differentiation and diversification depend in large part on the ability and propensity of organisms to successfully disperse. However, our understanding of these processes in organisms with high dispersal ability is biased by the limited genetic resolution offered by traditional genotypic markers. Many neustonic animals disperse not only as pelagic larvae, but also as juveniles and adults while drifting or rafting at the surface of the open ocean. In theory, the heightened dispersal ability of these animals should limit opportunities for species diversification and population differentiation. To test these predictions, we used next-generation sequencing of genomewide restriction-site-associated DNA tags (RADseq) and traditional mitochondrial DNA sequencing, to investigate the species-level relationships and global population structure of Planes crabs collected from oceanic flotsam and sea turtles. Our results indicate that species diversity in this clade is low-likely three closely related species-with no evidence of cryptic or undescribed species. Moreover, our results indicate weak population differentiation among widely separated aggregations with genetic indices showing only subtle genetic discontinuities across all oceans of the world (RADseq F ST = 0.08-0.16). The results of this study provide unprecedented resolution of the systematics and global biogeography of this group and contribute valuable information to our understanding of how theoretical dispersal potential relates to actual population differentiation and diversification among marine organisms. Moreover, these results demonstrate the limitations of single gene analyses and the value of genomic-level resolution for estimating contemporary population structure in organisms with large, highly connected populations.

Terrestrial species adapted to sea dispersal: Differences in propagule dispersal of two Caribbean mangroves.

A central goal of comparative phylogeography is to understand how species-specific traits interact with geomorphological history to govern the geographic distribution of genetic variation within species. One key biotic trait with an immense impact on the spatial patterns of intraspecific genetic differentiation is dispersal. Here, we quantify how species-specific traits directly related to dispersal affect genetic variation in terrestrial organisms with adaptations for dispersal by sea, not land-the mangroves of the Caribbean. We investigate the phylogeography of white mangroves (Laguncularia racemosa, Combretaceae) and red mangroves (Rhizophora mangle, Rhizophoraceae) using chloroplast genomes and nuclear markers (thousands of RAD-Seq loci) from individuals throughout the Caribbean. Both coastal tree species have viviparous propagules that can float in salt water for months, meaning they are capable of dispersing long distances. Spatially explicit tests of the role of ocean currents on patterning genetic diversity revealed that ocean currents act as a mechanism for facilitating dispersal, but other means of moving genetic material are also important. We measured pollen- vs. propagule-mediated gene flow and discovered that in white mangroves, seeds were more important for promoting genetic connectivity between populations, but in red mangroves, the opposite was true: pollen contributed more. This result challenges our concept of the importance of both proximity to ocean currents for moving mangrove seeds and the extent of long-distance pollen dispersal. This study also highlights the importance of spatially explicit quantification of both abiotic (ocean currents) and biotic (dispersal) factors contributing to gene flow to understand fully the phylogeographic histories of species.

Novel bacterial lineages associated with boreal moss species.

Mosses are critical components of boreal ecosystems where they typically account for a large proportion of net primary productivity and harbour diverse bacterial communities that can be the major source of biologically-fixed nitrogen in these ecosystems. Despite their ecological importance, we have limited understanding of how microbial communities vary across boreal moss species and the extent to which local site conditions may influence the composition of these bacterial communities. We used marker gene sequencing to analyze bacterial communities associated with seven boreal moss species collected near Fairbanks, AK, USA. We found that host identity was more important than site in determining bacterial community composition and that mosses harbour diverse lineages of potential N2 -fixers as well as an abundance of novel taxa assigned to understudied bacterial phyla (including candidate phylum WPS-2). We performed shotgun metagenomic sequencing to assemble genomes from the WPS-2 candidate phylum and found that these moss-associated bacteria are likely anoxygenic phototrophs capable of carbon fixation via RuBisCo with an ability to utilize byproducts of photorespiration from hosts via a glyoxylate shunt. These results give new insights into the metabolic capabilities of understudied bacterial lineages that associate with mosses and the importance of plant hosts in shaping their microbiomes.

Peripatric speciation associated with genome expansion and female-biased sex ratios in the moss genus Ceratodon.

PREMISE OF THE STUDY: A period of allopatry is widely believed to be essential for the evolution of reproductive isolation. However, strict allopatry may be difficult to achieve in some cosmopolitan, spore-dispersed groups, like mosses. We examined the genetic and genome size diversity in Mediterranean populations of the moss Ceratodon purpureus s.l. to evaluate the role of allopatry and ploidy change in population divergence. METHODS: We sampled populations of the genus Ceratodon from mountainous areas and lowlands of the Mediterranean region, and from Western and Central Europe. We performed phylogenetic and coalescent analyses on sequences from five nuclear introns and a chloroplast locus to reconstruct their evolutionary history. We also estimated genome size using flow cytometry (employing propidium iodide) and determined the sex of samples using a sex-linked PCR marker. KEY RESULTS: Two well-differentiated clades were resolved, discriminating two homogeneous groups: the widespread C. purpureus and a local group mostly restricted to the mountains in Southern Spain. The latter also possessed a genome size 25% larger than the widespread C. purpureus, and the samples of this group consist entirely of females. We also found hybrids, and some of them had a genome size equivalent to the sum of the C. purpureus and Spanish genome, suggesting that they arose by allopolyploidy. CONCLUSIONS: These data suggest that a new species of Ceratodon arose via peripatric speciation, potentially involving a genome size change and a strong female-biased sex ratio. The new species has hybridized in the past with C. purpureus.