Global leaf sulfur stoichiometry and the relationships with nitrogen and phosphorus: phylogeny, growth form and environmental controls.

Sulfur (S) is an essential bioelement with vital roles in serving regulatory and catalytic functions and tightly coupled with N and P in plants. However, globally stoichiometric patterns of leaf S and its relationships to leaf N and P are less well studied. We compiled 31 939 records of leaf-based data for 2600 plant species across 6652 sites worldwide. All plant species were divided into different phylogenetic taxa and growth forms. Standard major axis analysis was employed to fit the bivariate element relationships. A phylogenetic linear mixed-effect model and a multiple-regression model were used to partition the variations of bioelements into phylogeny and environments, and then to estimate the importance of environmental variables. Global geometric mean leaf S, N and P concentrations were 1.44, 15.70 and 1.27 mg g-1, respectively, with significant differences among plant groups. Leaf S-N-P positively correlated with each other, ignoring plant groups. The scaling exponents of LN-LS, LP-LS and LN-LP were 0.64, 0.76 and 0.79, respectively, for all species, but differed among plant groups. Both phylogeny and environments regulated the bioelements. The variability, rather than mean temperature, controlled the bioelements. Phylogeny explained more for the concentrations of all the three bioelements than environments, of which S was the one most affected by phylogenetic taxa.

Allometry and phylogenetic divergence: Correspondence or incongruence?

The potential connection between trends of within species variation, such as those of allometric change in morphology, and phylogenetic divergence has been a central topic in evolutionary biology for more than a century, including in the context of human evolution. In this study, I focus on size-related shape change in craniofacial proportions using a sample of more than 3200 adult Old World monkeys belonging to 78 species, of which 2942 specimens of 51 species are selected for the analysis. Using geometric morphometrics, I assess whether the divergence in the direction of static allometries increases in relation to phyletic differences. Because both small samples and taxonomic sampling may bias the results, I explore the sensitivity of the main analyses to the inclusion of more or less taxa depending on the choice of a threshold for the minimum sample size of a species. To better understand the impact of sampling error, I also use randomized subsampling experiments in the largest species samples. The study shows that static allometries vary broadly in directions without any evident phylogenetic signal. This variation is much larger than previously found in ontogenetic trajectories of Old World monkeys, but the conclusion of no congruence with phylogenetic divergence is the same. Yet, the effect of sampling error clearly contributes to inaccuracies and tends to magnify the differences in allometric change. Thus, morphometric research at the boundary between micro- and macro-evolution in primates, and more generally in mammals, critically needs very large and representative samples. Besides sampling error, I suggest other non-mutually exclusive explanations for the lack of correspondence between allometric and phylogenetic divergence in Old World monkeys, and also discuss why directions might be more variable in static compared to ontogenetic trajectories. Even if allometric variation may be a poor source of information in relation to phylogeny, the evolution of allometry is a fascinating subject and the study of size-related shape changes remains a fundamental piece of the puzzle to understand morphological variation within and between species in primates and other animals.

Variation of carbon, nitrogen and phosphorus content in fungi reflects their ecology and phylogeny.

Fungi are an integral part of the nitrogen and phosphorus cycling in trophic networks, as they participate in biomass decomposition and facilitate plant nutrition through root symbioses. Nutrient content varies considerably between the main fungal habitats, such as soil, plant litter or decomposing dead wood, but there are also large differences within habitats. While some soils are heavily loaded with N, others are limited by N or P. One way in which nutrient availability can be reflected in fungi is their content in biomass. In this study, we determined the C, N, and P content (in dry mass) of fruiting bodies of 214 fungal species to inspect how phylogeny and membership in ecological guilds (soil saprotrophs, wood saprotrophs, and ectomycorrhizal fungi) affect the nutrient content of fungal biomass. The C content of fruiting bodies (415 ± 25 mg g-1) showed little variation (324-494 mg g-1), while the range of N (46 ± 20 mg g-1) and P (5.5 ± 3.0 mg g-1) contents was within one order of magnitude (8-103 mg g-1 and 1.0-18.9 mg g-1, respectively). Importantly, the N and P contents were significantly higher in the biomass of soil saprotrophic fungi compared to wood saprotrophic and ectomycorrhizal fungi. While the average C/N ratio in fungal biomass was 11.2, values exceeding 40 were recorded for some fungi living on dead wood, typically characterized by low N content. The N and P content of fungal mycelium also showed a significant phylogenetic signal, with differences in nutrient content being relatively low within species and genera of fungi. A strong correlation was found between N and P content in fungal biomass, while the correlation of N content and the N-containing fungal cell wall biopolymer-chitin showed only weak significance. The content of macronutrients in fungal biomass is influenced by the fungal life style and nutrient availability and is also limited by phylogeny.

Diverse Mycena Fungi and Their Potential for Gastrodia elata Germination.

It remains to be determined whether there is a geographical distribution pattern and phylogenetic signals for the Mycena strains with seed germination of the orchid plant Gastrodia elata. This study analyzed the community composition and phylogenetics of 72 Mycena strains associated with G. elata varieties (G. elata. f. glauca and G. elata. f. viridis) using multiple gene fragments (ITS+nLSU+SSU). We found that (1) these diverse Mycena phylogenetically belong to the Basidiospore amyloid group. (2) There is a phylogenetic signal of Mycena for germination of G. elata. Those strains phylogenetically close to M. abramsii, M. polygramma, and an unclassified Mycena had significantly higher germination rates than those to M. citrinomarginata. (3) The Mycena distribution depends on geographic site and G. elata variety. Both unclassified Mycena group 1 and the M. abramsii group were dominant for the two varieties of G. elata; in contrast, the M. citrinomarginata group was dominant in G. elata f. glauca but absent in G. elata f. viridis. Our results indicate that the community composition of numerous Mycena resources in the Zhaotong area varies by geographical location and G. elata variety. Importantly, our results also indicate that Mycena's phylogenetic status is correlated with its germination rate.

Shaft structure of the first metatarsal contains a strong phylogenetic signal in apes and humans.

OBJECTIVES: Metatarsal bones constitute a key functional unit of the foot in primates. While the form-function relationships of metatarsals have been extensively studied, particularly in relation to the loss of the grasping ability of the foot in humans in contrast to apes, the effect of phyletic history on the metatarsal morphology and its variability remains largely unknown. MATERIALS AND METHODS: Here, we evaluate how the strength of the phylogenetic signal varies from the first to the fifth metatarsal in humans, chimpanzees, gorillas, orangutans, gibbons, and Japanese macaques. We use computed tomography imaging and morphometric mapping to quantify the second moment of area around and along the metatarsal shaft and evaluate the strength of the phylogenetic signal with multivariate K-statistics. RESULTS: The shaft structure of the first metatarsal, but not the others, correlates well with the phylogeny of apes and humans. DISCUSSION: Given the importance of the first metatarsal for grasping and bipedal/quadrupedal locomotion, the strong phylogenetic but weak functional signal in its structure is unexpected. These findings suggest that the evolutionary diversification of hominoid locomotor behaviors, including human bipedality, is only partly reflected in form-function relationships of key skeletal elements, and that phylogenetic history acted as a major evolutionary constraint.

Local climate change velocities and evolutionary history explain multidirectional range shifts in a North American butterfly assemblage.

Species are often expected to shift their distributions either poleward or upslope to evade warming climates and colonise new suitable climatic niches. However, from 18-years of fixed transect monitoring data on 88 species of butterfly in the midwestern United States, we show that butterflies are shifting their centroids in all directions, except towards regions that are warming the fastest (southeast). Butterflies shifted their centroids at a mean rate of 4.87 km year-1. The rate of centroid shift was significantly associated with local climate change velocity (temperature by precipitation interaction), but not with mean climate change velocity throughout the species' ranges. Species tended to shift their centroids at a faster rate towards regions that are warming at slower velocities but increasing in precipitation velocity. Surprisingly, species' thermal niche breadth (range of climates butterflies experience throughout their distribution) and wingspan (often used as metric for dispersal capability) were not correlated with the rate at which species shifted their ranges. We observed high phylogenetic signal in the direction species shifted their centroids. However, we found no phylogenetic signal in the rate species shifted their centroids, suggesting less conserved processes determine the rate of range shift than the direction species shift their ranges. This research shows important signatures of multidirectional range shifts (latitudinal and longitudinal) and uniquely shows that local climate change velocities are more important in driving range shifts than the mean climate change velocity throughout a species' entire range.

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.

Climate of origin shapes variations in wood anatomical properties of 17 Picea species.

BACKGROUND: Variations in hydraulic conductivity may arise from species-specific differences in the anatomical structure and function of the xylem, reflecting a spectrum of plant strategies along a slow-fast resource economy continuum. Spruce (Picea spp.), a widely distributed and highly adaptable tree species, is crucial in preventing soil erosion and enabling climate regulation. However, a comprehensive understanding of the variability in anatomical traits of stems and their underlying drivers in the Picea genus is currently lacking especially in a common garden. RESULTS: We assessed 19 stem economic properties and hydraulic characteristics of 17 Picea species grown in a common garden in Tianshui, Gansu Province, China. Significant interspecific differences in growth and anatomical characteristics were observed among the species. Specifically, xylem hydraulic conductivity (Ks) and hydraulic diameter exhibited a significant negative correlation with the thickness to span ratio (TSR), cell wall ratio, and tracheid density and a significant positive correlation with fiber length, and size of the radial tracheid. PCA revealed that the first two axes accounted for 64.40% of the variance, with PC1 reflecting the trade-off between hydraulic efficiency and mechanical support and PC2 representing the trade-off between high embolism resistance and strong pit flexibility. Regression analysis and structural equation modelling further confirmed that tracheid size positively influenced Ks, whereas the traits DWT, D_r, and TSR have influenced Ks indirectly. All traits failed to show significant phylogenetic associations. Pearson's correlation analysis demonstrated strong correlations between most traits and longitude, with the notable influence of the mean temperature during the driest quarter, annual precipitation, precipitation during the wettest quarter, and aridity index. CONCLUSIONS: Our results showed that xylem anatomical traits demonstrated considerable variability across phylogenies, consistent with the pattern of parallel sympatric radiation evolution and global diversity in spruce. By integrating the anatomical structure of the stem xylem as well as environmental factors of origin and evolutionary relationships, our findings provide novel insights into the ecological adaptations of the Picea genus.

An evolutionary optimum amid moderate heritability in prokaryotic cell size.

We investigate the distribution and evolution of prokaryotic cell size based on a compilation of 5,380 species. Size spans four orders of magnitude, from 100 nm (Mycoplasma) to more than 1 cm (Thiomargarita); however, most species congregate heavily around the mean. The distribution approximates but is distinct from log normality. Comparative phylogenetics suggests that size is heritable, yet the phylogenetic signal is moderate, and the degree of heritability is independent of taxonomic scale (i.e., fractal). Evolutionary modeling indicates the presence of an optimal cell size to which most species gravitate. The size is equivalent to a coccus of 0.70 µm in diameter. Analyses of 1,361 species with sequenced genomes show that genomic traits contribute to size evolution moderately and synergistically. Given our results, scaling theory, and empirical evidence, we discuss potential drivers that may expand or shrink cells around the optimum and propose a stability landscape model for prokaryotic cell size.

Evolutionary trends of reproductive phenotype in Cycadales: an analysis of morphological evolution in Ceratozamia.

BACKGROUND AND AIMS: The size and shape of reproductive structures is especially relevant in evolution because these characters are directly related to the capacity of pollination and seed dispersal, a process that plays a basic role in evolutionary patterns. The evolutionary trajectories of reproductive phenotypes in gymnosperms have received special attention in terms of pollination and innovations related to the emergence of the Spermatophytes. However, variability of reproductive structures, evolutionary trends and the role of environment in the evolution of cycad species have not been well documented and explored. This study considered this topic under an explicitly phylogenetic and evolutionary approach that included a broad sampling of reproductive structures in the genus Ceratozamia. METHODS: We sampled 1400 individuals of 36 Ceratozamia species to explore evolutionary pattern and identify and evaluate factors that potentially drove their evolution. We analyzed characters for both pollen and ovulate strobili within a phylogenetic framework using different methods and characters (i. e., molecular and both quantitative and qualitative morphological) to infer phylogenetic relationships. Using this phylogenetic framework, evolutionary models of trait evolution for strobilar size were evaluated. In addition, quantitative morphological variation and its relation to environmental variables across species were analyzed. KEY RESULTS: We found contrasting phylogenetic signals between characters of pollen and ovulate strobili. These structures exhibited high morphological disparity in several characters related to size. Results of analyses of evolutionary trajectories suggested a stabilizing selection model. In regards to phenotype-environment, the analysis produced mixed results and differences for groups in the vegetation type where the species occur; however, a positive relationship with climatic variables was found. CONCLUSIONS: The integrated approach synthesized reproductive phenotypic variation with current phylogenetic hypotheses and provided explicit statements of character evolution. The characters of volume for ovulate strobili were the most informative, which could provide a reference for further study of the evolutionary complexity in Ceratozamia. Finally, heterogeneous environments, which are under changing weather conditions, promote variability of reproductive structures.

Evolution of wing shape in geometrid moths: phylogenetic effects dominate over ecology.

Locomotory performance is an important determinant of fitness in most animals, including flying insects. Strong selective pressures on wing morphology are therefore expected. Previous studies on wing shape in Lepidoptera have found some support for hypotheses relating wing shape to environment-specific selective pressures on aerodynamic performance. Here, we present a phylogenetic comparative study on wing shape in the lepidopteran family Geometridae, covering 374 species of the northern European fauna. We focused on 11 wing traits including aspect ratio, wing roundness, and the pointedness of the apex, as well as the ratio of forewing and hindwing areas. All measures were taken from images available on the internet, using a combination of tools available in Fiji software and R. We found that wing shape demonstrates a phylogenetically conservative pattern of evolution in Geometridae, showing similar or stronger phylogenetic signal than many of its potential predictors. Several wing traits showed statistically significant associations with predictors such as body size, phenology, and preference for forest habitats. Overall, however, all of these associations remained notably weak, with no wing shape being excluded for any value of the predictors, including body size. We conclude that, in geometrids, wing traits do not readily respond to selective pressures optimizing aerodynamic performance of the moths in different environments. Selection on wing shape may nevertheless operate through other functions of the wings, with the effectiveness of crypsis at rest being a promising candidate for further studies.

Geometric morphometrics casts light on phylogenetic relevance of cephalopod beak morphological.

The feeding organ of cephalopod species, the beak, can be used to reveal important ecological information. In this study, geometric morphometric approaches were employed to investigate the phylogenetic relevance and classification effect of beak lateral profile shape. The two-dimensional beak morphologies of 1164 pairs of 24 species from 13 genera and five families were constructed, and their evolutionary relationships and taxonomic status were confirmed using geometric morphometrics and molecular biology approaches. We also assessed the phylogenetic signals of beak shape. The analysis results show shape variation in the beak mainly in the rostrum, hood, and lateral wall. The overall shape parameters (all PCs) of the upper and lower beak are more useful for species identification. The shapes of the upper and lower beak show a strong phylogenetic signal, and the phenogram based on the beak shape basically reflected the families' taxonomic positions. We also hypothesized that the shape variation in the beaks of cephalopods may be ascribed to genetic and environmental differences. In summary, beaks are a reliable material for the classification of cephalopod species. Geometric morphometric approaches are a powerful tool to reveal the identification, phylogenetic relevance and phenotypic diversity of beak shape in cephalopods.

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.

New horizons for comparative studies and meta-analyses.

Comparative analyses and meta-analyses are key tools to elucidate broad biological principles, yet the two approaches often appear different in purpose. We propose an integrated approach that can generate deeper insights into ecoevolutionary processes. Marrying comparative and meta-analytic approaches will allow for (i) a more accurate investigation of drivers of biological variation, (ii) a greater ability to account for sources of non-independence in experimental data, (iii) more effective control of publication bias, and (iv) improved transparency and reproducibility. Stronger integration of meta-analytic and comparative studies can also broaden the scope from species-centric investigations to community-level responses and function-valued traits (e.g., reaction norms). We illuminate commonalities, differences, and the transformative potential of combining these methodologies for advancing ecology and evolutionary biology.

Testing the Tropical Niche Conservatism Hypothesis: Climatic Niche Evolution of Escallonia Mutis ex L. F. (Escalloniaceae).

We assess the Tropical Niche Conservatism Hypothesis in the genus Escallonia in South America using phylogeny, paleoclimate estimation and current niche modelling. We tested four predictions: (1) the climatic condition where the ancestor of Escallonia grew is megathermal; (2) the temperate niche is a derived condition from tropical clades; (3) the most closely related species have a similar current climate niche (conservation of the phylogenetic niche); and (4) there is a range expansion from the northern Andes to high latitudes during warm times. Our phylogenetic hypothesis shows that Escallonia originated 52.17 ± 0.85 My, in the early Eocene, with an annual mean temperature of 13.8 °C and annual precipitation of 1081 mm, corresponding to a microthermal to mesothermal climate; the species of the northern and central tropical Andes would be the ancestral ones, and the temperate species evolved between 32 and 20 My in a microthermal climate. The predominant evolutionary models were Brownian and Ornstein-Uhlenbeck. There was phylogenetic signal in 7 of the 9 variables, indicating conservation of the climatic niche. Escallonia would have originated in the central and southern Andes and reached the other environments by dispersion.

Effects of UV-B radiation on pollen germination and tube growth: A global meta-analysis.

Despite widespread recognition of pollen's potential sensitivity to ultraviolet-B (UV-B) radiation (280-315 nm), there remains ongoing debate surrounding the extent and mechanisms of this effect. In this study, using published data on pollen germination and tube growth including 377 pair-wise comparisons from 77 species in 30 families, we present the first global quantification of the effects of UV-B radiation on pollen germination and tube growth, along with its underlying mechanisms. Our results showed a substantial reduction in both pollen germination and tube growth in response to UV-B radiation, affecting 90.9 % and 84.2 % of species, respectively. Notably, these reductions exhibited phylogenetic constraints, highlighting the role of evolutionary history in shaping the sensitivity of pollen germination and tube growth to UV-B radiation. A negative correlation between elevation and the sensitivity of pollen tube growth was detected, suggesting that pollens from plants at higher elevations exhibit greater resistance to UV-B radiation. Our investigation also revealed that the effects of UV-B radiation on pollen germination and tube growth were influenced by a range of abiotic and biotic factors. Nevertheless, the intensity and duration of UV-B radiation exposure exhibited the highest explanatory power for the effects on both pollen germination and tube growth. This suggests that the responses of pollens to UV-B radiation are profoundly influenced by its dose, a critical consideration within the context of global change. In conclusion, our study provides valuable insights into the diverse responses of pollen germination and tube growth to UV-B radiation, highlighting the environment and species-dependent nature of pollen's susceptibility to UV-B radiation, with substantial implications for our understanding of the ecological and agricultural consequences of ongoing changes in UV-B radiation.

Phylogenetic signal in gut microbial community rather than in rodent metabolic traits.

Host phylogeny and environment have all been implicated in shaping the gut microbiota and host metabolic traits of mammals. However, few studies have evaluated phylogeny-associated microbial assembly and host metabolic plasticity concurrently, and their relationships on both short-term and evolutionary timescales. We report that the branching order of a gut microbial dendrogram was nearly congruent with phylogenetic relationships of seven rodent species, and this pattern of phylosymbiosis was intact after diverse laboratory manipulations. Laboratory rearing, diet or air temperature (Ta) acclimation induced alterations in gut microbial communities, but could not override host phylogeny in shaping microbial community assembly. A simulative heatwave reduced core microbiota diversity by 26% in these species, and led to an unmatched relationship between the microbiota and host metabolic phenotypes in desert species. Moreover, the similarity of metabolic traits across species at different Tas was not correlated with phylogenetic distance. These data demonstrated that the gut microbial assembly showed strong concordance with host phylogeny and may be shaped by environmental variables, whereas host metabolic traits did not seem to be linked with phylogeny.

Low interspecific variation and no phylogenetic signal in additive genetic variance in wild bird and mammal populations.

Evolutionary adaptation through genetic change requires genetic variation and is a key mechanism enabling species to persist in changing environments. Although a substantial body of work has focused on understanding how and why additive genetic variance (V A) differs among traits within species, we still know little about how they vary among species. Here we make a first attempt at testing for interspecific variation in two complementary measures of V A and the role of phylogeny in shaping this variation. To this end, we performed a phylogenetic comparative analysis using 1822 narrow-sense heritability (h 2) for 68 species of birds and mammals and 378 coefficients of additive genetic variance (CV A) estimates for 23 species. Controlling for within-species variation attributable to estimation method and trait type, we found some interspecific variation in h 2 (~15%) but not CV A. Although suggestive of interspecific variation in the importance of non-(additive) genetic sources of variance, sample sizes were insufficient to test this hypothesis directly. Additionally, although power was low, no phylogenetic signal was detected for either measure. Hence, while this suggests interspecific variation in V A is probably small, our understanding of interspecific variation in the adaptive potential of wild vertebrate populations is currently hampered by data limitations, a scarcity of CV A estimates and a measure of their uncertainty in particular.

Claw shape variation in oribatid mites of the genera Carabodes and Caleremaeus: exploring the interplay of habitat, ecology and phylogenetics.

Background: Claws are a commonly observed biological adaptation across a wide range of animal groups. They serve different functions and their link to evolution is challenging to analyze. While there are many studies on the comparative anatomy and morphology of claws in reptiles, birds and several arthropods, knowledge about claws of soil-living oribatid mites, is still limited. Recent research on intertidal oribatid mites has shown that claw shape is strongly correlated with microhabitat and is subject to ecological selective pressures. However, the selective constraints shaping claws in terrestrial oribatid mites are still unknown. Methods: In this study, 300 specimens from 12 different species and two genera were examined. Geometric morphometrics were used to quantify claw length and curvature, and to analyze two-dimensional claw shape. In combination with molecular phylogenetic analyses of investigated populations phylogenetic signal was quantified within genera using Blomberg's K and random replicates. Additionally, ecological information on the investigated species was gathered from previous studies and compiled into tables. Results: The claw shapes of Carabodes species vary moderately, with the three species C. reticulatus, C. rugosior and C. tenuis deviating the most from the others. These three species are only found in a small number of habitats, which may require a more specialized claw shape. Our results show that there is a phylogenetic influence on claw shape in Carabodes but not in Caleremaeus. Additionally, habitat specificity and lifestyle were found to have ecological impact on claw shape in both genera. The present results demonstrate that characteristics of the claws of terrestrial oribatid mites are correlated with ecology, but this correlation is apparently weaker than in intertidal oribatid mites that are prone to strong external forces.

Evolution of the Mutation Spectrum Across a Mammalian Phylogeny.

Although evolutionary biologists have long theorized that variation in DNA repair efficacy might explain some of the diversity of lifespan and cancer incidence across species, we have little data on the variability of normal germline mutagenesis outside of humans. Here, we shed light on the spectrum and etiology of mutagenesis across mammals by quantifying mutational sequence context biases using polymorphism data from thirteen species of mice, apes, bears, wolves, and cetaceans. After normalizing the mutation spectrum for reference genome accessibility and k-mer content, we use the Mantel test to deduce that mutation spectrum divergence is highly correlated with genetic divergence between species, whereas life history traits like reproductive age are weaker predictors of mutation spectrum divergence. Potential bioinformatic confounders are only weakly related to a small set of mutation spectrum features. We find that clock-like mutational signatures previously inferred from human cancers cannot explain the phylogenetic signal exhibited by the mammalian mutation spectrum, despite the ability of these signatures to fit each species' 3-mer spectrum with high cosine similarity. In contrast, parental aging signatures inferred from human de novo mutation data appear to explain much of the 1-mer spectrum's phylogenetic signal in combination with a novel mutational signature. We posit that future models purporting to explain the etiology of mammalian mutagenesis need to capture the fact that more closely related species have more similar mutation spectra; a model that fits each marginal spectrum with high cosine similarity is not guaranteed to capture this hierarchy of mutation spectrum variation among species.