Liverwort oil bodies: diversity, biochemistry, and molecular cell biology of the earliest secretory structure of land plants.

Liverworts are known for their large chemical diversity. Much of this diversity is synthesized and enclosed within oil bodies (OBs), a synapomorphy of the lineage. OBs contain the enzymes to biosynthesize and store large quantities of sesquiterpenoids and other compounds while limiting their cytotoxicity. Recent important biochemical and molecular discoveries related to OB formation, diversity, and biochemistry allow comparison with other secretory structures of land plants from an evo-devo perspective. This review addresses and discusses the most recent advances in OB origin, development, and function towards understanding the importance of these organelles in liverwort physiology and adaptation to changing environments. Our mapping of OB types and chemical compounds to the current liverwort phylogeny suggests that OBs were present in the most recent common ancestor of liverworts, supporting that OBs evolved as the first secretory structures in land plants. Yet, we require better sampling to define the macroevolutionary pattern governing the ancestral type of OB. We conclude that current efforts to find molecular mechanisms responsible for the morphological and chemical diversity of secretory structures will help understand the evolution of each major group of land plants, and open new avenues in biochemical research on bioactive compounds in bryophytes and vascular plants.

Defying death: incorporating fossils into the phylogeny of the complex thalloid liverworts (Marchantiidae, Marchantiophyta) confirms high order clades but reveals discrepancies in family-level relationships.

In recent years, the use of extensive molecular and morphological datasets has clarified the phylogenetic relationships among the orders of complex thalloid liverworts (Marchantiidae). However, previous studies excluded extinct taxa; thereby, undersampling the actual taxonomic diversity of the group. Here, we conducted a total-evidence analysis of Marchantiidae incorporating fossils. The combined dataset consisted of 11 genes-sampled from the nuclear, mitochondrial and plastid genomes-and 128 morphological characters. Sixty-two species, representing all classes and orders within Marchantiophyta and genera within Marchantiidae were included in the analyses. Six fossils were scored from literature: two assigned to the outgroup (Metzgeriothallus sharonae and Pallaviciniites sandaolingensis) and four to the ingroup (Marchantites cyathodoides, M. huolinhensis, Ricciopsis ferganica and R. sandaolingensis). Tree searches were conducted using parsimony as the optimality criterion. Clade sensitivity was assessed across a wide range of weighting regimes. Also, we evaluated the influence of fossils on the inferred topologies and branch support. Our results were congruent with previously inferred clades above the order level: Neohodgsoniales was sister to a clade formed by Sphaerocarpales and Marchantiales. However, relationships among families within Marchantiales contradicted recent studies. For instance, a clade consisting of Monosoleniaceae, Wiesnerellaceae and Targioniaceae was sister to the morphologically simple taxa instead of being nested within them as in previous studies. Novel synapomorphies were found for several clades within Marchantiales. Outgroup fossils were more influential than Marchantiidae fossils on overall topologies and branch support values. Except for a single weighting scheme, sampling continuous characters and down-weighting characters improved fossil stability. Ultimately, our results challenge the widespread notion that bryophyte fossils are problematic for phylogenetic inference.

Phylogenomics of a new fungal phylum reveals multiple waves of reductive evolution across Holomycota.

Compared to multicellular fungi and unicellular yeasts, unicellular fungi with free-living flagellated stages (zoospores) remain poorly known and their phylogenetic position is often unresolved. Recently, rRNA gene phylogenetic analyses of two atypical parasitic fungi with amoeboid zoospores and long kinetosomes, the sanchytrids Amoeboradix gromovi and Sanchytrium tribonematis, showed that they formed a monophyletic group without close affinity with known fungal clades. Here, we sequence single-cell genomes for both species to assess their phylogenetic position and evolution. Phylogenomic analyses using different protein datasets and a comprehensive taxon sampling result in an almost fully-resolved fungal tree, with Chytridiomycota as sister to all other fungi, and sanchytrids forming a well-supported, fast-evolving clade sister to Blastocladiomycota. Comparative genomic analyses across fungi and their allies (Holomycota) reveal an atypically reduced metabolic repertoire for sanchytrids. We infer three main independent flagellum losses from the distribution of over 60 flagellum-specific proteins across Holomycota. Based on sanchytrids' phylogenetic position and unique traits, we propose the designation of a novel phylum, Sanchytriomycota. In addition, our results indicate that most of the hyphal morphogenesis gene repertoire of multicellular fungi had already evolved in early holomycotan lineages.

Chemotaxonomy and cytotoxicity of the liverwort Porella Viridissima.

The first chemotaxonomic study based on volatile components of Porella viridissima (Mitt.) Grolle is reported. The GC-MS analysis of ether extract was performed; ten santalane and five pinguisane-type sesquiterpenes were identified together with perrottetianal A as major diterpene. Most of detected santalane-type sesquiterpenes are reported for the first time in liverwort. P. viridissima was found to belong to the chemotype III (pinguisane/sacculatane) and shared chemical similarities with P. navicularis. Perrotettianal A was isolated and has shown strong cytotoxicity against ovarian cancer.

An ancient tropical origin, dispersals via land bridges and Miocene diversification explain the subcosmopolitan disjunctions of the liverwort genus Lejeunea.

Understanding the biogeographical and diversification processes explaining current diversity patterns of subcosmopolitan-distributed groups is challenging. We aimed at disentangling the historical biogeography of the subcosmopolitan liverwort genus Lejeunea with estimation of ancestral areas of origin and testing if sexual system and palaeotemperature variations can be factors of diversification. We assembled a dense taxon sampling for 120 species sampled throughout the geographical distribution of the genus. Lejeunea diverged from its sister group after the Paleocene-Eocene boundary (52.2 Ma, 95% credibility intervals 50.1-54.2 Ma), and the initial diversification of the crown group occurred in the early to middle Eocene (44.5 Ma, 95% credibility intervals 38.5-50.8 Ma). The DEC model indicated that (1) Lejeunea likely originated in an area composed of the Neotropics and the Nearctic, (2) dispersals through terrestrial land bridges in the late Oligocene and Miocene allowed Lejeunea to colonize the Old World, (3) the Boreotropical forest covering the northern regions until the late Eocene did not facilitate Lejeunea dispersals, and (4) a single long-distance dispersal event was inferred between the Neotropics and Africa. Biogeographical and diversification analyses show the Miocene was an important period when Lejeunea diversified globally. We found slight support for higher diversification rates of species with both male and female reproductive organs on the same individual (monoicy), and a moderate positive influence of palaeotemperatures on diversification. Our study shows that an ancient origin associated with a dispersal history facilitated by terrestrial land bridges and not long-distance dispersals are likely to explain the subcosmopolitan distribution of Lejeunea. By enhancing the diversification rates, monoicy likely favoured the colonisations of new areas, especially in the Miocene that was a key epoch shaping the worldwide distribution.

Molecular delimitation of European leafy liverworts of the genus Calypogeia based on plastid super-barcodes.

BACKGROUND: Molecular research revealed that some of the European Calypogeia species described on the basis of morphological criteria are genetically heterogeneous and, in fact, are species complexes. DNA barcoding is already commonly used for correct identification of difficult to determine species, to disclose cryptic species, or detecting new taxa. Among liverworts, some DNA fragments, recommend as universal plant DNA barcodes, cause problems in amplification. Super-barcoding based on genomic data, makes new opportunities in a species identification. RESULTS: On the basis of 22 individuals, representing 10 Calypogeia species, plastid genome was tested as a super-barcode. It is not effective in 100%, nonetheless its success of species discrimination (95.45%) is still conspicuous. It is not excluded that the above outcome may have been upset by cryptic speciation in C. suecica, as our results indicate. Having the sequences of entire plastomes of European Calypogeia species, we also discovered that the ndhB and ndhH genes and the trnT-trnL spacer identify species in 100%. CONCLUSIONS: This study shows that even if a super-barcoding is not effective in 100%, this method does not close the door to a traditional single- or multi-locus barcoding. Moreover, it avoids many complication resulting from the need to amplify selected DNA fragments. It seems that a good solution for species discrimination is a development of so-called "specific barcodes" for a given taxonomic group, based on plastome data.

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

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

The identification and functional characterization of three liverwort class I O-methyltransferases.

Previously it has been shown that the caffeoyl coenzyme A O-methyltransferase (CCoAOMT) type enzyme PaF6OMT, synthesized by the liverwort Plagiochasma appendiculatum Lehm. & Lindenb., (Aytoniaceae), interacts preferentially with 6-OH flavones. To clarify the biochemistry and evolution of liverwort OMTs, a comparison was made between the nucleotide sequence and biological activity of PaF6OMT and those of three of its homologs MpOMT1 (from Marchantia paleacea Bertol., (Marchantiaceae)), MeOMT1 (Marchantia emarginata Reinw et al., (Marchantiaceae)) and HmOMT1 (Haplomitrium mnioides (Lindb.) Schust., (Haplomitriaceae)). The four genes shared >60% level of sequence identity with one another but a <20% level of similarity with typical CCoAOMT or CCoAOMT-like sequences; they clustered with genes encoding animal catechol methyltransferases. The recombinant OMTs recognized phenylpropanoids, flavonoids and coumarins as substrates, but not catechol. MpOMT1 and PaF6OMT exhibited some differences with respect to their substrate preference, and the key residues underlying this preference were identified using site-directed mutagenesis. The co-expression of MpOMT1 and the Arabidopsis thaliana gene encoding S-adenosyl-L-methionine synthase in Escherichia coli was shown to be an effective means of enhancing the production of the pharmacologically active compounds scopoletin and oroxylin A. Liverwort OMTs are thought likely to represent an ancestral out-group of bona fide higher plant CCoAOMTs in evolution and have the potential to be exploited for the production of methylated flavones and coumarins.

Does Calypogeia azurea (Calypogeiaceae, Marchantiophyta) occur outside Europe? Molecular and morphological evidence.

Oil bodies are the unique feature of most liverworts. Their shape, color and distribution pattern in leaf and underleaf cells are important taxonomic features of the genus Calypogeia. Most species of the genus Calypogeia have pellucid and colorless oil bodies, whereas colored, including gray to pale brown, purple-brown or blue oil bodies, are rare. To date, C. azurea was the only species with blue oil bodies to have been considered as a species of the Holarctic range. This species has been noted in various parts of the northern hemisphere-from North America, through Europe to the Far East. The aim of this study was to determine the genetic diversity of C. azurea from different parts of its distribution range and to ascertain whether blue oil bodies appeared once or several times in the evolution of the genus Calypogeia. The phylogenetic analyses based on four plastid regions (rbcL, trnG, trnL, trnH-psbA) and one nuclear region (ITS2) revealed that C. azurea is presently a paraphyletic taxon, with other Calypogeia species nested among C. azurea accessions that were clustered into four different clades. Based on the level of genetic divergence (1.03-2.17%) and the observed morphological, ecological and geographical differences, the evaluated clades could be regarded as previously unrecognized species. Four species were identified: C. azurea Stotler & Crotz (a European species corresponding to the holotype), two new species from Pacific Asia-C. orientalis Buczkowska & Bakalin and C. sinensis Bakalin & Buczkowska, and a North American species which, due to the lack of identifiable morphological features, must be regarded as the cryptic species of C. azurea with a provisional name of C. azurea species NA.

Diversity of Secondary Metabolites in the Liverwort Syzygiella rubricaulis (Nees) Stephani (Jamesoniellaceae, Marchantiophyta) from Neotropical High Mountains.

In order to evaluate the chemical diversity of Syzygiella rubricaulis (Nees) Stephani, a species with a disjunct distribution in the neotropical high mountains, a phytochemical study was carried out with samples from 12 different populations of different altitudes from four South American countries. The chemical profiles of lipophilic extracts were analyzed by GC/MS for each population and 50 different compounds were found with the predominance and richness of sesquiterpenes. The majority of the compounds were found only in one population and the total number of substances ranged from 1 to 15 among the populations, but these numbers were not correlated with altitude, and characterize each population as distinct, based on similarity analysis. The qualitative and quantitative variations of metabolites found are a response to different conditions, under which they live, mostly likely altitudinal conditions. Further studies on the quantification of these chemicals may provide information on their ecological roles and importance for the distribution of S. rubricaulis at different altitudes. Despite the known richness of secondary metabolites produced by bryophytes, they are still poorly explored in the context of the ecological expressions.

Chemotypes and Biomarkers of Seven Species of New Caledonian Liverworts from the Bazzanioideae Subfamily.

Volatile components of seven species of the Bazzanioideae sub-family (Lepidoziaceae) native to New Caledonia, including three endemic species (Bazzania marginata, Acromastigum caledonicum and A. tenax), were analyzed by GC-FID-MS in order to index these plants to known or new chemotypes. Detected volatile constituents in studied species were constituted mainly by sesquiterpene, as well as diterpene compounds. All so-established compositions cannot successfully index some of them to known chemotypes but afforded the discovery of new chemotypes such as cuparane/fusicoccane. The major component of B. francana was isolated and characterized as a new zierane-type sesquiterpene called ziera-12(13),10(14)-dien-5-ol (23). In addition, qualitative intraspecies variations of chemical composition were very important particularly for B. francana which possessed three clearly defined different compositions. We report here also the first phytochemical investigation of Acromastigum species. Moreover, crude diethyl ether extract of B. vitatta afforded a new bis(bibenzyl) called vittatin (51), for which a putative biosynthesis was suggested.

DNA barcoding, ecology and geography of the cryptic species of Aneura pinguis and their relationships with Aneura maxima and Aneura mirabilis (Metzgeriales, Marchantiophyta).

Aneura pinguis is a thalloid liverwort species with broad geographical distribution. It is composed of cryptic species, however, the number of cryptic species within A. pinguis is not known. Five cpDNA regions (matK, rbcL, rpoC1, trnH-psbA and trnL-trnF) and the entire nuclear ITS region were studied in 130 samples of A. pinguis from different geographical regions. The relationships between the cryptic species of A. pinguis, A. maxima and A. mirabilis were analyzed. All of the examined samples were clustered into 10 clades corresponding to 10 cryptic species of A. pinguis (marked A to J). Aneura mirabilis and A. maxima were nested among different cryptic species of A. pinguis, which indicates that A. pinguis is a paraphyletic taxon. Subgroups were found in cryptic species A, B, C and E. As single barcodes, all tested DNA regions had 100% discriminant power and fulfilled DNA barcode criteria for species identification; however, the only combination detected in all subgroups was trnL-trnF with trnH-psbA or ITS2. The distances between cryptic species were 11- to 35-fold higher than intraspecific distances. In all analyzed DNA regions, the distances between most pairs of cryptic A. pinguis species were higher than between A. maxima and A. mirabilis. All cryptic species of A. pinguis clearly differed in their habitat preferences, which suggests that habitat adaptation could be the main driving force behind cryptic speciation within this taxon.

The extraordinary variation of the organellar genomes of the Aneura pinguis revealed advanced cryptic speciation of the early land plants.

Aneura pinguis is known as a species complex with several morphologically indiscernible species, which are often reproductively isolated from each other and show distinguishable genetic differences. Genetic dissimilarity of cryptic species may be detected by genomes comparison. This study presents the first complete sequences of chloroplast and mitochondrial genomes of six cryptic species of A. pinguis complex: A. pinguis A, B, C, E, F, J. These genomes have been compared to each other in order to reconstruct phylogenetic relationships and to gain better understanding of the evolutionary process of cryptic speciation in this complex. The chloroplast genome with the nucleotide diversity 0.05111 and 1537 indels is by far more variable than mitogenome with π value 0.00233 and number of indels 1526. Tests of selection evidenced that on about 36% of chloroplast genes and on 10% of mitochondrial genes of A. pinguis acts positive selection. It suggests an advanced speciation of species. The phylogenetic analyses based on genomes show that A. pinguis is differentiated and forms three distinct clades. Moreover, on the cpDNA trees, Aneura mirabilis is nested among the cryptic species of A. pinguis. This indicates that the A. pinguis cryptic species do not derive directly from one common ancestor.

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.

Conserved Organisation of 45S rDNA Sites and rDNA Gene Copy Number among Major Clades of Early Land Plants.

Genes encoding ribosomal RNA (rDNA) are universal key constituents of eukaryotic genomes, and the nuclear genome harbours hundreds to several thousand copies of each species. Knowledge about the number of rDNA loci and gene copy number provides information for comparative studies of organismal and molecular evolution at various phylogenetic levels. With the exception of seed plants, the range of 45S rDNA locus (encoding 18S, 5.8S and 26S rRNA) and gene copy number variation within key evolutionary plant groups is largely unknown. This is especially true for the three earliest land plant lineages Marchantiophyta (liverworts), Bryophyta (mosses), and Anthocerotophyta (hornworts). In this work, we report the extent of rDNA variation in early land plants, assessing the number of 45S rDNA loci and gene copy number in 106 species and 25 species, respectively, of mosses, liverworts and hornworts. Unexpectedly, the results show a narrow range of ribosomal locus variation (one or two 45S rDNA loci) and gene copies not present in vascular plant lineages, where a wide spectrum is recorded. Mutation analysis of whole genomic reads showed higher (3-fold) intragenomic heterogeneity of Marchantia polymorpha (Marchantiophyta) rDNA compared to Physcomitrella patens (Bryophyta) and two angiosperms (Arabidopsis thaliana and Nicotiana tomentosifomis) suggesting the presence of rDNA pseudogenes in its genome. No association between phylogenetic position, taxonomic adscription and the number of rDNA loci and gene copy number was found. Our results suggest a likely evolutionary rDNA stasis during land colonisation and diversification across 480 myr of bryophyte evolution. We hypothesise that strong selection forces may be acting against ribosomal gene locus amplification. Despite showing a predominant haploid phase and infrequent meiosis, overall rDNA homogeneity is not severely compromised in bryophytes.

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.

[Conserved motifs in the primary and secondary ITS1 structures in bryophytes].

A study of the ITS1 nucleotide sequences of 1000 moss species of 62 families, 11 liverwort species from five orders, and one hornwort Anthoceros agrestis identified five highly conserved motifs (CM1-CM5), which are presumably involved in pre-rRNA processing. Although the ITS1 sequences substantially differ in length and the extent of divergence, the conserved motifs are found in all of them. ITS1 secondary structures were constructed for 76 mosses, and main regularities at conserved motif positioning were observed. The positions of processing sites in the ITS1 secondary structure of the yeast Saccharomyces cerevisiae were found to be similar to the positions of the conserved motifs in the ITS1 secondary structures of mosses and liverworts. In addition, a potential hairpin formation in the putative secondary structure of a pre-rRNA fragment was considered for the region between ITS1 CM4-CM5 and a highly conserved region between hairpins 49 and 50 (H49 and H50) of the 18S rRNA.

Sesquiterpene hydrocarbons from the liverwort Treubia isignensis var. isignensis with chemotaxonomic significance.

The New Caledonian endemic Treubia isignensis var. isignensis, which is known as a morphologically primitive liverwort, was extracted with diethyl ether and the crude extract analyzed by TLC and GC/MS. The species is chemically very primitive because it produces only maaliane-, eudesmane-, aristolane and gorgonane sesquiterpene hydrocarbons, which are significant chemical marker of the species; neither oxygenated terpenoids nor aromatic compounds were detected.

Epiphytic leafy liverworts diversified in angiosperm-dominated forests.

Recent studies have provided evidence for pulses in the diversification of angiosperms, ferns, gymnosperms, and mosses as well as various groups of animals during the Cretaceous revolution of terrestrial ecosystems. However, evidence for such pulses has not been reported so far for liverworts. Here we provide new insight into liverwort evolution by integrating a comprehensive molecular dataset with a set of 20 fossil age constraints. We found evidence for a relative constant diversification rate of generalistic liverworts (Jungermanniales) since the Palaeozoic, whereas epiphytic liverworts (Porellales) show a sudden increase of lineage accumulation in the Cretaceous. This difference is likely caused by the pronounced response of Porellales to the ecological opportunities provided by humid, megathermal forests, which were increasingly available as a result of the rise of the angiosperms.