Phylogenetic Models Research Articles

Exploring the Macroevolutionary Signature of Asymmetric Inheritance at Speciation.

Popular comparative phylogenetic models such as Brownian Motion, Ornstein-Ulhenbeck, and their extensions, assume that, at speciation, a trait value is inherited identically by two descendant species. This assumption contrasts with models of speciation at a micro-evolutionary scale where descendants' phenotypic distributions are sub-samples of the ancestral distribution. Different speciation mechanisms can lead to a displacement of the ancestral phenotypic mean among descendants and an asymmetric inheritance of the ancestral phenotypic variance. In contrast, even macro-evolutionary models that account for intraspecific variance assume symmetrically conserved inheritance of ancestral phenotypic distribution at speciation. Here we develop an Asymmetric Brownian Motion model (ABM) that relaxes the assumption of symmetric and conserved inheritance of the ancestral distribution at the time of speciation. The ABM jointly models the evolution of both intra- and inter-specific phenotypic variation. It also infers the mode of phenotypic inheritance at speciation, which can range from a symmetric and conserved inheritance, where descendants inherit the ancestral distribution, to an asymmetric and displaced inheritance, where descendants inherit divergent phenotypic means and variances. To demonstrate this model, we analyze the evolution of beak morphology in Darwin finches, finding evidence of displacement at speciation. The ABM model helps to bridge micro- and macro-evolutionary models of trait evolution by providing a more robust framework for testing the effects of ecological speciation, character displacement, and niche partitioning on trait evolution at the macro-evolutionary scale.

Maximum Likelihood Estimation for Unrooted 3-Leaf Trees: An Analytic Solution for the CFN Model

Maximum likelihood estimation is among the most widely-used methods for inferring phylogenetic trees from sequence data. This paper solves the problem of computing solutions to the maximum likelihood problem for 3-leaf trees under the 2-state symmetric mutation model (CFN model). Our main result is a closed-form solution to the maximum likelihood problem for unrooted 3-leaf trees, given generic data; this result characterizes all of the ways that a maximum likelihood estimate can fail to exist for generic data and provides theoretical validation for predictions made in Parks and Goldman (Syst Biol 63(5):798–811, 2014). Our proof makes use of both classical tools for studying group-based phylogenetic models such as Hadamard conjugation and reparameterization in terms of Fourier coordinates, as well as more recent results concerning the semi-algebraic constraints of the CFN model. To be able to put these into practice, we also give a complete characterization to test genericity.

Fast mvSLOUCH: Multivariate Ornstein–Uhlenbeck‐based models of trait evolution on large phylogenies

Abstract The PCMBase R package is a powerful computational tool that enables efficient calculations of likelihoods for a wide range of phylogenetic Gaussian models. Taking advantage of it, we redesigned the R package mvSLOUCH. Here, we demonstrate how the new version of the package can be used to thoroughly examine the evolution and adaptation of traits in a large dataset of 1252 vascular plants through the use of multivariate Ornstein–Uhlenbeck processes. The results of our analysis demonstrate the ability of the modelling framework to distinguish between various alternative hypotheses regarding the evolution of functional traits in angiosperms.

Body mass, diet, foraging strata and foraging behaviour predict variation in the propensity to join mixed‐species flocks across birds

Mixed‐species flocking is common in birds and is associated with increased foraging efficiency and reduced predation pressure. The propensity to join mixed‐species flocks within a community varies across species, possibly reflecting species‐specific needs for increased protection from predators, opportunities to increase foraging efficiency or both. Previous studies on predictors of mixed‐species flocking propensity across species have typically focused on local communities, with or without accounting for phylogenetic relatedness. Here, we examined predictors of mixed‐species flocking propensity across the published literature using a phylogenetic linear mixed model framework. We obtained 724 mixed‐species flocking propensity estimates from 31 different study sites located primarily in South America and East Asia. Flocking propensity was associated with diet, foraging strata and foraging behaviour, and was inversely related to body mass. Controlling for phylogeny and body mass, flocking propensity was significantly lower for nectarivores than for species with other diets, lower for species foraging on the ground than in higher strata, and lower for sallying species than for bark and near‐perch foragers. These results suggest that variation across species in the tendency to join mixed‐species flocks can be predicted from knowledge about various species‐specific ecological traits in birds.

Polymorphism-Aware Models in RevBayes: Species Trees, Disentangling Balancing Selection, and GC-Biased Gene Conversion.

The role of balancing selection is a long-standing evolutionary puzzle. Balancing selection is a crucial evolutionary process that maintains genetic variation (polymorphism) over extended periods of time; however, detecting it poses a significant challenge. Building upon the Polymorphism-aware phylogenetic Models (PoMos) framework rooted in the Moran model, we introduce a PoMoBalance model. This novel approach is designed to disentangle the interplay of mutation, genetic drift, and directional selection (GC-biased gene conversion), along with the previously unexplored balancing selection pressures on ultra-long timescales comparable with species divergence times by analyzing multi-individual genomic and phylogenetic divergence data. Implemented in the open-source RevBayes Bayesian framework, PoMoBalance offers a versatile tool for inferring phylogenetic trees as well as quantifying various selective pressures. The novel aspect of our approach in studying balancing selection lies in polymorphism-aware phylogenetic models' ability to account for ancestral polymorphisms and incorporate parameters that measure frequency-dependent selection, allowing us to determine the strength of the effect and exact frequencies under selection. We implemented validation tests and assessed the model on the data simulated with SLiM and a custom Moran model simulator. Real sequence analysis of Drosophila populations reveals insights into the evolutionary dynamics of regions subject to frequency-dependent balancing selection, particularly in the context of sex-limited color dimorphism in Drosophila erecta.

Exploring Paleogene Tibet's warm temperate environments through target enrichment and phylogenetic niche modelling of Himalayan spiny frogs (Paini, Dicroglossidae).

The Cenozoic topographic development of the Himalaya-Tibet orogen (HTO) substantially affected the paleoenvironment and biodiversity patterns of High Asia. However, concepts on the evolution and paleoenvironmental history of the HTO differ massively in timing, elevational increase and sequence of surface uplift of the different elements of the orogen. Using target enrichment of a large set of transcriptome-derived markers, ancestral range estimation and paleoclimatic niche modelling, we assess a recently proposed concept of a warm temperate paleo-Tibet in Asian spiny frogs of the tribe Paini and reconstruct their historical biogeography. That concept was previously developed in invertebrates. Because of their early evolutionary origin, low dispersal capacity, high degree of local endemism, and strict dependence on temperature and humidity, the cladogenesis of spiny frogs may echo the evolution of the HTO paleoenvironment. We show that diversification of main lineages occurred during the early to Mid-Miocene, while the evolution of alpine taxa started during the late Miocene/early Pliocene. Our distribution and niche modelling results indicate range shifts and niche stability that may explain the modern disjunct distributions of spiny frogs. They probably maintained their (sub)tropical or (warm)temperate preferences and moved out of the ancestral paleo-Tibetan area into the Himalaya as the climate shifted, as opposed to adapting insitu. Based on ancestral range estimation, we assume the existence of low-elevation, climatically suitable corridors across paleo-Tibet during the Miocene along the Kunlun, Qiangtang and/or Gangdese Shan. Our results contribute to a deeper understanding of the mechanisms and processes of faunal evolution in the HTO.

Pathogen genomics and One Health: A scoping review of current practices in zoonotic disease research

ObjectivesWhole-genome sequencing has revolutionized the field of infectious disease surveillance, enabling near real-time detection of pathogens and tracking how infections may spread. Our study aimed to characterize genomic applications to cross-domain zoonotic pathogen transmission at the human-animal and/or human-environment interfaces. MethodsWe performed a scoping review of studies that have applied genomic epidemiology to zoonotic disease transmission across One Health domains (human, animal, and environment). We identified 114 records published between 2005 and 2022 which reported multi-domain genomic data of zoonotic pathogens integrated into phylogenetic models. ResultsMost studies investigated bacterial pathogens, highlighting key knowledge gaps for other zoonotic agents, particularly arboviruses. Sampling and sequencing efforts varied greatly across domains: the median number and range of pathogen genomes analyzed were highest for humans (23; 1–29,586) and lowest for the environment domain (13; 1–956). Genomics was used to track zoonotic disease outbreaks and cross-domain transmission, to improve pathogen surveillance, and to disentangle evolutionary dynamics driving lineage diversification and virulence. ConclusionsOur study highlights current practices and knowledge gaps to guide future study designs and genomic applications to multi-domain and cross-species transmission of zoonoses, with the potential to identify key infection sources and inform interventions for local and global health security.

Prevalence and composition of haemosporidians in an avian community from a World Heritage area: Associations with host foraging strata and forest regeneration

Forest regeneration is becoming a powerful tool to combat land conversion which covers 30 % of the Neotropical territory. However, little is known about the effect of forest regeneration on vector-borne diseases. Here, we describe the haemosporidian lineage composition across a successional gradient within an Atlantic Forest bird community. We test whether forest successional stages, in addition to host life history traits affect haemosporidian infection probability. We sampled birds at 16 sampling units with different successional stages between 2017 and 2018 within a forest remnant located in Antonina, Paraná, Brazil. We captured bird individuals using mist-nets, identified them to the species level, and collected blood samples to detect and identify Plasmodium and Haemoproteus lineages based on molecular analysis. We used a Bayesian phylogenetic linear model with a Bernoulli distribution to test whether the haemosporidian infection probability is affected by nest type, foraging stratum, and forest successional stage. We captured 322 bird individuals belonging to 52 species and 21 families. We found 31 parasite lineages and an overall haemosporidian prevalence of 23.9 %, with most infections being caused by Plasmodium (21.7 % of prevalence). The Plasmodium probability of infection was associated with forest successional stage and bird foraging stratum. Birds from the secondary forest in an intermediate stage of succession are more likely to be infected by the parasites than birds from the primary forests (β = 1.21, 95 % CI = 0.11 – 2.43), birds from upper strata exhibit a lower probability of infection than birds from lower foraging strata (β = -1.81, 95 % CI = -3.80 – -0.08). Nest type did not affect the Plasmodium probability of infection. Our results highlight the relevance of forest succession on haemosporidian infection dynamics, which is particularly relevant in a world where natural regeneration is the main tool used in forest restoration.

Deterministic processes dominate microbial assembly mechanisms in the gut microbiota of cold-water fish between summer and winter.

Exploring the effects of seasonal variation on the gut microbiota of cold-water fish plays an important role in understanding the relationship between seasonal variation and cold-water fish. Gut samples of cold-water fish and environmental samples were collected during summer and winter from the lower reaches of the Yalong River. The results of the 16S rRNA sequencing showed that significant differences were identified in the composition and diversity of gut bacteria of cold-water fish. Co-occurrence network complexity of the gut bacteria of cold-water fish was higher in summer compared to winter (Sum: nodes: 256; edges: 20,450; Win: nodes: 580; edges: 16,725). Furthermore, from summer to winter, the contribution of sediment bacteria (Sum: 5.3%; Win: 23.7%) decreased in the gut bacteria of cold-water fish, while the contribution of water bacteria (Sum: 0%; Win: 27.7%) increased. The normalized stochastic ratio (NST) and infer community assembly mechanisms by phylogenetic bin-based null model analysis (iCAMP) showed that deterministic processes played a more important role than stochastic processes in the microbial assembly mechanism of gut bacteria of cold-water fish. From summer to winter, the contribution of deterministic processes to gut bacteria community assembly mechanisms decreased, while the contribution of stochastic processes increased. Overall, these results demonstrated that seasonal variation influenced the gut bacteria of cold-water fish and served as a potential reference for future research to understand the adaptation of fish to varying environments.

Measuring complexity for hierarchical models using effective degrees of freedom

AbstractHierarchical models can express ecological dynamics using a combination of fixed and random effects, and measurement of their complexity (effective degrees of freedom, EDF) requires estimating how much random effects are shrunk toward a shared mean. Estimating EDF is helpful to (1) penalize complexity during model selection and (2) to improve understanding of model behavior. I applied the conditional Akaike Information Criterion (cAIC) to estimate EDF from the finite‐difference approximation to the gradient of model predictions with respect to each datum. I confirmed that this has similar behavior to widely used Bayesian criteria, and I illustrated ecological applications using three case studies. The first compared model parsimony with or without time‐varying parameters when predicting density‐dependent survival, where cAIC favors time‐varying demographic parameters more than conventional Akaike Information Criterion. The second estimates EDF in a phylogenetic structural equation model, and identifies a larger EDF when predicting longevity than mortality rates in fishes. The third compares EDF for a species distribution model fitted for 20 bird species and identifies those species requiring more model complexity. These highlight the ecological and statistical insight from comparing EDF among experimental units, models, and data partitions, using an approach that can be broadly adopted for nonlinear ecological models.

Selection Shapes the Genomic Landscape of Introgressed Ancestry in a Pair of Sympatric Sea Urchin Species.

A growing number of recent studies have demonstrated that introgression is common across the tree of life. However, we still have a limited understanding of the fate and fitness consequence of introgressed variation at the whole-genome scale across diverse taxonomic groups. Here, we implemented a phylogenetic hidden Markov model to identify and characterize introgressed genomic regions in a pair of well-diverged, nonsister sea urchin species: Strongylocentrotus pallidus and Strongylocentrotus droebachiensis. Despite the old age of introgression, a sizable fraction of the genome (1% to 5%) exhibited introgressed ancestry, including numerous genes showing signals of historical positive selection that may represent cases of adaptive introgression. One striking result was the overrepresentation of hyalin genes in the identified introgressed regions despite observing considerable overall evidence of selection against introgression. There was a negative correlation between introgression and chromosome gene density, and two chromosomes were observed with considerably reduced introgression. Relative to the nonintrogressed genome-wide background, introgressed regions had significantly reduced nucleotide divergence (dXY) and overlapped fewer protein-coding genes, coding bases, and genes with a history of positive selection. Additionally, genes residing within introgressed regions showed slower rates of evolution (dN, dS, dN/dS) than random samples of genes without introgressed ancestry. Overall, our findings are consistent with widespread selection against introgressed ancestry across the genome and suggest that slowly evolving, low-divergence genomic regions are more likely to move between species and avoid negative selection following hybridization and introgression.

Evolutionary patterns of variations in chromosome counts and genome sizes show positive correlations with taxonomic diversity in tropical gingers.

Cytogenetic traits such as an organism's chromosome number and genome size are taxonomically critical as they are instrumental in defining angiosperm diversity. Variations in these traits can be traced to evolutionary processes such as polyploidization, although geographic variations across cytogenetic traits remain underexplored. In the pantropical monocot family Zingiberaceae (~1500 species), cytogenetic traits have been well documented; however, the role of these traits in shaping taxonomic diversity and biogeographic patterns of gingers is not known. A time-calibrated Bayesian phylogenetic tree was constructed for 290 taxa covering three of the four subfamilies in Zingiberaceae. We tested models of chromosome number and genome size evolution within the family and whether lineage age, taxonomic diversity, and distributional range explain the variations in the cytogenetic traits. Tests were carried out at two taxonomic ranks: within Zingiberaceae and within genus Hedychium using correlations, generalized linear models and phylogenetic least square models. The most frequent changes in chromosome number within Zingiberaceae were noted to be demi-polyploidization and polyploidization (~57% of the time), followed by ascending dysploidy (~27%). The subfamily Zingiberoideae showed descending dysploidy at its base, while Alpinioideae showed polyploidization at its internal nodes. Although chromosome counts and genome sizes did not corroborate with each other, suggesting that they are not equivalent; higher chromosome number variations and higher genome size variations were associated with higher taxonomic diversity and wider biogeographic distribution. Within Zingiberaceae, multiple incidences of polyploidization were discovered, and cytogenetic events appear to have reduced the genome sizes and increased taxonomic diversity, distributional ranges and invasiveness.

Minimal and hybrid hydrogenases are active from archaea

Microbial hydrogen (H2) cycling underpins the diversity and functionality of diverse anoxic ecosystems. Among the three evolutionarily distinct hydrogenase superfamilies responsible, [FeFe] hydrogenases were thought to be restricted to bacteria and eukaryotes. Here, we show that anaerobic archaea encode diverse, active, and ancient lineages of [FeFe] hydrogenases through combining analysis of existing and new genomes with extensive biochemical experiments. [FeFe] hydrogenases are encoded by genomes of nine archaeal phyla and expressed by H2-producing Asgard archaeon cultures. We report an ultraminimal hydrogenase in DPANN archaea that binds the catalytic H-cluster and produces H2. Moreover, we identify and characterize remarkable hybrid complexes formed through the fusion of [FeFe] and [NiFe] hydrogenases in ten other archaeal orders. Phylogenetic analysis and structural modeling suggest a deep evolutionary history of hybrid hydrogenases. These findings reveal new metabolic adaptations of archaea, streamlined H2 catalysts for biotechnological development, and a surprisingly intertwined evolutionary history between the two major H2-metabolizing enzymes.

Bacterial sepsis triggers stronger transcriptomic immune responses in larger primates.

Empirical data relating body mass to immune defence against infections remain limited. Although the metabolic theory of ecology predicts that larger organisms would have weaker immune responses, recent studies have suggested that the opposite may be true. These discoveries have led to the safety factor hypothesis, which proposes that larger organisms have evolved stronger immune defences because they carry greater risks of exposure to pathogens and parasites. In this study, we simulated sepsis by exposing blood from nine primate species to a bacterial lipopolysaccharide (LPS), measured the relative expression of immune and other genes using RNAseq, and fitted phylogenetic models to determine how gene expression was related to body mass. In contrast to non-immune-annotated genes, we discovered hypermetric scaling in the LPS-induced expression of innate immune genes, such that large primates had a disproportionately greater increase in gene expression of immune genes compared to small primates. Hypermetric immune gene expression appears to support the safety factor hypothesis, though this pattern may represent a balanced evolutionary mechanism to compensate for lower per-transcript immunological effectiveness. This study contributes to the growing body of immune allometry research, highlighting its importance in understanding the complex interplay between body size and immunity over evolutionary timescales.

Non-structural carbohydrate concentrations in tree organs vary across biomes and leaf habits, but are independent of the fast-slow plant economic spectrum.

Carbohydrate reserves play a vital role in plant survival during periods of negative carbon balance. Under a carbon-limited scenario, we expect a trade-offs between carbon allocation to growth, reserves, and defense. A resulting hypothesis is that carbon allocation to reserves exhibits a coordinated variation with functional traits associated with the 'fast-slow' plant economics spectrum. We tested the relationship between non-structural carbohydrates (NSC) of tree organs and functional traits using 61 angiosperm tree species from temperate and tropical forests with phylogenetic hierarchical Bayesian models. Our results provide evidence that NSC concentrations in stems and branches are decoupled from plant functional traits. while those in roots are weakly coupled with plant functional traits. In contrast, we found that variation between NSC concentrations in leaves and the fast-slow trait spectrum was coordinated, as species with higher leaf NSC had trait values associated with resource conservative species, such as lower SLA, leaf N, and leaf P. We also detected a small effect of leaf habit on the variation of NSC concentrations in branches and roots. Efforts to predict the response of ecosystems to global change will need to integrate a suite of plant traits, such as NSC concentrations in woody organs, that are independent of the 'fast-slow' plant economics spectrum and that capture how species respond to a broad range of global change drivers.

Temporal changes in spatial scale and autocorrelation structure of forest openings based on taxonomic, phylogenetic, and functional turnover

ContextDuring European settlement, Illinois grasslands were converted for agricultural purposes. Remaining natural areas in southern Illinois include natural xeric forest openings, with communities representative of remnant grasslands and adjacent hardwood forests. Previous research in these openings shows plant communities are driven by edaphic conditions. ObjectivesThe first objective aimed to characterize spatial scale and autocorrelation structure of these openings based on climatic, environmental, and diversity variables. The second objective was to predict taxonomic, phylogenetic, and functional turnover between 1988 and 2019, using climatic and environmental variables. MethodsSurveys were conducted to calculate taxonomic, phylogenetics and functional trait metrics and analyses of these dimensions of diversity. Randomization tests were used to assess phylogenetic and functional clustering and over-dispersion at each site. Spatially-explicit climatic and environmental variables were included from earlier surveys and data repositories. Global Moran's I and spatial autocorrelograms were used to assess spatial structure of climatic, environmental, and diversity variables and generalized dissimilarity modeling was used to characterize taxonomic, phylogenetic, and functional turnover based on environmental variables. ResultsSoil depth was the only environmental variable which exhibited significant global spatial autocorrelation. Overall, sandstone sites were phylogenetically over-dispersed while loess sites were phylogenetically clustered. Climate variables and diversity metrics exhibited significant spatial structure during surveys. Generalized dissimilarity models showed that geographic distance between openings was the most influential driver of turnover across surveys. ConclusionsPrevious glacial events explained the spatial structure of soil depth across sites, due to Quaternary loess deposition in loess sites. High diversity values were clustered in the southeastern portions of the study area. Functional generalized dissimilarity models best predicted turnover in these openings compared to taxonomic and phylogenetic models.

Quantification of microbial community assembly processes during degradation on diverse plastispheres based on physicochemical characters and phylogenetic bin-based null model analysis

To understand the differences in degradation processes depending on the chemical properties of polymers, it is necessary to both quantify the microbiome composition and evaluate the process of microbial turnover (i.e., community assembly processes) in a variety of polymer materials. In this study, using a phylogenetic bin-based null model analysis (i.e., iCAMP), we evaluated community assembly processes from original estuary water to 37 types of polymers, which provide overwhelmingly diverse niches for microbes, in 14-day incubation experiments. First, we evaluated the polymer properties related to degradation rates. Polymers with higher adipic acid (AdA) monomer exhibited higher motility, hydrophilicity, and degradation rates, whereas those with higher aromatic monomer exhibited the opposite trends. Second, microbiome composition analysis was performed, and the microbiomes were significantly changed by the AdA or aromatic content. This was consistent with the polymer properties, suggesting that polymer motility and hydrophilicity attributable to the first-order structure modify the accessibility of the enzyme to the reaction site and hence the degradation rate, resulting in differences in microbiome community composition. Finally, we determined community assembly processes from estuary water to plastics using a phylogenetic bin-based null model analysis. The importance of heterogeneous selection was higher in mobile, hydrophilic, and fast-degrading polymers, while that of homogeneous selection was lower. This suggests that the environmental difference between before and after incubation becomes significant under rapid degradation, which select microbes adapted to biofilm environments. In addition, the more stochastic turnover prevailed, the more variation in the communities (i.e., β-diversity) increased. This suggests that turnover processes not dictated by the environment lead to instability in community compositions.

Evolutionary analyses of intrinsically disordered regions reveal widespread signals of conservation.

Intrinsically disordered regions (IDRs) are segments of proteins without stable three-dimensional structures. As this flexibility allows them to interact with diverse binding partners, IDRs play key roles in cell signaling and gene expression. Despite the prevalence and importance of IDRs in eukaryotic proteomes and various biological processes, associating them with specific molecular functions remains a significant challenge due to their high rates of sequence evolution. However, by comparing the observed values of various IDR-associated properties against those generated under a simulated model of evolution, a recent study found most IDRs across the entire yeast proteome contain conserved features. Furthermore, it showed clusters of IDRs with common "evolutionary signatures," i.e. patterns of conserved features, were associated with specific biological functions. To determine if similar patterns of conservation are found in the IDRs of other systems, in this work we applied a series of phylogenetic models to over 7,500 orthologous IDRs identified in the Drosophila genome to dissect the forces driving their evolution. By comparing models of constrained and unconstrained continuous trait evolution using the Brownian motion and Ornstein-Uhlenbeck models, respectively, we identified signals of widespread constraint, indicating conservation of distributed features is mechanism of IDR evolution common to multiple biological systems. In contrast to the previous study in yeast, however, we observed limited evidence of IDR clusters with specific biological functions, which suggests a more complex relationship between evolutionary constraints and function in the IDRs of multicellular organisms.

Evolutionary characteristics of the mitochondrial NADH dehydrogenase subunit 6 gene in some populations of four sympatric Mustela species (Mustelidae, Mammalia) from central Europe.

Selection on or reticulate evolution of mtDNA is documented in various mammalian taxa and could lead to misleading phylogenetic conclusions if not recognized. We sequenced the MT-ND6 gene of four sympatric Mustelid species of the genus Mustela from some central European populations. We hypothesised positive selection on MT-ND6, given its functional importance and the different body sizes and life histories of the species, even though climatic differences may be unimportant for adaptation in sympatry. MT-ND6 genes were sequenced in 187 sympatric specimens of weasels, Mustela nivalis, stoats, M. erminea, polecats, M. putorius, and steppe polecats, M. eversmannii, from eastern Austria and of fourteen allopatric polecats from eastern-central Germany. Median joining networks, neighbour joining and maximum likelihood analyses as well as Bayesian inference grouped all species according to earlier published phylogenetic models. However, polecats and steppe polecats, two very closely related species, shared the same two haplotypes. We found only negative selection within the Mustela sequences, including 131 downloaded ones covering thirteen species. Positive selection was observed on three MT-ND6 codons of other mustelid genera retrieved from GenBank. Negative selection for MT-ND6 within the genus Mustela suggests absence of both environmental and species-specific effects of cellular energy metabolism despite large species-specific differences in body size. The presently found shared polymorphism in European polecats and steppe polecats may result from ancestral polymorphism before speciation and historical or recent introgressive hybridization; it may indicate mtDNA capture of steppe polecats by M. putorius in Europe.

Binding Evolution of the Dengue Virus Envelope Against DC-SIGN: A Combined Approach of Phylogenetics and Molecular Dynamics Analyses Over 30 Years of Dengue Virus in Brazil

The Red Queen Hypothesis (RQH), derived from Lewis Carroll's “Through the Looking-Glass”, postulates that organisms must continually adapt in response to each other to maintain relative fitness. Within the context of host-pathogen interactions, the RQH implies an evolutionary arms race, wherein viruses evolve to exploit hosts and hosts evolve to resist viral invasion. This study delves into the dynamics of the RQH in the context of virus-cell interactions, specifically focusing on virus receptors and cell receptors. We observed multiple virus-host systems and noted patterns of co-evolution. As viruses evolved receptor-binding proteins to effectively engage with cell receptors, cells countered by altering their receptor genes. This ongoing mutual adaptation cycle has influenced the molecular intricacies of receptor-ligand interactions. Our data supports the RQH as a driving force behind the diversification and specialization of both viral and host cell receptors. Understanding this co-evolutionary dance offers insights into the unpredictability of emerging viral diseases and potential therapeutic interventions. Future research is crucial to dissect the nuanced molecular changes and the broader ecological consequences of this ever-evolving battle. Here, we combine phylogenetic inferences, structural modeling, and molecular dynamics analyses to describe the epidemiological characteristics of major Brazilian DENV strains that circulated from 1990 to 2022 from a combined perspective, thus providing us with a more detailed picture on the dynamics of such interactions over time.