Measuring phylogenetic signal between categorical traits and phylogenies

Anther modes influence diversification rates in the animal-pollinated species-rich Didymocarpoideae

Stamen traits significantly influence pollen presentation patterns and play a critical role in mating systems, floral evolution, and the diversification of angiosperms. Anthers within a flower can either develop freely or be united, resulting in 3 primary modes: separated anthers, synandry, and paired-united anthers. The impact of these anther modes on species diversification remains inadequately understood. In this study, we employed 14 plastid and nuclear ribosomal markers from 789 species to reconstruct the phylogenetic relationships of the Old World Didymocarpoideae, which are predominantly bee-pollinated. We then investigated the evolutionary patterns of anther modes and assessed their potential role in species diversification. Our findings indicate that the evolution of anther modes and associated floral traits exhibited significant trait conservatism. Paired-united anthers likely represent the ancestral form of Didymocarpoideae, characterized by 4 fertile stamens with exserted locations for the anthers. Notably, derived synandry was associated with relatively high rates of species diversification, particularly in the species-rich Cyrtandra, Primulina, and Streptocarpus, which could be due to enhanced pollination precision facilitated by aggregations of anthers and pollen grains. This study elucidates the evolutionary transitions of different anther modes while highlighting their influence on diversification rates within Didymocarpoideae.

The Increased Environmental Niche of Dual-Mycorrhizal Woody Species

The presence and distribution of mycorrhizal symbionts can influence plant distribution through specific host-mycorrhiza symbiosis interactions. However, generalist hosts also exist, such as dual-mycorrhizal plants that form symbiotic associations with both ectomycorrhizal fungi (EM) and arbuscular mycorrhizal fungi (AM). Little is known about the effect of dual mycorrhization status on the hosts' global distribution and acclimation to specific environments. This study investigates the potential advantage of dual associations of more than 400 woody genera spread at a global scale. We found that dual-host woody species occupy a broader geographical range and environmental niche space compared to those associating exclusively with either AM or EM. We show that the increased geographic range and expanded environmental niche space are independent of the phylogenetic architecture and evolutionary history of the woody genera. Our results highlight the advantage of generalist host-microbe symbioses between woody species and fungi to expand their range, and their potential role in colonising dry climates.

Conserved genetic basis for microbial colonization of the gut

Despite the fundamental importance of gut microbes, the genetic basis of their colonization remains largely unexplored. Here, by applying cross-species genotype-habitat association at the tree-of-life scale, we identify conserved microbial gene modules associated with gut colonization. Across thousands of species, we discovered 79 taxonomically diverse putative colonization factors organized into operonic and non-operonic modules. They include previously characterized colonization pathways such as autoinducer-2 biosynthesis and novel processes including tRNA modification and translation. In vivo functional validation revealed YigZ (IMPACT family) and tRNA hydroxylation protein-P (TrhP) are required for E. coli intestinal colonization. Overexpressing YigZ alone is sufficient to enhance colonization of the poorly colonizing MG1655 E. coli by >100-fold. Moreover, natural allelic variations in YigZ impact inter-strain colonization efficiency. Our findings highlight the power of large-scale comparative genomics in revealing the genetic basis of microbial adaptations. These broadly conserved colonization factors may prove critical for understanding gastrointestinal (GI) dysbiosis and developing therapeutics.

A Unified Method for Detecting Phylogenetic Signals in Continuous, Discrete, and Multiple Trait Combinations

Phylogenetic signals are widely used in ecological and evolutionary studies. Trait data used to detect phylogenetic signals can be continuous or discrete, but existing indices are designed for either type, not both. Moreover, most existing methods can only perform phylogenetic detection of individual traits, despite the fact that biological functions are often the result of interactions among multiple traits. Some attempts to detect phylogenetic signals across multiple trait combinations have employed alternative indicators, which may not align perfectly with the rigorous criteria for defining phylogenetic signals. In this study, we developed a new index (the M statistic) to detect phylogenetic signals for continuous traits, discrete traits, and multiple trait combinations. This capability is inherited from Gower's distance, which is used in the calculation of the M statistic to convert various types of traits into distances. The M statistic strictly adheres to the definition of phylogenetic signals and detects them by comparing these distances from phylogenies and traits. Using simulated data, we compared the performance of our new approach with that of existing commonly used indices. The results show that our method is not inferior to the existing methods. It performs well in handling continuous variables, discrete variables, and multiple trait combinations. We used trait data of turtles (Testudines) to demonstrate the utility of our new method. We suggest this new index as an original method for the detection of phylogenetic signals across various variable types. We provide an R package called "phylosignalDB" to facilitate all calculations.

Species that require long-day conditions to flower are not advancing their flowering phenology as fast as species without photoperiod requirements

BACKGROUND AND AIMS

Over the last few decades, many plant species have shown changes in phenology, such as the date on which they germinate, bud or flower. However, some species are changing more slowly than others, potentially owing to daylength (photoperiod) requirements.

METHODS

We combined data on flowering-time advancement with published records of photoperiod sensitivity to try to predict which species are advancing their flowering time. Data availability limited us to the Northern Hemisphere.

KEY RESULTS

Cross-species analyses showed that short-day plants advanced their flowering time by 1.4 days per decade and day-neutral plants by 0.9 days per decade, but long-day plants delayed their flowering by 0.2 days per decade. However, photoperiod-sensitivity status exhibited moderate phylogenetic conservation, and the differences in flowering-time advancement were not significant after phylogeny was accounted for. Both annual and perennial herbs were more likely to have long-day photoperiod cues than woody species, which were more likely to have short-day photoperiod cues.

CONCLUSIONS

Short-day plants are keeping up with plants that do not have photoperiod requirements, suggesting that daylength requirements do not hinder changes in phenology. However, long-day plants are not changing their phenology and might risk falling behind as competitors and pollinators adapt to climate change.

Dynamic quinone repertoire accompanied the diversification of energy metabolism in Pseudomonadota

It is currently unclear how Pseudomonadota, a phylum that originated around the time of the Great Oxidation Event, became one of the most abundant and diverse bacterial phyla on Earth, with metabolically versatile members colonizing a wide range of environments with different O2 concentrations. Here, we address this question by studying isoprenoid quinones, which are central components of energy metabolism covering a wide range of redox potentials. We demonstrate that a dynamic repertoire of quinone biosynthetic pathways accompanied the diversification of Pseudomonadota. The low potential menaquinone (MK) was lost in an ancestor of Pseudomonadota while the high potential ubiquinone (UQ) emerged. We show that the O2-dependent and O2-independent UQ pathways were both present in the last common ancestor of Pseudomonadota, and transmitted vertically. The O2-independent pathway has a conserved genetic organization and displays signs of positive regulation by the master regulator "fumarate and nitrate reductase" (FNR), suggesting a conserved role for UQ in anaerobiosis across Pseudomonadota. The O2-independent pathway was lost in some lineages but maintained in others, where it favoured a secondary reacquisition of low potential quinones (MK or rhodoquinone), which promoted diversification towards aerobic facultative and anaerobic metabolisms. Our results support that the ecological success of Pseudomonadota is linked to the acquisition of the largest known repertoire of quinones, which allowed adaptation to oxic niches as O2 levels increased on Earth, and subsequent diversification into anoxic or O2-fluctuating environments.

Biome evolution in subfamily Cercidoideae (Leguminosae): a tropical arborescent clade with a relictual depauperate temperate lineage

Some plant lineages remain within the same biome over time (biome conservatism), whereas others seem to adapt more easily to new biomes. The c. 398 species (14 genera) of subfamily Cercidoideae (Leguminosae or Fabaceae) are found in many biomes around the world, particularly in the tropical regions of South America, Asia and Africa, and display a variety of growth forms (small trees, shrubs, lianas and herbaceous perennials). Species distribution maps derived from cleaned occurrence records were compiled and compared with existing biome maps and with the literature to assign species to biomes. Rainforest (144 species), succulent (44 species), savanna (36 species), and temperate (10 species) biomes were found to be important in describing the global distribution of Cercidoideae, with many species occurring in more than one biome. Two phylogenetically isolated species-poor temperate (Cercis) and succulent (Adenolobus) biome lineages are sister to two broadly distributed species-rich tropical clades. Ancestral state reconstructions on a time-calibrated phylogeny suggest biome shifts occurred throughout the evolutionary history of the subfamily, with shifts between the succulent and rainforest biomes, from the rainforest to savanna, from the succulent to savanna biome, and one early occurring shift into (or from) the temperate biome. Of the 26 inferred shifts in biome, three are closely associated with a shift from the ancestral tree/shrub growth form to a liana or herbaceous perennial habit. Only three of the 13 inferred transcontinental dispersal events are associated with biome shifts. Overall, we find that biome shifts tend to occur within the same continent and that dispersals to new continents tend to occur within the same biome, but that nonetheless the biome-conserved and biogeographically structured Cercidoideae have been able to adapt to different environments through time.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40415-024-01058-z.

RERconverge Expansion: Using Relative Evolutionary Rates to Study Complex Categorical Trait Evolution

Comparative genomics approaches seek to associate molecular evolution with the evolution of phenotypes across a phylogeny. Many of these methods lack the ability to analyze non-ordinal categorical traits with more than two categories. To address this limitation, we introduce an expansion to RERconverge that associates shifts in evolutionary rates with the convergent evolution of categorical traits. The categorical RERconverge expansion includes methods for performing categorical ancestral state reconstruction, statistical tests for associating relative evolutionary rates with categorical variables, and a new method for performing phylogeny-aware permutations, "permulations", on categorical traits. We demonstrate our new method on a three-category diet phenotype, and we compare its performance to binary RERconverge analyses and two existing methods for comparative genomic analyses of categorical traits: phylogenetic simulations and a phylogenetic signal based method. We present an analysis of how the categorical permulations scale with the number of species and the number of categories included in the analysis. Our results show that our new categorical method outperforms phylogenetic simulations at identifying genes and enriched pathways significantly associated with the diet phenotypes and that the categorical ancestral state reconstruction drives an improvement in our ability to capture diet-related enriched pathways compared to binary RERconverge when implemented without user input on phenotype evolution. The categorical expansion to RERconverge will provide a strong foundation for applying the comparative method to categorical traits on larger data sets with more species and more complex trait evolution than have previously been analyzed.

Viral diversity and co-evolutionary dynamics across the ant phylogeny

Knowledge of viral biodiversity within insects, particularly within ants, is extremely limited with only a few environmental viruses from invasive ant species identified to date. This study documents and explores the viral communities in ants. We comprehensively profile the metagenomes of a phylogenetically broad group of 35 ant species with varied ecological traits and report the discovery of 3710 novel and unique ant-associated viral genomes. These previously unknown viruses discovered within this study constitute over 95% of all currently described ant viruses, significantly increasing our knowledge of the ant virosphere. The identified RNA and DNA viruses fill gaps in insect-associated viral phylogenies and uncover evolutionary histories characterized by both frequent host switching and co-divergence. Many ants also host diverse bacterial communities, and we discovered that approximately one-third of these new ant-associated viruses are bacteriophages. Two ecological categories, bacterial abundance in the host and habitat degradation are both correlated with ant viral diversity and help to structure viral communities within ants. These data demonstrate that the ant virosphere is remarkably diverse phylogenetically and genomically and provide a substantial foundation for studies in virus ecology and evolution within eukaryotes. We highlight the importance of studying insect-associated viruses in natural ecosystems in order to more thoroughly and effectively understand host-microbe evolutionary dynamics.

Plastome phylogenomics and new fossil evidence from Dominican amber shed light on the evolutionary history of the Neotropical fern genus Pecluma

PREMISE

Molecular studies based on chloroplast markers have questioned the monophyly of the fern genus Pecluma (Polypodioideae, Polypodiaceae), which has several species of Polypodium nested within it. We explored the delimitation of Pecluma and its biogeographic pattern by evaluating the phylogenetic position of four Polypodium species not sequenced thus far and integrating the first fossil evidence of Pecluma.

METHODS

Using herbarium material, we applied a genome-skimming approach to obtain a phylogenetic hypothesis of Polypodioideae; assessed the combination of character states observed in the fossil from Miocene Dominican amber using a previously published phylogeny of Polypodioideae based on four plastid markers as framework; calculated divergence times; and conducted an ancestral area estimation.

RESULTS

Within Polypodioideae, Pecluma was recovered as sister to Phlebodium. Three of the newly sequenced species-Polypodium otites, P. pinnatissimum, and P. ursipes-were recovered with maximum support within the Pecluma clade, whereas P. christensenii remained within Polypodium. The closest combination of character states of the fossil was found within Pecluma. Our biogeographic analyses suggest an Eocene origin of the genus in South America, with several subsequent Oligocene and Miocene colonization events to Mexico-Central America and to the West Indies.

CONCLUSIONS

Although the circumscription of Pecluma is still challenging, our results elucidate the origin and age of the genus. The newly described fossil, Pecluma hispaniolae sp. nov., supports the hypothesis that the epiphytic communities of the Greater Antilles exhibit a constant generic composition since the Miocene. We propose new combinations (Pecluma otites, Pecluma pinnatissima, and Pecluma ursipes) to accommodate three species previously classified in Polypodium.

Effects of individual traits vs. trait syndromes on assemblages of various herbivore guilds associated with central European Salix

Plants employ diverse anti-herbivore defences that can covary to form syndromes consisting of multiple traits. Such syndromes are hypothesized to impact herbivores more than individual defences. We studied 16 species of lowland willows occurring in central Europe and explored if their chemical and physical traits form detectable syndromes. We tested for phylogenetic trends in the syndromes and explored whether three herbivore guilds (i.e., generalist leaf-chewers, specialist leaf-chewers, and gallers) are affected more by the detected syndromes or individual traits. The recovered syndromes showed low phylogenetic signal and were mainly defined by investment in concentration, richness, or uniqueness of structurally related phenolic metabolites. Resource acquisition traits or inducible volatile organic compounds exhibited a limited correlation with the syndromes. Individual traits composing the syndromes showed various correlations to the assemblages of herbivores from the three studied guilds. In turn, we found some support for the hypothesis that defence syndromes are composed of traits that provide defence against various herbivores. However, individual traits rather than trait syndromes explained more variation for all studied herbivore assemblages. The detected negative correlations between various phenolics suggest that investment trade-offs may occur primarily among plant metabolites with shared metabolic pathways that may compete for their precursors. Moreover, several traits characterizing the recovered syndromes play additional roles in willows other than defence from herbivory. Taken together, our findings suggest that the detected syndromes did not solely evolve as an anti-herbivore defence.

Annual temperature, body size, and sexual size dimorphism in the evolution of Pyrgomorphidae

In many animal species, larger body size is positively correlated with male mating success and female fecundity. However, in the case of insects, in high seasonality environments, natural selection favors a faster maturation that decreases the risk of pre-reproductive death. However, this advantageous adaptation comes at a tradeoff, resulting in a reduction in body size. Maturation time is influenced by environmental factors, such as temperature, season length, and food availability during the rains. The geographic variation in these parameters provides an opportunity to study their impact on the adaptive evolution of body size in Pyrgomorphidae grasshoppers. These grasshoppers exhibit remarkable variation in body size and wing development and can be found in diverse plant communities across Africa, Asia, Australia, and tropical America. In this study, we utilized a phylogenetic approach to examine the evolution of body size, considering climatic factors, and the influence of sexual selection on size differences between males and females. We found a positive correlation between mean annual temperature and sexual size dimorphism (SSD). Remarkably, species exhibiting a strong bias toward larger females were found to be adapted to regions with higher temperatures. In the Pyrgomorphidae family, an intermediate body size was identified as the ancestral trait. Additionally, winged male and female grasshoppers were observed to be larger than their wingless counterparts. Despite the potential conflicting pressures on body size in males and females, these grasshoppers adhere to Rench's Rule, suggesting that sexual selection on males' body size may explain the evolution of SSD.

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

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

Decoupled evolution of ventral and dorsal scales in agamid lizards: ventral keels are associated with arboreality

Arboreality has evolved in all major vertebrate lineages and is often associated with morphological adaptations and increased diversification concomitant with accessing novel niche space. In squamate reptiles, foot, claw, and tail morphology are well-studied adaptations shown to be associated with transitions to arboreality. Here, we examined a less well understood trait-the keeled scale-in relation to microhabitat, climate, and diversification dynamics across a diverse lizard radiation, Agamidae. We found that the ancestral agamid had keeled dorsal but not ventral scales; further, dorsal and ventral keels are evolutionarily decoupled. Ventral keeled scales evolved repeatedly in association with arboreality and may be advantageous in reducing wear or by promoting interlocking when climbing. We did not find an association between keeled scales and diversification, suggesting keels do not allow finer-scale microhabitat partitioning observed in other arboreal-associated traits. We additionally found a relationship between keeled ventral scales and precipitation in terrestrial species where we posit that the keels may function to reduce scale degradation. Our results suggest that keeled ventral scales facilitated transitions to arboreality across agamid lizards, and highlight a need for future studies that explore their biomechanical function in relation to microhabitat and climate.

Resting posture drives the evolution of agonistic displays in bats

Agonistic displays are one of the most diverse social behaviors that have important functions in animal's life history. However, their origin and driving factors have largely been unexplored. Here, we evaluated agonistic displays of 71 bat species across 10 families and classified these displays into two categories: (a) boxing displays where a bat attacks its opponent with its wrist and thumb and (b) pushing displays where a bat uses its head or body to hit a rival. We estimated the strength of the phylogenetic signal of the agonistic displays, revealed their origin, and tested the potential evolutionary relationships between agonistic behaviors and body size or resting posture (free hanging vs. contact hanging where the bat is in contact with some surface). We found that agonistic displays were phylogenetically conserved and that boxing displays are the ancestral state. Moreover, we found that bats with a free-hanging resting posture were more likely to exhibit boxing displays than pushing displays. In addition, bats with longer forearms do not have a higher propensity for boxing displays. This study expands our limited knowledge of the evolution of agonistic displays and highlights the importance of resting posture as a driving force in the diversity of agonistic displays.

Bilirubin reductase shows host-specific associations in animal large intestines

Animal gastrointestinal tracts contain diverse metabolites, including various host-derived compounds that gut-associated microbes interact with. Here, we explore the diversity and evolution of bilirubin reductase, a bacterial enzyme that metabolizes the host-derived tetrapyrrole bilirubin, performing a key role in the animal heme degradation pathway. Through an analysis of the bilirubin reductase phylogeny and predicted structures, we found that the enzyme family can be divided into three distinct clades with different structural features. Using these clade definitions, we analyzed metagenomic sequencing data from multiple animal species, finding that bilirubin reductase is significantly enriched in the large intestines of animals and that the clades exhibit differences in distribution among animals. Combined with phylogenetic signal analysis, we find that the bilirubin reductase clades exhibit significant associations with specific animals and animal physiological traits like gastrointestinal anatomy and diet. These patterns demonstrate that bilirubin reductase is specifically adapted to the anoxic lower gut environment of animals and that its evolutionary history is complex, involving adaptation to a diverse collection of animals harboring bilirubin-reducing microbes. The findings suggest that bilirubin reductase evolution has been shaped by the host environment, providing a new perspective on heme metabolism in animals and highlighting the importance of the microbiome in animal physiology and evolution.

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.

RERconverge Expansion: Using Relative Evolutionary Rates to Study Complex Categorical Trait Evolution

Comparative genomics approaches seek to associate evolutionary genetic changes with the evolution of phenotypes across a phylogeny. Many of these methods, including our evolutionary rates based method, RERconverge, lack the capability of analyzing non-ordinal, multicategorical traits. To address this limitation, we introduce an expansion to RERconverge that associates shifts in evolutionary rates with the convergent evolution of multi-categorical traits. The categorical RERconverge expansion includes methods for performing categorical ancestral state reconstruction, statistical tests for associating relative evolutionary rates with categorical variables, and a new method for performing phylogenetic permulations on multi-categorical traits. In addition to demonstrating our new method on a three-category diet phenotype, we compare its performance to naive pairwise binary RERconverge analyses and two existing methods for comparative genomic analyses of categorical traits: phylogenetic simulations and a phylogenetic signal based method. We also present a diagnostic analysis of the new permulations approach demonstrating how the method scales with the number of species and the number of categories included in the analysis. Our results show that our new categorical method outperforms phylogenetic simulations at identifying genes and enriched pathways significantly associated with the diet phenotype and that the new ancestral reconstruction drives an improvement in our ability to capture diet-related enriched pathways. Our categorical permulations were able to account for non-uniform null distributions and correct for non-independence in gene rank during pathway enrichment analysis. The categorical expansion to RERconverge will provide a strong foundation for applying the comparative method to categorical traits on larger data sets with more species and more complex trait evolution.

Why are so many enteric pathogen infections asymptomatic? Pathogen and gut microbiome characteristics associated with diarrhea symptoms and carriage of diarrheagenic E. coli in northern Ecuador

A high proportion of enteric infections, including those caused by diarrheagenic Escherichia coli (DEC), are asymptomatic for diarrhea. The factors responsible for the development of diarrhea symptoms, or lack thereof, remain unclear. Here, we used DEC isolate genome and whole stool microbiome data from a case-control study of diarrhea in Ecuador to examine factors associated with diarrhea symptoms accompanying DEC carriage. We investigated i) pathogen abundance, ii) gut microbiome characteristics, and iii) strain-level pathogen characteristics from DEC infections with diarrhea symptoms (symptomatic infections) and without diarrhea symptoms (asymptomatic infections). We also included data from individuals with and without diarrhea who were not infected with DEC (uninfected cases and controls). i) E. coli relative abundance in the gut microbiome was highly variable, but higher on-average in individuals with symptomatic compared to asymptomatic DEC infections. Similarly, the number and relative abundances of virulence genes in the gut were higher in symptomatic than asymptomatic DEC infections. ii) Measures of microbiome diversity were similar regardless of diarrhea symptoms or DEC carriage. Proteobacterial families that have been described as pathobionts were enriched in symptomatic infections and uninfected cases, whereas potentially beneficial taxa, including the Bacteroidaceae and Bifidobacteriaceae, were more abundant in individuals without diarrhea. An analysis of high-level gene functions recovered in metagenomes revealed that genes that were differentially abundant by diarrhea and DEC infection status were more abundant in symptomatic than asymptomatic DEC infections. iii) DEC isolates from symptomatic versus asymptomatic individuals showed no significant differences in virulence or accessory gene content, and there was no phylogenetic signal associated with diarrhea symptoms. Together, these data suggest signals that distinguish symptomatic from asymptomatic DEC infections. In particular, the abundance of E. coli, the virulence gene content of the gut microbiome, and the taxa present in the gut microbiome have an apparent role.

Impact of the climatic changes in the Pliocene-Pleistocene transition on Irano-Turanian species. The radiation of genus Jurinea (Compositae)

The Irano-Turanian region is one of the world's richest floristic regions and the centre of diversity for numerous xerophytic plant lineages. However, we still have limited knowledge on the timing of evolution and biogeographic history of its flora, and potential drivers of diversification remain underexplored. To fill this knowledge gap, we focus on the Eurasian genus Jurinea (ca. 200 species), one of the largest plant radiations that diversified in the region. We applied a macroevolutionary integrative approach to explicitly test diversification hypotheses and investigate the relative roles of geography vs. ecology and niche conservatism vs. niche lability in speciation processes. To do so, we gathered a sample comprising 77% of total genus richness and obtained data about (1) its phylogenetic history, recovering 502 nuclear loci sequences; (2) growth forms; (3) ecological niche, compiling data of 21 variables for more than 2500 occurrences; and (4) paleoclimatic conditions, to estimate climatic stability. Our results revealed that climate was a key factor in the evolutionary dynamics of Jurinea. The main diversification and biogeographic events that occurred during past climate changes, which led to colder and drier conditions, are the following: (1) the origin of the genus (10.7 Ma); (2) long-distance dispersals from the Iranian Plateau to adjacent regions (∼7-4 Ma); and (3) the diversification shift during Pliocene-Pleistocene Transition (ca. 3 Ma), when net diversification rate almost doubled. Our results supported the pre-adaptation hypothesis, i.e., the evolutionary success of Jurinea was linked to the retention of the ancestral niche adapted to aridity. Interestingly, the paleoclimatic analyses revealed that in the Iranian Plateau long-term climatic stability favoured old-lineage persistence, resulting in current high species richness of semi-arid and cold adapted clades; whereas moderate climate oscillations stimulated allopatric diversification in the lineages distributed in the Circumboreal region. In contrast, growth form lability and high niche disparity among closely related species in the Central Asian clade suggest adaptive radiation to mountain habitats. In sum, the radiation of Jurinea is the result of both adaptive and non-adaptive processes influenced by climatic, orogenic and ecological factors.

Environmental modulation of plant mycorrhizal traits in the global flora

Mycorrhizal symbioses are known to strongly influence plant performance, structure plant communities and shape ecosystem dynamics. Plant mycorrhizal traits, such as those characterising mycorrhizal type (arbuscular (AM), ecto-, ericoid or orchid mycorrhiza) and status (obligately (OM), facultatively (FM) or non-mycorrhizal) offer valuable insight into plant belowground functionality. Here, we compile available plant mycorrhizal trait information and global occurrence data (∼ 100 million records) for 11,770 vascular plant species. Using a plant phylogenetic mega-tree and high-resolution climatic and edaphic data layers, we assess phylogenetic and environmental correlates of plant mycorrhizal traits. We find that plant mycorrhizal type is more phylogenetically conserved than plant mycorrhizal status, while environmental variables (both climatic and edaphic; notably soil texture) explain more variation in mycorrhizal status, especially FM. The previously underestimated role of environmental conditions has far-reaching implications for our understanding of ecosystem functioning under changing climatic and soil conditions.

Population genomics of diarrheagenic Escherichia coli uncovers high connectivity between urban and rural communities in Ecuador

Human movement may be an important driver of transmission dynamics for enteric pathogens but has largely been underappreciated except for international 'travelers' diarrhea or cholera. Phylodynamic methods, which combine genomic and epidemiological data, are used to examine rates and dynamics of disease matching underlying evolutionary history and biogeographic distributions, but these methods often are not applied to enteric bacterial pathogens. We used phylodynamics to explore the phylogeographic and evolutionary patterns of diarrheagenic E. coli in northern Ecuador to investigate the role of human travel in the geographic distribution of strains across the country. Using whole genome sequences of diarrheagenic E. coli isolates, we built a core genome phylogeny, reconstructed discrete ancestral states across urban and rural sites, and estimated migration rates between E. coli populations. We found minimal structuring based on site locations, urban vs. rural locality, pathotype, or clinical status. Ancestral states of phylogenomic nodes and tips were inferred to have 51% urban ancestry and 49% rural ancestry. Lack of structuring by location or pathotype E. coli isolates imply highly connected communities and extensive sharing of genomic characteristics across isolates. Using an approximate structured coalescent model, we estimated rates of migration among circulating isolates were 6.7 times larger for urban towards rural populations compared to rural towards urban populations. This suggests increased inferred migration rates of diarrheagenic E. coli from urban populations towards rural populations. Our results indicate that investments in water and sanitation prevention in urban areas could limit the spread of enteric bacterial pathogens among rural populations.

Human Amygdala Volumetric Patterns Convergently Evolved in Cooperatively Breeding and Domesticated Species

The amygdala is a hub in brain networks that supports social life and fear processing. Compared with other apes, humans have a relatively larger lateral nucleus of the amygdala, which is consistent with both the self-domestication and the cooperative breeding hypotheses of human evolution. Here, we take a comparative approach to the evolutionary origin of the relatively larger lateral amygdala nucleus in humans. We carry out phylogenetic analysis on a sample of 17 mammalian species for which we acquired single amygdala nuclei volumetric data. Our results indicate that there has been convergent evolution toward larger lateral amygdala nuclei in both domesticated and cooperatively breeding mammals. These results suggest that changes in processing fearful stimuli to reduce fear-induced aggression, which are necessary for domesticated and cooperatively breeding species alike, tap into the same neurobiological proximate mechanism. However, humans show changes not only in processing fearful stimuli but also in proactive prosociality. Since cooperative breeding, but not domestication, is also associated with increased proactive prosociality, a prominent role of the former during human evolution is more parsimonious, whereas self-domestication may have been involved as an additional stepping stone.

Evolution of posture in amniotes-Diving into the trabecular architecture of the femoral head

Extant amniotes show remarkable postural diversity. Broadly speaking, limbs with erect (strongly adducted, more vertically oriented) posture are found in mammals that are particularly heavy (graviportal) or show good running skills (cursorial), while crouched (highly flexed) limbs are found in taxa with more generalized locomotion. In Reptilia, crocodylians have a "semi-erect" (somewhat adducted) posture, birds have more crouched limbs and lepidosaurs have sprawling (well-abducted) limbs. Both synapsids and reptiles underwent a postural transition from sprawling to more erect limbs during the Mesozoic Era. In Reptilia, this postural change is prominent among archosauriforms in the Triassic Period. However, limb posture in many key Triassic taxa remains poorly known. In Synapsida, the chronology of this transition is less clear, and competing hypotheses exist. On land, the limb bones are subject to various stresses related to body support that partly shape their external and internal morphology. Indeed, bone trabeculae (lattice-like bony struts that form the spongy bone tissue) tend to orient themselves along lines of force. Here, we study the link between femoral posture and the femoral trabecular architecture using phylogenetic generalized least squares. We show that microanatomical parameters measured on bone cubes extracted from the femoral head of a sample of amniote femora depend strongly on body mass, but not on femoral posture or lifestyle. We reconstruct ancestral states of femoral posture and various microanatomical parameters to study the "sprawling-to-erect" transition in reptiles and synapsids, and obtain conflicting results. We tentatively infer femoral posture in several hypothetical ancestors using phylogenetic flexible discriminant analysis from maximum likelihood estimates of the microanatomical parameters. In general, the trabecular network of the femoral head is not a good indicator of femoral posture. However, ancestral state reconstruction methods hold great promise for advancing our understanding of the evolution of posture in amniotes.

Profile of a flower: How rates of morphological evolution drive floral diversification in Ericales and angiosperms

PREMISE

Recent studies of floral disparity in the asterid order Ericales have shown that flowers vary strongly among families and that disparity is unequally distributed between the three flower modules (perianth, androecium, gynoecium). However, it remains unknown whether these patterns are driven by heterogeneous rates of morphological evolution or other factors.

METHODS

Here, we compiled a data set of 33 floral characters scored for 414 species of Ericales sampled from 346 genera and all 22 families. We conducted ancestral state reconstructions using an equal-rates Markov model for each character. We estimated rates of morphological evolution for Ericales and for a separate angiosperm-wide data set of 19 characters and 792 species, creating "rate profiles" for Ericales, angiosperms, and major angiosperm subclades. We compared morphological rates among flower modules within each data set separately and between data sets, and we compared rates among angiosperm subclades using the angiosperm data set.

RESULTS

The androecium exhibits the highest evolutionary rates across most characters, whereas most perianth and gynoecium characters evolve more slowly in both Ericales and angiosperms. Both high and low rates of morphological evolution can result in high floral disparity in Ericales. Analyses of an angiosperm-wide floral data set reveal that this pattern appears to be conserved across most major angiosperm clades.

CONCLUSIONS

Elevated rates of morphological evolution in the androecium of Ericales may explain the higher disparity reported for this floral module. Comparing rates of morphological evolution through rate profiles proves to be a powerful tool in understanding floral evolution.

Parallel evolution of fish bi-modal breathing and expansion of olfactory receptor (OR) genes: toward a universal ORs nomenclature

Olfactory receptors (ORs) play a key role in the prime sensorial perception, being highly relevant for intra/interspecific interactions. ORs are a subgroup of G-protein coupled receptors that exhibit highly complex subgenomes in vertebrates. However, OR repertoires remain poorly studied in fish lineages, precluding finely retracing their origin, evolution, and diversification, especially in the most basal groups. Here, we conduct an exhaustive gene screening upon 43 high-quality fish genomes exhibiting varied gene repertoires (2-583 genes). While the early vertebrates performed gas exchange through gills, we hypothesize that the emergence of new breathing structures (swim bladder and paired lungs) in early osteichthyans may be associated with expansions in the ORs gene families sensitive to airborne molecules. Additionally, we verify that the OR repertoire of moderns actinopterygians has not increased as expected following a whole genome duplication, likely due to regulatory mechanisms compensating the gene load excess. Finally, we identify 25 distinct OR families, allowing us to propose an updated universal nomenclature for the fish ORs.

Ultra-deep sequencing of Hadza hunter-gatherers recovers vanishing gut microbes

The gut microbiome modulates immune and metabolic health. Human microbiome data are biased toward industrialized populations, limiting our understanding of non-industrialized microbiomes. Here, we performed ultra-deep metagenomic sequencing on 351 fecal samples from the Hadza hunter-gatherers of Tanzania and comparative populations in Nepal and California. We recovered 91,662 genomes of bacteria, archaea, bacteriophages, and eukaryotes, 44% of which are absent from existing unified datasets. We identified 124 gut-resident species vanishing in industrialized populations and highlighted distinct aspects of the Hadza gut microbiome related to in situ replication rates, signatures of selection, and strain sharing. Industrialized gut microbes were found to be enriched in genes associated with oxidative stress, possibly a result of microbiome adaptation to inflammatory processes. This unparalleled view of the Hadza gut microbiome provides a valuable resource, expands our understanding of microbes capable of colonizing the human gut, and clarifies the extensive perturbation induced by the industrialized lifestyle.

Evolutionary history explains foliar spectral differences between arbuscular and ectomycorrhizal plant species

Leaf spectra are integrated foliar phenotypes that capture a range of traits and can provide insight into ecological processes. Leaf traits, and therefore leaf spectra, may reflect belowground processes such as mycorrhizal associations. However, evidence for the relationship between leaf traits and mycorrhizal association is mixed, and few studies account for shared evolutionary history. We conduct partial least squares discriminant analysis to assess the ability of spectra to predict mycorrhizal type. We model the evolution of leaf spectra for 92 vascular plant species and use phylogenetic comparative methods to assess differences in spectral properties between arbuscular mycorrhizal and ectomycorrhizal plant species. Partial least squares discriminant analysis classified spectra by mycorrhizal type with 90% (arbuscular) and 85% (ectomycorrhizal) accuracy. Univariate models of principal components identified multiple spectral optima corresponding with mycorrhizal type due to the close relationship between mycorrhizal type and phylogeny. Importantly, we found that spectra of arbuscular mycorrhizal and ectomycorrhizal species do not statistically differ from each other after accounting for phylogeny. While mycorrhizal type can be predicted from spectra, enabling the use of spectra to identify belowground traits using remote sensing, this is due to evolutionary history and not because of fundamental differences in leaf spectra due to mycorrhizal type.

Geometric regularity in webs of non-orb-weaving spiders

Geometric regularity of spider webs has been intensively studied in orb-weaving spiders, although it is not exclusive of orb weavers. Here, we document the geometrically regular, repetitive elements in the webs of the non-orb-weaving groups Leptonetidae and Telemidae for the first time. Similar to orb weavers, we found areas with regularly spaced parallel lines in the webs of Calileptoneta helferi, Sulcia sp., and cf. Pinelema sp. Furthermore, we provide a detailed account of the regular webs of Ochyrocera (Ochyroceratidae). The sections of the web with regularly disposed parallel lines are built as U-shaped modules reminiscent of orb webs. It has been suggested that the regularly spaced parallel lines in the webs of Ochyroceratidae and Psilodercidae may be produced in a single sweep of their posterior lateral spinnerets, which have regularly spaced aciniform gland spigots, perhaps involving expansion of the spinnerets. To test this hypothesis, we compared the spacing between parallel lines with the spacing between spigots, searched for expansible membranes in the spinnerets, and examined the junctions of regularly spaced lines. The distance between parallel lines was 10-20 times the distance between spigots, and we found no expansible membranes, and the intersection of parallel lines are cemented, which opposes the single sweep hypothesis. Furthermore, we found cues of viscid silk in the parallel lines of the psilodercid Althepus and broadened piriform gland spigots that may be responsible of its production. Finally, we evaluated the presence or absence of geometrically regular web elements across the spider tree of life. We found reports of regular webs in 31 spider families, including 20 families that are not orb weavers and hypothesize that the two basic aspects of regularity (parallel lines spaced at regular intervals, and radial lines spaced at regular angles) probably appeared many times in the evolution of spiders.

Locomotor and postural diversity among reptiles viewed through the prism of femoral microanatomy: Palaeobiological implications for some Permian and Mesozoic taxa

The water-to-land transition by the first tetrapod vertebrates represents a key stage in their evolution. Selection pressures exerted by this new environment on animals led to the emergence of new locomotor and postural strategies that favoured access to different ecological niches and contributed to their evolutionary success. Today, amniotes show great locomotor and postural diversity, particularly among Reptilia, whose extant representatives include parasagittally locomoting erect and crouched bipeds (birds), sub-parasagittal 'semi-erect' quadrupeds (crocodylians) and sprawling quadrupeds (squamates and turtles). But the different steps leading to such diversity remain enigmatic and the type of locomotion adopted by many extinct species raises questions. This is notably the case of certain Triassic taxa such as Euparkeria and Marasuchus. The exploration of the bone microanatomy in reptiles could help to overcome these uncertainties. Indeed, this locomotor and postural diversity is accompanied by great microanatomical disparity. On land, the bones of the appendicular skeleton support the weight of the body and are subject to multiple constraints that partly shape their external and internal morphology. Here we show how microanatomical parameters measured in cross-section, such as bone compactness or the position of the medullocortical transition, can be related to locomotion. We hypothesised that this could be due to variations in cortical thickness. Using statistical methods that take phylogeny into account (phylogenetic flexible discriminant analyses), we develop different models of locomotion from a sample of femur cross-sections from 51 reptile species. We use these models to infer locomotion and posture in 7 extinct reptile taxa for which they remain debated or not fully clear. Our models produced reliable inferences for taxa that preceded and followed the quadruped/biped and sprawling/erect transitions, notably within the Captorhinidae and Dinosauria. For taxa contemporary with these transitions, such as Terrestrisuchus and Marasuchus, the inferences are more questionable. We use linear models to investigate the effect of body mass and functional ecology on our inference models. We show that body mass seems to significantly impact our model predictions in most cases, unlike the functional ecology. Finally, we illustrate how taphonomic processes can impact certain microanatomical parameters, especially the eccentricity of the section, while addressing some other potential limitations of our methods. Our study provides insight into the evolution of enigmatic locomotion in various early reptiles. Our models and methods could be used by palaeontologists to infer the locomotion and posture in other extinct reptile taxa, especially when considered in combination with other lines of evidence.

Phylotranscriptomics of Swertiinae (Gentianaceae) reveals that key floral traits are not phylogenetically correlated

Establishing how lineages with similar traits are phylogenetically related remains critical for understanding the origin of biodiversity on Earth. Floral traits in plants are widely used to explore phylogenetic relationships and to delineate taxonomic groups. The subtribe Swertiinae (Gentianaceae) comprises more than 350 species with high floral diversity ranging from rotate to tubular corollas and possessing diverse nectaries. Here we performed phylogenetic analysis of 60 species from all 15 genera of the subtribe Swertiinae sensu Ho and Liu, representing the range of floral diversity, using data from the nuclear and plastid genomes. Extensive topological conflicts were present between the nuclear and plastome trees. Three of the 15 genera represented by multiple species are polyphyletic in both trees. Key floral traits including corolla type, absence or presence of lobe scales, nectary type, nectary position, and stigma type are randomly distributed in the nuclear and plastome trees without phylogenetic correlation. We also revealed the likely ancient hybrid origin of one large clade comprising 10 genera with diverse floral traits. These results highlight the complex evolutionary history of this subtribe. The phylogenies constructed here provide a basic framework for further exploring the ecological and genetic mechanisms underlying both species diversification and floral diversity.

A comprehensive list of bacterial pathogens infecting humans

There exists an enormous diversity of bacteria capable of human infection, but no up-to-date, publicly accessible list is available. Combining a pragmatic definition of pathogenicity with an extensive search strategy, we report 1513 bacterial pathogens known to infect humans described pre-2021. Of these, 73 % were regarded as established (have infected at least three persons in three or more references) and 27 % as putative (fewer than three known cases). Pathogen species belong to 10 phyla and 24 classes scattered throughout the bacterial phylogeny. We show that new human pathogens are discovered at a rapid rate. Finally, we discuss how our results could be expanded to a database, which could provide a useful resource for microbiologists. Our list is freely available and archived on GitHub and Zenodo and we have provided walkthroughs to facilitate access and use.

Data collected by citizen scientists reveal the role of climate and phylogeny on the frequency of shelter types used by frogs across the Americas

Shelters are microhabitats where animals rest and hide. These microhabitats can be used from short daily periods to long-term estivation or hibernation. Environmental conditions and the phenotypical characteristics of the animal drive habitat selection in relation to shelters. Based on this, climate regions and phylogeny are expected to affect the use of different shelter types. Although shelters are yet to be described for most anuran species, a variety of microhabitats have already been reported as shelter-sites, including dense vegetation, rock crevices, and holes in the ground. In this study, we evaluated photos of frogs for sheltering behaviour from 29 countries in the Americas deposited on the popular citizen-science platform, iNaturalist. We compared the frequency of use of different shelter types identified on the photos among different climate regions and anuran families, also testing possible phylogenetic signals. We identified 11,133 photographs of 378 frog species showing individuals hiding in shelters or in a resting position. We classified observations into 10 shelter types, with live vegetation (24.7 %) being the most commonly recorded natural shelter, followed by hole in the ground (11.4 %) and tree trunk (11.1 %). The use of different shelter types varied between arid and humid climates, and also among different anuran families. We found strong phylogenetic signal for three shelter types (hole in the ground, live vegetation, and water) and the differences in shelter use among taxa suggest a relation with body characteristics. Approximately 47 % of observations of threatened and near threatened species were in hole in the ground, while artificial habitat represented only 3.6 % of the observations in this group. The daily pattern of shelter use corroborated the nocturnal activity of most species. Our findings also expanded the description of shelter sites for 330 species that had no published information on this behaviour. This study contributes to our current knowledge about animal behaviour and highlights the use of citizen science as an effective approach to understand the natural history of amphibians at a large scale.

Bacterial composition in Swedish raw drinking water reveals three major interacting ubiquitous metacommunities

BACKGROUND

Surface raw water used as a source for drinking water production is a critical resource, sensitive to contamination. We conducted a study on Swedish raw water sources, aiming to identify mutually co-occurring metacommunities of bacteria, and environmental factors driving such patterns.

METHODS

The water sources were different regarding nutrient composition, water quality, and climate characteristics, and displayed various degrees of anthropogenic impact. Water inlet samples were collected at six drinking water treatment plants over 3 years, totaling 230 samples. The bacterial communities of DNA sequenced samples (n = 175), obtained by 16S metabarcoding, were analyzed using a joint model for taxa abundance.

RESULTS

Two major groups of well-defined metacommunities of microorganisms were identified, in addition to a third, less distinct, and taxonomically more diverse group. These three metacommunities showed various associations to the measured environmental data. Predictions for the well-defined metacommunities revealed differing sets of favored metabolic pathways and life strategies. In one community, taxa with methanogenic metabolism were common, while a second community was dominated by taxa with carbohydrate and lipid-focused metabolism.

CONCLUSION

The identification of ubiquitous persistent co-occurring bacterial metacommunities in freshwater habitats could potentially facilitate microbial source tracking analysis of contamination issues in freshwater sources.

Evolution of the rostral scale and mimicry in the genus Xenodon Boie, 1826 (Serpentes: Dipsadidae: Xenodontinae)

Snakes are a stimulating life form from an evolutionary perspective. Despite the basic morphological body shape (limbless, with a tubular body), these vertebrates are extremely diverse. The Neotropical region is one of the most diverse regions for snakes in the world, with >650 known species. Within this great diversity, the genus Xenodon includes 12 species with interesting adaptations to terrestrial and semi-fossorial habitats. Members of this genus are mostly diurnal and terrestrial, feed mainly on anurans and exhibit Batesian mimicry of venomous snakes of the genera Bothrops or Micrurus. Here, through phylogenetic analysis and ancestral state estimation, we explore the evolution of the rostral scale and mimicry within the genus Xenodon. Our results suggest that the ancestral lineage of Xenodon had a rounded rostral scale and exhibited Bothrops mimicry. The evolution of the rostral scale in Xenodon might be related to abiotic factors, as an adaptation for open and forested habitats, and mimicry is likely to be related to biotic factors, as a defensive strategy resembling those of venomous snakes.

Floral evolution and pollinator diversification in Hedychium: Revisiting Darwin's predictions using an integrative taxonomic approach

PREMISE

Hedychium J. Koenig (Zingiberaceae) is endemic to the Indo-Malayan Realm and is known for its colorful and fragrant flowers. Historically, two different pollination syndromes characterize Hedychium: diurnal or bird pollination, and nocturnal or moth pollination. In this study, we aim to understand the evolution of nocturnal and diurnal flowers, and to test its putative association with lineage diversification in Hedychium.

METHODS

A molecular tree of Hedychium was used as a scaffold upon which we estimated ancestral character states, phylogenetic signals, and correlations for certain categorical and continuous floral traits. Furthermore, we used phylomorphospace and trait-dependent diversification rate estimation analyses to understand phenotypic evolution and associated lineage diversification in Hedychium.

RESULTS

Although floral color and size lacked any association with specific pollinators, white or pale flowers were most common in the early branching clades when compared to bright-colored flowers, which were more widely represented in the most-derived clade IV. Five categorical and two continuous characters were identified to have informative evolutionary patterns, which also emphasized that ecology may have played a critical role in the diversification of Hedychium.

CONCLUSIONS

From our phylogenetic analyses and ecological observations, we conclude that specializations in pollinator interactions are rare in the hyperdiverse clade IV, thus challenging the role of both moth-specialization and bird-specialization as central factors in the diversification of Hedychium. However, our results also suggest that clade III (predominantly island clade) may show specializations, and future studies should investigate ecological and pollinator interactions, along with inclusion of new traits such as floral fragrance and anthesis time.

Genome diploidization associates with cladogenesis, trait disparity, and plastid gene evolution

Angiosperm genome evolution was marked by many clade-specific whole-genome duplication events. The Microlepidieae is one of the monophyletic clades in the mustard family (Brassicaceae) formed after an ancient allotetraploidization. Postpolyploid cladogenesis has resulted in the extant c. 17 genera and 60 species endemic to Australia and New Zealand (10 species). As postpolyploid genome diploidization is a trial-and-error process under natural selection, it may proceed with different intensity and be associated with speciation events. In Microlepidieae, different extents of homoeologous recombination between the two parental subgenomes generated clades marked by slow ("cold") versus fast ("hot") genome diploidization. To gain a deeper understanding of postpolyploid genome evolution in Microlepidieae, we analyzed phylogenetic relationships in this tribe using complete chloroplast sequences, entire 35S rDNA units, and abundant repetitive sequences. The four recovered intra-tribal clades mirror the varied diploidization of Microlepidieae genomes, suggesting that the intrinsic genomic features underlying the extent of diploidization are shared among genera and species within one clade. Nevertheless, even congeneric species may exert considerable morphological disparity (e.g. in fruit shape), whereas some species within different clades experience extensive morphological convergence despite the different pace of their genome diploidization. We showed that faster genome diploidization is positively associated with mean morphological disparity and evolution of chloroplast genes (plastid-nuclear genome coevolution). Higher speciation rates in perennials than in annual species were observed. Altogether, our results confirm the potential of Microlepidieae as a promising subject for the analysis of postpolyploid genome diploidization in Brassicaceae.

Quantifying research interests in 7,521 mammalian species with h-index: a case study

BACKGROUND

Taxonomic bias is a known issue within the field of biology, causing scientific knowledge to be unevenly distributed across species. However, a systematic quantification of the research interest that the scientific community has allocated to individual species remains a big data problem. Scalable approaches are needed to integrate biodiversity data sets and bibliometric methods across large numbers of species. The outputs of these analyses are important for identifying understudied species and directing future research to fill these gaps.

FINDINGS

In this study, we used the species h-index to quantity the research interest in 7,521 species of mammals. We tested factors potentially driving species h-index, by using a Bayesian phylogenetic generalized linear mixed model (GLMM). We found that a third of the mammals had a species h-index of zero, while a select few had inflated research interest. Further, mammals with higher species h-index had larger body masses; were found in temperate latitudes; had their humans uses documented, including domestication; and were in lower-risk International Union for Conservation of Nature Red List categories. These results surprisingly suggested that critically endangered mammals are understudied. A higher interest in domesticated species suggested that human use is a major driver and focus in mammalian scientific literature.

CONCLUSIONS

Our study has demonstrated a scalable workflow and systematically identified understudied species of mammals, as well as identified the likely drivers of this taxonomic bias in the literature. This case study can become a benchmark for future research that asks similar biological and meta-research questions for other taxa.

Macroevolutionary dynamics in the transition of angiosperms to aquatic environments

Angiosperm lineages in aquatic environments are characterized by high structural and functional diversity, and wide distributions. A long-standing evolutionary riddle is what processes have caused the relatively low diversity of aquatic angiosperms compared to their terrestrial relatives. We use diversification and ancestral reconstruction models with a comprehensive > 10 000 genus angiosperm phylogeny to elucidate the macroevolutionary dynamics associated with transitions of terrestrial plants to water. Our study reveals that net diversification rates are significantly lower in aquatic than in terrestrial angiosperms due to lower speciation and higher extinction. Shifts from land to water started early in angiosperm evolution, but most events were concentrated during the last c. 25 million years. Reversals to a terrestrial habitat started only 40 million years ago, but occurred at much higher rates. Within aquatic angiosperms, the estimated pattern is one of gradual accumulation of lineages, and relatively low and constant diversification rates throughout the Cenozoic. Low diversification rates, together with infrequent water transitions, account for the low diversity of aquatic angiosperms today. The stressful conditions and small global surface of the aquatic habitat available for angiosperms are hypothesized to explain this pattern.

Resident birds are more behaviourally plastic than migrants

Species subjected to more variable environments should have greater phenotypic plasticity than those that are more restricted to specific habitat types leading to the expectation that migratory birds should be relatively more plastic than resident birds. We tested this comparatively by studying variation in flight initiation distance (FID), a well-studied antipredator behaviour. We predicted that variation in FID would be greater for migratory species because they encountered a variety of locations during their lives and therefore had less predictable assessments of risk compared to more sedentary species. Contrary to our prediction, we found that non-migratory species (sedentary) had greater variation in FID than migratory ones. Migratory and partially migratory birds had greater average FIDs than sedentary birds, suggesting that they were generally more wary. These results suggest that the predictability associated with not migrating permits more nuanced risk assessment which was seen in the greater variation in FID of sedentary bird species.

Ultra-deep Sequencing of Hadza Hunter-Gatherers Recovers Vanishing Microbes.. [PMID: 36238714 PMCID: PMC9558438 DOI: 10.1101/2022.03.30.486478]

The gut microbiome is a key modulator of immune and metabolic health. Human microbiome data is biased towards industrialized populations, providing limited understanding of the distinct and diverse non-industrialized microbiomes. Here, we performed ultra-deep metagenomic sequencing and strain cultivation on 351 fecal samples from the Hadza, hunter-gatherers in Tanzania, and comparative populations in Nepal and California. We recover 94,971 total genomes of bacteria, archaea, bacteriophages, and eukaryotes, 43% of which are absent from existing unified datasets. Analysis of in situ growth rates, genetic pN/pS signatures, high-resolution strain tracking, and 124 gut-resident species vanishing in industrialized populations reveals differentiating dynamics of the Hadza gut microbiome. Industrialized gut microbes are enriched in genes associated with oxidative stress, possibly a result of microbiome adaptation to inflammatory processes. This unparalleled view of the Hadza gut microbiome provides a valuable resource that expands our understanding of microbes capable of colonizing the human gut and clarifies the extensive perturbation brought on by the industrialized lifestyle.

Assessing Bayesian Phylogenetic Information Content of Morphological Data Using Knowledge from Anatomy Ontologies

Morphology remains a primary source of phylogenetic information for many groups of organisms, and the only one for most fossil taxa. Organismal anatomy is not a collection of randomly assembled and independent “parts”, but instead a set of dependent and hierarchically nested entities resulting from ontogeny and phylogeny. How do we make sense of these dependent and at times redundant characters? One promising approach is using ontologies—structured controlled vocabularies that summarize knowledge about different properties of anatomical entities, including developmental and structural dependencies. Here, we assess whether evolutionary patterns can explain the proximity of ontology-annotated characters within an ontology. To do so, we measure phylogenetic information across characters and evaluate if it matches the hierarchical structure given by ontological knowledge—in much the same way as across-species diversity structure is given by phylogeny. We implement an approach to evaluate the Bayesian phylogenetic information (BPI) content and phylogenetic dissonance among ontology-annotated anatomical data subsets. We applied this to data sets representing two disparate animal groups: bees (Hexapoda: Hymenoptera: Apoidea, 209 chars) and characiform fishes (Actinopterygii: Ostariophysi: Characiformes, 463 chars). For bees, we find that BPI is not substantially explained by anatomy since dissonance is often high among morphologically related anatomical entities. For fishes, we find substantial information for two clusters of anatomical entities instantiating concepts from the jaws and branchial arch bones, but among-subset information decreases and dissonance increases substantially moving to higher-level subsets in the ontology. We further applied our approach to address particular evolutionary hypotheses with an example of morphological evolution in miniature fishes. While we show that phylogenetic information does match ontology structure for some anatomical entities, additional relationships and processes, such as convergence, likely play a substantial role in explaining BPI and dissonance, and merit future investigation. Our work demonstrates how complex morphological data sets can be interrogated with ontologies by allowing one to access how information is spread hierarchically across anatomical concepts, how congruent this information is, and what sorts of processes may play a role in explaining it: phylogeny, development, or convergence. [Apidae; Bayesian phylogenetic information; Ostariophysi; Phenoscape; phylogenetic dissonance; semantic similarity.]

A milk-sharing economy allows placental mammals to overcome their metabolic limits

Here, we demonstrate that a naturally evolving behavior (allonursing) has greater effect on reproductive power (mass per unit of time) and output (litter mass at birth) than does artificial selection (domestication). Additionally, we demonstrate the importance of resource optimization afforded by sociality (rather than resource abundance per se) in shaping a species’ life history profile and its ability to overcome its own physiological constraints.

Maternal resource availability and metabolism have a strong limiting effect on reproductive output. Allomaternal care and domestication increase the energy available to the mother and should correlate with an increase in reproductive output. Here, we take a comparative approach to understand how this increase is accomplished (e.g., litter mass, reproductive frequency, etc.) and the strength of the effect among different forms of external energetic supplementation. We find that domestication and all forms of allocare correlate with increased fertility. All forms of provisioning correlate with larger litters without compromising offspring size. The greatest increase we observe in reproductive power is in species that practice allonursing. Our results suggest that the ultimate factor limiting reproductive output in placental mammals is maternal metabolic power rather than resource availability.

Mimicry can drive convergence in structural and light transmission features of transparent wings in Lepidoptera

Müllerian mimicry is a positive interspecific interaction, whereby co-occurring defended prey species share a common aposematic signal. In Lepidoptera, aposematic species typically harbour conspicuous opaque wing colour patterns with convergent optical properties among co-mimetic species. Surprisingly, some aposematic mimetic species have partially transparent wings, raising the questions of whether optical properties of transparent patches are also convergent, and of how transparency is achieved. Here, we conducted a comparative study of wing optics, micro and nanostructures in neotropical mimetic clearwing Lepidoptera, using spectrophotometry and microscopy imaging. We show that transparency, as perceived by predators, is convergent among co-mimics in some mimicry rings. Underlying micro- and nanostructures are also sometimes convergent despite a large structural diversity. We reveal that while transparency is primarily produced by microstructure modifications, nanostructures largely influence light transmission, potentially enabling additional fine-tuning in transmission properties. This study shows that transparency might not only enable camouflage but can also be part of aposematic signals.

Prevalence and Specificity of Chemoreceptor Profiles in Plant-Associated Bacteria

Chemosensory pathways are among the most abundant prokaryotic signal transduction systems, allowing bacteria to sense and respond to environmental stimuli. Signaling is typically initiated by the binding of specific molecules to the ligand binding domain (LBD) of chemoreceptor proteins (CRs). Although CRs play a central role in plant-microbiome interactions such as colonization and infection, little is known about their phylogenetic and ecological specificity. Here, we analyzed 82,277 CR sequences from 11,806 representative microbial species covering the whole prokaryotic phylogeny, and we classified them according to their LBD type using a de novo homology clustering method. Through phylogenomic analysis, we identified hundreds of LBDs that are found predominantly in plant-associated bacteria, including several LBDs specific to phytopathogens and plant symbionts. Functional annotation of our catalogue showed that many of the LBD clusters identified might constitute unknown types of LBDs. Moreover, we found that the taxonomic distribution of most LBD types that are specific to plant-associated bacteria is only partially explained by phylogeny, suggesting that lifestyle and niche adaptation are important factors in their selection. Finally, our results show that the profile of LBD types in a given genome is related to the lifestyle specialization, with plant symbionts and phytopathogens showing the highest number of niche-specific LBDs. The LBD catalogue and information on how to profile novel genomes are available at https://github.com/compgenomicslab/CRs.

IMPORTANCE Considering the enormous variety of LBDs at sensor proteins, an important question resides in establishing the forces that have driven their evolution and selection. We present here the first clear demonstration that environmental factors play an important role in the selection and evolution of LBDs. We were able to demonstrate the existence of LBD families that are highly enriched in plant-associated bacteria but show a wide phylogenetic spread. These findings offer a number of research opportunities in the field of single transduction, such as the exploration of similar relationships in chemoreceptors of bacteria with a different lifestyle, like those inhabiting or infecting the human intestine. Similarly, our results raise the question whether similar LBD types might be shared by members of different sensor protein families. Lastly, we provide a comprehensive catalogue of CRs classified by their LBD region that includes a large number of putative new LBD types.

Horizontal Transfer of Microbial Toxin Genes to Gall Midge Genomes

A growing body of evidence has underscored the role of horizontal gene transfer (HGT) in animal evolution. Previously, we discovered the horizontal transfer of the gene encoding the eukaryotic genotoxin cytolethal distending toxin B (cdtB) from the pea aphid Acyrthosiphon pisum secondary endosymbiont (APSE) phages to drosophilid and aphid nuclear genomes. Here, we report cdtB in the nuclear genome of the gall-forming "swede midge" Contarinia nasturtii (Diptera: Cecidomyiidae) via HGT. We searched all available gall midge genome sequences for evidence of APSE-to-insect HGT events and found five toxin genes (aip56, cdtB, lysozyme, rhs, and sltxB) transferred horizontally to cecidomyiid nuclear genomes. Surprisingly, phylogenetic analyses of HGT candidates indicated APSE phages were often not the ancestral donor lineage of the toxin gene to cecidomyiids. We used a phylogenetic signal statistic to test a transfer-by-proximity hypothesis for animal HGT, which suggested that microbe-to-insect HGT was more likely between taxa that share environments than those from different environments. Many of the toxins we found in midge genomes target eukaryotic cells, and catalytic residues important for toxin function are conserved in insect copies. This class of horizontally transferred, eukaryotic cell-targeting genes is potentially important in insect adaptation.