Phylogenetic Regression Research Articles

Cutting the sap: First molecular phylogeny of twig‐girdler longhorn beetles (Coleoptera: Cerambycidae: Lamiinae: Onciderini) suggests shifts in host plant attack behaviours contributed to morphological evolution

AbstractUnderstanding how novel adaptive traits arise, evolve and impact other aspects of an organism's phenotype is a foundational question in evolutionary biology. We explore this by focusing on Onciderini (Coleoptera: Cerambycidae: Lamiinae), a tribe of longhorn beetles commonly referred to as twig girdlers because the females of some species girdle live trees to deposit their eggs. We reconstructed the first time‐calibrated phylogeny of the Onciderini, based on three genetic markers (cytochrome c oxidase subunit 1 ‘cox1’, Wingless ‘Wg’ and carbamoyl‐phosphate synthase domain of the CAD locus ‘CPS’), and used morphometric data, employing multiple models of trait evolution and phylogenetic regressions, to assess whether girdling behaviour predicts head size or head size sexual dimorphism. Our results indicate that onciderines originated c. 49 million years ago ‘Ma’ (95% highest posterior densities ‘HPD’: 44.1–54.57 Ma) and consist of two major lineages, which we formalize as subtribes Hypsiomina and Onciderina. Additionally, our analyses revealed several taxonomic inconsistencies within the tribe, which we rectify by proposing new synonymies, including Jamesia Jekel as a synonym of Hypselomus Perty, Cipriscola Dillon & Dillon of Hypsioma Audinet‐Serville and Psyllotoxus Thomson and Taricanus Thomson as a synonym of Oncideres Lacordaire, and a new combination for Periergates kenjii Nearns & Swift as Oncideres kenjii (Nearns & Swift) comb. nov. Using this new phylogenetic framework, we identified five independent unidirectional origins of girdling behaviour within the tribe. Overall, our results suggest that girdling influences head size evolution in Onciderini, as it is associated with an increase in both male and female head size of girdler species. Despite this, girdling and non‐girdling lineages do not consistently differ in head size sexual dimorphism. This study refines the classification of Onciderini and marks a significant step in understanding the evolutionary dynamics shaping the diversity of twig‐girdler beetles, a group with notable ecological and economic importance.

Evolutionary origins, macroevolutionary dynamics, and climatic niche space of the succulent plant syndrome in the Caryophyllales.

The succulent plant syndrome is defined by the coordination of traits that enhance internal water storage within plant tissues. Although distributed globally in different habitats, succulent plants are thought to have evolved to avoid drought in arid regions, due to trait modifications that decrease tissue water deficits. We evaluated the evolution and the ecological significance of the succulent strategy at a global scale by comparing the climatic niche of species displaying succulence within the Core Caryophyllales with their non-succulent relatives. We assembled and curated a worldwide dataset of 201K georeferenced records belonging to 5447 species within 28 families, and analyzed the climatic niche of species along with their origin and evolutionary trajectories using ecological niche modeling, phylogenetic regression, divergence dates and ancestral state estimation. Results indicate the Core Caryophyllales have inhabited drylands since their origin in the Early Cretaceous. However, the succulent syndrome appeared and diversified during later geologic periods. The climatic niche space of succulents is narrower than non-succulent relatives, but no climate niche separation was detected between groups. Our results support alternative interpretations on the environmental and ecological forces that spurred the origin and diversification of the succulent plant syndrome and the radiation of rich succulent lineages.

Unifying approaches from statistical genetics and phylogenetics for mapping phenotypes in structured populations.

In both statistical genetics and phylogenetics, a major goal is to identify correlations between genetic loci or other aspects of the phenotype or environment and a focal trait. In these 2 fields, there are sophisticated but disparate statistical traditions aimed at these tasks. The disconnect between their respective approaches is becoming untenable as questions in medicine, conservation biology, and evolutionary biology increasingly rely on integrating data from within and among species, and once-clear conceptual divisions are becoming increasingly blurred. To help bridge this divide, we lay out a general model describing the covariance between the genetic contributions to the quantitative phenotypes of different individuals. Taking this approach shows that standard models in both statistical genetics (e.g., genome-wide association studies; GWAS) and phylogenetic comparative biology (e.g., phylogenetic regression) can be interpreted as special cases of this more general quantitative-genetic model. The fact that these models share the same core architecture means that we can build a unified understanding of the strengths and limitations of different methods for controlling for genetic structure when testing for associations. We develop intuition for why and when spurious correlations may occur analytically and conduct population-genetic and phylogenetic simulations of quantitative traits. The structural similarity of problems in statistical genetics and phylogenetics enables us to take methodological advances from one field and apply them in the other. We demonstrate by showing how a standard GWAS technique-including both the genetic relatedness matrix (GRM) as well as its leading eigenvectors, corresponding to the principal components of the genotype matrix, in a regression model-can mitigate spurious correlations in phylogenetic analyses. As a case study, we re-examine an analysis testing for coevolution of expression levels between genes across a fungal phylogeny and show that including eigenvectors of the covariance matrix as covariates decreases the false positive rate while simultaneously increasing the true positive rate. More generally, this work provides a foundation for more integrative approaches for understanding the genetic architecture of phenotypes and how evolutionary processes shape it.

Drivers of interlineage variability in mitogenomic evolutionary rates in Platyhelminthes.

Studies of forces driving interlineage variability in the evolutionary rates (both sequence and architecture) of mitochondrial genomes often produce contradictory results. Flatworms (Platyhelminthes) exhibit the fastest-evolving mitogenomic sequences among all bilaterian phyla. To test the effects of multiple factors previously associated with different aspects of mitogenomic evolution, we used mitogenomes of 223 flatworm species, phylogenetic multilevel regression models, and causal inference. Thermic host environment (endothermic vs. ectothermic) had nonsignificant impacts on both sequence evolution and mitogenomic size. Mitogenomic gene order rearrangements (GORR) were mostly positively correlated with mitogenomic size (R2 ≈ 20-30%). Longevity was not (negatively) correlated with sequence evolution in flatworms. The predominantly free-living "turbellaria" exhibited much shorter branches and faster-evolving mitogenomic architecture than parasitic Neodermata. As a result, "parasitism" had a strong explanatory power on the branch length variability (>90%), and there was a negative correlation between GORR and branch length. However, the stem branch of Neodermata comprised 63.6% of the total average branch length. This evolutionary period was also marked by a high rate of gene order rearrangements in the ancestral Neodermata. We discuss how this period of rapid evolution deep in the evolutionary history may have decoupled sequence evolution rates from longevity and GORR, and overestimated the explanatory power of "parasitism". This study shows that impacts of variables often vary across lineages, and stresses the importance accounting for the episodic nature of evolutionary patterns in studies of mitogenomic evolution.

Centromere drive may propel the evolution of chromosome and genome size in plants.

Genome size is influenced by natural selection and genetic drift acting on variations from polyploidy and repetitive DNA sequences. We hypothesized that centromere drive, where centromeres compete for inclusion in the functional gamete during meiosis, may also affect genome and chromosome size. This competition occurs in asymmetric meiosis, where only one of the four meiotic products becomes a gamete. If centromere drive influences chromosome size evolution, it may also impact post-polyploid diploidization, where a polyploid genome is restructured to function more like a diploid through chromosomal rearrangements, including fusions. We tested if plant lineages with asymmetric meiosis exhibit faster chromosome size evolution compared to those with only symmetric meiosis, which lack centromere drive as all four meiotic products become gametes. We also examined if positive selection on centromeric histone H3 (CENH3), a protein that can suppress centromere drive, is more frequent in these asymmetric lineages. We analyzed plant groups with different meiotic modes: asymmetric in gymnosperms and angiosperms, and symmetric in bryophytes, lycophytes, and ferns. We selected species based on available CENH3 gene sequences and chromosome size data. Using Ornstein-Uhlenbeck evolutionary models and phylogenetic regressions, we assessed the rates of chromosome size evolution and the frequency of positive selection on CENH3 in these clades. Our analyses showed that clades with asymmetric meiosis have a higher frequency of positive selection on CENH3 and increased rates of chromosome size evolution compared to symmetric clades. Our findings support the hypothesis that centromere drive accelerates chromosome and genome size evolution, potenatially also influencing the process of post-polyploid diploidization. We propose a model which in a single famework helps explain the stability of chromosome size in symmetric lineages (bryophytes, lycophytes, and ferns) and its variability in asymmetric lineages (gymnosperms and angiosperms), providing a foundation for future research in plant genome evolution.

Non‐linear phylogenetic regression using regularised kernels

Abstract Phylogenetic regression is a type of generalised least squares (GLS) method that incorporates a modelled covariance matrix based on the evolutionary relationships between species (i.e. phylogenetic relationships). While this method has found widespread use in hypothesis testing via phylogenetic comparative methods, such as phylogenetic ANOVA, its ability to account for non‐linear relationships has received little attention. To address this, here we implement a phylogenetic Kernel Ridge Regression (phyloKRR) method that utilises GLS in a high‐dimensional feature space, employing linear combinations of phylogenetically weighted data to account for non‐linearity. We analysed two biological datasets using the Radial Basis Function and linear kernel function. The first dataset contained morphometric data, while the second dataset comprised discrete trait data and diversification rates as response variable. Hyperparameter tuning of the model was achieved through cross‐validation rounds in the training set. In the tested biological datasets, phyloKRR reduced the error rate (as measured by RMSE) by around 20% compared to linear‐based regression when data did not exhibit linear relationships. In simulated datasets, the error rate decreased almost exponentially with the level of non‐linearity. These results show that introducing kernels into phylogenetic regression analysis presents a novel and promising tool for complementing phylogenetic comparative methods. We have integrated this method into Python package named phyloKRR, which is freely available at: https://github.com/ulises‐rosas/phylokrr.

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

Diversity in pappus shapes and size in Asteraceae suggests an adaptive response to dispersion challenges adjusting diaspores to optimal phenotypic configurations. Here, by analysing the relationship among pappus-cypsela size relationships, flight performance and pappus types in an evolutionary context, we evaluate the role of natural selection acting on the evolution of diaspore configuration at a macro-ecological scale in the daisy family. To link pappus-cypsela size relationships with flight performance we collected published data on these traits from 82 species. This allowed us to translate morphometric traits in flight performance for 150 species represented in a fully resolved backbone phylogeny of the daisy family. Through ancestral reconstructions and evolutionary model selection, we assessed whether flight performance was associated with and constrained by different pappus types. Additionally, we evaluated, through phylogenetic regressions, whether species with different pappus types exhibited evolutionary allometric pappus-cypsela size relationships. The setose pappus type had the highest flight performances and represented the most probable ancestral state in the family. Stepwise changes in pappus types independently led from setose to multiple instances of pappus loss with associated reduction in flight performance. Flight performance evolution was best modelled as constrained by five adaptive regimes represented by specific pappus types which correspond to specific optimal diaspore configurations that are distinct in pappus-cypsela allometric relationships. Evolutionary modelling suggests natural selection as the main factor of diaspore configuration changes which proceeded towards five optima, often overcoming constraints imposed by allometric relationships and favouring evolution in certain directions. With the perspective that natural selection is the main process driving the observed patterns, various biotic and abiotic are suggested as principal drivers of transitions in diaspore configurations along space and time in the daisy family history. The results also allow discussion of evolutionary changes in a historical context.

Reduced palatability, fast flight, and tails: decoding the defence arsenal of Eudaminae skipper butterflies in a Neotropical locality.

Prey often rely on multiple defences against predators, such as flight speed, attack deflection from vital body parts, or unpleasant taste, but our understanding on how often and why they are co-exhibited remains limited. Eudaminae skipper butterflies use fast flight and mechanical defences (hindwing tails), but whether they use other defences like unpalatability (consumption deterrence) and how these defences interact have not been assessed. We tested the palatability of 12 abundant Eudaminae species in Peru, using training and feeding experiments with domestic chicks. Further, we approximated the difficulty of capture based on flight speed and quantified it by wing loading. We performed phylogenetic regressions to find any association between multiple defences, body size, and habitat preference. We found a broad range of palatability in Eudaminae, within and among species. Contrary to current understanding, palatability was negatively correlated with wing loading, suggesting that faster butterflies tend to have lower palatability. The relative length of hindwing tails did not explain the level of butterfly palatability, showing that attack deflection and consumption deterrence are not mutually exclusive. Habitat preference (open or forested environments) did not explain the level of palatability either, although butterflies with high wing loading tended to occupy semi-closed or closed habitats. Finally, the level of unpalatability in Eudaminae is size dependent. Larger butterflies are less palatable, perhaps because of higher detectability/preference by predators. Altogether, our findings shed light on the contexts favouring the prevalence of single versus multiple defensive strategies in prey.

The smallest angiosperm genomes may be the price for effective traps of bladderworts.

Species of the carnivorous family Lentibulariaceae exhibit the smallest genomes in flowering plants. We explored the hypothesis that their minute genomes result from the unique mitochondrial cytochrome c oxidase (COX) mutation. The mutation may boost mitochondrial efficiency, which is especially useful for suction-bladder traps of Utricularia, but also increase DNA-damaging reactive oxygen species, leading to genome shrinkage through deletion-biased DNA repair. We aimed to explore this mutation's impact on genome size, providing insights into genetic mutation roles in plant genome evolution under environmental pressures. We compiled and measured genome and mean chromosome sizes for 127 and 67 species, respectively, representing all three genera (Genlisea, Pinguicula, and Utricularia) of Lentibulariaceae. We also isolated and analyzed COX sequences to detect the mutation. Through phylogenetic regressions and Ornstein-Uhlenbeck models of trait evolution, we assessed the impact of the COX mutation on the genome and chromosome sizes across the family. Our findings reveal significant correlations between the COX mutations and smaller genome and chromosome sizes. Specifically, species carrying the ancestral COX sequence exhibited larger genomes and chromosomes than those with the mutation. This evidence supports the notion that the COX mutation contributes to genome downsizing, with statistical analyses confirming a directional evolution towards smaller genomes in species harboring these mutations. Our study confirms that the COX mutation in Lentibulariaceae is associated with genome downsizing, likely driven by increased reactive oxygen species production and subsequent DNA damage requiring deletion-biased repair mechanisms. While boosting mitochondrial energy output, this genetic mutation compromises genome integrity and may potentially affect recombination rates, illustrating a complex trade-off between evolutionary advantages and disadvantages. Our results highlight the intricate processes by which genetic mutations and environmental pressures shape genome size evolution in carnivorous plants.

Disentangling ecological drivers of interspecific achromatic plumage variation in birds

AbstractAimUnderstanding the ecological determinants of interspecific achromatic (light‐to‐dark) plumage variation in birds is crucial yet challenging due to the complex interplay of climatic, habitat‐related, and morphological influences. This study aimed to disentangle the effects of temperature, precipitation, habitat openness, body mass and hand‐wing index (HWI, a widely used single‐parameter proxy for the extent to which a species relies on flight) on shaping achromatic plumage variation among bird species.LocationGlobal.Time PeriodContemporary.Major Taxa StudiedBirds.MethodsBased on data from over 8000 sessile bird species globally, we employed phylogenetic linear regressions to account for achromatic plumage colour in relation to temperature, precipitation, habitat openness, body mass and HWI, while correcting for phylogenetic non‐independence between species. Furthermore, we conducted phylogenetic path analyses to decompose direct from indirect effects.ResultsWe found that temperature, precipitation, habitat openness and body mass exerted separate but interactive effects on the variation in achromatic colour across species. Species inhabiting cold, wet or densely vegetated environments were darker coloured, while smaller species were lighter. Darker plumage was more strongly related to higher precipitation in colder regions for nocturnal species. For diurnal species, darker plumage was more closely associated with higher precipitation in more open habitats, whereas lighter plumage was more linked to lower mass in denser habitats. Noteworthy was the identification of a substantial correlation between achromatic colour and HWI. Diurnal species that are more aerial were lighter. Conversely, nocturnal flyers, especially females, tended to be darker.Main ConclusionsThe findings highlight the multifaceted nature of plumage coloration evolution, with adaptations for thermal efficiency, crypsis, signalling, waterproofing or protection against bacteria. However, the variable relative importance of these factors among groups emphasizes the significance of each factor in different contexts.

The Relationship between Grinnellian and Eltonian Niche Characteristics and Passerine Distribution across a Latitudinal Gradient

The degree to which Grinnellian and Eltonian niche characteristics influence species distribution may depend on latitude. Tropical regions are environmentally stable and resource-rich, whereas temperate regions are comparatively less environmentally stable (e.g., environmental filtering). Moreover, phylogenetic niche conservatism could influence distributions by inhibiting the ability for species to colonize environmentally different locations. Herein, we examine relationships between niche characteristics, passerine distributions, and phylogenetic niche conservatism across the latitudinal gradient. We used environmental and climatic variables to characterize Grinnellian niches and diets to characterize Eltonian niches. We conducted variation partitioning with retained components from ordination methods to evaluate the degree of association of Grinnellian and Eltonian niche characteristics with passerine distribution across latitudes. We examined the relationship between phylogenetic signal and niche characteristics with a phylogenetic regression. Passerine distributions were more related to environmental gradients than resources across latitudes. While niche conservatism was prevalent in Eltonian niche characteristics, phylogeny was related to Grinnellian niche characteristics in only 46% of biomes. There was no latitudinal gradient in phylogenetic niche conservatism or the degree to which Eltonian and Grinnellian niche characteristics relate to passerine distribution. Niche conservatism, albeit weak, was present for Grinnellian niche characteristics, thus potentially influencing the expansion of passerine distributions into the northern hemisphere.

Gene Protein Sequence Evolution Can Predict the Rapid Divergence of Ovariole Numbers in the Drosophila melanogaster Subgroup.

Ovaries play key roles in fitness and evolution: they are essential female reproductive structures that develop and house the eggs in sexually reproducing animals. In Drosophila, the mature ovary contains multiple tubular egg-producing structures known as ovarioles. Ovarioles arise from somatic cellular structures in the larval ovary called terminal filaments (TFs), formed by TF cells and subsequently enclosed by sheath (SH) cells. As in many other insects, ovariole number per female varies extensively in Drosophila. At present, however, there is a striking gap of information on genetic mechanisms and evolutionary forces that shape the well-documented rapid interspecies divergence of ovariole numbers. To address this gap, here we studied genes associated with Drosophila melanogaster ovariole number or functions based on recent experimental and transcriptional datasets from larval ovaries, including TFs and SH cells, and assessed their rates and patterns of molecular evolution in five closely related species of the melanogaster subgroup that exhibit species-specific differences in ovariole numbers. From comprehensive analyses of protein sequence evolution (dN/dS), branch-site positive selection, expression specificity (tau), and phylogenetic regressions (phylogenetic generalized least squares), we report evidence of 42 genes that showed signs of playing roles in the genetic basis of interspecies evolutionary change of Drosophila ovariole number. These included the signaling genes upd2 and Ilp5 and extracellular matrix genes vkg and Col4a1, whose dN/dS predicted ovariole numbers among species. Together, we propose a model whereby a set of ovariole-involved gene proteins have an enhanced evolvability, including adaptive evolution, facilitating rapid shifts in ovariole number among Drosophila species.

Evolutionary allometry of the canid baculum (Carnivora: Mammalia)

Abstract While the mammalian baculum shows enormous morphological variability, the baculum of canids is highly conserved, with most variation restricted to size. Here, we explore the allometric relationship between baculum length and body size in extant and extinct canids. Examination of 26 species in the extant subfamily Caninae using standard linear regression revealed isometry. Phylogenetic regression also revealed an allometric slope indistinguishable from isometry. This pattern differs from the substantially negative slopes seen in other mammalian clades. The strength of the canid allometric relationship (r2) is also greater than in other clades, suggesting functional constraints on their baculum size. The constraints may be related to the copulatory tie that is characteristic of canids, and/or their monogamous mating system. Complete bacula are known from just four extinct species. The two complete bacula from the extinct subfamily Borophaginae (Aelurodon ferox and Aelurodon stirtoni) fall on the same allometric relationship as the living canids. However, the baculum of the extinct dire wolf (Aenocyon dirus, from the extant subfamily Caninae) and from the extinct subfamily Herperocyoninae, Hesperocyon gregarius, are significantly longer than expected based on their body sizes, suggesting that they may have had a different reproductive biology from that of extant canines.

Range and climate niche shifts in European and North American breedingbirds.

Species respond dynamically to climate change and exhibit time lags. Consequently, species may not occupy their full climatic niche during range shifting. Here, we assessed climate niche tracking during recent range shifts of European and UnitedStates (US) birds. Using data from two European bird atlases and from the North American Breeding Bird Survey betweenthe 1980s and 2010s, we analysed range overlap and climate niche overlap based on kernel density estimation. Phylogenetic multiple regression was used to assess the effect of species morphological, ecological and biogeographic traits on range and niche metrics. European birds shifted their ranges north and north-eastwards, US birds westwards. Range unfilling was lower than expected by null models, and niche expansion was more common than niche unfilling. Also, climate niche tracking was generally lower in US birds and poorly explained by species traits. Overall, our results suggest that dispersal limitations were minor in range shifting birds in Europe and theUSA while delayed extinctions from unfavourable areas seem more important. Regional differences could be related to differences in land use history and monitoring schemes. Comparative analyses of range and niche shifts provide a useful screening approach for identifying the importance of transient dynamics and time-lagged responses to climate change. This article is part of the theme issue 'Ecological novelty and planetary stewardship: biodiversity dynamics in a transforming biosphere'.

Strong habitat and seasonal phenology effects on the evolution of self-compatibility, clonality and pollinator shifts in Lachenalia (Asparagaceae: Scilloideae).

Plants employ a diversity of reproductive safeguarding strategies to circumvent the challenge of pollen limitation. Focusing on southern African Lachenalia (Asparagaceae: Scilloideae), we test the hypothesis that the evolution of reproductive safeguarding traits (self-compatibility, autonomous selfing, bird pollination and clonal propagation) is favoured in species occupying conditions of low insect abundance imposed by critically infertile fynbos heathland vegetation and by flowering outside the austral spring insect abundance peak. We trace the evolution of these traits and selective regimes on a dated, multi-locus phylogeny of Lachenalia and assess their evolutionary associations using ordinary and phylogenetic regression. Ancestral state reconstructions identify an association with non-fynbos vegetation and spring flowering as ancestral in Lachenalia, the transition to fynbos vegetation and non-spring flowering taking place multiple times. They also show that self-compatibility, autofertility, bird pollination and production of multiple clonal offsets have evolved repeatedly. Regression models suggest that bird pollination and self-compatibility are selected for in fynbos and in non-spring flowering lineages, with autofertility being positively associated with non-spring flowering. These patterns support the interpretation of these traits as reproductive safeguarding adaptations under reduced insect pollinator abundance. We find no evidence to support the interpretation of clonal propagation as a reproductive safeguarding strategy.

Gymnosperms demonstrate patterns of fine‐root trait coordination consistent with the global root economics space

Abstract Gymnosperms encompass a diverse group of mostly woody plants with high ecological and economic value, yet little is known about the scope and organization of fine‐root trait diversity among gymnosperms due to the undersampling of most gymnosperm families and the dominance of angiosperm groups in recent syntheses. New and existing data were compiled for morphological traits (root diameter, length, tissue density, specific root length [SRL] and specific root area [SRA]), the architectural trait branching ratio, root nitrogen content [N] and mycorrhizal colonization. We used phylogenetic least squares regression and principal component analysis to determine trait–trait relationships and coordination across 66 species, representing 11 of the 12 extant gymnosperm families from boreal, temperate, subtropical and tropical biomes. Finally, we compared the relationship between family divergence time and mean trait values to determine whether evolutionary history structured variation in fine‐root traits within the gymnosperm phylogeny. Wide variation in gymnosperm root traits could be largely captured by two primary axes of variation defined by SRL and diameter, and root tissue density and root nitrogen, respectively. However, individual root length and SRA each had significant correlations with traits defining both main axes of variation. Neither mycorrhizal colonization nor root branching ratio were closely related to other traits. We did not observe a directional evolution of mean trait values from older to more recently diverged gymnosperm families. Synthesis. Despite their unique evolutionary history, gymnosperms display a root economic space similar to that identified in angiosperms, likely reflecting common constraints on plants adapting to diverse environments in both groups. These findings provide greater confidence that patterns observed in broad syntheses justly capture patterns of trait diversity among multiple, distinct lineages. Additionally, independence between root architecture and other traits may support greater diversity in below‐ground resource acquisition strategies. Unlike angiosperms, there were no clear trends towards increasingly thin roots over evolutionary time, possibly because of lower diversification rates or unique biogeographic history among gymnosperms, though additional observations are needed to more richly test evolutionary trends among gymnosperms.

Proteomic diversification of spermatostyles among six species of whirligig beetles.

Seminal fluid protein composition is complex and commonly assumed to be rapidly divergent due to functional interactions with both sperm and the female reproductive tract (FRT), both of which evolve rapidly. In addition to sperm, seminal fluid may contain structures, such as mating plugs and spermatophores. Here, we investigate the evolutionary diversification of a lesser-known ejaculate structure: the spermatostyle, which has independently arisen in several families of beetles and true bugs. We characterized the spermatostyle proteome, in addition to spermatostyle and FRT morphology, in six species of whirligig beetles (family Gyrinidae). Spermatostyles were enriched for proteolytic enzymes, and assays confirmed they possess proteolytic activity. Sperm-leucylaminopeptidases (S-LAPs) were particularly abundant, and their localization to spermatostyles was confirmed by immunohistochemistry. Although there was evidence for functional conservation of spermatostyle proteomes across species, phylogenetic regressions suggest evolutionary covariation between protein composition and the morphology of both spermatostyles and FRTs. We postulate that S-LAPs (and other proteases) have evolved a novel structural role in spermatostyles and discuss spermatostyles as adaptations for delivering male-derived materials to females.

The environmental conditions of endemism hotspots shape the functional traits of mammalian assemblages.

Endemic (small-ranged) species are distributed non-randomly across the globe. Regions of high topography and stable climates have higher endemism than flat, climatically unstable regions. However, it is unclear how these environmental conditions interact with and filter mammalian traits. Here, we characterize the functional traits of highly endemic mammalian assemblages in multiple ways, testing the hypothesis that these assemblages are trait-filtered (less functionally diverse) and dominated by species with traits associated with small range sizes. Compiling trait data for more than 5000 mammal species, we calculated assemblage means and multidimensional functional metrics to evaluate the distribution of traits across each assemblage. We then related these metrics to the endemism of global World Wildlife Fund ecoregions using linear models and phylogenetic fourth-corner regression. Highly endemic mammalian assemblages had small average body masses, low fecundity, short lifespans and specialized habitats. These traits relate to the stable climate and rough topography of endemism hotspots and to mammals' ability to expand their ranges, suggesting that the environmental conditions of endemism hotspots allowed their survival. Furthermore, species living in endemism hotspots clustered near the edges of their communities' functional spaces, indicating that abiotic trait filtering and biotic interactions act in tandem to shape these communities.

Niche position and niche breadth effects on population abundances: A case study of New World Warblers (Parulidae).

Species abundance patterns are influenced by a myriad of factors, including habitat availability and ecological niche characteristics. However, the evidence concerning the specific impact factors such as niche position and niche breadth on mean and maximum abundances in vertebrates at a broad geographical scale remains inconclusive. In this study, we investigated the influence of niche position and breadth on the abundance of 47 species of birds belonging to the Parulidae family, commonly known as New World Warblers. We obtained data on abundance and presence records spanning the reproductive distribution of these species and employed the outlying mean index analysis to calculate niche position and niche breadth. We assessed the relationship between abundance metrics and niche descriptors using phylogenetic regressions to account for the non-independence resulting from phylogenetic ancestry. Initially, we developed individual models for each predictor and subsequently formulated a multi-predictor model encompassing niche position, niche breadth, and their interaction. Our findings revealed a negative relationship between niche position and both mean and maximum abundance, while niche breadth exhibited a positive relationship with these niche characteristics. Notably, the results of the multi-predictor models indicated that niche position exerted the most substantial influence on both mean and maximum abundance. Additionally, the interaction between niche position and niche breadth had the most positive and significant contribution to mean population abundance. This study underscores the need for future research in other vertebrates to delve into the mechanisms underlying these patterns. Such endeavors will not only enhance our understanding of ecological dynamics but also equip us with predictive capabilities to anticipate population responses to environmental changes effectively.

Clonal alien plants in the mountains spread upward more extensively and faster than non-clonal

Alien species are colonizing mountain ecosystems and increasing their elevation ranges in response to ongoing climate change and anthropogenic disturbances, posing increasing threats to native species. However, how quickly alien species spread upward and what drives their invasion remains insufficiently understood. Here, using 26,952 occurrence records of 58 alien plant species collected over two centuries in the Czech Republic, we explored the elevation range and invasion speed of each alien species and the underlying factors driving these variables. We collected species traits relevant for invasion (e.g., clonality, flowering time, life span, invasion status, height, mycorrhizal type, native range, naturalized range, monoploid genome size, and Ellenberg-type indicator values for light, temperature, and nitrogen), human-associated factors (e.g., introduction pathways and the sum of economic use types), and minimum residence time. We explored the relationships between these factors and species’ elevation range and invasion speed using phylogenetic regressions. Our results showed that 58 alien species have been expanding upward along mountain elevations in the Czech Republic over the past two centuries. A stronger effect of species’ traits than human-associated factors has been revealed, e.g., clonality was a key trait supporting the invasion of alien species into the mountains, while human-associated factors showed no effect. Our findings highlight that the characteristics associated with rapid reproduction and spread are crucial for alien species’ invasion into montane regions. Identifying key drivers of this process is important for predicting the spatiotemporal dynamics of alien species in high-altitude ecosystems and thus employing apposite measures to reduce the threat to native plant species.