Computational Modeling of Anthocyanin Pathway Evolution: Biases, Hotspots, and Trade-offs

Floral pigments and their perception by avian pollinators in three Chilean Puya species

The Chilean Puya species, Puya coerulea var. violacea and P. chilensis bear blue and pale-yellow flowers, respectively, while P. alpestris considered to be their hybrid-derived species has unique turquoise flowers. In this study, the chemical basis underlying the different coloration of the three Puya species was explored. We first isolated and identified three anthocyanins: delphinidin 3,3',5'-tri-O-glucoside, delphinidin 3,3'-di-O-glucoside and delphinidin 3-O-glucoside; seven flavonols: quercetin 3-O-rutinoside-3'-O-glucoside, quercetin 3,3'-di-O-glucoside, quercetin 3-O-rutinoside, isorhamnetin 3-O-rutinoside, myricetin 3,3',5'-tri-O-glucoside, myricetin 3,3'-di-O-glucoside and laricitrin 3,5'-di-O-glucoside; and six flavones: luteolin 4'-O-glucoside, apigenin 4'-O-glucoside, tricetin 4'-O-glucoside, tricetin 3',5'-di-O-glucoside, tricetin 3'-O-glucoside and selagin 5'-O-glucoside, which is a previously undescribed flavone, from their petals. We also compared compositions of floral flavonoid and their aglycone among these species, which suggested that the turquoise species P. alpestris has an essentially intermediate composition between the blue and pale-yellow species. The vacuolar pH was relatively higher in the turquoise (pH 6.2) and pale-yellow (pH 6.2) flower species, while that of blue flower species was usual (pH 5.2). The flower color was reconstructed in vitro using isolated anthocyanin, flavonol and flavone at neutral and acidic pH, and its color was analyzed by reflectance spectra and the visual modeling of their avian pollinators. The modeling demonstrated that the higher pH of the turquoise and pale-yellow species enhances the chromatic contrast and spectral purity. The precise regulation of flower color by flavonoid composition and vacuolar pH may be adapted to the visual perception of their avian pollinator vision.

Model-assisted analysis for tuning anthocyanin composition in grape berries

Anthocyanin composition is responsible for the red colour of grape berries and wines, and contributes to their organoleptic quality. However, anthocyanin biosynthesis is under genetic, developmental and environmental regulation, making its targeted fine-tuning challenging. We constructed a mechanistic model to simulate the dynamics of anthocyanin composition throughout grape ripening in Vitis vinifera, employing a consensus anthocyanin biosynthesis pathway. The model was calibrated and validated using six datasets from eight cultivars and 37 growth conditions. Tuning the transformation and degradation parameters allowed us to accurately simulate the accumulation process of each individual anthocyanin under different environmental conditions. The model parameters were robust across environments for each genotype. The coefficients of determination (R2) for the simulated versus observed values for the six datasets ranged from 0.92 to 0.99, while the relative root mean square errors (RRMSEs) were between 16.8 and 42.1 %. The leave-one-out cross-validation for three datasets showed R2 values of 0.99, 0.96 and 0.91, and RRMSE values of 28.8, 32.9 and 26.4 %, respectively, suggesting a high prediction quality of the model. Model analysis showed that the anthocyanin profiles of diverse genotypes are relatively stable in response to parameter perturbations. Virtual experiments further suggested that targeted anthocyanin profiles may be reached by manipulating a minimum of three parameters, in a genotype-dependent manner. This model presents a promising methodology for characterizing the temporal progression of anthocyanin composition, while also offering a logical foundation for bioengineering endeavours focused on precisely adjusting the anthocyanin composition of grapes.

Plant morpho evo-devo

Evo-devo is often thought of as being the study of which genes underlie which phenotypes. However, evo-devo is much more than this, especially in plant science. In leaf scars along stems, cell changes across wood growth rings, or flowers along inflorescences, plants trace a record of their own development. Plant morpho evo-devo provides data that genes could never furnish on themes such as heterochrony, the evolution of temporal phenotypes, modularity, and phenotype-first evolution. As plant science surges into increasingly -omic realms, it is essential to keep plant morpho evo-devo in full view as an honored member of the evo-devo canon, ensuring that plant scientists can, wherever they are, generate fundamental insights at the appropriate level of biological organization.

Insight into the molecular mechanisms of leaf coloration in Cymbidium ensifolium

Cymbidiumensifolium L. is a significant ornamental plant in Orchidaceae. Aside from its attractive flowers, its leaf coloration is also an important ornamental trait. However, there is an apparent lack of studies concerning the intricate mechanism of leaf coloration in C. ensifolium. In this study, we report a systematic evaluation of leaf coloration utilizing transcriptome and metabolome profiles of purple, yellow, and green leaves. In total, 40 anthocyanins and 67 flavonoids were quantified along with chlorophyll content. The tissue–transcriptome profile identified 26,499 differentially expressed genes (DEGs). The highest chlorophyll contents were identified in green leaves, followed by yellow and purple leaves. We identified key anthocyanins and flavonoids associated with leaf coloration, including cyanidin-3-O-sophoroside, naringenin-7-O-glucoside, delphinidin, cyanidin, petunidin, and quercetin, diosmetin, sinensetin, and naringenin chalcone. Moreover, genes encoding UDP-glucoronosyl, UDP-glucosyl transferase, chalcone synthesis, flavodoxin, cytochrome P450, and AMP-binding enzyme were identified as key structural genes affecting leaf coloration in C. ensifolium. In summary, copigmentation resulting from several key metabolites modulated by structural genes was identified as governing leaf coloration in C. ensifolium. Further functional verification of the identified DEGs and co-accumulation of metabolites can provide a tool to modify leaf color and improve the aesthetic value of C. ensifolium.

The Common Bean V Gene Encodes Flavonoid 3'5' Hydroxylase: A Major Mutational Target for Flavonoid Diversity in Angiosperms

The classic V (violet, purple) gene of common bean (Phaseolus vulgaris) functions in a complex genetic network that controls seed coat and flower color and flavonoid content. V was cloned to understand its role in the network and the evolution of its orthologs in the Viridiplantae. V mapped genetically to a narrow interval on chromosome Pv06. A candidate gene was selected based on flavonoid analysis and confirmed by recombinational mapping. Protein and domain modeling determined V encodes flavonoid 3'5' hydroxylase (F3'5'H), a P450 enzyme required for the expression of dihydromyricetin-derived flavonoids in the flavonoid pathway. Eight recessive haplotypes, defined by mutations of key functional domains required for P450 activities, evolved independently in the two bean gene pools from a common ancestral gene. V homologs were identified in Viridiplantae orders by functional domain searches. A phylogenetic analysis determined F3'5'H first appeared in the Streptophyta and is present in only 41% of Angiosperm reference genomes. The evolutionarily related flavonoid pathway gene flavonoid 3' hydroxylase (F3'H) is found nearly universally in all Angiosperms. F3'H may be conserved because of its role in abiotic stress, while F3'5'H evolved as a major target gene for the evolution of flower and seed coat color in plants.

Transcription factors evolve faster than their structural gene targets in the flavonoid pigment pathway

Dissecting the relationship between gene function and substitution rates is key to understanding genome-wide patterns of molecular evolution. Biochemical pathways provide powerful systems for investigating this relationship because the functional role of each gene is often well characterized. Here, we investigate the evolution of the flavonoid pigment pathway in the colorful Petunieae clade of the tomato family (Solanaceae). This pathway is broadly conserved in plants, both in terms of its structural elements and its MYB, basic helix–loop–helix, and WD40 transcriptional regulators, and its function has been extensively studied, particularly in model species of petunia. We built a phylotranscriptomic data set for 69 species of Petunieae to infer patterns of molecular evolution across pathway genes and across lineages. We found that transcription factors exhibit faster rates of molecular evolution (dN/dS) than their targets, with the highly specialized MYB genes evolving fastest. Using the largest comparative data set to date, we recovered little support for the hypothesis that upstream enzymes evolve slower than those occupying more downstream positions, although expression levels do predict molecular evolutionary rates. Although shifts in floral pigmentation were only weakly related to changes affecting coding regions, we found a strong relationship with the presence/absence patterns of MYB transcripts. Intensely pigmented species express all three main MYB anthocyanin activators in petals, whereas pale or white species express few or none. Our findings reinforce the notion that pathway regulators have a dynamic history, involving higher rates of molecular evolution than structural components, along with frequent changes in expression during color transitions.

Elevational divergence in pigmentation plasticity is associated with selection and pigment biochemistry

Phenotypic plasticity is predicted to evolve in environmentally variable habitats, or those experiencing a high frequency of strong selection. However, the evolution of plasticity may be constrained by costs or physiological constraints. In flowers, UV-absorbing pigmentation ameliorates UV damage to pollen, and is linked with elevated UV exposure. Whether plasticity contributes to this pattern remains unclear. Petals of Argentina anserina have larger UV-absorbing petal areas at high elevations where they experience higher and more variable UV exposure compared to low elevations. We measured UV-induced pigmentation plasticity in high- and low-elevation populations (hereafter, 'high, 'low'), and selection on pigmentation via male fitness. We dissected UV pigment biochemistry using metabolomics to explore biochemical mechanisms underlying plasticity. High displayed positive UV-induced pigmentation plasticity but low lacked plasticity. Selection favored elevated pigmentation under UV in high, supporting adaptive plasticity. In high, UV-absorption was conferred by flavonoids produced in one flavonoid pathway branch. However, in low, UV-absorption was associated with many compounds spanning many branches. Elevated plasticity was thus associated with reduced pigment diversity. The results are consistent with adaptive floral pigmentation plasticity in more extreme and variable environments. We discuss how biochemical underpinnings of pigmentation may permit or constrain the evolution of pigmentation plasticity. This article is protected by copyright. All rights reserved.

TbMYC4A Is a Candidate Gene Controlling the Blue Aleurone Trait in a Wheat-Triticum boeoticum Substitution Line

Triticum boeoticum Boiss (A^bA^b, 2n = 2x = 14) is one of the sources of the blue grain trait controlled by blue aleurone layer 2 (Ba2). However, the underlying genes have not been cloned. In this study, a transcriptomic comparison between a blue-grained wheat-T. boeoticum substitution line and its wheat parent identified 41 unigenes related to anthocyanin biosynthesis and 29 unigenes related to transport. The bHLH transcription factor gene TbMYC4A showed a higher expression level in the blue-grained substitution line. TbMYC4A contained the three characteristic bHLH transcription factor domains (bHLH-MYC_N, HLH and ACT-like) and clustered with genes identified from other wheat lines with the blue grain trait derived from other Triticeae species. TbMYC4A overexpression confirmed that it was a functional bHLH transcription factor. The analysis of a TbMYC4A-specific marker showed that the gene was also present in T. boeoticum and T. monococcum with blue aleurone but absent in other Triticeae materials with white aleurone. These results indicate that TbMYC4A is a candidate gene of Ba2 controlling the blue aleurone trait. The isolation of TbMYC4A is helpful for further clarifying the genetic mechanism of the blue aleurone trait and is of great significance for breeding blue-grained wheat varieties.

Quantitative dissection of color patterning in the foliar ornamental coleus

Coleus (Coleus scutellarioides) is a popular ornamental plant that exhibits a diverse array of foliar color patterns. New cultivars are currently hand selected by both amateur and experienced plant breeders. In this study, we reimagine breeding for color patterning using a quantitative color analysis framework. Despite impressive advances in high-throughput data collection and processing, complex color patterns remain challenging to extract from image datasets. Using a phenotyping approach called "ColourQuant," we extract and analyze pigmentation patterns from one of the largest coleus breeding populations in the world. Working with this massive dataset, we can analyze quantitative relationships between maternal plants and their progeny, identify features that underlie breeder-selections, and collect and compare public input on trait preferences. This study is one of the most comprehensive explorations into complex color patterning in plant biology and provides insights and tools for exploring the color pallet of the plant kingdom.

Computational modeling and evolutionary implications of biochemical reactions in bacterial microcompartments

Bacterial microcompartments (BMCs) are protein-encapsulated compartments found across at least 23 bacterial phyla. BMCs contain a variety of metabolic processes that share the commonality of toxic or volatile intermediates, oxygen-sensitive enzymes and cofactors, or increased substrate concentration for magnified reaction rates. These compartmentalized reactions have been computationally modeled to explore the encapsulated dynamics, ask evolutionary-based questions, and develop a more systematic understanding required for the engineering of novel BMCs. Many crucial aspects of these systems remain unknown or unmeasured, such as substrate permeabilities across the protein shell, feasibility of pH gradients, and transport rates of associated substrates into the cell. This review explores existing BMC models, dominated in the literature by cyanobacterial carboxysomes, and highlights potentially important areas for exploration.

Convergent Biochemical Pathways for Xanthine Alkaloid Production in Plants Evolved from Ancestral Enzymes with Different Catalytic Properties

Convergent evolution is widespread but the extent to which common ancestral conditions are necessary to facilitate the independent acquisition of similar traits remains unclear. In order to better understand how ancestral biosynthetic catalytic capabilities might lead to convergent evolution of similar modern-day biochemical pathways, we resurrected ancient enzymes of the caffeine synthase (CS) methyltransferases that are responsible for theobromine and caffeine production in flowering plants. Ancestral CS enzymes of Theobroma, Paullinia, and Camellia exhibited similar substrate preferences but these resulted in the formation of different sets of products. From these ancestral enzymes, descendants with similar substrate preference and product formation independently evolved after gene duplication events in Theobroma and Paullinia. Thus, it appears that the convergent modern-day pathways likely originated from ancestral pathways with different inferred flux. Subsequently, the modern-day enzymes originated independently via gene duplication and their convergent catalytic characteristics evolved to partition the multiple ancestral activities by different mutations that occurred in homologous regions of the ancestral proteins. These results show that even when modern-day pathways and recruited genes are similar, the antecedent conditions may be distinctive such that different evolutionary steps are required to generate convergence.

Changes at a Critical Branchpoint in the Anthocyanin Biosynthetic Pathway Underlie the Blue to Orange Flower Color Transition in Lysimachia arvensis

Anthocyanins are the primary pigments contributing to the variety of flower colors among angiosperms and are considered essential for survival and reproduction. Anthocyanins are members of the flavonoids, a broader class of secondary metabolites, of which there are numerous structural genes and regulators thereof. In western European populations of Lysimachia arvensis, there are blue- and orange-petaled individuals. The proportion of blue-flowered plants increases with temperature and daylength yet decreases with precipitation. Here, we performed a transcriptome analysis to characterize the coding sequences of a large group of flavonoid biosynthetic genes, examine their expression and compare our results to flavonoid biochemical analysis for blue and orange petals. Among a set of 140 structural and regulatory genes broadly representing the flavonoid biosynthetic pathway, we found 39 genes with significant differential expression including some that have previously been reported to be involved in similar flower color transitions. In particular, F3'5'H and DFR, two genes at a critical branchpoint in the ABP for determining flower color, showed differential expression. The expression results were complemented by careful examination of the SNPs that differentiate the two color types for these two critical genes. The decreased expression of F3'5'H in orange petals and differential expression of two distinct copies of DFR, which also exhibit amino acid changes in the color-determining substrate specificity region, strongly correlate with the blue to orange transition. Our biochemical analysis was consistent with the transcriptome data indicating that the shift from blue to orange petals is caused by a change from primarily malvidin to largely pelargonidin forms of anthocyanins. Overall, we have identified several flavonoid biosynthetic pathway loci likely involved in the shift in flower color in L. arvensis and even more loci that may represent the complex network of genetic and physiological consequences of this flower color polymorphism.

Structure and contingency determine mutational hotspots for flower color evolution

Evolutionary genetic studies have uncovered abundant evidence for genomic hotspots of phenotypic evolution, as well as biased patterns of mutations at those loci. However, the theoretical basis for this concentration of particular types of mutations at particular loci remains largely unexplored. In addition, historical contingency is known to play a major role in evolutionary trajectories, but has not been reconciled with the existence of such hotspots. For example, do the appearance of hotspots and the fixation of different types of mutations at those loci depend on the starting state and/or on the nature and direction of selection? Here, we use a computational approach to examine these questions, focusing the anthocyanin pigmentation pathway, which has been extensively studied in the context of flower color transitions. We investigate two transitions that are common in nature, the transition from blue to purple pigmentation and from purple to red pigmentation. Both sets of simulated transitions occur with a small number of mutations at just four loci and show strikingly similar peaked shapes of evolutionary trajectories, with the mutations of the largest effect occurring early but not first. Nevertheless, the types of mutations (biochemical vs. regulatory) as well as their direction and magnitude are contingent on the particular transition. These simulated color transitions largely mirror findings from natural flower color transitions, which are known to occur via repeated changes at a few hotspot loci. Still, some types of mutations observed in our simulated color evolution are rarely observed in nature, suggesting that pleiotropic effects further limit the trajectories between color phenotypes. Overall, our results indicate that the branching structure of the pathway leads to a predictable concentration of evolutionary change at the hotspot loci, but the types of mutations at these loci and their order is contingent on the evolutionary context.

Celebrating a New Division of Botany at SICB: An Introduction to the Integrative Plant Biology Symposium

The Society for Integrative and Comparative Biology (SICB) should, in theory, be a home for scientists working across the entire Tree of Life. In practice, SICB has remained principally a society that supports integrative zoological research. Here we highlight a broad collection of what we consider to the best in integrative and comparative plant biology, gathered together for a special symposium at the 2019 SICB meeting. This symposium and special issue mark the initiation of a new Division of Botany within SICB, which we hope will usher in a new era of SICB where botanists and zoologists engage, collaborate, and celebrate together in this especially creative period of integrative and comparative biology.

Chemical Basis of Floral Color Signals in Gesneriaceae: The Effect of Alternative Anthocyanin Pathways

Changes in floral pigmentation can have dramatic effects on angiosperm evolution by making flowers either attractive or inconspicuous to different pollinator groups. Flower color largely depends on the type and abundance of pigments produced in the petals, but it is still unclear whether similar color signals rely on same biosynthetic pathways and to which extent the activation of certain pathways influences the course of floral color evolution. To address these questions, we investigated the physical and chemical aspects of floral color in the Neotropical Gesnerioideae (ca. 1,200 spp.), in which two types of anthocyanins, hydroxyanthocyanins, and deoxyanthocyanins, have been recorded as floral pigments. Using spectrophotometry, we measured flower reflectance for over 150 species representing different clades and pollination syndromes. We analyzed these reflectance data to estimate how the Gesnerioideae flowers are perceived by bees and hummingbirds using the visual system models of these pollinators. Floral anthocyanins were further identified using high performance liquid chromatography coupled to mass spectrometry. We found that orange/red floral colors in Gesnerioideae are produced either by deoxyanthocyanins (e.g., apigenidin, luteolinidin) or hydroxyanthocyanins (e.g., pelargonidin). The presence of deoxyanthocyanins in several lineages suggests that the activation of the deoxyanthocyanin pathway has evolved multiple times in the Gesnerioideae. The hydroxyanthocyanin-producing flowers span a wide range of colors, which enables them to be discriminated by hummingbirds or bees. By contrast, color diversity among the deoxyanthocyanin-producing species is lower and mainly represented at longer wavelengths, which is in line with the hue discrimination optima for hummingbirds. These results indicate that Gesnerioideae have evolved two different biochemical mechanisms to generate orange/red flowers, which is associated with hummingbird pollination. Our findings also suggest that the activation of the deoxyanthocyanin pathway has restricted flower color diversification to orange/red hues, supporting the potential constraining role of this alternative biosynthetic pathway on the evolutionary outcome of phenotypical and ecological diversification.

A bHLH transcription factor TsMYC2 is associated with the blue grain character in triticale (Triticum × Secale)

RNA-Seq was employed to compare the transcriptome differences between the triticale lines and to identify the key gene responsible for the blue aleurone trait. The accumulation of anthocyanins in the aleurone of triticale results in the formation of the blue-grained trait, but the identity of the genes associated with anthocyanin biosynthesis in the aleurone has not yet been reported. In this manuscript, RNA-Seq was employed to compare the transcriptome differences between the triticale lines HM13 (blue aleurone) and HM5 (white aleurone), and to identify the key genes responsible for the blue aleurone trait. There were 32,406 differentially expressed genes between HM13 and HM5. Seventy-three unigenes were homologous to the structural genes related to anthocyanin biosynthesis, and the average transcript level of the structural genes was higher in HM13 than in HM5, so that quantitative differences between the two lines in transcription rates could be the cause of the blue aleurone. The MYB and bHLH transcription factors had two homologous unigenes, but contained only one differentially expressed unigene each. The relative transcript level of bHLH Unigene5672_All (TsMYC2) in HM13 was 42.71 times that in HM5, while the relative transcript level of the MYB transcription factor Unigene12228_All in HM13 was 2.20 times that in HM5. qPCR experiments determined the relative transcript level of TsMYC2 in developing grain, with the expression of TsMYC2 in grain being the highest compared with that in root, stem or leaf tissue. TsMYC2 was homologous to the bHLH transcription factor regulating anthocyanin biosynthesis and contained three entire functional domains: bHLH-MYC_N, HLH and ACT-like, which were important for exercising regulation of anthocyanin biosynthesis as a bHLH transcription factor. Transient expression of ZmC1 and TsMYC2 could induce anthocyanin biosynthesis in white wheat coleoptile cells, demonstrating that TsMYC2 was a functional bHLH transcription factor. These results indicated that TsMYC2 was associated with the blue aleurone trait and could prove to be a valuable gene with which to breed new triticale cultivars with the blue aleurone trait.