Divergent selection and climate adaptation fuel genomic differentiation between sister species of Sphagnum (peat moss).

BACKGROUND AND AIMS: New plant species can evolve through the reinforcement of reproductive isolation via local adaptation along habitat gradients. Peat mosses (Sphagnaceae) are an emerging model system for the study of evolutionary genomics and have well-documented niche differentiation among species. Recent molecular studies have demonstrated that the globally distributed species Sphagnum magellanicum is a complex of morphologically cryptic lineages that are phylogenetically and ecologically distinct. Here, we describe the architecture of genomic differentiation between two sister species in this complex known from eastern North America: the northern S. diabolicum and the largely southern S. magniae. METHODS: We sampled plant populations from across a latitudinal gradient in eastern North America and performed whole genome and restriction-site associated DNA sequencing. These sequencing data were then analyzed computationally. KEY RESULTS: Using sliding-window population genetic analyses we find that differentiation is concentrated within 'islands' of the genome spanning up to 400 kb that are characterized by elevated genetic divergence, suppressed recombination, reduced nucleotide diversity and increased rates of non-synonymous substitution. Sequence variants that are significantly associated with genetic structure and bioclimatic variables occur within genes that have functional enrichment for biological processes including abiotic stress response, photoperiodism and hormone-mediated signalling. Demographic modelling demonstrates that these two species diverged no more than 225 000 generations ago with secondary contact occurring where their ranges overlap. CONCLUSIONS: We suggest that this heterogeneity of genomic differentiation is a result of linked selection and reflects the role of local adaptation to contrasting climatic zones in driving speciation. This research provides insight into the process of speciation in a group of ecologically important plants and strengthens our predictive understanding of how plant populations will respond as Earth's climate rapidly changes.

Clonality, local population structure and gametophyte sex ratios in cryptic species of the Sphagnum magellanicum complex.

BACKGROUND AND AIMS: Sphagnum (peatmoss) comprises a moss (Bryophyta) clade with ~300-500 species. The genus has unparalleled ecological importance because Sphagnum-dominated peatlands store almost a third of the terrestrial carbon pool and peatmosses engineer the formation and microtopography of peatlands. Genomic resources for Sphagnum are being actively expanded, but many aspects of their biology are still poorly known. Among these are the degree to which Sphagnum species reproduce asexually, and the relative frequencies of male and female gametophytes in these haploid-dominant plants. We assess clonality and gametophyte sex ratios and test hypotheses about the local-scale distribution of clones and sexes in four North American species of the S. magellanicum complex. These four species are difficult to distinguish morphologically and are very closely related. We also assess microbial communities associated with Sphagnum host plant clones and sexes at two sites. METHODS: Four hundred and five samples of the four species, representing 57 populations, were subjected to restriction site-associated DNA sequencing (RADseq). Analyses of population structure and clonality based on the molecular data utilized both phylogenetic and phenetic approaches. Multi-locus genotypes (genets) were identified using the RADseq data. Sexes of sampled ramets were determined using a molecular approach that utilized coverage of loci on the sex chromosomes after the method was validated using a sample of plants that expressed sex phenotypically. Sex ratios were estimated for each species, and populations within species. Difference in fitness between genets was estimated as the numbers of ramets each genet comprised. Degrees of clonality [numbers of genets/numbers of ramets (samples)] within species, among sites, and between gametophyte sexes were estimated. Sex ratios were estimated for each species, and populations within species. Sphagnum-associated microbial communities were assessed at two sites in relation to Sphagnum clonality and sex. KEY RESULTS: All four species appear to engage in a mixture of sexual and asexual (clonal) reproduction. A single ramet represents most genets but two to eight ramets were dsumbers ansd text etected for some genets. Only one genet is represented by ramets in multiple populations; all other genets are restricted to a single population. Within populations ramets of individual genets are spatially clustered, suggesting limited dispersal even within peatlands. Sex ratios are male-biased in S. diabolicum but female-biased in the other three species, although significantly so only in S. divinum. Neither species nor males/females differ in levels of clonal propagation. At St Regis Lake (NY) and Franklin Bog (VT), microbial community composition is strongly differentiated between the sites, but differences between species, genets and sexes were not detected. Within S. divinum, however, female gametophytes harboured two to three times the number of microbial taxa as males. CONCLUSIONS: These four Sphagnum species all exhibit similar reproductive patterns that result from a mixture of sexual and asexual reproduction. The spatial patterns of clonally replicated ramets of genets suggest that these species fall between the so-called phalanx patterns, where genets abut one another but do not extensively mix because of limited ramet fragmentation, and the guerrilla patterns, where extensive genet fragmentation and dispersal result in greater mixing of different genets. Although sex ratios in bryophytes are most often female-biased, both male and female biases occur in this complex of closely related species. The association of far greater microbial diversity for female gametophytes in S. divinum, which has a female-biased sex ratio, suggests additional research to determine if levels of microbial diversity are consistently correlated with differing patterns of sex ratio biases.

Extensive Genome-Wide Phylogenetic Discordance Is Due to Incomplete Lineage Sorting and Not Ongoing Introgression in a Rapidly Radiated Bryophyte Genus.

The relative importance of introgression for diversification has long been a highly disputed topic in speciation research and remains an open question despite the great attention it has received over the past decade. Gene flow leaves traces in the genome similar to those created by incomplete lineage sorting (ILS), and identification and quantification of gene flow in the presence of ILS is challenging and requires knowledge about the true phylogenetic relationship among the species. We use whole nuclear, plastid, and organellar genomes from 12 species in the rapidly radiated, ecologically diverse, actively hybridizing genus of peatmoss (Sphagnum) to reconstruct the species phylogeny and quantify introgression using a suite of phylogenomic methods. We found extensive phylogenetic discordance among nuclear and organellar phylogenies, as well as across the nuclear genome and the nodes in the species tree, best explained by extensive ILS following the rapid radiation of the genus rather than by postspeciation introgression. Our analyses support the idea of ancient introgression among the ancestral lineages followed by ILS, whereas recent gene flow among the species is highly restricted despite widespread interspecific hybridization known in the group. Our results contribute to phylogenomic understanding of how speciation proceeds in rapidly radiated, actively hybridizing species groups, and demonstrate that employing a combination of diverse phylogenomic methods can facilitate untangling complex phylogenetic patterns created by ILS and introgression.

Habitat-adapted microbial communities mediate Sphagnum peatmoss resilience to warming.

Sphagnum peatmosses are fundamental members of peatland ecosystems, where they contribute to the uptake and long-term storage of atmospheric carbon. Warming threatens Sphagnum mosses and is known to alter the composition of their associated microbiome. Here, we use a microbiome transfer approach to test if microbiome thermal origin influences host plant thermotolerance. We leveraged an experimental whole-ecosystem warming study to collect field-grown Sphagnum, mechanically separate the associated microbiome and then transfer onto germ-free laboratory Sphagnum for temperature experiments. Host and microbiome dynamics were assessed with growth analysis, Chla fluorescence imaging, metagenomics, metatranscriptomics and 16S rDNA profiling. Microbiomes originating from warming field conditions imparted enhanced thermotolerance and growth recovery at elevated temperatures. Metagenome and metatranscriptome analyses revealed that warming altered microbial community structure in a manner that induced the plant heat shock response, especially the HSP70 family and jasmonic acid production. The heat shock response was induced even without warming treatment in the laboratory, suggesting that the warm-microbiome isolated from the field provided the host plant with thermal preconditioning. Our results demonstrate that microbes, which respond rapidly to temperature alterations, can play key roles in host plant growth response to rapidly changing environments.

Phylogenomic structure and speciation in an emerging model: the Sphagnum magellanicum complex (Bryophyta).

Sphagnum magellanicum is one of two Sphagnum species for which a reference-quality genome exists to facilitate research in ecological genomics. Phylogenetic and comparative genomic analyses were conducted based on resequencing data from 48 samples and RADseq analyses based on 187 samples. We report herein that there are four clades/species within the S. magellanicum complex in eastern North America and that the reference genome belongs to Sphagnum divinum. The species exhibit tens of thousands (RADseq) to millions (resequencing) of fixed nucleotide differences. Two species, however, referred to informally as S. diabolicum and S. magni because they have not been formally described, are differentiated by only 100 (RADseq) to 1000 (resequencing) of differences. Introgression among species in the complex is demonstrated using D-statistics and f4 ratios. One ecologically important functional trait, tissue decomposability, which underlies peat (carbon) accumulation, does not differ between segregates in the S. magellanicum complex, although previous research showed that many closely related Sphagnum species have evolved differences in decomposability/carbon sequestration. Phylogenetic resolution and more accurate species delimitation in the S. magellanicum complex substantially increase the value of this group for studying the early evolutionary stages of climate adaptation and ecological evolution more broadly.

Natural selection on a carbon cycling trait drives ecosystem engineering by Sphagnum (peat moss).

Sphagnum peat mosses have an extraordinary impact on the global carbon cycle as they control long-term carbon sequestration in boreal peatland ecosystems. Sphagnum species engineer peatlands, which harbour roughly a quarter of all terrestrial carbon, through peat accumulation by constructing their own niche that allows them to outcompete other plants. Interspecific variation in peat production, largely resulting from differences in tissue decomposability, is hypothesized to drive niche differentiation along microhabitat gradients thereby alleviating competitive pressure. However, little empirical evidence exists for the role of selection in the creation and maintenance of such gradients. In order to document how niche construction and differentiation evolved in Sphagnum, we quantified decomposability for 54 species under natural conditions and used phylogenetic comparative methods to model the evolution of this carbon cycling trait. We show that decomposability tracks the phylogenetic diversification of peat mosses, that natural selection favours different levels of decomposability corresponding to optimum niche and that divergence in this trait occurred early in the evolution of the genus prior to the divergence of most extant species. Our results demonstrate the evolution of ecosystem engineering via natural selection on an extended phenotype, of a fundamental ecosystem process, and one of the Earth's largest soil carbon pools.

Phylogenomics reveals convergent evolution of red-violet coloration in land plants and the origins of the anthocyanin biosynthetic pathway.

The flavonoids, one of the largest classes of plant secondary metabolites, are found in lineages that span the land plant phylogeny and play important roles in stress responses and as pigments. Perhaps the most well-studied flavonoids are the anthocyanins that have human health benefits and help plants attract pollinators, regulate hormone production, and confer resistance to abiotic and biotic stresses. The canonical biochemical pathway responsible for the production of these pigments is well-characterized for flowering plants yet its conservation across deep divergences in land plants remains debated and poorly understood. Many early land plants such as mosses, liverworts, and ferns produce flavonoid pigments, but their biosynthetic origins and homologies to the anthocyanin pathway remain uncertain. We conducted phylogenetic analyses using full genome sequences representing nearly all major green plant lineages to reconstruct the evolutionary history of the anthocyanin biosynthetic pathway then test the hypothesis that genes in this pathway are present in early land plants. We found that the entire pathway was not intact until the most recent common ancestor of seed plants and that orthologs of many downstream enzymes are absent from seedless plants including mosses, liverworts, and ferns. Our results also highlight the utility of phylogenetic inference, as compared to pairwise sequence similarity, in orthology assessment within large gene families that have complex duplication-loss histories. We suggest that the production of red-violet flavonoid pigments widespread in seedless plants, including the 3-deoxyanthocyanins, requires the activity of novel, as-yet discovered enzymes, and represents convergent evolution of red-violet coloration across land plants.

Phylogenetic structure in the Sphagnum recurvum complex (Bryophyta) in relation to taxonomy and geography.

PREMISE: The Sphagnum recurvum complex comprises a group of closely related peat mosses that are dominant components of many northern wetland ecosystems. Taxonomic hypotheses for the group range from interpreting the whole complex as one polymorphic species to distinguishing 6-10 species. The complex occurs throughout the Northern Hemisphere, and some of the putative species have intercontinental ranges. Our goals were to delimit the complex and assess its phylogenetic structure in relation to morphologically defined species and intercontinental geography. METHODS: RADseq analyses were applied to a sample of 384 collections from Europe, North America, and Asia. The data were subjected to maximum likelihood phylogenetic analyses and analyses of genetic structure using the software STRUCTURE and multivariate ordination approaches. RESULTS: The S. recurvum complex includes S. angustifolium, S. fallax, S. flexuosum, S. pacificum, and S. recurvum as clades with little evidence of admixture. We also resolved an unnamed clade that is referred to here as S. "pseudopacificum." We confirm that S. balticum and S. obtusum are nested within the complex. Species with bluntly acute to obtuse stem leaf apices are sister to those with acute to apiculate leaves. Most of the species exhibit some differentiation between intraspecific population systems disjunct on different continents. CONCLUSIONS: We recognize seven species in the amended S. recurvum complex, including S. balticum and S. obtusum, in addition to the informal clade S. "pseudopacificum." Although we detected some geographically correlated phylogenetic structure within widespread morphospecies, our RADseq data support the interpretation that these species have intercontinental geographic ranges.

Organellomic data sets confirm a cryptic consensus on (unrooted) land-plant relationships and provide new insights into bryophyte molecular evolution.

PREMISE: Phylogenetic trees of bryophytes provide important evolutionary context for land plants. However, published inferences of overall embryophyte relationships vary considerably. We performed phylogenomic analyses of bryophytes and relatives using both mitochondrial and plastid gene sets, and investigated bryophyte plastome evolution. METHODS: We employed diverse likelihood-based analyses to infer large-scale bryophyte phylogeny for mitochondrial and plastid data sets. We tested for changes in purifying selection in plastid genes of a mycoheterotrophic liverwort (Aneura mirabilis) and a putatively mycoheterotrophic moss (Buxbaumia), and compared 15 bryophyte plastomes for major structural rearrangements. RESULTS: Overall land-plant relationships conflict across analyses, generally weakly. However, an underlying (unrooted) four-taxon tree is consistent across most analyses and published studies. Despite gene coverage patchiness, relationships within mosses, liverworts, and hornworts are largely congruent with previous studies, with plastid results generally better supported. Exclusion of RNA edit sites restores cases of unexpected non-monophyly to monophyly for Takakia and two hornwort genera. Relaxed purifying selection affects multiple plastid genes in mycoheterotrophic Aneura but not Buxbaumia. Plastid genome structure is nearly invariant across bryophytes, but the tufA locus, presumed lost in embryophytes, is unexpectedly retained in several mosses. CONCLUSIONS: A common unrooted tree underlies embryophyte phylogeny, [(liverworts, mosses), (hornworts, vascular plants)]; rooting inconsistency across studies likely reflects substantial distance to algal outgroups. Analyses combining genomic and transcriptomic data may be misled locally for heavily RNA-edited taxa. The Buxbaumia plastome lacks hallmarks of relaxed selection found in mycoheterotrophic Aneura. Autotrophic bryophyte plastomes, including Buxbaumia, hardly vary in overall structure.

Functional trait evolution in Sphagnum peat mosses and its relationship to niche construction.

Species in the genus Sphagnum create, maintain, and dominate boreal peatlands through 'extended phenotypes' that allow these organisms to engineer peatland ecosystems and thereby impact global biogeochemical cycles. One such phenotype is the production of peat, or incompletely decomposed biomass, that accumulates when rates of growth exceed decomposition. Interspecific variation in peat production is thought to be responsible for the establishment and maintenance of ecological gradients such as the microtopographic hummock-hollow gradient, along which sympatric species sort within communities. This study investigated the mode and tempo of functional trait evolution across 15 species of Sphagnum using data from the most extensive studies of Sphagnum functional traits to date and phylogenetic comparative methods. We found evidence for phylogenetic conservatism of the niche descriptor height-above-water-table and of traits related to growth, decay and litter quality. However, we failed to detect the influence of phylogeny on interspecific variation in other traits such as shoot density and suggest that environmental context can obscure phylogenetic signal. Trait correlations indicate possible adaptive syndromes that may relate to niche and its construction. This study is the first to formally test the extent to which functional trait variation among Sphagnum species is a result of shared evolutionary history.

Range change evolution of peat mosses (Sphagnum) within and between climate zones.

Peat mosses (Sphagnum) hold exceptional importance in the control of global carbon fluxes and climate because of the vast stores of carbon bound up in partially decomposed biomass (peat). This study tests the hypothesis that the early diversification of Sphagnum was in the Northern Hemisphere, with subsequent range expansions to tropical latitudes and the Southern Hemisphere. A phylogenetic analysis of 192 accessions representing the moss class Sphagnopsida based on four plastid loci was conducted in conjunction with biogeographic analyses using BioGeoBEARS to investigate the tempo and mode of geographic range evolution. Analyses support the hypothesis that the major intrageneric clades of peat-forming species accounting for >90% of peat moss diversity originated and diversified at northern latitudes. The genus underwent multiple range expansions into tropical and Southern Hemisphere regions. Range evolution in peat mosses was most common within latitudinal zones, attesting to the relative difficulty of successfully invading new climate zones. Allopolyploidy in Sphagnum (inferred from microsatellite heterozygosity) does not appear to be biased with regard to geographic region nor intrageneric clade. The inference that Sphagnum diversified in cool-or cold-climate regions and repeatedly expanded its range into tropical regions makes the genus an excellent model for studying morphological, physiological, and genomic traits associated with adaptation to warming climates.

A novel experimental system using the liverwort Marchantia polymorpha and its fungal endophytes reveals diverse and context-dependent effects.

Fungal symbioses are ubiquitous in plants, but their effects have mostly been studied in seed plants. This study aimed to assess the diversity of fungal endophyte effects in a bryophyte and identify factors contributing to the variability of outcomes in these interactions. Fungal endophyte cultures and axenic liverwort clones were isolated from wild populations of the liverwort, Marchantia polymorpha. These collections were combined in a gnotobiotic system to test the effects of fungal isolates on the growth rates of hosts under laboratory conditions. Under the experimental conditions, fungi isolated from M. polymorpha ranged from aggressively pathogenic to strongly growth-promoting, but the majority of isolates caused no detectable change in host growth. Growth promotion by selected fungi depended on nutrient concentrations and was inhibited by coinoculation with multiple fungi. The M. polymorpha endophyte system expands the resources for this model liverwort. The experiments presented here demonstrate a wealth of diversity in fungal interactions even in a host reported to lack standard mycorrhizal symbiosis. In addition, they show that some known pathogens of vascular plants live in M. polymorpha and can confer benefits to this nonvascular host. This highlights the importance of studying endophyte effects across the plant tree of life.

The Sphagnome Project: enabling ecological and evolutionary insights through a genus-level sequencing project.

Considerable progress has been made in ecological and evolutionary genetics with studies demonstrating how genes underlying plant and microbial traits can influence adaptation and even 'extend' to influence community structure and ecosystem level processes. Progress in this area is limited to model systems with deep genetic and genomic resources that often have negligible ecological impact or interest. Thus, important linkages between genetic adaptations and their consequences at organismal and ecological scales are often lacking. Here we introduce the Sphagnome Project, which incorporates genomics into a long-running history of Sphagnum research that has documented unparalleled contributions to peatland ecology, carbon sequestration, biogeochemistry, microbiome research, niche construction, and ecosystem engineering. The Sphagnome Project encompasses a genus-level sequencing effort that represents a new type of model system driven not only by genetic tractability, but by ecologically relevant questions and hypotheses.

Evolutionary origin of the latitudinal diversity gradient in liverworts.

A latitudinal diversity gradient towards the tropics appears as one most recurrent patterns in ecology, but the mechanisms underlying this pattern remain an area of controversy. In angiosperms, the tropical conservatism hypothesis proposes that most groups originated in the tropics and are adapted to a tropical climatic regime, and that relatively few species have evolved physiological adaptations to cold, dry or unpredictable climates. This mechanism is, however, unlikely to apply across land plants, and in particular, to liverworts, a group of about 7500 species, whose ability to withstand cold much better than their tracheophyte counterparts is at odds with the tropical conservatism hypothesis. Molecular dating, diversification rate analyses and ancestral area reconstructions were employed to explore the evolutionary mechanisms that account for the latitudinal diversity gradient in liverworts. As opposed to angiosperms, tropical liverwort genera are not older than their extra-tropical counterparts (median stem age of tropical and extra-tropical liverwort genera of 24.35 ±â€¯39.65 Ma and 39.57 ±â€¯49.07 Ma, respectively), weakening the 'time for speciation hypothesis'. Models of ancestral area reconstructions with equal migration rates between tropical and extra-tropical regions outperformed models with asymmetrical migration rates in either direction. The symmetry and intensity of migrations between tropical and extra-tropical regions suggested by the lack of resolution in ancestral area reconstructions towards the deepest nodes are at odds with the tropical niche conservatism hypothesis. In turn, tropical genera exhibited significantly higher net diversification rates than extra-tropical ones, suggesting that the observed latitudinal diversity gradient results from either higher extinction rates in extra-tropical lineages or higher speciation rates in the tropics. We discuss a series of experiments to help deciphering the underlying evolutionary mechanisms.

Species delimitation and biogeography of a southern hemisphere liverwort clade, Frullania subgenus Microfrullania (Frullaniaceae, Marchantiophyta).

Frullania subgenus Microfrullania is a clade of ca. 15 liverwort species occurring in Australasia, Malesia, and southern South America. We used combined nuclear and chloroplast sequence data from 265 ingroup accessions to test species circumscriptions and estimate the biogeographic history of the subgenus. With dense infra-specific sampling, we document an important role of long-distance dispersal in establishing phylogeographic patterns of extant species. At deeper time scales, a combination of phylogenetic analyses, divergence time estimation and ancestral range estimation were used to reject vicariance and to document the role of long-distance dispersal in explaining the evolution and biogeography of the clade across the southern Hemisphere. A backbone phylogeny for the subgenus is proposed, providing insight into evolution of morphological patterns and establishing the basis for an improved sectional classification of species within Microfrullania. Several species complexes are identified, the presence of two undescribed but genetically and morphologically distinct species is noted, and previously neglected names are discussed.

Range size heritability and diversification patterns in the liverwort genus Radula.

Why some species exhibit larger geographical ranges than others, and to what extent does variation in range size affect diversification rates, remains a fundamental, but largely unanswered question in ecology and evolution. Here, we implement phylogenetic comparative analyses and ancestral area estimations in Radula, a liverwort genus of Cretaceous origin, to investigate the mechanisms that explain differences in geographical range size and diversification rates among lineages. Range size was phylogenetically constrained in the two sub-genera characterized by their almost complete Australasian and Neotropical endemicity, respectively. The congruence between the divergence time of these lineages and continental split suggests that plate tectonics could have played a major role in their present distribution, suggesting that a strong imprint of vicariance can still be found in extant distribution patterns in these highly mobile organisms. Amentuloradula, Volutoradula and Metaradula species did not appear to exhibit losses of dispersal capacities in terms of dispersal life-history traits, but evidence for significant phylogenetic signal in macroecological niche traits suggests that niche conservatism accounts for their restricted geographic ranges. Despite their greatly restricted distribution to Australasia and Neotropics respectively, Amentuloradula and Volutoradula did not exhibit significantly lower diversification rates than more widespread lineages, in contrast with the hypothesis that the probability of speciation increases with range size by promoting geographic isolation and increasing the rate at which novel habitats are encountered. We suggest that stochastic long-distance dispersal events may balance allele frequencies across large spatial scales, leading to low genetic structure among geographically distant areas or even continents, ultimately decreasing the diversification rates in highly mobile, widespread lineages.

Divergent evolution and niche differentiation within the common peatmoss Sphagnum magellanicum.

PREMISE OF THE STUDY: Populations with phenotypic polymorphism in discrete characters may be good models for investigating genome evolution and speciation. Sphagnum magellanicum Brid. is found throughout the northern hemisphere, and despite considerable variation in morphological characters, it is considered one of the least taxonomically controversial peatmoss species. We have observed two main morphs of the species associated with different microhabitats. Here we investigated the genomic and environmental basis of this intraspecific morphological variation. METHODS: We conducted transplant and common garden experiments to test whether the two morphs are genetically differentiated. We then used RAD-sequencing to quantify the genomic divergence between the morphs and approximate Bayesian computation (ABC) to infer the most likely demographic scenario explaining the genome-wide differentiation of the two morphs. KEY RESULTS: We found that genomic differentiation between the two morphs is unexpectedly high and that several of the differentiated morphological characters have a genetic basis. Using simulation approaches, we found support for a scenario of ancient divergence followed by recent secondary contact. CONCLUSIONS: We show that the two morphs represent the two main genetic clusters previously found worldwide. Our results demonstrate that relatively minor morphological differentiation in a presumed phenotypically plastic peatmoss may be associated with massive divergence across the genome.

Phylotranscriptomic analysis of the origin and early diversification of land plants.

Reconstructing the origin and evolution of land plants and their algal relatives is a fundamental problem in plant phylogenetics, and is essential for understanding how critical adaptations arose, including the embryo, vascular tissue, seeds, and flowers. Despite advances in molecular systematics, some hypotheses of relationships remain weakly resolved. Inferring deep phylogenies with bouts of rapid diversification can be problematic; however, genome-scale data should significantly increase the number of informative characters for analyses. Recent phylogenomic reconstructions focused on the major divergences of plants have resulted in promising but inconsistent results. One limitation is sparse taxon sampling, likely resulting from the difficulty and cost of data generation. To address this limitation, transcriptome data for 92 streptophyte taxa were generated and analyzed along with 11 published plant genome sequences. Phylogenetic reconstructions were conducted using up to 852 nuclear genes and 1,701,170 aligned sites. Sixty-nine analyses were performed to test the robustness of phylogenetic inferences to permutations of the data matrix or to phylogenetic method, including supermatrix, supertree, and coalescent-based approaches, maximum-likelihood and Bayesian methods, partitioned and unpartitioned analyses, and amino acid versus DNA alignments. Among other results, we find robust support for a sister-group relationship between land plants and one group of streptophyte green algae, the Zygnematophyceae. Strong and robust support for a clade comprising liverworts and mosses is inconsistent with a widely accepted view of early land plant evolution, and suggests that phylogenetic hypotheses used to understand the evolution of fundamental plant traits should be reevaluated.

Horizontal transfer of an adaptive chimeric photoreceptor from bryophytes to ferns.

Ferns are well known for their shade-dwelling habits. Their ability to thrive under low-light conditions has been linked to the evolution of a novel chimeric photoreceptor--neochrome--that fuses red-sensing phytochrome and blue-sensing phototropin modules into a single gene, thereby optimizing phototropic responses. Despite being implicated in facilitating the diversification of modern ferns, the origin of neochrome has remained a mystery. We present evidence for neochrome in hornworts (a bryophyte lineage) and demonstrate that ferns acquired neochrome from hornworts via horizontal gene transfer (HGT). Fern neochromes are nested within hornwort neochromes in our large-scale phylogenetic reconstructions of phototropin and phytochrome gene families. Divergence date estimates further support the HGT hypothesis, with fern and hornwort neochromes diverging 179 Mya, long after the split between the two plant lineages (at least 400 Mya). By analyzing the draft genome of the hornwort Anthoceros punctatus, we also discovered a previously unidentified phototropin gene that likely represents the ancestral lineage of the neochrome phototropin module. Thus, a neochrome originating in hornworts was transferred horizontally to ferns, where it may have played a significant role in the diversification of modern ferns.

Analyses of transcriptome sequences reveal multiple ancient large-scale duplication events in the ancestor of Sphagnopsida (Bryophyta).

The goal of this research was to investigate whether there has been a whole-genome duplication (WGD) in the ancestry of Sphagnum (peatmoss) or the class Sphagnopsida, and to determine if the timing of any such duplication(s) and patterns of paralog retention could help explain the rapid radiation and current ecological dominance of peatmosses. RNA sequencing (RNA-seq) data were generated for nine taxa in Sphagnopsida (Bryophyta). Analyses of frequency plots for synonymous substitutions per synonymous site (Ks ) between paralogous gene pairs and reconciliation of 578 gene trees were conducted to assess evidence of large-scale or genome-wide duplication events in each transcriptome. Both Ks frequency plots and gene tree-based analyses indicate multiple duplication events in the history of the Sphagnopsida. The most recent WGD event predates divergence of Sphagnum from the two other genera of Sphagnopsida. Duplicate retention is highly variable across species, which might be best explained by local adaptation. Our analyses indicate that the last WGD could have been an important factor underlying the diversification of peatmosses and facilitated their rise to ecological dominance in peatlands. The timing of the duplication events and their significance in the evolutionary history of peat mosses are discussed.