Bacterial lifestyle shapes pangenomes.

Pangenomes vary across bacteria. Some species have fluid pangenomes, with a high proportion of genes varying between individual genomes. Other species have less fluid pangenomes, with different genomes tending to contain the same genes. Two main hypotheses have been suggested to explain this variation: differences in species' bacterial lifestyle and effective population size. However, previous studies have not been able to test between these hypotheses because the different features of lifestyle and effective population size are highly correlated with each other, and phylogenetically conserved, making it hard to disentangle their relative importance. We used phylogeny-based analyses, across 126 bacterial species, to tease apart the causal role of different factors. We found that pangenome fluidity was lower in i) host-associated compared with free-living species and ii) host-associated species that are obligately dependent on a host, live inside cells, and are more pathogenic and less motile. In contrast, we found no support for the competing hypothesis that larger effective population sizes lead to more fluid pangenomes. Effective population size appears to correlate with pangenome variation because it is also driven by bacterial lifestyle, rather than because of a causal relationship.

The metabolic domestication syndrome of budding yeast.

Cellular metabolism evolves through changes in the structure and quantitative states of metabolic networks. Here, we explore the evolutionary dynamics of metabolic states by focusing on the collection of metabolite levels, the metabolome, which captures key aspects of cellular physiology. Using a phylogenetic framework, we profiled metabolites in 27 populations of nine budding yeast species, providing a graduated view of metabolic variation across multiple evolutionary time scales. Metabolite levels evolve more rapidly and independently of changes in the metabolic network's structure, providing complementary information to enzyme repertoire. Although metabolome variation accumulates mainly gradually over time, it is profoundly affected by domestication. We found pervasive signatures of convergent evolution in the metabolomes of independently domesticated clades of Saccharomyces cerevisiae. Such recurring metabolite differences between wild and domesticated populations affect a substantial part of the metabolome, including rewiring of the TCA cycle and several amino acids that influence aroma production, likely reflecting adaptation to human niches. Overall, our work reveals previously unrecognized diversity in central metabolism and the pervasive influence of human-driven selection on metabolite levels in yeasts.

Phylogenomic comparative methods: Accurate evolutionary inferences in the presence of gene tree discordance.

Phylogenetic comparative methods have long been a mainstay of evolutionary biology, allowing for the study of trait evolution across species while accounting for their common ancestry. These analyses typically assume a single, bifurcating phylogenetic tree describing the shared history among species. However, modern phylogenomic analyses have shown that genomes are often composed of mosaic histories that can disagree both with the species tree and with each other-so-called discordant gene trees. These gene trees describe shared histories that are not captured by the species tree, and therefore that are unaccounted for in classic comparative approaches. The application of standard comparative methods to species histories containing discordance leads to incorrect inferences about the timing, direction, and rate of evolution. Here, we develop two approaches for incorporating gene tree histories into comparative methods: one that constructs an updated phylogenetic variance-covariance matrix from gene trees, and another that applies Felsenstein's pruning algorithm over a set of gene trees to calculate trait histories and likelihoods. Using simulation, we demonstrate that our approaches generate much more accurate estimates of tree-wide rates of trait evolution than standard methods. We apply our methods to two clades of the wild tomato genus Solanum with varying rates of discordance, demonstrating the contribution of gene tree discordance to variation in a set of floral traits. Our approaches have the potential to be applied to a broad range of classic inference problems in phylogenetics, including ancestral state reconstruction and the inference of lineage-specific rate shifts.

Wiring Between Close Nodes in Molecular Networks Evolves More Quickly Than Between Distant Nodes.

As species diverge, a wide range of evolutionary processes lead to changes in protein-protein interaction (PPI) networks and metabolic networks. The rate at which molecular networks evolve is an important question in evolutionary biology. Previous empirical work has focused on interactomes from model organisms to calculate rewiring rates, but this is limited by the relatively small number of species and sparse nature of network data across species. We present a proxy for variation in network topology: variation in drug-drug interactions (DDIs), obtained by studying drug combinations (DCs) across taxa. Here, we propose the rate at which DDIs change across species as an estimate of the rate at which the underlying molecular network changes as species diverge. We computed the evolutionary rates of DDIs using previously published data from a high-throughput study in gram-negative bacteria. Using phylogenetic comparative methods, we found that DDIs diverge rapidly over short evolutionary time periods, but that divergence saturates over longer time periods. In parallel, we mapped drugs with known targets in PPI and cofunctional networks. We found that the targets of synergistic DDIs are closer in these networks than other types of DCs and that synergistic interactions have a higher evolutionary rate, meaning that nodes that are closer evolve at a faster rate. Future studies of network evolution may use DC data to gain larger-scale perspectives on the details of network evolution within and between species.

Blouch: Bayesian Linear Ornstein-Uhlenbeck Models for Comparative Hypotheses.

Relationships among species in the tree of life can complicate comparative methods and testing adaptive hypotheses. Models based on the Ornstein-Uhlenbeck process permit hypotheses about adaptation to be tested by allowing traits to either evolve towards fixed adaptive optima (e.g., regimes or niches) or track continuously changing optima that can be influenced by other traits. These models allow estimation of the effects of both adaptation and phylogenetic inertia - resistance to adaptation due to any source - on trait evolution, an approach known as the "adaptation-inertia" framework. However, previous applications of this framework, and most approaches suggested to deal with the issue of species non-independence, are based on a maximum likelihood approach and thus it is difficult to include information based on prior biological knowledge in the analysis, which can affect resulting inferences. Here I present Blouch, (Bayesian Linear Ornstein-Uhlenbeck Models for Comparative Hypotheses), which fits allometric and adaptive models of continuous trait evolution in a Bayesian framework based on fixed or continuous predictors and incorporates measurement error. I first briefly discuss the models implemented in Blouch, and then the new applications for these models provided by a Bayesian framework. This includes the advantages of assigning biologically meaningful priors when compared to non-Bayesian approaches, allowing for varying effects (intercepts and slopes), and multilevel modeling. Validations on simulated data show good performance in recovering the true evolutionary parameters for all models. To demonstrate the workflow of Blouch on an empirical dataset, I test the hypothesis that the relatively larger antlers of larger bodied deer are the result of more intense sexual selection that comes along with their tendency to live in larger breeding groups. While results show that larger bodied deer that live in larger breeding groups have relatively larger antlers, deer living in the smallest groups appear to have a different and steeper scaling pattern of antler size to body size than other groups. These results are contrary to previous findings and may argue that a different type of sexual selection or other selective pressures govern optimum antler size in the smallest breeding groups.

Exploring the Macroevolutionary Signature of Asymmetric Inheritance at Speciation.

Popular comparative phylogenetic models such as Brownian Motion, Ornstein-Ulhenbeck, and their extensions, assume that, at speciation, a trait value is inherited identically by two descendant species. This assumption contrasts with models of speciation at a micro-evolutionary scale where descendants' phenotypic distributions are sub-samples of the ancestral distribution. Different speciation mechanisms can lead to a displacement of the ancestral phenotypic mean among descendants and an asymmetric inheritance of the ancestral phenotypic variance. In contrast, even macro-evolutionary models that account for intraspecific variance assume symmetrically conserved inheritance of ancestral phenotypic distribution at speciation. Here we develop an Asymmetric Brownian Motion model (ABM) that relaxes the assumption of symmetric and conserved inheritance of the ancestral distribution at the time of speciation. The ABM jointly models the evolution of both intra- and inter-specific phenotypic variation. It also infers the mode of phenotypic inheritance at speciation, which can range from a symmetric and conserved inheritance, where descendants inherit the ancestral distribution, to an asymmetric and displaced inheritance, where descendants inherit divergent phenotypic means and variances. To demonstrate this model, we analyze the evolution of beak morphology in Darwin finches, finding evidence of displacement at speciation. The ABM model helps to bridge micro- and macro-evolutionary models of trait evolution by providing a more robust framework for testing the effects of ecological speciation, character displacement, and niche partitioning on trait evolution at the macro-evolutionary scale.

Ecological Predictors of Organelle Genome Evolution: Phylogenetic Correlations with Taxonomically Broad, Sparse, Unsystematized Data.

Comparative analysis of variables across phylogenetically linked observations can reveal mechanisms and insights in evolutionary biology. As the taxonomic breadth of the sample of interest increases, challenges of data sparsity, poor phylogenetic resolution, and complicated evolutionary dynamics emerge. Here, we investigate a cross-eukaryotic question where all these problems exist: which organismal ecology features are correlated with gene retention in mitochondrial and chloroplast DNA (organelle DNA or oDNA). Through a wide palette of synthetic control studies, we first characterize the specificity and sensitivity of a collection of parametric and non-parametric phylogenetic comparative approaches to identify relationships in the face of such sparse and awkward datasets. This analysis is not directly focused on oDNA, and so provides generalizable insights into comparative approaches with challenging data. We then combine and curate ecological data coupled to oDNA genome information across eukaryotes, including a new semi-automated approach for gathering data on organismal traits from less systematized open-access resources including encyclopedia articles on species and taxa. The curation process also involved resolving several issues with existing datasets, including enforcing the clade-specificity of several ecological features and fixing incorrect annotations. Combining this unique dataset with our benchmarked comparative approaches, we confirm support for several known links between organismal ecology and organelle gene retention, identify several previously unidentified relationships constituting possible ecological contributors to oDNA genome evolution, and provide support for a recently hypothesized link between environmental demand and oDNA retention. We, with caution, discuss the implications of these findings for organelle evolution and of this pipeline for broad comparative analyses in other fields.

Mutualisms drive plant trait evolution beyond interaction-related traits.

Mutualisms have driven the evolution of extraordinary structures and behavioural traits, but their impact on traits beyond those directly involved in the interaction remains unclear. We addressed this gap using a highly evolutionarily replicated system - epiphytes in the Rubiaceae forming symbioses with ants. We employed models that allow us to test the influence of discrete mutualistic traits on continuous non-mutualistic traits. Our findings are consistent with mutualism shaping the pace of morphological evolution, strength of selection and long-term mean of non-mutualistic traits in function of mutualistic dependency. While specialised and obligate mutualisms are associated with slower trait change, less intimate, facultative and generalist mutualistic interactions - which are the most common - have a greater impact on non-mutualistic trait evolution. These results challenge the prevailing notion that mutualisms solely affect the evolution of interaction-related traits via stabilizing selection and instead demonstrate a broader role for mutualisms in shaping trait evolution.

CAGEE: Computational Analysis of Gene Expression Evolution.

Despite the increasing abundance of whole transcriptome data, few methods are available to analyze global gene expression across phylogenies. Here, we present a new software package (Computational Analysis of Gene Expression Evolution [CAGEE]) for inferring patterns of increases and decreases in gene expression across a phylogenetic tree, as well as the rate at which these changes occur. In contrast to previous methods that treat each gene independently, CAGEE can calculate genome-wide rates of gene expression, along with ancestral states for each gene. The statistical approach developed here makes it possible to infer lineage-specific shifts in rates of evolution across the genome, in addition to possible differences in rates among multiple tissues sampled from the same species. We demonstrate the accuracy and robustness of our method on simulated data and apply it to a data set of ovule gene expression collected from multiple self-compatible and self-incompatible species in the genus Solanum to test hypotheses about the evolutionary forces acting during mating system shifts. These comparisons allow us to highlight the power of CAGEE, demonstrating its utility for use in any empirical system and for the analysis of most morphological traits. Our software is available at https://github.com/hahnlab/CAGEE/.

A Phylogenetic Framework to Simulate Synthetic Interspecies RNA-Seq Data.

Interspecies RNA-Seq datasets are increasingly common, and have the potential to answer new questions about the evolution of gene expression. Single-species differential expression analysis is now a well-studied problem that benefits from sound statistical methods. Extensive reviews on biological or synthetic datasets have provided the community with a clear picture on the relative performances of the available methods in various settings. However, synthetic dataset simulation tools are still missing in the interspecies gene expression context. In this work, we develop and implement a new simulation framework. This tool builds on both the RNA-Seq and the phylogenetic comparative methods literatures to generate realistic count datasets, while taking into account the phylogenetic relationships between the samples. We illustrate the usefulness of this new framework through a targeted simulation study, that reproduces the features of a recently published dataset, containing gene expression data in adult eye tissue across blind and sighted freshwater crayfish species. Using our simulated datasets, we perform a fair comparison of several approaches used for differential expression analysis. This benchmark reveals some of the strengths and weaknesses of both the classical and phylogenetic approaches for interspecies differential expression analysis, and allows for a reanalysis of the crayfish dataset. The tool has been integrated in the R package compcodeR, freely available on Bioconductor.

Phylogenetic Comparative Approach Reveals Evolutionary Conservatism, Ancestral Composition, and Integration of Vertebrate Gut Microbiota.

How host-associated microbial communities evolve as their hosts diversify remains equivocal: how conserved is their composition? What was the composition of ancestral microbiota? Do microbial taxa covary in abundance over millions of years? Multivariate phylogenetic models of trait evolution are key to answering similar questions for complex host phenotypes, yet they are not directly applicable to relative abundances, which usually characterize microbiota. Here, we extend these models in this context, thereby providing a powerful approach for estimating phylosymbiosis (the extent to which closely related host species harbor similar microbiota), ancestral microbiota composition, and integration (evolutionary covariations in bacterial abundances). We apply our model to the gut microbiota of mammals and birds. We find significant phylosymbiosis that is not entirely explained by diet and geographic location, indicating that other evolutionary-conserved traits shape microbiota composition. We identify main shifts in microbiota composition during the evolution of the two groups and infer an ancestral mammalian microbiota consistent with an insectivorous diet. We also find remarkably consistent evolutionary covariations among bacterial orders in mammals and birds. Surprisingly, despite the substantial variability of present-day gut microbiota, some aspects of their composition are conserved over millions of years of host evolutionary history.

The Evolution of Rodent Tail Morphology.

AbstractPopulation-level variation in rodent tail structures has been variously attributed to facilitating social communication, locomotion, thermoregulation, and predator avoidance. Little is known, however, about the applicability of these ecological and social correlates to explaining the tremendous interspecific diversity of this appendage. To investigate the potential drivers of rodent tail morphology at a macroevolutionary level, we first carefully reviewed the literature and constructed a list of major hypotheses regarding this variation. We then compiled a database of 11 different tail traits related to length, color, texture, and ecological characteristics for 2,101 species of rodents (order Rodentia) and examined their key evolutionary correlates. Using Bayesian phylogenetic mixed models across the entire order and additionally within the five rodent suborders, we found that tail length is correlated with both temperature (Allen's rule) and locomotory mode, that black tips are more common in brightly lit environments, that naked tails are often found in warmer climates, that fluffy-tipped tails are more common in smaller and/or arboreal species, that prehensility is predominant in arboreal species and/or species with longer tails, and that tail autotomy is more common in open environments. Most of our tested predictions, largely drawn from population-level studies, are not recapitulated across the entire order, potentially indicating a role of local ecological context in shaping tail morphology.

Genes for cooperation are not more likely to be carried by plasmids.

Cooperation is prevalent across bacteria, but risks being exploited by non-cooperative cheats. Horizontal gene transfer, particularly via plasmids, has been suggested as a mechanism to stabilize cooperation. A key prediction of this hypothesis is that genes which are more likely to be transferred, such as those on plasmids, should be more likely to code for cooperative traits. Testing this prediction requires identifying all genes for cooperation in bacterial genomes. However, previous studies used a method which likely misses some of these genes for cooperation. To solve this, we used a new genomics tool, SOCfinder, which uses three distinct modules to identify all kinds of genes for cooperation. We compared where these genes were located across 4648 genomes from 146 bacterial species. In contrast to the prediction of the hypothesis, we found no evidence that plasmid genes are more likely to code for cooperative traits. Instead, we found the opposite-that genes for cooperation were more likely to be carried on chromosomes. Overall, the vast majority of genes for cooperation are not located on plasmids, suggesting that the more general mechanism of kin selection is sufficient to explain the prevalence of cooperation across bacteria.

Diving into the diversity of colour patterns in reef fishes.

Colours and associated patterns are probably some of the most obvious phenotypic traits in animals and reef teleost fishes are often cited as a textbook example for illustrating this type of diversity. Even if it is well established that colour patterns play a central role in the ecology and evolution of reef fishes, we still lack the necessary toolkits to fully grasp the mechanisms driving the diversification of this obvious phenotypic trait. On the one hand, genotyping power seems now limitless thanks to current DNA sequencing technologies. Today, entire genomes of fishes can be easily produced for large sets of species. On the other hand, the description of colour patterns and the quantification of their variation across reef fishes might be highly challenging. In a cover manuscript in this issue of Molecular Ecology, Coulmance et al. (2023) introduced an innovative approach for extracting and quantifying the major colour pattern elements present in the hamlets (Hypoplectrus spp., Serranidae), a recent reef fish radiation from the Caribbean. Then, they intelligently used the quantified colour pattern variation as a phenotypic trait for a genome-wide association study (GWAS). Interestingly, using a method that required no a priori knowledge, they were able to recover well-established marks (e.g., vertical bars) and to highlight less expected colour pattern elements (e.g., dark to light gradient on ventral part as well as caudal and anal fins), which show strong association peaks on linkage group (LG) 12 and 04. Beyond the demonstration of the potential of their new quantitative analysis of colour pattern variation in reef fishes combined with GWAS, their findings offer new perspectives on our understanding of the intrinsic and extrinsic factors generating this outstanding diversity of the fish world.

Body size evolution in otters distinguished from terrestrial mustelids.

Some taxa of mammals live in water, all of which evolved from land-dwelling ancestors. In the family Mustelidae (Mammalia: Carnivora), most species live on land, while otters, comprising the subfamily Lutrinae, inhabit aquatic environments, which include the almost exclusively aquatic sea otters (Enhydra lutris). Thus, the transition from a terrestrial to an aquatic lifestyle has occurred within this family. Despite potentially different selection pressures on body size in aquatic and terrestrial habitats, no divergence in the evolutionary pattern of body size between otters and other mustelids has previously been shown using models of trait evolution on a phylogeny. We applied models that explicitly incorporated lineage-specific directional selection to the evolution of body mass in living mustelids. Using a simulation-based likelihood and approximate Bayesian computation approach, we demonstrated lineage-specific directional selection for larger body mass in otters, which is distinct from other mustelids. There was no evidence of a difference between sea otters and other otters in the strength of directional selection for larger body mass. Additionally, our analyses supported no difference in the rate at which body mass evolves in both directions between otters and other mustelids. These findings suggest that the evolution of body mass in otters is associated with selective advantages of larger size rather than the relaxation of constraints on body size in aquatic habitats, like other aquatic mammals such as sirenians, cetaceans, and pinnipeds.

The Non-Dereliction in Evolution: Trophic Specialisation Drives Convergence in the Radiation of Red Devil Spiders (Araneae: Dysderidae) in the Canary Islands.

Natural selection plays a key role in deterministic evolution, as clearly illustrated by the multiple cases of repeated evolution of ecomorphological characters observed in adaptive radiations. Unlike most spiders, Dysdera species display a high variability of cheliceral morphologies, which has been suggested to reflect different levels of specialization to feed on isopods. In this study, we integrate geometric morphometrics and experimental trials with a fully resolved phylogeny of the highly diverse endemic species from the Canary Islands to 1) quantitatively delimit the different cheliceral morphotypes present in the archipelago, 2) test their association with trophic specialization, as reported for continental species, 3) reconstruct the evolution of these ecomorphs throughout the diversification of the group, 4) test the hypothesis of convergent evolution of the different morphotypes, and 5) examine whether specialization constitutes a case of evolutionary irreversibility in this group. We show the existence of 9 cheliceral morphotypes and uncovered their significance for trophic ecology. Further, we demonstrate that similar ecomorphs evolved multiple times in the archipelago, providing a novel study system to explain how convergent evolution and irreversibility due to specialization may be combined to shape phenotypic diversification in adaptive radiations.

The Cauchy Process on Phylogenies: A Tractable Model for Pulsed Evolution.

Phylogenetic comparative methods use random processes, such as the Brownian Motion, to model the evolution of continuous traits on phylogenetic trees. Growing evidence for non-gradual evolution motivated the development of complex models, often based on Lévy processes. However, their statistical inference is computationally intensive and currently relies on approximations, high-dimensional sampling, or numerical integration. We consider here the Cauchy Process (CP), a particular pure-jump Lévy process in which the trait increment along each branch follows a centered Cauchy distribution with a dispersion proportional to its length. In this work, we derive an exact algorithm to compute both the joint probability density of the tip trait values of a phylogeny under a CP and the ancestral trait values and branch increments posterior densities in quadratic time. A simulation study shows that the CP generates patterns in comparative data that are distinct from any Gaussian process, and that restricted maximum likelihood parameter estimates and root trait reconstruction are unbiased and accurate for trees with 200 tips or less. The CP has only two parameters but is rich enough to capture complex-pulsed evolution. It can reconstruct posterior ancestral trait distributions that are multimodal, reflecting the uncertainty associated with the inference of the evolutionary history of a trait from extant taxa only. Applied on empirical datasets taken from the Evolutionary Ecology and Virology literature, the CP suggests nuanced scenarios for the body size evolution of Greater Antilles Lizards and for the geographical spread of the West Nile Virus epidemics in North America, both consistent with previous studies using more complex models. The method is efficiently implemented in C with an R interface in package cauphy, which is open source and freely available online.

How Important Is Budding Speciation for Comparative Studies?

The acknowledgment of evolutionary dependence among species has fundamentally changed how we ask biological questions. Phylogenetic models became the standard approach for studies with 3 or more lineages, in particular those using extant species. Most phylogenetic comparative methods (PCMs) translate relatedness into covariance, meaning that evolutionary changes before lineages split should be interpreted together whereas after the split lineages are expected to change independently. This clever realization has shaped decades of research. Here, we discuss one element of the comparative method often ignored or assumed as unimportant: if nodes of a phylogeny represent the dissolution of the ancestral lineage into two new ones or if the ancestral lineage can survive speciation events (i.e., budding). Budding speciation is often reported in paleontological studies, due to the nature of the evidence for budding in the fossil record, but it is surprisingly absent in comparative methods. Here, we show that many PCMs assume that divergence happens as a symmetric split, even if these methods do not explicitly mention this assumption. We discuss the properties of trait evolution models for continuous and discrete traits and their adequacy under a scenario of budding speciation. We discuss the effects of budding speciation under a series of plausible evolutionary scenarios and show when and how these can influence our estimates. We also propose that long-lived lineages that have survived through a series of budding speciation events and given birth to multiple new lineages can produce evolutionary patterns that challenge our intuition about the most parsimonious history of trait changes in a clade. We hope our discussion can help bridge comparative approaches in paleontology and neontology as well as foster awareness about the assumptions we make when we use phylogenetic trees.

Palms of the past: can morphometric phytolith analysis inform deep time evolution and palaeoecology of Arecaceae?

BACKGROUND AND AIMS: Palm fossils are often used as evidence for warm and wet palaeoenvironments, reflecting the affinities of most modern palms. However, several extant palm lineages tolerate cool and/or arid climates, making a clear understanding of the taxonomic composition of ancient palm communities important for reliable palaeoenvironmental inference. However, taxonomically identifiable palm fossils are rare and often confined to specific facies. Although the resolution of taxonomic information they provide remains unclear, phytoliths (microscopic silica bodies) provide a possible solution because of their high preservation potential under conditions where other plant fossils are scarce. We thus evaluate the taxonomic and palaeoenvironmental utility of palm phytoliths. METHODS: We quantified phytolith morphology of 97 modern palm and other monocot species. Using this dataset, we tested the ability of five common discriminant methods to identify nine major palm clades. We then compiled a dataset of species' climate preferences and tested if they were correlated with phytolith morphology using a phylogenetic comparative approach. Finally, we reconstructed palm communities and palaeoenvironmental conditions at six fossil sites. KEY RESULTS: Best-performing models correctly identified phytoliths to their clade of origin only 59 % of the time. Although palms were generally distinguished from non-palms, few palm clades were highly distinct, and phytolith morphology was weakly correlated with species' environmental preferences. Reconstructions at all fossil sites suggested that palm communities were dominated by Trachycarpeae and Areceae, with warm, equable climates and high, potentially seasonal rainfall. However, fossil site reconstructions had high uncertainty and often conflicted with other climate proxies. CONCLUSIONS: While phytolith morphology provides some distinction among palm clades, caution is warranted. Unlike prior spatially restricted studies, our geographically and phylogenetically broad study indicates phytolith morphology may not reliably differentiate most palm taxa in deep time. Nevertheless, it reveals distinct clades, including some likely to be palaeoenvironmentally informative.

Traditional medicinal use is linked with apparency, not specialized metabolite profiles in the order Caryophyllales.

PREMISE: Better understanding of the relationship between plant specialized metabolism and traditional medicine has the potential to aid in bioprospecting and untangling of cross-cultural use patterns. However, given the limited information available for metabolites in most plant species, understanding medicinal use-metabolite relationships can be difficult. The order Caryophyllales has a unique pattern of lineages of tyrosine- or phenylalanine-dominated specialized metabolism, represented by mutually exclusive anthocyanin and betalain pigments, making Caryophyllales a compelling system to explore the relationship between medicine and metabolites by using pigment as a proxy for dominant metabolism. METHODS: We compiled a list of medicinal species in select tyrosine- or phenylalanine-dominant families of Caryophyllales (Nepenthaceae, Polygonaceae, Simmondsiaceae, Microteaceae, Caryophyllaceae, Amaranthaceae, Limeaceae, Molluginaceae, Portulacaceae, Cactaceae, and Nyctaginaceae) by searching scientific literature until no new uses were recovered. We then tested for phylogenetic clustering of uses using a "hot nodes" approach. To test potential non-metabolite drivers of medicinal use, like how often humans encounter a species (apparency), we repeated the analysis using only North American species across the entire order and performed phylogenetic generalized least squares regression (PGLS) with occurrence data from the Global Biodiversity Information Facility (GBIF). RESULTS: We hypothesized families with tyrosine-enriched metabolism would show clustering of different types of medicinal use compared to phenylalanine-enriched metabolism. Instead, wide-ranging, apparent clades in Polygonaceae and Amaranthaceae are overrepresented across nearly all types of medicinal use. CONCLUSIONS: Our results suggest that apparency is a better predictor of medicinal use than metabolism, although metabolism type may still be a contributing factor.