Chemical Composition of Essential Oils from Natural Populations of Artemisia scoparia Collected at Different Altitudes: Antibacterial, Mosquito Repellent, and Larvicidal Effects.

The current study aimed to evaluate the presence of chemical variations in essential oils (EOs) extracted from Artemisia scoparia growing at different altitudes and to reveal their antibacterial, mosquito larvicidal, and repellent activity. The gas chromatographic-mass spectrometric analysis of A. scoparia EOs revealed that the major compounds were capillene (9.6-31.8%), methyleugenol (0.2-26.6%), ß-myrcene (1.9-21.4%), γ-terpinene (1.5-19.4%), trans-ß-caryophyllene (0.8-12.4%), and eugenol (0.1-9.1%). The EO of A. scoparia collected from the city of Attock at low elevation was the most active against Staphylococcus aureus, Bacillus subtilis, Escherichia coli, and Pseudomonas aeruginosa bacteria (minimum inhibitory concentration of 156-1250 µg/mL) and showed the best mosquito larvicidal activity (LC50, 55.3 mg/L). The EOs of A. scoparia collected from the high-altitude areas of Abbottabad and Swat were the most repellent for females of Ae. aegypti and exhibited repellency for 120 min and 165 min, respectively. The results of the study reveal that different climatic conditions and altitudes have significant effects on the chemical compositions and the biological activity of essential oils extracted from the same species.

Phylogenomics reveals patterns of ancient hybridization and differential diversification that contribute to phylogenetic conflict in willows, poplars, and close relatives.

Despite the economic, ecological, and scientific importance of the genera Salix L. (willows) and Populus L. (poplars, cottonwoods, and aspens) Salicaceae, we know little about the sources of differences in species diversity between the genera and of the phylogenetic conflict that often confounds estimating phylogenetic trees. Salix subgenera and sections, in particular, have been difficult to classify, with one recent attempt termed a "spectacular failure" due to a speculated radiation of the subgenera Vetrix and Chamaetia. Here, we use targeted sequence capture to understand the evolutionary history of this portion of the Salicaceae plant family. Our phylogenetic hypothesis was based on 787 gene regions and identified extensive phylogenetic conflict among genes. Our analysis supported some previously described subgeneric relationships and confirmed the polyphyly of others. Using an fbranch analysis, we identified several cases of hybridization in deep branches of the phylogeny, which likely contributed to discordance among gene trees. In addition, we identified a rapid increase in diversification rate near the origination of the Vetrix-Chamaetia clade in Salix. This region of the tree coincided with several nodes that lacked strong statistical support, indicating a possible increase in incomplete lineage sorting due to rapid diversification. The extraordinary level of both recent and ancient hybridization in both Salix and Populus have played important roles in the diversification and diversity in these two genera.

The chemical biogeography of a widespread aromatic plant species shows both spatial and temporal variation.

Plants produce a wide variety of secondary metabolites, but intraspecific variation in space and time can alter the ecological interactions these compounds mediate. The aim of this work was to document the spatial and temporal chemical biogeography of Monarda fistulosa. I collected leaves from 1587 M. fistulosa individuals from 86 populations from Colorado to Manitoba, extracted and analyzed their terpenes with gas chromatography, mapped monoterpene chemotypes, and analyzed chemical variation with principal component analysis. I also measured the amounts of terpenes in different plant tissues to examine intraplant variation and monitored leaf terpene chemistry over a single growing season to examine temporal patterns. Finally, I extracted terpenes from herbarium samples up to 125 years old and compared the chemotypes with recent samples from the same sites. M. fistulosa populations consisted mostly of thymol (T) and carvacrol (C) chemotypes, but geraniol (G) and (R)-(-)-linalool (L), a chemotype new to this species, were also present. A principal component analysis showed that most of the chemical variation across populations was due to the amounts of the dominant terpene in plants. Intraplant tissue chemistry revealed that leaves mostly had the greatest amounts of terpenes, followed by floral structures, stems, and roots. Short-term temporal variation in leaf chemistry of T and C plants over a growing season showed that plants produced the highest levels of biosynthetic precursors early in the season and their dominant monoterpenes peaked in mid-summer. Plant chemotype was discernable in the oldest herbarium samples, and 15 of 18 historic samples matched the majority chemotype currently at the site, indicating that population chemotype ratios may remain stable over longer time scales. Overall, the results show that plant species' secondary chemistry can vary both spatially and temporally, which may alter the biological interactions that these compounds mediate over space and time.

Repeated turnovers keep sex chromosomes young in willows.

BACKGROUND: Salicaceae species have diverse sex determination systems and frequent sex chromosome turnovers. However, compared with poplars, the diversity of sex determination in willows is poorly understood, and little is known about the evolutionary forces driving their turnover. Here, we characterized the sex determination in two Salix species, S. chaenomeloides and S. arbutifolia, which have an XY system on chromosome 7 and 15, respectively. RESULTS: Based on the assemblies of their sex determination regions, we found that the sex determination mechanism of willows may have underlying similarities with poplars, both involving intact and/or partial homologs of a type A cytokinin response regulator (RR) gene. Comparative analyses suggested that at least two sex turnover events have occurred in Salix, one preserving the ancestral pattern of male heterogamety, and the other changing heterogametic sex from XY to ZW, which could be partly explained by the "deleterious mutation load" and "sexually antagonistic selection" theoretical models. We hypothesize that these repeated turnovers keep sex chromosomes of willow species in a perpetually young state, leading to limited degeneration. CONCLUSIONS: Our findings further improve the evolutionary trajectory of sex chromosomes in Salicaceae species, explore the evolutionary forces driving the repeated turnovers of their sex chromosomes, and provide a valuable reference for the study of sex chromosomes in other species.

Genetic mapping of sexually dimorphic volatile and non-volatile floral secondary chemistry of a dioecious willow.

Secondary chemistry often differs between sexes in dioecious plant species, a pattern attributed to its possible role in the evolution and/or maintenance of dioecy. We used GC-MS to measure floral volatiles emitted from, and LC-MS to quantitate non-volatile secondary compounds contained in, female and male Salix purpurea willow catkins from an F2 family. Using the abundance of these chemicals, we then performed quantitative trait locus (QTL) mapping to locate them on the genome, identified biosynthetic candidate genes in the QTL intervals, and examined expression patterns of candidate genes using RNA-seq. Male flowers emitted more total terpenoids than females, but females produced more benzenoids. Male tissue contained greater amounts of phenolic glycosides, but females had more chalcones and flavonoids. A flavonoid pigment and a spermidine derivative were found only in males. Male catkins were almost twice the mass of females. Forty-two QTL were mapped for 25 chemical traits and catkin mass across 16 of the 19 S. purpurea chromosomes. Several candidate genes were identified, including a chalcone isomerase associated with seven compounds. A better understanding of the genetic basis of the sexually dimorphic chemistry of a dioecious species may shed light on how chemically mediated ecological interactions may have helped in the evolution and maintenance of dioecy.

Population-specific responses of floral volatiles to abiotic factors in changing environments.

PREMISE: Shifts in abiotic factors can affect many plant traits, including floral volatiles. This study examined the response of floral volatiles to water availability and whether phenotypic plasticity to water availability differs among populations. It also investigated genetic differentiation in floral volatiles, determined the effect of temperature on phenotypic plasticity to water availability, and assessed temporal variation in floral scent emission between day and evening, since pollinator visitation differs at those times. METHODS: Rocky Mountain columbine plants (Aquilegia coerulea), started from seeds collected in three wild populations in Colorado, Utah, and Arizona, were grown under two water treatments in a greenhouse in Madison, Wisconsin, United States. One population was also grown under the two water treatments, at two temperatures. Air samples were collected from enclosed flowers using dynamic headspace methods and floral volatiles were identified and quantified by gas chromatography (GC) with mass spectrometry (MS). RESULTS: Emission of three floral volatiles increased in the wetter environment, indicating phenotypic plasticity. The response of six floral volatiles to water availability differed among populations, suggesting genetic differentiation in phenotypic plasticity. Five floral volatiles varied among populations, and emission of most floral volatiles was greater during the day. CONCLUSIONS: Phenotypic plasticity to water availability permits a quick response of floral volatiles in changing environments. The genetic differentiation in phenotypic plasticity suggests that phenotypic plasticity can evolve but complicates predictions of the effects of environmental changes on a plant and its pollinators.

Combined drought and bark beetle attacks deplete non-structural carbohydrates and promote death of mature pine trees.

How carbohydrate reserves in conifers respond to drought and bark beetle attacks are poorly understood. We investigated changes in carbohydrate reserves and carbon-dependent diterpene defences in ponderosa pine trees that were experimentally subjected to two levels of drought stress (via root trenching) and two types of biotic challenge treatments (pheromone-induced bark beetle attacks or inoculations with crushed beetles that include beetle-associated fungi) for two consecutive years. Our results showed that trenching did not influence carbohydrates, whereas both biotic challenges reduced amounts of starch and sugars of trees. However, only the combined trenched-bark beetle attacked trees depleted carbohydrates and died during the first year of attacks. While live trees contained higher carbohydrates than dying trees, amounts of constitutive and induced diterpenes produced did not vary between live and beetle-attacked dying trees, respectively. Based on these results we propose that reallocation of carbohydrates to diterpenes during the early stages of beetle attacks is limited in drought-stricken trees, and that the combination of biotic and abiotic stress leads to tree death. The process of tree death is subsequently aggravated by beetle girdling of phloem, occlusion of vascular tissue by bark beetle-vectored fungi, and potential exploitation of host carbohydrates by bark beetle symbionts as nutrients.

Growth-defense trade-offs shape population genetic composition in an iconic forest tree species.

All organisms experience fundamental conflicts between divergent metabolic processes. In plants, a pivotal conflict occurs between allocation to growth, which accelerates resource acquisition, and to defense, which protects existing tissue against herbivory. Trade-offs between growth and defense traits are not universally observed, and a central prediction of plant evolutionary ecology is that context-dependence of these trade-offs contributes to the maintenance of intraspecific variation in defense [Züst and Agrawal, Annu. Rev. Plant Biol., 68, 513-534 (2017)]. This prediction has rarely been tested, however, and the evolutionary consequences of growth-defense trade-offs in different environments are poorly understood, especially in long-lived species [Cipollini et al., Annual Plant Reviews (Wiley, 2014), pp. 263-307]. Here we show that intraspecific trait trade-offs, even when fixed across divergent environments, interact with competition to drive natural selection of tree genotypes corresponding to their growth-defense phenotypes. Our results show that a functional trait trade-off, when coupled with environmental variation, causes real-time divergence in the genetic architecture of tree populations in an experimental setting. Specifically, competitive selection for faster growth resulted in dominance by fast-growing tree genotypes that were poorly defended against natural enemies. This outcome is a signature example of eco-evolutionary dynamics: Competitive interactions affected microevolutionary trajectories on a timescale relevant to subsequent ecological interactions [Brunner et al., Funct. Ecol. 33, 7-12 (2019)]. Eco-evolutionary drivers of tree growth and defense are thus critical to stand-level trait variation, which structures communities and ecosystems over expansive spatiotemporal scales.

Characterization of Salix nigra floral insect community and activity of three native Andrena bees.

Salix nigra (black willow) is a widespread tree that hosts many species of polylectic hymenopterans and oligolectic bees of the genus Andrena. The early flowering of S. nigra makes it an important nutritive resource for arthropods emerging from hibernation. However, since S. nigra is dioecious, not all insect visits will lead to successful pollination. Using both visual observation and pan-trapping, we characterized the community of arthropods that visited S. nigra flowers and assessed differences among male and female trees as well as the chemical and visual drivers that influenced community composition across 3 years. We found that male trees consistently supported higher diversity of insects than female trees and only three insect species, all Andrena spp., consistently visited both sexes. Additionally, Andrena nigrae, which was the only insect that occurred more on female than male flowers, correlated strongly to volatile cues. This suggests that cross-pollinators cue into specific aspects of floral scent, but diversity of floral visitors is driven strongly by visual cues of yellow male pollen. Through time, the floral activity of two Andrena species remained stable, but A. nigrae visited less in 2017 when flowers bloomed earlier than other years. When native bee emergence does not synchronize with bloom, activity appears to be diminished which could threaten species that subsist on a single host. Despite the community diversity of S. nigra flowers, its productivity depends on a small fraction of species that are not threatened by competition, but rather rapidly changing conditions that lead to host-insect asynchrony.

Evolutionary history and ecology shape the diversity and abundance of phytochemical arsenals across monkeyflowers.

We examine the extent to which phylogenetic effects and ecology are associated with macroevolutionary patterns of phytochemical defence production across the Mimulus phylogeny. We grew plants from 21 species representing the five major sections of the Mimulus phylogeny in a common garden to assess how the arsenals (NMDS groupings) and abundances (concentrations) of a phytochemical defence, phenylpropanoid glycosides (PPGs), vary across the phylogeny. Very few PPGs are widespread across the genus, but many are common to multiple sections of the genus. Phytochemical arsenals cluster among sections in an NMDS and are not associated with total concentration of PPGs. There is a strong phylogenetic signal for phytochemical arsenal composition across the Mimulus genus, whereas ecological variables such as growing season length, latitude, and elevation do not significantly influence arsenal. In contrast, there is little phylogenetic signal for total PPG concentration, and this trait is significantly influenced by several ecological factors. Phytochemical arsenals and abundances are influenced by plant life history form. Both phylogenetic effects and ecology are related to phytochemical patterns across species, albeit in different ways. The independence of phytochemical defence concentrations from arsenal compositions indicates that these aspects of defence may continue to evolve independently of one another.

Sex determination through X-Y heterogamety in Salix nigra.

The development of non-recombining sex chromosomes has radical effects on the evolution of discrete sexes and sexual dimorphism. Although dioecy is rare in plants, sex chromosomes have evolved repeatedly throughout the diversification of angiosperms, and many of these sex chromosomes are relatively young compared to those found in vertebrates. In this study, we designed and used a sequence capture array to identify a novel sex-linked region (SLR) in Salix nigra, a basal species in the willow clade, and demonstrated that this species has XY heterogamety. We did not detect any genetic overlap with the previously characterized ZW SLRs in willows, which map to a different chromosome. The S. nigra SLR is characterized by strong recombination suppression across a 2 MB region and an excess of low-frequency alleles, resulting in a low Tajima's D compared to the remainder of the genome. We speculate that either a recent bottleneck in population size or factors related to positive or background selection generated this differential pattern of Tajima's D on the X and autosomes. This discovery provides insights into factors that may influence the evolution of sex chromosomes in plants and contributes to a large number of recent observations that underscore their dynamic nature.

Trait plasticity and trade-offs shape intra-specific variation in competitive response in a foundation tree species.

The ability to tolerate neighboring plants (i.e. degree of competitive response) is a key determinant of plant success in high-competition environments. Plant genotypes adjust their functional trait expression under high levels of competition, which may help explain intra-specific variation in competitive response. However, the relationships between traits and competitive response are not well understood, especially in trees. In this study, we investigated among-genotype associations between tree trait plasticity and competitive response. We manipulated competition intensity in experimental stands of trembling aspen (Populus tremuloides) to address the covariance between competition-induced changes in functional trait expression and aspects of competitive ability at the genotype level. Genotypic variation in the direction and magnitude of functional trait responses, especially those of crown foliar mass, phytochemistry, and leaf physiology, was associated with genotypic variation in competitive response. Traits exhibited distinct plastic responses to competition, with varying degrees of genotypic variation and covariance with other trait responses. The combination of genotypic diversity and covariance among functional traits led to tree responses to competition that were coordinated among traits yet variable among genotypes. Such relationships between tree traits and competitive success have the potential to shape stand-level trait distributions over space and time.

A General Model to Explain Repeated Turnovers of Sex Determination in the Salicaceae.

Dioecy, the presence of separate sexes on distinct individuals, has evolved repeatedly in multiple plant lineages. However, the specific mechanisms by which sex systems evolve and their commonalities among plant species remain poorly understood. With both XY and ZW sex systems, the family Salicaceae provides a system to uncover the evolutionary forces driving sex chromosome turnovers. In this study, we performed a genome-wide association study to characterize sex determination in two Populus species, P. euphratica and P. alba. Our results reveal an XY system of sex determination on chromosome 14 of P. euphratica, and a ZW system on chromosome 19 of P. alba. We further assembled the corresponding sex-determination regions, and found that their sex chromosome turnovers may be driven by the repeated translocations of a Helitron-like transposon. During the translocation, this factor may have captured partial or intact sequences that are orthologous to a type-A cytokinin response regulator gene. Based on results from this and other recently published studies, we hypothesize that this gene may act as a master regulator of sex determination for the entire family. We propose a general model to explain how the XY and ZW sex systems in this family can be determined by the same RR gene. Our study provides new insights into the diversification of incipient sex chromosomes in flowering plants by showing how transposition and rearrangement of a single gene can control sex in both XY and ZW systems.

Metagenomics Reveals Diet-Specific Specialization of Bacterial Communities in Fungus Gardens of Grass- and Dicot-Cutter Ants.

Leaf-cutter ants in the genus Atta are dominant herbivores in the Neotropics. While most species of Atta cut dicots to incorporate into their fungus gardens, some species specialize on grasses. Here we examine the bacterial community associated with the fungus gardens of grass- and dicot-cutter ants to examine how changes in substrate input affect the bacterial community. We sequenced the metagenomes of 12 Atta fungus gardens, across four species of ants, with a total of 5.316 Gbp of sequence data. We show significant differences in the fungus garden bacterial community composition between dicot- and grass-cutter ants, with grass-cutter ants having lower diversity. Reflecting this difference in community composition, the bacterial functional profiles between the fungus gardens are significantly different. Specifically, grass-cutter ant fungus garden metagenomes are particularly enriched for genes responsible for amino acid, siderophore, and terpenoid biosynthesis while dicot-cutter ant fungus gardens metagenomes are enriched in genes involved in membrane transport. Differences between community composition and functional capacity of the bacteria in the two types of fungus gardens reflect differences in the substrates that the ants incorporated. These results show that different substrate inputs matter for fungus garden bacteria and shed light on the potential role of bacteria in mediating the ants' transition to the use of a novel substrate.

Bacteria Contribute to Plant Secondary Compound Degradation in a Generalist Herbivore System.

Herbivores must overcome a variety of plant defenses, including coping with plant secondary compounds (PSCs). To help detoxify these defensive chemicals, several insect herbivores are known to harbor gut microbiota with the metabolic capacity to degrade PSCs. Leaf-cutter ants are generalist herbivores, obtaining sustenance from specialized fungus gardens that act as external digestive systems and which degrade the diverse collection of plants foraged by the ants. There is in vitro evidence that certain PSCs harm Leucoagaricus gongylophorus, the fungal cultivar of leaf-cutter ants, suggesting a role for the Proteobacteria-dominant bacterial community present within fungus gardens. In this study, we investigated the ability of symbiotic bacteria present within fungus gardens of leaf-cutter ants to degrade PSCs. We cultured fungus garden bacteria, sequenced the genomes of 42 isolates, and identified genes involved in PSC degradation, including genes encoding cytochrome P450 enzymes and genes in geraniol, cumate, cinnamate, and α-pinene/limonene degradation pathways. Using metatranscriptomic analysis, we showed that some of these degradation genes are expressed in situ Most of the bacterial isolates grew unhindered in the presence of PSCs and, using gas chromatography-mass spectrometry (GC-MS), we determined that isolates from the genera Bacillus, Burkholderia, Enterobacter, Klebsiella, and Pseudomonas degrade α-pinene, ß-caryophyllene, or linalool. Using a headspace sampler, we show that subcolonies of fungus gardens reduced α-pinene and linalool over a 36-h period, while L. gongylophorus strains alone reduced only linalool. Overall, our results reveal that the bacterial communities in fungus gardens play a pivotal role in alleviating the effect of PSCs on the leaf-cutter ant system.IMPORTANCE Leaf-cutter ants are dominant neotropical herbivores capable of deriving energy from a wide range of plant substrates. The success of leaf-cutter ants is largely due to their external gut, composed of key microbial symbionts, specifically, the fungal mutualist L. gongylophorus and a consistent bacterial community. Both symbionts are known to have critical roles in extracting energy from plant material, yet comparatively little is known about their roles in the detoxification of plant secondary compounds. In this study, we assessed if the bacterial communities associated with leaf-cutter ant fungus gardens can degrade harmful plant chemicals. We identify plant secondary compound detoxification in leaf-cutter ant gardens as a process that depends on the degradative potential of both the bacterial community and L. gongylophorus Our findings suggest that the fungus garden and its associated microbial community influence the generalist foraging abilities of the ants, underscoring the importance of microbial symbionts in plant substrate suitability for herbivores.

Host plant genetic control of associated fungal and insect species in a Populus hybrid cross.

Plants employ a diverse set of defense mechanisms to mediate interactions with insects and fungi. These relationships can leave lasting impacts on host plant genome structure such as rapid expansion of gene families through tandem duplication. These genomic signatures provide important clues about the complexities of plant/biotic stress interactions and evolution. We used a pseudo-backcross hybrid family to identify quantitative trait loci (QTL) controlling associations between Populus trees and several common Populus diseases and insects. Using whole-genome sequences from each parent, we identified candidate genes that may mediate these interactions. Candidates were partially validated using mass spectrometry to identify corresponding QTL for defensive compounds. We detected significant QTL for two interacting fungal pathogens and three insects. The QTL intervals contained candidate genes potentially involved in physical and chemical mechanisms of host-plant resistance and susceptibility. In particular, we identified adjoining QTLs for a phenolic glycoside and Phyllocolpa sawfly abundance. There was also significant enrichment of recent tandem duplications in the genomic intervals of the native parent, but not the exotic parent. Tandem gene duplication may be an important mechanism for rapid response to biotic stressors, enabling trees with long juvenile periods to reach maturity despite many coevolving biotic stressors.

Pathways to sex determination in plants: how many roads lead to Rome?

The presence of thousands of independent origins of dioecy in angiosperms provides a unique opportunity to address the parallel evolution of the molecular pathways underlying unisexual flowers. Recent progress towards identifying sex determination genes has identified hormone response pathways, mainly associated with cytokinin and ethylene response pathways, as having been recruited multiple times independently to control unisexuality. Moreover, transcriptomics has begun to identify commonalities among intermediate sections of signal transduction pathways. These recent advances set the stage for development of a comparative evolutionary development research program to identify the shared and unique aspects of the genetic pathways of unisexual flower development in angiosperms.

To compete or defend: linking functional trait variation with life-history tradeoffs in a foundation tree species.

Although chemical deterrents to herbivory often exact costs in terms of plant growth, the manner in which those costs arise, and their physiological relationship to other functional traits, remain unclear. In the absence of appreciable herbivory, we examined interrelationships among chemical defense levels and other foliar functional traits (e.g., light-saturated photosynthesis, specific leaf area, nitrogen concentration) as co-determinants of tree growth and, by extension, competitive ability in high-density populations comprising 16 genotypes of Populus tremuloides. Across genotypes, concentrations of chemical defenses were not significantly related to other leaf functional traits, but levels of the salicinoid phenolic glycosides (SPGs) salicin, salicortin and tremulacin were each negatively correlated with relative mass growth (RMG) of aboveground woody tissue (P ≤ 0.001). RMG, in turn, underpinned 77% of the genotypic variation in relative height growth (our index of competitive ability). RMG was also positively related to light-saturated photosynthesis (P ≤ 0.001), which, together with the three SPGs, explained 86% of genotypic RMG variation (P ≤ 0.001). Moreover, results of a carbon balance simulation indicated that costs of resource allocation to SPGs, reaching nearly a third of annual crown photosynthesis, were likely mediated by substantial metabolic turnover, particularly for salicin. The lack of discernible links between foliar defense allocation and other (measured) functional traits, and the illustrated potential of metabolic turnover to reconcile influences of SPG allocation on RMG, shed additional light on fundamental physiological mechanisms underlying evolutionary tradeoffs between chemical defense investment and competitive ability in a foundation tree species.

Phylogenomics of the genus Populus reveals extensive interspecific gene flow and balancing selection.

Phylogenetic analysis is complicated by interspecific gene flow and the presence of shared ancestral polymorphisms, particularly those maintained by balancing selection. In this study, we aimed to examine the prevalence of these factors during the diversification of Populus, a model tree genus in the Northern Hemisphere. We constructed phylogenetic trees of 29 Populus taxa using 80 individuals based on re-sequenced genomes. Our species tree analyses recovered four main clades in the genus based on consensus nuclear phylogenies, but in conflict with the plastome phylogeny. A few interspecific relationships remained unresolved within the multiple-species clade because of inconsistent gene trees. Our results indicated that gene flow has been widespread within each clade and also occurred among the four clades during their early divergence. We identified 45 candidate genes with ancient polymorphisms maintained by balancing selection. These genes were mainly associated with mating compatibility, growth and stress resistance. Both gene flow and selection-mediated ancient polymorphisms are prevalent in the genus Populus. These are potentially important contributors to adaptive variation. Our results provide a framework for the diversification of model tree genus that will facilitate future comparative studies.