Epiphytic leafy liverworts diversified in angiosperm-dominated forests

The global pattern of epiphytic liverwort disparity: insights from Frullania

The geographical and ecological patterns of morphological disparity are crucial to understand how species are assembled within communities in the context of the evolutionary history, morphological evolution and ecological interactions. However, with limited exceptions, rather few studies have been conducted on the global pattern of disparity, particularly in early land plants. Here we explored the spatial accumulation of disparity in a morphologically variable and species rich liverwort genus Frullania in order to test the hypothesis of latitude disparity gradient. We compiled a morphological data set consisting of eight continuous traits for 244 currently accepted species, and scored the species distribution into 19 floristic regions worldwide. By reconstructing the morphospace of all defined regions and comparisons, we identified a general Gondwana-Laurasia pattern of disparity in Frullania. This likely results from an increase of ecological opportunities and / or relaxed constraints towards low latitudes. The lowest disparity occurred in arid tropical regions, largely due to a high extinction rate as a consequence of paleoaridification. There was weak correlation between species diversity and disparity at different spatial scales. Furthermore, long-distance dispersal may have partially shaped the present-day distribution of Frullania disparity, given its frequency and the great contribution of widely distributed species to local morphospace. This study not only highlighted the crucial roles of paleoenvironmental changes, ecological opportunities, and efficient dispersal on the global pattern of plant disparity, but also implied its dependence on the ecological and physiological function of traits.

The Neotropical endemic liverwort subfamily Micropterygioideae had circum-Antarctic links to the rest of the Lepidoziaceae during the early Cretaceous

Lepidoziaceae are the third-largest family of liverworts, with about 860 species distributed on all continents. The evolutionary history of this family has not been satisfactorily resolved, with taxa such as Micropterygioideae yet to be included in phylogenetic analyses. We inferred a dated phylogeny of Lepidoziaceae using a data set consisting of 13 genetic markers, sampled from 147 species. Based on our phylogenetic estimate, we used statistical dispersal-vicariance analysis to reconstruct the biogeographic history of the family. We inferred a crown age of 197 Ma (95% credible interval 157-240 Ma) for the family in the Australian region, with most major lineages also originating in the same region. Micropterygioideae are placed as the sister group to Lembidioideae, with these two lineages diverging from each other about 132 Ma in the South American-Australian region. With South America and Australia being connected through Antarctica at the time, our results suggest a circum-Antarctic link between Micropterygioideae and the rest of the family. Crown Micropterygioideae were inferred to have arisen 45 Ma in South America before the continent separated from Antarctica. Extinction from southern temperate regions might explain the present-day restriction of Micropterygioideae to the Neotropics. Our study reveals the influence of past geological events, such as continental drift, on the evolution and distribution of a widespread and diverse family of liverworts.

Dating the evolution of the complex thalloid liverworts (Marchantiopsida): total-evidence dating analysis supports a Late Silurian-Early Devonian origin and post-Mesozoic morphological stasis

Divergence times based on molecular clock analyses often differ from those derived from total-evidence dating (TED) approaches. For bryophytes, fossils have been excluded from previous assessments of divergence times, and thus, their utility in dating analyses remains unexplored. Here, we conduct the first TED analyses of the complex thalloid liverworts (Marchantiopsida) that include fossils and evaluate macroevolutionary trends in morphological 'diversity' (disparity) and rates. Phylogenetic analyses were performed on a combined dataset of 130 discrete characters and 11 molecular markers (sampled from nuclear, plastid and mitochondrial genomes). Taxon sampling spanned 56 extant species - representing all the orders within Marchantiophyta and extant genera within Marchantiales - and eight fossil taxa. Total-evidence dating analyses support the radiation of Marchantiopsida during Late Silurian-Early Devonian (or Middle Ordovician when the outgroup is excluded) and that of Ricciaceae in the Middle Jurassic. Morphological change rate was high early in the history of the group, but it barely increased after Late Cretaceous. Disparity-through-time analyses support a fast increase in diversity until the Middle Triassic (c. 250 Ma), after which phenotypic evolution slows down considerably. Incorporating fossils in analyses challenges previous assumptions on the affinities of extinct taxa and indicates that complex thalloid liverworts radiated c. 125 Ma earlier than previously inferred.

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.

Recent origin and diversification accompanied by repeated host shifts of thallus-mining flies (Diptera: Agromyzidae) on liverworts and hornworts

Despite the vast diversity of phytophagous insects that feed on vascular plants (tracheophytes), insects that feed on bryophytes remain understudied. Agromyzidae, one of the most species-rich phytophagous clades in Diptera, consists mainly of leaf-mining species that feed on tracheophytes. However, a recent discovery of thallus-mining species on liverworts and hornworts within the Liriomyza group of Phytomyzinae provides an opportunity to study host shifts between tracheophytes and bryophytes. This study aimed to explore the origin and diversification of thallus-miners and estimate the pattern and timing of host shifts. Phylogenetic analysis of Phytomyzinae has revealed that the thallus-mining agromyzids formed a separate clade, which was sister to a fern pinnule-miner. The diversification of bryophyte-associated agromyzids since the Oligocene involved multiple host shifts across various bryophyte taxa. The diversification of the thallus-mining Phytoliriomyza may have occurred at the same time as the leaf-mining agromyzid flies on herbaceous plants, indicating a dynamic history of interactions between bryophytes and herbivores in angiosperms-dominated ecosystems.

A complete genus-level phylogeny reveals the Cretaceous biogeographic diversification of the poppy family

Angiosperms, a trigger for the Cretaceous Terrestrial Revolution (KTR), underwent a rapid expansion and occupied all the environments during the Mid-Upper Cretaceous. Yet, Cretaceous biogeographic patterns and processes underlying the distribution of angiosperm diversity in the Northern Hemisphere are still poorly known. Here, we elucidated the biogeographic diversification of the angiosperm family Papaveraceae, an ancient Northern Hemisphere clade characterized by poor dispersal ability and high level of regional endemism. Based on both plastome and multi-locus datasets, we reconstructed a robust time-calibrated phylogeny that includes all currently recognized 45 genera of this family. Within the time-calibrated phylogenetic framework, we conducted 72 biogeographic analyses by testing the sensitivity of uncertainties of area delimitation, maxarea constraints, and the parameters of the model, i.e., j (describing jump-dispersal events) and w (modifying dispersal multiplier matrices), to ancestral range estimations. We also inferred ancestral habitat and ecological niches. Phylogenetic analyses strongly support Papaveraceae as monophyletic. Pteridophylloideae is strongly supported as sister to Hypecoideae-Fumarioideae. Our results indicate that the j parameter and number of predefined areas strongly affect ancestral range estimates, generating questionable ancestral ranges, whereas maxarea constraint and w parameter have no effect and improve model fit. After accounting for these uncertainties, our results indicate that Papaveraceae differentiated in Asian wet forests during the Lower Cretaceous and subsequently occupied the Asian and western North American arid and open areas. Three dispersals from Asia to western North America via the Bering land bridge occurred in the Mid-Upper Cretaceous, largely in agreement with the KTR. Habitat shift and ecological niche divergence resulted in the subsequent disjunctions between Asia and western North America. These findings suggest that the interplay of range expansion and niche divergence-driven vicariance might have shaped Cretaceous biogeographic patterns of angiosperms with Papaveraceae-like ecological requirements and dispersal abilities in the Northern Hemisphere, hence contributing to the knowledge on the geographic expansion of angiosperms during the KTR.

Filmy Ferns (Hymenophyllaceae) and Associated Spike-Mosses (Selaginellaceae) from the Mid-Cretaceous Kachin Amber, Myanmar

Filmy ferns (Hymenophyllaceae) are the most diverse lineage of the early-diverging leptosporangiate ferns with ca. 430 species widely distributed around the world but with the highest diversity in the humid tropics. However, their fossil record is scarce because of the low preservation potential of the delicate, membranous laminae. So far, no Hymenophyllaceae fossils have been reported from tropical Asia. Here, we describe some fern remains and their syninclusions (spike-mosses) in four pieces of Kachin amber from the mid-Cretaceous of Hukawng Valley, Northern Myanmar, as Hymenophyllites angustus sp. nov., H. kachinensis sp. nov., H. setosus sp. nov. (Hymenophyllaceae) and Selaginella alata sp. nov. (Selaginellaceae), respectively. These fern remains are assigned to Hymenophyllaceae based on the filmy, one-cell thick, decompound pinnatifid laminae and dichotomous venation. They represent the first fossil record of Hymenophyllaceae in tropical Asia. The growth habits of these ferns and associated spike-mosses and their implication for paleoenvironment are discussed. Our study expands the diversity of the cryptogams in mid-Cretaceous Kachin amber. Together with other contemporaneous findings, the present fossils indicate that Hymenophyllaceae have already accumulated some notable diversity in the Cretaceous.

The renaissance and enlightenment of Marchantia as a model system

The liverwort Marchantia polymorpha has been utilized as a model for biological studies since the 18th century. In the past few decades, there has been a Renaissance in its utilization in genomic and genetic approaches to investigating physiological, developmental, and evolutionary aspects of land plant biology. The reasons for its adoption are similar to those of other genetic models, e.g. simple cultivation, ready access via its worldwide distribution, ease of crossing, facile genetics, and more recently, efficient transformation, genome editing, and genomic resources. The haploid gametophyte dominant life cycle of M. polymorpha is conducive to forward genetic approaches. The lack of ancient whole-genome duplications within liverworts facilitates reverse genetic approaches, and possibly related to this genomic stability, liverworts possess sex chromosomes that evolved in the ancestral liverwort. As a representative of one of the three bryophyte lineages, its phylogenetic position allows comparative approaches to provide insights into ancestral land plants. Given the karyotype and genome stability within liverworts, the resources developed for M. polymorpha have facilitated the development of related species as models for biological processes lacking in M. polymorpha.

The role of paleontological data in bryophyte systematics

Systematics reconstructs tempo and mode in biological evolution by resolving the phylogenetic fabric of biodiversity. The staggering duration and complexity of evolution, coupled with loss of information (extinction), render exhaustive reconstruction of the evolutionary history of life unattainable. Instead, we sample its products-phenotypes and genotypes-to generate phylogenetic hypotheses, which we sequentially reassess and update against new data. Current consensus in evolutionary biology emphasizes fossil integration in total-evidence analyses, requiring in-depth understanding of fossils-age, phenotypes, and systematic affinities-and a detailed morphological framework uniting fossil and extant taxa. Bryophytes present a special case: deep evolutionary history but sparse fossil record and phenotypic diversity encompassing small dimensional scales. We review how these peculiarities shape fossil inclusion in bryophyte systematics. Paucity of the bryophyte fossil record, driven primarily by phenotypic (small plant size) and ecological constraints (patchy substrate-hugging populations), and incomplete exploration, results in many morphologically isolated, taxonomically ambiguous fossil taxa. Nevertheless, instances of exquisite preservation and pioneering studies demonstrate the feasibility of including bryophyte fossils in evolutionary inference. Further progress will arise from developing extensive morphological matrices for bryophytes, continued exploration of the fossil record, re-evaluation of previously described fossils, and training specialists in identification and characterization of bryophyte fossils, and in bryophyte morphology.

The Angiosperm Terrestrial Revolution and the origins of modern biodiversity

Biodiversity today has the unusual property that 85% of plant and animal species live on land rather than in the sea, and half of these live in tropical rainforests. An explosive boost to terrestrial diversity occurred from c. 100-50 million years ago, the Late Cretaceous and early Palaeogene. During this interval, the Earth-life system on land was reset, and the biosphere expanded to a new level of productivity, enhancing the capacity and species diversity of terrestrial environments. This boost in terrestrial biodiversity coincided with innovations in flowering plant biology and evolutionary ecology, including their flowers and efficiencies in reproduction; coevolution with animals, especially pollinators and herbivores; photosynthetic capacities; adaptability; and ability to modify habitats. The rise of angiosperms triggered a macroecological revolution on land and drove modern biodiversity in a secular, prolonged shift to new, high levels, a series of processes we name here the Angiosperm Terrestrial Revolution.

Deep time diversity and the early radiations of birds

Reconstructing the history of biodiversity has been hindered by often-separate analyses of stem and crown groups of the clades in question that are not easily understood within the same unified evolutionary framework. Here, we investigate the evolutionary history of birds by analyzing three supertrees that combine published phylogenies of both stem and crown birds. Our analyses reveal three distinct large-scale increases in the diversification rate across bird evolutionary history. The first increase, which began between 160 and 170 Ma and reached its peak between 130 and 135 Ma, corresponds to an accelerated morphological evolutionary rate associated with the locomotory systems among early stem birds. This radiation resulted in morphospace occupation that is larger and different from their close dinosaurian relatives, demonstrating the occurrence of a radiation among early stem birds. The second increase, which started ∼90 Ma and reached its peak between 65 and 55 Ma, is associated with rapid evolution of the cranial skeleton among early crown birds, driven differently from the first radiation. The third increase, which occurred after ∼40 to 45 Ma, has yet to be supported by quantitative morphological data but gains some support from the fossil record. Our analyses indicate that the bird biodiversity evolution was influenced mainly by long-term climatic changes and also by major paleobiological events such as the Cretaceous-Paleogene (K-Pg) extinction.

Insights Into the Evolutionary History of the Subfamily Orthotrichoideae (Orthotrichaceae, Bryophyta): New and Former Supra-Specific Taxa So Far Obscured by Prevailing Homoplasy

Mosses of the subfamily Orthotrichoideae represent one of the main components of the cryptogam epiphytic communities in temperate areas. During the last two decades, this taxonomical group has undergone an extensive revision that has led to its rearrangement at the generic level. However, their phylogenetic relationships and inferences on the evolutionary patterns that have driven the present diversity have little advanced. In this study, we present a dated molecular phylogenetic reconstruction at the subfamily level, including 130 samples that represent the 12 genera currently recognized within the subfamily, and the analysis of four molecular markers: ITS2, rps4, trnG, and trnL-F. We also analyze 13 morphological characters of systematic value to infer their origin and diagnostic utility within the subfamily. The phylogenetic reconstruction yields three main clades within the subfamily, two of which correspond to the tribe Zygodonteae, and one to Orthotricheae. Within Zygodonteae, the genus Zygodon results to be a polyphyletic artificial assembly, and we propose to separate a new genus named Australoria. Conversely, our results do not support the separation of Pentastichella and Pleurorthotrichum at the genus level and we therefore propose to include Pleurorthotrichum in Pentastichella. Regarding Orthotricheae, our analyses clearly allow the distinction of two subtribes: Orthotrichinae and Lewinskyinae. Within the latter, Ulota results a polyphyletic entity, and therefore we propose the segregation of a separate new genus named Atlantichella. Dating analyses allow us to conclude that the split of the tribes within Orthotrichoideae dates from the Middle Jurassic, while the diversification of Orthotrichum and Zygodon probably started during the Late Cretaceous. However, most of the extant genera of this subfamily seem to be younger, and apparently its highest diversification burst took place during the Oligocene. Finally, the analysis of the morphological traits reveals that most of the characters previously used to separate genera and here tested are homoplastic, which has hindered the taxonomical and systematic proposals for decades. However, even if there are no exclusive characters, all of the genera can be defined by the combination of a few characters.

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

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

Life in the canopy: community trait assessments reveal substantial functional diversity among fern epiphytes

The expansion of angiosperm-dominated forests in the Cretaceous and early Cenozoic had a profound effect on terrestrial biota by creating novel ecological niches. The majority of modern fern lineages are hypothesized to have arisen in response to this expansion, particularly fern epiphytes that radiated into the canopy. Recent evidence, however, suggests that epiphytism does not correlate with increased diversification rates in ferns, calling into question the role of the canopy habitat in fern evolution. To understand the role of the canopy in structuring fern community diversity, we investigated functional traits of fern sporophytes and gametophytes across a broad phylogenetic sampling on the island of Moorea, French Polynesia, including > 120 species and representatives of multiple epiphytic radiations. While epiphytes showed convergence in small size and a higher frequency of noncordate gametophytes, they showed greater functional diversity at the community level relative to terrestrial ferns. These results suggest previously overlooked functional diversity among fern epiphytes, and raise the hypothesis that while the angiosperm canopy acted as a complex filter that restricted plant size, it also facilitated diversification into finely partitioned niches. Characterizing these niche axes and adaptations of epiphytic ferns occupying them should be a priority for future pteridological research.

Recent accelerated diversification in rosids occurred outside the tropics

Conflicting relationships have been found between diversification rate and temperature across disparate clades of life. Here, we use a supermatrix comprising nearly 20,000 species of rosids-a clade of ~25% of all angiosperm species-to understand global patterns of diversification and its climatic association. Our approach incorporates historical global temperature, assessment of species' temperature niche, and two broad-scale characterizations of tropical versus non-tropical niche occupancy. We find the diversification rates of most subclades dramatically increased over the last 15 million years (Myr) during cooling associated with global expansion of temperate habitats. Climatic niche is negatively associated with diversification rates, with tropical rosids forming older communities and experiencing speciation rates ~2-fold below rosids in cooler climates. Our results suggest long-term cooling had a disproportionate effect on non-tropical diversification rates, leading to dynamic young communities outside of the tropics, while relative stability in tropical climes led to older, slower-evolving but still species-rich communities.

Estimating rates and patterns of diversification with incomplete sampling: a case study in the rosids

PREMISE

Recent advances in generating large-scale phylogenies enable broad-scale estimation of species diversification. These now common approaches typically are characterized by (1) incomplete species coverage without explicit sampling methodologies and/or (2) sparse backbone representation, and usually rely on presumed phylogenetic placements to account for species without molecular data. We used empirical examples to examine the effects of incomplete sampling on diversification estimation and provide constructive suggestions to ecologists and evolutionary biologists based on those results.

METHODS

We used a supermatrix for rosids and one well-sampled subclade (Cucurbitaceae) as empirical case studies. We compared results using these large phylogenies with those based on a previously inferred, smaller supermatrix and on a synthetic tree resource with complete taxonomic coverage. Finally, we simulated random and representative taxon sampling and explored the impact of sampling on three commonly used methods, both parametric (RPANDA and BAMM) and semiparametric (DR).

RESULTS

We found that the impact of sampling on diversification estimates was idiosyncratic and often strong. Compared to full empirical sampling, representative and random sampling schemes either depressed or inflated speciation rates, depending on methods and sampling schemes. No method was entirely robust to poor sampling, but BAMM was least sensitive to moderate levels of missing taxa.

CONCLUSIONS

We suggest caution against uncritical modeling of missing taxa using taxonomic data for poorly sampled trees and in the use of summary backbone trees and other data sets with high representative bias, and we stress the importance of explicit sampling methodologies in macroevolutionary studies.

Fossil fern rhizomes as a model system for exploring epiphyte community structure across geologic time: evidence from Patagonia

Background

In extant ecosystems, complex networks of ecological interactions between organisms can be readily studied. In contrast, understanding of such interactions in ecosystems of the geologic past is incomplete. Specifically, in past terrestrial ecosystems we know comparatively little about plant biotic interactions besides saprotrophy, herbivory, mycorrhizal associations, and oviposition. Due to taphonomic biases, epiphyte communities are particularly rare in the plant-fossil record, despite their prominence in modern ecosystems. Accordingly, little is known about how terrestrial epiphyte communities have changed across geologic time. Here, we describe a tiny in situ fossil epiphyte community that sheds light on plant-animal and plant-plant interactions more than 50 million years ago.

Methods

A single silicified Todea (Osmundaceae) rhizome from a new locality of the early Eocene (ca. 52 Ma) Tufolitas Laguna del Hunco (Patagonia, Argentina) was studied in serial thin sections using light microscopy. The community of organisms colonizing the tissues of the rhizome was characterized by identifying the organisms and mapping and quantifying their distribution. A 200 × 200 µm grid was superimposed onto the rhizome cross section, and the colonizers present at each node of the grid were tallied.

Results

Preserved in situ, this community offers a rare window onto aspects of ancient ecosystems usually lost to time and taphonomic processes. The community is surprisingly diverse and includes the first fossilized leafy liverworts in South America, also marking the only fossil record of leafy bryophyte epiphytes outside of amber deposits; as well as several types of fungal hyphae and spores; microsclerotia with possible affinities in several ascomycete families; and evidence for oribatid mites.

Discussion

The community associated with the Patagonian rhizome enriches our understanding of terrestrial epiphyte communities in the distant past and adds to a growing body of literature on osmundaceous rhizomes as important hosts for component communities in ancient ecosystems, just as they are today. Because osmundaceous rhizomes represent an ecological niche that has remained virtually unchanged over time and space and are abundant in the fossil record, they provide a paleoecological model system that could be used to explore epiphyte community structure through time.

Accelerated diversifications in three diverse families of morphologically complex lichen-forming fungi link to major historical events

Historical mass extinction events had major impacts on biodiversity patterns. The most recent and intensively studied event is the Cretaceous - Paleogene (K-Pg) boundary (ca. 66 million years ago [MYA]). However, the factors that may have impacted diversification dynamics vary across lineages. We investigated the macroevolutionary dynamics with a specific focus on the impact of major historical events such as the K-Pg mass extinction event on two major subclasses - Lecanoromycetidae and Ostropomycetidae - of lichen-forming fungi and tested whether variation in the rate of diversification can be associated with the evolution of a specific trait state - macrolichen. Our results reveal accelerated diversification events in three families of morphologically complex lichen-forming fungi - Cladoniaceae, Parmeliaceae, and Peltigeraceae - which are from the subclass Lecanoromycetidae and mostly composed of macrolichens, those that form three dimensional structures. Our RTT plot result for the subclass Lecanoromycetidae also reveals accelerated diversification. Changes in diversification rates occurred around the transition between Mesozoic and Cenozoic eras and was likely related to the K-Pg mass extinction event. The phylogenetic positions for rate increases estimated based on marginal shift probability are, however, scattered from 100 to 40 MYA preventing us from making explicit inference. Although we reveal that the phenotypic state of macrolichens is associated with a higher diversification rate than microlichens, we also show that the evolution of macrolichens predated the K-Pg event. Furthermore, the association between macrolichens and increased diversification is not universal and can be explained, in part, by phylogenetic relatedness. By investigating the macroevolutionary dynamics of lichen-forming fungi our study provides a new empirical system suitable to test the effect of major historical event on shaping biodiversity patterns and to investigate why changes in biodiversity patterns are not in concordance across clades. Our results imply that multiple historical events during the transition from Mesozoic to Cenozoic eras, including the K-Pg mass extinction event, impacted the evolutionary dynamics in lichen-forming fungi. However, future studies focusing on individual lichen-forming fungal families are required to ascertain whether diversification rates are associated with growth form and certain geological events.

Multiple historical processes obscure phylogenetic relationships in a taxonomically difficult group (Lobariaceae, Ascomycota)

In the age of next-generation sequencing, the number of loci available for phylogenetic analyses has increased by orders of magnitude. But despite this dramatic increase in the amount of data, some phylogenomic studies have revealed rampant gene-tree discordance that can be caused by many historical processes, such as rapid diversification, gene duplication, or reticulate evolution. We used a target enrichment approach to sample 400 single-copy nuclear genes and estimate the phylogenetic relationships of 13 genera in the lichen-forming family Lobariaceae to address the effect of data type (nucleotides and amino acids) and phylogenetic reconstruction method (concatenation and species tree approaches). Furthermore, we examined datasets for evidence of historical processes, such as rapid diversification and reticulate evolution. We found incongruence associated with sequence data types (nucleotide vs. amino acid sequences) and with different methods of phylogenetic reconstruction (species tree vs. concatenation). The resulting phylogenetic trees provided evidence for rapid and reticulate evolution based on extremely short branches in the backbone of the phylogenies. The observed rapid and reticulate diversifications may explain conflicts among gene trees and the challenges to resolving evolutionary relationships. Based on divergence times, the diversification at the backbone occurred near the Cretaceous-Paleogene (K-Pg) boundary (65 Mya) which is consistent with other rapid diversifications in the tree of life. Although some phylogenetic relationships within the Lobariaceae family remain with low support, even with our powerful phylogenomic dataset of up to 376 genes, our use of target-capturing data allowed for the novel exploration of the mechanisms underlying phylogenetic and systematic incongruence.

MicroRNAs in Marchantia polymorpha

Contents Summary 409 I. Introduction 409 II. RNA silencing machinery in Marchantia polymorpha 410 III. miRNA prediction by integrating omics approach 410 IV. miRNAs and their targets in Marchantia polymorpha 410 V. Mpo-miR390-mediated MpTAS3 tasiRNA biogenesis and potential tasiARF target MpARF2 414 VI. Artificial miRNA and CRISPR-CAS9 edited MIR gene in Marchantia polymorpha 414 VII. Conclusions 415 Acknowledgements 415 References 415 SUMMARY: The liverwort Marchantia polymorpha occupies an important phylogenetic position for comparative studies of land plant gene regulation. Multiple gene regulatory pathways mediated by small RNAs, including microRNAs (miRNAs), trans-acting short-interfering RNAs, and heterochromatic siRNAs often associated with RNA-dependent DNA methylation, have been characterized in flowering plants. Genes for essential components for all of these small RNA-mediated gene silencing pathways are found in M. polymorpha as well as the moss Phsycomitrella patens, indicating that these pathways existed in the ancestral land plant. However, only seven miRNAs are conserved across land plants, with both ancestral and novel targets identified in M. polymorpha. There is little or no evidence that any of these conserved miRNAs are present in algae. As with other plants investigated, most miRNAs in M. polypmorpha exhibit lineage-specific evolution. Application of artificial miRNA and CRISPR-Cas9 technologies in genetic studies of M. polymorpha provide avenues to further investigate miRNA biology.

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.

Constraining uncertainty in the timescale of angiosperm evolution and the veracity of a Cretaceous Terrestrial Revolution

Through the lens of the fossil record, angiosperm diversification precipitated a Cretaceous Terrestrial Revolution (KTR) in which pollinators, herbivores and predators underwent explosive co-diversification. Molecular dating studies imply that early angiosperm evolution is not documented in the fossil record. This mismatch remains controversial. We used a Bayesian molecular dating method to analyse a dataset of 83 genes from 644 taxa and 52 fossil calibrations to explore the effect of different interpretations of the fossil record, molecular clock models, data partitioning, among other factors, on angiosperm divergence time estimation. Controlling for different sources of uncertainty indicates that the timescale of angiosperm diversification is much less certain than previous molecular dating studies have suggested. Discord between molecular clock and purely fossil-based interpretations of angiosperm diversification may be a consequence of false precision on both sides. We reject a post-Jurassic origin of angiosperms, supporting the notion of a cryptic early history of angiosperms, but this history may be as much as 121 Myr, or as little as 23 Myr. These conclusions remain compatible with palaeobotanical evidence and a more general KTR in which major groups of angiosperms diverged later within the Cretaceous, alongside the diversification of pollinators, herbivores and their predators.

Challenges of comprehensive taxon sampling in comparative biology: Wrestling with rosids

Using phylogenetic approaches to test hypotheses on a large scale, in terms of both species sampling and associated species traits and occurrence data-and doing this with rigor despite all the attendant challenges-is critical for addressing many broad questions in evolution and ecology. However, application of such approaches to empirical systems is hampered by a lingering series of theoretical and practical bottlenecks. The community is still wrestling with the challenges of how to develop species-level, comprehensively sampled phylogenies and associated geographic and phenotypic resources that enable global-scale analyses. We illustrate difficulties and opportunities using the rosids as a case study, arguing that assembly of biodiversity data that is scale-appropriate-and therefore comprehensive and global in scope-is required to test global-scale hypotheses. Synthesizing comprehensive biodiversity data sets in clades such as the rosids will be key to understanding the origin and present-day evolutionary and ecological dynamics of the angiosperms.

Novel climates reverse carbon uptake of atmospherically dependent epiphytes: Climatic constraints on the iconic boreal forest lichen Evernia mesomorpha

PREMISE OF THE STUDY

Changing climates are expected to affect the abundance and distribution of global vegetation, especially plants and lichens with an epiphytic lifestyle and direct exposure to atmospheric variation. The study of epiphytes could improve understanding of biological responses to climatic changes, but only if the conditions that elicit physiological performance changes are clearly defined.

METHODS

We evaluated individual growth performance of the epiphytic lichen Evernia mesomorpha, an iconic boreal forest indicator species, in the first year of a decade-long experiment featuring whole-ecosystem warming and drying. Field experimental enclosures were located near the southern edge of the species' range.

KEY RESULTS

Mean annual biomass growth of Evernia significantly declined 6 percentage points for every +1°C of experimental warming after accounting for interactions with atmospheric drying. Mean annual biomass growth was 14% in ambient treatments, 2% in unheated control treatments, and -9% to -19% (decreases) in energy-added treatments ranging from +2.25 to +9.00°C above ambient temperatures. Warming-induced biomass losses among persistent individuals were suggestive evidence of an extinction debt that could precede further local mortality events.

CONCLUSIONS

Changing patterns of warming and drying would decrease or reverse Evernia growth at its southern range margins, with potential consequences for the maintenance of local and regional populations. Negative carbon balances among persisting individuals could physiologically commit these epiphytes to local extinction. Our findings illuminate the processes underlying local extinctions of epiphytes and suggest broader consequences for range shrinkage if dispersal and recruitment rates cannot keep pace.

Rates of ecological divergence and body size evolution are correlated with species diversification in scaly tree ferns

Variation in species richness across regions and between different groups of organisms is a major feature of evolution. Several factors have been proposed to explain these differences, including heterogeneity in the rates of species diversification and the age of clades. It has been frequently assumed that rapid rates of diversification are coupled to high rates of ecological and morphological evolution, leading to a prediction that remains poorly explored for most species: the positive association between ecological niche divergence, morphological evolution and species diversification. We combined a time-calibrated phylogeny with distribution, ecological and body size data for scaly tree ferns (Cyatheaceae) to test whether rates of species diversification are predicted by the rates at which clades have evolved distinct ecological niches and body sizes. We found that rates of species diversification are positively correlated with rates of ecological and morphological evolution, with rapidly diversifying clades also showing rapidly evolving ecological niches and body sizes. Our results show that rapid diversification of scaly tree ferns is associated with the evolution of species with comparable morphologies that diversified into similar, yet distinct, environments. This suggests parallel evolutionary pathways opening in different tropical regions whenever ecological and geographical opportunities arise. Accordingly, rates of ecological niche and body size evolution are relevant to explain the current patterns of species richness in this 'ancient' fern lineage across the tropics.

Species richness and composition of epiphytic bryophytes in flooded forests of Caxiuanã National Forest, Eastern Amazon, Brazil

This study aimed to compare the richness and composition of the epiphytic bryoflora between várzea and igapó forests in Caxiuanã National Forest, Brazilian Amazon. Bryophytes were collected on 502 phorophytes of Virola surinamensis. Average richness per phorophyte and composition between forests and between dry and rainy periods was tested by two-way analysis and by cluster analysis, respectively. In total, 54 species of 13 families were identified. Richness was greater in igapó forest (44 species) compared to várzea forest (38 species). There was no significant difference in the number of species between the studied periods. Cluster analysis showed the bryoflora composition was different between várzea and igapó, but not between dry and rainy periods. Results did not corroborate the hypothesis that várzea forests harbor higher species richness than igapó forests.

Are diversification rates and chromosome evolution in the temperate grasses (Pooideae) associated with major environmental changes in the Oligocene-Miocene?

The Pooideae are a highly diverse C3 grass subfamily that includes some of the most economically important crops, nested within the highly speciose core-pooid clade. Here, we build and explore the phylogeny of the Pooideae within a temporal framework, assessing its patterns of diversification and its chromosomal evolutionary changes in the light of past environmental transformations. We sequenced five plastid DNA loci, two coding (ndhF, matk) and three non-coding (trnH-psbA, trnT-L and trnL-F), in 163 Poaceae taxa, including representatives for all subfamilies of the grasses and all but four ingroup Pooideae tribes. Parsimony and Bayesian phylogenetic analyses were conducted and divergence times were inferred in BEAST using a relaxed molecular clock. Diversification rates were assessed using the MEDUSA approach, and chromosome evolution was analyzed using the chromEvol software. Diversification of the Pooideae started in the Late-Eocene and was especially intense during the Oligocene-Miocene. The background diversification rate increased significantly at the time of the origin of the Poodae + Triticodae clade. This shift in diversification occurred in a context of falling temperatures that potentially increased ecological opportunities for grasses adapted to open areas around the world. The base haploid chromosome number n = 7 has remained stable throughout the phylogenetic history of the core pooids and we found no link between chromosome transitions and major diversification events in the Pooideae.

Geographical structure, narrow species ranges, and Cenozoic diversification in a pantropical clade of epiphyllous leafy liverworts

The evolutionary history and classification of epiphyllous cryptogams are still poorly known. Leptolejeunea is a largely epiphyllous pantropical liverwort genus with about 25 species characterized by deeply bilobed underleaves, elliptic to narrowly obovate leaf lobes, the presence of ocelli, and vegetative reproduction by cladia. Sequences of three chloroplast regions (rbcL, trnL-F, psbA) and the nuclear ribosomal ITS region were obtained for 66 accessions of Leptolejeunea and six outgroup species to explore the phylogeny, divergence times, and ancestral areas of this genus. The phylogeny was estimated using maximum-likelihood and Bayesian inference approaches, and divergence times were estimated with a Bayesian relaxed clock method. Leptolejeunea likely originated in Asia or the Neotropics within a time interval from the Early Eocene to the Late Cretaceous (67.9 Ma, 95% highest posterior density [HPD]: 47.9-93.7). Diversification of the crown group initiated in the Eocene or early Oligocene (38.4 Ma, 95% HPD: 27.2-52.6). Most species clades were established in the Miocene. Leptolejeunea epiphylla and L. schiffneri originated in Asia and colonized African islands during the Plio-Pleistocene. Accessions of supposedly pantropical species are placed in different main clades. Several monophyletic morphospecies exhibit considerable sequence variation related to a geographical pattern. The clear geographic structure of the Leptolejeunea crown group points to evolutionary processes including rare long-distance dispersal and subsequent speciation. Leptolejeunea may have benefitted from the large-scale distribution of humid tropical angiosperm forests in the Eocene.

By how much do we underestimate species diversity of liverworts using morphological evidence? An example from Australasian Plagiochila (Plagiochilaceae: Jungermanniopsida)

As a framework for revisionary study of the leafy liverwort Plagiochila in Australia, two methods for species delimitation on molecular sequence data, General Mixed Yule Coalescence model (GMYC) and Automatic Barcode Gap Discovery (ABGD) were applied to a dataset including 265 individuals from Australia, New Zealand, and the Pacific. Groups returned by GMYC and ABGD were incongruent in some lineages, and ABGD tended to lump groups. This may reflect underlying heterogeneity in the history of diversification within different lineages of Plagiochila. GMYC from trees calculated using three different molecular clocks were compared, in some lineages different primary species hypotheses were returned by analyses of trees estimated under different clock models, suggesting clock model selection should be a routine component of phylogeny reconstruction for tree-based species delimitation methods, such as GMYC. Our results suggest that a minimum of 71 Plagiochilaceae species occur in Australasia, 16 more than currently accepted for the region, comprising 8 undetermined species and 8 synonyms requiring reinstatement. Despite modern taxonomic investigation over a four decade period, (1) real diversity is 29% higher than currently recognized; and (2) 12 of 33, or 36%, of currently accepted and previously untested Australasian species have circumscription issues, including polyphyly, paraphyly, internal phylogenetic structure, or combinations of two or more of these issues. These both reflect the many challenges associated with grouping decisions based solely on morphological data in morphologically simple yet polymorphic plant lineages. Our results highlight again the critical need for combined molecular-morphological datasets as a basis for resolving robust species hypotheses in species-rich bryophyte lineages.

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.

Evolution of Epiphytism and Fruit Traits Act Unevenly on the Diversification of the Species-Rich Genus Peperomia (Piperaceae)

The species-rich genus Peperomia (Black Pepper relatives) is the only genus among early diverging angiosperms where epiphytism evolved. The majority of fruits of Peperomia release sticky secretions or exhibit hook-shaped appendages indicative of epizoochorous dispersal, which is in contrast to other flowering plants, where epiphytes are generally characterized by fruit morphological adaptations for anemochory or endozoochory. We investigate fruit characters using Cryo-SEM. Comparative phylogenetic analyses are applied for the first time to include life form and fruit character information to study diversification in Peperomia. Likelihood ratio tests uncover correlated character evolution. We demonstrate that diversification within Peperomia is not homogenous across its phylogeny, and that net diversification rates increase by twofold within the most species-rich subgenus. In contrast to former land plant studies that provide general evidence for increased diversification in epiphytic lineages, we demonstrate that the evolution of epiphytism within Peperomia predates the diversification shift. An epiphytic-dependent diversification is only observed for the background phylogeny. An elevated frequency of life form transitions between epiphytes and terrestrials and thus evolutionary flexibility of life forms is uncovered to coincide with the diversification shift. The evolution of fruits showing dispersal related structures is key to diversification in the foreground region of the phylogeny and postdates the evolution of epiphytism. We conclude that the success of Peperomia, measured in species numbers, is likely the result of enhanced vertical and horizontal dispersal ability and life form flexibility but not the evolution of epiphytism itself.

The rise of angiosperm-dominated herbaceous floras: Insights from Ranunculaceae

The rise of angiosperms has been regarded as a trigger for the Cretaceous revolution of terrestrial ecosystems. However, the timeframe of the rise angiosperm-dominated herbaceous floras (ADHFs) is lacking. Here, we used the buttercup family (Ranunculaceae) as a proxy to provide insights into the rise of ADHFs. An integration of phylogenetic, molecular dating, ancestral state inferring, and diversification analytical methods was used to infer the early evolutionary history of Ranunculaceae. We found that Ranunculaceae became differentiated in forests between about 108–90 Ma. Diversification rates markedly elevated during the Campanian, mainly resulted from the rapid divergence of the non-forest lineages, but did not change across the Cretaceous-Paleogene boundary. Our data for Ranunculaceae indicate that forest-dwelling ADHFs may have appeared almost simultaneously with angiosperm-dominated forests during the mid-Cretaceous, whereas non-forest ADHFs arose later, by the end of the Cretaceous terrestrial revolution. Furthermore, ADHFs were relatively unaffected by the Cretaceous-Paleogene mass extinction.

Crown Group Lejeuneaceae and Pleurocarpous Mosses in Early Eocene (Ypresian) Indian Amber

Cambay amber originates from the warmest period of the Eocene, which is also well known for the appearance of early angiosperm-dominated megathermal forests. The humid climate of these forests may have triggered the evolution of epiphytic lineages of bryophytes; however, early Eocene fossils of bryophytes are rare. Here, we present evidence for lejeuneoid liverworts and pleurocarpous mosses in Cambay amber. The preserved morphology of the moss fossil is inconclusive for a detailed taxonomic treatment. The liverwort fossil is, however, distinctive; its zig-zagged stems, suberect complicate-bilobed leaves, large leaf lobules, and small, deeply bifid underleaves suggest a member of Lejeuneaceae subtribe Lejeuneinae (Harpalejeunea, Lejeunea, Microlejeunea). We tested alternative classification possibilities by conducting divergence time estimates based on DNA sequence variation of Lejeuneinae using the age of the fossil for corresponding age constraints. Consideration of the fossil as a stem group member of Microlejeunea or Lejeunea resulted in an Eocene to Late Cretaceous age of the Lejeuneinae crown group. This reconstruction is in good accordance with published divergence time estimates generated without the newly presented fossil evidence. Balancing available evidence, we describe the liverwort fossil as the extinct species Microlejeunea nyiahae, representing the oldest crown group fossil of Lejeuneaceae.

Divergence Times and Phylogenetic Patterns of Sebacinales, a Highly Diverse and Widespread Fungal Lineage

Patterns of geographic distribution and composition of fungal communities are still poorly understood. Widespread occurrence in terrestrial ecosystems and the unique richness of interactions of Sebacinales with plants make them a target group to study evolutionary events in the light of nutritional lifestyle. We inferred diversity patterns, phylogenetic structures and divergence times of Sebacinales with respect to their nutritional lifestyles by integrating data from fossil-calibrated phylogenetic analyses. Relaxed molecular clock analyses indicated that Sebacinales originated late Permian within Basidiomycota, and their split into Sebacinaceae and Serendipitaceae nom. prov. likely occurred during the late Jurassic and the early Cretaceous, coinciding with major diversifications of land plants. In Sebacinaceae, diversification of species with ectomycorrhizal lifestyle presumably started during the Paleocene. Lineage radiations of the core group of ericoid and cavendishioid mycorrhizal Sebacinales started probably in the Eocene, coinciding with diversification events of their hosts. The diversification of Sebacinales with jungermannioid interactions started during the Oligocene, and occurred much later than the diversification of their hosts. Sebacinales communities associated either with ectomycorrhizal plants, achlorophyllous orchids, ericoid and cavendishioid Ericaceae or liverworts were phylogenetically clustered and globally distributed. Major Sebacinales lineage diversifications started after the continents had drifted apart. We also briefly discuss dispersal patterns of extant Sebacinales.

Divergence times and the evolution of morphological complexity in an early land plant lineage (Marchantiopsida) with a slow molecular rate

We present a complete generic-level phylogeny of the complex thalloid liverworts, a lineage that includes the model system Marchantia polymorpha. The complex thalloids are remarkable for their slow rate of molecular evolution and for being the only extant plant lineage to differentiate gas exchange tissues in the gametophyte generation. We estimated the divergence times and analyzed the evolutionary trends of morphological traits, including air chambers, rhizoids and specialized reproductive structures. A multilocus dataset was analyzed using maximum likelihood and Bayesian approaches. Relative rates were estimated using local clocks. Our phylogeny cements the early branching in complex thalloids. Marchantia is supported in one of the earliest divergent lineages. The rate of evolution in organellar loci is slower than for other liverwort lineages, except for two annual lineages. Most genera diverged in the Cretaceous. Marchantia polymorpha diversified in the Late Miocene, giving a minimum age estimate for the evolution of its sex chromosomes. The complex thalloid ancestor, excluding Blasiales, is reconstructed as a plant with a carpocephalum, with filament-less air chambers opening via compound pores, and without pegged rhizoids. Our comprehensive study of the group provides a temporal framework for the analysis of the evolution of critical traits essential for plants during land colonization.

Identification of miRNAs and Their Targets in the Liverwort Marchantia polymorpha by Integrating RNA-Seq and Degradome Analyses

Bryophytes (liverworts, hornworts and mosses) comprise the three earliest diverging lineages of land plants (embryophytes). Marchantia polymorpha, a complex thalloid Marchantiopsida liverwort that has been developed into a model genetic system, occupies a key phylogenetic position. Therefore, M. polymorpha is useful in studies aiming to elucidate the evolution of gene regulation mechanisms in plants. In this study, we used computational, transcriptomic, small RNA and degradome analyses to characterize microRNA (miRNA)-mediated pathways of gene regulation in M. polymorpha. The data have been integrated into the open access ContigViews-miRNA platform for further reference. In addition to core components of the miRNA pathway, 129 unique miRNA sequences, 11 of which could be classified into seven miRNA families that are conserved in embryophytes (miR166a, miR390, miR529c, miR171-3p, miR408a, miR160 and miR319a), were identified. A combination of computational and degradome analyses allowed us to identify and experimentally validate 249 targets. In some cases, the target genes are orthologous to those of other embryophytes, but in other cases, the conserved miRNAs target either paralogs or members of different gene families. In addition, the newly discovered Mpo-miR11707.1 and Mpo-miR11707.2 are generated from a common precursor and target MpARGONAUTE1 (LW1759). Two other newly discovered miRNAs, Mpo-miR11687.1 and Mpo-miR11681.1, target the MADS-box transcription factors MpMADS1 and MpMADS2, respectively. Interestingly, one of the pentatricopeptide repeat (PPR) gene family members, MpPPR_66 (LW9825), the protein products of which are generally involved in various steps of RNA metabolism, has a long stem-loop transcript that can generate Mpo-miR11692.1 to autoregulate MpPPR_66 (LW9825) mRNA. This study provides a foundation for further investigations of the RNA-mediated silencing mechanism in M. polymorpha as well as of the evolution of this gene silencing pathway in embryophytes.

Molecular and Morphological Evidence Challenges the Records of the Extant Liverwort Ptilidium pulcherrimum in Eocene Baltic Amber

Preservation of liverworts in amber, a fossilized tree resin, is often exquisite. Twenty-three fossil species of liverworts have been described to date from Eocene (35–50 Ma) Baltic amber. In addition, two inclusions have been assigned to the extant species Ptilidium pulcherrimum (Ptilidiales or Porellales). However, the presence of the boreal P. pulcherrimum in the subtropical or warm-temperate Baltic amber forest challenges the phytogeographical interpretation of the Eocene flora. A re-investigation of one of the fossils believed to be P. pulcherrimum reveals that this specimen in fact represents the first fossil evidence of the genus Tetralophozia, and thus is re-described here as Tetralophozia groehnii sp. nov. A second fossil initially assigned to P. pulcherrimum is apparently lost, and can be reassessed only based on the original description and illustrations. This fossil is morphologically similar to the extant North Pacific endemic Ptilidium californicum, rather than P. pulcherrimum. Divergence time estimates based on chloroplast DNA sequences provide evidence of a Miocene origin of P. pulcherrimum, and thus also argue against the presence of this taxon in the Eocene. Ptilidium californicum originated 25–43 Ma ago. As a result, we cannot rule out that the Eocene fossil belongs to P. californicum. Alternatively, the fossil might represent a stem lineage element of Ptilidium or an early crown group species with morphological similarities to P. californicum.

The resurrection of Neohattoria Kamim. (Jubulaceae, Marchantiophyta): a six decade systematic conflict resolved through a molecular perspective

The systematic placement of Frullaniaherzogii has been contentious since its description six decades ago. Over the years it has been interpreted as either a member of the genus Frullania or segregated into its own genus, Neohattoria, due to morphological similarities with both Frullania and Jubula. Here we provide molecular evidence that supports the recognition of the genus Neohattoria and its inclusion within the Jubulaceae, together with Jubula and Nipponolejeunea. Jubulaceae are placed sister to Lejeuneaceae rather than to the monogeneric Frullaniaceae.

A new species of Brevianthus (Brevianthaceae, Marchantiophyta) from New Caledonia with unusual underleaf production

Brevianthus is a distinctive genus of leafy liverwort in its succubously inserted, entire leaves, lack of underleaves, restriction of sexual organs to lateral-intercalary branches, scattered rhizoids and dense leaf-surface ornamentation. The sole species, Brevianthusflavus, is divided into two subspecies, one in Tasmania the other in New Zealand. A second species, Brevianthushypocanthidium, is described as new and is the first record of the genus for New Caledonia. Among its distinguishing characters are its shallowly bilobed leaves, and triangular underleaves present on small to medium-sized shoot sectors, the lack of a hyaline leaf margin, and the crenulate leaf margin formed by heavily thickened external cell walls. The most unusual features of the new species are the presence of underleaves between lateral leaf insertion lines that reach the ventral stem mid-line, and the absence of underleaves from larger shoots. To explain these features we propose a competitive model of shoot formation wherein the ventral merophyte progressively loses vigor as its relative stature decreases, and its derivative cells become discontinuous and isolated along the ventral stem surface, with intervening areas occupied by derivatives of the more vigorous lateral merophytes.