Phylogenomic analysis of Stylops reveals the evolutionary history of a Holarctic Strepsiptera radiation parasitizing wild bees

Holarctic Stylops is the largest genus of the enigmatic insect order Strepsiptera, twisted winged parasites. Members of Stylops are obligate endoparasites of Andrena mining bees and exhibit extreme sexual dimorphism typical of Strepsiptera. So far, molecular studies on Stylops have focused on questions on species delimitation. Here, we utilize the power of whole genome sequencing to infer the phylogeny of this morphologically challenging genus from thousands of loci. We use a species tree method, concatenated maximum likelihood analysis and Bayesian analysis with a relaxed clock model to reconstruct the phylogeny of 46 Stylops species, estimate divergence times, evaluate topological consistency across methods and infer the root position. Furthermore, the biogeographical history and coevolutionary patterns with host species are assessed. All methods recovered a well resolved topology with close to all nodes maximally supported and only a handful of minor topological variations. Based on the result, we find that included species can be divided into 12 species groups, seven of them including only Palaearctic species, three Nearctic and two were geographically mixed. We find a strongly supported root position between a clade formed by the spreta, thwaitesi and gwynanae species groups and the remaining species and that the sister group of Stylops is Eurystylops or Eurystylops + Kinzelbachus. Our results indicate that Stylops originated in the Western Palaearctic or Western Palaearctic and Nearctic in the early Neogene or late Paleogene, with four independent dispersal events to the Nearctic. Cophylogenetic analyses indicate that the diversification of Stylops has been shaped by both significant coevolution with the mining bee hosts and host-shifting. The well resolved and strongly supported phylogeny will provide a valuable phylogenetic basis for further studies into the fascinating world of Strepsipterans.

The role of rhizosphere phages in soil health

While the One Health framework has emphasized the importance of soil microbiomes for plant and human health, one of the most diverse and abundant groups-bacterial viruses, i.e. phages-has been mostly neglected. This perspective reviews the significance of phages for plant health in rhizosphere and explores their ecological and evolutionary impacts on soil ecosystems. We first summarize our current understanding of the diversity and ecological roles of phages in soil microbiomes in terms of nutrient cycling, top-down density regulation and pathogen suppression. We then consider how phages drive bacterial evolution in soils by promoting horizontal gene transfer, encoding auxiliary metabolic genes that increase host bacterial fitness and selecting for phage-resistant mutants with altered ecology due to trade-offs with pathogen competitiveness and virulence. Finally, we consider challenges and avenues for phage research in soil ecosystems and how to elucidate the significance of phages for microbial ecology and evolution and soil ecosystem functioning in the future. We conclude that similar to bacteria, phages likely play important roles in connecting different One Health compartments, affecting microbiome diversity and functions in soils. From the applied perspective, phages could offer novel approaches to modulate and optimize microbial and microbe-plant interactions to enhance soil health.

A positive selection at binding site 501 in the B.1 lineage might have triggered the highly infectious sub-lineages of SARS-CoV-2

The descendants of the B lineage are the most predominant variants among the SARS-CoV-2 virus due to the incorporation of new mutations augmenting the infectivity of the virus. There is a substantial increase in the transition transversion bias, nucleotide diversity and purifying selection on the spike protein in the descendants of the B lineage of the SARS-CoV-2 virus on a temporal scale. A strong bias for C-to-U substitutions is found in the genes encoding spike protein in this lineage. The positive selection has operated on the spike gene of B lineages and its sub-lineages. The B.1 lineage has undergone positive selection on site 501 of the receptor binding domain ultimately reflected in a key substitution N501Y in its three descendant lineages namely B.1.1.7, B.1.351 and P.1. The intensity of purifying selection on the multiple sites of the spike gene has increased substantially in the sub-lineages of B.1 in a timescale. The binding site 501 on the spike protein in B lineage is found to coevolve with other amino acid sites. This study sheds light on the evolutionary trajectory of the B lineage into highly infectious descendants in the recent past under the influence of positive and purifying selection exerted by natural immunity and vaccination of the host.

Autorepression of yeast Hsp70 cochaperones by intramolecular interactions involving their J-domains

The 70 kDa heat shock protein (Hsp70) chaperones control protein homeostasis in all ATP-containing cellular compartments. J-domain proteins (JDPs) coevolved with Hsp70s to trigger ATP hydrolysis and catalytically upload various substrate polypeptides in need to be structurally modified by the chaperone. Here, we measured the protein disaggregation and refolding activities of the main yeast cytosolic Hsp70, Ssa1, in the presence of its most abundant JDPs, Sis1 and Ydj1, and two swap mutants, in which the J-domains have been interchanged. The observed differences by which the four constructs differently cooperate with Ssa1 and cooperate with each other, as well as their observed intrinsic ability to bind misfolded substrates and trigger Ssa1's ATPase, indicate the presence of yet uncharacterized intramolecular dynamic interactions between the J-domains and the remaining C-terminal segments of these proteins. Taken together, the data suggest an autoregulatory role to these intramolecular interactions within both type A and B JDPs, which might have evolved to reduce energy-costly ATPase cycles by the Ssa1-4 chaperones that are the most abundant Hsp70s in the yeast cytosol.

Metamorphic proteins and how to find them

In the last two decades, our existing notion that most foldable proteins have a unique native state has been challenged by the discovery of metamorphic proteins, which reversibly interconvert between multiple, sometimes highly dissimilar, native states. As the number of known metamorphic proteins increases, several computational and experimental strategies have emerged for gaining insights about their refolding processes and identifying unknown metamorphic proteins amongst the known proteome. In this review, we describe the current advances in biophysically and functionally ascertaining the structural interconversions of metamorphic proteins and how coevolution can be harnessed to identify novel metamorphic proteins from sequence information. We also discuss the challenges and ongoing efforts in using artificial intelligence-based protein structure prediction methods to discover metamorphic proteins and predict their corresponding three-dimensional structures.

Personality and cognition: shoal size discrimination performance is related to boldness and sociability among ten freshwater fish species

Several studies have reported that animals' personalities are often correlated with individual differences in cognition. Here, we tested whether personality is related to cognition across species, focusing on 10 freshwater fishes and a task relevant for fitness, the ability to discriminate shoal size. Bolder species exhibited more 'shuttle' behavior for information sampling during shoal selection and showed high performance (HP) in the numerical discrimination than shyer species, i.e., low performance (LP) species. Species at both the high and low ends of sociability showed LP, possibly due to loosened selection pressure because of either no need to perform shoal size discrimination tasks frequently in nature for very high sociability species or decreased willingness and motivation to join and stay within shoals for very low sociability species. Notably, the numerical discrimination was sensitive to the numerical contrast ratio in LP species but not in HP species, suggesting that the numerical system used for size discrimination also varied between species. Overall, we demonstrated the interspecies relationship between personality and shoal size discrimination across fish species, suggesting an evolutionary link between numerical abilities and behavior.

A plethora of rodents: Rattlesnake predators generate unanticipated patterns of venom resistance in a grassland ecosystem

Predation has the potential to impart strong selective pressures on organisms within their environments, resulting in adaptive changes in prey that minimize risk of predation. Pressures from venomous snakes present an exceptional challenge to prey, as venom represents a unique chemical arsenal evolutionarily tailored to incapacitate prey. In response, venom resistance has been detected in various snake prey species, and to varying degrees. This study analyzes venom resistance in an eastern Colorado grassland habitat, where the Prairie Rattlesnake (Crotalus viridis) and Desert Massasauga Rattlesnake (Sistrurus tergeminus edwardsii) co-occur with a suite of grassland rodents. We test for venom resistance across rodent and snake pairings using two geographically distant field sites to determine the role of 1) predation pressure and trophic ecology, and 2) sympatric and allopatric patterns of venom resistance. Resistance was measured using serum-based metalloproteinase inhibition assays to determine potential inhibition of proteolytic activity, augmented by median lethal dose (LD50) assays on rodent species to assess toxicity of crude venoms. Resistance is present in several rodent species, with strong resistance present in populations of Eastern Woodrat (Neotoma floridana), Ord's Kangaroo Rat (Dipodomys ordii), and Northern Grasshopper Mouse (Onychomys leucogaster). Resistance is less developed in other species, including the House Mouse (Mus musculus) and Plains Pocket Mouse (Perognathus flavescens). An unexpected differential is present, where Lincoln County Kangaroo Rats are highly resistant to venom of co-occurring Prairie Rattlesnakes yet are sensitive to an allopatric population of Prairie Rattlesnakes in Weld County. Lincoln Co. Northern Grasshopper Mice also demonstrate extremely elevated resistance to Weld Co. Prairie Rattlesnake venoms, and they may possess resistance mechanisms for myotoxin a, an abundant component of Weld Co. C. v viridis venoms. This study illustrates the complexity of venom resistance in biological communities that can exist when incorporating multiple species interactions. Future studies aimed at characterizing resistance mechanisms at the molecular level will provide a more detailed physiological context for understanding mechanisms by which resistance to venoms occurs.

Separation of evolutionary timescales in coevolving species

Many coevolutionary processes, including host-parasite and host-symbiont interactions, involve one species or trait which evolves much faster than the other. Whether or not a coevolutionary trajectory converges depends on the relative rates of evolutionary change in the two species, and so current adaptive dynamics approaches generally either determine convergence stability by considering arbitrary (often comparable) rates of evolutionary change or else rely on necessary or sufficient conditions for convergence stability. We propose a method for determining convergence stability in the case where one species is expected to evolve much faster than the other. This requires a second separation of timescales, which assumes that the faster evolving species will reach its evolutionary equilibrium (if one exists) before a new mutation arises in the more slowly evolving species. This method, which is likely to be a reasonable approximation for many coevolving species, both provides straightforward conditions for convergence stability and is less computationally expensive than traditional analysis of coevolution models, as it reduces the trait space from a two-dimensional plane to a one-dimensional manifold. In this paper, we present the theory underlying this new separation of timescales and provide examples of how it could be used to determine coevolutionary outcomes from models.

Evaluating the role of rhizosphere microbial home-field advantage in Betula luminifera adaptation to antimony mining areas

It has been established that the coevolution of plants and the rhizosphere microbiome in response to abiotic stress can result in the recruitment of specific functional microbiomes. However, the potential of inoculated rhizosphere microbiomes to enhance plant fitness and the inheritance of adaptive traits in subsequent generations remains unclear. In this study, cross-inoculation trials were conducted using seeds, rhizosphere microbiome, and in situ soil collected from areas of Betula luminifera grown in both antimony mining and control sites. Compared to the control site, plants originating from mining areas exhibited stronger adaptive traits, specifically manifested as significant increases in hundred-seed weight, specific surface area, and germination rate, as well as markedly enhanced seedling survival rate and biomass. Inoculation with mining microbiomes could enhance the fitness of plants in mining sites through a "home-field advantage" while also improving the fitness of plants originating from control sites. During the initial phase of seedling development, bacteria play a crucial role in promoting growth, primarily due to their mechanisms of metal resistance and nutrient cycling. This study provided evidence that the outcomes of long-term coevolution between plants and the rhizosphere microbiome in mining areas can be passed on to future generations. Moreover, it has been demonstrated that transgenerational inheritance and rhizosphere microbiome inoculation are important factors in improving the adaptability of plants in mining areas. The findings have important implications for vegetation restoration and ecological environment improvement in mining areas.

Disentangling the mechanisms underlying phylosymbiosis in mammals

Mammalian gut microbial communities are frequently found to be host-specific-microbial community compositions are more similar within than between host species-and some individual microbial taxa consistently associate with a single or small set of host species. The ecoevolutionary dynamics that result in this pattern of phylosymbiosis or host specificity have been proposed, but robust tests of the mechanisms driving these relationships are lacking. In this issue of Molecular Ecology, Mazel et al. (2023) combine large amplicon sequencing data sets with bacterial phenotypic traits to test whether microbial dispersal patterns contribute to the host specificity of the gut microbiome. They find that both transmission mode and oxygen tolerance are predictive of how specialized a microbe is. Horizontally transmitted, oxygen-tolerant microbes are more likely to be generalists, and vertically transmitted anaerobes are more likely to be limited to a few host species. This creative use of publicly available data provides a roadmap for testing hypotheses about the mechanisms underlying phylosymbiosis.

Structure-oriented engineering of nitrile hydratase: Reshaping of substrate access tunnel and binding pocket for efficient synthesis of cinnamamide

Cinnamamide and its derivatives are the most common and important building blocks widely present in natural products. Currently, nitrile hydratase (NHase, EC 4.2.1.84) has been widely used in large-scale industrial production of nicotinamide and acrylamide, while its catalytic activity is extremely low or inactive for bulky nitrile substrates such as cinnamonitrile. Therefore, beneficial variant βF37P/L48P/F51N were obtained from PtNHase of Pseudonocardia thermophila JCM3095 by reshaping of substrate access tunnel and binding pocket, which exhibited 14.88-fold improved catalytic efficiency compared to the wild-type PtNHase. Structure analysis, molecular dynamics simulations and dynamical cross-correlation matrix (DCCM) analysis revealed that the introduced mutations enlarged the substrate access tunnel and binding pocket, enhanced overall anti-correlated movements of enzymes, which would promote product release during the dynamic process of catalysis. In a hydration process, the complete conversion of 5 mM cinnamonitrile was achieved by βF37P/L48P/F51N in a 50 mL reaction, with cinnamamide yield of almost 100 % and productivity of 0.736 g L-1 h-1. The study demonstrates the co-evolution of substrate access tunnel and binding pocket is an effective strategy, and provides a valuable reference for future research. Furthermore, NHases have huge potential for catalyzing bulky nitriles to form corresponding amides in large-scale industrial production.

Coevolution-based computational approach to detect resistance mechanism of epidermal growth factor receptor

Tyrosine kinase epidermal growth factor receptor (EGFR) correlates the neoplastic cell metastasis, angiogenesis, neoplastic incursion, and apoptosis. Due to the involvement of EGFR in these biological processes, it becomes a most potent target for treating non-small cell lung cancer (NSCLC). The tyrosine kinase inhibitors (TKI) have endorsed high efficacy and anticipation to patients but unfortunately, within a year of treatment, drug targets develop resistance due to mutations. The present study detected the compensatory mutations in EGFR to know the evolutionary mechanism of drug resistance. The results of this study demonstrate that compensatory mutations enlarge the drug-binding pocket which may lead to the altered orientation of the ligand (gefitinib and erlotinib) causing drug resistance. This indicates that coevolutionary forces play a significant role in fine-tuning the structure of EGFR protein against the drugs. The analysis provides insight into the evolution-induced structural aspects of drug resistance changes in EGFR which in turn be useful in designing drugs with better efficacy.

A New Termitophilous Rove Beetle from Brazil Suggests an Undersampling Scenario and Cospeciation between Termites and their Termitophiles in the Neotropical region (Coleoptera: Staphylinidae: Aleocharinae)

Numerous staphylinid species can be found in termite nests. While the distributional range of termite hosts has been well studied, the same cannot be said for their associated fauna. There is a lack of comprehensive sampling even in areas where termite fauna has been well-sampled. In recent surveys in Southeastern Brazil, we collected a new species of termitophilous rove beetle from the tribe Termitohospitini (Aleocharinae): Termitosocius boraceiensis sp. nov., is herein described. The new species is the first termitophile found in Boracéia, SP, Brazil, even with decades of termite sampling in the place. Additionally, we discuss the diversity and description rate of neotropical Termitohospitini, and investigate the potential of cospeciation between termites and their termitophilous beetles on the Emerson-Kistner principle. Our findings suggest that the neotropical Termitohospitini may be more diverse than previously believed and that Termitosocius species are present in host species that are geographically separated, implying a cospeciation scenario. Additional sampling is required to fully comprehend the relationships between termitophiles and their termite hosts. The outcomes of this study provide new insights into the undersampling scenario and possibly cospeciation between termites and their termitophiles in the Neotropical region.

Genetic Covariation Between the Vertically Transmitted Endophyte Epichloë canadensis and Its Host Canada Wildrye

Symbiotic mutualisms are thought to be stabilized by correlations between the interacting genotypes which may be strengthened via vertical transmission and/or reduced genetic variability within each species. Vertical transmission, however, may weaken interactions over time as the endosymbionts would acquire mutations that could not be purged. Additionally, temporal variation in a conditional mutualism could create genetic variation and increased variation in the interaction outcome. In this study, we assessed genetic variation in both members of a symbiosis, the endosymbiotic fungal endophyte Epichloë canadensis and its grass host Canada wildrye (Elymus canadensis). Both species exhibited comparable levels of diversity, mostly within populations rather than between. There were significant differences between populations, although not in the same pattern for the two species, and the differences were not correlated with geographic distance for either species. Interindividual genetic distance matrices for the two species were significantly correlated, although all combinations of discriminant analysis of principle components (DAPC) defined multilocus genotype groups were found suggesting that strict genotype matching is not necessary. Variation in interaction outcome is common in grass/endophyte interactions, and our results suggest that the accumulation of mutations overtime combined with temporal variation in selection pressures increasing genetic variation in the symbiosis may be the cause.

DO PARASITIC LICE EXHIBIT ENDEMISM IN PARALLEL WITH THEIR AVIAN HOSTS? A COMPARISON ACROSS NORTHERN AMAZONIAN AREAS OF ENDEMISM

Areas of endemism are the smallest units in biogeography and can be defined as biologically unique areas comprising taxa with common geographic limits to their distributions. High beta diversity within Amazonia is often related to turnover among these areas. For decades, evolutionary biologists have tried to comprehend the mechanisms generating and maintaining the spatial structure and high diversity of free-living Amazonian organisms, particularly birds. However, few studies have tried to analyze these patterns among their parasites. Host and parasite associations involve shared history that may allow us to better understand the fine-scale evolutionary history of the host. Here we compare the coevolutionary patterns among 2 avian host species with distinct patterns of genetic structure in northern Amazonia, Dendrocincla fuliginosa (Aves: Dendrocolaptidae) and Dixiphia pipra (Aves: Pipridae), and their ectoparasitic lice (Insecta: Phthiraptera), Furnaricola sp. ex Dendrocincla fuliginosa, Myrsidea sp. ex Dixiphia pipra, and Tyranniphilopterus sp. ex Dixiphia pipra. We obtained sequences of the mitochondrial gene cytochrome oxidase subunit I from hosts and parasites collected on opposite banks of the Negro and Japurá rivers, which delimit 3 areas of endemism in northern Amazonia: Napo, Jau, and Guiana. Our results demonstrate that the Negro River is a geographical barrier for both Furnaricola sp. and its avian host, Dendrocincla fuliginosa. Phylogenies of both hosts, Dendrocincla fuliginosa, and the parasites, Furnaricola sp., show monophyletic clades on opposite margins of the river that are not sister taxa. These clades have a mean uncorrected p-distance of 17.8% for Furnaricola sp. and 6.0% for Dendrocincla fuliginosa. Thus, these parasite clades constitute distinct evolutionary lineages and may even be distinct species. In contrast, Dixiphia pipra has no population structure associated with either river. Accordingly, data from their lice Myrsidea sp. indicate weak support for different clades on opposite margins of the Negro River, whereas data from their lice Tyranniphilopterus sp. indicate weak structure across the Japurá. This study is a first step toward understanding the effects of biogeographic history on permanent ectoparasites and suggests that host biogeographic history is to some extent a determinant of the parasite's history. Furthermore, the parasite's evolutionary history is an additional source of information about their hosts' evolution in this highly diverse region of northern Amazonia.

Coevolution based immunoinformatics approach considering variability of epitopes to combat different strains: A case study using spike protein of SARS-CoV-2

In the recent past several vaccines were developed to combat the COVID-19 disease. Unfortunately, the protective efficacy of the current vaccines has been reduced due to the high mutation rate in SARS-CoV-2. Here, we successfully implemented a coevolution based immunoinformatics approach to design an epitope-based peptide vaccine considering variability in spike protein of SARS-CoV-2. The spike glycoprotein was investigated for B- and T-cell epitope prediction. Identified T-cell epitopes were mapped on previously reported coevolving amino acids in the spike protein to introduce mutation. The non-mutated and mutated vaccine components were constructed by selecting epitopes showing overlapping with the predicted B-cell epitopes and highest antigenicity. Selected epitopes were linked with the help of a linker to construct a single vaccine component. Non-mutated and mutated vaccine component sequences were modelled and validated. The in-silico expression level of the vaccine constructs (non-mutated and mutated) in E. coli K12 shows promising results. The molecular docking analysis of vaccine components with toll-like receptor 5 (TLR5) demonstrated strong binding affinity. The time series calculations including root mean square deviation (RMSD), radius of gyration (RGYR), and energy of the system over 100 ns trajectory obtained from all atom molecular dynamics simulation showed stability of the system. The combined coevolutionary and immunoinformatics approach used in this study will certainly help to design an effective peptide vaccine that may work against different strains of SARS-CoV-2. Moreover, the strategy used in this study can be implemented on other pathogens.

Facultative symbiont virulence determines horizontal transmission rate without host strain specificity

In facultative symbioses, only a fraction of hosts are associated with a symbiont. Understanding why specific host and symbiont strains are associated can inform us of the evolutionary forces affecting facultative symbioses. Possibilities include ongoing host-symbiont coevolution driven by reciprocal selection, or priority effects that are neutral in respect to the host-symbiont interaction. We hypothesized that ongoing host-symbiont coevolution would lead to higher fitness estimates for naturally co-occurring (native) host and symbiont combinations compared to nonnative combinations. We used the Dictyostelium discoideum - Paraburkholderia bonniea system to test this hypothesis. P. bonniea features a reduced genome size relative to another Paraburkholderia symbiont of D. discoideum, indicating a significant history of coevolution with its host. Facultative symbionts may experience continued genome reduction if coevolution is ongoing, or their genome size may have reached a stable state if the symbiosis has also stabilized. Our work demonstrates that ongoing coevolution is unlikely for D. discoideum and P. bonniea. The system instead represents a stable facultative symbiosis. Specifically associated host and symbiont strains in this system are the result of priority effects, and presently unassociated hosts are simply uncolonized. We find evidence for a virulence-transmission trade-off without host strain specificity, and identify candidate virulence factors in the genomes of P. bonniea strains that may contribute to variation in benevolence.

Uncovering the genomic basis of infection through co-genomic sequencing of hosts and parasites

Understanding the genomic basis of infectious disease is fundamental objective in coevolutionary theory with relevance to healthcare, agriculture, and epidemiology. Models of host-parasite coevolution often assume that infection requires specific combinations of host and parasite genotypes. Coevolving host and parasite loci are therefor expected to show associations that reflects an underlying infection/resistance allele matrix, yet little evidence for such genome-to-genome interactions has been observed among natural populations. We conducted a study to search for this genomic signature across 258 linked host (Daphnia magna) and parasite (Pasteuria ramosa) genomes. Our results show a clear signal of genomic association between multiple epistatically-interacting loci in the host genome, and a family of genes encoding for collagen-like protein in the parasite genome. These findings are supported by laboratory-based infection trials, which show strong correspondence between phenotype and genotype at the identified loci. Our study provides clear genomic evidence of antagonistic coevolution among wild populations.

CastNet: a systems-level sequence evolution simulator

BACKGROUND

Simulating DNA evolution has been done through coevolution-agnostic probabilistic frameworks for the past 3 decades. The most common implementation is by using the converse of the probabilistic approach used to infer phylogenies which, in the simplest form, simulates a single sequence at a time. However, biological systems are multi-genic, and gene products can affect each other's evolutionary paths through coevolution. These crucial evolutionary dynamics still remain to be simulated, and we believe that modelling them can lead to profound insights for comparative genomics.

RESULTS

Here we present CastNet, a genome evolution simulator that assumes each genome is a collection of genes with constantly evolving regulatory interactions in between them. The regulatory interactions produce a phenotype in the form of gene expression profiles, upon which fitness is calculated. A genetic algorithm is then used to evolve a population of such entities through a user-defined phylogeny. Importantly, the regulatory mutations are a response to sequence mutations, thus making a 1-1 relationship between the rate of evolution of sequences and of regulatory parameters. This is, to our knowledge, the first time the evolution of sequences and regulation have been explicitly linked in a simulation, despite there being a multitude of sequence evolution simulators, and a handful of models to simulate Gene Regulatory Network (GRN) evolution. In our test runs, we see a coevolutionary signal among genes that are active in the GRN, and neutral evolution in genes that are not included in the network, showing that selective pressures imposed on the regulatory output of the genes are reflected in their sequences.

CONCLUSION

We believe that CastNet represents a substantial step for developing new tools to study genome evolution, and more broadly, coevolutionary webs and complex evolving systems. This simulator also provides a new framework to study molecular evolution where sequence coevolution has a leading role.

Naïve Pine Terpene Response to the Mountain Pine Beetle (Dendroctonus ponderosae) through the Seasons

Insect herbivores must contend with constitutive and induced plant defenses. The mountain pine beetle, Dendroctonus ponderosae Hopkins (Coleoptera: Curculionidae, Scolytinae) has expanded its range east of the Rocky Mountains into the western boreal forest and is encountering evolutionarily naïve lodgepole pines (Pinus contorta) and jack pines (Pinus banksiana). Pinus contorta and P. banksiana in the expanded range have different constitutive and induced defenses in response to wounding and inoculation with fungal associates of D. ponderosae. In the historic range, previous studies have examined phloem terpene content prior to and just after D. ponderosae mass attack, but the terpene profile of attacked trees post-overwintering is unknown. We examined the response of mature P. contorta and P. banksiana trees to experimentally-applied mass attack by D. ponderosae and quantified phloem terpenes at three time points, pre-attack, post-attack (same season), and the following spring, post-overwintering. Phloem content of total terpenes as well as many individual terpenes increased after D. ponderosae attack but were only significantly higher than pre-attack levels at the post-overwintering time point in both P. contorta and P. banksiana. The absence of a significant increase in phloem terpenes in the month following attack in naïve pines is a potential cause for increased D. ponderosae offspring production reported in naïve P. contorta. Beetle attack density did not influence the phloem terpene profiles of either species and there was no significant interaction between attack density and sampling time on terpene content. High phloem terpenes in trees that are attacked at low densities could prime these trees for defense against attacks in the following season but it could also make these trees more apparent to early-foraging beetles and facilitate efficient mass attack at low D. ponderosae population densities in the expanded range.

Eco-evolutionary cyclic dominance among predators, prey, and parasites

Predator-prey interactions are one of ecology's central research themes, but with many interdisciplinary implications across the social and natural sciences. Here we consider an often-overlooked species in these interactions, namely parasites. We first show that a simple predator-prey-parasite model, inspired by the classical Lotka-Volterra equations, fails to produce a stable coexistence of all three species, thus failing to provide a biologically realistic outcome. To improve this, we introduce free space as a relevant eco-evolutionary component in a new mathematical model that uses a game-theoretical payoff matrix to describe a more realistic setup. We then show that the consideration of free space stabilizes the dynamics by means of cyclic dominance that emerges between the three species. We determine the parameter regions of coexistence as well as the types of bifurcations leading to it by means of analytical derivations as well as by means of numerical simulations. We conclude that the consideration of free space as a finite resource reveals the limits of biodiversity in predator-prey-parasite interactions, and it may also help us in the determination of factors that promote a healthy biota.

The coevolution of the spread of a disease and competing opinions in multiplex networks

The COVID-19 pandemic has resulted in a proliferation of conflicting opinions on physical distancing across various media platforms, which has had a significant impact on human behavior and the transmission dynamics of the disease. Inspired by this social phenomenon, we present a novel UAP-SIS model to study the interaction between conflicting opinions and epidemic spreading in multiplex networks, in which individual behavior is based on diverse opinions. We distinguish susceptibility and infectivity among individuals who are unaware, pro-physical distancing and anti-physical distancing, and we incorporate three kinds of mechanisms for generating individual awareness. The coupled dynamics are analyzed in terms of a microscopic Markov chain approach that encompasses the aforementioned elements. With this model, we derive the epidemic threshold which is related to the diffusion of competing opinions and their coupling configuration. Our findings demonstrate that the transmission of the disease is shaped in a significant manner by conflicting opinions, due to the complex interaction between such opinions and the disease itself. Furthermore, the implementation of awareness-generating mechanisms can help to mitigate the overall prevalence of the epidemic, and global awareness and self-awareness can be interchangeable in certain instances. To effectively curb the spread of epidemics, policymakers should take steps to regulate social media and promote physical distancing as the mainstream opinion.

To explore strange new worlds - The diversification in Tremella caloplacae was linked to the adaptive radiation of the Teloschistaceae

Lichenicolous fungi are a heterogeneous group of organisms that grow exclusively on lichens, forming obligate associations with them. It has often been assumed that cospeciation has occurred between lichens and lichenicolous fungi, but this has been seldom analysed from a macroevolutionary perspective. Many lichenicolous species are rare or are rarely observed, which results in frequent and large gaps in the knowledge of the diversity of many groups. This, in turn, hampers evolutionary studies that necessarily are based on a reasonable knowledge of this diversity. Tremella caloplacae is a heterobasidiomycete growing on various hosts from the lichen-forming family Teloschistaceae, and evidence suggests that it may represent a species complex. We combine an exhaustive sampling with molecular and ecological data to study species delimitation, cophylogenetic events and temporal concordance of this association. Tremella caloplacae is here shown to include at least six distinct host-specific lineages (=putative species). Host switch is the dominant and most plausible event influencing diversification and explaining the coupled evolutionary history in this system, although cospeciation cannot be discarded. Speciation in T. caloplacae would therefore have occurred coinciding with the rapid diversification - by an adaptive radiation starting in the late Cretaceous - of their hosts. New species in T. caloplacae would have developed as a result of specialization on diversifying lichen hosts that suddenly offered abundant new ecological niches to explore or adapt to.

Phage-bacterial evolutionary interactions: experimental models and complications

Phage treatment of bacterial infections has offered some hope even as the crisis of antimicrobial resistance continues to be on the rise. However, bacterial resistance to phage is another looming challenge capable of undermining the effectiveness of phage therapy. Moreover, the consideration of including phage therapy in modern medicine calls for more careful research around every aspect of phage study. In an attempt to adequately prepare for the events of phage resistance, many studies have attempted to experimentally evolve phage resistance in different bacterial strains, as well as train phages to evolve counter-infectivity of resistant bacterial mutants, in view of answering such questions as coevolutionary dynamics between phage and bacteria, mechanisms of phage resistance, fitness costs of phage resistance on bacteria, etc. In this review, we summarised many such studies and by careful examination, highlighted critical issues to the outcome of phage therapy. We also discuss the insufficiency of many of these in vitro studies to represent actual disease conditions during phage application, alongside other complications that exist in phage-bacterial evolutionary interactions. Conclusively, we present the exploitation of phage-bacterial interactions for successful infection managements, as well as some future perspectives to direct phage research.

Draft genome sequencing data of a feather mite, Amerodectes protonotaria Hernandes 2018 (Acariformes: Proctophyllodidae)

Feather mites are ubiquitous, permanent, obligate ectosymbionts of avian hosts and are a valuable natural system for studying host-symbiont evolutionary and ecological dynamics at multiple levels of biological organization. However, a lack of a sequenced genome impedes molecular studies using this system. Therefore, we present the first draft genome of a symbiotic feather mite, Amerodectes protonotaria Hernandes 2018. The genome sequence data presented here were derived from an individual female mite that was collected in the field from Protonotaria citrea, its only known host species. Short read sequence data were obtained using an Illumina NovaSeq 6000 platform. From these data, we assembled a 59,665,063 bp draft genome consisting of 2,399 contigs. Raw short reads and the assembled genome sequence are available at the National Center for Biotechnology Information (NCBI)'s Sequence Read Archive (SRA) under BioProject PRJNA884722. The data presented here are beneficial for future research on the biology and evolution of closely related mites and the genomics of host-symbiont interactions.

Degree of Dietary Specialization on Furanocoumarin-Containing Hostplants in a Newly Invasive Web Building Caterpillar

The genus Depressaria (Lepidoptera: Depressariidae) mostly comprises specialist herbivores with varying capacity for detoxification of defensive phytochemistry. Depressaria depressana, a Eurasian moth recently introduced into North America, is a family-level specialist of the Apiaceae, whose hosts include more than a dozen species in multiple tribes; Depressaria radiella is a super-specialist of Eurasian origin that feeds exclusively on species in the genera Pastinaca and Heracleum throughout its native and introduced range. In eastern North America, it feeds upon Pastinaca sativa, an invasive European species, and Heracleum maximum, a native species. We determined whether differences in furanocoumarin metabolism exist between D. depressana and two isolated populations of D. radiella, feeding exclusively on either P. sativa or H. maximum. We also compared gravimetric estimates of feeding efficiency to assess D. depressana larval performance on different diets. Both populations of D. radiella metabolized furanocoumarins at a greater rate than D. depressana. Although there was no difference in rates of metabolism of linear furanocoumarins in the two populations of D. radiella, individuals collected from H. maximum metabolized angular furanocoumarins more rapidly. The gravimetric assessments of feeding efficiencies revealed that D. depressana exhibited highest efficiencies consuming Daucus carota; moreover, this species survived to pupation consuming fruits of Zizia aurea, an apiaceous species native to North America. Our preliminary phylogenetic analysis, building on an earlier morphological analysis, incorporates mitochondrial cytochrome oxidase subunit 1 data from the BOLD database and revealed that the presence or absence of furanocoumarins is not a strong predictor of species-level evolution in Depressaria.

Fruit scent as an honest signal for fruit quality

BACKGROUND

Fleshy fruits evolved to be attractive to seed dispersers through various signals such as color and scent. Signals can evolve through different trajectories and have various degrees of reliability. The strongest substrate on which reliable signals can evolve is when there is an inherent link between signal and reward, rendering cheating costly or impossible. It was recently proposed that aliphatic esters in fruit scent may be predictive of sugar content due to their synthesis from products of sugar fermentation. We test this hypothesis on a case study of wild fig species (Ficus tiliifolia) from Madagascar, which relies on seed dispersal by lemurs.

RESULTS

We found a strong positive correlation between signal (esters) and reward (sugar). We also found that non-esters, including direct fermentation products, in fruit scent do not indicate sugar levels, which implies that this relationship is not simply a product of fruit maturation wherein more mature fruits emit more scent and contain more sugar.

CONCLUSIONS

While based on a single taxon, these results strongly support the hypothesis that a biochemical link between ester synthesis and sugar may render the ester fraction of fruit scent an honest signal for fruit quality, with consequences for animal sensory and feeding ecology, and the evolution of plants in the context of seed dispersal.

Evolutionary consequences of vector-borne transmission: how using vectors shapes host, vector and pathogen evolution

Transmission mode is a key factor that influences host–parasite coevolution. Vector-borne pathogens are among the most important disease agents for humans and wildlife due to their broad distribution, high diversity, prevalence and lethality. They comprise some of the most important and widespread human pathogens, such as yellow fever, leishmania and malaria. Vector-borne parasites (in this review, those transmitted by blood-feeding Diptera) follow unique transmission routes towards their vertebrate hosts. Consequently, each part of this tri-partite (i.e. parasite, vector and host) interaction can influence co- and counter-evolutionary pressures among antagonists. This mode of transmission may favour the evolution of greater virulence to the vertebrate host; however, pathogen–vector interactions can also have a broad spectrum of fitness costs to the insect vector. To complete their life cycle, vector-borne pathogens must overcome immune responses from 2 unrelated organisms, since they can activate responses in both vertebrate and invertebrate hosts, possibly creating a trade-off between investments against both types of immunity. Here, we assess how dipteran vector-borne transmission shapes the evolution of hosts, vectors and the pathogens themselves. Hosts, vectors and pathogens co-evolve together in a constant antagonistic arms race with each participant's primary goal being to maximize its performance and fitness.

Global diversity and distribution of prophages are lineage-specific within the Ralstonia solanacearum species complex

Background

Ralstonia solanacearum species complex (RSSC) strains are destructive plant pathogenic bacteria and the causative agents of bacterial wilt disease, infecting over 200 plant species worldwide. In addition to chromosomal genes, their virulence is mediated by mobile genetic elements including integrated DNA of bacteriophages, i.e., prophages, which may carry fitness-associated auxiliary genes or modulate host gene expression. Although experimental studies have characterised several prophages that shape RSSC virulence, the global diversity, distribution, and wider functional gene content of RSSC prophages are unknown. In this study, prophages were identified in a diverse collection of 192 RSSC draft genome assemblies originating from six continents.

Results

Prophages were identified bioinformatically and their diversity investigated using genetic distance measures, gene content, GC, and total length. Prophage distributions were characterised using metadata on RSSC strain geographic origin and lineage classification (phylotypes), and their functional gene content was assessed by identifying putative prophage-encoded auxiliary genes. In total, 313 intact prophages were identified, forming ten genetically distinct clusters. These included six prophage clusters with similarity to the Inoviridae, Myoviridae, and Siphoviridae phage families, and four uncharacterised clusters, possibly representing novel, previously undescribed phages. The prophages had broad geographical distributions, being present across multiple continents. However, they were generally host phylogenetic lineage-specific, and overall, prophage diversity was proportional to the genetic diversity of their hosts. The prophages contained many auxiliary genes involved in metabolism and virulence of both phage and bacteria.

Conclusions

Our results show that while RSSC prophages are highly diverse globally, they make lineage-specific contributions to the RSSC accessory genome, which could have resulted from shared coevolutionary history.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08909-7.

Genetic diversification of an invasive honey bee ectoparasite across sympatric and allopatric host populations

Invasive parasites are major threats to biodiversity. The honey bee ectoparasite, Varroa destructor, has shifted host and spread almost globally several decades ago. This pest is generally considered to be the main global threat to Western honey bees, Apis mellifera, although the damages it causes are not equivalent in all its new host's populations. Due to the high virulence of this parasite and the viruses it vectors, beekeepers generally rely on acaricide treatments to keep their colonies alive. However, some populations of A. mellifera can survive without anthropogenic mite control, through the expression of diverse resistance and tolerance traits. Such surviving colonies are currently found throughout the globe, with the biggest populations being found in Sub-Saharan Africa and Latin America. Recently, genetic differences between mite populations infesting surviving and treated A. mellifera colonies in Europe were found, suggesting that adaptations of honey bees drive mite evolution. Yet, the prevalence of such co-evolutionary adaptations in other invasive populations of V. destructor remain unknown. Using the previous data from Europe and novel genetic data from V. destructor populations in South America and Africa, we here investigated whether mites display signs of adaptations to different host populations of diverse origins and undergoing differing management. Our results show that, contrary to the differences previously documented in Europe, mites infesting treated and untreated honey bee populations in Africa and South America are genetically similar. However, strong levels of genetic differentiation were found when comparing mites across continents, suggesting ongoing allopatric speciation despite a recent spread from genetically homogenous lineages. This study provides novel insights into the co-evolution of V. destructor and A. mellifera, and confirms that these species are ideal to investigate coevolution in newly established host-parasite systems.

Evidence for saponin diversity-mycobiome links and conservatism of plant-fungi interaction patterns across Holarctic disjunct Panax species

Although interplays between plant and coevolved microorganisms are believed to drive landscape formation and ecosystem services, the relationships between the mycobiome and phytochemical evolution and the evolutionary characteristics of plant-mycobiome interaction patterns are still unclear. The present study explored fungal communities from 405 multiniche samples of three Holarctic disjunct Panax species. The overall mycobiomes showed compartment-dominated variations and dynamic universality. Neutral models were fitted for each compartment at the Panax genus (I) and species (II) levels to infer the community assembly mechanism and identify fungal subgroups potentially representing different plant-fungi interaction results, i.e., the potentially selected, opposed, and neutral taxa. Selection contributed more to the endosphere than to external compartments. The nonneutral taxa showed significant phylogenetic clustering. In Model I, the opposed subgroups could best reflect Panax saponin diversities (r = 0.69), and genera with highly positive correlations to specific saponins were identified using machine learning. Although mycobiomes in the three species differed significantly, subgroups in Model II were phylogenetically clustered based on potential interaction type rather than plant species, indicating potentially conservative plant-fungi interactions. In summary, the finding of strong links between invaders and saponin diversity can help explore the underlying mechanisms of saponin biosynthesis evolution from microbial insights, which is important to understanding the formation of the current landscape. The potential conservatism of plant-fungi interaction patterns suggests that the related genetic modules and selection pressures were convergent across Panax species, advancing our understanding of plant interplay with biotic environments.

Epochal coevolution of minor capsid protein in norovirus GII.4 variants with major capsid protein based on their interactions over the last five decades

Norovirus is a leading cause of viral gastroenteritis outbreaks worldwide, with GII.4 responsible for the majority of infections. Minor capsid protein VP2 has been found to have functions such as stabilizing virus particles, and VP2 is one of the highly variable proteins of norovirus, similar to major capsid protein VP1. However, whether the variation of VP2 is functionally driven still remains unclear. In this study, VP2 showed a higher evolutionary rate (2.642×10^-3 substitutions/site/year) than VP1 (1.587×10^-3 substitutions/site/year), and a hypervariable region in VP2 in a serial of norovirus GII.4 over the past 50 years had been observed. Notably, the high variation of VP2 was not haphazard. The evolutionary process of VP2 is similar to that of VP1 with comparable topologies when the phylogenetic trees were constructed. Moreover, VP2 was found to interact with VP1 among epidemic variants of GII.4 using the yeast two-hybrid experiments. The results of interactions were grouped into time-adjacent (e.g. Ancestral-VP1 plus US95-VP2) and non-adjacent (e.g. Ancestral-VP1 plus Sydney-VP2) according to the epochal chronologically based prevalence of GII.4 norovirus. Interestingly, the interaction of the former group was significantly stronger than that of the latter group (P=0.0001). Furthermore, the interaction regions on VP2 (residues 131-160 and 171-180) were mapped to the hypervariable region. And these interaction regions did show an important role in the evolutionary process of VP2, which was consistent with that of VP1. In summary, the minor capsid protein VP2 of GII.4 noroviruses had shown the epochal coevolution with VP1 based on their interactions over the past 50 years. The findings of this study provided valuable information for further understanding and completing the evolutionary mechanism of norovirus.

Ehrlich and Raven escape and radiate coevolution hypothesis at different levels of organization: Past and future perspectives

The classic paper by Ehrlich and Raven on coevolution will soon be 60 years old. Although they were not the first to develop the idea of coevolution, their thought-provoking paper certainly popularized this idea and inspired several generations of scientists interested in coevolution. Here, we describe some of their main contributions, quantitatively measure the impact of their seminal paper on different fields of research, and discuss how ideas related to their original paper might push the study of coevolution forward. To guide our discussion, we explore their original hypothesis into three research fields that are associated with distinct scales/levels of organization: (1) the genetic mechanisms underlying coevolutionary interactions; (2) the potential association between coevolutionary diversification and the organization of ecological networks; and (3) the micro- and macroevolutionary mechanisms and expected patterns under their hypothesis. By doing so, we discuss potentially overlooked aspects and future directions for the study of coevolutionary dynamics and diversification.

Human parasites as tracers of the evolution of their hosts

Because parasite data reveal essential information about the behavior and history of their hosts, it is possible to use them as tracers of host evolution. A table built from the analysis of the data contained in the book by Ashford and Crewe "The Parasites of Homo sapiens" allows counting and cross comparing the parasites according to the main descriptors used by the authors: Taxonomic groups, for each group number of parasites species identified in humans; Status, numbers of reported human cases and their dispersion; Geographic distribution, parasite specific richness recorded in biogeographic regions; Habitat, parasite location in or on the human body; Transmission, contamination pathways to man; Hosts, non-human hosts, which have a role in the maintenance of a parasite; Host-specificity status, relative role of man or other hosts in the maintenance of parasite populations. A strong positive correlation is observed between the number of parasites species recorded in humans and the global parasite species richness for each taxonomic group. About 74% of the parasites recorded in humans are rare, sporadic or nowhere common; 10% only are common or abundant worldwide. The Palearctic exhibits the highest parasite species diversity; the Oriental, Nearctic, Neotropical and Aethiopian regions have roughly similar richness values; the Australian Region is the poorest. Earliest domesticated animals, such as dog, cat, cattle or pig, share more parasite species with Humans than tardily domesticated as horse, rabbit or camel. More than one third of our parasites have elected our alimentary canal as a home and about two third are using the digestive tract path for contamination. Time of occupancy of new territories, diversity in feeding habits and commensalism with other animals, widely explain Human particular parasite richness. As suggested by the authors: "There must be few parasitic species which have never had the opportunity to infect a human".

Interacting phenotypes and the coevolutionary process: Interspecific indirect genetic effects alter coevolutionary dynamics

Coevolution occurs when species interact to influence one another's fitness, resulting in reciprocal evolutionary change. In many coevolving lineages, trait expression in one species is modified by the genotypes and phenotypes of the other, forming feedback loops reminiscent of models of intraspecific social evolution. Here, we adapt the theory of within-species social evolution, characterized by indirect genetic effects and social selection imposed by interacting individuals, to the case of interspecific interactions. In a trait-based model, we derive general expressions for multivariate evolutionary change in two species and the expected between-species covariance in evolutionary change when selection varies across space. We show that reciprocal interspecific indirect genetic effects can dominate the coevolutionary process and drive patterns of correlated evolution beyond what is expected from direct selection alone. In extreme cases, interspecific indirect genetic effects can lead to coevolution when selection does not covary between species or even when one species lacks genetic variance. Moreover, our model indicates that interspecific indirect genetic effects may interact in complex ways with cross-species selection to determine the course of coevolution. Importantly, our model makes empirically testable predictions for how different forms of reciprocal interactions contribute to the coevolutionary process.

The S100A7 nuclear interactors in autoimmune diseases: a coevolutionary study in mammals

S100A7, a member of the S100A family of Ca²⁺-binding proteins, is considered a key effector in immune response. In particular, S100A7 dysregulation has been associated with several diseases, including autoimmune disorders. At the nuclear level, S100A7 interacts with several protein-binding partners which are involved in transcriptional regulation and DNA repair. By using the BioGRID and GAAD databases, S100A7 nuclear interactors with a putative involvement in autoimmune diseases were retrieved. We selected fatty acid-binding protein 5 (FABP5), autoimmune regulator (AIRE), cystic fibrosis transmembrane conductance regulator (CFTR), chromodomain helicase DNA-binding protein 4 (CHD4), epidermal growth factor receptor (EGFR), estrogen receptor 1 (ESR1), histone deacetylase 2 (HDAC2), v-myc avian myelocytomatosis viral oncogene homolog (MYC), protection of telomeres protein 1 (POT1), telomeric repeat-binding factor (NIMA-interacting) 1 (TERF1), telomeric repeat-binding factor 2 (TERF2), and Zic family member 1 (ZIC1). Linear correlation coefficients between interprotein distances were calculated with MirrorTree. Coevolution clusters were also identified with the use of a recent version of the Blocks in Sequences (BIS2) algorithm implemented in the BIS2Analyzer web server. Analysis of pair positions identified interprotein coevolving clusters between S100A7 and the binding partners CFTR and TERF1. Such findings could guide further analysis to better elucidate the function of S100A7 and its binding partners and to design drugs targeting for these molecules in autoimmune diseases.

Chemical Ecology of the North American Newt Genera Taricha and Notophthalmus

The North American newt genera Taricha and Notophthalmus (order Caudata) are well known for the combination of potent toxicity, aposematic coloration, and striking defense postures that protects these animals from predation. This suite of traits is centered around the neurotoxin tetrodotoxin, which causes paralysis and death in metazoans by disrupting the initiation and propagation of electrical signals in the nerves and muscles. Tetrodotoxin defends newts from predation across multiple life history stages and its role in generating arms-race coevolution between Taricha newts and garter snake (genus Thamnophis) predators is well studied. However, understanding the broader picture of chemical defenses in Taricha and Notophthalmus requires an expanded comprehension of the defensive chemical ecology of tetrodotoxin that includes possible coevolutionary interactions with insect egg predators, protection against parasites, as well as mimicry complexes associated with tetrodotoxin and aposematic coloration in both genera. Herein the authors review what is known about the structure, function, and pharmacology of tetrodotoxin to explore its evolution and chemical ecology in the North American newt. Focus is made specifically on the origin and possible biosynthesis of tetrodotoxin in these taxa as well as providing an expanded picture of the web of interactions that contribute to landscape level patterns of toxicity and defense in Taricha and Notophthalmus.

The genus Entomophthora: bringing the insect destroyers into the twenty-first century

The fungal genus Entomophthora consists of highly host-specific pathogens that cause deadly epizootics in their various insect hosts. The most well-known among these is the "zombie fly" fungus E. muscae, which, like other Entomophthora species, elicits a series of dramatic behaviors in infected hosts to promote optimal spore dispersal. Despite having been first described more than 160 years ago, there are still many open questions about Entomophthora biology, including the molecular underpinnings of host behavior manipulation and host specificity. This review provides a comprehensive overview of our current understanding of the biology of Entomophthora fungi and enumerates the most pressing outstanding questions that should be addressed in the field. We briefly review the discovery of Entomophthora and provide a summary of the 21 recognized Entomophthora species, including their type hosts, methods of transmission (ejection of spores after or before host death), and for which molecular data are available. Further, we argue that this genus is globally distributed, based on a compilation of Entomophthora records in the literature and in online naturalist databases, and likely to contain additional species. Evidence for strain-level specificity of hosts is summarized and directly compared to phylogenies of Entomophthora and the class Insecta. A detailed description of Entomophthora's life-cycle and observed manipulated behaviors is provided and used to summarize a consensus for ideal growth conditions. We discuss evidence for Entomophthora's adaptation to growth exclusively inside insects, such as producing wall-less hyphal bodies and a unique set of subtilisin-like proteases to penetrate the insect cuticle. However, we are only starting to understand the functions of unusual molecular and genomic characteristics, such as having large > 1 Gb genomes full of repetitive elements and potential functional diploidy. We argue that the high host-specificity and obligate life-style of most Entomophthora species provides ample scope for having been shaped by close coevolution with insects despite the current general lack of such evidence. Finally, we propose six major directions for future Entomophthora research and in doing so hope to provide a foundation for future studies of these fungi and their interaction with insects.

A geographic mosaic of coevolution between Eurosta solidaginis (Fitch) and its host plant tall goldenrod Solidago altissima (L.)

A geographic mosaic of coevolution has produced local reciprocal adaptation in tall goldenrod, Solidago altissima (L.), and the goldenrod ball-gall fly, Eurosta solidaginis (Fitch 1855). The fly is selected to induce gall diameters that minimize mortality from natural enemies, and the plant is selected to limit gall growth that reduces plant fitness. We conducted a double reciprocal transplant experiment where S. altissima and E. solidaginis from three sites were grown in gardens at each site to partition the gall morphology variation into fly genotype, plant genotype, and the environment components. The host plant gall diameter induced by each E. solidaginis population was adapted to inhibit local natural enemies from ovipositing on or consuming enclosed larvae. Reciprocally, increasing the gall size induced by the local fly population increased the resistance of the local plant host population to gall growth. Differences among sites in natural enemies produced a mosaic of hotspots of coevolutionary arms races between flies selecting for greater gall diameter and plants for smaller diameters, and coldspots where there is no selection on plant or fly for a change in gall diameter. In contrast, the geographic variations of gall length and gall shape were not due to coevolutionary interactions.

Co-structure analysis and genetic associations reveal insights into pinworms (Trypanoxyuris) and primates (Alouatta palliata) microevolutionary dynamics

BACKGROUND

In parasitism arm race processes and red queen dynamics between host and parasites reciprocally mold many aspects of their genetics and evolution. We performed a parallel assessment of population genetics and demography of two species of pinworms with different degrees of host specificity (Trypanoxyuris multilabiatus, species-specific; and T. minutus, genus-specific) and their host, the mantled howler monkey (Alouatta palliata), based on mitochondrial DNA sequences and microsatellite loci (these only for the host). Given that pinworms and primates have a close co-evolutionary history, covariation in several genetic aspects of their populations is expected.

RESULTS

Mitochondrial DNA revealed two genetic clusters (West and East) in both pinworm species and howler monkeys, although population structure and genetic differentiation were stronger in the host, while genetic diversity was higher in pinworms than howler populations. Co-divergence tests showed no congruence between host and parasite phylogenies; nonetheless, a significant correlation was found between both pinworms and A. palliata genetic pairwise distances suggesting that the parasites' gene flow is mediated by the host dispersal. Moreover, the parasite most infective and the host most susceptible haplotypes were also the most frequent, whereas the less divergent haplotypes tended to be either more infective (for pinworms) or more susceptible (for howlers). Finally, a positive correlation was found between pairwise p-distance of host haplotypes and that of their associated pinworm haplotypes.

CONCLUSION

The genetic configuration of pinworm populations appears to be molded by their own demography and life history traits in conjunction with the biology and evolutionary history of their hosts, including host genetic variation, social interactions, dispersal and biogeography. Similarity in patterns of genetic structure, differentiation and diversity is higher between howler monkeys and T. multilabiatus in comparison with T. minutus, highlighting the role of host-specificity in coevolving processes. Trypanoxyuris minutus exhibits genetic specificity towards the most frequent host haplotype as well as geographic specificity. Results suggest signals of potential local adaptation in pinworms and further support the notion of correlated evolution between pinworms and their primate hosts.

Novel ACE2 protein interactions relevant to COVID-19 predicted by evolutionary rate correlations

Angiotensin-converting enzyme 2 (ACE2) is the cell receptor that the coronavirus SARS-CoV-2 binds to and uses to enter and infect human cells. COVID-19, the pandemic disease caused by the coronavirus, involves diverse pathologies beyond those of a respiratory disease, including micro-thrombosis (micro-clotting), cytokine storms, and inflammatory responses affecting many organ systems. Longer-term chronic illness can persist for many months, often well after the pathogen is no longer detected. A better understanding of the proteins that ACE2 interacts with can reveal information relevant to these disease manifestations and possible avenues for treatment. We have undertaken an approach to predict candidate ACE2 interacting proteins which uses evolutionary inference to identify a set of mammalian proteins that “coevolve” with ACE2. The approach, called evolutionary rate correlation (ERC), detects proteins that show highly correlated evolutionary rates during mammalian evolution. Such proteins are candidates for biological interactions with the ACE2 receptor. The approach has uncovered a number of key ACE2 protein interactions of potential relevance to COVID-19 pathologies. Some proteins have previously been reported to be associated with severe COVID-19, but are not currently known to interact with ACE2, while additional predicted novel ACE2 interactors are of potential relevance to the disease. Using reciprocal rankings of protein ERCs, we have identified strongly interconnected ACE2 associated protein networks relevant to COVID-19 pathologies. ACE2 has clear connections to coagulation pathway proteins, such as Coagulation Factor V and fibrinogen components FGA, FGB, and FGG, the latter possibly mediated through ACE2 connections to Clusterin (which clears misfolded extracellular proteins) and GPR141 (whose functions are relatively unknown). ACE2 also connects to proteins involved in cytokine signaling and immune response (e.g. XCR1, IFNAR2 and TLR8), and to Androgen Receptor (AR). The ERC prescreening approach has elucidated possible functions for relatively uncharacterized proteins and possible new functions for well-characterized ones. Suggestions are made for the validation of ERC-predicted ACE2 protein interactions. We propose that ACE2 has novel protein interactions that are disrupted during SARS-CoV-2 infection, contributing to the spectrum of COVID-19 pathologies.

Molecular Phylogenetic Position of Microjoenia (Parabasalia: Spirotrichonymphea) from Reticulitermes and Hodotermopsis Termite Hosts

Microjoenia are obligate symbionts of termites. The genus was erected in 1892 for small cells with many flagella that insert near, but not directly from, the cell apex, and an axostyle that can protrude from the cell posterior. Although ultrastructural studies have been carried out on three Microjoenia species to date, no molecular data have been directly attributed to any species. Microjoenia are classified within the parabasalian class Spirotrichonymphea, which is characterized by flagellar bands that emerge near the cell apex and proceed posteriorly in a right-handed helix. In Microjoenia, however, the flagellar bands are very short and proceed longitudinally or with a weakly observable helix. In this study, we have amplified and sequenced the 18S ribosomal RNA gene from individually isolated Microjoenia cells from Reticulitermes and Hodotermopsis hosts as part of an ongoing effort to understand the phylogeny of Spirotrichonymphea and their coevolution with termites. In our 18S rRNA gene phylogeny, Microjoenia forms the sister lineage to Spirotrichonympha, though many other evolutionary relationships within Spirotrichonymphea remain unresolved.

The coevolution of large and small terminases of bacteriophages is a result of purifying selection leading to phenotypic stabilization

Genome packaging in many dsDNA phages requires a series of precisely coordinated actions of two phage-coded proteins, namely, large terminase (TerL) and small terminase (TerS) with DNA and ATP, and with each other. Despite the strict functional conservation, TerL and TerS homologs exhibit large sequence variations. We investigated the sequence variability across eight phage types and observed a coevolutionary framework wherein the genealogy of TerL homologs mirrored that of the corresponding TerS homologs. Furthermore, a high purifying selection observed (dN/dS«1) indicated strong structural constraints on both TerL and TerS, and identify coevolving residues in TerL and TerS of phage T4 and lambda. Using the highly coevolving (correlation coefficient of 0.99) TerL and TerS of phage N4, we show that their biochemical features are similar to the phylogenetically divergent phage λ terminases. We also demonstrate using the Surface Plasma Resonance (SPR) technique that phage N4 TerL transiently interacts with TerS.

Evolution of nuptial gifts and its coevolutionary dynamics with male-like persistence traits of females for multiple mating

BACKGROUND

Many male animals donate nutritive materials during courtship or mating to their female mates. Donation of large-sized gifts, though costly to prepare, can result in increased sperm transfer during mating and delayed remating of the females, resulting in higher paternity. Nuptial gifting sometimes causes severe female-female competition for obtaining gifts (i.e., sex-role reversal in mate competition) and selection on females to increase their mating rate, changing the intensity of sperm competition and the resultant paternity gains. We built a theoretical model to simulate such coevolutionary feedbacks between nuptial gift size (male trait) and propensity for multiple mating (female trait). Donation of nuptial gifts sometimes causes development of female persistence trait for gift acquisition. We also analyzed the causes and consequences of this type of traits, taking double receptacles for nutritious seminal gifts, which are known to occur in an insect group with a "female penis" (Neotrogla spp.), as an illustrative example.

RESULTS

Our individual-based simulations demonstrated that female-female competition for male-derived nutrients always occur when the environment is oligotrophic and mating costs are low for females. However, a positive correlation between donated gift size and the resultant paternity gain was a requisite for the co-occurrence of large gifts and females' competitive multiple mating for the gifts. When gift donation satisfied female demands and thus resulted in monandry, exaggeration of nuptial gift size also occurred under the assumption that the last male monopolizes paternity. The evolution of double slots for gift acquisition and digestion (female persistence trait) always occurred when males could not satisfy the demands of females for gifts. However, through coevolutionary reduction in male gift size, fixation of this trait in a population drastically reduced the average female fitness.

CONCLUSION

Sperm usage patterns, which have rarely been examined for animals with nuptial gifts, can be a critical factor for determining the extent of exaggeration in nuptial gifting. Sex-role reversals in mate competition, as a result of donation of nuptial gifts from males to females, can involve the evolution of male-like, persistent traits in females that reduce population productivity, as is the case with persistence traits in males.

Symbiosis and the Anthropocene

Recent human activity has profoundly transformed Earth biomes on a scale and at rates that are unprecedented. Given the central role of symbioses in ecosystem processes, functions, and services throughout the Earth biosphere, the impacts of human-driven change on symbioses are critical to understand. Symbioses are not merely collections of organisms, but co-evolved partners that arise from the synergistic combination and action of different genetic programs. They function with varying degrees of permanence and selection as emergent units with substantial potential for combinatorial and evolutionary innovation in both structure and function. Following an articulation of operational definitions of symbiosis and related concepts and characteristics of the Anthropocene, we outline a basic typology of anthropogenic change (AC) and a conceptual framework for how AC might mechanistically impact symbioses with select case examples to highlight our perspective. We discuss surprising connections between symbiosis and the Anthropocene, suggesting ways in which new symbioses could arise due to AC, how symbioses could be agents of ecosystem change, and how symbioses, broadly defined, of humans and “farmed” organisms may have launched the Anthropocene. We conclude with reflections on the robustness of symbioses to AC and our perspective on the importance of symbioses as ecosystem keystones and the need to tackle anthropogenic challenges as wise and humble stewards embedded within the system.

Interspecific interactions that affect ageing: Age-distorters manipulate host ageing to their own evolutionary benefits

Genetic causes for ageing are traditionally investigated within a species. Yet, the lifecycles of many organisms intersect. Additional evolutionary and genetic causes of ageing, external to a focal species/organism, may thus be overlooked. Here, we introduce the phrase and concept of age-distorters and its evidence. Age-distorters carry ageing interfering genes, used to manipulate the biological age of other entities upon which the reproduction of age-distorters relies, e.g. age-distorters bias the reproduction/maintenance trade-offs of cells/organisms for their own evolutionary interests. Candidate age-distorters include viruses, parasites and symbionts, operating through specific, genetically encoded interferences resulting from co-evolution and arms race between manipulative non-kins and manipulable species. This interference results in organismal ageing when age-distorters prompt manipulated organisms to favor their reproduction at the expense of their maintenance, turning these hosts into expanded disposable soma. By relying on reproduction/maintenance trade-offs affecting disposable entities, which are left ageing to the reproductive benefit of other physically connected lineages with conflicting evolutionary interests, the concept of age-distorters expands the logic of the Disposable Soma theory beyond species with fixed germen/soma distinctions. Moreover, acknowledging age-distorters as external sources of mutation accumulation and antagonistic pleiotropic genes expands the scope of the mutation accumulation and of the antagonistic pleiotropy theories.

The interplay between host biogeography and phylogeny in structuring diversification of the feather louse genus Penenirmus

Parasite diversification is influenced by many of the same factors that affect speciation of free-living organisms, such as biogeographic barriers. However, the ecology and evolution of the host lineage also has a major impact on parasite speciation. Here we explore the interplay between biogeography and host-association on the pattern of diversification in a group of ectoparasitic lice (Insecta: Phthiraptera: Penenirmus) that feeds on the feathers of woodpeckers, barbets, and honeyguides (Piciformes) and some songbirds (Passeriformes). We use whole genome sequencing of 41 ingroup and 12 outgroup samples to develop a phylogenomic dataset of DNA sequences from a reference set of 2395 single copy ortholog genes, for a total of nearly four million aligned base positions. The phylogenetic trees resulting from both concatenated and gene-tree/species-tree coalescent analyses were nearly identical and highly supported. These trees recovered the genus Penenirmus as monophyletic and identified several major clades, which tended to be associated with one major host group. However, cophylogenetic analysis revealed that host-switching was a prominent process in the diversification of this group. This host-switching generally occurred within single major biogeographic regions. We did, however, find one case in which it appears that a rare dispersal event by a woodpecker lineage from North America to Africa allowed its associated louse to colonize a woodpecker in Africa, even though the woodpecker lineage from North America never became established there.

The rise of ecological parasitology: twelve landmark advances that changed its history

In the five decades since the first publication of the International Journal for Parasitology, ecological parasitology has grown from modest beginnings to become a modern discipline with a strong theoretical foundation, a diverse toolkit, and a multidisciplinary approach. In this review, I highlight 12 advances in the field that have spurred its growth over the past 50 years. Where relevant, I identify pivotal contributions that have altered the course of research, as well as the influence of developments in other fields such as mainstream ecology and molecular biology. The 12 key advances discussed are in areas including parasite population dynamics and community assembly, the regulation of host population abundance and food web structure, parasites as agents of natural selection, the impacts of biodiversity and anthropogenic changes on host-parasite interactions, the biogeography of parasite diversity, and the evolutionary genetics of parasites. I conclude by identifying some challenges and opportunities lying ahead, which need to be met for the future growth of ecological research on host-parasite interactions.

Bog ecosystems as a playground for plant-microbe coevolution: bryophytes and vascular plants harbour functionally adapted bacteria

BACKGROUND

Bogs are unique ecosystems inhabited by distinctive, coevolved assemblages of organisms, which play a global role for carbon storage, climate stability, water quality and biodiversity. To understand ecology and plant-microbe co-occurrence in bogs, we selected 12 representative species of bryophytes and vascular plants and subjected them to a shotgun metagenomic sequencing approach. We explored specific plant-microbe associations as well as functional implications of the respective communities on their host plants and the bog ecosystem.

RESULTS

Microbial communities were shown to be functionally adapted to their plant hosts; a higher colonization specificity was found for vascular plants. Bryophytes that commonly constitute the predominant Sphagnum layer in bogs were characterized by a higher bacterial richness and diversity. Each plant group showed an enrichment of distinct phylogenetic and functional bacterial lineages. Detailed analyses of the metabolic potential of 28 metagenome-assembled genomes (MAGs) supported the observed functional specification of prevalent bacteria. We found that novel lineages of Betaproteobacteria and Actinobacteria in the bog environment harboured genes required for carbon fixation via RuBisCo. Interestingly, several of the highly abundant bacteria in both plant types harboured pathogenicity potential and carried similar virulence factors as found with corresponding human pathogens.

CONCLUSIONS

The unexpectedly high specificity of the plant microbiota reflects intimate plant-microbe interactions and coevolution in bog environments. We assume that the detected pathogenicity factors might be involved in coevolution processes, but the finding also reinforces the role of the natural plant microbiota as a potential reservoir for human pathogens. Overall, the study demonstrates how plant-microbe assemblages can ensure stability, functioning and ecosystem health in bogs. It also highlights the role of bog ecosystems as a playground for plant-microbe coevolution. Video abstract.

Cospeciation of coronavirus and paramyxovirus with their bat hosts in the same geographical areas

Background

Bat-borne viruses are relatively host specific. We hypothesize that this host specificity is due to coevolution of the viruses with their hosts. To test this hypothesis, we investigated the coevolution of coronavirus and paramyxovirus with their bat hosts. Published nucleotide sequences of the RNA-dependent RNA polymerase (RdRp) gene of 60 coronavirus strains identified from 37 bat species, the RNA polymerase large (L) gene of 36 paramyxovirus strains from 29 bat species, and the cytochrome B (cytB) gene of 35 bat species were analyzed for coevolution signals. Each coevolution signal detected was tested and verified by global-fit cophylogenic analysis using software ParaFit, PACo, and eMPRess.

Results

Significant coevolution signals were detected in coronaviruses and paramyxoviruses and their bat hosts, and closely related bat hosts were found to carry closely related viruses.

Conclusions

Our results suggest that paramyxovirus and coronavirus coevolve with their hosts.