- 1The New Zealand Institute for Plant and Food Research Limited, Palmerston North, New Zealand
- 2Faculty of Agriculture and Life Sciences, Lincoln University, Christchurch, New Zealand
- 3School of Biological Sciences, Monash University, Melbourne, VIC, Australia
The flavonoid pathway is one of the best characterized specialized metabolite pathways of plants. In angiosperms, the flavonoids have varied roles in assisting with tolerance to abiotic stress and are also key for signaling to pollinators and seed dispersal agents. The pathway is thought to be specific to land plants and to have arisen during the period of land colonization around 550–470 million years ago. In this review we consider current knowledge of the flavonoid pathway in the bryophytes, consisting of the liverworts, hornworts, and mosses. The pathway is less characterized for bryophytes than angiosperms, and the first genetic and molecular studies on bryophytes are finding both commonalities and significant differences in flavonoid biosynthesis and pathway regulation between angiosperms and bryophytes. This includes biosynthetic pathway branches specific to each plant group and the apparent complete absence of flavonoids from the hornworts.
Opening Comments
The flavonoid pathway is a core component enabling land plants to interact with their environment. Flavonoids have been demonstrated to assist in tolerance to both abiotic and biotic stresses, and their production can be induced by cold, UV-B light (UVB) or strong white light, nutrient deprivation, desiccation, salinity, metal toxicity, and pest and pathogen attack (Agati and Tattini, 2010; Cheynier et al., 2013; Landi et al., 2015; Davies et al., 2018). In angiosperms, flavonoids are also key for signaling to pollinators and seed dispersal agents. Flavonoids are formed within the larger phenylpropanoid pathway, starting with chalcones as the first flavonoids. During plant evolution, flavonoid pathway diversity has increased greatly, with more than 8,000 different structures reported from the relatively small number of plant species studied to date (Andersen and Markham, 2006). This multiplicity of structures and functions is thought to have assisted land plants to colonize the wide range of environments they now occupy. The over 8,000 compounds are grouped into relatively few classes of flavonoids, based on the core structure and/or biosynthetic origin. The major flavonoid classes are the flavones, flavonols, isoflavonoids, aurones, 3-deoxyanthocyanins, anthocyanins, proanthocyanidins (condensed tannins), and the recently reported “auronidins” (Berland et al., 2019) (Figures 1 and 2). However, there are also notable groups of related non-flavonoid compounds produced with the same starting precursors as used for chalcone formation, such as the stilbenes and bibenzyls. Most flavonoids are targeted to the vacuole as water-soluble glycosides, although some are transported to the cell wall or are released from the plant to the environment. Many can absorb light in the UV-spectrum, while anthocyanins and auronidins provide colored pigments that can screen in the visible part of the light spectrum (Lee and Gould, 2002; Landi et al., 2015; Berland et al., 2019).
Figure 1 Examples of compounds from phenylpropanoid pathway branches discussed in the text. Bibenzyls, stilbenes, and lignans are phenylpropanoids formed from branches prior to the start of the specific flavonoid pathway branch. The example compounds are from liverworts, angiosperms, and hornworts, respectively. All the other compounds shown are flavonoids.
Figure 2 Schematic of a section of the phenylpropanoid pathway leading to the production of major compound groups discussed in the article. Enzyme abbreviations are: PAL, phenylalanine ammonia-lyase; C4H, cinnamate 4-hydroxylase; 4CL, 4-coumarate:CoA ligase; CHS, chalcone synthase; CHI, chalcone isomerase; CHIL, chalcone isomerase-like; F3H, flavanone 3-hydroxylase; DFR, dihydroflavonol 4-reductase; FNR, flavanone reductase; ANS, anthocyanidin synthase; F3GT, flavonoid 3-O-glucosyltransferase; A5GT, anthocyanidin 5-O-lucosyltransferase; AUS/AS, aureusidin/aurone synthase; FNS, flavone synthase (2OGD = I, Cyp450 = II); FLS, flavonol synthase.
Flavonoid biosynthesis is frequently considered unique to land plants. Some charophycean algae do tolerate terrestrial conditions (Karsten and Holzinger, 2014; Holzinger and Pichrtová, 2016) and can produce polyphenolics in response to abiotic stresses such as UVB, salinity, or dehydration. Indeed, one hypothesis is that land plants arose from algal ancestors that were already terrestrial (Harholt et al., 2016) and therefore may have had some of the biosynthetic pathways characteristic of land plants. However, there are as yet no convincing reports of the flavonoid pathway branch existing outside of land plants, with the possible exception of some fungi. A common adaptation for UVB tolerance in extant algae is the production of mycosporine-like amino acids (MAA). However, while MAA production is frequently reported for red algae (Rhodophyta) and other marine organisms, there are few reports for the Chlorophyta (a green algal clade that with Streptophyta form the Viridiplantae) and none for the charophyte algae that is thought to be most closely related to the land plant ancestor. Although genes for the initial steps of the phenylpropanoid pathway may be present in the genomes of extant Rhodophyta, Glaucophyta, Chlorophyta, and charophytes (Labeeuw et al., 2015; de Vries et al., 2017; Davies et al., 2020), there are no substantiated examples of flavonoid-specific genes being identified. Indeed, many metabolite studies on algae report only amounts of “total flavonoids” measured using general assays. Nevertheless, there are some reports detailing specific flavonoid structures in algal preparations, including compounds such as chalcones, flavones, flavonols, isoflavonoids, and proanthocyanidins (Klejdus et al., 2010; Goiris et al., 2014; El Shoubaky et al., 2016; Agregán et al., 2017; Ben Saad et al., 2017). Although the concentrations of flavonoid compounds reported in most examples are extremely low compared with those commonly found in land plants (ng gDW−1 amounts compared with mg gDW−1), these reports mean that the presence of a biosynthetic pathway to flavonoids in algae cannot be ruled out.
The phenolics most commonly produced by algae are the phlorotannins, diverse oligomers derived from phloroglucinol that are found in brown algae (Heterokonts) (Imbs and Zvyagintseva, 2018). In Ectocarpus siliculosus, the phloroglucinol precursor is formed by the condensation of malonyl-CoA by a polyketide synthase (PKS) (Meslet-Cladière et al., 2013). This is analogous to the action of CHALCONE SYNTHASE (CHS) in the condensation of malonyl-CoA with a p-coumaroyl-CoA starter molecule for flavonoid biosynthesis. However, given the very large phylogenetic distance between brown algae and land plants this may be an example of parallel evolution. Notably, in fungal species there have been recent reports of both the production of flavonoids and the presence of genes with significant sequence similarity to those of the phenylpropanoid pathway of land plants. Findings include the detection of a range of phenylpropanoids, including flavonoids, in Fusarium (Bilska et al., 2018) and detection of flavonoids and candidate genes for stilbene production in Alternaria (Lu et al., 2019). The significance of phenylpropanoid biosynthesis being present in fungi for the current proposals on the evolutionary origins of flavonoid biosynthesis has yet to be addressed. The separation of the fungi and the algae/plant ancestors is thought to be an ancient event, preceding the divergence fungi and animals (Burki, 2014; Burki et al., 2020).
Origins and Vegetative Functions of the Flavonoid Pathway
Regardless of whether the ancestral genes for flavonoid biosynthesis were present in algal ancestors, the flavonoid pathway we see in extant land plants is hypothesized to have arisen when the land plant ancestors were first colonizing the land about 550–470 million years ago (MYA) (Markham, 1988; Stafford, 1991; Jorgensen, 1994; Koes et al., 1994; Kenrick and Crane, 1997; Rozema et al., 2002). Two major hypotheses have been proposed for the initial role of flavonoids. Firstly, that flavonoids may have helped in coping with the additional abiotic stresses resulting from a terrestrial lifestyle, in particular increased exposure to UVB, but potentially also drought and extreme temperature fluctuations (Markham, 1988; Jorgensen, 1994; Kenrick and Crane, 1997; Cockell and Knowland, 1999; Rozema et al., 2002; Ligrone et al., 2012; Mouradov and Spangenberg, 2014; Demarsy et al., 2017; Davies et al., 2018; de Vries and Archibald, 2018; Rensing, 2018). The alternative proposal is that flavonoids arose as physiological regulators or chemical messengers. This was outlined in Stafford (1991) for the regulation of auxin action, with signaling to mycorrhizal and symbiotic fungi proposed as possible additional communication functions. It was argued that flavonoids would probably have been present at only low concentrations when the pathway first evolved, limiting their efficacy as UVB-screening compounds. This was in the context of other arguments against the need for flavonoids as UVB-screens, such as the effective UVB-screening properties of non-flavonoid phenylpropanoids like the hydroxycinnamic acids (HCAs). More recently, arguments for the early functions of flavonoids being other than UVB-screening have been extended by discoveries on their antioxidant properties and possible signaling actions through the redox pathway or by affecting H2O2 retrograde signals between the chloroplast and nucleus (Taylor and Grotewold, 2005; Agati and Tattini, 2010; Pollastri and Tattini, 2011; Agati et al., 2012; Brunetti et al., 2018; Foyer, 2018; Muhlemann et al., 2018; Brunetti et al., 2019).
That flavonoids affect auxin transport has now been demonstrated in several angiosperm species, including Arabidopsis, apple, and tomato, by analysis of mutants or transgenic lines with reduced flavonoid biosynthesis (Brown et al., 2001; Buer and Muday, 2004; Taylor and Grotewold, 2005; Peer and Murphy, 2007; Dare and Hellens, 2013; Maloney et al., 2014). Altered developmental traits in such plants include dwarfing, loss of pollen fertility and altered root development and gravitropic responses (Van der Meer et al., 1992; Napoli et al., 1999; Brown et al., 2001; Dare and Hellens, 2013; Maloney et al., 2014; Muhlemann et al., 2018). The phenotypes observed vary between species, for example the complete loss of flavonoids in the Arabidopsis chs mutant affects root patterning but not pollen viability (Burbulis et al., 1996; Ylstra et al., 1996). To date, the great majority of data are for angiosperms, and studies on other plant groups are required to determine whether these flavonoid functions are shared across land plants and so may have a common evolutionary origin in the early land plant ancestor. Differential distribution of flavonols and auxin has been observed accompanying stem reorientation in a gymnosperm (Ramos et al., 2016), supporting a conserved function in auxin transport within seed plants. If the genetic mechanisms involved are also conserved, then that would support an evolutionary origin before 350–300 MYA. However, while definitive experiments on flavonoids and hormone function have not been conducted in bryophytes, indications are that flavonoids are not necessary for normal development in this plant group. No flavonoids have been detected to date in hornworts, and a genetic mutant of the liverwort Marchantia polymorpha lacking flavonoids has normal developmental patterns (Clayton et al., 2018). Addition of phenylpropanoids or phenylpropanoid pathway inhibitors can alter bryophyte development in culture, but whether this is because of altered hormone action has not been tested (Chattopadhyay et al., 2018). Even with additional data it may be difficult to determine the most probable option between developmental roles for flavonoids having been acquired in seed plants since the last common ancestor, or having being present but then lost during subsequent evolution of the bryophytes. Nevertheless, establishing whether flavonoids regulate auxin action in bryophytes is an important goal.
Fossils, such as those found in the Rhynie Chert in Scotland (about 410 MYA) (Ligrone et al., 2012), provide detail of the structure of early land plants but little information on their specialized/secondary metabolism. The presence of specialized biosynthetic or storage structures in fossils, such as possible equivalents to the terpenoid-accumulating oil bodies of extant liverworts (Labandeira and Currano, 2013), can support the presence of specialized metabolite pathways but not provide details of the specific compounds produced. To generate hypotheses on the origins and subsequent evolution of the flavonoid pathway we need to compare the genetics and biochemistry of the pathway across diverse extant plant groups, as this can identify conserved pathway components that may have originated with the last common ancestor. In this respect, bryophytes are of key importance (Figure 3). “Bryophytes” is the collective name for non-vascular land plants, comprising the liverworts (Marchantiophyta, approximately 9,000 species), hornworts (Anthocerotophyta, approximately 300 species) and mosses (Bryophyta, approximately 12,000 species). Evidence such as morphological comparison with the fossil record has placed liverworts as the “sister” group to extant land plants—that is, at the base of the land plant evolutionary tree, making bryophytes paraphyletic. However, DNA sequencing data have suggested alternatives, in particular, either paraphyletic bryophytes with hornworts as a sister group to all other land plants, or a single monophyletic bryophyte clade that is sister to the vascular plants (Wickett et al., 2014; Puttick et al., 2018). It is generally accepted that land plants evolved from an ancestral charophyte, and the extant algal sister group is probably the order Zygnematales or a clade of the Zygnematales and Coleochaetales together (Zhong et al., 2014; Delwiche and Cooper, 2015; de Vries and Archibald, 2018).
Figure 3 Phylogenetic context of land plants. The phylogenetic tree on the left shows the relationships of the animals, fungi, and plants, including the major divisions within the plants. The tree on the right shows the relationships among the land plants (Embryophyta). The phylogeny of the bryophytes is unresolved, but the current proposal of a sister relationship between liverworts and mosses is shown. MYA, million years ago.
Flavonoids and Tolerance to Ultraviolet B Light
The origin of the flavonoid pathway for providing tolerance to UVB has been supported by recent studies on the liverwort species Marchantia. (“Marchantia” is used in this article to refer to M. polymorpha subsp. ruderalis, which is the model experimental species.) Marchantia is an excellent research model: it is small; has a rapid growth rate; can asexually reproduce in large numbers through single-cell-derived clonal gemmae; and, has a small genome (approximately 220 Mb) which, although larger than that of Arabidopsis (at 135 Mb), contains significantly fewer genes (around 19,000 gene models compared with around 28,000 protein-coding genes in Arabidopsis) (Ishizaki et al., 2015; Bowman et al., 2016; Shimamura, 2016; Bowman et al., 2017). It also offers efficient CRISPR/Cas9 mutagenesis in a dominant haploid gametophytic generation (Sugano et al., 2018).
Significantly, the UVB response of Marchantia has many components in common with that of Arabidopsis. Flavonol O-glycosides are key for UVB tolerance of Arabidopsis (Kusano et al., 2011; Morales et al., 2013; Yin and Ulm, 2017), while the related flavone O-glycosides contribute to Marchantia UVB tolerance (Clayton et al., 2018). In both Arabidopsis and Marchantia, mutants with reduced flavonoid production are more easily damaged by UVB, while mutants or transgenics with increased flavonoid content have increased UVB tolerance. Moreover, the signaling pathway for flavonoid pathway induction through the UV RESISTANCE LOCUS8 (UVR8) photoreceptor, the bZIP transcription factor (TF) ELONGATED HYPOCOTYL5 (HY5), and the modifiers of protein stability such as CONSTITUTIVELY PHOTOMORPHOGENIC1 (COP1) and REPRESSOR OF UVB PHOTOMORPHOGENESIS1 (RUP1) is also conserved between the species (Clayton et al., 2018; Kondou et al., 2019). Equivalent functional studies are lacking on other major basal plant groups, such as mosses and lycophytes. However, phenolics do seem to be important for UVB tolerance in mosses (Clarke and Robinson, 2008; Wolf et al., 2010; Waterman et al., 2017; Soriano et al., 2018b); and UVB exposure of the Antarctic moss Pohlia nutans increased transcript abundance for genes of the UVR8 and flavonoid pathways (Li et al., 2019). This suggests that the core UVB protection mechanism of UVR8-induced flavonoid production may already have been established in the last common ancestor of bryophytes and angiosperms. While some of the same genetic components, such as UVR8 and HY5-like genes have been identified in distantly related groups within the Viridiplantae (Allorent et al., 2016; Bowman et al., 2017; Clayton et al., 2018; Soriano et al., 2018a; Kondou et al., 2019), and algae can produce purple phenolic pigments in response to abiotic stress (Aigner et al., 2013; Holzinger and Pichrtová, 2016), UVB-induced flavonoid production has not been characterized in algae. Thus, the UVR8-induction pathway for UVB-absorbing flavonoids (potentially flavone glycosides as the first compounds) may have been a character rapidly acquired during the water-to-land transition.
There are variations among land plants to the UVR8/flavonoid system for providing tolerance to UVB exposure, and further research would be beneficial to establish whether these may also have an early evolutionary origin. In Arabidopsis, the HCA compounds sinapate esters have an important role in UVB-screening. Comparison of Arabidopsis mutants at different biosynthetic steps found sinapate esters have a protective role comparable to, or perhaps more important than, that of flavonols (Li et al., 1993; Landry et al., 1995). Furthermore, a genetic screen for UVB-tolerance genes identified the transcriptional repressor AtMYB4 as being downregulated in response to UVB to facilitate increased sinapate ester production (Jin et al., 2000). The principal absorption maxima from 290 to 330 nm makes HCAs particularly effective UVB screening compounds, and would be a more carbon-efficient screen than the flavonoids. In addition to flavonol glycosides, mosses have been shown to produce biflavonoids or cell wall-bound phenolics in response to UVB exposure (Clarke and Robinson, 2008; Wolf et al., 2010; Waterman et al., 2017; Soriano et al., 2018b). Biflavonoid induction is a notable feature of the moss Ceratodon purpureus, a species found from Antarctica to hot desert environments (Waterman et al., 2017). The production of cell wall-localized phenolics as part of the UVB screening capacity may be a common feature of bryophytes, including for the flavonoid-lacking hornworts (Monforte et al., 2018; Soriano et al., 2018c).
Pigmented Flavonoids and Tolerance to Abiotic Stress
The other major group of flavonoids shown to be involved in tolerance to abiotic stress is the 3-hydroxyanthocyanins (typically cyanidin-derivatives) that are found in gymnosperms and almost all angiosperms (Lee and Gould, 2002; Landi et al., 2015). Under the physiological conditions commonly found in vegetative tissues, these provide red pigmentation. The structurally similar 3-deoxyanthocyanins have been extensively characterized in ferns, and also reported for lycophytes and mosses (Andersen and Markham, 2006). Although the color of 3-deoxyanthocyanins is shifted toward orange compared to the equivalent 3-hydroxyanthocyanins, in vegetative tissues similar red pigmentation typically occurs from both compound types. Additionally, cell wall-bound red flavonoids have been reported from liverworts and mosses: riccionidin (an auronidin) and sphagnorubin, respectively (Vowinkel, 1975; Kunz et al., 1993; Berland et al., 2019). Until recently, riccionidin and sphagnorubin were considered anthocyanidins (the non-glycosylated anthocyanin core molecule) with additional rings. Thus, an evolutionary path could be envisioned of cell-wall bound anthocyanidin being the basal state that could have been present in the last common ancestor, then a progression to 3-deoxyanthocyanins, and then 3-hydroxyanthocyanins. However, now we know that riccionidins represent a separate flavonoid class unrelated to anthocyanins (see later sections for biosynthetic details) it is less clear as to what the class of red flavonoid pigment (if any) was present in the early land plants. Auronidin and anthocyanin biosynthesis may be derived characters in each lineage, or the missing pathway in each linage may have been lost during evolution. Furthermore, no proanthocyanidins have been reported from bryophytes or lycophytes, suggesting this branch of flavonoid biosynthesis probably arose later in vascular plant evolution, since proanthocyanidins are found in extant ferns, gymnosperms, and angiosperms.
With regard to function, given the lack of flowers or enclosed seeds in the early land plants, any initial functions of pigmented flavonoids were probably unrelated to animal interactions. In vascular plants anthocyanins have been linked to improving tolerance to a range of abiotic stresses, as well as some biotic challenges. However, how anthocyanins achieve this in the different stress situations, and whether there is a single mechanism or specific functional variations, is the subject of much research and debate. Anthocyanins, auronidins and sphagnorubins can screen out white light to reduce photooxidative damage, and it has been estimated that anthocyanins could absorb over 40% of photosynthetically active radiation in the range containing the most damaging wavelengths for photoinhibition (Merzlyak and Chivkunova, 2000; Pietrini et al., 2002). However, anthocyanins can simultaneously reduce the cellular stress associated with photooxidation through quenching reactive oxygen species (ROS). The relative importance of these two mechanisms is unresolved, even for extensively studied situations such as the appearance of red anthocyanins during autumn senescence of leaves of deciduous trees. Localization of the pigmented flavonoids in the cell wall and vacuole could be thought to argue for a screening mechanism, as ROS are principally generated in plastids and mitochondria. However, it has been suggested that this does not rule out an antioxidant primary function (Agati et al., 2012). Many of the same arguments for or against light screening versus antioxidant primary functions can also be applied to the function of flavones/flavonols in UVB tolerance (Agati and Tattini, 2010; Agati et al., 2012; Davies et al., 2018). Further complicating the development of a unifying theory for anthocyanin abiotic stress function are alternative hypotheses that involve neither light screening nor ROS scavenging, such as drought tolerance through decreased osmotic potential, increasing light absorption to help warm leaves, providing camouflage against insect herbivores, “honest” signaling to herbivores that leaves contain anti-feedant compounds and/or are about to be shed, and making leaves more noticeable to insect predators (anti-crypsis) (Gould et al., 1995; Lee and Gould, 2002; Manetas, 2006; Archetti, 2009; Hughes, 2011; Agati et al., 2012; Landi et al., 2015; Davies et al., 2018). Additionally, as mentioned earlier, there is also evidence supporting flavonoid roles as signaling molecules (Taylor and Grotewold, 2005; Agati and Tattini, 2010; Agati et al., 2012; Foyer, 2018).
The cell wall-bound nature of auronidins and sphagnorubins complicates the theory of an antioxidant or signaling role for red flavonoids. The cell wall-localization of the red pigments of mosses and liverworts, termed “tissue fixed”, has been known from early studies on Sphagnum and various liverworts (Nagai, 1915; Rudolph et al., 1977). Nevertheless, that the environmental stimuli inducing cell wall-bound red pigments in liverworts and mosses appear similar to those that trigger anthocyanin production in vascular plants, was also noted as far back as 1915 (Nagai, 1915). Kny (1890) and Cesares Gil (1902) noted that M. polymorpha or Reboulia hemisphaerica, respectively, grown in sunny locations produced more pigment than those growing in shady environments. Nagai (1915) was able to show that limiting nitrogen and phosphorus supply intensified the pigmentation of M. polymorpha and M. paleacea, but upon transfer to nutrient-rich media the newly developed tissue lacked pigmentation. Moreover, combinations of stresses that can cause oxidative stress are strong signals for red flavonoid biosynthesis in many bryophytes. Thus, as with angiosperms, in liverworts and mosses cold and light can individually induce reddening, but strong sunlight in cold conditions induces much stronger pigmentation, whether this be at altitude, in the Antarctic, or during the cold nights and bright days of autumn (Gerdol, 1996; Gerdol et al., 1998; Newsham et al., 2005; Hooijmaijers and Gould, 2007; Glime, 2007; Bonnett et al., 2010).
Detailed studies on Marchantia (Albert et al., 2018; Kubo et al., 2018) and Ricciocarpos natans (Kunz and Becker, 1995) found nitrogen deprivation and increased white light exposure induced auronidin accumulation, as has also been shown for anthocyanin accumulation in Arabidopsis and apple (Rubin et al., 2009; Wang et al., 2018). The signaling pathways for anthocyanin induction by nitrogen and phosphorus deficiency are well-characterized for Arabidopsis, with R2R3MYBs being the key activating transcription factors (Lillo et al., 2008; Rubin et al., 2009). Induction of auronidin in Marchantia by nitrogen and phosphorus also requires an R2R3MYB (Kubo et al., 2018), suggesting signaling components may be conserved. For Antarctic liverworts and mosses UVB exposure also induced production of red flavonoids, which most commonly were cell wall-bound (Newsham et al., 2005; Waterman et al., 2018). UVB induces anthocyanin production in some angiosperms, but it is much less common a response than induction of flavones/flavonols. Flavones and flavonols are more effective at screening UVB than anthocyanins, although aromatic acylation can give anthocyanins absorbance maxima in the UV range. The induction of anthocyanins by UVB has thus been suggested to be more for ROS scavenging and/or screening of white light than for UVB screening. In the case of the non-acylated cell wall-bound flavonoid pigments of mosses and liverworts, it seems probable that production is induced to screen white light and prevent further ROS generation, especially as the summer conditions in the Antarctic present a combination of stresses from continuous white light, cold, and drought.
There are other red/purple plant pigments besides the flavonoids able to screen in photosynthetically active wavelengths. Notable among these are the betacyanins, which are produced in the many species of the core Caryophyllales that do not produce anthocyanins (Polturak and Aharoni, 2019), and the phenolic pigments of algae. Zygogonium ericetorum is a charophyte green alga that can grow in alpine environments and when exposed to abiotic stress produces vacuolar-localized purple pigments, thought to be polymers of glucose and gallic acid, which can absorb in both UVB and photosynthetically active wavelengths (Aigner et al., 2013). In brown algae, phlorotannins can accumulate to more than 15% of dry weight (Imbs and Zvyagintseva, 2018). Phlorotannins are highly hydrophilic polymers, and may be cell wall-bound, stored intracellularly in vesicles, or exported.
Progress on determining the biological roles of cell wall-bound pigments in bryophytes has been limited by the lack of genetic systems, and the difficulty of extracting the pigments. However, genetic tools are now available in Marchantia that will allow tests of the functions of the pigments in abiotic or biotic stress tolerance. Mutants are available that have loss of auronidin pigmentation but retain flavone production (Albert et al., 2018; Kubo et al., 2018), have loss of flavone production but retain auronidin pigmentation, or have reduced amounts of both compounds. These have been used for physiological studies with respect to UVB tolerance (Clayton et al., 2018) and pathogen attack (Carella et al., 2019). In angiosperms, a range of flavonoids are localized to the cell wall (Agati et al., 2012), including rare examples of cell wall-bound anthocyanins (Philpott et al., 2009), although the physiological roles of these are generally unclear. The cell wall localization of other phenylpropanoids, in particular HCA derivatives, is common in angiosperms. These may contribute to lignin formation or be accumulated as monomers or dimers in the wall. Besides having structural roles these polymers may also contribute to physical barriers to pathogens (Zhao and Dixon, 2014). Although lignin is thought to be absent from non-vascular plants, cinnamic acid derivatives such as rosmarinic acid and (neo)lignans are common in bryophytes (Asakawa et al., 2013; Asakawa, 2017), and may be cell wall-localized (Wang et al., 2013). In Sphagnum moss, oxidative derivatives of sphagnum acid, p-hydroxyacetophenone, hydroxybutenolide, and p-hydroxybenzoic acid, as well as the phenolics p-coumaric acid and trans-cinnamic acid, were predominantly bound to the cell wall (Verhoeven and Liefveld, 1997). It seems a strong possibility that the red flavonoid pigments of bryophytes contribute, along with the cinnamic acid derivatives, to forming a physical barrier against pathogens. The recent study of Carella et al. (2019) demonstrated that the production of auronidin in Marchantia greatly enhanced resistance to Phytophthora palmivora infection, with a lack of hyphae penetration into the highly pigmented regions of plants. In relation to the mechanism of action, it would be of much interest to determine the nature of the incorporation of auronidin and sphagnorubins into the wall and whether polymerization occurs. Dimers of auronidin/riccionidin A have been isolated (termed riccionidin B) (Kunz et al., 1993), providing a basis for polymerization.
Several thalloid liverwort genera, and many moss genera, have species with considerable drought tolerance, with examples in both plant groups of individuals withstanding continuous desiccation for more than 20 years (Breuil-Sée, 1993; Stark et al., 2017). As the plants of liverwort genera such as Riccia and Targionia dry out, the sides of the thallus roll over the dorsal surface so that it is covered by the darkly pigmented ventral scales and rhizoids [see Figure 4 and Reeb et al. (2018) for examples]. This forms a “capsule” that can recover and renew growth even after extended periods without additional water. The function of the very strong pigmentation of the ventral scales, presumably by cell wall-bound auronidin, is not known. It may provide protection of the DNA against UVB damage during a period when DNA repair mechanisms are not active, given that auronidin accumulation is induced by UVB in some Antarctic species. Alternatively, the modification of the cell wall could prevent pathogen ingress, as demonstrated for Marchantia (Carella et al., 2019), or reduce water loss. Plants of desiccation-tolerant species in leafy liverwort genera such as Herbertus and Cephaloziella also often have dark red pigmentation (Vitt et al., 2012). A related but little studied example first described in the 1890s (Campbell, 1896) is the formation of “tubers” by some liverwort species, notably Geothallus tuberosus. G. tuberosus can form thickened inner regions of the thallus that are presumed to store carbohydrates. As the tubers form, the associated cells become strongly dark red pigmented with “thick walls”. The tubers can become buried in the soil and, although the surrounding plant may die, the thallus and associated meristem survive the long dry season of the Southern Californian regions to which the species is native.
Figure 4 Illustrations of some of the variations in form of the purple-pigmented ventral scales of thalloid liverworts. Left The aquatic form of Ricciocarpos natans (previously Riccia natans) has long ventral scales that extend into the water below the plant. Right Riccia squamata is one of the drought-adapted thalloid liverworts. As the environment dries out it curls over so that the ventral scales encase the thallus. Illustrations are from Lindenberg (1836).
The ventral scales of many thalloid liverworts, frequently strongly pigmented by auronidin, often extend around the apex of the thallus to provide a barrier layer between the meristem and the soil (Figure 5). Protection of the meristem from physical damage and pathogen ingress could explain this pigmentation. This suggestion could be extended to include protection against herbivory. The extended ventral scales of the aquatic form of R. natans also have strong auronidin-based purple pigmentation (Figure 4), so perhaps auronidins contribute to aquatic herbivore deterrence in these cells.
Figure 5 Examples of pigmentation of the thallus and antheridiophores of Marchantia foliacea. Ventral (A, C) and dorsal (B, D) views of a thallus branch. The ventral surface contains the strongly pigmented scales. These extend past the meristematic notch and can be seen from the dorsal surface. The thallus in image (A) shows the strong pigmentation that can also occur in the non-scale cells. Image (D) shows a close up of the meristem region from the dorsal side, with the corresponding ventral view (C) “flipped” vertically to present the same orientation. Image (E) shows the strong pigmentation common for the surface of the antheridiophores.
There are few studies on the biological functions of the 3-deoxyanthocyanins that are common in mosses and ferns, but there is evidence they also are involved in plant defense. Greater amounts of 3-deoxyanthocyanins in fronds of the aquatic fern Azolla correlated with increased feeding deterrence to snails and the tadpoles of frogs (Cohen et al., 2002a). In the same species, 3-deoxyanthocyanins also may promote the establishment of the symbiosis with the cyanobacterium Nostoc (Cohen et al., 2002b). This suggests there are specific biological functions of the flavonoid pigments in different bryophyte and fern species, although the induction of 3-deoxyanthocyanin production in ferns by general abiotic stresses also indicates a general light screening/antioxidant function in common with that of 3-hydroxyanthocyanins in angiosperms. Perhaps the cell wall-bound pigments of the bryophytes have evolved to have elegant multi-functionality, providing abiotic stress tolerance through antioxidant and/or light screening actions and altering physical properties of the cell for biotic stress resistance.
The Evolutionary Significance of the Occurrence of Different Flavonoid Structural Groups
The identification of different flavonoid groups across land plants has been conducted for many years, both to further understand the evolutionary significance of flavonoid distribution by chemotaxonomy and for the discovery of novel bioactives (Markham, 1988; Asakawa, 2017; Jiang et al., 2016; de Vries et al., 2017; Yonekura-Sakakibara et al., 2019). Flavones and/or flavonols are almost ubiquitous across land plants (Berim, 2016) but variations in the specific types of flavonols or flavones produced have occurred during evolution, for example resulting in the rarity of flavone O-glycosides in leafy liverworts (Markham, 1988) or of polymethoxylated flavones in gymnosperms (Berim, 2016). Overall, it is probable that their biosynthesis was acquired very early during land plant evolution as an important stress adaptation. The remarkable exception is the hornworts. Hornworts produce polyphenolics, notably rosmarinic acid and lignan-like compounds (e.g., anthocerotonic acid and megacerotonic acid) (Petersen and Simmonds, 2003; Soriano et al., 2018b) (Figure 1), but there is no report of any flavonoid being found (Markham, 1988). Thus, either hornworts diverged from the last common land plant ancestor before the evolution of the flavonoid pathway, or the ability to make flavonoids was subsequently lost in this lineage. The completion of a hornwort genome sequence (Szövényi et al., 2015) and transcriptomic studies examining land plant evolution (Wickett et al., 2014; Puttick et al., 2018) may provide the data to help in resolving this question. Analysis of the transcriptomic (SRA PRJEB21674) and genomic (SRA ERR771108 and SRR1278954 for Anthoceros agrestis and Anthoceros punctatus, respectively) data currently available on GenBank can identify with confidence hornwort deduced sequences corresponding to the early steps of the phenylpropanoid and flavonoid pathway, including for PHENYLALANINE AMMONIA LYASE (PAL), CINNAMATE 4-HYDROXYLASE (C4H), 4-COUMARATE-COA LIGASE (4CL), CHS, and CHALCONE ISOMERASE-LIKE (CHIL) (Supplementary Table 1). However, additional analysis is required to show whether these produce functional enzymes. No clear CHALCONE ISOMERASE (CHI)-encoding sequence is present in the data, but mosses can make flavonoids without a gene corresponding to the typical CHI (Cheng et al., 2018). Rosmarinic acid is also found in some algae (Agregán et al., 2017), but comparison of the biosynthetic pathways between land plants and algae has not been made.
Until the recent clarification of the riccionidin structures as auronidins, anthocyanins were thought to be present in all extant lineages of land plants except hornworts. A progression in anthocyanin complexity was suggested, with liverworts producing “primitive” anthocyanidins (the non-glycosylated anthocyanin core), mosses and ferns 3-deoxyanthocyanins, seed plants the 3-hydroxyanthocyanins, and angiosperms a great range of substituted anthocyanins (including 5'-hydroxylation and variation in glycosylation, acylation, and methylation). However, as mentioned earlier, it is now more difficult to speculate on the possible red pigments present in the last common ancestor of land plants. Riccionidin A has been reported from the root cultures of the angiosperm Rhus chinensis (syn. Rhus javanica) (Taniguchi et al., 2000), but it has not been examined whether this is synthesized via an aurone intermediate route.
In addition to the core flavonoid pathway found across most land plants, there are groups of flavonoids prevalent in specific taxonomic groups, such as the isoflavonoids typical of legumes. There are also flavonoid types that occur sporadically, such as aurones that are found in liverworts and some angiosperms. For aurones this may well represent convergent evolution, as even within angiosperms there are alternative biosynthetic mechanisms (Boucherle et al., 2017). New metabolomic technologies combined with additional genome sequences for non-angiosperm species should help clarify the distribution across the land plants of different flavonoid types and the associated biosynthetic genes (Yonekura-Sakakibara et al., 2019).
The Phenylpropanoid Biosynthetic Pathway in Bryophytes
The core steps of the phenylpropanoid pathway through to the first flavonoids (the chalcones) are conserved across land plants (Tohge et al., 2013), including the presence of PAL, C4H, 4CL, CHS, and CHIL gene sequences in hornworts. Sequences relating to some of these genes are present in the genome sequences of charophyte and chlorophyte algae (Labeeuw et al., 2015; de Vries et al., 2017), but without functional assays the conclusions that can be drawn are limited. Most phenylpropanoid pathway enzymes are thought to have evolved from primary metabolism enzymes (Tohge et al., 2013; Yonekura-Sakakibara et al., 2019), and so related sequences might be expected to be present. For PAL, whether it arose during land plant evolution or is an ancestral gene from algae has yet to be resolved (de Vries et al., 2017). It was suggested that PAL was acquired by the land plant ancestor via a horizontal gene transfer event (Emiliani et al., 2009), but genes related to PAL are present in the charophyte Klebsormidium flaccidum and could have been acquired by endosymbiotic gene transfer from cyanobacteria to algal ancestors of land plants (de Vries et al., 2017). C4H, which belongs to the CYP73A sub-family of cytochrome P450 monooxygenases (Cyp450s), shows strong sequence conservation across land plants, including characteristic motifs and residues, but no authentic gene sequences are apparent in chlorophyte genomes (Tohge et al., 2013; Davies et al., 2020). In contrast, sequences with similarity to 4CL do occur in rhodophyte and chlorophyte genomes (Labeeuw et al., 2015; de Vries et al., 2017), suggesting the existence of this enzyme in a shared ancestor of land plants and algae before the ancestral divergence of the red algae (Labeeuw et al., 2015). A further aspect yet to be addressed is the presence in fungi of genes with significant sequence similarity to those of the phenylpropanoid pathway (Bilska et al., 2018; Lu et al., 2019). As the separation of fungi and plants is thought to have occurred during the early stages of eukaryote divergence (Burki, 2014; Burki et al., 2020), it is possible that these may represent cases of convergent evolution.
The type III PKS superfamily that contains CHS is present in all plant genomes examined to date (Pandith et al., 2020). PKS genes are found also in fungi, and some bacteria and algae, and the plant PKS genes contain conserved structural elements with the bacterial PKS genes involved in primary metabolism. Across plants there is a wide variety of PKS enzymes with close sequence similarity to CHS but which either use alternative substrates (such as acridone synthases and pyrone synthases) or catalyze different cyclisation reactions using the same starter molecules (notably, STILBENE SYNTHASE, STS). It is thought that STS has independently evolved from CHS several times in the course of evolution (Yonekura-Sakakibara et al., 2019; Pandith et al., 2020). It is probable that there are many novel PKS activities still to be discovered in plants, including bryophytes. This may include steps in bibenzyl biosynthesis, a group of liverwort phenylpropanoid compounds related to plant defense that includes cannabinoid-like structures (Hussain et al., 2018). The presence of at least 24 PKS genes in the Marchantia genome suggests potential biosynthetic diversity, and at least one gene (Mapoly0014s0122) is closely related to the anther-specific chalcone synthase-like enzymes (ASCLs) involved in the biosynthesis of sporopollenin in angiosperms (Bowman et al., 2017). However, the majority of the annotated MpPKS genes appear to have resulted from a strongly conserved duplication of a CHS/PAL gene pair (Bowman et al., 2017).
The occurrence and function of CHI and CHIL in basal plants is proving to be an interesting question. Liverworts have both types of gene, and knockout chi mutants of Marchantia completely lose production of flavones (Clayton et al., 2018). Thus, in Marchantia, as in angiosperms examined, CHI is an essential in planta activity for flavanone production. However, no gene sequences for CHI have been found in moss or hornwort genome sequences or transcriptomes (Ngaki et al., 2012; Cheng et al., 2018; Berland et al., 2019). Although spontaneous closure to form the C-ring to produce flavanones from chalcones has been shown to occur in vitro, comparative studies on the spontaneous and enzyme catalyzed reactions suggest this is unlikely to be significant in planta (Jez and Noel, 2002). Studies on mutants for chi in Arabidopsis (tt5), carnation (i), and rice (gh1) found that flavonoid biosynthesis was not fully prevented (Stich et al., 1992; Hong et al., 2012; Jiang et al., 2015), suggesting some spontaneous conversion. However, in the case of carnation at least, the residual production of flavanones in the chi mutant has been found to be due to a second, weakly expressed, CHI gene (Miyahara et al., 2018). Thus, how flavonoid biosynthesis occurs in mosses is an open question. CHI and CHIL are thought to be examples of the rare event of catalytic activity arising in a noncatalytic scaffold protein (Kaltenbach et al., 2018). The mechanism of action of CHIL is unclear, and it may have differing activities across land plants, perhaps based on the promotion of activity of different biosynthetic enzymes through protein-protein interaction. In hop (Humulus lupulus), HlCHIL2 enhances the activities of CHS and an aromatic prenyltransferase (HlPT1L) through protein–protein interaction (Ban et al., 2018), and the promotion of flavonol and proanthocyanidin biosynthesis in Arabidopsis is proposed to be through direct interaction of CHIL and CHI (Jiang et al., 2015). In Marchantia, CHIL may interact with CHS or more than one phenylpropanoid pathway enzyme, since the production of both flavones and auronidins in chil mutants is only about 10% of wild-type amounts (Clayton et al., 2018). Thus, one possibility is that in mosses and hornworts CHIL can replace CHI. However, the moss CHIL genes assayed to date do not have CHI activity (Cheng et al., 2018), making this less probable.
Two major hydroxylase groups, the Cyp450s and 2-oxoglutarate dioxygenases (2OGDs, divided into the three classes DOXA, B, and C), contribute several enzymes to the phenylpropanoid pathway of angiosperms. Cyp450s include C4H, FLAVONOID 3'-HYDROXYLASE (F3’H), and FLAVONE SYNTHASE II (FNSII). 2OGDs include the FLAVANONE 2-HYDROXYLASE (F2H), FLAVANONE 3-HYDROXYLASE (F3H), FLAVONOL SYNTHASE (FLS), FLAVONE SYNTHASE I (FNSI), and ANTHOCYANIDIN SYNTHASE/LEUCOANTHOCYANIDIN DIOXYGENASE (ANS). The evolutionary aspects of these gene families with regard to flavonoid biosynthesis were recently reviewed by Yonekura-Sakakibara et al. (2019). C4H is conserved in bryophytes, and the presence of all the other enzymes in liverworts and/or mosses would be expected based on the compounds produced. However, the close similarity of the sequences within the Cyp450 and 2OGD enzyme groups means that assignments based only on sequence similarity to the angiosperm genes should be treated with caution, and conclusive identification of other genes requires functional analysis. FNSI, F2H, and F3H have high sequence similarity and are in the DOXC28 clade and FLS and ANS are close in sequence and in the DOXC47 clade. A review of the two clades, and possible evolutionary timing of the origin of each, is given in Yonekura-Sakakibara et al. (2019).
Based on the occurrence of flavones in liverworts and mosses, it is expected that F3H and FNS activities evolved early in land plants, and two DOXC28 genes have increased transcript abundance during UVB-induced flavone production in Marchantia (Clayton et al., 2018). However, to date, the only functional characterization is for a F2H that may contribute to flavone biosynthesis in the liverwort Plagiochasma appendiculatum (Han et al., 2014). The biosynthesis of flavones illustrates the difficulties of making assumptions about gene function, as a variety of alternative routes to flavone O- and C-glycosides have evolved in angiosperms (Jiang et al., 2016). A further complication is that the 2OGD enzymes (particularly the FLS and ANS) show promiscuous and sometimes overlapping activities when assayed in vitro (reviewed in Martens et al., 2010). Studies with Arabidopsis have shown that these “secondary” activities can also be present in planta, as ANS can contribute to (relatively weak) flavonol biosynthesis in the Arabidopsis fls-1 mutant (Martens et al., 2010).
As yet, it is not clear precisely what phenylpropanoid biosynthetic activities may be present in bryophytes but not found in other plant groups. There are certainly some major pathway branches prevalent in bryophytes that are absent or rare in other groups, such as those for bibenzyls, auronidins, and sphagnorubins. Corresponding evolutionary divergence of specialized metabolic pathways would be expected to underpin the occurrence of the differing compound types. Phylogenetic analysis of the 148 Cyp450, 38 2OGD, and 41 Family-1 UDP-glycosyltransferase (the UGT family containing the “plant secondary product glycosyltransferase” motif) genes of Marchantia found that the majority formed individual clades that also suggested substantial lineage-specific diversification of specialized metabolism (Bowman et al., 2017). Moreover, the emerging transcriptome and genome sequence information from bryophytes is suggesting expanded functionality may have occurred for other classes of enzymes involved in phenylpropanoid biosynthesis. In the next section we examine two specialized metabolism gene families that show unexpectedly large gene family sizes in the Marchantia genome: those for POLYPHENOL OXIDASE (PPO) and DIRIGENT (DIR) proteins.
Liverworts May Have Expanded Functional Roles in Specialized Metabolism for Polyphenol Oxidase and Dirigent Proteins
PPO genes are found throughout land plants, as well as in bacteria, fungi, and animals, but are absent from algae. PPOs are type-III-copper proteins and the name PPO covers two major enzyme types: tyrosinases, which hydroxylate para-substituted monophenols to ortho-diphenols (monophenolase activity) and use molecular oxygen to oxidize ortho-diphenols to ortho-quinones (diphenolase activity); and the catechol oxidases, which have only the diphenolase activity. However, it has been recently proposed that monophenolase activity could be a widespread feature of PPOs, but that the activity has remained cryptic because activity assays usually use tyrosine rather than the natural substrates, which are often not known (Molitor et al., 2016). PPOs are commonly thought of as plant defense enzymes that oxidize and/or polymerize a range of phenolic substrates with which they come into contact during cell disruption, resulting in the familiar browning reactions following tissue damage, for example in cut apples or potatoes. However, in addition to these general activities, some PPOs can conduct cross-linking reactions in biosynthetic pathways, such as latex formation; and new specific roles for PPOs have emerged in recent years (Figure 6). The published PPO gene family size in plants varies from zero (e.g., Arabidopsis) to 13 in Physcomitrella patens (Tran et al., 2012). Several angiosperm species examined have only a single PPO gene, but 11 genes have been found in genome sequences of Glycine max (the legume soybean), Populus trichocarpa (poplar), and Selaginella moellendorffii (a lycophyte) (Tran et al., 2012). However, we found a much larger PPO gene family in the Marchantia genome: there are 64 candidate PPO genes (including gene fragments and unresolved gene models). Excluding those having partial gene models, 46 of the 64 PPO genes were represented in the RNA-seq data of Berland et al. (2019) and so are actively transcribed. Given the relatively small number of total gene models in the draft genome sequence of Marchantia, this represents a significant gene family, larger than the annotated 2OGD and UGT families.
Figure 6 Examples of plant POLYPHENOL OXIDASES with specific activities in specialized metabolism. The biosynthetic activities of three PPOs are shown: (+) Larreatricin hydroxylase, which is part of lignan biosynthesis in the creosote bush (Larrea tridentata), and the AURONE SYNTHASE from Coreopsis grandiflora (CgAUS) and Antirrhinum majus (AmAUS).
Plant PPOs characterized to date are produced in a latent state as proteins of about 64–68 kDa. Besides the N-terminal targeting peptide (usually for plastid localization), PPOs contain a catalytically active domain of about 40 kDa, and a C-terminal domain of about 19 kDa that shields the active site and is later cleaved off to release the active protein. The C-terminal domain is ubiquitous in plant PPOs examined to date. Based on predicted amino acid sequences, PPOs with this typical structure are found in Marchantia; however, there are also members of the PPO family that lack this C-terminal domain (Figure 7), including the auronidin-related Mapoly0021s0041. These “short” type PPOs have also been found in fungi and bacteria (Huber et al., 1985; Shuster and Fishman, 2009; Gasparetti et al., 2010). Only a few of this short type, from the bacteria Streptomyces and Bacillus, have been extensively studied. The Streptomyces PPO is thought to be initially in an inactive form that is bound with a “caddie” protein. The caddie protein subsequently transfers copper to the PPO and disassociates to release an active PPO (Chen et al., 1992; Matoba et al., 2006). In contrast, the PPO from Bacillus does not need a caddie protein (Sendovski et al., 2011).
Figure 7 Amino acid sequence features of polyphenol oxidases (PPOs) in Marchantia polymorpha. (A) Number of expressed PPO gene models found and their PPO type. RNA-seq data from Berland et al. (2019) were used to check for expression. The top structure (not to scale) is standard for plant PPOs. Transit/signal peptide (P), copper binding domain TYR : PFAM 00264 (CuA/CuB), PPO1_DWL: PFAM 12142 domain (DWL), tyrosine motif (YxY), PPO1_KFDV : PFAM12143 domain (KFDV). (B) Weblogo display of conserved residues of DWL and KFDV domains based on all gene models having an intact region (generated using https://weblogo.berkeley.edu/logo.cgi). DWL, tyrosine, and KFDV core motifs are underlined in black. Residues for the regions identified by Tran et al. (2012) with core motifs of EEEVLV (left) and EFAGSF (right) are underlined in blue.
The first PPO found to have an unexpected role in plant specialized metabolism was the AUREUSIDIN SYNTHASE (AUS) that converts chalcone 4'-O-glucosides to aurone 6-O-glucosides in Antirrhinum majus (Nakayama et al., 2000; Davies et al., 2006; Ono et al., 2006; Elumalai and Liu, 2011). AmAUS differed from previously characterized PPOs in three important aspects: it was vacuole localized (Ono et al., 2006), it was a glycoprotein, and it lacked activity against common PPO substrates such as tyrosine or 3,4-dihydroxyl L-phenylalanine (L-DOPA). AmAUS conducts oxygenation of the B-ring of the chalcone, which is followed by cyclization into the aurone (Nakayama et al., 2000) (Figure 6). Although it can use chalcone aglycones in vitro, in planta aurone production in A. majus requires the activity of the CHALCONE 4'-O-GLUCOSYLTRANSFERASE (C4’GT) since only the glucoside is transported into the vacuole (Ono et al., 2006; Bradley et al., 2017). Subsequently, PPOs that form aurones were identified in other species, with the AURONE SYNTHASE of Coreopsis grandiflora that makes 4-deoxyaurones being studied in detail (Molitor et al., 2015; Molitor et al., 2016). In contrast to AmAUS, the CgAUS has the N-terminal chloroplast transit peptide and thylakoid transfer domain characteristic of plastid-localized PPOs involved in browning reactions, and uses chalcone aglycones to make aurone aglycones (Kaintz et al., 2014; Molitor et al., 2015). A PPO (Mapoly0021s0041) is strongly up-regulated by MpMYB14 in association with auronidin production in Marchantia (Berland et al., 2019). Loss-of-function mapoly0021s0041 mutants have greatly reduced amounts of auronidin, suggesting it too may encode an aurone biosynthetic activity, or is involved in later steps of auronidin biosynthesis and/or polymerization. Aurones have been found across land plant groups, but with sporadic occurrence. This suggests that their biosynthesis may have arisen independently on a number of occasions. Although both the aurone biosynthetic enzymes characterized to date are PPOs, the sequences are phylogenetically distinct and have differing activities and sub-cellular localization. Additionally, the biosynthesis of the aurone hispidol in Medicago truncatula may be conducted by a peroxidase rather than a PPO (Farag et al., 2009). Besides AUS, PPOs have been implicated in tyrosine or phenylpropanoid biosynthetic pathways of walnut (Araji et al., 2014) and creosote bush (Larrea tridentata) (Cho et al., 2003).
Forty of the Marchantia PPO genes (66%) occur as tandem repeats or small gene clusters (local tandemly arrayed genes; TAGs). Although this figure may be either an under- or over-estimate as it is based on the initial scaffold assembly of the genome (Bowman et al., 2017), it is nevertheless a much higher value than the overall percentage of Marchantia TAGs estimated on the same basis, which at 5.9% is near the lower end of the range observed in flowering plants (Bowman et al., 2017). The TAG percentage is also relatively high for some of the other characterized specialized metabolite gene families of Marchantia. For example, there are 18 occurrences of neighboring PAL and/or CHS genes. TAGs are notable in some angiosperm species that have prominent specialized metabolic characteristics—such as the terpenoid pathways of the tree species Eucalyptus grandis and teak (Tectona grandis). Teak has at least 14 TAGs for the terpene synthase gene family (Zhao et al., 2019). E. grandis has the largest number of genes in tandem repeats reported among sequenced plant genomes, at 34% of total genes (Myburg et al., 2014). For the Marchantia phenylpropanoid biosynthetic pathway, 10 multigene families have expanded, mostly through tandem duplication, to result in a total of 174 genes. In angiosperms, gene diversification is a result of a combination of local duplication events and whole genome duplications, but it is probable that no whole-genome duplication events have occurred during liverwort evolution (Bowman et al., 2017). Therefore, although there is no overall increase in the frequency of TAGs in liverworts (at least for Marchantia), local gene duplication events are likely to have been a common mechanism for generating gene neofunctionalization in specialized metabolism. Whether this is typical of other liverworts requires the completion of further genome sequences. However, BLAST analysis of the Lunularia cruciata transcriptome (www.polebio.lrsv.ups-tlse.fr/Luc_v1/Luc_v1.fa) identified more than 20 sequences with the conserved features of PPOs (data not shown), suggesting a large gene family in this species also.
Dirigent proteins (DIR) are small (~16–18 kDa) cell wall-localized proteins that may control the regio- and stereospecific outcome of phenoxy radical coupling in lignin and lignan polymerization reactions (Davin et al., 1997; Gang et al., 1999). The polymerization reactions also require the activity of laccase or peroxidase to produce electron oxidative capacity to generate the phenoxy radical. In vascular plants, lignins are complex, amorphous heteropolymers involved in wall strengthening and pathogen resistance, with species-specific composition produced by polymerization of coniferyl, sinapyl, and p-coumaryl alcohols. A role for DIRs in directing the reactions leading to the formation of lignin has been proposed but not definitively established, although there is strong genetic evidence in support of some specific cases (Hosmani et al., 2013).
The role of DIRs in determining stereospecificity has been best described in the formation of lignans, a class of 8-8’ linked C6C3 phenylpropanoid dimers involved in pathogen resistance, for example in the production of (+)- or (−)-pinoresinol compounds in flax (Linum usitatissimum) and pea (Pisum sativum) (Kim et al., 2015; Corbin et al., 2018). The X-ray crystal structure of PsDRR206 involved in (+)-pinoresinol formation suggested that the active protein had a trimeric structure (Kim et al., 2015). Recent work has suggested that at least some DIRs may do more than the hypothesized positioning of phenoxy radicals prior to coupling, and may themselves possess enzymatic activity. The crystal structure of Arabidopsis AtDIR6 identified potentially catalytic residues including aspartic acids that were essential for activity, and it was proposed that this protein catalyzed the cyclization of the bis-quinone methide intermediate during (+)- or (–)-pinoresinol formation (Gasper et al., 2016). Also, a recombinantly expressed DIR from Glycyrrhiza echinata was found to possess isoflavanol dehydratase activity and carry out the final ring-closure step of the biosynthesis of the anti-microbial phytoalexin (–)-pterocarpan (Uchida et al., 2017).
DIR gene families can be quite large, with 26 genes in Arabidopsis (Paniagua et al., 2017) and 44 genes in flax (of which seven appeared to be gene fragments or result from chromosomal rearrangements; Corbin et al., 2018). Of the 37 genes with classical DIR structure in flax, 15 paralogous gene pairs were identified. Kubo et al. (2018) identified 52 dirigent-like predicted proteins in the Marchantia genome sequence. Our analysis for this article found that at least 35 of these occur as TAGs. However, the deduced protein sequences of the family members are diverse, and the functionality of the proteins has not yet been established. Our BLAST analysis of the Marchantia genome and transcript resources with the 24 annotated Arabidopsis DIR genes gave us 60 initial candidate gene models, with strong evidence of some very recent gene duplications giving groups of adjacent genes with highly similar or identical deduced amino acid sequences.
Phenylpropanoid biosynthesis and lignification are common plant responses to biotic and abiotic stress (Zhao and Dixon, 2014; Paniagua et al., 2017) and consequently, as a component of lignification, DIRs have been implicated in responses to pathogen and drought stress (e.g., Thamil Arasan et al., 2013; Paniagua et al., 2017). In P. patens, fungal infection resulted in increased incorporation of phenolic compounds into the wall and up-regulation of a DIR gene (Reboledo et al., 2015). In Marchantia, abiotic stresses such as UVB irradiation, N deficiency and salinity (Albert et al., 2018; Kubo et al., 2018), and pathogen attack (Carella et al., 2019) increased the expression of MpMYB14. MpMYB14 promotes auronidin production and up-regulates transcript abundance for at least three DIR genes (Mapoly0006s0216, Mapoly0006s0217, Mapoly0078s0058) (Albert et al., 2018; Kubo et al., 2018). The deduced protein products of these genes possess predicted signal peptides (SignalP 5.0, http://www.cbs.dtu.dk/services/SignalP/), indicative of secretion to the vacuole or, extracellularly, to the cell wall. However, prediction of subcellular localization using WoLFPSORT (https://wolfpsort.hgc.jp) indicated with low confidence different compartments for the three predicted protein products of the genes: extracellular (Mapoly0006s0217), vacuolar (Mapoly0006s0216), and cytoplasmic (Mapoly0078s0058). Intracellular coupling of monolignol radicals has been described in Arabidopsis (Dima et al., 2015). There have been no studies on DIR genes in hornworts. However, as lignans are prominent specialized metabolites of hornworts, and DIRs have roles in lignan biosynthesis in angiosperms, this could be a worthwhile area to investigate.
Evolution of the Transcriptional Regulation of the Phenylpropanoid Pathway
In angiosperms and gymnosperms, the key regulatory complex consists of R2R3MYB and bHLH TFs joined with a WD-Repeat (WDR) protein, a composition of proteins known as an “MYB-bHLH-WD repeat (MBW)” complex. The MBW complex that activates anthocyanin and proanthocyanidin production contains R2R3MYB proteins from sub-group (SG) 5 or 6, and commonly promotes transcription of the biosynthetic genes throughout the pathway. The action of the MBW complex is modified by a WRKY class activator TF and a series of proteins with repressor actions (Lloyd et al., 2017). In particular, R2R3MYBs from SG4 can join an activating MBW complex and turn it into one that represses target gene transcription, and R3MYBs can bind the bHLH to prevent it from forming the MBW complex, thus competitively inhibiting activation (Albert et al., 2014). The SG4 R2R3MYBs are characterized by the presence of an ethylene response factor (ERF)-associated amphiphilic repression (EAR) motif (LxLxL or DLNxxP) or a TLLLFR motif in the C terminus that mediates transcriptional repression (Chen et al., 2019a; Chen et al., 2019b; Ma and Constabel, 2019). For the activation of the flavonol and flavone branches, a SG7 R2R3MYB acts without being part of the complex. There is also regulation upstream of the flavonoid pathway, as HY5 activates the production of the SG7 R2R3MYB. Additionally, in Arabidopsis it has been shown that HY5 directly activates transcription of some flavonoid biosynthetic genes, such as CHS. The conservation of HY5 function in the UVB responses of both bryophytes and angiosperms was mentioned earlier, although its target gene set has yet to be resolved.
The expansion of TF families during evolution has been a driver of diversity in land plants, as a consequence of multicellularity and increased organismal complexity and/or for coping with the increased stress of a sessile land-based lifestyle. The MYB gene family is one of the largest TF families in plants, with Arabidopsis having 137 R2R3MYB genes (Feller et al., 2011). This includes one SG5, four SG6, and three SG7 genes in Arabidopsis for proanthocyanidin, anthocyanin, and flavonol production, respectively. The presence of small gene families for sub-groups regulating specialized metabolic pathways is common for angiosperms, and has enabled sub-functionalization and diversification of flavonoid temporal and spatial regulation in flowers, seeds, and vegetative tissues. The bHLH and WDR components are less specific in their regulatory targets, and can regulate other characters as well as flavonoid biosynthesis, such as epidermal cell differentiation.
The great majority of information on the transcriptional regulation of specialized metabolite pathways is available from studies on angiosperms, with only a small number of studies on gymnosperm, fern, or bryophyte species. Identifying the genetic components for flavonoid pathway regulation in these other plant groups will help establish a model for how regulation of specialized metabolism may have changed during evolution. For bryophytes, notable questions relating to flavonoid pathway regulation include: are R2R3MYB and bHLH genes the key direct activators? If MYBs are the direct activators, which SGs are present in bryophytes and do small gene families occur for each SG? Does a MBW complex form in bryophytes? Do repressor TFs modify pathway regulation? Characterizing these aspects in species such as M. polymorpha, P. patens, and the lycophyte model S. moellendorffii should indicate which aspects of flavonoid regulation are conserved across land plants, and thus may have been present in the early land plant ancestor, and which aspects may have arisen as part of evolutionary diversification of the different land plant groups.
Compared with angiosperms, the characterized bryophytes and lycophytes have small TF families. There are only 22, 49, and 62 R2R3MYBs in the genomes of M. polymorpha, P. patens, and S. moellendorffii, respectively (Feller et al., 2011; Bowman et al., 2017). For Marchantia genes, a phylogenetic comparison of this gene family shows that MpMYB02 and MpMYB14 fall basal to a clade that contains all the phenylpropanoid-related R2R3MYB genes of Arabidopsis (SGs 4, 5, 6, 7, 15, and 44) (Bowman et al., 2017). Concluding whether these correspond to descendants of the flavonoid regulatory R2R3MYBs of the ancestral land plant requires further study, although both MpMYB02 and MpMYB14 activate phenylpropanoid biosynthetic genes. MpMYB02 is required for production of bibenzyls (Kubo et al., 2018) while MpMYB14 is essential for auronidin production and promotes the production of flavone O-glycosides (Albert et al., 2018; Clayton et al., 2018; Kubo et al., 2018). The profiles of transcripts up-regulated by MpMYB02 and MpMYB14 include DIR genes (Albert et al., 2018; Kubo et al., 2018; Berland et al., 2019). For MpMYB14, this includes the three DIR genes discussed earlier as well as other DIR genes that have been shown to be direct targets (Kubo et al., 2018). Co-expression analysis in flax found that MYB TFs were up-regulated along with DIR genes during secondary wall biosynthesis (Corbin et al., 2018), suggesting that MYB proteins could control DIR expression in both angiosperms and bryophytes.
MpMYB14 must act redundantly with other uncharacterized TFs for flavone production, as Mpmyb14 mutants still show the induction of flavones in response to UVB (Clayton et al., 2018), nutrient stress, or high-irradiance white light (Albert et al., 2018). Flavone production is reduced in Mphy5 mutants, so it is possible that HY5 is a direct activator of flavonoid biosynthetic genes as in Arabidopsis, but there may also be HY5-independent activation pathways for flavone production (Kondou et al., 2019). Analysis of changes in transcriptomes in response to UVB treatment does not present any alternative R2R3MYB candidate for flavone regulation (Clayton et al., 2018). Thus, Marchantia may lack the equivalent of the angiosperm SG7 activators of flavonol and flavone biosynthesis.
The PabHLH gene of the liverwort P. appendiculatum is a probable activator of bibenzyl biosynthesis (Wu et al., 2018). Over-expression of PabHLH in P. appendiculatum increased bibenzyl concentration and up-regulated transcript abundance from known phenylpropanoid biosynthetic genes (PAL, 4CL) and candidate bibenzyl biosynthetic genes, whereas RNA interference-induced suppression down-regulated the same genes and reduced bibenzyl accumulation. Phylogenetically, PabHLH falls within clades containing the flavonoid MBW bHLH sequences of angiosperms (within bHLH subgroup IIIf), suggesting it may be homologous to them. In Marchantia, MpbHLH12 is the gene with the highest sequence identity to PabHLH and the flavonoid-related bHLHs of angiosperms, and transcriptomic analysis of MpBHLH12 overexpression transgenics suggests it may also be involved in flavonoid regulation (Arai et al., 2019). However, although R2R3MYB and bHLH genes do regulate flavonoid biosynthesis in liverworts, and there are conserved WDR sequences in the genome (Bowman et al., 2017), there is no answer yet on whether the MBW complex exists in bryophytes. A flavonoid-related MBW complex has been characterized in the gymnosperm Norway spruce (Picea abies) (Nemesio-Gorriz et al., 2017), supporting an origin for the MBW complex in the plant lineage prior to the last common ancestor of gymnosperms and angiosperms, around 350–300 MYA. However, although the conserved amino acid motif ([D/E]Lx2[R/K]x3Lx6Lx3R) identified as necessary for R2R3MYB proteins to bind the bHLH partners (Zimmermann et al., 2004) is present in the S. moellendorffii sequence SmXP002978781, it is lacking in bryophyte R2R3MYBs studied to date. The closest matches in P. patens (PpXP001752936) and Marchantia (MpMYB02 and MpMYB14) lack one and two deduced amino acid residues, respectively.
Whether bryophytes possess MYB genes with a repressive action in phenylpropanoid regulation, either the R2R3MYB active repressors that form part of the MBW complex or the R3MYBs that are thought to “compete” for the bHLH proteins, is also an open question. We were unable to identify (known) repression motifs in any of the Marchantia R2R3MYB sequences. Both P. patens and S. moellendorffii have R2R3MYB genes with putative EAR motif sequences (LxLxL), but the possible function of these in regulating phenylpropanoid biosynthesis has not been examined. Analysis of the auxin signaling pathway of Marchantia has identified an orthologue of TOPLESS, which in angiosperms interacts with the EAR motif to mediate transcriptional repression (Flores-Sandoval et al., 2015). The Marchantia genome contains an expanded R3MYB gene family (Bowman et al., 2017), but no analysis of these with regard to flavonoid biosynthesis has been published.
In summary, based on the evidence from Marchantia, it seems probable that the ancestral R2R3MYB regulators of phenylpropanoid metabolism were activators acting outside of an MBW complex. R2R3MYB-repressive TFs and the MBW complex probably evolved after the last common ancestor of liverworts and gymnosperms/angiosperms. As the flavone pathway probably evolved prior to anthocyanin biosynthesis, it could be expected that R2R3MYBs most similar to SG7 might be the ancestral type. However, the specific flavone activators of Marchantia have yet to be identified. MpMYB02 and MpMYB14 may correspond to the ancestral phenylpropanoid pathway activators, and like SG7 probably act outside the MBW complex, but it is difficult to state which is the most closely related SG because of the extent of sequence divergence, with no conservation of sequence outside the MYB domains themselves. Furthermore, additional data are required from other bryophyte species, as the evolutionary path to Marchantia will have resulted in extensive genetic changes and the loss of characters that were present in the last common ancestor.
Concluding Comments
The commonality of phenylpropanoid biosynthetic genes between bryophytes and angiosperms, and the conserved functions of flavonoids in assisting in tolerance to stresses such as UVB and pathogen attack, support the proposal that the pathway arose before the last common ancestor of these land plant groups, relatively early during the process of land colonization. The exception to this is the hornworts, which lack flavonoids. Unless the divergence of hornworts occurred before the pathway arose, the hornwort ancestor must have acquired mutations that caused loss of the biosynthetic or regulatory capacity. This may be analogous to the loss of anthocyanin biosynthesis in some lineages of the Caryophyllales, where they are replaced by betalains. As the main red pigments of angiosperms (soluble anthocyanins) and bryophytes (cell wall-bound auronidins and sphagnorubins) differ in structure and cellular properties, it is difficult to suggest what the original common ancestor may have possessed with regard to red pigments. Establishing which components of anthocyanin biosynthesis are present or lacking in bryophytes may help in this regard.
The diversification of both specialized metabolite biosynthesis and the transcription factors that regulate the pathways are thought to be important contributors to the evolution of plants to occupy the varied ecological niches offered on land (Pichersky and Gang, 2000). To date, much of our understanding of the genetic basis of the diversification process has been based on studies of flowering plants. However, the completion of the first genome sequences for a moss (P. patens) and liverwort (M. polymorpha) has started to reveal the details of the specialized metabolite gene families, such as for phenylpropanoid biosynthesis. Notably, the TF families thought to regulate the phenylpropanoid pathway are much smaller in Marchantia than in flowering plants. However, there are relatively large Marchantia gene families for enzymes that are often involved in specialized metabolism, such as the Cyp450, 2OGD, and UGT families. Moreover, Marchantia has large PPO and DIR gene families compared to angiosperms, suggesting these enzyme groups may make a greater contribution than previously anticipated to phenylpropanoid and other specialized metabolite biosynthesis in the liverworts. Thus, in liverworts some of the gene families involved in the biosynthesis of specialized metabolites appear to have undergone more gene duplication (allowing consequent sub- and neofunctionalization for particular family members) than the TFs that regulate the same pathways. Expansion of the regulatory TF families through duplication and sub/neofunctionalization is seen in the angiosperms, probably reflecting increased organismal complexity.
Author Contributions
KD, RJ, YZ, NA, DB, BJ, JB and KS reviewed literature, formulated ideas and wrote the manuscript. KD and KS prepared the figures. KD, RJ, YZ, DB and KS conducted bioinformatic analysis.
Funding
Financial support was provided by the Marsden Fund of New Zealand Grant PAF1701.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2020.00007/full#supplementary-material
References
Agati, G., Tattini, M. (2010). Multiple functional roles of flavonoids in photoprotection. New Phytol. 186, 786–793. doi: 10.1111/j.1469-8137.2010.03269.x
Agati, G., Azzarello, E., Pollastri, S., Tattini, M. (2012). Flavonoids as antioxidants in plants: Location and functional significance. Plant Sci. 196, 67–76. doi: 10.1016/j.plantsci.2012.07.014
Agregán, R., Munekata, P. E. S., Franco, D., Dominguez, R., Carballo, J., Lorenzo, J. M. (2017). Phenolic compounds from three brown seaweed species using LC-DAD–ESI-MS/MS. Food Res. Int. 99, 979–985. doi: 10.1016/j.foodres.2017.03.043
Aigner, S., Remias, D., Karsten, U., Holzinger, A. (2013). Unusual phenolic compounds contribute to ecophysiological performance in the purple-colored green alga Zygogonium ericetorum (zygnematophyceae, streptophyta) from a high-alpine habitat. J. Phycol. 49, 648–660. doi: 10.1111/jpy.12075
Albert, N. W., Davies, K. M., Lewis, D. H., Zhang, H., Montefiori, M., Brendolise, C., et al. (2014). A conserved network of transcriptional activators and repressors regulates anthocyanin pigmentation in eudicots. Plant Cell 26, 962–980. doi: 10.1105/tpc.113.122069
Albert, N. W., Thirmawithana, A. H., Mcghie, T. K., Clayton, W., Deroles, S. C., Schwinn, K. E., et al. (2018). Genetic analysis of the liverwort Marchantia polymorpha reveals that R2R3MYB activation of flavonoid production in response to abiotic stress is an ancient character in land plants. New Phytol. 218, 554–566. doi: 10.1111/nph.15002
Allorent, G., Lefebvre-Legendre, L., Chappuis, R., Kuntz, M., Truong, T. B., Niyogi, K. K., et al. (2016). UV-B photoreceptor-mediated protection of the photosynthetic machinery in Chlamydomonas reinhardtii. Proc. Natl. Acad. Sci. U. S. A. 113, 14864–14869. doi: 10.1073/pnas.1607695114
Andersen, O. M., Markham, K. R. (Eds.) (2006). The Flavonoids: Chemistry, Biochemistry and Applications. (Boca Raton, USA: CRC Press).
Arai, H., Yanagiura, K., Toyama, Y., Morohashi, K. (2019). Genome-wide analysis of MpBHLH12, a IIIf basic helix-loop-helix transcription factor of Marchantia polymorpha. J. Plant Res. 132, 197–209. doi: 10.1007/s10265-019-01095-w
Araji, S., Grammer, T. A., Gertzen, R., Anderson, S. D., Mikulic-Petkovsek, M., Veberic, R., et al. (2014). Novel roles for the polyphenol oxidase enzyme in secondary metabolism and the regulation of cell death in walnut. Plant Physiol. 164, 1191–1203. doi: 10.1104/pp.113.228593
Archetti, M. (2009). Classification of hypotheses for the evolution of autumn colours. Oikos 118, 328–333. doi: 10.1111/j.1600-0706.2008.17164.x
Asakawa, Y., Ludwiczuk, A., Nagashima, F. (2013). “Chemical Constituents of Bryophytes; Bio- and Chemical Diversity, Biological Activity, and Chemosystematics,” in Progress in the Chemistry of Organic Natural Products. Eds. Kinghorn, A. D., Falk, H., Kobayashi, J. (Vienna: Springer), 1–796. doi: 10.1007/978-3-7091-1084-3
Asakawa, Y. (2017). “Polyphenols in Bryophytes,” in Recent Advances in Polyphenol Research. Eds. Yoshida, K., Cheynier, V., Quideau, S., Wiley-Blackwell: Oxford, UK. 36–66. doi: 10.1002/9781118883303.ch2
Ban, Z., Qin, H., Mitchell, A. J., Liu, B., Zhang, F., Weng, J.-K., et al. (2018). Noncatalytic chalcone isomerase-fold proteins in Humulus lupulus are auxiliary components in prenylated flavonoid biosynthesis. Proc. Natl. Acad. Sci. U. S. A. 115, E5223–E5232. doi: 10.1073/pnas.1802223115
Ben Saad, H., Gargouri, M., Kallel, F., Chaabene, R., Boudawara, T., Jamoussi, K., et al. (2017). Flavonoid compounds from the red marine alga Alsidium corallinum protect against potassium bromate-induced nephrotoxicity in adult mice. Environ. Toxicol. 32, 1475–1486. doi: 10.1002/tox.22368
Berim, A. (2016). Methoxylated flavones: occurrence, importance, biosynthesis. Phytochem. Rev. 15, 363–390. doi: 10.1007/s11101-015-9426-0
Berland, H., Albert, N. W., Stavland, A., Jordheim, M., McGhie, T. M., Zhou, Y., et al. (2019). Auronidins are a previously unreported class of flavonoid pigments that challenges when anthocyanin biosynthesis evolved in plants. Proc. Natl. Acad. Sci. U. S. A. 116, 20232–20239. doi: 10.1073/pnas.1912741116
Bilska, K., Stuper-Szablewska, K., Kulik, T., Buśko, M., Załuski, D., Jurczak, S., et al. (2018). Changes in phenylpropanoid and trichothecene production by Fusarium culmorum and F. graminearum Sensu Stricto via exposure to flavonoids. Toxins 10, 110. doi: 10.3390/toxins10030110
Bonnett, S. A. F., Ostle, N., Freeman, C. (2010). Short-term effect of deep shade and enhanced nitrogen supply on Sphagnum capillifolium morphophysiology. Plant Ecol. 207, 347–358. doi: 10.1007/s11258-009-9678-0
Boucherle, B., Peuchmaur, M., Boumendjel, A., Haudecoeur, R. (2017). Occurrences, biosynthesis and properties of aurones as high-end evolutionary products. Phytochemistry 142, 92–111. doi: 10.1016/j.phytochem.2017.06.017
Bowman, J. L., Araki, T., Kohchi, T. (2016). Marchantia: Past, Present and Future. Plant Cell Physiol. 57, 205–209. doi: 10.1093/pcp/pcw023
Bowman, J. L., Kohchi, T., Yamato, K. T., Jenkins, J., Shu, S., Ishizaki, K., et al. (2017). Insights into land plant evolution garnered from the Marchantia polymorpha genome. Cell 171, 287–304. doi: 10.1016/j.cell.2017.09.030
Bradley, D., Xu, P., Mohorianu, I.-I., Whibley, A., Field, D., Tavares, H., et al. (2017). Evolution of flower color pattern through selection on regulatory small RNAs. Science 358, 925–928. doi: 10.1126/science.aao3526
Breuil-Sée, A. (1993). Recorded desiccation-survival times in bryophytes. J. Bryol. 17, 679–680. doi: 10.1179/jbr.1993.17.4.679
Brown, D. E., Rashotte, A. M., Murphy, A. S., Normanly, J., Tague, B. W., Peer, W. A., et al. (2001). Flavonoids act as negative regulators of auxin transport in vivo in arabidopsis. Plant Physiol. 126, 524–535. doi: 10.1104/pp.126.2.524
Brunetti, C., Fini, A., Sebastiani, F., Gori, A., Tattini, M. (2018). Modulation of phytohormone signaling: a primary function of flavonoids in plant–environment interactions. Front. Plant Sci. 9, 1042. doi: 10.3389/fpls.2018.01042
Brunetti, C., Sebastiani, F., Tattini, M. (2019). Review: ABA, flavonols, and the evolvability of land plants. Plant Sci. 280, 448–454. doi: 10.1016/j.plantsci.2018.12.010
Buer, C. S., Muday, G. K. (2004). The transparent testa4 mutation prevents flavonoid synthesis and alters auxin transport and the responses of Arabidopsis roots to gravity and light. Plant Cell 16, 1191–1205. doi: 10.1105/tpc.020313
Burbulis, I. E., Iacobucci, M., Shirley, B. W. (1996). A null mutation in the first enzyme of flavonoid biosynthesis does not affect male fertility in Arabidopsis. Plant Cell 8, 1013–1025. doi: 10.1105/tpc.8.6.1013
Burki, F., Roger, A. J., Brown, M. W., Simpson, A. G. B. (2020). The new tree of eukaryotes. Trends Ecol. Evol. in press. 43–55.
Burki, F. (2014). The eukaryotic tree of life from a global phylogenomic perspective. Cold Spring Harb. Perspect. Biol. 6, a016147. doi: 10.1101/cshperspect.a016147
Campbell, D. H. (1896). The development of Geothallus tuberosus, Campbell. Ann. Bot. 10, 489–510. doi: 10.1093/oxfordjournals.aob.a088624
Carella, P., Gogleva, A., Hoey, D. J., Bridgen, A. J., Stolze, S. C., Nakagami, H., et al. (2019). Conserved biochemical defenses underpin host responses to oomycete infection in an early divergent land plant lineage. Curr. Biol. 29, 1–13. doi: 10.1016/j.cub.2019.05.078
Cesares Gil, D. A. (1902). Algunas observaciones sobre la coloracion rojiza de ciertas hepáticas. Bol. Soc Española Hist. Nat. 2, 217–220.
Chattopadhyay, A., Erland, L. A. E., Jones, A. M. P., Saxena, P. K. (2018). Indoleamines and phenylpropanoids modify development in the bryophyte Plagiomnium cuspidatum (Hedw.) T.J. Kop. In Vitro Cell. Dev. Biol. Plant 54, 454–464. doi: 10.1007/s11627-018-9904-3
Chen, L.-Y., Leu, W.-M., Wang, K.-T., Lee, Y.-H. W. (1992). Copper transfer and activation of the Streptomyces apotyrosinase are mediated through a complex formation between apotyrosinase and its trans-activator MelC1. J. Biol. Chem. 267, 20100–20107.
Chen, C., Zhang, K., Khurshid, M., Li, J., He, M., Georgiev, M. I., et al. (2019a). MYB transcription repressors regulate plant secondary metabolism. Crit. Rev. Plant Sci. 38, 159–170. doi: 10.1080/07352689.2019.1632542
Chen, L., Hu, B., Qin, Y., Hu, G., Zhoa, J. (2019b). Advance of the negative regulation of anthocyanin biosynthesis by MYB transcription factors. Plant Physiol. Biochem. 136, 178–187. doi: 10.1016/j.plaphy.2019.01.024
Cheng, A.-X., Zhang, X., Han, X.-J., Zhang, Y.-Y., Gao, S., Liu, C.-J., et al. (2018). Identification of chalcone isomerase in the basal land plants reveals an ancient evolution of enzymatic cyclization activity for synthesis of flavonoids. New Phytol. 217, 909–924. doi: 10.1111/nph.14852
Cheynier, V., Comte, G., Davies, K. M., Lattanzio, V., Martens, S. (2013). Plant phenolics: Recent advances on their biosynthesis, genetics, and ecophysiology. Plant Physiol. Biochem. 72, 1–20. doi: 10.1016/j.plaphy.2013.05.009
Cho, M. H., Moinuddin, S. G., Helms, G. L., Hishiyama, S., Eichinger, D., Davin, L. B., et al. (2003). (+)-Larreatricin hydroxylase, an enantio-specific polyphenol oxidase from the creosote bush (Larrea tridentata). Proc. Natl. Acad. Sci. U. S. A 100, 10641–10646. doi: 10.1073/pnas.1934562100
Clarke, L. J., Robinson, S. A. (2008). Cell wall-bound ultraviolet-screening compounds explain the high ultraviolet tolerance of the Antarctic moss, Ceratodon purpureus. New Phytol. 179, 776–783. doi: 10.1111/j.1469-8137.2008.02499.x
Clayton, W. A., Albert, N. W., Thrimawithana, A. H., Mcghie, T. K., Deroles, S. C., Schwinn, K. E., et al. (2018). UVR8-mediated induction of flavonoid biosynthesis for UVB tolerance is conserved between the liverwort Marchantia polymorpha and flowering plants. Plant J. 96, 503–517. doi: 10.1111/tpj.14044
Cockell, M. M., Knowland, J. (1999). Ultraviolet radiation screening compounds. Biol. Rev. 74, 311–345. doi: 10.1111/j.1469-185X.1999.tb00189.x
Cohen, M. F., Meziane, T., Tsuchiya, M., Yamasaki, H. (2002a). Feeding deterrence of Azolla in relation to deoxyanthocyanin and fatty acid composition. Aquat. Bot. 74, 181–187. doi: 10.1016/S0304-3770(02)00077-3
Cohen, M. F., Sakihama, Y., Takagi, Y. C., Ichiba, T., Yamasaki, H. (2002b). Synergistic effect of deoxyanthocyanins from symbiotic fern Azolla spp. on hrmA gene induction in the cyanobacterium Nostoc punctiforme. Mol. Plant Microbe Interact. 15, 875–882. doi: 10.1094/MPMI.2002.15.9.875
Corbin, C., Drouet, S., Markulin, L., Auguin, D., Laine, E., Davin, L. B., et al. (2018). A genome-wide analysis of the flax (Linum usitatissimum L.) dirigent protein family: from gene identification and evolution to differential regulation. Plant Mol. Biol. 97, 73–101. doi: 10.1007/s11103-018-0725-x
Dare, A. P., Hellens, R. P. (2013). RNA interference silencing of CHS greatly alters the growth pattern of apple (Malus x domestica). Plant Signal. Behav. 8, e25033. doi: 10.4161/psb.25033
Davies, K. M., Marshall, G. B., Bradley, J. M., Schwinn, K. E., Bloor, S. J., Winefield, C. S., et al. (2006). Characterisation of aurone biosynthesis in Antirrhinum majus. Physiol. Plant 128, 593–603. doi: 10.1111/j.1399-3054.2006.00774.x
Davies, K., Albert, N., Zhou, Y., Schwinn, K. (2018). Functions of flavonoid and betalain pigments in abiotic stress tolerance in plants. Annu. Plant Rev. Online 1, 1–41. doi: 10.1002/9781119312994.apr0131
Davies, K. M., Jibran, R., Albert, N. W., Zhou, Y., Schwinn, K. E. (2020). “Conservation and divergence between bryophytes and angiosperms in the biosynthesis and regulation of flavonoid production,” in Recent Advances in Polyphenols Research. Eds. Quideau, S., Reed, J., De Freitas, V. (Oxford, UK: Wiley-Blackwell). in press.
Davin, L. B., Wang, H. B., Crowell, A. L., Bedgar, D. L., Martin, D. M., Sarkanen, S., et al. (1997). Stereoselective bimolecular phenoxy radical coupling by an auxiliary (dirigent) protein without an active center. Science 275, 362–366. doi: 10.1126/science.275.5298.362
de Vries, J., Archibald, J. M. (2018). Plant evolution: landmarks on the path to terrestrial life. New Phytol. 217, 1428–1434. doi: 10.1111/nph.14975
de Vries, J., de Vries, S., Slamovits, C. H., Rose, L. E., Archibald, J. M. (2017). How embryophytic is the biosynthesis of phenylpropanoids and their derivatives in streptophyte algae? Plant Cell Physiol. 58, 934–945. doi: 10.1093/pcp/pcx037
Delwiche, C. F., Cooper, E. D. (2015). The evolutionary origin of a terrestrial flora. Curr. Biol. 25, R899–R910. doi: 10.1016/j.cub.2015.08.029
Demarsy, E., Goldschmidt-Clermont, M., Ulm, R. (2017). Coping with ‘Dark Sides of the Sun’ through photoreceptor signaling. Trends Plant Sci. 23, 260–271. doi: 10.1016/j.tplants.2017.11.007
Dima, O., Morreel, K., Vanholme, B., Kim, H., Ralph, J., Boerjan, W. (2015). Small glycosylated lignin oligomers are stored in Arabidopsis leaf vacuoles. Plant Cell 27, 695–710. doi: 10.1105/tpc.114.134643
El Shoubaky, G. A., Abdel-Daim, M. M., Mansour, M. H., Salem, E. A. (2016). Isolation and identification of a flavone apigenin from marine red alga Acanthophora spicifera with antinociceptive and anti-Inflammatory activities. J. Exp. Neurosci. 10, 21–29. doi: 10.4137/JEN.S25096
Elumalai, P., Liu, H. L. (2011). Homology modeling and dynamics study of aureusidin synthase - an important enzyme in aurone biosynthesis of snapdragon flower. Int. J. Biol. Macromol. 49, 134–142. doi: 10.1016/j.ijbiomac.2011.03.017
Emiliani, G., Fondi, M., Fani, R., Gribaldo, S. (2009). A horizontal gene transfer at the origin of phenylpropanoid metabolism: a key adaptation of plants to land. Biol. Direct 4, 7. doi: 10.1186/1745-6150-4-7
Farag, M. A., Deavours, B. E., De Fatima, A., Naoumkina, M., Dixon, R. A., Sumner, L. W. (2009). Integrated metabolite and transcript profiling identify a biosynthetic mechanism for hispidol in Medicago truncatula cell cultures. Plant Physiol. 151, 1096–1113. doi: 10.1104/pp.109.141481
Feller, A., Machemer, K., Braun, E. L., Grotewold, E. (2011). Evolutionary and comparative analysis of MYB and bHLH plant transcription factors. Plant J. 66, 94–116. doi: 10.1111/j.1365-313X.2010.04459.x
Flores-Sandoval, E., Eklund, D. M., Bowman, J. L. (2015). A simple auxin transcriptional response system regulates multiple morphogenetic processes in the liverwort Marchantia polymorpha. PloS Genet. 11, e1005207. doi: 10.1371/journal.pgen.1005207
Foyer, C. H. (2018). Reactive oxygen species, oxidative signaling and the regulation of photosynthesis. Environ. Exp. Bot. 154, 134–142. doi: 10.1016/j.envexpbot.2018.05.003
Gang, D. R., Costa, M. A., Fujita, M., Dinkova-Kostova, A. T., Wang, H. B., Burlat, V., et al. (1999). Regiochemical control of monolignol radical coupling: a new paradigm for lignin and lignan biosynthesis. Chem. Biol. 6, 143–151. doi: 10.1016/S1074-5521(99)89006-1
Gasparetti, C., Faccio, G., Arvas, M., Buchert, J., Saloheimo, M., Kruus, K. (2010). Discovery of a new tyrosinase-like enzyme family lacking a C-terminally processed domain: production and characterization of an Aspergillus oryzae catechol oxidase. Appl. Microbiol. Biotechnol. 86, 213–226. doi: 10.1007/s00253-009-2258-3
Gasper, R., Effenberger, I., Kolesinski, P., Terlecka, B., Hofmann, E., Schaller, A. (2016). Dirigent protein mode of action revealed by the crystal structure of AtDIR6. Plant Physiol. 172, 2165–2175. doi: 10.1104/pp.16.01281
Gerdol, R., Bonora, A., Marchesini, R., Gualandri, R., Pancaldi, S. (1998). Growth response of Sphagnum capillifolium to nighttime temperature and nutrient level: Mechanisms and implications for global change. Arct. Alpine Res. 30, 388–395. doi: 10.2307/1552011
Gerdol, R. (1996). The seasonal growth pattern of Sphagnum magellanicum Brid. in different microhabitats on a mire in the southern Alps (Italy). Oecologia Montana 5, 13–20.
Glime, J. M. (2007). Bryophyte Ecology. Volume 1. Physiological Ecology (Michigan Technological University and the International Association of Bryologists).
Goiris, K., Muylaert, K., Voorspoels, S., Noten, B., Baart, G. J. E., De Cooman, L. (2014). Detection of flavonoids in microalgae from different evolutionary lineages. J. Phycol. 50, 483–492. doi: 10.1111/jpy.12180
Gould, K. S., Kuhn, D. N., Lee, D. W., Oberbauer, S. F. (1995). Why leaves are sometimes red. Nature 378, 241–242. doi: 10.1038/378241b0
Han, X. J., Wu, Y. F., Gao, S., Yu, H. N., Xu, R. X., Lou, H. X., et al. (2014). Functional characterization of a Plagiochasma appendiculatum flavone synthase I showing flavanone 2-hydroxylase activity. FEBS Lett. 588, 2307–2314. doi: 10.1016/j.febslet.2014.05.023
Harholt, J., Moestrup, O., Ulvskov, P. (2016). Why plants were terrestrial from the beginning. Trends Plant Sci. 21, 96–101. doi: 10.1016/j.tplants.2015.11.010
Holzinger, A., Pichrtová, M. (2016). Abiotic stress tolerance of Charophyte green algae: new challenges for Omics techniques. Front. Plant Sci. 7, 678. doi: 10.3389/fpls.2016.00678
Hong, L., Qian, Q., Tang, D., Wang, K., Li, M., Cheng, Z. (2012). A mutation in the rice chalcone isomerase gene causes the golden hull and internode 1 phenotype. Planta 236, 141–151. doi: 10.1007/s00425-012-1598-x
Hooijmaijers, C.a.M., Gould, K. S. (2007). Photoprotective pigments in red and green gametophytes of two New Zealand liverworts. N. Z. J. Bot. 45, 451–461. doi: 10.1080/00288250709509728
Hosmani, P. S., Kamiya, T., Danku, J., Naseer, S., Geldner, N., Guerinot, M. L., et al. (2013). Dirigent domain-containing protein is part of the machinery required for formation of the lignin-based Casparian strip in the root. Proc. Natl. Acad. Sci. U.S.A. 110, 14498–14503. doi: 10.1073/pnas.1308412110
Huber, M., Hintermann, G., Lerch, K. (1985). Primary structure of tyrosinase from Streptomyces glaucescens. Biochemistry 24, 6038–6044. doi: 10.1021/bi00343a003
Hughes, N. M. (2011). Winter leaf reddening in ‘evergreen’ species. New Phytol. 190, 573–581. doi: 10.1111/j.1469-8137.2011.03662.x
Hussain, T., Plunkett, B., Ejaz, M., Espley, R. V., Kayser, O. (2018). Identification of putative precursor genes for the biosynthesis of cannabinoid-like compound in Radula marginata. Front. Plant Sci. 9, 537. doi: 10.3389/fpls.2018.00537
Imbs, T. I., Zvyagintseva, T. N. (2018). Phlorotannins are polyphenolic metabolites of brown algae. Russ. J. Mar. Biol. 44, 263–273. doi: 10.1134/S106307401804003X
Ishizaki, K., Nishihama, R., Yamato, K. T., Kohchi, T. (2015). Molecular genetic tools and techniques for Marchantia polymorpha research. Plant Cell Physiol. 57, 262–270. doi: 10.1093/pcp/pcv097
Jez, J. M., Noel, J. P. (2002). Reaction Mechanism of chalcone isomerase: pH dependence, diffusion control, and product binding differences. J. Biol. Chem. 277, 1361–1369. doi: 10.1074/jbc.M109224200
Jiang, W., Yin, Q., Wu, R., Zheng, G., Liu, J., Dixon, R. A., et al. (2015). Role of a chalcone isomerase-like protein in flavonoid biosynthesis in Arabidopsis thaliana. J. Exp. Bot. 66, 7165–7179. doi: 10.1093/jxb/erv413
Jiang, N., Doseff, A. I., Grotewold, E. (2016). Flavones: From biosynthesis to health benefits. Plants 5, E27. doi: 10.3390/plants5020027
Jin, H., Cominelli, E., Bailey, P., Parr, A., Mehrtens, F., Jones, J., et al. (2000). Transcriptional repression by AtMYB4 controls production of UV-protecting sunscreens in Arabidopsis. EMBO J. 19, 6150–6161. doi: 10.1093/emboj/19.22.6150
Jorgensen, R. (1994). “The genetic origins of biosynthesis and light-responsive control of the chemical UV screen of land plants,” in Recent Advances in Phytochemistry, vol. 28 . Eds. Ellis, B. E., Kuroki, G. W., Stafford, H. A. (Boston, USA: Springer), 179–192. doi: 10.1007/978-1-4615-2544-8_7
Kaintz, C., Molitor, C., Thill, J., Kampatsikas, I., Michael, C., Halbwirth, H., et al. (2014). Cloning and functional expression in E. coli of a polyphenol oxidase transcript from Coreopsis grandiflora involved in aurone formation. FEBS Lett. 588, 3417–3426. doi: 10.1016/j.febslet.2014.07.034
Kaltenbach, M., Burke, J. R., Dindo, M., Pabis, A., Munsberg, F. S., Rabin, A., et al. (2018). Evolution of chalcone isomerase from a noncatalytic ancestor. Nat. Chem. Biol. 14, 988. doi: 10.1038/s41589-018-0042-3
Karsten, U., Holzinger, A. (2014). Green algae in alpine biological soil crust communities: acclimation strategies against ultraviolet radiation and dehydration. Biodivers. Conserv. 23, 1845–1858. doi: 10.1007/s10531-014-0653-2
Kenrick, P., Crane, P. R. (1997). The origin and early evolution of plants on land. Nature 389, 33–39. doi: 10.1038/37918
Kim, K. W., Smith, C. A., Daily, M. D., Cort, J. R., Davin, L. B., Lewis, N. G. (2015). Trimeric structure of (+)-pinoresinol-forming dirigent protein at 1.95 Å resolution with three isolated active sites. J. Biol. Chem. 290, 1308–1318. doi: 10.1074/jbc.M114.611780
Klejdus, B., Lojkova, L., Plaza, M., Snoblova, M., Sterbova, D. (2010). Hyphenated technique for the extraction and determination of isoflavones in algae: ultrasound-assisted supercritical fluid extraction followed by fast chromatography with tandem mass spectrometry. J. Chromatogr. A 1217, 7956–7965. doi: 10.1016/j.chroma.2010.07.020
Kny, L. (1890). “Bau und Entwickelung von Marchantia polymorpha,” in Botanische Wandtafeln, Vol VIII Abtheilung (Berlin, Germany: Paul Parey), 364–401.
Koes, R. E., Quattrocchio, F., Mol, J. N. M. (1994). The favonoid biosynthetic pathway in plants: function and evolution. Bioessays 16, 123–132. doi: 10.1002/bies.950160209
Kondou, Y., Miyagi, Y., Morito, T., Fujihira, K., Miyauchi, W., Moriyama, A., et al. (2019). Physiological function of photoreceptor UVR8 in UV-B tolerance in the liverwort Marchantia polymorpha. Planta 249, 1349–1364. doi: 10.1007/s00425-019-03090-w
Kubo, H., Nozawa, S., Hiwatashi, T., Kondou, Y., Nakabayashi, R., Mori, T., et al. (2018). Biosynthesis of riccionidins and marchantins is regulated by R2R3-MYB transcription factors in Marchantia polymorpha. J. Plant Res. 131, 849–864. doi: 10.1007/s10265-018-1044-7
Kunz, S., Becker, H. (1995). Cell wall pigment formation of in vitro cultures of the liverwort Ricciocarpos natans. Z.ür Naturforschung C 50, 235–240. doi: 10.1515/znc-1995-3-412
Kunz, S., Burkhardt, G., Becker, H. (1993). Riccionidins A and B, anthocyanidins from the cell walls of the liverwort Ricciocarpos natans. Phytochemistry 35, 233–235. doi: 10.1016/S0031-9422(00)90540-5
Kusano, M., Tohge, T., Fukushima, A., Kobayashi, M., Hayashi, N., Otsuki, H., et al. (2011). Metabolomics reveals comprehensive reprogramming involving two independent metabolic responses of Arabidopsis to UV-B light. Plant J. 67, 354–369. doi: 10.1111/j.1365-313X.2011.04599.x
Labandeira, C., Currano, E. (2013). The fossil record of plant-insect dynamics. Ann. Rev. Earth Planet. Sci. 41, 287–311. doi: 10.1146/annurev-earth-050212-124139
Labeeuw, L., Martone, P. T., Boucher, Y., Case, R. J. (2015). Ancient origin of the biosynthesis of lignin precursors. Biol. Direct 10, 23. doi: 10.1186/s13062-015-0052-y
Landi, M., Tattini, M., Gould, K. S. (2015). Multiple functional roles of anthocyanins in plant-environment interactions. Environ. Exp. Bot. 119, 4–17. doi: 10.1016/j.envexpbot.2015.05.012
Landry, L. G., Chapple, C. C., Last, R. L. (1995). Arabidopsis mutants lacking phenolic sunscreens exhibit enhanced ultraviolet-B injury and oxidative damage. Plant Physiol. 109, 1159–1166. doi: 10.1104/pp.109.4.1159
Lee, D. W., Gould, K. S. (2002). Anthocyanins in leaves and other vegetative organs: An introduction. Adv. Bot. Res. 37, 1–16. doi: 10.1016/S0065-2296(02)37040-X
Li, J., Ou-Lee, M., Amundson, R. G., Last, R. L. (1993). Arabidopsis flavonoid mutants are hypersensitive to UV-B irradiation. Plant Cell 5, 171–179. doi:10.1105/tpc.5.2.171
Li, C., Liu, S., Zhang, W., Chen, K., Zhang, P. (2019). Transcriptional profiling and physiological analysis reveal the critical roles of ROS-scavenging system in the Antarctic moss Pohlia nutans under Ultraviolet-B radiation. Plant Physiol. Biochem. 134, 113–122. doi: 10.1016/j.plaphy.2018.10.034
Ligrone, R., Duckett, J. G., Renzaglia, K. S. (2012). Major transitions in the evolution of early land plants: a bryological perspective. Ann. Bot. 109, 851–871. doi: 10.1093/aob/mcs017
Lillo, C., Lea, U. S., Ruoff, P. (2008). Nutrient depletion as a key factor for manipulating gene expression and product formation in different branches of the flavonoid pathway. Plant Cell Environ. 31, 587–601. doi: 10.1111/j.1365-3040.2007.01748.x
Lindenberg, J. B. W. (1836). Monographie der Riccieen. Nova Acta Physico-Medica Academiae Caesareae Leopoldino - Carolinae Naturae Curiosum XVIII, 361–504.
Lloyd, A., Brockman, A., Aguirre, L., Campbell, A., Bean, A., Cantero, A., et al. (2017). Advances in the MYB-bHLH-WD Repeat (MBW) pigment regulatory model: Addition of a WRKY factor and co-option of an anthocyanin MYB for betalain regulation. Plant Cell Physiol. 58, 1431–1441. doi: 10.1093/pcp/pcx075
Lu, Y., Ye, C., Che, J., Xu, X., Shao, D., Jiang, C., et al. (2019). Genomic sequencing, genome-scale metabolic network reconstruction, and in silico flux analysis of the grape endophytic fungus Alternaria sp. MG1. Microbial Cell Fact. 18, 13. doi: 10.1186/s12934-019-1063-7
Ma, D., Constabel, C. P. (2019). MYB repressors as regulators of phenylpropanoid metabolism in plants. Trends Plant Sci. 24, 275–289. doi: 10.1016/j.tplants.2018.12.003
Maloney, G. S., DiNapoli, K. T., Muday, G. K. (2014). The anthocyanin reduced tomato mutant demonstrates the role of flavonols in tomato lateral root and root hair development. Plant Physiol. 166, 614–631. doi: 10.1104/pp.114.240507
Manetas, Y. (2006). Why some leaves are anthocyanic and why most anthocyanic leaves are red? Flora 201, 163–177. doi: 10.1016/j.flora.2005.06.010
Markham, K. R. (1988). “Distribution of flavonoids in the lower plants and its evolutionary significance,” in The Flavonoids. Ed. Harborne, J. B. (London, UK: Chapman and Hall), 427–468. doi: 10.1007/978-1-4899-2913-6_12
Martens, S., Preuß, A., Matern, U. (2010). Multifunctional flavonoid dioxygenases: Flavonol and anthocyanin biosynthesis in Arabidopsis thaliana L. Phytochemistry 71, 1040–1049. doi: 10.1016/j.phytochem.2010.04.016
Matoba, Y., Kumagai, T., Yamamoto, A., Yoshitsu, H., Sugiyama, M. (2006). Crystallographic evidence that the dinuclear copper center of tyrosinase is flexible during catalysis. J. Biol. Chem. 281, 8981–8990. doi: 10.1074/jbc.M509785200
Merzlyak, M. N., Chivkunova, O. B. (2000). Light-stress-induced pigment changes and evidence for anthocyanin photoprotection in apples. J. Photochem. Photobiol. B Biol. 55, 155–163. doi: 10.1016/S1011-1344(00)00042-7
Meslet-Cladière, L., Delage, L., Leroux, C. J.-J., Goulitquer, S., Leblanc, C., Creis, E., et al. (2013). Structure/function analysis of a Type III polyketide synthase in the brown alga Ectocarpus siliculosus reveals a biochemical pathway in phlorotannin monomer biosynthesis. Plant Cell 25, 3089–3103. doi: 10.1105/tpc.113.111336
Miyahara, T., Sugishita, N., Ishida-Dei, M., Okamoto, E., Kouno, T., Cano, E. A., et al. (2018). Carnation I locus contains two chalcone isomerase genes involved in orange flower coloration. Breed. Sci. 68, 481–487. doi: 10.1270/jsbbs.18029
Molitor, C., Mauracher, S. G., Pargan, S., Mayer, R. L., Halbwirth, H., Rompel, A. (2015). Latent and active aurone synthase from petals of C. grandiflora: a polyphenol oxidase with unique characteristics. Planta 242, 519–537. doi: 10.1007/s00425-015-2261-0
Molitor, C., Mauracher, S. G., Rompel, A. (2016). Aurone synthase is a catechol oxidase with hydroxylase activity and provides insights into the mechanism of plant polyphenol oxidases. Proc. Nat. Acad. Sci. U.S.A. 113, E1806–E1815. doi: 10.1073/pnas.1523575113
Monforte, L., Soriano, G., Núñez-Olivera, E., Martínez-Abaigar, J. (2018). Cell compartmentation of ultraviolet-absorbing compounds: An underexplored tool related to bryophyte ecology, phylogeny and evolution. Funct. Ecol. 32, 882–893. doi: 10.1111/1365-2435.13048
Morales, L. O., Brosché, M., Vainonen, J., Jenkins, G. I., Wargent, J. J., Sipari, N., et al. (2013). Multiple roles for UV RESISTANCE LOCUS8 in regulating gene expression and metabolite accumulation in Arabidopsis under solar ultraviolet radiation. Plant Physiol. 161, 744–759. doi: 10.1104/pp.112.211375
Mouradov, A., Spangenberg, G. (2014). Flavonoids: a metabolic network mediating plants adaptation to their real estate. Front. Plant Sci. 5, 620. doi: 10.3389/fpls.2014.00620
Muhlemann, J. K., Younts, T. L. B., Muday, G. K. (2018). Flavonols control pollen tube growth and integrity by regulating ROS homeostasis during high-temperature stress. Proc. Nat. Acad. Sci. U.S.A. 115, E11188–EE1119. doi: 10.1073/pnas.1811492115
Myburg, A. A., Grattapaglia, D., Tuskan, G. A., Hellsten, U., Hayes, R. D., Grimwood, J., et al. (2014). The genome of Eucalyptus grandis. Nature 510, 356–362. doi: 10.1038/nature13308
Nagai, I. (1915). Ueber roten pigmentbildung beu einigen Marchantia-Arten. Bot. Mag. 29, 90–98. doi: 10.15281/jplantres1887.29.342_90
Nakayama, T., Yonekura-Sakakibara, K., Sato, T., Kikuchi, S., Fukui, Y., Fukuchi-Mizutani, M., et al. (2000). Aureusidin synthase: a polyphenol oxidase homolog responsible for flower coloration. Science 290, 1163–1166. doi: 10.1126/science.290.5494.1163
Napoli, C. A., Fahy, D., Wang, H. Y., Taylor, L. P. (1999). white anther: a Petunia mutant that abolishes pollen flavonol accumulation, induces male sterility, and is complemented by a chalcone synthase transgene. Plant Physiol. 120, 615–622. doi: 10.1104/pp.120.2.615
Nemesio-Gorriz, M., Blair, P. B., Dalman, K., Hammerbacher, A., Arnerup, J., Stenlid, J., et al. (2017). Identification of Norway Spruce MYB-bHLH-WDR transcription factor complex members linked to regulation of the flavonoid pathway. Front. Plant Sci. 8, 305–305. doi: 10.3389/fpls.2017.00305
Newsham, K. K., Geissler, P. A., Nicolson, M. J., Peat, H. J., Lewis-Smith, R. I. (2005). Sequential reduction of UV-B radiation in the field alters the pigmentation of an Antarctic leafy liverwort. Environ. Exp. Bot. 54, 22–32. doi: 10.1016/j.envexpbot.2004.05.006
Ngaki, M. N., Louie, G. V., Philippe, R. N., Manning, G., Pojer, F., Bowman, M. E., et al. (2012). Evolution of the chalcone-isomerase fold from fatty-acid binding to stereospecific catalysis. Nature 485, 530–533. doi: 10.1038/nature11009
Ono, E., Hatayama, M., Isono, Y., Sato, T., Watanabe, R., Yonekura-Sakakibara, K., et al. (2006). Localization of a flavonoid biosynthetic polyphenol oxidase in vacuoles. Plant J. 45, 133–143. doi: 10.1111/j.1365-313X.2005.02625.x
Pandith, S. A., Ramazan, S., Khan, M. I., Reshi, Z. A., Shah, M. A. (2020). Chalcone synthases (CHSs): the symbolic type III polyketide synthases. Planta 251, 15. doi: 10.1007/s00425-019-03307-y
Paniagua, C., Bilkova, A., Jackson, P., Dabravolski, S., Riber, W., Didi, V., et al. (2017). Dirigent proteins in plants: modulating cell wall metabolism during abiotic and biotic stress exposure. J. Exp. Bot. 68, 3287–3301. doi: 10.1093/jxb/erx141
Peer, W. A., Murphy, A. S. (2007). Flavonoids and auxin transport: modulators or regulators? Trends Plant Sci. 12, 556–563. doi: 10.1016/j.tplants.2007.10.003
Petersen, M., Simmonds, M. S. (2003). Rosmarinic acid. Phytochemistry 62, 121–125. doi: 10.1016/S0031-9422(02)00513-7
Philpott, M., Ferguson, L. R., Gould, K. S., Harris, P. J. (2009). Anthocyanidin-containing compounds occur in the periderm cell walls of the storage roots of sweet potato (Ipomoea batatas). J. Plant Physiol. 166, 1112–1117. doi: 10.1016/j.jplph.2008.12.007
Pichersky, E., Gang, D. R. (2000). Genetics and biochemistry of secondary metabolites in plants: an evolutionary perspective. Trends Plant Sci. 5, 439–445. doi: 10.1016/S1360-1385(00)01741-6
Pietrini, F., Iannelli, M. A., Massacci, A. (2002). Anthocyanin accumulation in the illuminated surface of maize leaves enhances protection from photo-inhibitory risks at low temperature, without further limitation to photosynthesis. Plant Cell Environ. 25, 1251–1259. doi: 10.1046/j.1365-3040.2002.00917.x
Pollastri, S., Tattini, M. (2011). Flavonols: old compounds for old roles. Ann. Bot. 108, 1225–1233. doi: 10.1093/aob/mcr234
Polturak, G., Aharoni, A. (2019). Advances and future directions in betalain metabolic engineering. New Phytol. 224, 1472–1478. doi: 10.1111/nph.15973
Puttick, M. N., Morris, J. L., Williams, T. A., Cox, C. J., Edwards, D., Kenrick, P., et al. (2018). The interrelationships of land plants and the nature of the ancestral embryophyte. Curr. Biol. 28, 733–745. doi: 10.1016/j.cub.2018.01.063
Ramos, P., Guajardo, J., Moya-León, M. A., Herrera, R. (2016). A differential distribution of auxin and flavonols in radiata pine stem seedlings exposed to inclination. Tree Genet. Genomes 12, 42. doi: 10.1007/s11295-016-1003-1
Reboledo, G., Del Campo, R., Alvarez, A., Montesano, M., Mara, H., Ponce De Leon, I. (2015). Physcomitrella patens activates defense responses against the pathogen Colletotrichum gloeosporioides. Int. J. Mol. Sci. 16, 22280–22298. doi: 10.3390/ijms160922280
Reeb, C., Marline, L., Rabeau, L., Andriamanantena, A., Andriamiarisoa, R. L., Ranarijaona, H., et al. (2018). A survey of Marchantiales from Madagascar. Acta Biol. Plant Agriensis 6, 3–72. doi: 10.21406/abpa.2018.6.3
Rensing, S. A. (2018). Great moments in evolution: the conquest of land by plants. Curr. Opin. Plant Biol. 42, 49–54. doi: 10.1016/j.pbi.2018.02.006
Rozema, J., Bjorn, L. O., Bornman, J. F., Gaberscik, A., Hader, D. P., Trost, T., et al. (2002). The role of UV-B radiation in aquatic and terrestrial ecosystems–an experimental and functional analysis of the evolution of UV-absorbing compounds. J. Photochem. Photobiol. B: Biol. 66, 2–12. doi: 10.1016/S1011-1344(01)00269-X
Rubin, G., Tohge, T., Matsuda, F., Saito, K., Scheible, W. R. (2009). Members of the LBD family of transcription factors repress anthocyanin synthesis and affect additional nitrogen responses in Arabidopsis. Plant Cell 21, 3567–3584. doi: 10.1105/tpc.109.067041
Rudolph, H., Kabsch, U., Schmidt-Stohn, G. (1977). Änderungen des Chloroplastenpigment-Spiegels bei Sphagnum magellanicum im Verlauf der Synthese von Sphagnorubin und anderer membranochromer Pigmente. Z. für Pflanzenphysiologie 82, 107–116. doi: 10.1016/S0044-328X(77)80219-5
Sendovski, M., Kanteev, M., Ben-Yosef, V. S., Adir, N., Fishman, A. (2011). First structures of an active bacterial tyrosinase reveal copper plasticity. J. Mol. Biol. 405, 227–237. doi: 10.1016/j.jmb.2010.10.048
Shimamura, M. (2016). Marchantia polymorpha: Taxonomy, phylogeny and morphology of a model system. Plant Cell Physiol. 57, 230–256. doi: 10.1093/pcp/pcv192
Shuster, V., Fishman, A. (2009). Isolation, cloning and characterization of a tyrosinase with improved activity in organic solvents from Bacillus megaterium. J. Mol. Microbiol. Biotechnol. 17, 188–200. doi: 10.1159/000233506
Soriano, G., Cloix, C., Heilmann, M., Nunez-Olivera, E., Martinez-Abaigar, J., Jenkins, G. I. (2018a). Evolutionary conservation of structure and function of the UVR8 photoreceptor from the liverwort Marchantia polymorpha and the moss Physcomitrella patens. New Phytol. 217, 151–162. doi: 10.1111/nph.14767
Soriano, G., Del-Castillo-Alonso, M. A., Monforte, L., Nunez-Olivera, E., Martinez-Abaigar, J. (2018b). Phenolic compounds from different bryophyte species and cell compartments respond specifically to ultraviolet radiation, but not particularly quickly. Plant Physiol. Biochem. 134, 137–144. doi: 10.1016/j.plaphy.2018.07.020
Soriano, G., Del-Castillo-Alonso, M. A., Monforte, L., Nunez-Olivera, E., Martinez-Abaigar, J. (2018c). First data on the effects of ultraviolet radiation on phenolic compounds in the model hornwort Anthoceros agrestis. Cryptogamie Bryologie 39, 201–211. doi: 10.7872/cryb/v39.iss2.2018.201
Stafford, H. A. (1991). Flavonoid evolution: An enzymic approach. Plant Physiol. 96, 680–685. doi: 10.1104/pp.96.3.680
Stark, L. R., Greenwood, J. L., Brinda, J. C. (2017). Desiccated Syntrichia ruralis shoots regenerate after 20 years in the herbarium. J. Bryol. 39, 85–93. doi: 10.1080/03736687.2016.1176307
Stich, K., Eidenberger, T., Wurst, F., Forkmann, G. (1992). Enzymatic conversion of dihydroflavonols to flavan-3,4-diols using flower extracts of Dianthus caryophyllus L. (carnation). Planta 187, 103–108. doi: 10.1007/BF00201630
Sugano, S. S., Nishihama, R., Shirakawa, M., Takagi, J., Matsuda, Y., Ishida, S., et al. (2018). Efficient CRISPR/Cas9-based genome editing and its application to conditional genetic analysis in Marchantia polymorpha. PloS One 13, e0205117. doi: 10.1371/journal.pone.0205117
Szövényi, P., Frangedakis, E., Ricca, M., Quandt, D., Wicke, S., Langdale, J. A. (2015). Establishment of Anthoceros agrestis as a model species for studying the biology of hornworts. BMC Plant Biol. 15, 98. doi: 10.1186/s12870-015-0481-x
Taniguchi, S., Yazaki, K., Yabuuchi, R., Kawakami, K., Ito, H., Hatano, T., et al. (2000). Galloylglucoses and riccionidin A in Rhus javanica adventitious root cultures. Phytochemistry 53, 357–363. doi: 10.1016/S0031-9422(99)00399-4
Taylor, L. P., Grotewold, E. (2005). Flavonoids as developmental regulators. Curr. Opin. Plant Biol. 8, 317–323. doi: 10.1016/j.pbi.2005.03.005
Thamil Arasan, S. K., Park, J. I., Ahmed, N. U., Jung, H. J., Hur, Y., Kang, K. K., et al. (2013). Characterization and expression analysis of dirigent family genes related to stresses in Brassica. Plant Physiol. Biochem. 67, 144–153. doi: 10.1016/j.plaphy.2013.02.030
Tohge, T., Watanabe, M., Hoefgen, R., Fernie, A. R. (2013). The evolution of phenylpropanoid metabolism in the green lineage. Crit. Rev. Biochem. Mol. Biol. 48, 123–152. doi: 10.3109/10409238.2012.758083
Tran, L. T., Taylor, J. S., Constabel, C. P. (2012). The polyphenol oxidase gene family in land plants: Lineage-specific duplication and expansion. BMC Genomics 13, 395. doi: 10.1186/1471-2164-13-395
Uchida, K., Akashi, T., Aoki, T. (2017). The missing link in Leguminous pterocarpan biosynthesis is a dirigent domain-containing protein with isoflavanol dehydratase activity. Plant Cell Physiol. 58, 398–408. doi: 10.1093/pcp/pcw213
Van der Meer, I. M., Stam, M. E., Van Tunen, A. J., Mol, J. N., Stuitje, A. R. (1992). Antisense inhibition of flavonoid biosynthesis in Petunia anthers results in male sterility. Plant Cell 4, 253–262. doi: 10.1105/tpc.4.3.253
Verhoeven, J. T. A., Liefveld, W. M. (1997). The ecological significance of organochemical compounds in Sphagnum. Acta Botanica Neerlandica 46, 117–130. doi: 10.1111/plb.1997.46.2.117
Vitt, D. H., Crandall-Stotler, B. J., Wood, A. (2012). “Bryophytes: survival in a dry world through tolerance and avoidance,” in Plant Ecology and Evolution in Harsh Environments. Eds. Rajakaruna, N., Boyd, R. S., Harris, T. B. (New York, USA: Nova Science), 267–295.
Vowinkel, E. (1975). Torfmoosmembranochrome, 2. Die struktur des sphagnorubins. Chemische Berichte-Recueil 108, 1166–1181. doi: 10.1002/cber.19751080423
Wang, Y., Chantreau, M., Sibout, R., Hawkins, S. (2013). Plant cell wall lignification and monolignol metabolism. Front. Plant Sci. 4, 220. doi: 10.3389/fpls.2013.00220
Wang, N., Liu, W., Zhang, T., Jiang, S., Xu, H., Wang, Y., et al. (2018). Transcriptomic analysis of red-fleshed apples reveals the novel role of MdWRKY11 in flavonoid and anthocyanin biosynthesis. J. Agric. Food Chem. 66, 7076–7086. doi: 10.1021/acs.jafc.8b01273
Waterman, M. J., Nugraha, A. S., Hendra, R., Ball, G. E., Robinson, S. A., Keller, P. A. (2017). Antarctic moss biflavonoids show high antioxidant and ultraviolet-screening activity. J. Nat. Prod. 80, 2224–2231. doi: 10.1021/acs.jnatprod.7b00085
Waterman, M. J., Bramley-Alves, J., Miller, R. E., Keller, P. A., Robinson, S. A. (2018). Photoprotection enhanced by red cell wall pigments in three East Antarctic mosses. Biol. Res. 51, 49. doi: 10.1186/s40659-018-0196-1
Wickett, N. J., Mirarab, S., Nguyen, N., Warnow, T., Carpenter, E., Matasci, N. (2014). Phylotranscriptomic analysis of the origin and early diversification of land plants. Proc. Natl. Acad. Sci. U.S.A. 111, E4859–E4868. doi: 10.1073/pnas.1323926111
Wolf, L., Rizzini, L., Stracke, R., Ulm, R., Rensing, S. A. (2010). The molecular and physiological responses of Physcomitrella patens to Ultraviolet-B radiation. Plant Physiol. 153, 1123–1134. doi: 10.1104/pp.110.154658
Wu, Y. F., Zhao, Y., Liu, X. Y., Gao, S., Cheng, A. X., Lou, H. X. (2018). A bHLH transcription factor regulates bisbibenzyl biosynthesis in the liverwort Plagiochasma appendiculatum. Plant Cell Physiol. 59, 1187–1199. doi: 10.1093/pcp/pcy053
Yin, R., Ulm, R. (2017). How plants cope with UV-B: from perception to response. Curr. Opin. Plant Biol. 37, 42–48. doi: 10.1016/j.pbi.2017.03.013
Ylstra, B., Muskens, M., Van Tunen, A. J. (1996). Flavonols are not essential for fertilization in Arabidopsis thaliana. Plant Mol. Biol. 32, 1155–1158. doi: 10.1007/BF00041399
Yonekura-Sakakibara, K., Higashi, Y., Nakabayashi, R. (2019). The origin and evolution of plant flavonoid metabolism. Front. Plant Sci. 10, 943. doi: 10.3389/fpls.2019.00943
Zhao, Q., Dixon, R. A. (2014). Altering the cell wall and its impact on plant disease: from forage to bioenergy. Ann. Rev. Phytopathol. 52, 69–91. doi: 10.1146/annurev-phyto-082712-102237
Zhao, D., Hamilton, J. P., Bhat, W. W., Johnson, S. R., Godden, G. T., Kinser, T. J., et al. (2019). A chromosomal-scale genome assembly of Tectona grandis reveals the importance of tandem gene duplication and enables discovery of genes in natural product biosynthetic pathways. GigaScience 8, giz005. doi: 10.1093/gigascience/giz005
Zhong, B., Xi, Z., Goremykin, V. V., Fong, R., Mclenachan, P. A., Novis, P. M., et al. (2014). Streptophyte algae and the origin of land plants revisited using heterogeneous models with three new algal chloroplast genomes. Mol. Biol. Evol. 31, 177–183. doi: 10.1093/molbev/mst200
Keywords: anthocyanin, auronidin, dirigent, polyphenol oxidase, transcription factor
Citation: Davies KM, Jibran R, Zhou Y, Albert NW, Brummell DA, Jordan BR, Bowman JL and Schwinn KE (2020) The Evolution of Flavonoid Biosynthesis: A Bryophyte Perspective. Front. Plant Sci. 11:7. doi: 10.3389/fpls.2020.00007
Received: 29 October 2019; Accepted: 07 January 2020;
Published: 04 February 2020.
Edited by:
Kazuki Saito, RIKEN Center for Sustainable Resource Science (CSRS), JapanReviewed by:
Stefan Martens, Fondazione Edmund Mach, ItalyBrenda S.J. Winkel, Virginia Tech, United States
Copyright © 2020 Davies, Jibran, Zhou, Albert, Brummell, Jordan, Bowman and Schwinn. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Kevin M. Davies, kevin.davies@plantandfood.co.nz