Establishing CRISPR-Cas9 in the sexually dimorphic moss, Ceratodon purpureus.

The development of CRISPR technologies provides a powerful tool for understanding the evolution and functionality of essential biological processes. Here we demonstrate successful CRISPR-Cas9 genome editing in the dioecious moss species, Ceratodon purpureus. Using an existing selection system from the distantly related hermaphroditic moss, Physcomitrium patens, we generated knock-outs of the APT reporter gene by employing CRISPR-targeted mutagenesis under expression of native U6 snRNA promoters. Next, we used the native homology-directed repair (HDR) pathway, combined with CRISPR-Cas9, to knock in two reporter genes under expression of an endogenous RPS5A promoter in a newly developed landing site in C. purpureus. Our results show that the molecular tools developed in P. patens can be extended to other mosses across this ecologically important and developmentally variable group. These findings pave the way for precise and powerful experiments aimed at identifying the genetic basis of key functional variation within the bryophytes and between the bryophytes and other land plants.

Local adaptation, recombination, and the fate of neopolyploids.

Polyploidy is widely recognized as an important speciation mechanism because it isolates tetraploids from their diploid progenitors. Polyploidy also provides new genetic material that may facilitate adaptive evolution. However, new mutations are more likely to arise after a neopolyploid has already successfully invaded a population. Thus, the role of adaptive forces in establishing a polyploid remains unclear. One solution to this apparent paradox may lie in the capacity of polyploids to suppress recombination among preexisting locally adapted alleles. The local adaptation mechanism requires that spatially heterogeneous selection acts on multiple loci and that gene flow introduces maladapted alleles to the population where the polyploid forms. The mechanism requires neither strong genetic drift nor any intrinsic benefit of genome doubling and can accommodate any mode of gene action. A unique prediction of the mechanism is that adaptive alleles should predate polyploidization, a pattern consistent with observations from a few well-studied polyploids. The mechanism is also consistent with the coexistence of both diploid and tetraploid cytotypes, fitness heterogeneity among independently derived polyploids, and the prevalence of outcrossing among older polyploids. The local adaptation mechanism also makes novel predictions about circumstances favoring polyploid invasions that can be tested using molecular genetic or comparative approaches.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Multiple factors influence population sex ratios in the Mojave Desert moss Syntrichia caninervis.

PREMISE OF RESEARCH: Natural populations of many mosses appear highly female-biased based on the presence of reproductive structures. This bias could be caused by increased male mortality, lower male growth rate, or a higher threshold for achieving sexual maturity in males. Here we test these hypotheses using samples from two populations of the Mojave Desert moss Syntrichia caninervis. METHODS: We used double-digest restriction-site associated DNA (RAD) sequencing to identify candidate sex-associated loci in a panel of sex-expressing plants. Next, we used putative sex-associated markers to identify the sex of individuals without sex structures. KEY RESULTS: We found a 17:1 patch-level phenotypic female to male sex ratio in the higher elevation site (Wrightwood) and no sex expression at the low elevation site (Phelan). In contrast, on the basis of genetic data, we found a 2:1 female bias at the Wrightwood site and only females at the Phelan site. The relative area occupied by male and female genets was indistinguishable, but males were less genetically diverse. CONCLUSIONS: Our data suggest that both male-biased mortality and sexual dimorphism in thresholds for sex expression could explain genetic and phenotypic sex ratio biases and that phenotypic sex expression alone over-estimates the extent of actual sex ratio bias present in these two populations of S. caninervis.

High migration rates shape the postglacial history of amphi-Atlantic bryophytes.

Paleontological evidence and current patterns of angiosperm species richness suggest that European biota experienced more severe bottlenecks than North American ones during the last glacial maximum. How well this pattern fits other plant species is less clear. Bryophytes offer a unique opportunity to contrast the impact of the last glacial maximum in North America and Europe because about 60% of the European bryoflora is shared with North America. Here, we use population genetic analyses based on approximate Bayesian computation on eight amphi-Atlantic species to test the hypothesis that North American populations were less impacted by the last glacial maximum, exhibiting higher levels of genetic diversity than European ones and ultimately serving as a refugium for the postglacial recolonization of Europe. In contrast with this hypothesis, the best-fit demographic model involved similar patterns of population size contractions, comparable levels of genetic diversity and balanced migration rates between European and North American populations. Our results thus suggest that bryophytes have experienced comparable demographic glacial histories on both sides of the Atlantic. Although a weak, but significant genetic structure was systematically recovered between European and North American populations, evidence for migration from and towards both continents suggests that amphi-Atlantic bryophyte population may function as a metapopulation network. Reconstructing the biogeographic history of either North American or European bryophyte populations therefore requires a large, trans-Atlantic geographic framework.

Molecular evidence for convergent evolution and allopolyploid speciation within the Physcomitrium-Physcomitrella species complex.

BACKGROUND: The moss Physcomitrella patens (Hedw.) Bruch & Schimp. is an important experimental model system for evolutionary-developmental studies. In order to shed light on the evolutionary history of Physcomitrella and related species within the Funariaceae, we analyzed the natural genetic diversity of the Physcomitrium-Physcomitrella species complex. RESULTS: Molecular analysis of the nuclear single copy gene BRK1 reveals that three Physcomitrium species feature larger genome sizes than Physcomitrella patens and encode two expressed BRK1 homeologs (polyploidization-derived paralogs), indicating that they may be allopolyploid hybrids. Phylogenetic analyses of BRK1 as well as microsatellite simple sequence repeat (SSR) data confirm a polyphyletic origin for three Physcomitrella lineages. Differences in the conservation of mitochondrial editing sites further support hybridization and cryptic speciation within the Physcomitrium-Physcomitrella species complex. CONCLUSIONS: We propose a revised classification of the previously described four subspecies of Physcomitrella patens into three distinct species, namely Physcomitrella patens, Physcomitrella readeri and Physcomitrella magdalenae. We argue that secondary reduction of sporophyte complexity in these species is due to the establishment of an ecological niche, namely spores resting in mud and possible spore dispersal by migratory birds. Besides the Physcomitrium-Physcomitrella species complex, the Funariaceae are host to their type species, Funaria hygrometrica, featuring a sporophyte morphology which is more complex. Their considerable developmental variation among closely related lineages and remarkable trait evolution render the Funariaceae an interesting group for evolutionary and genetic research.

Are all sex chromosomes created equal?

Three principal types of chromosomal sex determination are found in nature: male heterogamety (XY systems, as in mammals), female heterogamety (ZW systems, as in birds), and haploid phase determination (UV systems, as in some algae and bryophytes). Although these systems share many common features, there are important biological differences between them that have broad evolutionary and genomic implications. Here we combine theoretical predictions with empirical observations to discuss how differences in selection, genetic properties and transmission uniquely shape each system. We elucidate how the differences among these systems can be exploited to gain insights about general evolutionary processes, genome structure, and gene expression. We suggest directions for research that will greatly increase our general understanding of the forces driving sex-chromosome evolution in diverse organisms.

Meiotic sex ratio variation in natural populations of Ceratodon purpureus (Ditrichaceae).

UNLABELLED: • PREMISE OF STUDY: Sex ratio variation is a common but often unexplained phenomenon in species across the tree of life. Here we evaluate the hypothesis that meiotic sex ratio variation can contribute to the biased sex ratios found in natural populations of the moss Ceratodon purpureus.• METHODS: We obtained sporophytes from several populations of C. purpureus from eastern North America. From each sporophyte, we estimated the mean spore viability by germinating replicate samples on agar plates. We estimated the meiotic sex ratio of each sporophyte by inferring the sex of a random sample of germinated spores (mean = 77) using a PCR-RFLP test. We tested for among-sporophyte variation in viability using an ANOVA and for deviations from 1:1 sex ratio using a χ(2)-test and evaluated the relationship between these quantities using a linear regression.• KEY RESULTS: We found among-sporophyte variation in spore viability and meiotic sex ratio, suggesting that genetic variants that contribute to variation in both of these traits segregate within populations of this species. However, we found no relationship between these quantities, suggesting that factors other than sex ratio distorters contribute to variation in spore viability within populations.• CONCLUSIONS: These results demonstrate that sex ratio distortion may partially explain the population sex ratio variation seen in C. purpureus, but more generally that genetic conflict over meiotic segregation may contribute to fitness variation in this species. Overall, this study lays the groundwork for future studies on the genetic basis of meiotic sex ratio variation.

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.

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.