{'fraction_non_alphanumeric': 0.03397190995794739, 'fraction_numerical': 0.014484786851001841, 'mean_word_length': 5.166779661016949, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 3, 'https://': 2, 'lorem ipsum': 0, 'www.': 2, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Phylogeny, classification, and character evolution of tribe Citharexyleae (Verbenaceae)

Premise of the study Tribe Citharexyleae comprises three genera: Baillonia, Citharexylum, and Rehdera. While there is good support for these genera as a clade, relationships between genera remain unresolved due to low sampling of the largest genus, Citharexylum. A molecular phylogenetic approach was taken to resolve intergeneric relationships in Citharexyleae and infrageneric relationships in Citharexylum. Methods Seven chloroplast regions, two nuclear ribosomal spacers, and six low-copy nuclear loci were analyzed for 64 species of Citharexyleae. Phylogenetic analyses were conducted using maximum likelihood, Bayesian inference, and Bayesian multi-species coalescent approaches. Habit, presence/absence of thorns, inflorescence architecture, flower color, fruit color, and geography were examined to identify diagnostic characters for clades within Citharexylum. Key results Intergeneric relationships resolved Rehdera as sister to Citharexylum and Baillonia nested within Citharexylum. Two species, C. oleinum and C. tetramerum, fell outside of Citharexyleae close to tribe Duranteae. There is strong support for seven clades within Citharexylum, each characterized by a unique combination of geography, fruit color/maturation, and inflorescence architecture. Conclusions Baillonia is included in Citharexylum; Rehdera is retained as a distinct genus. A subgeneric classification for Citharexylum is proposed.

uni-carpellate, two-seeded fruit (Briquet, 1895). Baillonia was traditionally classified in tribe Lantaneae, where all unicarpellate members of Verbenaceae were assigned before molecular phylogenies demonstrated the multiple origins of this trait. In addition to carpel number, these genera differ in whether their fruits are fleshy or dry: while Citharexylum and Baillonia both produce fleshy, drupaceous fruits, those of Rehdera are dry and dehiscent. Due to overall similarity in non-fruit traits, the dry-fruited species initially were classified as Citharexylum (Moore, 1895;Junell, 1934), but were segregated into Rehdera by Moldenke (1935).
Estimated to comprise 130 species (Atkins 2004), Citharexylum is by far the largest genus of Citharexyleae, as well as one of the largest genera of Verbenaceae. This would be true even if the suggestions that this estimate is inflated (Sanders, 2001;O'Leary and Frost, in review) are validated. In addition to the fleshy, bi-carpellate fruits that define the genus within the tribe, Citharexylum exhibits cryptic dioecy (Tomlinson and Fawcett, 1972;Rocca and Sazima, 2006;Rueda and Hammel, in prep.), with flowers that contain both stamens and pistils, but one reproductive whorl is sterile so that individuals are functionally staminate or pistillate.
Dioecy is rare in Verbenaceae, known only in Citharexylum and Lippia section Dioicolippia (Múlgura, 2010). It is unknown whether dioecy is universal in Citharexylum. Despite the size of the genus, its taxonomy has been poorly studied and described. Moldenke published a series of notes (1957)(1958)(1959) preparatory to a monograph that was never completed and Walpers (1845) proposed infrageneric groupings that were never validly published. These sections, Spinosa and Inermia, were defined based on presence or absence of spines, respectively; section Inermia was further divided by calyx morphology. Section Spinosa included five taxa that are recognized today as C. flexuosum, C. montevidense, and Rhaphithamnus spinosus (now distinct from Citharexylum); section Inermia encompassed the remaining 29 taxa in the treatment (Walpers, 1845). Moldenke (1957Moldenke ( -1959 mirrored Walpers' sections in a dichotomous key, updating the sections to include the many taxa described (mostly by Moldenke) since the informal treatment by Walpers. Moldenke placed C. andinum, C. herrerae, and C. weberbaueri in section Spinosa, leaving over 100 taxa-the vast majority of species in Citharexylum-listed in section Inermia.
Due to the lack of structure within Inermia, this division provides no hypothesis of relationships between species, which vary in other traits including habit, flower color, inflorescence architecture, and fruit color.
Further phylogenetic study is necessary to fully understand the generic relationships within Citharexyleae and their implications for reclassification. The first phylogeny to include all three genera and sample multiple species of Citharexylum and Rehdera, supported inclusion of Baillonia within Citharexylum and a sister relationship of Rehdera to Citharexylum + Baillonia (Marx et al., 2010). However, because only 11 species in Citharexyleae were sampled and only chloroplast data were used, greater taxonomic sampling (especially of the large genus Citharexylum) and more loci are needed to fully resolve generic relationships. A more recent study expanded taxonomic sampling in Citharexylum and included one nuclear locus (Frost et al., 2017). This provided a framework phylogeny of Citharexylum, but did not test intergeneric relationships in the tribe.
We fill this knowledge gap by significantly expanding taxonomic and genetic sampling in a molecular phylogeny of Citharexyleae. Microfluidic PCR was used to target seven chloroplast regions, two nuclear rDNA spacers, and five low-copy nuclear genes. These data were analyzed in multiple frameworks, including a species-tree approach. With this improved sampling, we (1) resolve relationships between genera of Citharexyleae, (2) identify major clades of Citharexylum and the characters that unite them, (3) infer trait evolution within the tribe, and (4) establish a subgeneric classification for Citharexylum.

MATERIALS AND METHODS
Taxon sampling-One hundred eighty-two accessions representing 67 species of Citharexyleae were included in this study: one accession of monotypic Baillonia, 176 accessions of 64 Citharexylum species, and five accessions of two Rehdera species. Nine species from across Verbenaceae were included as outgroups. Vouchers and GenBank accession numbers may be found in Appendix S1 (see the Supplementary Data with this article). Herbarium acronyms follow Index Herbariorum (Thiers, constantly updated: http://sweetgum.nybg.org/science/ih/).
Primer design and validation followed Uribe-Convers et al. (2016), with the exception that primer sequences were designed by eye, instead of using an R script, from alignments of published data (ncbi.nlm.nih.gov). Each primer pair was designed to amplify 400-600 bp.
Validation reactions were performed with four experimental samples: two Citharexylum species (C. schottii and C. myrianthum) to test performance at the genus level, one taxon from the outgroup (Duranta coriacea) to test performance at the family level, and a negative control.
Sixteen primer pairs were selected for seven chloroplast regions. One primer pair was used for each of the nrDNA and low-copy nuclear loci except for waxy, which was amplified with two primer pairs. Primer pair sequences can be found in Appendix S2 (see the Supplementary Data with this article).

DNA extraction, Microfluidic PCR, and sequencing-Genomic DNA was extracted
from field-collected, silica-gel-dried tissues using a modified CTAB method (Doyle and Doyle, 1987) or from herbarium specimens using a DNeasy plant mini kit (Qiagen, Valencia, California, USA). A total of 7603 bp of cpDNA were targeted for amplification and sequencing, including: Microfluidic PCR was performed in two separate runs on an Access Array System (Fluidigm), one with a 48.48 Access Array integrated fluidic circuit (IFC) and another with a Juno 192.24 IFC. The same primer pairs were submitted for each IFC. Amplicons were harvested and pooled as described in Uribe-Convers et al. (2016). For each IFC, the resulting pools were multiplexed in an Illumina Miseq using the Reagent Kit v3 600 cycles. Microfluidic PCR in the Access Array System, downstream quality control, and Illumina sequencing were performed in the University of Idaho (Moscow, ID) Institute for Bioinformatics and Evolutionary Studies (ibest) Genomic Resources Core facility.
Data processing-Reads from the Illumina Miseq runs were processed using the Fluidigm2PURC pipeline for processing paired-end Illumina data generated using the Fluidigm platform (Blischak et al., 2018). Reads were demultiplexed using the program dbcAmplicons (Uribe-Convers et al., 2016). Paired reads were separated by locus, with the fluidigm2purc script (github.com/pblischak/fluidigm2PURC). Each locus was processed using PURC's purc_recluster.py script (github.com/pblischak/fluidigm2PURC; Edgar, 2004Edgar, , 2010Camacho et al., 2009;Edgar et al., 2011;Martin 2011;Rothfels et al., 2016) to iteratively run chimera detection and sequence clustering to produce a set of putative haplotypes. Reads were processed with clustering thresholds of 0.975, 0.99, 0.995, 0.997 and a size threshold of 5 (i.e. only clusters with ≥ 5 identical reads were retained). The clusters were run through the crunch-clusters script to infer haplotypes (github.com/pblischak/fluidigm2PURC). Chloroplast regions were processed as haploid. Reads from nuclear loci were processed with "unknown" ploidy, meaning all possible haplotypes were kept and haplotype number could vary by sample. Gene trees were generated with RAxML ver. 8.2.11 (Stamatakis, 2014) to process the haplotypes for nuclear loci. For individuals with haplotypes that nested together in gene trees (e.g., haplotypes that exhibited sequence similarity, but differed in amplicon length or a few heterogeneous sites) either the longest read or largest cluster was selected to represent the individual. For individuals with haplotypes that nested in different clades throughout the tree, all haplotypes were kept.
Raw reads for each locus were parsed into separate files by the purc_recluster.py script.
In the instances in which raw data were available, but with insufficient coverage to meet the filters set in purc_recluster.py, the raw data were assessed and incorporated manually. A consensus sequence was generated from the available reads for regions with a single haplotype.
For regions with multiple haplotypes, reads were analyzed individually in a phylogenetic context using RAxML ver. 8.2.11 (Stamatakis, 2014); consensus sequences were generated for reads that formed clades.
Four datasets were assembled for each region with different tiers of quality: 1) the original output from data processing, 2) the former plus individuals/haplotypes with ≥ 5 raw reads, 3) the former plus individuals/haplotypes with 2-4 raw reads, and 4) the former plus individuals/haplotypes with only a single read. Phylogenetic analyses were performed for each tier with RAxML ver. 8.2.11 (Stamatakis, 2014) to assess effects of the supplemented data on support values and topology. This supplementation technique was used only for individuals that had sufficient data to be included in the study, but had missing data for some loci due to low coverage. No additional individuals were added to the study as a result of these low-coverage reads.

Phylogenetic analyses
Gene tree analyses-Sequences were aligned using AliView ver. 1.18.1 (Larsson, 2014); alignments were inspected and manually adjusted where necessary. The best-fit model of nucleotide substitution for each alignment was determined with jModelTest 2.1.7 (Darriba et al., 2012;Guindon and Gascuel, 2003); the GTR + G model was selected and applied to all alignments. Loci were concatenated using Geneious R7 7.0.6 (Biomatters, Auckland, New Zealand) and treated as a single dataset for regions with shared evolutionary histories: ndhF, rpoC2, trnT-trnL-trnF, rbcL, matK, and rpl32 are part of the haploid chloroplast genome (cpDNA), and ETS and ITS are tightly linked in the nuclear rDNA repeat (nrDNA). Each lowcopy nuclear gene was treated as an independent locus. When multiple haplotypes were recovered for a region, haplotypes were separated and treated as independent loci. Eight independent loci were analyzed: cp, nr, PPR24, PPR62, PPR70, PPR123, waxy haplotype1, and waxy haplotype 2.
Phylogenetic analyses were conducted for individual datasets, low-copy nuclear genes combined, and all data combined using maximum likelihood and Bayesian approaches, as implemented in RAxML and MrBayes ver. 3.2.6 ( Ronquist & Huelsenbeck, 2003), respectively.
Maximum likelihood analyses consisted of bootstrap analyses with 500 replicates for each dataset. Datasets were partitioned into individual gene regions with individual parameters unlinked in Bayesian analyses. Analyses consisted of two replicate runs, each with four chains, which were run for at least one million generations, sampled every 200 generations. Runs were assessed for convergence and stationarity in Tracer v1.6.0 (Rambaut et al., 2014); our benchmark for assessing convergence were effective sample sizes (ESS) above 200 for all parameter values across combined files. The first 20% of sampled trees were discarded as burnin from each run; the remaining trees were summarized into a majority rule consensus tree with all compatible groups.
Species tree analyses-Species tree estimation was performed on the ingroup using the multispecies coalescent (MSC) model implemented in *BEAST (BEAST v1.8.4;Drummond et al. 2012). Substitution models and clock models were unlinked for the fourteen individual alignments. The GTR+G substitution model, empirical base frequencies, and a random local clock (Drummond and Suchard, 2010) were applied to each. Trees were linked based on shared evolutionary history, in the case of cp and nr datasets, or based on gene tree topology for lowcopy nuclear genes. Two partitions were set for low-copy nuclear genes: (1) those where Mesoamerican species formed a clade (PPR24, PPR123, waxy haplotype 2) and (2) those where Mesoamerican species were paraphyletic with respect to South American species (PPR62, PPR70, waxy haplotype 1). Chloroplast data were treated as an organellar (haploid) locus with half the effective population size of a bi-parentally inherited locus, and the nuclear loci were treated as autosomal nuclear. The birth-death process was used as the species tree prior. Two replicate runs were performed for 500 million generations, sampling every 25,000. Convergence between runs and stationarity within were assessed in Tracer v1.6.0 (Rambaut et al., 2014); the duplicate runs were considered to have converged when effective sample sizes (ESS) for parameter values were above 200 for combined files. The first 10% of sampled tree states were discarded from each run as burn-in. Accessions of C. bourgeauianum, C. flabellifolium, C.
guatemalense, and C. reticulatum were excluded from the species tree analyses due a high proportion of missing data. Additionally, Citharexylum oleinum and C. tetramerum were excluded because they were inferred to belong outside of the ingroup.

Morphology-In order to understand character evolution and describe clades in
Citharexylum, morphological characters were examined. Targeted characters included habit, presence/absence of thorns, inflorescence architecture, flower number, corolla color, fruit color, and geography. These traits, which are easily observable in the field, from photographs, and from herbarium specimens, were chosen because they are practical, likely informative in recognizing large clades, and have the potential to fill gaps in natural history knowledge. For habit, species were scored as shrubs (i.e., less than six m tall and with several main stems) or trees (i.e., a single trunk and/or taller than 6 m). To best describe variation in inflorescence architecture terminal and axillary inflorescences were scored separately; these were scored as absent, simple, or compound. Average flower number per inflorescence was categorized as >20, 10-20, or <10.
Corolla color was categorized as solid white, solid yellow, solid purple, or white/yellow with purple/orange on the abaxial surface of the petal lobes. Fruit color/maturation was categorized as red synchronous (green to red), black synchronous (green to black), black sequential (green to red to black) in Citharexylum. Because fruit type does not change within Citharexylum and fruit color and mode of maturation in Citharexylum was the focus of this analysis, the dry fruits of sister genus Rehdera were coded as brown. Geographic range was divided into three regions: Mesoamerica, Caribbean, and South America. Character data were pooled from Moldenke (1957)(1958)(1959); field observations; herbarium records via personal visits to A/GH, F, MEXU, MO, and USM; and online databases GBIF (https://www.gbif.org/), CV Starr Virtual Herbarium (NYBG; sweetgum.nybg.org/science/vh/), SEINET (http://swbiodiversity.org/seinet/), and Tropicos (MOBOT; https://www.tropicos.org/).

Ancestral State Reconstructions -Ancestral character reconstructions were performed
for all traits except geography using parsimony and maximum likelihood estimation. Parsimony reconstructions were performed using the "Trace Character History" > "Parsimony Ancestral States" function in Mesquite v3.51 (Maddison and Maddison, 2018). Maximum likelihood estimations were conducted using the R package phytools (Revell, 2012) using the ace function and an equal rates (ER) model. The biogeography of Citharexylum will be addressed in a separate, more detailed study using more appropriate methods for biogeographic reconstruction.
Because the presence and/or complexity of terminal and/or axillary inflorescences varied within species, each inflorescence type was simply scored as present or absent for ancestral character reconstructions.
Tests for correlated evolution-The BayesDiscrete module in BayesTraits (Pagel and Meade, 2006) was used to test for correlation between morphological characters. Maximum likelihood analyses were conducted for each possible pair of characters coded as binary variables. Traits that were coded as multi-state in ancestral state reconstructions-flower number, flower color and fruit color/maturation-were simplified to binary variables or excluded. Flower number was scored as ≤ 20 flowers per inflorescence vs. ≥ 20 flowers per inflorescence. Flower color was scored pale or purple, including species that have purple lobes.
Fruit maturation was simplified to synchronous or sequential. Fruit color was not included in correlation tests due to the difficulty of scoring red, black, and red-to-black fruits in a binary and biologically accurate way.
For each pair of characters, two maximum likelihood models were fit to the phylogeny in BayesTraits: an independent model (I) and a dependent model (D). The independent model assessed the likelihood of the data under the null hypothesis that the traits evolved independently, whereas the dependent model assumed the transition rates for one character were dependent upon the state of the other character. To determine if traits were correlated, the likelihood estimate of the independent model was compared to the likelihood estimate of the dependent model via a likelihood ratio test. The likelihood ratio was tested against a χ 2 distribution, with four degrees of freedom (following Pagel 1994). Probabilities were calculated from χ 2 values.

Dataset characteristics-
The total concatenated dataset consisted of 12927 aligned bp and 176 accessions, including 167 accessions in Citharexyleae. Among the characters, 4140 were variable, of which 2398 were potentially parsimony-informative. One hundred thirty-nine accessions had sequence data for at least five of the eight loci, with 82 of these having data for at least seven loci. All accessions included in the concatenated dataset had data for both cpDNA and nrDNA. Nine accessions of Citharexylum with insufficient data to be included in the concatenated analyses, but representing species not otherwise included, were retained in individual gene tree analyses. Characteristics of individual loci are shown in Table 1.

Phylogenetic analyses
Gene tree analyses-All gene trees support the monophyly of Citharexyleae, except for two species (C. oleinum and C. tetramerum), which nest among the outgroups. As in previous studies (Marx et al, 2010), Rehdera is inferred to be sister to Citharexylum, while Baillonia is Clade A includes clades I, II, and III in the chloroplast, low-copy nuclear combined, and all data combined trees; the placement of clades II and III varied in individual gene trees. In all trees, clade B included clades IV, V, and VI (species in the latter often form a grade in individual gene trees), though relationships among taxa in clade VI differ among datasets (Appendices S3-S12).
Due to the geographic distribution of most of their species, clade A (including clades I, II, and III, as supported by combined analyses) will be referred to as the Mesoamerican clade and clade B as the South American clade.
Citharexylum spinosum, the type species for the genus, nested within different clades for different datasets. All three accessions nested within the Mesoamerican clade in cpDNA, PPR70, and PPR123 trees. Accessions were split, with one or more nesting in the Mesoamerican clade and one or more in the South American clade in nrDNA, PPR24, and PPR62 trees. Two accessions of C. spinosum had an additional copy for each of the waxy of haplotypes (i.e. four total copies of the waxy gene, two copies of haplotype 1 and two copies of haplotype 2). For each haplotype, one copy nested within the Mesoamerican clade and another nested within the South American clade (Appendices S9, S10). One accession also showed this pattern in PPR62 (C. spinosum_JCS5727; Appendix S6).

Species tree analyses-
The species tree topology was consistent with the topology of the all-combined dataset (Fig. 4). The split between the Mesoamerican clade and the South American clade was supported in the species tree. The same seven clades were inferred in the species tree as in the individual and combined gene trees with high support. Relationships within subclade VI varied between replicate runs in the multispecies coalescent and differed from relationships recovered by the all combined dataset. Three clades within clade VI were consistently inferred in the species tree and all combined analyses: a clade of C. flexuosum, C. kobuskianum, C. peruvianum, and C. weberbaueri, a clade of C. andinum, C. joergensenii, and C. montevidense, and a clade including C. argutedentatum, C. herrerae, and C. sulcatum.

Morphology-Observations
were available for all characters for all species except fruit color in C. racemosum (Fig. 4). Variation within species was observed for traits related to inflorescence architecture. Terminal inflorescences were present in all species, but often absent in individuals of C. subflavescens and C. montanum. Terminal inflorescence complexity (i.e., simple or compound) varied within all species of Citharexylum for which compound terminal inflorescences were observed; only Rehdera penninervia was found to consistently have compound terminal inflorescences. Similarly, there was variation within species for presence/absence of axillary inflorescences. Individuals of 27 of the species-nearly half the taxonomic diversity of the tribe-had axillary inflorescences. Of those 27, 12 species were observed to consistently have axillary inflorescences, and six included individuals with compound axillary inflorescences. As with the terminal infloresence, Rehdera penninervia was the only species for which members consistently had compound axillary inflorescences.
Ancestral State Reconstructions-Habit type shifted repeatedly across the phylogeny; and proved difficult to reconstruct at deeper nodes (Appendix S13). Terminal inflorescences, which were present in all species but often absent in C. montanum and C. subflavescens, either became labile in the most recent common ancestor of C. montanum, C. subflavescens, and C. rimbachii and that lability was lost in C. rimbachii or lability was independently gained in the two species (Appendix S14). Axillary inflorescences were lost in the ancestor of the two major clades of Citharexylum and regained independently nine times (Appendix S15). Reduction of flower number per inflorescence from >20 to <20 was estimated to have occurred in C. altamiranum, C.
alanii and C. microphyllum, clade II, and clade VI (Appendix S16). Flowers with at least partial purple coloration arose six times throughout the evolutionary history genus (Appendix S17).
Reconstructions of the evolutionary history of fruit color were inconsistent. Red synchronous fruits were estimated to be the ancestral fruit color in Citharexylum, with black synchronous fruits evolving once in C. altamiranum and one to three times in clade VI of the South American clade by maximum likelihood analyses (Fig. 5, Appendix S18A). Maximum parsimony reconstructions, however, suggested that black synchronous fruits were the ancestral state for the genus and the South American clade, with red synchronous fruits arising three times in clade VI (Appendix S18B). Likewise, maximum likelihood reconstructions had higher support for multiple origins of black sequential fruits in the Mesoamerican clade, whereas maximum parsimony reconstructions supported a single origin and multiple losses (Appendix S18).
The variable topologies of subclade VI have implications for character evolution of thorns (Appendix S19). Thorns are exclusively found in members of subclade VI, though not all species bear thorns. There was variable support for a single evolution of thorns in the ancestor of subclade VI followed by multiple losses (sometimes with secondary regains), or multiple independence evolutions of thorns (Appendix S20). Support for these reconstructions approached equivocal.

Intergeneric relationships-Phylogenetic reconstructions infer Rehdera to be sister to
Citharexylum and support the inclusion of Baillonia in Citharexylum. This is consistent with previous studies of Verbenaceae and is bolstered by increased sampling of species in Citharexylum. Two species of Citharexylum, C. oleinum and C. tetramerum, are found to occur among the outgroup (Fig. 3). Preliminary results suggest these species may belong to tribe Duranteae.
One additional species in subgenus Southamericanum is tentatively placed in section Sylvaticum.
Citharexylum reticulatum had too much missing data to include the MSC analysis, but was strongly supported as being related to the other members of section Sylvacticum. This species would be difficult to place by morphology alone: it has the long, many-flowered inflorescences with white flowers that are typical of section Sylvaticum, but differs in many ways, including its small stature, reportedly black fruits (Moldenke, 1957(Moldenke, -1959, and its distribution at higher elevations, all of which would be consistent with its placement in section Armatum, the Andean shrub clade. Section Armatum is a particularly difficult group to parse further taxonomically. Its species all have synchronous fruit maturation, but the ripe fruit color can be either red or black, and fruit color does not define lineages within this section. Similarly, many species of this section also produce axillary shoots that are modified into thorns-the first character applied in . This specific hypothesis is not supported in the exclusively tropical Citharexylum, nor does there seem to be a connection between seasonality in dry habitats and the loss of sequential maturation in subgenus Citharexylum. Fruit color could also be connected to the degree to which a plant species is able to invest in costly anthocyanins; darker colors of fleshy fruits are perceived as an honest signal of food quality in birds, and thus would be beneficial in conditions that allow this costly investment (Schaefer et al., 2008).

Flower color evolution-Purple floral pigmentation has arisen at least four times within
Citharexylum (Fig. S6). This includes in the common ancestor of section Pluriflorum and subgenus Purpuratum, twice within each section Citharexlyum and Armatum, in which the flower buds of several species are tipped with purple. These species also have dark purple to black fruits. Though not all black-fruited species have pigmented flowers, there may be some connection between fruit color and expression of anthocyanins in the flowers, as has been found in some plant communities (Renoult et al., 2014). Species with black fruits also occur at high elevations in the Andes, where UV radiation is higher (Steyn, 2008). In these species, darker pigments in flowers and fruits might offer some protection to developing reproductive structures and seeds by absorbing excess UV radiation. This hypothesis does not apply, however, to the purple flowered species in Mesoamerica, which do not occur at high elevations. Flower color in these species may be linked to other factors such as pollination, but little is known about this aspect of their natural history. However, the large, showy corollas of some purple flowers (e.g.,

Infloresence evolution and dioecy-Dioecy needs further investigation throughout
Cithareyxlum. Several species have been confirmed as dioecious in the literature, including C.
spinosum (Tomlinson and Fawcett, 1972), C. myrianthum (Rocca and Sazima, 2006), and all Costa Rican species (Rueda and Hammil, in prep.). Dioecy has been observed in the field for many others (L. Frost, personal observation); it can be inferred in Citharexylum from observation of post-anthesis inflorescences. In staminate individuals, the entire flower disarticulates from the pedicel at senescence, leaving naked racemes. In pistillate individuals, the calyx becomes lignified as the fruits mature and is retained on the pedicel, even as fruits are eaten or fall off. individuals could increase flower production either through the development of axillary inflorescences, or compound inflorescences, or both. Variation between sexes is easily observed for species with large, many-flowered inflorescences, but is difficult to detect in shrub species with reduced inflorescences. These species have simple, few-flowered inflorescences that terminate short, axillary branches. One individual does not obviously produce more flowers or more inflorescences than another. It is also more difficult in general to infer the sex of individuals in these few-flowered species using the observations described above and dioecy has not been postulated for any of these species. Further investigation of dioecy in the Andean and arid-adapted Mexican shrubs is warranted.
Cryptic dieocy is not always readily observable and may be more widespread in the in Verbenaceae is in Lippia section Dioicolippia (Múlgura, 2010), although gynodioecy also has been documented in Rhaphithamnus venustus (Crawford et al., 1993).

Evolution of dry-adaptive traits-Of the morphological traits examined, habit, axillary
inflorescences, and flower number were found to be significantly correlated with each other. In Citharexylum combinations of these traits were typically found in shrubs with simple, exclusively terminal inflorescences (i.e., no axillary inflorescences) and few flowers per inflorescence, or trees with many-flowered inflorescences that often include axillary inflorescences, compound inflorescences or both. These alternate suites of traits are found repeatedly across the phylogeny of Citharexylum. For example, shrubs with few-flowered, simple, terminal racemes occur in at least four distantly related lineages: C. roxanae, C.
flabellifolium (phylogenetic placement uncertain), section Armatum, section Caribe, and section Mexicanum. While not closely related, these clades share similar environmental constraints: limited water availability in their high elevation and/or arid to semi-arid biomes. These environmental conditions are known to constrain above-ground growth (Shantz, 1927), and woody plants are often reduced to shrubs in these habitats. Plants growing at high elevations may be additionally subject to high winds, cloud cover, and low temperatures, which can also limit the height and reproductive effort of woody species (Howard, 1970). We suggest that this combination of characters-short habit and reduced investment in reproductive growth-also reflects a high/dry habitat syndrome in Citharexylum. herrerae, which is more closely related to the high-elevation grassland/cloud forest species that lack thorns. Though not technically bearing thorns, some species outside of the inter-Andean dry valley shrubs also have sharp defenses. Citharexylum amabilis has small, sharp projections borne out of leaf scars (referred to a "spiny sterigmata" by Moldenke, 1941), and many arid-adapted Mexican shrubs (C. brachyanthum, C. flabeillifolium, C. lycioides, C. racemosum, and C. rosei) have relatively short, pointed, branches (Moldenke 1957(Moldenke -1959.
The evolution of distinctive leaf traits is also associated with the dry-adaptive syndrome in Citharexylum. Though they share some traits related to adaptation to arid environments C.
roxanae, C. flabellifolium, and section Mexicanum have evolved different vegetative strategies for coping with life in the desert. Whereas many Andean species of section Armatus have tough, leathery leaves along with their thorns, Citharexylum roxanae has chartaceous, ephemeral leaves and an ephedroid habit in which the persistently green stems likely contribute to photosynthesis when most of the leaves have fallen. Cithareyxlum flabellifolium, which occurs with C. roxanae in Baja California Sur, but extends further into Sonora, has fan-shaped leaves that are densely hairy, presumably protecting leaf tissue from the sun and maintaining humidity around stomatal pores (Ehleringer, 1982). Finally, the four species of section Mexicanum found in the Chihuahuan desert have small, narrow leaves that dissipate heat better than the large, broad leaves that characterize other subclades, reducing the need for evaporative cooling and water loss from the process (Solbrig and Orians, 1977 Future work is aimed at understanding patterns of movement and evolution into (and out of) biomes that led to the present distribution of Citharexylum and likely drove the repeated evolution of suites of environmentally-relevant characters.

Conclusions
We present the first phylogeny of Citharexyleae to include near complete taxonomic sampling, confirming earlier hypotheses of intergeneric relationships. The first robust phylogeny of Citharexylum provides a framework to identify major clades, map morphological characters which define those clades, and understand patterns of character evolution. Three subgenera are described based on fruit maturation and geography. Six sections are described based on habit and inflorescence architecture, some by flower color and fruit color as well. Morphology at lower ranks, sections and subclades often corresponds with environment. Citharexylum is a widespread in the Neotropics and inhabits multiple biomes. We demonstrate a dynamic pattern of trait evolution, including the correlated evolution of multiple traits associated with colonization of dry habitats. This labile trait evolution likely facilitated the rich biogeographic history of the genus, reflected here in the morphological diversity. Future studies will examine in detail the biogeographic patterns underlying diversity in Citharexylum, as well as continue to clarify evolutionary relationships and refine taxonomic limits in parts of Citharexylum characterized by rapid radiations.

DATA AVAILABILITY STATEMENT
Submission of Illumina MiSeq reads to the NCBI Sequence Read Archive (SRA) is in progress under submission number SUB7989566 and will be made available upon publication. All other data formats (tree files, alignments, character matrices, etc.) will be uploaded on dryad and made available upon publication.

SUPPORTING INFORMATION
Additional Supporting Information may be found online in the supporting information section at the end of the article.  A  C  A  C  T  G  A  C  G  A  C  A  T  G  G  T  T  C  T  A  C  A   T  a  r  g  e  t  N  a  m  e  T  a  r  g  e  t  -S  p  e  c  i  f  i  c  S  e  q  u  e  n  c  e  (  5  '  -3  '  )  O  l  i  g  o  n  u  c  l  e  o  t  i  d  e  S  e  q  u  e  n  c  e  f  o  r  m  i  c  r  o  f  l