How to measure and test phylogenetic signal

Evolutionary history and root trait coordination predict nutrient strategy in tropical legume trees

Plants express diverse nutrient use and acquisition traits, but it is unclear how trait combinations at the species level are constrained by phylogeny, trait coordination, or trade-offs in resource investment. One trait - nitrogen (N) fixation - is assumed to correlate with other traits and used to define plant functional groups, despite potential confounding effects of phylogeny. We quantified growth, carbon metabolism, fixation rate, root phosphatase activity (RPA), mycorrhizal colonization, and leaf and root morphology/chemistry across 22 species of fixing and nonfixing tropical Fabaceae trees under common conditions. Belowground trait variation was high even among closely related species, and most traits displayed a phylogenetic signal, including N-fixation rate and nodule biomass. Across species, we observed strong positive correlations between physiological traits such as RPA and root respiration. RPA increased ~ fourfold per unit increase in fixation, supporting the debated hypothesis that N-fixers 'trade' N for phosphatases to enhance phosphorus acquisition. Specific root length and root N differed between functional groups, though for other traits, apparent differences became nonsignificant after accounting for phylogenetic nonindependence. We conclude that evolutionary history, trait coordination, and fixation ability contribute to nutrient trait expression at the species level, and recommend explicitly considering phylogeny in analyses of functional groupings.

Comparative muscle anatomy of the anuran pelvis and hindlimb in relation to locomotor mode

Frogs have a highly conserved body plan, yet they employ a diverse array of locomotor modes, making them ideal organisms for investigating the relationships between morphology and locomotor function, in particular whether anatomical complexity is a prerequisite for functional complexity. We use diffusible iodine contrast-enhanced microCT (diceCT) imaging to digitally dissect the gross muscle anatomy of the pelvis and hindlimbs for 30 species of frogs representing five primary locomotor modes, including the first known detailed dissection for some of the world's smallest frogs, forming the largest digital comparative analysis of musculoskeletal structure in any vertebrate clade to date. By linking musculoskeletal dissections and phylogenetic comparative methods, we then quantify and compare relationships between anatomy and function across over 160 million years of anuran evolution. In summary, we have found that bone lengths and pelvic crest sizes are generally not reliable predictors of muscle sizes, which highlights important implications for future palaeontological studies. Our investigation also presents previously unreported differences in muscle anatomy between frogs specialising in different locomotor modes, including several of the smallest frog hindlimb muscles, which are extremely difficult to extract and measure using traditional approaches. Furthermore, we find evidence of many-to-one and one-to-many mapping of form to function across the phylogeny. Additionally, we perform the first quantitative analysis of how the degree of muscle separation can differ between frogs. We find evidence that phylogenetic history is the key contributing factor to muscle separation in the pelvis and thigh, while the separation of shank muscles is influenced more strongly by locomotor mode. Finally, our anatomical 3D reconstructions are published alongside this manuscript to contribute towards future research and serve as educational materials.

Evolutionary drivers of sex-specific parasite prevalence in wild birds

Males and females often differ in ecology, behaviour and lifestyle, and these differences are expected to lead to sex differences in parasite susceptibility. However, neither the sex differences in parasite prevalence, nor their ecological and evolutionary drivers have been investigated across a broad range of taxa using phylogenetically corrected analyses. Using the most extensive dataset yet that includes 755 prevalence estimates from 151 wild bird species in a meta-analytic framework, here we compare sex differences in blood and gastrointestinal parasites. We show that despite sex differences in parasite infection being frequently reported in the literature, only Haemoproteus infections were more prevalent in females than in males. Notably, only seasonality was strongly associated with the sex-specific parasite prevalence of both Leucocytozoon and Haemoproteus, where birds showed greater female bias in prevalence during breeding periods compared to the non-breeding period. No other ecological or sexual selection variables were associated with sex-specific prevalence of parasite prevalence. We suggest that much of the variation in sex-biased prevalence could be idiosyncratic, and driven by local ecology and behavioural differences of the parasite and the host. Therefore, breeding ecology and sexual selection may only have a modest influence on sex-different parasite prevalence across wild birds.

Phylogeny predicts sensitivity in aquatic animals for only a minority of chemicals

There are substantial gaps in our empirical knowledge of the effects of chemical exposure on aquatic life that are unlikely to be filled by traditional laboratory toxicity testing alone. One possible alternative of generating new toxicity data is cross-species extrapolation (CSE), a statistical approach in which existing data are used to predict the effect of a chemical on untested species. Some CSE models use relatedness as a predictor of chemical sensitivity, but relatively little is known about how strongly shared evolutionary history influences sensitivity across all chemicals. To address this question, we conducted a survey of phylogenetic signal in the toxicity data from aquatic animal species for a large set of chemicals using a phylogeny inferred from taxonomy. Strong phylogenetic signal was present in just nine of thirty-six toxicity datasets, and there were no clear shared properties among those datasets with strong signal. Strong signal was rare even among chemicals specifically developed to target insects, meaning that these chemicals may be equally lethal to non-target taxa, including chordates. When signal was strong, distinct patterns of sensitivity were evident in the data, which may be informative when assembling toxicity datasets for regulatory use. Although strong signal does not appear to manifest in aquatic toxicity data for most chemicals, we encourage additional phylogenetic evaluations of toxicity data in order to guide the selection of CSE tools and as a means to explore the patterns of chemical sensitivity across the broad diversity of life.

Comparative histology of abscission zones reveals the extent of convergence and divergence in seed shattering in weedy and cultivated rice

The modification of seed shattering has been a recurring theme in rice evolution. The wild ancestor of cultivated rice disperses its seeds, but reduced shattering was selected during multiple domestication events to facilitate harvesting. Conversely, selection for increased shattering occurred during the evolution of weedy rice, a weed invading cultivated rice fields that has originated multiple times from domesticated ancestors. Shattering requires formation of a tissue known as the abscission zone (AZ), but how the AZ has been modified throughout rice evolution is unclear. We quantitatively characterized the AZ characteristics of relative length, discontinuity, and intensity in 86 cultivated and weedy rice accessions. We reconstructed AZ evolutionary trajectories and determined the degree of convergence among different cultivated varieties and among independent weedy rice populations. AZ relative length emerged as the best feature to distinguish high and low shattering rice. Cultivated varieties differed in average AZ morphology, revealing lack of convergence in how shattering reduction was achieved during domestication. In contrast, weedy rice populations typically converged on complete AZs, irrespective of origin. By examining AZ population-level morphology, our study reveals its evolutionary plasticity, and suggests that the genetic potential to modify the ecologically and agronomically important trait of shattering is plentiful in rice lineages.

Insights into the evolution of photoreceptor oil droplets in frogs and toads

Photoreceptor oil droplets (ODs) are spherical organelles placed most commonly within the inner segment of the cone photoreceptors. Comprising neutral lipids, ODs can be either non-pigmented or pigmented and have been considered optically functional in various studies. Among living amphibians, ODs were only reported to occur in frogs and toads (Anura), while they are absent in salamanders and caecilians. Nonetheless, the limited understanding of their taxonomic distribution in anurans impedes a comprehensive assessment of their evolution and relationship with visual ecology. We studied the retinae of 134 anuran species, extending the knowledge of the distribution of ODs to 46 of the 58 currently recognized families, and providing a new perspective on this group that complements the available information from other vertebrates. The occurrence of ODs in anurans shows a strong phylogenetic signal, and our findings revealed that ODs evolved at least six times during the evolutionary history of the group, independently from other vertebrates. Although no evident correlation was found between OD occurrence, adult habits and diel activity, it is inferred that each independent origin involves distinct scenarios in the evolution of ODs concerning photic habits. Furthermore, our results revealed significant differences in the size of the ODs between nocturnal and arrhythmic anurans relative to the length of the cones' outer segment.

shinyTempSignal: an R shiny application for exploring temporal and other phylogenetic signals

The molecular clock model is fundamental for inferring species divergence times from molecular sequences. However, its direct application may introduce significant biases due to sequencing errors, recombination events, and inaccurately labeled sampling times. Improving accuracy necessitates rigorous quality control measures to identify and remove potentially erroneous sequences. Furthermore, while not all branches of a phylogenetic tree may exhibit a clear temporal signal, specific branches may still adhere to the assumptions, with varying evolutionary rates. Supporting a relaxed molecular clock model better aligns with the complexities of evolution. The root-to-tip regression method has been widely used to analyze the temporal signal in phylogenetic studies and can be generalized for detecting other phylogenetic signals. Despite its utility, there remains a lack of corresponding software implementations for broader applications. To address this gap, we present shinyTempSignal, an interactive web application implemented with the shiny framework, available as an R package and publicly accessible at https://github.com/YuLab-SMU/shinyTempSignal. This tool facilitates the analysis of temporal and other phylogenetic signals under both strict and relaxed models. By extending the root-to-tip regression method to diverse signals, shinyTempSignal helps in the detection of evolving features or traits, thereby laying the foundation for deeper insights and subsequent analyses.

Male and female contributions to diversity among birdwing butterfly images

Machine learning (ML) newly enables tests for higher inter-species diversity in visible phenotype (disparity) among males versus females, predictions made from Darwinian sexual selection versus Wallacean natural selection, respectively. Here, we use ML to quantify variation across a sample of > 16,000 dorsal and ventral photographs of the sexually dimorphic birdwing butterflies (Lepidoptera: Papilionidae). Validation of image embedding distances, learnt by a triplet-trained, deep convolutional neural network, shows ML can be used for automated reconstruction of phenotypic evolution achieving measures of phylogenetic congruence to genetic species trees within a range sampled among genetic trees themselves. Quantification of sexual disparity difference (male versus female embedding distance), shows sexually and phylogenetically variable inter-species disparity. Ornithoptera exemplify high embedded male image disparity, diversification of selective optima in fitted multi-peak OU models and accelerated divergence, with cases of extreme divergence in allopatry and sympatry. However, genus Troides shows inverted patterns, including comparatively static male embedded phenotype, and higher female than male disparity - though within an inferred selective regime common to these females. Birdwing shapes and colour patterns that are most phenotypically distinctive in ML similarity are generally those of males. However, either sex can contribute majoritively to observed phenotypic diversity among species.

A trait-based root acquisition-defence-decomposition framework in angiosperm tree species

To adapt to the complex belowground environment, plants make trade-offs between root resource acquisition and defence ability. This includes forming partnerships with different types of root associating microorganisms, such as arbuscular mycorrhizal and ectomycorrhizal fungi. These trade-offs, by mediating root chemistry, exert legacy effects on nutrient release during decomposition, which may, in turn, affect the ability of new roots to re-acquire resources, thereby generating a feedback loop. However, the linkages at the basis of this potential feedback loop remain largely unquantified. Here, we propose a trait-based root 'acquisition-defence-decomposition' conceptual framework and test the strength of relevant linkages across 90 angiosperm tree species. We show that, at the plant species level, the root-fungal symbiosis gradient within the root economics space, root chemical defence (condensed tannins), and root decomposition rate are closely linked, providing support to this framework. Beyond the dichotomy between arbuscular mycorrhizal-dominated versus ectomycorrhizal-dominated systems, we suggest a continuous shift in feedback loops, from 'high arbuscular mycorrhizal symbiosis-low defence-fast decomposition-inorganic nutrition' by evolutionarily ancient taxa to 'high ectomycorrhizal symbiosis-high defence-slow decomposition-organic nutrition' by more modern taxa. This 'acquisition-defence-decomposition' framework provides a foundation for testable hypotheses on multidimensional linkages between species' belowground strategies and ecosystem nutrient cycling in an evolutionary context.

Rapid Targeted Assembly of the Proteome Reveals Evolutionary Variation of GC Content in Avian Lice

Nucleotide base composition plays an influential role in the molecular mechanisms involved in gene function, phenotype, and amino acid composition. GC content (proportion of guanine and cytosine in DNA sequences) shows a high level of variation within and among species. Many studies measure GC content in a small number of genes, which may not be representative of genome-wide GC variation. One challenge when assembling extensive genomic data sets for these studies is the significant amount of resources (monetary and computational) associated with data processing, and many bioinformatic tools have not been optimized for resource efficiency. Using a high-performance computing (HPC) cluster, we manipulated resources provided to the targeted gene assembly program, automated target restricted assembly method (aTRAM), to determine an optimum way to run the program to maximize resource use. Using our optimum assembly approach, we assembled and measured GC content of all of the protein-coding genes of a diverse group of parasitic feather lice. Of the 499 426 genes assembled across 57 species, feather lice were GC-poor (mean GC = 42.96%) with a significant amount of variation within and between species (GC range = 19.57%-73.33%). We found a significant correlation between GC content and standard deviation per taxon for overall GC and GC3, which could indicate selection for G and C nucleotides in some species. Phylogenetic signal of GC content was detected in both GC and GC3. This research provides a large-scale investigation of GC content in parasitic lice laying the foundation for understanding the basis of variation in base composition across species.

Interconnected Codons: Unravelling the Epigenetic Significance of Flanking Sequences in CpG Dyads

Hypothesizing that CpG codon dyads, formed by consecutive codons containing a cytosine-guanine pair (NNC-GNN), may play a crucial role in gene function, we conducted an extensive analysis to investigate their distribution and conservation within mammalian genes. Our findings reveal that genes characterized by a high density of CpG codon dyads are notably associated with homeobox domains and RNA polymerase II transcription factors. Conversely, genes exhibiting low CpG codon dyad density have links to DNA damage repair and mitosis. Importantly, our study identifies a re-markable increase in expressed genes that harbor CpG during embryonic development, suggesting their potential involvement in gene regulation at these developmental stages. These results under-score the functional significance of CpG codon dyads in DNA methylation and gene expression, fur-ther demonstrating the coevolution of consecutive codons and their contribution to codon usage bias.

Environmental filtering mediates desert ant community assembly at two spatial scales

Understanding the mechanisms that maintain species coexistence and determine patterns of community assembly are fundamental goals of ecology. Quantifying the relationship between species traits and stress gradients is a necessary step to disentangle assembly processes and to be able to predict the outcome of environmental change. We examined the hypothesis that desert ant communities are assembled by niche-based processes i.e., environmental filtering and limiting similarity. First, we used population-level morphological trait measurements to study the functional structure of ant communities along a dryland environmental stress gradient. Second, we developed species distribution models for each species to quantify large-scale climatic niche overlap between species. Body, femur, antennal scape, and head lengths were correlated with environmental gradients. Regionally, the ant community was significantly and functionally overdispersed in terms of morphological traits which suggests the importance of competition to ant community structure. Ant community assembly was also strongly influenced by environmental factors as the degree of functional trait divergence, but not phylogenetic divergence, decreased with increasing environmental stress. Thus, environmental stress likely mediates limiting similarity in these desert ecosystems. Species with lower climatic niche overlap were more dissimilar in morphological traits. This suggests that environmental filtering on ant functional traits is important at the scale of species distributions in addition to regional scales. This study shows that environmental and biotic filtering (i.e., niche-based assembly mechanisms) are jointly and non-independently structuring the ant community.

Climatic niche evolution and niche conservatism of Nymphaea species in Africa, South America, and Australia

BACKGROUND

Interest in the evolution of climatic niches, particularly in understanding the potential adaptive responses of species under climate change, has increased both theoretically and within macroecological studies. These studies have provided valuable insights into how climatic traits of species influence their niche evolution. In this study, we aim to investigate whether niche conservatism plays a role in the species diversification of Nymphaea, a group of aquatic plants with a cosmopolitan distribution that is facing severe habitat loss. We will use climatic models and phylogenetic data for 23 species to reconstruct Nymphaea's niche evolution, measure niche overlap, and assess disparity through time while testing for evolutionary models.

RESULTS

There was a lot of overlap in niches both within and between groups, especially for species that can be found in many places. The breadth and peaks of the niche profile varied depending on the bioclimatic variables, which suggested that the species evolved differently to cope with changes in climate. The analysis also showed that evolutionary changes happened across the phylogeny, with weak to moderate signals. The morphological disparity index (MDI) values indicated that there were disparities within subclades over time but not between or among them. Niche reconstruction and evolution analysis revealed both convergent and divergent evolution among various variables. For example, N. immutabilis, N. atrans, N. violancea, and N. nouchali evolved towards intermediate temperatures for bio2 and bio3 (isothermity) while moving towards extreme temperatures for bio8 and bio9 (wettest and driest average quarterly temperatures).

CONCLUSION

Our study will improve our understanding of how changes in climatic niches are potentially driving the evolution of Nymphaea. It has significant scientific implications for the limits, assemblages, evolution, and diversification of species. This information is crucial for the ongoing efforts of conservation and management, particularly considering the inevitable effects of climate change.

Fast Bayesian Inference of Phylogenies from Multiple Continuous Characters

Time-scaled phylogenetic trees are an ultimate goal of evolutionary biology and a necessary ingredient in comparative studies. The accumulation of genomic data has resolved the tree of life to a great extent, yet timing evolutionary events remain challenging if not impossible without external information such as fossil ages and morphological characters. Methods for incorporating morphology in tree estimation have lagged behind their molecular counterparts, especially in the case of continuous characters. Despite recent advances, such tools are still direly needed as we approach the limits of what molecules can teach us. Here, we implement a suite of state-of-the-art methods for leveraging continuous morphology in phylogenetics, and by conducting extensive simulation studies we thoroughly validate and explore our methods' properties. While retaining model generality and scalability, we make it possible to estimate absolute and relative divergence times from multiple continuous characters while accounting for uncertainty. We compile and analyze one of the most data-type diverse data sets to date, comprised of contemporaneous and ancient molecular sequences, and discrete and continuous morphological characters from living and extinct Carnivora taxa. We conclude by synthesizing lessons about our method's behavior, and suggest future research venues.

Habitat differentiation and environmental adaptability contribute to leaf size variations globally in C3 and C4 grasses

The grass family (Poaceae) dominates ~43 % of Earth's land area and contributes 33 % of terrestrial primary productivity that is critical to naturally regulating atmosphere CO2 concentration and global climate change. Currently grasses comprise ~11,780 species and ~50 % of them (~6000 species) utilize C4 photosynthetic pathway. Generally, grass species have smaller leaves under colder and drier environments, but it is unclear whether the primary drivers of leaf size differ between C3 and C4 grasses on a global scale. Here, we analyzed 34 environmental variables, such as latitude, elevation, mean annual temperature, mean annual precipitation, and solar radiation etc., through a comparatively comprehensive database of ~3.0 million occurrence records from 1380 C3 and 978 C4 grass species (2358 species in total). Results from this study confirm that C4 grasses have occupied habitats with lower latitudes and elevations, characterized by warmer, sunnier, drier and less fertile environmental conditions. Grass leaf size correlates positively with mean annual temperature and precipitation as expected. Our results also demonstrate that the mean temperature of the wettest quarter of the year is the primary control for C3 leaf size, whereas C4 leaf size is negatively correlated with the difference between summer and winter temperatures. For C4 grasses, phylogeny exerts a significant effect on leaf size but is less important than environmental factors. Our findings highlight the importance of evolutionarily contrasting variations in leaf size between C3 and C4 grasses for shaping their geographical distribution and habitat suitability at the global scale.

Behavioral and phylogenetic correlates of limb length proportions in extant apes and monkeys: Implications for interpreting hominin fossils

The body proportions of extant animals help inform inferences about the behaviors of their extinct relatives, but relationships between body proportions, behavior, and phylogeny in extant primates remain unclear. Advances in behavioral data, molecular phylogenies, and multivariate analytical tools make it an opportune time to perform comprehensive comparative analyses of primate traditional limb length proportions (e.g., intermembral, humerofemoral, brachial, and crural indices), body size-adjusted long bone proportions, and principal components. In this study we used a mix of newly-collected and published data to investigate whether and how the limb length proportions of a diverse sample of primates, including monkeys, apes, and modern humans, are influenced by behavior and phylogeny. We reconfirm that the intermembral index, followed by the first principal component of traditional limb length proportions, is the single most effective variable distinguishing hominoids and other anthropoids. Combined limb length proportions and positional behaviors are strongly correlated in extant anthropoid groups, but phylogeny is a better predictor of limb length proportion variation than of behavior. We confirm convergences between members of the Atelidae and extant apes (especially Pan), members of the Hylobatidae and Pongo, and a potential divergence of Presbytis limb proportions from some other cercopithecoids, which correlate with adaptations for forelimb-dominated behaviors in some colobines. Collectively, these results substantiate hypotheses indicating that extinct hominins and other hominoid taxa can be distinguished by analyzing combinations of their limb length proportions at different taxonomic levels. From these results, we hypothesize that fossil skeletons characterized by notably disparate limb length proportions are unlikely to have exhibited similar behavioral patterns.

Evolutionary history shapes variation of wood density of tree species across the world

The effect of evolutionary history on wood density variation may play an important role in shaping variation in wood density, but this has largely not been tested. Using a comprehensive global dataset including 27,297 measurements of wood density from 2621 tree species worldwide, we test the hypothesis that the legacy of evolutionary history plays an important role in driving the variation of wood density among tree species. We assessed phylogenetic signal in different taxonomic (e.g., angiosperms and gymnosperms) and ecological (e.g., tropical, temperate, and boreal) groups of tree species, explored the biogeographical and phylogenetic patterns of wood density, and quantified the relative importance of current environmental factors (e.g., climatic and soil variables) and evolutionary history (i.e., phylogenetic relatedness among species and lineages) in driving global wood density variation. We found that wood density displayed a significant phylogenetic signal. Wood density differed among different biomes and climatic zones, with higher mean values of wood density in relatively drier regions (highest in subtropical desert). Our study revealed that at a global scale, for angiosperms and gymnosperms combined, phylogeny and species (representing the variance explained by taxonomy and not direct explained by long-term evolution process) explained 84.3% and 7.7% of total wood density variation, respectively, whereas current environment explained 2.7% of total wood density variation when phylogeny and species were taken into account. When angiosperms and gymnosperms were considered separately, the three proportions of explained variation are, respectively, 84.2%, 7.5% and 6.7% for angiosperms, and 45.7%, 21.3% and 18.6% for gymnosperms. Our study shows that evolutionary history outpaced current environmental factors in shaping global variation in wood density.

Linking fine root lifespan to root chemical and morphological traits-A global analysis

Fine root lifespan is a critical trait associated with contrasting root strategies of resource acquisition and protection. Yet, its position within the multidimensional "root economics space" synthesizing global root economics strategies is largely uncertain, and it is rarely represented in frameworks integrating plant trait variations. Here, we compiled the most comprehensive dataset of absorptive median root lifespan (MRL) data including 98 observations from 79 woody species using (mini-)rhizotrons across 40 sites and linked MRL to other plant traits to address questions of the regulators of MRL at large spatial scales. We demonstrate that MRL not only decreases with plant investment in root nitrogen (associated with more metabolically active tissues) but also increases with construction of larger diameter roots which is often associated with greater plant reliance on mycorrhizal symbionts. Although theories linking organ structure and function suggest that root traits should play a role in modulating MRL, we found no correlation between root traits associated with structural defense (root tissue density and specific root length) and MRL. Moreover, fine root and leaf lifespan were globally unrelated, except among evergreen species, suggesting contrasting evolutionary selection between leaves and roots facing contrasting environmental influences above vs. belowground. At large geographic scales, MRL was typically longer at sites with lower mean annual temperature and higher mean annual precipitation. Overall, this synthesis uncovered several key ecophysiological covariates and environmental drivers of MRL, highlighting broad avenues for accurate parametrization of global biogeochemical models and the understanding of ecosystem response to global climate change.

Fungal traits help to understand the decomposition of simple and complex plant litter

Litter decomposition is a key ecosystem process, relevant for the release and storage of nutrients and carbon in soil. Soil fungi are one of the dominant drivers of organic matter decomposition, but fungal taxa differ substantially in their functional ability to decompose plant litter. Knowledge is mostly based on observational data and subsequent molecular analyses and in vitro studies have been limited to forest ecosystems. In order to better understand functional traits of saprotrophic soil fungi in grassland ecosystems, we isolated 31 fungi from a natural grassland and performed several in vitro studies testing for i) leaf and wood litter decomposition, ii) the ability to use carbon sources of differing complexity, iii) the enzyme repertoire. Decomposition strongly varied among phyla and isolates, with Ascomycota decomposing the most and Mucoromycota decomposing the least. The phylogeny of the fungi and their ability to use complex carbon were the most important predictors for decomposition. Our findings show that it is crucial to understand the role of individual members and functional groups within the microbial community. This is an important way forward to understand the role of microbial community composition for the prediction of litter decomposition and subsequent potential carbon storage in grassland soils.

Phylogenetic conservation in plant phenological traits varies between temperate and subtropical climates in China

Phenological traits, such as leaf and flowering dates, are proven to be phylogenetically conserved. The relationship between phylogenetic conservation, plant phenology, and climatic factors remains unknown. Here, we assessed phenological features among flowering plants as evidence for phylogenetic conservatism, the tendency for closely related species to share similar ecological and biological attributes. We use spring phenological traits data from 1968-2018 of 65 trees and 49 shrubs in Xi'an (temperate climate) and Guiyang (subtropical climate) to understand plant phenological traits' relationship with phylogeny. Molecular datasets are employed in evolutionary models to test the phylogenetic conservatism in spring phenological characteristics in response to climate-sensitive phenological features. Significant phylogenetic conservation was found in the Xi'an plant's phenological traits, while there was a non-significant conservation in the Guiyang plant species. Phylogenetic generalized least squares (PGLS) models correlate with phenological features significantly in Xi'an while non-significantly in Guiyang. Based on the findings of molecular dating, it was suggested that the Guiyang species split off from their relatives around 46.0 mya during the middle Eocene of the Tertiary Cenozoic Era, while Xi'an species showed a long evolutionary history and diverged from their relatives around 95 mya during the late Cretaceous Mesozoic Era. First leaf dates (FLD) indicative of spring phenology, show that Xi'an adjourned the case later than Guiyang. Unlike FLD, first flower dates (FFD) yield different results as Guiyang flowers appear later than Xi'an's. Our research revealed that various factors, including phylogeny, growth form, and functional features, influenced the diversity of flowering phenology within species in conjunction with local climate circumstances. These results are conducive to understanding evolutionary conservation mechanisms in plant phenology concerning evolutionary processes in different geographical and climate zones.

Foliar spectra accurately distinguish most temperate tree species and show strong phylogenetic signal

PREMISE

Spectroscopy is a powerful remote sensing tool for monitoring plant biodiversity over broad geographic areas. Increasing evidence suggests that foliar spectral reflectance can be used to identify trees at the species level. However, most studies have focused on only a limited number of species at a time, and few studies have explored the underlying phylogenetic structure of leaf spectra. Accurate species identifications are important for reliable estimations of biodiversity from spectral data.

METHODS

Using over 3500 leaf-level spectral measurements, we evaluated whether foliar reflectance spectra (400-2400 nm) can accurately differentiate most tree species from a regional species pool in eastern North America. We explored relationships between spectral, phylogenetic, and leaf functional trait variation as well as their influence on species classification using a hurdle regression model.

RESULTS

Spectral reflectance accurately differentiated tree species (κ = 0.736, ±0.005). Foliar spectra showed strong phylogenetic signal, and classification errors from foliar spectra, although present at higher taxonomic levels, were found predominantly between closely related species, often of the same genus. In addition, we find functional and phylogenetic distance broadly control the occurrence and frequency of spectral classification mistakes among species.

CONCLUSIONS

Our results further support the link between leaf spectral diversity, taxonomic hierarchy, and phylogenetic and functional diversity, and highlight the potential of spectroscopy to remotely sense plant biodiversity and vegetation response to global change.

Optimum ratio of dietary protein and carbohydrate that maximises lifespan is shared among related insect species

Animals often regulate the intake and quantity of nutrients to maximise fitness through life-history traits such as lifespan, but we still lack a proper understanding of how specific nutrients influence these traits. Here, I developed an algorithm which allowed me to create a nutrient-specific database from literature data, and investigated how the requirements of protein (P) and carbohydrate (C) needed to maximise lifespan evolved across nine insect species. I found moderate evidence of a phylogenetic signal on the optimal ratio of protein to carbohydrate ratio (PC ratio) that maximised lifespan, suggesting that optimal PC ratio for lifespan could have evolved non-independently among related species. I also found evidence for weak-to-strong sex-specific optimal PC ratios for lifespan, suggesting that sex-specific nutritional needs to maximise lifespan can emerge and persist in some species. Although limited in the number of species, the approach adopted here is portable to experiments withn number of nutrients and, thus, can be used in complex comparative precision nutrition studies for insights into the evolution of animal nutrition.

Ulnar shape of extant primates: Functional signals and covariation with triquetrum shape

OBJECTIVES

In this study, we investigated the shape differences of the distal ulna in a phylogenetic context among a broad range of primate taxa. Furthermore, we evaluated covariation between ulnar and triquetrum shape and a possible association between ulnar shape and locomotor behavior.

MATERIALS AND METHODS

We applied 3D geometric morphometrics on a large dataset comprising the distal ulna of 124 anthropoid primate specimens belonging to 12 different genera. For each species, a mean shape was calculated using 11 Procrustes-aligned surface landmarks on the distal ulna. These mean shapes are used in a bgPCA, pPCA, and PACA and 3D morphs were used to visualize more subtle differences between taxa. A p2B-PLS analysis was performed to test the covariance between distal ulnar and triquetrum shape.

RESULTS

The results show that more closely related species exhibit a similar distal ulnar shape. Overall, extant hominid ulnae show a shape shift compared to those of extant monkeys and hylobatids. This includes a shortening of the ulnar styloid process and dorspalmarly widening of the ulnar head, shape characteristics that are independent of phylogeny. Within the hominids, Pongo pygmaeus seem to possess the most plesiomorphic distal ulnar shape, while Gorilla and Homo sapiens display the most derived distal ulna. Cercopithecoids, hylobatids, and P. pygmaeus are characterized by a relatively deep ECU groove, which is a shape trait dependent of phylogeny. Although there was no significant covariation between distal ulnar shape and triquetrum shape, the shape differences of the distal ulna between the different primate taxa reveal a possible link with locomotor behavior.

CONCLUSIONS

The comparative analyses of this study reveal different shape trends in a phylogenetic context. Highly arboreal primates, such as hylobatids and Ateles fusciceps, show a distal ulnar morphology that appears to be adapted to tensile and torsional forces. In primates that use their wrist under more compressive conditions, such as quadrupedal cercopithecoids and great apes, the distal ulnar morphology seems to reflect increased compressive forces. In modern humans, the distal ulnar shape can be associated to enhanced manipulative skills and power grips. There was no significant covariation between distal ulnar shape and triquetrum shape, probably due to the variation in the amount of contact between the triquetrum and ulna. In combination with future research on wrist mobility in diverse primate taxa, the results of this study will allow us to establish form-function relationships of the primate wrist and contribute towards an evidence-based interpretation of fossil remains.

Convergent and environmentally associated chromatic polymorphism in Bryconops Kner, 1858 (Ostariophysi: Characiformes: Iguanodectidae)

Bryconops Kner, 1858, includes two well defined subgenera based on morphological evidence, with each containing at least one species (B. (Bryconops) caudomaculatus and B. (Creatochanes) melanurus) with a very wide distribution, within which regional populations present color variations. To test if phenotypic variation is related to cladogenetic events, we performed tests for phylogenetic independence and determined the strength of convergence for color characters in relation to water type, as the variation between clear, black and white waters is considered to be one of the major driving forces in the evolution of Amazonian fishes. Color characters for fins above the median line of the body were generally found to be independent from phylogeny and the Wheatsheaf test strongly supports convergence of the dorsal fin color between populations of species in the same type of water, with a similar trend suggested for the color of the dorsal lobe of the caudal fin. This means that simple color characters cannot necessarily be relied upon for taxonomic revisions of the genus as local phenotypic variants may represent environmentally determined plasticity or convergent evolution. Further studies are required to determine the validity of these characters.

Investigating the evolutionary dynamics of diploid number variation in Ctenomys (Ctenomyidae, Rodentia)

Contrary to predictions from classical hybrid sterility models of chromosomal speciation, some organisms display high rates of karyotype variation. Ctenomys are the current mammals with the greatest interspecific and intraspecific chromosomal variation. A large number of species have been studied cytogenetically. The diploid numbers of chromosomes range from 2n = 10 to 2n = 70. Here, we analyzed karyotype evolution in Ctenomys using comparative phylogenetic methods. We found a strong phylogenetic signal with chromosome number. This refutes the chromosomal megaevolution model, which proposes the independent accumulation of multiple chromosomal rearrangements in each closely related species. We found that Brownian motion (BM) described the observed characteristic changes more thoroughly than the Ornstein-Uhlenbeck and Early-Burst models. This suggests that the evolution of chromosome numbers occurs by a random walk along phylogenetic clades. However, our data indicate that the BM model alone does not fully characterize the chromosomal evolution of Ctenomys.

A phylogenomic perspective on interspecific competition

Evolutionary processes may have substantial impacts on community assembly, but evidence for phylogenetic relatedness as a determinant of interspecific interaction strength remains mixed. In this perspective, we consider a possible role for discordance between gene trees and species trees in the interpretation of phylogenetic signal in studies of community ecology. Modern genomic data show that the evolutionary histories of many taxa are better described by a patchwork of histories that vary along the genome rather than a single species tree. If a subset of genomic loci harbour trait-related genetic variation, then the phylogeny at these loci may be more informative of interspecific trait differences than the genome background. We develop a simple method to detect loci harbouring phylogenetic signal and demonstrate its application through a proof-of-principle analysis of Penicillium genomes and pairwise interaction strength. Our results show that phylogenetic signal that may be masked genome-wide could be detectable using phylogenomic techniques and may provide a window into the genetic basis for interspecific interactions.

Machine learning and phylogenetic analysis allow for predicting antibiotic resistance in M. tuberculosis

BACKGROUND

Antimicrobial resistance (AMR) poses a significant global health threat, and an accurate prediction of bacterial resistance patterns is critical for effective treatment and control strategies. In recent years, machine learning (ML) approaches have emerged as powerful tools for analyzing large-scale bacterial AMR data. However, ML methods often ignore evolutionary relationships among bacterial strains, which can greatly impact performance of the ML methods, especially if resistance-associated features are attempted to be detected. Genome-wide association studies (GWAS) methods like linear mixed models accounts for the evolutionary relationships in bacteria, but they uncover only highly significant variants which have already been reported in literature.

RESULTS

In this work, we introduce a novel phylogeny-related parallelism score (PRPS), which measures whether a certain feature is correlated with the population structure of a set of samples. We demonstrate that PRPS can be used, in combination with SVM- and random forest-based models, to reduce the number of features in the analysis, while simultaneously increasing models' performance. We applied our pipeline to publicly available AMR data from PATRIC database for Mycobacterium tuberculosis against six common antibiotics.

CONCLUSIONS

Using our pipeline, we re-discovered known resistance-associated mutations as well as new candidate mutations which can be related to resistance and not previously reported in the literature. We demonstrated that taking into account phylogenetic relationships not only improves the model performance, but also yields more biologically relevant predicted most contributing resistance markers.

Genome-scale phylogeny and comparative genomics of the fungal order Sordariales

The order Sordariales is taxonomically diverse, and harbours many species with different lifestyles and large economic importance. Despite its importance, a robust genome-scale phylogeny, and associated comparative genomic analysis of the order is lacking. In this study, we examined whole-genome data from 99 Sordariales, including 52 newly sequenced genomes, and seven outgroup taxa. We inferred a comprehensive phylogeny that resolved several contentious relationships amongst families in the order, and cleared-up intrafamily relationships within the Podosporaceae. Extensive comparative genomics showed that genomes from the three largest families in the dataset (Chaetomiaceae, Podosporaceae and Sordariaceae) differ greatly in GC content, genome size, gene number, repeat percentage, evolutionary rate, and genome content affected by repeat-induced point mutations (RIP). All genomic traits showed phylogenetic signal, and ancestral state reconstruction revealed that the variation of the properties stems primarily from within-family evolution. Together, the results provide a thorough framework for understanding genome evolution in this important group of fungi.

Evaluating the Performance of Widely Used Phylogenetic Models for Gene Expression Evolution

Phylogenetic comparative methods are increasingly used to test hypotheses about the evolutionary processes that drive divergence in gene expression among species. However, it is unknown whether the distributional assumptions of phylogenetic models designed for quantitative phenotypic traits are realistic for expression data and importantly, the reliability of conclusions of phylogenetic comparative studies of gene expression may depend on whether the data is well described by the chosen model. To evaluate this, we first fit several phylogenetic models of trait evolution to 8 previously published comparative expression datasets, comprising a total of 54,774 genes with 145,927 unique gene-tissue combinations. Using a previously developed approach, we then assessed how well the best model of the set described the data in an absolute (not just relative) sense. First, we find that Ornstein-Uhlenbeck models, in which expression values are constrained around an optimum, were the preferred models for 66% of gene-tissue combinations. Second, we find that for 61% of gene-tissue combinations, the best-fit model of the set was found to perform well; the rest were found to be performing poorly by at least one of the test statistics we examined. Third, we find that when simple models do not perform well, this appears to be typically a consequence of failing to fully account for heterogeneity in the rate of the evolution. We advocate that assessment of model performance should become a routine component of phylogenetic comparative expression studies; doing so can improve the reliability of inferences and inspire the development of novel models.

Variability in seed salinity tolerance in an island coastal community

BACKGROUND AND AIMS

Islands, with their long coastlines and increased vulnerability to sea level rise, offer compelling opportunities to investigate the salinity tolerance of coastal plants. Seeds are generally more vulnerable than other plant stages to increased stressors. The aim of this study was to characterize salinity tolerance during germination across a diverse pool of 21 species from 14 plant families found in coastal communities throughout the Hawaiian Islands in order to increase our general understanding of coastal plant ecology for conservation and restoration.

METHODS

Seeds of each species were exposed to unfiltered/untreated seawater (35 ppt total salinity) and two salinity treatments (10 and 20 ppt) in which the seawater was diluted with distilled water, and germination percent and timing were compared to seeds in a distilled water control. Non-germinated seeds were then tested for recovery germination. We quantified and compared germination percent, time and recovery among species and across salinity levels and tested for heterogeneity related to seed size, dormancy class, habit and threatened status.

KEY RESULTS

Although salinity tolerance varied considerably among species, salinity exposure generally reduced and delayed germination. The greatest effects were detected at higher salinity levels. Recovery germination overall was higher for seeds that had been exposed to higher salinity. None of the factors we explored emerged as predictors of salinity tolerance except seed mass, which tended to enhance germination at higher salinity.

CONCLUSIONS

Species responses to salinity exposure indicate high vulnerability of coastal systems to increased salinity stress, and variability among species could lead to shifts in community assembly and composition under sea level rise. These results can help guide coastal ecosystem conservation and restoration management decisions in the face of climate change.

Foliar elementome and functional traits relationships identify tree species niche in French Guiana rainforests

Biogeochemical niche (BN) hypothesis aims to relate species/genotype elemental composition with its niche based on the fact that different elements are involved differentially in distinct plant functions. We here test the BN hypothesis through the analysis of the 10 foliar elemental concentrations and 20 functional-morphological of 60 tree species in a French Guiana tropical forest. We observed strong legacy (phylogenic + species) signals in the species-specific foliar elemental composition (elementome) and, for the first time, provide empirical evidence for a relationship between species-specific foliar elementome and functional traits. Our study thus supports the BN hypothesis and confirms the general niche segregation process through which the species-specific use of bio-elements drives the high levels of α-diversity in this tropical forest. We show that the simple analysis of foliar elementomes may be used to test for BNs of co-occurring species in highly diverse ecosystems, such as tropical rainforests. Although cause and effect mechanisms of leaf functional and morphological traits in species-specific use of bio-elements require confirmation, we posit the hypothesis that divergences in functional-morphological niches and species-specific biogeochemical use are likely to have co-evolved.

The reconstructed cranium of Pierolapithecus and the evolution of the great ape face

Pierolapithecus catalaunicus (~12 million years ago, northeastern Spain) is key to understanding the mosaic nature of hominid (great ape and human) evolution. Notably, its skeleton indicates that an orthograde (upright) body plan preceded suspensory adaptations in hominid evolution. However, there is ongoing debate about this species, partly because the sole known cranium, preserving a nearly complete face, suffers from taphonomic damage. We 1) carried out a micro computerized tomography (CT) based virtual reconstruction of the Pierolapithecus cranium, 2) assessed its morphological affinities using a series of two-dimensional (2D) and three-dimensional (3D) morphometric analyses, and 3) modeled the evolution of key aspects of ape face form. The reconstruction clarifies many aspects of the facial morphology of Pierolapithecus. Our results indicate that it is most similar to great apes (fossil and extant) in overall face shape and size and is morphologically distinct from other Middle Miocene apes. Crown great apes can be distinguished from other taxa in several facial metrics (e.g., low midfacial prognathism, relatively tall faces) and only some of these features are found in Pierolapithecus, which is most consistent with a stem (basal) hominid position. The inferred morphology at all ancestral nodes within the hominoid (ape and human) tree is closer to great apes than to hylobatids (gibbons and siamangs), which are convergent with other smaller anthropoids. Our analyses support a hominid ancestor that was distinct from all extant and fossil hominids in overall facial shape and shared many features with Pierolapithecus. This reconstructed ancestral morphotype represents a testable hypothesis that can be reevaluated as new fossils are discovered.

Plant parentage influences the type of timber use by traditional peoples of the Brazilian Caatinga

Local populations select different plants to meet their demands, so that morphologically similar species can be more used for a given use. Herein, we seek to understand whether plant species that are phylogenetically closer together are used more similarly than distant species in the phylogeny. Ethnobotanical data were collected in five rural communities in a semi-arid region of Brazil. A total of 120 local experts were selected and interviewed using semi-structured questionnaires. The people's knowledge of plants was organized into usage subcategories. We estimated the redundancy values for the mentioned species, and we compiled data from the literature on the wood density values of the cited species. We constructed our phylogenetic hypothesis of useful plants and used comparative phylogenetic methods to estimate the signal. Our results showed a strong phylogenetic grouping for both tool handle and craft uses. We observed a moderate phylogenetic grouping in which related cited plants exhibit similar redundancy and a weak grouping in which cited plants present similar wood density values. Our results revealed the importance of using phylogeny for useful plants. We conclude the phylogenetic proximity of useful plants and the lower redundancy for some species in our study may suggest greater use pressure, given that few species fulfill the same function.

Protective Geometry and Reproductive Anatomy as Candidate Determinants of Clutch Size Variation in Pentatomid Bugs

AbstractMany animals lay their eggs in clusters. Eggs on the periphery of clusters can be at higher risk of mortality. We asked whether the most commonly occurring clutch sizes in pentatomid bugs could result from geometrical arrangements that maximize the proportion of eggs in the cluster's interior. Although the most common clutch sizes do not correspond with geometric optimality, stink bugs do tend to lay clusters of eggs in shapes that protect increasing proportions of their offspring as clutch sizes increase. We also considered whether ovariole number, an aspect of reproductive anatomy that may be a fixed trait across many pentatomids, could explain observed distributions of clutch sizes. The most common clutch sizes across many species correspond with multiples of ovariole number. However, there are species with the same number of ovarioles that lay clutches of widely varying size, among which multiples of ovariole number are not overrepresented. In pentatomid bugs, reproductive anatomy appears to be more important than egg mass geometry in determining clutch size uniformity. In addition, our analysis demonstrates that groups of animals with little variation in ovariole number may nonetheless lay a broad range of clutch shapes and sizes.

Consistent changes in muscle metabolism underlie dive performance across multiple lineages of diving ducks

Diving animals must sustain high activity with limited O2 stores to successfully capture prey. Studies suggest that increasing body O2 stores supports breath-hold diving, but less is known about metabolic specializations that underlie underwater locomotion. We measured maximal activities of 10 key enzymes in locomotory muscles (gastrocnemius and pectoralis) to identify biochemical changes associated with diving in pathways of oxidative and substrate-level phosphorylation and compared them across three groups of ducks-the longest diving sea ducks (eight spp.), the mid-tier diving pochards (three spp.) and the non-diving dabblers (five spp.). Relative to dabblers, both diving groups had increased activities of succinate dehydrogenase and cytochrome c oxidase, and sea ducks further showed increases in citrate synthase (CS) and hydroxyacyl-CoA dehydrogenase (HOAD). Both diving groups had relative decreases in capacity for anaerobic metabolism (lower ratio of lactate dehydrogenase to CS), with sea ducks also showing a greater capacity for oxidative phosphorylation and lipid oxidation (lower ratio of pyruvate kinase to CS, higher ratio of HOAD to hexokinase). These data suggest that the locomotory muscles of diving ducks are specialized for sustaining high rates of aerobic metabolism, emphasizing the importance of body O2 stores for dive performance in these species.

Community-wide patterns in pollen and ovule production, their ratio (P/O), and other floral traits along an elevation gradient in southwestern China

BACKGROUND

As the male and female gametophytes of flowering plants, pollen and ovules largely determine the upper and lower boundaries of plant reproductive success. It is commonly predicted that pollen and ovule number per flower should increase, and pollen-ovule ratio (P/O) per flower should decrease with increasing elevation in response to a more stochastic pollination environment. Here, we aimed to determine the response of pollen number, ovule number, and P/O to other floral traits and elevation gradients for 84 insect-pollinated herbaceous flowering plant species in five sub-alpine and alpine communities (2709 to 3896 m a.s.l.) on Yulong Snow Mountain, southwestern China.

RESULTS

Six floral traits, including P/O, floral display area, flower number, tube depth, flower shape, and pollen presentation, were highly correlated with pollen and ovule number per flower. With increasing elevation, pollen number and P/O per flower increased marginally and significantly, respectively; ovule number per individual, flower number per individual, stigma stamen separation, and inflorescence height decreased significantly. However, ovule number per flower and other floral traits (i.e., floral display area, tube depth, stigma height, stamen height, and pollen and P/O per individual) did not change with elevation. We detected significant phylogenetic signals for pollen number, ovule number, and P/O, suggesting that these traits may be highly conserved and with limited response to changing environmental conditions.

CONCLUSIONS

Results revealed patterns of plant reproductive character evolution along elevation gradients and the potential factors governing their spatial variation in high-elevation environments. Plant species at high elevations are more likely adapted to cross-pollination, indicated by increased P/O per flower at high elevations on Yulong Mountain. Combined effects of phylogenetic history and plant-pollinator interactions should determine plant trait evolution.

Mixed effects of a national protected area network on terrestrial and freshwater biodiversity

Protected areas are considered fundamental to counter biodiversity loss. However, evidence for their effectiveness in averting local extinctions remains scarce and taxonomically biased. We employ a robust counterfactual multi-taxon approach to compare occupancy patterns of 638 species, including birds (150), mammals (23), plants (39) and phytoplankton (426) between protected and unprotected sites across four decades in Finland. We find mixed impacts of protected areas, with only a small proportion of species explicitly benefiting from protection-mainly through slower rates of decline inside protected areas. The benefits of protection are enhanced for larger protected areas and are traceable to when the sites were protected, but are mostly unrelated to species conservation status or traits (size, climatic niche and threat status). Our results suggest that the current protected area network can partly contribute to slow down declines in occupancy rates, but alone will not suffice to halt the biodiversity crisis. Efforts aimed at improving coverage, connectivity and management will be key to enhance the effectiveness of protected areas towards bending the curve of biodiversity loss.

Predictability of the community-function landscape in wine yeast ecosystems

Predictively linking taxonomic composition and quantitative ecosystem functions is a major aspiration in microbial ecology, which must be resolved if we wish to engineer microbial consortia. Here, we have addressed this open question for an ecological function of major biotechnological relevance: alcoholic fermentation in wine yeast communities. By exhaustively phenotyping an extensive collection of naturally occurring wine yeast strains, we find that most ecologically and industrially relevant traits exhibit phylogenetic signal, allowing functional traits in wine yeast communities to be predicted from taxonomy. Furthermore, we demonstrate that the quantitative contributions of individual wine yeast strains to the function of complex communities followed simple quantitative rules. These regularities can be integrated to quantitatively predict the function of newly assembled consortia. Besides addressing theoretical questions in functional ecology, our results and methodologies can provide a blueprint for rationally managing microbial processes of biotechnological relevance.

Cell wall polysaccharide and glycoprotein content tracks growth-form diversity and an aridity gradient in the leaf-succulent genus Crassula

Cell wall traits are believed to be a key component of the succulent syndrome, an adaptive syndrome to drought, yet the variability of such traits remains largely unknown. In this study, we surveyed the leaf polysaccharide and glycoprotein composition in a wide sampling of Crassula species that occur naturally along an aridity gradient in southern Africa, and we interpreted its adaptive significance in relation to growth form and arid adaptation. To study the glycomic diversity, we sampled leaf material from 56 Crassula taxa and performed comprehensive microarray polymer profiling to obtain the relative content of cell wall polysaccharides and glycoproteins. This analysis was complemented by the determination of monosaccharide composition and immunolocalization in leaf sections using glycan-targeting antibodies. We found that compact and non-compact Crassula species occupy distinct phenotypic spaces in terms of leaf glycomics, particularly in regard to rhamnogalacturonan I, its arabinan side chains, and arabinogalactan proteins (AGPs). Moreover, these cell wall components also correlated positively with increasing aridity, which suggests that they are likely advantageous in terms of arid adaptation. These differences point to compact Crassula species having more elastic cell walls with plasticizing properties, which can be interpreted as an adaptation toward increased drought resistance. Furthermore, we report an intracellular pool of AGPs associated with oil bodies and calcium oxalate crystals, which could be a peculiarity of Crassula and could be linked to increased drought resistance. Our results indicate that glycomics may be underlying arid adaptation and drought resistance in succulent plants.

On the Origins of Phenotypic Parallelism in Benthic and Limnetic Stickleback

Rapid evolution of similar phenotypes in similar environments, giving rise to in situ parallel adaptation, is an important hallmark of ecological speciation. However, what appears to be in situ adaptation can also arise by dispersal of divergent lineages from elsewhere. We test whether two contrasting phenotypes repeatedly evolved in parallel, or have a single origin, in an archetypal example of ecological adaptive radiation: benthic-limnetic three-spined stickleback (Gasterosteus aculeatus) across species pair and solitary lakes in British Columbia. We identify two genomic clusters across freshwater populations, which differ in benthic-limnetic divergent phenotypic traits and separate benthic from limnetic individuals in species pair lakes. Phylogenetic reconstruction and niche evolution modeling both suggest a single evolutionary origin for each of these clusters. We detected strong phylogenetic signal in benthic-limnetic divergent traits, suggesting that they are ancestrally retained. Accounting for ancestral state retention, we identify local adaptation of body armor due to the presence of an intraguild predator, the sculpin (Cottus asper), and environmental effects of lake depth and pH on body size. Taken together, our results imply a predominant role for retention of ancestral characteristics in driving trait distribution, with further selection imposed on some traits by environmental factors.

Variation in enamel mechanical properties throughout the crown in catarrhine primates

Enamel mechanical properties vary across molar crowns, but the relationship among mechanical properties, tooth function, and phylogeny are not well understood. Fifteen primate lower molars representing fourteen taxa (catarrhine, n = 13; platyrrhine, n = 1) were sectioned in the lingual-buccal plane through the mesial cusps. Gradients of enamel mechanical properties, specifically hardness and elastic modulus, were quantified using nanoindentation from inner (near the enamel-dentine junction), through middle, to outer enamel (near the outer enamel surface) at five positions (buccal lateral, buccal cuspal, occlusal middle, lingual cuspal, lingual lateral). Cuspal positions had higher mechanical property values than lateral positions. Middle enamel had higher mean hardness and elastic modulus values than inner and outer locations in all five crown positions. Functionally, the thicker-enameled buccal cusps of lower molars did not show evidence of increased resistance to failure; instead, lingual cusps-which show higher rates of fracture-had higher average mechanical property values, with no significant differences observed between sides. Preliminary phylogenetic results suggest there is relatively little phylogenetic signal in gradients of mechanical properties through the enamel or across the crown. There appears to be common mechanical property patterns across molar crowns in Catarrhini and potentially among primates more broadly. These results may allow more precise interpretations of dental biomechanics and processes resulting in mechanical failure of enamel in primates, such as wear and fracture.

Evaluating the Performance of Widely Used Phylogenetic Models for Gene Expression Evolution

Phylogenetic comparative methods are increasingly used to test hypotheses about the evolutionary processes that drive divergence in gene expression among species. However, it is unknown whether the distributional assumptions of phylogenetic models designed for quantitative phenotypic traits are realistic for expression data and importantly, the reliability of conclusions of phylogenetic comparative studies of gene expression may depend on whether the data is well-described by the chosen model. To evaluate this, we first fit several phylogenetic models of trait evolution to 8 previously published comparative expression datasets, comprising a total of 54,774 genes with 145,927 unique gene-tissue combinations. Using a previously developed approach, we then assessed how well the best model of the set described the data in an absolute (not just relative) sense. First, we find that Ornstein-Uhlenbeck models, in which expression values are constrained around an optimum, were the preferred model for 66% of gene-tissue combinations. Second, we find that for 61% of gene-tissue combinations, the best fit model of the set was found to perform well; the rest were found to be performing poorly by at least one of the test statistics we examined. Third, we find that when simple models do not perform well, this appears to be typically a consequence of failing to fully account for heterogeneity in the rate of the evolution. We advocate that assessment of model performance should become a routine component of phylogenetic comparative expression studies; doing so can improve the reliability of inferences and inspire the development of novel models.

Reproductive and environmental traits explain the variation in egg size among Medusozoa (Cnidaria)

Medusozoa (Cnidaria) are characterized by diverse life cycles, with different semaphoronts (medusa, medusoid, fixed gonophore, polyp) representing the sexual phase and carrying the gametes. Although egg size is often considered a proxy to understand reproductive and developmental traits of medusozoans, understanding of the processes influencing egg size variation in the group under an evolutionary context is still limited. We carried out a comprehensive review of the variation of egg size in Medusozoa to test whether this variation is related to biological/sexual or environmental traits. Egg size presents a strong phylogenetic signal (λ = 0.79, K = 0.67), explaining why closely related species with different reproductive strategies and different individual sizes have similar egg sizes. However, variation in egg size is influenced by the number of eggs, depth and temperature, with larger eggs frequently present in species with few eggs (1-15), in deep-sea species and in cold-water species. Conversely, the production of small eggs among cold-water species of Staurozoa might be associated with the development of a small benthic larvae in this group. Our study reinforces that egg sizes respond to reproductive and environmental traits, although egg size is highly conserved within medusa classes.

Response to Tcherkez and Farquhar: Rubisco adaptation is more limited by phylogenetic constraint than by catalytic trade-off

Rubisco is the primary entry point for carbon into the biosphere. It has been widely proposed that rubisco is highly constrained by catalytic trade-offs due to correlations between the enzyme's kinetic traits across species. In previous work, we have shown that the strength of these correlations, and thus the strength of catalytic trade-offs, have been overestimated due to the presence of phylogenetic signal in the kinetic trait data (Bouvier et al., 2021). We demonstrated that only the trade-offs between the Michaelis constant for CO2 and carboxylase turnover, and between the Michaelis constants for CO2 and O2 were robust to phylogenetic effects. We further demonstrated that phylogenetic constraints have limited rubisco adaptation to a greater extent than the combined action of catalytic trade-offs. Recently, however, our claims have been contested by Tcherkez and Farquhar (2021), who have argued that the phylogenetic signal we detect in rubisco kinetic traits is an artefact of species sampling, the use of rbcL-based trees for phylogenetic inference, laboratory-to-laboratory variability in kinetic measurements, and homoplasy of the C4 trait. In the present article, we respond to these criticisms on a point-by-point basis and conclusively show that all are unfounded. As such, we stand by our original conclusions. Namely, although rubisco kinetic evolution has been limited by biochemical trade-offs, these are not absolute and have been previously overestimated due to phylogenetic biases. Instead, rubisco adaptation has in fact been more limited by phylogenetic constraint.

Intrinsic and extrinsic drivers of shape variation in the albatross compound bill

Albatross are the largest seabirds on Earth and have a suite of adaptations for their pelagic lifestyle. Rather than having a bill made of a single piece of keratin, Procellariiformes have a compound rhamphotheca, made of several joined plates. Drivers of the shape of the albatross bill have not been explored. Here we use three-dimensional scans of 61 upper bills from 12 species of albatross to understand whether intrinsic (species assignment & size) or extrinsic (diet) factors predict bill shape. Diet is a significant predictor of bill shape with coarse dietary categories providing higher R2 values than dietary proportion data. We also find that of the intrinsic factors, species assignment accounts for ten times more of the variation than size (72% versus 6.8%) and that there is a common allometric vector of shape change between all species. When considering species averages in a phylogenetic framework, there are significant Blomberg's K results for both shape and size (K = 0.29 & 1.10) with the first axis of variation having a much higher K value (K = 1.9), reflecting the split in shape at the root of the tree. The influence of size on bill shape is limited, with species assignment and diet predicting far more of the variation. The results show that both intrinsic and extrinsic factors are needed to understand morphological evolution.

Seed dispersal syndrome predicts ethanol concentration of fruits in a tropical dry forest

Studying fruit traits and their interactions with seed dispersers can improve how we interpret patterns of biodiversity, ecosystem function and evolution. Mounting evidence suggests that fruit ethanol is common and variable, and may exert selective pressures on seed dispersers. To test this, we comprehensively assess fruit ethanol content in a wild ecosystem and explore sources of variation. We hypothesize that both phylogeny and seed dispersal syndrome explain variation in ethanol levels, and we predict that fruits with mammalian dispersal traits will contain higher levels of ethanol than those with bird dispersal traits. We measured ripe fruit ethanol content in species with mammal- (n = 16), bird- (n = 14) or mixed-dispersal (n = 7) syndromes in a Costa Rican tropical dry forest. Seventy-eight per cent of fruit species yielded measurable ethanol concentrations. We detected a phylogenetic signal in maximum ethanol levels (Pagel's λ = 0.82). Controlling for phylogeny, we observed greater ethanol concentrations in mammal-dispersed fruits, indicating that dispersal syndrome helps explain variation in ethanol content, and that mammals may be more exposed to ethanol in their diets than birds. Our findings further our understanding of wild fruit ethanol and its potential role as a selective pressure on frugivore sensory systems and metabolism.

Repeated evolution of similar phenotypes: Integrating comparative methods with developmental pathways

Repeated phenotypes, often referred to as 'homoplasies' in cladistic analyses, may evolve through changes in developmental processes. Genetic bases of recurrent evolution gained attention and have been studied in the past years using approaches that combine modern analytical phylogenetic tools with the stunning assemblage of new information on developmental mechanisms. In this review, we evaluated the topic under an integrated perspective, revisiting the classical definitions of convergence and parallelism and detailing comparative methods used to evaluate evolution of repeated phenotypes, which include phylogenetic inference, estimates of evolutionary rates and reconstruction of ancestral states. We provide examples to illustrate how a given methodological approach can be used to identify evolutionary patterns and evaluate developmental mechanisms associated with the intermittent expression of a given trait along the phylogeny. Finally, we address why repeated trait loss challenges strict definitions of convergence and parallelism, discussing how changes in developmental pathways might explain the high frequency of repeated trait loss in specific lineages.

Patterns and Factors Influencing Parrot (Order: Psittaciformes) Success in Establishing Thriving Naturalized Populations within the Contiguous United States

Parrots (Order: Psittaciformes) represent one of the most striking and ecomorphologically diverse avian clades, spanning more than two orders of magnitude in body size with populations occupying six continents. The worldwide diaspora of parrots is largely due to the pet trade, driven by human desire for bright, colorful, and intelligent animals as companions. Some introduced species have aptly inserted themselves into the local ecosystem and established successful breeding colonies all around the globe. Notably, the United States is home to several thriving populations of introduced species including red-masked parakeets (Psittacara erythrogenys), monk parakeets (Myiopsitta monachus), nanday conures (Aratinga nenday), and red-crowned amazons (Amazona viridigenalis). Their incredible success globally begs the question as to how these birds adapt so readily to novel environments. In this commentary, we trace parrots through evolutionary history, contextualize existent naturalized parrot populations within the contiguous United States, and provide a phylogenetic regression analysis of body mass and brain size based on success in establishing breeding populations. The propensity for a parrot species to become established appears to be phylogenetically driven. Notably, parrots in the family Cacatuidae and Neotropical Pyrrhua appear to be poor at establishing themselves in the United States once released. Although brain size among Psittaciformes did not show a significant impact on successful breeding in the continental United States, we propose that the success of parrots can be attributed to their charismatic nature, significant intelligence relative to other avian lineages, and behavioral flexibility.

Whole-cell modeling of E. coli colonies enables quantification of single-cell heterogeneity in antibiotic responses

Antibiotic resistance poses mounting risks to human health, as current antibiotics are losing efficacy against increasingly resistant pathogenic bacteria. Of particular concern is the emergence of multidrug-resistant strains, which has been rapid among Gram-negative bacteria such as Escherichia coli. A large body of work has established that antibiotic resistance mechanisms depend on phenotypic heterogeneity, which may be mediated by stochastic expression of antibiotic resistance genes. The link between such molecular-level expression and the population levels that result is complex and multi-scale. Therefore, to better understand antibiotic resistance, what is needed are new mechanistic models that reflect single-cell phenotypic dynamics together with population-level heterogeneity, as an integrated whole. In this work, we sought to bridge single-cell and population-scale modeling by building upon our previous experience in "whole-cell" modeling, an approach which integrates mathematical and mechanistic descriptions of biological processes to recapitulate the experimentally observed behaviors of entire cells. To extend whole-cell modeling to the "whole-colony" scale, we embedded multiple instances of a whole-cell E. coli model within a model of a dynamic spatial environment, allowing us to run large, parallelized simulations on the cloud that contained all the molecular detail of the previous whole-cell model and many interactive effects of a colony growing in a shared environment. The resulting simulations were used to explore the response of E. coli to two antibiotics with different mechanisms of action, tetracycline and ampicillin, enabling us to identify sub-generationally-expressed genes, such as the beta-lactamase ampC, which contributed greatly to dramatic cellular differences in steady-state periplasmic ampicillin and was a significant factor in determining cell survival.

Spatial and temporal changes in the assembly mechanism and co-occurrence network of the chromophytic phytoplankton communities in coastal ecosystems under anthropogenic influences

As a typical semiclosed coastal sea area in China, the ecological environment of Bohai Bay has been significantly disturbed by human activities. As primary producers, the chromophytic phytoplankton are the basis of ecosystems, especially in coastal ecosystems, and changes in the chromophytic phytoplankton community can affect the stability of the entire ecosystem. In this study, we investigated the effects of the human activity-induced spatial and temporal environmental heterogeneity on the community composition, diversity, assembly mechanisms, and co-occurrence networks of chromophytic phytoplankton in Bohai Bay during the wet season and the dry season. The results showed that in both seasons, there was obvious environmental heterogeneity between the nearshore area and the offshore area, and the nearshore areas were more affected by human disturbance. Although higher diversity was supported by the abundance of nutrients in nearshore areas, co-occurrence network analysis revealed that the chromophytic phytoplankton were less closely connected to each other in nearshore areas than in offshore areas due to chemical oxygen demand (COD), eutrophication index (EI), and dissolved inorganic nitrogen (DIN). The nearshore network was less stable than the offshore co-occurrence network in both seasons, which was related to the concentration of dissolved oxygen and COD. Both stochastic and deterministic processes dominated the assembly of the chromophytic phytoplankton communities, with different importance rankings of stochastic and deterministic processes in the nearshore and offshore areas. Drift dominated the assembly of the communities in nearshore areas, while variable selection dominated the assembly of the communities in offshore areas. DIN, EI, and COD, rather than geographic distance, were the main environmental factors affecting the phylogenetic turnover of the chromophytic phytoplankton. Our study showed that environmental heterogeneity caused by human disturbance had a greater impact on the chromophytic phytoplankton communities in Bohai Bay than natural factors such as temperature and salinity.