Stochastic mapping of morphological characters

Allometry of sound reception structures and evidence for a mandibular middle ear in non-mammalian synapsids

The origin of sensory structures provides an excellent framework for studying how constraints and selective pressures affect the evolution of complex features. The evolution of the mammalian middle ear from the jaw hinge of non-mammalian synapsids offers a deep-time perspective on sensory evolution but is limited by a poor understanding of early synapsid hearing. This work tests the hypothesis that the size of the reflected lamina of the angular bone in non-mammalian synapsids followed a strict, negative allometric trend that may be expected for a sound receiver. Allometry is first investigated in the pterygoid bone of chameleons, which was co-opted for hearing in some species and represents a possible analog for the synapsid reflected lamina. Results indicate that chameleons with a pterygoid ear exhibit a similar allometric slope, while species without a pterygoid ear have variable slopes, suggesting an optimum allometric pattern in sound receivers. In the reflected lamina, we find reduced variation around the allometric trend in therocephalians and non-bidentalian anomodonts, and evolutionary modeling suggests constraint in these groups. These results are consistent with a mandibular middle ear in non-mammalian synapsids, adding valuable new insights to the hypothesis that selective pressures for hearing ability were present long before the evolution of the mammalian middle ear.

Parallel sensory compensation following independent subterranean colonization by groundwater salamanders (Eurycea)

Lineages that have invaded subterranean environments have repeatedly evolved remarkable adaptations to life in darkness. However, observational and experimental studies in additional natural systems are needed to further our understanding of repeated evolution and convergence. In Texas, a radiation of groundwater salamanders (genus Eurycea), with independent invasions of subterranean karstic environments, offers an opportunity to investigate phenotypic convergence, parallel evolution, and the enhancement and regression of sensory systems. Adaptations to a troglobitic life in this clade include morphological, behavioral, and physiological changes within and among species. Intraspecific and interspecific variation in morphology in response to the selective pressures of life underground allows for detailed examination of physical, behavioral, and physiological changes associated with subterranean adaptation within a comparative phylogenetic framework. We find a correlated change between two sensory systems repeated across multiple subterranean Eurycea lineages: the degeneration of the eye and the expansion of the mechanosensory lateral line. The increase in anterior neuromast organs in subterranean lineages was positively correlated with the expression of pax6 (Paired-box 6), a conserved transcription factor important for vertebrate neurogenesis. Our results show a decreasing trend of PAX6 labeling in the neuromasts of adult surface salamanders (Eurycea nana) relative to the maintained labeling in subterranean species (Eurycea rathbuni). These lateral line enhancements are correlated with reductions in the development of optic systems in subterranean salamander lineages. Altogether, our findings provide a starting point for future evolutionary developmental investigations examining the genetic underpinnings of adaptive, repeated evolution in a novel system.

Arcyria and allied genera: taxonomic backbone and character evolution

Arcyria, one of the most diverse genera of bright-spored myxomycetes, has recently been identified as polyphyletic within the paraphyletic family Arcyriaceae sensu lato. However, due to the high variability and inconsistency of morphological traits in Trichiales, no taxonomic decision has yet been proposed to rectify this situation. A thorough revision of the genus is necessary. To address this, we present results of the extensive taxonomic and geographic sampling of the genus by incorporating sequences of two marker genes from 192 specimens of Arcyria spp. from 19 countries and six continents into the existing three gene phylogenetic backbone of Trichiales. Our phylogenies decisively confirm the polyphyletic nature of Arcyria, with a considerable part of the genus forming a sister clade to Hemitrichia. Consequently, we erect here the new family Hemitrichiaceae and resurrect the genus Heterotrichia to accommodate the second major branch of former Arcyria species. Additionally, we describe the new genus Spiromyxa with intermediate phylogenetic position and unique combination of morphological traits. Furthermore, we show that most of the Arcyria morphospecies included into this study are either species complexes or para- or polyphyletic taxa. Traits, previously thought to delimit the genus Arcyria, show little consistency with the phylogeny. In particular, the expanding, net-like capillitium, attached to the cup or stalk, appears to have been present in the ancestor of three families of Trichiales, afterwards mostly persisted in two of them and evolved into something else in the third. Such traits as ornamentation of capillitium, cup and spores retain their taxonomic value primarily at the species level. However, the colour of the spore mass in fresh condition partially aligns with our phylogenetic results: most of the yellow "Arcyria" species cluster within Heterotrichia, while red and grey species sort mostly into separate subclades within the clade incorporating the remaining species of the genus Arcyria. Citation: Yatsiuk I, Leontyev D, Schnittler M, Ehlers T, Mikryukov V, Kõljalg U (2025). Arcyria and allied genera: taxonomic backbone and character evolution. Fungal Systematics and Evolution 15: 97-118. doi: 10.3114/fuse.2025.15.04.

Ecosystem relocation on Snowball Earth: Polar-alpine ancestry of the extant surface biosphere?

Geological observations informed by climate dynamics imply that the oceans were 99.9% covered by light-blocking ice shelves during two discrete, self-reversing Snowball Earth epochs spanning a combined 60 to 70 Myr of the Cryogenian Period (720 to 635 Ma). The timescale for initial ice advances across the tropical oceans is ~300 y in an ice-atmosphere-ocean general circulation model in Cryogenian paleogeography. Areas of optically thin oceanic ice are usually invoked to account for fossil marine phototrophs, including macroscopic multicellular eukaryotes, before and after each Snowball, but different taxa. Ecosystem relocation is a scenario that does not require thin marine ice. Assume that long before Cryogenian Snowballs, diverse supra- and periglacial biomes were established in polar-alpine regions. When the Snowball onsets occurred, those biomes migrated in step with their ice margins to the equatorial zone of net sublimation. There, they prospered and evolved, their habitat areas expanded, and the cruelty of winter reduced. Nutrients were supplied by dust (loess) derived from cozonal ablative lands where surface winds were strong. When each Snowball finally ended, those biomes were mostly inundated by the meltwater-dominated and rapidly warming lid of a nutrient-rich but depauperate ocean. Some taxa returned to the mountaintops while others restocked the oceans. This ecosystem relocation scenario makes testable predictions. The lineages required for post-Cryogenian biotic radiations should be present in modern polar-alpine biomes. Legacies of polar-alpine ancestry should be found in the genomes of living organisms. Examples of such tests are highlighted herein.

Most bacterial gene families are biased toward specific chromosomal positions

The arrangement of genes along bacterial chromosomes influences their expression through growth rate-dependent gene copy number changes during DNA replication. Although translation- and transcription-related genes often cluster near the origin of replication, the extent of positional biases across gene families remains unclear. We hypothesized that natural selection broadly favors specific chromosomal positions to optimize growth rate-dependent expression. Analyzing 910 bacterial species and proteomics data from Escherichia coli and Bacillus subtilis, we found that about two-thirds of bacterial gene families are positionally biased. Natural selection drives genes mainly toward the origin or terminus of replication, with the strongest selection in fast-growing species. Our findings reveal chromosomal positioning as a fundamental mechanism for coordinating gene expression with growth rate, highlighting evolutionary constraints on bacterial genome architecture.

Innovation in ant larval feeding facilitated queen-worker divergence and social complexity

Building differences between genetically equivalent units is a fundamental challenge for all multicellular organisms and superorganisms. In ants, reproductive or worker fate is typically determined during the larval stage, through feeding regimes managed by adult caretakers. However, the feeding care provided to larvae varies significantly across ants, as does phenotypic divergence between queen and worker castes. Here, we employed comparative phylogenetic methods and causal inference to investigate the relationships between larval feeding care, caste size dimorphism, and social complexity across ant diversity. We digitized the life's work of George and Jeanette Wheeler, cataloging the larval morphology of over 700 species, and we compiled data on species diets and larval feeding behaviors from the literature and our own observations. We measured queen-worker size dimorphism in 392 species and gathered data for colony size, worker polymorphism, and worker reproduction. Our analyses revealed that ancestral active-feeding larvae evolved passive morphologies when adults began feeding them individually, typically with processed material and often following a shift to nonpredatory diets. Greater queen-worker size dimorphism coevolved with larval passiveness, alongside traits indicative of increased social complexity, including larger colony sizes, worker subcastes, and a reduction in workers' reproductive potential. Likelihood comparisons of causal phylogenetic models support that extended alloparental care facilitated stronger caste dimorphism, which, in turn and along with increased colony sizes, promoted higher social complexity. Our results suggest that enhanced adult control over larval development enabled greater phenotypic specialization within colonies, with profound implications for social evolution.

When the sand blossoms: Phylogeny, trait evolution, and geography of speciation in Linanthus

PREMISE

Understanding how plants successfully diversified in novel environments is a central question in evolutionary biology. Linanthus occurs in arid areas of western North America and exhibits extensive floral trait variation, multiple color polymorphisms, differences in blooming time, and variation in life history strategies. We reconstructed the evolutionary history of this genus.

METHODS

We generated restriction-site associated (ddRAD) sequences for 180 individuals and target capture (TC) sequences for 63 individuals, with complete species sampling. Using maximum likelihood and pseudo-coalescent approaches, we inferred phylogenies of Linanthus and used them to model the evolution of phenotypic traits and investigate the genus's geographic speciation history.

RESULTS

Relationships are consistent and well supported with both ddRAD and TC data. Most species are monophyletic despite extensive local sympatry and range overlap, suggesting strong isolating barriers. The non-monophyly of the night-blooming and perennial species may be due to rapid speciation or cryptic diversity. Perenniality likely evolved from annuality, a rare shift in angiosperms. Night-blooming evolved three times independently. Flower color polymorphism is an evolutionarily labile trait that is likely ancestral. No single geographic mode of speciation characterizes this diversification, but most species overlap in range, which suggests that they evolved in parapatry.

CONCLUSIONS

Our results illustrate the complexity of phylogenetic inference for recent radiations, even with multiple sources of genomic data and extensive sampling. This analysis provides a foundation for understanding aridity adaptations, such as evolution of flower color polymorphisms, night-blooming, and perenniality, as well as speciation mechanisms.

Ecological change and conflict reduction led to a social circulatory system in ants

Behavioral innovations can be ecologically transformative for lineages that perform them and for their associated communities. Many ecologically dominant, superorganismal, and speciose ant lineages use mouth-to-mouth social regurgitation behavior - stomodeal trophallaxis - to share exogenous and endogenous materials within colonies. This behavior is less common in other species-poor, less cooperative ant lineages. How and why trophallaxis evolved and fixed in only some ant clades remains unclear, and whether this trait could be indicative of superorganismality has yet to be established. Here we show that trophallaxis evolved in two main events, in non-doryline formicoids around 130 Ma and in some ponerines around 90 Ma, lineages that today encompass 86% of all ant species. We found that trophallaxis evolved in lineages that began drinking sugary liquids and that had reduced intra-colonial conflict by constraining worker reproductive potential. Evolution of trophallaxis increased net diversification. Causal models indicate that trophallaxis required low reproductive conflict and contributed to the large colony sizes of the ants that use it. This suggests that the evolution of social regurgitation was enabled by both social conflict reduction and opportunistic inclusion of nectar and honeydew in the ant diet during the shifts in terrestrial ecosystems toward flowering plants.

Bioluminescence and environmental light drive the visual evolution of deep-sea shrimp (Oplophoroidea)

Light functions as the universal language in the deep sea (>200 m). Both bioluminescent emissions and downwelling light sources dimly illuminate the water column and can drive sensory system evolution. In pelagic environments, vertically migrating animals can experience drastic changes to their lighting environment across depth, subjecting them to unique selective pressures, possibly to distinguish between changes in ambient light and bioluminescent sources. Here we show that visual opsin diversity across a group of variable vertical migrators -bioluminescent deep-sea shrimp belonging to the Superfamily Oplophoroidea- is higher among species who migrate to shallower waters with more variable light conditions. Further, we provide evidence for adaptive visual evolution among species who have evolved an additional mode of bioluminescence (photophores), including positive selection for a putative mid-wavelength sensitive opsin that may facilitate light source discrimination. Diversification of this opsin appears to play an important role in the visual ecologies of photophore-bearing shrimp with its diversification in Oplophoroidea likely playing a critical role in the fitness and evolutionary success of this group.

Assessing the Adequacy of Morphological Models Using Posterior Predictive Simulations

Reconstructing the evolutionary history of different groups of organisms provides insight into how life originated and diversified on Earth. Phylogenetic trees are commonly used to estimate this evolutionary history. Within Bayesian phylogenetics a major step in estimating a tree is in choosing an appropriate model of character evolution. While the most common character data used is molecular sequence data, morphological data remains a vital source of information. The use of morphological characters allows for the incorporation fossil taxa, and despite advances in molecular sequencing, continues to play a significant role in neontology. Moreover, it is the main data source that allows us to unite extinct and extant taxa directly under the same generating process. We therefore require suitable models of morphological character evolution, the most common being the Mk Lewis model. While it is frequently used in both palaeobiology and neontology, it is not known whether the simple Mk substitution model, or any extensions to it, provide a sufficiently good description of the process of morphological evolution. In this study we investigate the impact of different morphological models on empirical tetrapod datasets. Specifically, we compare unpartitioned Mk models with those where characters are partitioned by the number of observed states, both with and without allowing for rate variation across sites and accounting for ascertainment bias. We show that the choice of substitution model has an impact on both topology and branch lengths, highlighting the importance of model choice. Through simulations, we validate the use of the model adequacy approach, posterior predictive simulations, for choosing an appropriate model. Additionally, we compare the performance of model adequacy with Bayesian model selection. We demonstrate how model selection approaches based on marginal likelihoods are not appropriate for choosing between models with partition schemes that vary in character state space (i.e., that vary in Q-matrix state size). Using posterior predictive simulations, we found that current variations of the Mk model are often performing adequately in capturing the evolutionary dynamics that generated our data. We do not find any preference for a particular model extension across multiple datasets, indicating that there is no "one size fits all" when it comes to morphological data and that careful consideration should be given to choosing models of discrete character evolution. By using suitable models of character evolution, we can increase our confidence in our phylogenetic estimates, which should in turn allow us to gain more accurate insights into the evolutionary history of both extinct and extant taxa.

Pandanus plastomes decoded: When climate mirrors morphology and phylogenetic relationships

PREMISE

Pandanus Parkinson (Pandanaceae) is a large genus of paleotropical tree-like monocots. Previous studies using small DNA regions questioned the monophyly of the seven Pandanus subgenera, but low phylogenetic branch support hindered further investigations. We aimed to (1) test Pandanus subgeneric monophyly, (2) identify clade morphological synapomorphies, (3) investigate correlations between leaf anatomy of water storage tissue and climatic differentiation across clades, and (4) construct hypotheses on the genus' spatiotemporal history.

METHODS

We sequenced 50 Pandanus species using genome skimming and reconstructed plastomes with MITObim. We inferred partitioned RAxML phylogenetic trees to test subgeneric monophyly using Shimodaira-Hasegawa tests. We inferred a partitioned dated BEAST phylogenetic tree used for ancestral state reconstructions of morphological traits. Phylogenetic clades were used to compare climatic (Bioclim) and soil (UNESCO Digital Soil Map) conditions using random forests. We correlated present morphology and climatic niche with past climate events.

RESULTS

Our phylogenetic analyses revealed two clades and four subclades. Only subgenus Coronata was monophyletic. Staminate synapomorphies were identified for three subclades. Hypertrophied and hyperplasic water-storage tissue was a synapomorphy for clade II, correlating with more seasonal temperature and precipitation regimes and more well-draining soil. Clades differentiated during the advent of the Southeast Asian monsoon in the early Miocene, whereas subclades differentiated during the Miocene Thermal Maximum.

CONCLUSIONS

Pandanus subgeneric classification needs to be revised. Hypertrophied hyperplasic water-storage tissue is a key trait in Pandanus evolution, possibly explaining climatic and biogeographic patterns because it is key to maintaining photosynthesis during periods of hydric stress.

Tempo and mode of winter diapause evolution in butterflies

Quantifying the tempo and mode via modern phylogenetic comparative methods can provide key insights into how selection and constraints shape trait evolution on a macroevolutionary time scale. Here, we elucidate the evolution of hibernation (winter) diapause, a complex and defining life-history trait that allows temporal escape from harsh winters in temperate regions for many insects, including our model system, butterflies. Butterflies can diapause in all major life stages, and the availability of global-scale phylogenies makes them an ideal model system for studying diapause evolution. First, using a thorough literature survey, we scored the developmental stage of hibernation diapause (egg, larva, pupa, adult) vs. absence of diapause. We find that larval diapause is most common, while pupal, egg, and adult diapause are relatively rare. Next, we determined that the loss of diapause occurred at a much higher rate and that gains primarily occurred from the non-diapause state. While ancestral state estimation at deeper nodes remained uncertain, we found consistent patterns for some families and strong evidence for extensive convergence in diapause evolution. Contrary to expectations, we find no support for increased gain of diapause during the Eocene-Oligocene glaciation (~35 million years ago). Overall, the evolution of diapause in butterflies has a complex history, has evolved convergently, and has likely predated the major glaciation event consistent with the deep history of diapause evolution in insects. This study advances our understanding of the evolution of a complex and important life-history trait and establishes a macroevolutionary foundation for future studies on the ultimate and proximate basis of diapause evolution.

Divergence time and environmental similarity predict the strength of morphological convergence in stick and leaf insects

Independent evolution of similar traits in lineages inhabiting similar environments (convergent or repeated evolution) is often taken as evidence for adaptation by natural selection, and used to illustrate the predictability of evolution. Yet convergence is rarely perfect for two reasons. First, environments may not be as similar as they appear. Second, responses to selection are contingent upon available genetic variation and independent lineages may differ in the alleles, genetic backgrounds, and even the developmental mechanisms responsible for the phenotypes in question. Both impediments to convergence are predicted to increase as the length of time separating two lineages increases, making it difficult to discern their relative importance. We quantified environmental similarity and the extent of convergence to show how habitat and divergence time each contribute to observed patterns of morphological evolution in 212 species of stick and leaf insects (order Phasmatodea). Dozens of phasmid lineages independently colonized similar habitats, repeatedly evolving in parallel directions on a 23-trait morphospace, though the magnitude and direction of these shifts varied. Lineages converging toward more similar environments ended up closer on the morphospace, as did closely related lineages, and closely related lineages followed more parallel evolutionary trajectories to arrive there than more distantly related ones. Remarkably, after accounting for habitat similarity, we show that divergence time reduced the extent of convergence at a constant rate across more than 100 My of separation, suggesting even the magnitude of contingency can be predictable, given sufficient spans of time.

Nuclear phylogenomics of Eperua (Leguminosae) highlights the role of habitat and morphological lability in dispersal and diversification across Amazonia and in the Caatinga-Cerrado ecotone

Eperua is a genus of Neotropical trees that forms a major component of tropical lowland forests in Amazonia, especially in the Guiana Shield and on white-sand forests. One species occurs in the Cerrado-Caatinga ecotone, and the genus also inhabits riverine and terra firme forests. Species in Eperua exhibit one of two drastically different floral architectures and inflorescence types, each associated with distinct pollinators. Prior phylogenetic studies of Eperua have revealed an unstructured topology concerning floral architectures and inflorescence types. In addition, no investigation has been conducted on how the evolution of these traits and habitat preferences influenced the dispersal and diversification of Eperua. Using target capture sequencing, we inferred the most comprehensive phylogeny for Eperua to date, sampling all 19 known species, five for the first time. We used coalescence, concatenation, and network methods to infer the Eperua phylogeny and investigate sources of incongruence impacting resolution and support. We reconstructed the biogeographic history and ancestral states for the flower architecture, inflorescence type, and habitat preference. Our phylogenomic analyses successfully resolved relationships within Eperua, attributing conflicts between the species tree and concatenated tree to gene tree discordance linked to reticulation events. Biogeographical analyses indicate that Eperua originated and initially diversified in the white-sand forests of the Guiana Shield. A subsequent adaptation to riverine and terra firme forests enabled Eperua to expand into new habitats and regions. Still, its historical preference for white-sand forests probably accounts for its absence in the southern and western parts of Amazonia. Ancestral geographic areas and corolla morphotype reconstructions suggest that speciation in Eperua has occurred in sympatry, likely driven by pollinator shifts mediated by drastic changes in floral architecture.

Allometric Constraint Predominates Over the Acoustic Adaptation Hypothesis in a Radiation of Neotropical Treefrogs

Male frogs emit stereotypical advertisement calls to attract mates and deter conspecific rivals. The evolution of these calls is thought to be linked to anatomical constraints and the acoustic characteristics of their surroundings. The acoustic adaptation hypothesis (AAH) posits that species evolve calls that maximize propagation distance and reduce signal degradation in the environment where they are emitted. We applied phylogenetic comparative analyses to study the association of body size, vegetation density, type of aquatic ecosystem, and calling site on the evolution of acoustic traits in Cophomantini, a large radiation of Neotropical treefrogs (Hylidae). We obtained and analyzed body size, acoustic, and habitat data from a total of 112 species (58% of Cophomantini), using the most inclusive available phylogeny. We found a significant negative correlation between peak frequency, body size, and calling site, but contrary to the predictions of the AAH, we did not find support for associations among call traits and environmental characteristics. Although spectral allometry is explained by an anatomical constraint, it could also be maintained by female choice. We recommend that future studies strive to incorporate factors such as female mate preferences, eavesdropping by predators or parasites, and genetic drift.

Social mating systems in birds: resource-defense polygamy-but not lekking-is a macroevolutionarily unstable trait

Our understanding of the evolution of social mating systems is largely based on an atemporal ecological framework, whereas macroevolutionary and phylogenetic perspectives looking at the causes of mating systems variation are less developed. Here, we present analyses of the evolution of social mating systems in birds at an unprecedented scale, including 66% of the world's birds and using trait-dependent speciation and extinction models. We found that lekking (no social bond between the sexes) is very rarely lost, in accordance with the hypothesis that a male shifting to investing in one rather than multiple mates would suffer a severe fitness cost. In contrast, resource-defense polygamous lineages (with a weak, transient socio-sexual bond) frequently revert back to monogamy (strong, durable socio-sexual bond) and have an elevated extinction fraction. We tentatively attribute this to the impossibility of females settling on an optimal parental care strategy under this system. Finally, we found that most gains of lekking have been directly from monogamy rather than through an intermediate stage of resource-defense polygamy.

Evolution of intraspecific floral variation in a generalist-specialist pollination system

Intraspecific processes impact macroevolutionary patterns through individual variation, selection, and ecological specialization. According to the niche variation hypothesis, the broader ecological niche of generalist species results in an increased morphological variation among individuals either because they are constituted of diversified specialized individuals each exploiting a fraction of the species' niche or because they are constituted of true generalist individuals that experience relaxed selection. To test this hypothesis, we surveyed the individual floral morphology of species of Antillean Gesneriaceae, a group that has transitioned between specialization for hummingbird pollination and generalization multiple times throughout its evolutionary history. We characterized the profiles of corollas using geometric morphometrics and compared the intraspecific shape variance of specialists and generalists in a phylogenetic context. We used three approaches that differently accounted for the high dimensionality of morphological traits, the ancestral reconstruction of pollination syndromes over time, and the error associated with the estimation of the intraspecific variance. Our findings provide partial support for the niche variation hypothesis. If considering the whole shape in the analysis corroborated this idea, decomposing the shape into principal components indicated that not all aspects of the corolla exhibit the same pattern of variation. Specifically, pollination generalists tend to display greater intraspecific variation than specialists in terms of tubularity, but not of curvature. Accounting for the error in the variance estimation also reduced the support for the hypothesis, suggesting that larger sample sizes may be required to reach stronger conclusions. This study emphasizes the reciprocal influence between plants and their pollinators on floral morphology at different biodiversity scales and suggests that ecological strategies of species can affect patterns of morphological variation at macroevolutionary scales.

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 toward 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.

Independent Transitions to Freshwater Environments Promote Phenotypic Divergence, Not Convergence, in Stingrays

Instances of convergent or parallel evolution provide a potent model system for exploring contingency and determinism in evolutionary biology. Likewise, the multiple, independent habitat transitions from saltwater to freshwater biomes offer opportunities for studying convergent evolution within and among different vertebrate lineages. For example, stingrays have invaded freshwater habitats multiple times across different continents, sometimes even several times within the same clade (e.g., Dasyatidae). We evaluated the frequency of saltwater-freshwater invasions in stingrays, compared ecological and phenotypic diversification among freshwater and saltwater lineages, and assessed the degree of convergence among freshwater species. Despite not being morphologically distinct from saltwater stingrays, freshwater stingrays do expand the margins of stingray morphological diversity. According to our data, trophic specialists occupied non-overlapping regions of morphospace, with piscivores and molluscivores being distinct from other diet guilds. Freshwater stingrays as a group did not strongly converge morphologically, neither did freshwater rays from different lineages, which shared similar niches. These findings could be explained by there not being enough time for convergence to occur among more ancient and more recent freshwater lineages. Alternatively, the different ancestral bauplans of various freshwater ray lineages and weak selection on optimal phenotypes could promote contingency in the form of evolution along paths of least resistance.

Phenotypic Convergence Is Stronger and More Frequent in Herbivorous Fishes

Constraints on phenotypic evolution can lead to patterns of convergent evolution, by limiting the "pool" of potential phenotypes in the face of endogenous (functional, developmental) or exogenous (competition, predation) selective pressures. Evaluation of convergence depends on integrating ecological and morphological data within a robust, comparative phylogenetic context. The staggering diversity of teleost fishes offers a multitude of lineages adapted for similar ecological roles and, therefore, offers numerous replicated evolutionary experiments for exploring phenotypic convergence. However, our understanding of fish feeding systems has been primarily shaped by marine species, with the monolithic exception of freshwater cichlids. Here we use piranhas and pacus (Serrasalmidae) to explore the evolution of different feeding ecologies and their morphological proxies in Neotropical freshwater environments. Specifically, we explore whether convergence is more widespread among plant-eating fishes, arising from strong constraints on phenotypic evolution in herbivores. Using osteological micro-computed tomographic imaging (μCT), we describe the major axes of morphological variation in pacus and piranhas, regarding their diet and feeding behaviors. Next, we evaluated whether herbivorous niches are less labile than other dietary guilds and whether herbivorous species' phenotypes evolve at a slower evolutionary rate than other taxa. We then assess how convergent herbivorous taxa are, using three different suites of morphological characters (dental, jaw, and abdominal morphometrics). Ecologically, herbivory is not a dead end, exhibiting similar observed transition rates as those between carnivores and omnivores. However, we documented widespread convergence in herbivores and that herbivores have slower rates of phenotypic evolution than carnivores. Most instances of convergence are found in herbivorous taxa, specifically in frugivores and folivores. Moreover, instances of "complete" convergence, indicated by positive convergence metrics observed in more than one morphometric dataset, were only found in herbivores. Herbivores do appear to evolve under constrained circumstances, but this has not limited their ecological ability.

Pappus phenotypes and flight performance across evolutionary history in the daisy family

BACKGROUND AND AIMS

Diversity in pappus shapes and size in Asteraceae suggests an adaptive response to dispersion challenges adjusting diaspores to optimal phenotypic configurations. Here, by analysing the relationship among pappus-cypsela size relationships, flight performance and pappus types in an evolutionary context, we evaluate the role of natural selection acting on the evolution of diaspore configuration at a macro-ecological scale in the daisy family.

METHODS

To link pappus-cypsela size relationships with flight performance we collected published data on these traits from 82 species. This allowed us to translate morphometric traits in flight performance for 150 species represented in a fully resolved backbone phylogeny of the daisy family. Through ancestral reconstructions and evolutionary model selection, we assessed whether flight performance was associated with and constrained by different pappus types. Additionally, we evaluated, through phylogenetic regressions, whether species with different pappus types exhibited evolutionary allometric pappus-cypsela size relationships.

RESULTS

The setose pappus type had the highest flight performances and represented the most probable ancestral state in the family. Stepwise changes in pappus types independently led from setose to multiple instances of pappus loss with associated reduction in flight performance. Flight performance evolution was best modelled as constrained by five adaptive regimes represented by specific pappus types which correspond to specific optimal diaspore configurations that are distinct in pappus-cypsela allometric relationships.

CONCLUSIONS

Evolutionary modelling suggests natural selection as the main factor of diaspore configuration changes which proceeded towards five optima, often overcoming constraints imposed by allometric relationships and favouring evolution in certain directions. With the perspective that natural selection is the main process driving the observed patterns, various biotic and abiotic are suggested as principal drivers of transitions in diaspore configurations along space and time in the daisy family history. The results also allow discussion of evolutionary changes in a historical context.

Complexity and weak integration promote the diversity of reef fish oral jaws

Major trade-offs often manifest as axes of diversity in organismal functional systems. Overarching trade-offs may result in high trait integration and restrict the trajectory of diversification to be along a single axis. Here, we explore the diversification of the feeding mechanism in coral reef fishes to establish the role of trade-offs and complexity in a spectacular ecological radiation. We show that the primary axis of variation in the measured musculo-skeletal traits is aligned with a trade-off between mobility and force transmission, spanning species that capture prey with suction and those that bite attached prey. We found weak or no covariation between about half the traits, reflecting deviations from the trade-off axis. The dramatic trophic range found among reef fishes occurs along the primary trade-off axis, with numerous departures that use a mosaic of trait combinations to adapt the feeding mechanism to diverse challenges. We suggest that morphological evolution both along and independent of a major axis of variation is a widespread mechanism of diversification in complex systems where a global trade-off shapes major patterns of diversity. Significant additional diversity emerges as systems use weak integration and complexity to assemble functional units with many trait combinations that meet varying ecological demands.

Inferring the locomotor ecology of two of the oldest fossil squirrels: influence of operationalization, trait, body size and machine learning method

Correlations between morphology and lifestyle of extant taxa are useful for predicting lifestyles of extinct relatives. Here, we infer the locomotor behaviour of Palaeosciurus goti from the middle Oligocene and Palaeosciurus feignouxi from the lower Miocene of France using their femoral morphology and different machine learning methods. We used two ways to operationalize morphology, in the form of a geometric morphometric shape dataset and a multivariate dataset of 11 femoral traits. The predictive models were built and tested using more than half (180) of the extant species of squirrel relatives. Both traditional models such as linear discriminant analysis and more sophisticated models like neural networks had the greatest predictive power. However, the predictive power also depended on the operationalization and the femoral traits used to build the model. We also found that predictive power tended to improve with increasing body size. Contrary to previous suggestions, the older species, P. goti, was most likely arboreal, whereas P. feignouxi was more likely terrestrial. This provides further evidence that arboreality was already the most common locomotor ecology among the earliest squirrels, while a predominantly terrestrial locomotor behaviour evolved shortly afterwards, before the vast establishment of grasslands in Europe.

Evolutionary innovation accelerates morphological diversification in pufferfishes and their relatives

Evolutionary innovations have played an important role in shaping the diversity of life on Earth. However, how these innovations arise and their downstream effects on patterns of morphological diversification remain poorly understood. Here, we examine the impact of evolutionary innovation on trait diversification in tetraodontiform fishes (pufferfishes, boxfishes, ocean sunfishes, and allies). This order provides an ideal model system for studying morphological diversification owing to their range of habitats and divergent morphologies, including the fusion of the teeth into a beak in several families. Using three-dimensional geometric morphometric data for 176 extant and fossil species, we examine the effect of skull integration and novel habitat association on the evolution of innovation. Strong integration may be a requirement for rapid trait evolution and facilitating the evolution of innovative structures, like the tetraodontiform beak. Our results show that the beak arose in the presence of highly conserved patterns of integration across the skull, suggesting that integration did not limit the range of available phenotypes to tetraodontiforms. Furthermore, we find that beaks have allowed tetraodontiforms to diversify into novel ecological niches, irrespective of habitat. Our results suggest that general rules pertaining to evolutionary innovation may be more nuanced than previously thought.

Stochastic Character Mapping, Bayesian Model Selection, and Biosynthetic Pathways Shed New Light on the Evolution of Habitat Preference in Cyanobacteria

Cyanobacteria are the only prokaryotes to have evolved oxygenic photosynthesis paving the way for complex life. Studying the evolution and ecological niche of cyanobacteria and their ancestors is crucial for understanding the intricate dynamics of biosphere evolution. These organisms frequently deal with environmental stressors such as salinity and drought, and they employ compatible solutes as a mechanism to cope with these challenges. Compatible solutes are small molecules that help maintain cellular osmotic balance in high-salinity environments, such as marine waters. Their production plays a crucial role in salt tolerance, which, in turn, influences habitat preference. Among the 5 known compatible solutes produced by cyanobacteria (sucrose, trehalose, glucosylglycerol, glucosylglycerate, and glycine betaine), their synthesis varies between individual strains. In this study, we work in a Bayesian stochastic mapping framework, integrating multiple sources of information about compatible solute biosynthesis in order to predict the ancestral habitat preference of Cyanobacteria. Through extensive model selection analyses and statistical tests for correlation, we identify glucosylglycerol and glucosylglycerate as the most significantly correlated with habitat preference, while trehalose exhibits the weakest correlation. Additionally, glucosylglycerol, glucosylglycerate, and glycine betaine show high loss/gain rate ratios, indicating their potential role in adaptability, while sucrose and trehalose are less likely to be lost due to their additional cellular functions. Contrary to previous findings, our analyses predict that the last common ancestor of Cyanobacteria (living at around 3180 Ma) had a 97% probability of a high salinity habitat preference and was likely able to synthesize glucosylglycerol and glucosylglycerate. Nevertheless, cyanobacteria likely colonized low-salinity environments shortly after their origin, with an 89% probability of the first cyanobacterium with low-salinity habitat preference arising prior to the Great Oxygenation Event (2460 Ma). Stochastic mapping analyses provide evidence of cyanobacteria inhabiting early marine habitats, aiding in the interpretation of the geological record. Our age estimate of ~2590 Ma for the divergence of 2 major cyanobacterial clades (Macro- and Microcyanobacteria) suggests that these were likely significant contributors to primary productivity in marine habitats in the lead-up to the Great Oxygenation Event, and thus played a pivotal role in triggering the sudden increase in atmospheric oxygen.

Females drive postmating reproductive trait evolution across Drosophila species, but not via remating rate

While traits that contribute to premating sexual interactions are known to be wildly diverse, much less is known about the diversity of postmating (especially female) reproductive traits and the mechanisms shaping this diversity. To assess the rate, pattern, and potential drivers of postmating reproductive trait evolution, we analyzed male and female traits across up to 30 Drosophila species within a phylogenetic comparative framework. In addition to postmating reproductive morphology (e.g., sperm length, reproductive tract length and mass), we also quantified mating behaviors including female remating rate-a common proxy for the strength of postmating sexual selection. We found evidence for strong coevolution between male and female postmating traits (specifically sperm length and sperm storage organ size). However, remating rate was not associated with the rate of evolution or exaggeration of either male or female postmating reproductive morphology, once phylogenetic relatedness was accounted for. We infer that female-mediated and intersexual selection predominantly drive the evolution of our postmating morphological traits, including via divergent male and female interests in controlling paternity. In comparison, remating rate has a complex and likely secondary role in shaping this evolution, in part because this trait can be both a driver and a product of postmating selection.

Independent repeated mutations within the alphaviruses Ross River virus and Barmah Forest virus indicates convergent evolution and past positive selection in ancestral populations despite ongoing purifying selection

Ross River virus (RRV) and Barmah Forest virus (BFV) are arthritogenic arthropod-borne viruses (arboviruses) that exhibit generalist host associations and share distributions in Australia and Papua New Guinea (PNG). Using stochastic mapping and discrete-trait phylogenetic analyses, we profiled the independent evolution of RRV and BFV signature mutations. Analysis of 186 RRV and 88 BFV genomes demonstrated their viral evolution trajectories have involved repeated selection of mutations, particularly in the nonstructural protein 1 (nsP1) and envelope 3 (E3) genes suggesting convergent evolution. Convergent mutations in the nsP1 genes of RRV (residues 248 and 441) and BFV (residues 297 and 447) may be involved with catalytic enzyme mechanisms and host membrane interactions during viral RNA replication and capping. Convergent E3 mutations (RRV site 59 and BFV site 57) may be associated with enzymatic furin activity and cleavage of E3 from protein precursors assisting viral maturation and infectivity. Given their requirement to replicate in disparate insect and vertebrate hosts, convergent evolution in RRV and BFV may represent a dynamic link between their requirement to selectively 'fine-tune' intracellular host interactions and viral replicative enzymatic processes. Despite evidence of evolutionary convergence, selection pressure analyses did not reveal any RRV or BFV amino acid sites under strong positive selection and only weak positive selection for nonstructural protein sites. These findings may indicate that their alphavirus ancestors were subject to positive selection events which predisposed ongoing pervasive convergent evolution, and this largely supports continued purifying selection in RRV and BFV populations during their replication in mosquito and vertebrate hosts.

Language follows a distinct mode of extra-genomic evolution

As one of the most specific, yet most diverse of human behaviors, language is shaped by both genomic and extra-genomic evolution. Sharing methods and models between these modes of evolution has significantly advanced our understanding of language and inspired generalized theories of its evolution. Progress is hampered, however, by the fact that the extra-genomic evolution of languages, i.e. linguistic evolution, maps only partially to other forms of evolution. Contrasting it with the biological evolution of eukaryotes and the cultural evolution of technology as the best understood models, we show that linguistic evolution is special by yielding a stationary dynamic rather than stable solutions, and that this dynamic allows the use of language change for social differentiation while maintaining its global adaptiveness. Linguistic evolution furthermore differs from technological evolution by requiring vertical transmission, allowing the reconstruction of phylogenies; and it differs from eukaryotic biological evolution by foregoing a genotype vs phenotype distinction, allowing deliberate and biased change. Recognising these differences will improve our empirical tools and open new avenues for analyzing how linguistic, cultural, and biological evolution interacted with each other when language emerged in the hominin lineage. Importantly, our framework will help to cope with unprecedented scientific and ethical challenges that presently arise from how rapid cultural evolution impacts language, most urgently from interventional clinical tools for language disorders, potential epigenetic effects of technology on language, artificial intelligence and linguistic communicators, and global losses of linguistic diversity and identity. Beyond language, the distinctions made here allow identifying variation in other forms of biological and cultural evolution, developing new perspectives for empirical research.

Practical guidelines for Bayesian phylogenetic inference using Markov chain Monte Carlo (MCMC)

Phylogenetic estimation is, and has always been, a complex endeavor. Estimating a phylogenetic tree involves evaluating many possible solutions and possible evolutionary histories that could explain a set of observed data, typically by using a model of evolution. Values for all model parameters need to be evaluated as well. Modern statistical methods involve not just the estimation of a tree, but also solutions to more complex models involving fossil record information and other data sources. Markov chain Monte Carlo (MCMC) is a leading method for approximating the posterior distribution of parameters in a mathematical model. It is deployed in all Bayesian phylogenetic tree estimation software. While many researchers use MCMC in phylogenetic analyses, interpreting results and diagnosing problems with MCMC remain vexing issues to many biologists. In this manuscript, we will offer an overview of how MCMC is used in Bayesian phylogenetic inference, with a particular emphasis on complex hierarchical models, such as the fossilized birth-death (FBD) model. We will discuss strategies to diagnose common MCMC problems and troubleshoot difficult analyses, in particular convergence issues. We will show how the study design, the choice of models and priors, but also technical features of the inference tools themselves can all be adjusted to obtain the best results. Finally, we will also discuss the unique challenges created by the incorporation of fossil information in phylogenetic inference, and present tips to address them.

From molecular variations to behavioral adaptations: Unveiling adaptive epistasis in primate oxytocin system

OBJECTIVE

Our primary objective was to investigate the variability of oxytocin (OT) and the GAMEN binding motif within the LNPEP oxytocinase in primates.

MATERIALS AND METHODS

We sequenced the LNPEP segment encompassing the GAMEN motif in 34 Platyrrhini species, with 21 of them also sequenced for the OT gene. Our dataset was supplemented with primate sequences of LNPEP, OT, and the oxytocin receptor (OTR) sourced from public databases. Evolutionary analysis and coevolution predictions were made followed by the macroevolution analysis of relevant amino acids associated with phenotypic traits, such as mating systems, parental care, and litter size. To account for phylogenetic structure, we utilized two distinct statistical tests. Additionally, we calculated binding energies focusing on the interaction between Callithtrix jacchus VAMEN and Pro8OT.

RESULTS

We identified two novel motifs (AAMEN and VAMEN), challenging the current knowledge of motif conservation in placental mammals. Coevolution analysis demonstrated a correlation between GAMEN, AAMEN, and VAMEN and their corresponding OTs and OTRs. Callithrix jacchus exhibited a higher binding energy between VAMEN and Pro8OT than orthologous molecules found in humans (GAMEN and Leu8OT).

DISCUSSION

The coevolution of AAMEN and VAMEN with their corresponding OTs and OTRs suggests a functional relationship that could have contributed to specific reproductive and adaptive behaviors, including paternal care, social monogamy, and twin births, prominent traits in Cebidae species, such as marmosets and tamarins. Our findings underscore the coevolution of taxon-specific amino acids among the three studied molecules, shedding light on the oxytocinergic system as an adaptive epistatic repertoire in primates.

Phylogenetic analysis of viviparity, matrotrophy, and other reproductive patterns in chondrichthyan fishes

The reproductive diversity of extant cartilaginous fishes (class Chondrichthyes) is extraordinarily broad, reflecting more than 400 million years of evolutionary history. Among their many notable reproductive specialisations are viviparity (live-bearing reproduction) and matrotrophy (maternal provision of nutrients during gestation). However, attempts to understand the evolution of these traits have yielded highly discrepant conclusions. Here, we compile and analyse the current knowledge on the evolution of reproductive diversity in Chondrichthyes with particular foci on the frequency, phylogenetic distribution, and directionality of evolutionary changes in their modes of reproduction. To characterise the evolutionary transformations, we amassed the largest empirical data set of reproductive parameters to date covering nearly 800 extant species and analysed it via a comprehensive molecular-based phylogeny. Our phylogenetic reconstructions indicated that the ancestral pattern for Chondrichthyes is 'short single oviparity' (as found in extant holocephalans) in which females lay successive clutches (broods) of one or two eggs. Viviparity has originated at least 12 times, with 10 origins among sharks, one in batoids, and (based on published evidence) another potential origin in a fossil holocephalan. Substantial matrotrophy has evolved at least six times, including one origin of placentotrophy, three separate origins of oophagy (egg ingestion), and two origins of histotrophy (uptake of uterine secretions). In two clades, placentation was replaced by histotrophy. Unlike past reconstructions, our analysis reveals no evidence that viviparity has ever reverted to oviparity in this group. Both viviparity and matrotrophy have arisen by a variety of evolutionary sequences. In addition, the ancestral pattern of oviparity has given rise to three distinct egg-laying patterns that increased clutch (brood) size and/or involved deposition of eggs at advanced stages of development. Geologically, the ancestral oviparous pattern arose in the Paleozoic. Most origins of viviparity and matrotrophy date to the Mesozoic, while a few that are represented at low taxonomic levels are of Cenozoic origin. Coupled with other recent work, this review points the way towards an emerging consensus on reproductive evolution in chondrichthyans while offering a basis for future functional and evolutionary analyses. This review also contributes to conservation efforts by highlighting taxa whose reproductive specialisations reflect distinctive evolutionary trajectories and that deserve special protection and further investigation.

Roosting ecology drives the evolution of diverse bat landing maneuvers

Bats and birds are the only living vertebrates capable of powered flight. However, bats differ from birds in that their flight required the evolution of ascending landing maneuvers that achieve their iconic head-under-heels roosting posture. We examined the evolution of landing flight in bats and tested its association with the physical properties of roosts. Bats performed four maneuvers, each correlated with patterns of peak impact force, impulse, and roosting ecology, a critical aspect of bat biology. Our findings indicate that the common ancestor of bats performed simple, four-limbed landings, similar to extant gliding mammals, and that rotationally complex landings enhancing control over impact forces coevolved multiple times with shifts to stiff, horizontal roosts. These results suggest landing biomechanics is central to bat biology: it was critical to flight adaptation in the past, mediates roost use in the present, and may affect bats' ability to respond to deforestation in the future.

A highly resolved nuclear phylogeny uncovers strong phylogenetic conservatism and correlated evolution of fruit color and size in Solanum L

Mutualisms between plants and fruit-eating animals were key to the radiation of angiosperms. Still, phylogenetic uncertainties limit our understanding of fleshy-fruit evolution, as in the case of Solanum, a genus with remarkable fleshy-fruit diversity, but with unresolved phylogenetic relationships. We used 1786 nuclear genes from 247 species, including 122 newly generated transcriptomes/genomes, to reconstruct the Solanum phylogeny and examine the tempo and mode of the evolution of fruit color and size. Our analysis resolved the backbone phylogeny of Solanum, providing high support for its clades. Our results pushed back the origin of Solanum to 53.1 million years ago (Ma), with most major clades diverging between 35 and 27 Ma. Evolution of Solanum fruit color and size revealed high levels of trait conservatism, where medium-sized berries that remain green when ripe are the likely ancestral form. Our analyses revealed that fruit size and color are evolutionary correlated, where dull-colored fruits are two times larger than black/purple and red fruits. We conclude that the strong phylogenetic conservatism shown in the color and size of Solanum fruits could limit the influences of fruit-eating animals on fleshy-fruit evolution. Our findings highlight the importance of phylogenetic constraints on the diversification of fleshy-fruit functional traits.

Four-bar Geometry is Shared among Ecologically DivergentFish Species

Understanding the factors that influence morphological evolution is a major goal in biology. One such factor is the ability to acquire and process prey. Prey hardness and evasiveness are important properties that can impact evolution of the jaws. Similar diets and biomechanical systems have repeatedly evolved among fish lineages, providing an opportunity to test for shared patterns of evolution across distantly related organisms. Four-bar linkages are structures often used by animals to transmit force and motion during feeding and that provide an excellent system to understand the impact of diet on morphological and biomechanical evolution. Here, we tested how diet influences the evolutionary dynamics of the oral four-bar linkage system in wrasses (Family: Labridae) and cichlids (Family: Cichlidae). We found that shifts in prey hardness/evasiveness are associated with limited modifications in four-bar geometry across these two distantly related fish lineages. Wrasse and cichlid four-bar systems largely exhibit many-to-one mapping in response to dietary shifts. Across two iconic adaptive radiations of fish, an optimal four-bar geometry has largely been co-opted for different dietary functions during their extensive ecological diversification. Given the exceptional jaw diversity of both lineages, many-to-one mapping of morphology to mechanical properties may be a core feature of fish adaptive radiation.

The respiratory system influences flight mechanics in soaring birds

The subpectoral diverticulum (SPD) is an extension of the respiratory system in birds that is located between the primary muscles responsible for flapping the wing1,2. Here we survey the pulmonary apparatus in 68 avian species, and show that the SPD was present in virtually all of the soaring taxa investigated but absent in non-soarers. We find that this structure evolved independently with soaring flight at least seven times, which indicates that the diverticulum might have a functional and adaptive relationship with this flight style. Using the soaring hawks Buteo jamaicensis and Buteo swainsoni as models, we show that the SPD is not integral for ventilation, that an inflated SPD can increase the moment arm of cranial parts of the pectoralis, and that pectoralis muscle fascicles are significantly shorter in soaring hawks than in non-soaring birds. This coupling of an SPD-mediated increase in pectoralis leverage with force-specialized muscle architecture produces a pneumatic system that is adapted for the isometric contractile conditions expected in soaring flight. The discovery of a mechanical role for the respiratory system in avian locomotion underscores the functional complexity and heterogeneity of this organ system, and suggests that pulmonary diverticula are likely to have other undiscovered secondary functions. These data provide a mechanistic explanation for the repeated appearance of the SPD in soaring lineages and show that the respiratory system can be co-opted to provide biomechanical solutions to the challenges of flight and thereby influence the evolution of avian volancy.

Estimating ancestral ranges and biogeographical processes in early hominins

Historical biogeography provides crucial insights into understanding the evolutionary history of hominins. We applied maximum-likelihood and biogeographical stochastic mapping to infer the ancestral ranges of hominins and estimate the frequency of biogeographical events. These events were inferred using two time-calibrated phylogenetic trees that differ in the position of Australopithecus sediba. Results suggest that regardless of which phylogeny was selected, Northcentral Africa was the preferred ancestral region for the ancestor of the Homo-Pan clade, as well as the ancestor of Sahelanthropus and later hominins. The northern and middle part of eastern Africa was the preferred ancestral region for several clades originating at subsequent deep nodes of the trees (∼5-4 Ma). The choice of tree topology had one important effect on results: whether hominin ancestors appearing after ∼4 Ma were widespread or endemic. These different patterns highlight the biogeographic significance of the phylogenetic relationships of A. sediba. Overall, the results showed that dispersal, local extinction, and sympatry played vital roles in creating the hominin distribution, whereas vicariance and jump dispersal were not as common. The results suggested symmetry in the directionality of dispersals. Distance probably influenced how rapidly taxa colonized a new region, and dispersals often followed the closest path. These findings are potentially impacted by the imperfection of the fossil record, suggesting that the results should be interpreted cautiously.

Molecular phylogeny and divergence time of Harpalyce (Leguminosae, Papilionoideae), a lineage with amphitropical diversification in seasonally dry forests and savannas

Our knowledge of the systematics of the papilionoid legume tribe Brongniartieae has greatly benefitted from recent advances in molecular phylogenetics. The tribe was initially described to include species marked by a strongly bilabiate calyx and an embryo with a straight radicle, but recent research has placed taxa from the distantly related core Sophoreae and Millettieae within it. Despite these advances, the most species-rich genera within the Brongniartieae are still not well studied, and their morphological and biogeographical evolution remains poorly understood. Comprising 35 species, Harpalyce is one of these poorly studied genera. In this study, we present a comprehensive, multi-locus molecular phylogeny of the Brongniartieae, with an increased sampling of Harpalyce, to investigate morphological and biogeographical evolution within the group. Our results confirm the monophyly of Harpalyce and indicate that peltate glandular trichomes and a strongly bilabiate calyx with a carinal lip and three fused lobes are synapomorphies for the genus, which is internally divided into three distinct ecologically and geographically divergent lineages, corresponding to the previously recognized sections. Our biogeographical reconstructions demonstrate that Brongniartieae originated in South America during the Eocene, with subsequent pulses of diversification in South America, Mesoamerica, and Australia. Harpalyce also originated in South America during the Miocene at around 20 Ma, with almost synchronous later diversification in South America and Mexico/Mesoamerica beginning 10 Ma, but mostly during the Pliocene. Migration of Harpalyce from South to North America was accompanied by a biome and ecological shift from savanna to seasonally dry forest.

Predicting photosynthetic pathway from anatomy using machine learning

Plants with Crassulacean acid metabolism (CAM) have long been associated with a specialized anatomy, including succulence and thick photosynthetic tissues. Firm, quantitative boundaries between non-CAM and CAM plants have yet to be established - if they indeed exist. Using novel computer vision software to measure anatomy, we combined new measurements with published data across flowering plants. We then used machine learning and phylogenetic comparative methods to investigate relationships between CAM and anatomy. We found significant differences in photosynthetic tissue anatomy between plants with differing CAM phenotypes. Machine learning-based classification was over 95% accurate in differentiating CAM from non-CAM anatomy, and had over 70% recall of distinct CAM phenotypes. Phylogenetic least squares regression and threshold analyses revealed that CAM evolution was significantly correlated with increased mesophyll cell size, thicker leaves, and decreased intercellular airspace. Our findings suggest that machine learning may be used to aid the discovery of new CAM species and that the evolutionary trajectory from non-CAM to strong, obligate CAM requires continual anatomical specialization.

Evolution of bioluminescence in Anthozoa with emphasis on Octocorallia

Bioluminescence is a widespread phenomenon that has evolved multiple times across the tree of life, converging among diverse fauna and habitat types. The ubiquity of bioluminescence, particularly in marine environments where it is commonly used for communication and defense, highlights the adaptive value of this trait, though the evolutionary origins and timing of emergence remain elusive for a majority of luminous organisms. Anthozoan cnidarians are a diverse group of animals with numerous bioluminescent species found throughout the world's oceans, from shallow waters to the light-limited deep sea where bioluminescence is particularly prominent. This study documents the presence of bioluminescent Anthozoa across depth and explores the diversity and evolutionary origins of bioluminescence among Octocorallia-a major anthozoan group of marine luminous organisms. Using a phylogenomic approach and ancestral state reconstruction, we provide evidence for a single origin of bioluminescence in Octocorallia and infer the age of occurrence to around the Cambrian era, approximately 540 Ma-setting a new record for the earliest timing of emergence of bioluminescence in the marine environment. Our results further suggest this trait was largely maintained in descendants of a deep-water ancestor and bioluminescent capabilities may have facilitated anthozoan diversification in the deep sea.

Systematics and character evolution of capitate hydrozoans

Capitate hydrozoans are a morphologically and ecologically diverse hydrozoan suborder, currently including about 200 species. Being grouped in two clades, Corynida and Zancleida, these hydrozoans still show a number of taxonomic uncertainties at the species, genus and family levels. Many Capitata species established symbiotic relationships with other benthic organisms, including bryozoans, other cnidarians, molluscs and poriferans, as well as with planktonic dinoflagellates for mixotrophic relationships and with bacteria for thiotrophic ectosymbioses. Our study aimed at providing an updated and comprehensive phylogeny reconstruction of the suborder, at modelling the evolution of selected morphological and ecological characters, and at testing evolutionary relationships between the symbiotic lifestyle and the other characters, by integrating taxonomic, ecological and evolutionary data. The phylogenetic hypotheses here presented shed light on the evolutionary relationships within Capitata, with most families and genera being recovered as monophyletic. The genus Zanclea and family Zancleidae, however, were divided into four divergent clades, requiring the establishment of the new genus Apatizanclea and the new combinations for species in Zanclea and Halocoryne genera. The ancestral state reconstructions revealed that symbiosis arose multiple times in the evolutionary history of the Capitata, and that homoplasy is a common phenomenon in the group. Correlations were found between the evolution of symbiosis and morphological characters, such as the perisarc. Overall, our results highlighted that the use of genetic data and a complete knowledge of the life cycles are strongly needed to disentangle taxonomic and systematic issues in capitate hydrozoans. Finally, the colonization of tropical habitat appears to have influenced the evolution of a symbiotic lifestyle, playing important roles in the evolution of the group.

On the origin of bird's nest fungi: Phylogenomic analyses of fungi in the Nidulariaceae (Agaricales, Basidiomycota)

Nidulariaceae, also known as bird's nest fungi, is an understudied group of mushroom-forming fungi. The common name is derived from their nest-like morphology. Bird's nest fungi are ubiquitous wood decomposers or saprobes on dung. Recent studies showed that species in the Nidulariaceae form a monophyletic group with five sub-clades. However, phylogenetic relationships among genera and placement of Nidulariaceae are still unclear. We present phylogenomic analyses of bird's nest fungi and related Agaricales fungi to gain insight into the evolution of Nidulariaceae. A species tree with 17 newly generated genomes of bird's nest fungi and representatives from all major clades of Agaricales was constructed using 1044 single-copy genes to explore the intergeneric relationships and pinpoint the placement of Nidulariaceae within Agaricales. We corroborated the hypothesis that bird's nest fungi are sister to Squamanitaceae, which includes mushroom-shaped fungi with a stipe and pileus that are saprobes and mycoparasites. Lastly, stochastic character mapping of discrete traits on phylogenies (SIMMAP) suggests that the ancestor of bird's nest fungi likely possessed an evanescent, globose peridium without strings attaching to the spore packets (funiculi). This analysis suggests that the funiculus was gained twice and that the persistent, cupulate peridium form was gained at least four times and lost once. However, alternative coding schemes and datasets with a wider array of Agaricales produced conflicting results during ancestral state reconstruction, indicating that there is some uncertainty in the number of peridium transitions and that taxon sampling may significantly alter ancestral state reconstructions. Overall, our results suggest that several key morphological characters of Nidulariaceae have been subject to homoplasy.

Hagfish genome elucidates vertebrate whole-genome duplication events and their evolutionary consequences

Polyploidy or whole-genome duplication (WGD) is a major event that drastically reshapes genome architecture and is often assumed to be causally associated with organismal innovations and radiations. The 2R hypothesis suggests that two WGD events (1R and 2R) occurred during early vertebrate evolution. However, the timing of the 2R event relative to the divergence of gnathostomes (jawed vertebrates) and cyclostomes (jawless hagfishes and lampreys) is unresolved and whether these WGD events underlie vertebrate phenotypic diversification remains elusive. Here we present the genome of the inshore hagfish, Eptatretus burgeri. Through comparative analysis with lamprey and gnathostome genomes, we reconstruct the early events in cyclostome genome evolution, leveraging insights into the ancestral vertebrate genome. Genome-wide synteny and phylogenetic analyses support a scenario in which 1R occurred in the vertebrate stem-lineage during the early Cambrian, and 2R occurred in the gnathostome stem-lineage, maximally in the late Cambrian-earliest Ordovician, after its divergence from cyclostomes. We find that the genome of stem-cyclostomes experienced an additional independent genome triplication. Functional genomic and morphospace analyses demonstrate that WGD events generally contribute to developmental evolution with similar changes in the regulatory genome of both vertebrate groups. However, appreciable morphological diversification occurred only in the gnathostome but not in the cyclostome lineage, calling into question the general expectation that WGDs lead to leaps of bodyplan complexity.

Exploring the Climatic Niche Evolution of the Genus Falco (Aves: Falconidae) in Europe

By integrating species distribution modeling techniques, phylogenetic comparative methods, and climatic data, we analyzed how European falcon climatic niches have changed over evolutionary time in order to understand their tempo and mode of evolution and gain phylogenetic insights related to the ecological context of falcon evolution. For this purpose, we tested the relative contributions of niche conservatism, convergent evolution, and divergent evolution in the evolutionary history of this group of species in Europe. The occupation of climatic niche spaces by falcon species in Europe was not similar, considering that their climatic niche evolution was characterized by heterotachy, especially after ca. 4 Mya. Our results indicate that convergent evolution and niche divergence played an important role in the evolutionary history of these species, with no significant evidence of closely related species retaining their fundamental niche over time (phylogenetic niche conservatism). In most analyses, less closely related falcon species occupied similar climatic environments. We found that speciation in the European genus Falco was influenced by climatic niche differentiation, more prevalent in the last 4 million years, with the main climatic niche shifts occurring between closely related falcon species.

Phylogeny, biogeography, and character evolution of Anaphalis (Gnaphalieae, Asteraceae)

The HAP clade, mainly including Helichrysum Mill, Anaphalis DC., and Pseudognaphalium Kirp., is a major component of tribe Gnaphalieae (Asteraceae). In this clade, Anaphalis represents the largest genus of Asian Gnaphalieae. The intergeneric relationships among Anaphalis and its related genera and the infrageneric taxonomy of this genus are complex and remain controversial. However, there are few studies that have focused on these issues. Herein, based on the current most comprehensive sampling of the HAP clade, especially Anaphalis, we conducted phylogenetic analyses using chloroplast (cp) genome and nuclear ribosomal DNA (nrDNA) to evaluate the relationships within HAP clade, test the monophyly of Anaphalis, and examine the infrageneric taxonomy of this genus. Meanwhile, the morphological characters were verified to determine the circumscription and infrageneric taxonomy system of Anaphalis. Additionally, the biogeographical history, diversification processes, and evolution of crucial morphological characters were estimated and inferred. Our phylogenetic analyses suggested that Anaphalis is polyphyletic because it nested with Helichrysum and Pseudognaphalium. Two and four main clades of Anaphalis were identified in cp genome and nrDNA trees, respectively. Compared with nrDNA trees, the cp genome trees were more effective for phylogenetic resolution. After comprehensively analyzing morphological and phylogenetic evidence, it was concluded that the achene surface ornamentation and leaf base showed less homoplasy and supported the two Anaphalis lineages that were inferred from cp genome. Our biogeographical analyses based on cp genome indicated that HAP clade underwent rapid diversification from late Miocene to Pliocene. The two Anaphalis lineages appeared to have originated in Africa, then spread to Western and Southern Asia, and subsequently moved into Southwestern China forming a diversity center. The dispersal patterns of the two Anaphalis lineages were different. One dispersed around the world, except in Africa and South America. The other one dispersed to Eastern and Southeastern Asia from the ancestral origin region.

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.

True grit? Comparative anatomy and evolution of gizzards in fishes

Gut morphology frequently reflects the food organisms digest. Gizzards are organs of the gut found in archosaurs and fishes that mechanically reduce food to aid digestion. Gizzards are thought to compensate for edentulism and/or provide an advantage when consuming small, tough food items (e.g., phytoplankton and algae). It is unknown how widespread gizzards are in fishes and how similar these structures are among different lineages. Here, we investigate the distribution of gizzards across bony fishes to (1) survey different fishes for gizzard presence, (2) compare the histological structure of gizzards in three species, (3) estimate how often gizzards have evolved in fishes, and (4) explore whether anatomical and ecological traits like edentulism and microphagy predict gizzard presence. According to our analyses, gizzards are rare across bony fishes, evolving only six times in a broad taxonomic sampling of 51 species, and gizzard presence is not clearly correlated with factors like gut length or dentition. We find that gizzard morphology varies among the lineages where one is present, both macroscopically (presence of a crop) and microscopically (varying tissue types). We conclude that gizzards likely aid in the mechanical reduction of food in fishes that have lost an oral dentition in their evolutionary past; however, the relative scarcity of gizzards suggests they are just one of many possible solutions for processing tough, nutrient-poor food items. Gizzards have long been present in the evolutionary history of fishes, can be found in a wide variety of marine and freshwater clades, and likely have been overlooked in many taxa.

Cryogenian Origins of Multicellularity in Archaeplastida

Earth was impacted by global glaciations during the Cryogenian (720 to 635 million years ago; Ma), events invoked to explain both the origins of multicellularity in Archaeplastida and radiation of the first land plants. However, the temporal relationship between these environmental and biological events is poorly established, due to a paucity of molecular and fossil data, precluding resolution of the phylogeny and timescale of archaeplastid evolution. We infer a time-calibrated phylogeny of early archaeplastid evolution based on a revised molecular dataset and reappraisal of the fossil record. Phylogenetic topology testing resolves deep archaeplastid relationships, identifying two clades of Viridiplantae and placing Bryopsidales as sister to the Chlorophyceae. Our molecular clock analysis infers an origin of Archaeplastida in the late-Paleoproterozoic to early-Mesoproterozoic (1712 to 1387 Ma). Ancestral state reconstruction of cytomorphological traits on this time-calibrated tree reveals many of the independent origins of multicellularity span the Cryogenian, consistent with the Cryogenian multicellularity hypothesis. Multicellular rhodophytes emerged 902 to 655 Ma while crown-Anydrophyta (Zygnematophyceae and Embryophyta) originated 796 to 671 Ma, broadly compatible with the Cryogenian plant terrestrialization hypothesis. Our analyses resolve the timetree of Archaeplastida with age estimates for ancestral multicellular archaeplastids coinciding with the Cryogenian, compatible with hypotheses that propose a role of Snowball Earth in plant evolution.

phytools 2.0: an updated R ecosystem for phylogenetic comparative methods (and other things)

Phylogenetic comparative methods comprise the general endeavor of using an estimated phylogenetic tree (or set of trees) to make secondary inferences: about trait evolution, diversification dynamics, biogeography, community ecology, and a wide range of other phenomena or processes. Over the past ten years or so, the phytools R package has grown to become an important research tool for phylogenetic comparative analysis. phytools is a diverse contributed R library now consisting of hundreds of different functions covering a variety of methods and purposes in phylogenetic biology. As of the time of writing, phytools included functionality for fitting models of trait evolution, for reconstructing ancestral states, for studying diversification on trees, and for visualizing phylogenies, comparative data, and fitted models, as well numerous other tasks related to phylogenetic biology. Here, I describe some significant features of and recent updates to phytools, while also illustrating several popular workflows of the phytools computational software.

Bare parts in the Galliformes: the evolution of a multifunctional structure

A morphological trait can have multiple functions shaped by varying selective forces. Bare parts in birds, such as wattles, casques and combs, are known to function in both signalling and thermoregulation. Studies have demonstrated such structures are targets of sexual selection via female choice in several species of Galliformes (junglefowl, turkeys and grouse), though other studies have shown some role in thermoregulation (guineafowl). Here, we tested fundamental hypotheses regarding the evolution and maintenance of bare parts in Galliformes. Using a phylogeny that included nearly 90% of species in the order, we evaluated the role of both sexual and natural selection in shaping the function of bare parts across different clades. We found a combination of both environmental and putative sexually selected traits strongly predicted the variation of bare parts for both males and females across Galliformes. When the analysis is restricted to the largest family, Phasianidae (pheasants, junglefowl and allies), sexually selected traits were the primary predictors of bare parts. Our results suggest that bare parts are important for both thermoregulation and sexual signalling across Galliformes but are primarily under strong sexual selection within the Phasianidae.

Integrative taxonomy of Serrasentis gibsoni n. sp. (Acanthocephala: Isthmosacanthidae) from flatfishes in the Gulf of Mexico

The Isthmosacanthidae acanthocephalan species of the genus Serrasentis are parasites of marine teleosts and an elasmobranch. In this study, Serrasentis gibsoni n. sp. is described from the intestines of four flatfish species (Paralichthyidae), namely Ancyclopsetta quadrocellata, Cyclopsetta chittendeni, Syacium gunteri, and S. papillosum from 10 oceanic sites in the Gulf of Mexico (GoM). Twenty sequences of the 'barcoding' region of cytochrome C oxidase subunit I gene were obtained from 20 adults of Serrasentis gibsoni n. sp. Additionally, five sequences of the barcoding region were obtained from five adults of rhadinorhynchid Gorgorhynchus lepidus from C. chittendeni, S. papillosum and one species of Haemulidae, Haemulom aurolineatum, from five oceanic sites from the GoM. Two phylogenetic approaches were followed: Bayesian inference and maximum likelihood. In both phylogenetic reconstructions, the sequences of Serrasentis gibsoni n. sp. were recovered as a monophyletic group within the genus Serrasentis and placed as a sister group to G. lepidus. However, due to the lack of molecular data for species of the Isthmosacanthidae and Rhadinorhynchidea, these phylogenetic inferences must be taken with caution. Serrasentis gibsoni n. sp. is the first species of Serrasentis described from Paralichthyidae flatfish species from marine waters of the Americas and from the GoM. Based on the barcoding data set analyzed, Serrasentis gibsoni n. sp. appears to have high intraspecific genetic variation; thus, it is necessary to continue exploring the genetic diversity of this species to infer its intraspecific evolutionary patterns.