Host-parasite co-evolution and its genomic signature

Long-term balancing selection for pathogen resistance maintains trans-species polymorphisms in a planktonic crustacean

Balancing selection is an evolutionary process that maintains genetic polymorphisms at selected loci and strongly reduces the likelihood of allele fixation. When allelic polymorphisms that predate speciation events are maintained independently in the resulting lineages, a pattern of trans-species polymorphisms may occur. Trans-species polymorphisms have been identified for loci related to mating systems and the MHC, but they are generally rare. Trans-species polymorphisms in disease loci are believed to be a consequence of long-term host-parasite coevolution by balancing selection, the so-called Red Queen dynamics. Here we scan the genomes of three crustaceans with a divergence of over 15 million years and identify 11 genes containing identical-by-descent trans-species polymorphisms with the same polymorphisms in all three species. Four of these genes display molecular footprints of balancing selection and have a function related to immunity. Three of them are located in or close to loci involved in resistance to a virulent bacterial pathogen, Pasteuria, with which the Daphnia host is known to coevolve. This provides rare evidence of trans-species polymorphisms for loci known to be functionally relevant in interactions with a widespread and highly specific parasite. These findings support the theory that specific antagonistic coevolution is able to maintain genetic diversity over millions of years.

Greater than the sum of their parts: an overview of the AvrRps4 effector family

Phytopathogenic microbes use secreted effector proteins to increase their virulence in planta. If these effectors or the results of their activity are detected by the plant cell, the plant will mount an immune response which applies evolutionary pressure by reducing growth and success of the pathogen. Bacterial effector proteins in the AvrRps4 family (AvrRps4, HopK1, and XopO) have commonly been used as tools to investigate plant immune components. At the same time, the in planta functions of this family of effectors have yet to be fully characterized. In this minireview we summarize current knowledge about the AvrRps4 effector family with emphasis on properties of the proteins themselves. We hypothesize that the HopK1 C-terminus and the AvrRps4 C-terminus, though unrelated in sequence and structure, are broadly related in functions that counteract plant defense responses.

A new dawn for the study of cell type evolution

Animal evolution is influenced by the emergence of new cell types, yet our understanding of this process remains elusive. This prompts the need for a broader exploration across diverse research organisms, facilitated by recent breakthroughs, such as gene editing tools and single-cell genomics. Essential to our understanding of cell type evolution is the accurate identification of homologous cells. We delve into the significance of considering developmental ontogeny and potential pitfalls when drawing conclusions about cell type homology. Additionally, we highlight recent discoveries in the study of cell type evolution through the application of single-cell transcriptomics and pinpoint areas ripe for further exploration.

Deep population structure linked to host vernalization requirement in the barley net blotch fungal pathogen

Invasive fungal pathogens pose a substantial threat to widely cultivated crop species, owing to their capacity to adapt to new hosts and new environmental conditions. Gaining insights into the demographic history of these pathogens and unravelling the mechanisms driving coevolutionary processes are crucial for developing durably effective disease management programmes. Pyrenophora teres is a significant fungal pathogen of barley, consisting of two lineages, Ptt and Ptm, with global distributions and demographic histories reflecting barley domestication and spread. However, the factors influencing the population structure of P. teres remain poorly understood, despite the varietal and environmental heterogeneity of barley agrosystems. Here, we report on the population genomic structure of P. teres in France and globally. We used genotyping-by-sequencing to show that Ptt and Ptm can coexist in the same area in France, with Ptt predominating. Furthermore, we showed that differences in the vernalization requirement of barley varieties were associated with population differentiation within Ptt in France and at a global scale, with one population cluster found on spring barley and another population cluster found on winter barley. Our results demonstrate how cultivation conditions, possibly associated with genetic differences between host populations, can be associated with the maintenance of divergent invasive pathogen populations coexisting over large geographic areas. This study not only advances our understanding of the coevolutionary dynamics of the Pt-barley pathosystem but also prompts further research on the relative contributions of adaptation to the host versus adaptation to abiotic conditions in shaping Ptt populations.

The germ theory revisited: A noncentric view on infection outcome

The germ theory states that pathogenic microorganisms are responsible for causing infectious diseases. The theory is inherently microbe-centric and does not account for variability in disease severity among individuals and asymptomatic carriership-two phenomena indicating an important role for host variability in infection outcome. The basic tenet of the germ theory was recently challenged, and a radically host-centric paradigm referred to as the "full-blown host theory" was proposed. According to this view, the pathogen is reduced to a passive environmental trigger, and the development of disease is instead due to pre-existing immunodeficiencies of the host. Here, we consider the factors that determine disease severity using established knowledge concerning evolutionary biology, microbial pathogenesis, and host-pathogen interactions. We note that the available data support a noncentric view that recognizes key roles for both the causative microbe and the host in dictating infection outcome.

Ectopical expression of bacterial collagen-like protein supports its role as adhesin in host-parasite coevolution

For a profound understanding of antagonistic coevolution, it is necessary to identify the coevolving genes. The bacterium Pasteuria and its host, the microcrustacean Daphnia, are a well-characterized paradigm for co-evolution, but the underlying genes remain largely unknown. A genome-wide association study suggested a Pasteuria collagen-like protein 7 (Pcl7) as a candidate mediating parasite attachment and driving its coevolution with the host. Since Pasteuria ramosa cannot currently be genetically manipulated, we used Bacillus thuringiensis to express a fusion protein of a Pcl7 carboxy-terminus from P. ramosa and the amino-terminal domain of a B. thuringiensis collagen-like protein (CLP). Mutant B. thuringiensis (Pcl7-Bt) spores but not wild-type B. thuringiensis (WT-Bt) spores attached to the same site of susceptible hosts as P. ramosa. Furthermore, Pcl7-Bt spores attached readily to susceptible host genotypes, but only slightly to resistant host genotypes. These findings indicated that the fusion protein was properly expressed and folded and demonstrated that indeed the C-terminus of Pcl7 mediates attachment in a host genotype-specific manner. These results provide strong evidence for the involvement of a CLP in the coevolution of Daphnia and P. ramosa and open new avenues for genetic epidemiological studies of host-parasite interactions.

Two-speed genomes of Epichloe fungal pathogens show contrasting signatures of selection between species and across populations

Antagonistic selection between pathogens and their hosts can drive rapid evolutionary change and leave distinct molecular footprints of past and ongoing selection in the genomes of the interacting species. Despite an increasing availability of tools able to identify signatures of selection, the genetic mechanisms underlying coevolutionary interactions and the specific genes involved are still poorly understood, especially in heterogeneous natural environments. We searched the genomes of two species of Epichloe plant pathogen for evidence of recent selection. The Epichloe genus includes highly host-specific species that can sterilize their grass hosts. We performed selection scans using genome-wide SNP data from seven natural populations of two co-occurring Epichloe sibling species specialized on different hosts. We found evidence of recent (and ongoing) selective sweeps across the genome in both species. However, selective sweeps were more abundant in the species with a larger effective population size. Sweep regions often overlapped with highly polymorphic AT-rich regions supporting the role of these genome compartments in adaptive evolution. Although most loci under selection were specific to individual populations, we could also identify several candidate genes targeted by selection in sweep regions shared among populations. The genes encoded small secreted proteins typical of fungal effectors and cell wall-degrading enzymes. By investigating the genomic signatures of selection across multiple populations and species, this study contributes to our understanding of complex adaptive processes in natural plant pathogen systems.

A synthesis of coevolution across levels of biological organization

In evolutionary ecology, coevolution is typically defined as reciprocal evolution of interacting species. However, outside the context of interacting species, the term "coevolution" is also used at levels of biological organization within species (e.g., between males and females, between cells, and between genes or proteins). Furthermore, although evolution is typically defined as "genetic change over time", coevolution need not involve genetic changes in the interacting parties, since cultures can also evolve. In this review, I propose that coevolution be defined more broadly as "reciprocal adaptive evolution at any level of biological organisation". The classification of reciprocal evolution at all levels of biological organization as coevolution would maintain consistency in terminology. More importantly, the broader definition should facilitate greater integration of coevolution research across disciplines. For example, principles usually discussed only in the context of coevolution between species or coevolution between genes (e.g., tight and diffuse coevolution, and compensatory coevolution, respectively) could be more readily applied to new fields. The application of coevolutionary principles to new contexts could also provide benefits to society, for instance in deducing the dynamics of coevolution between cancer cells and cells of the human immune system.

Whole-Genome Sequencing Analyses Reveal the Whip-like Tail Formation, Innate Immune Evolution, and DNA Repair Mechanisms of Eupleurogrammus muticus

Smallhead hairtail (Eupleurogrammus muticus) is an important marine economic fish distributed along the northern Indian Ocean and the northwest Pacific coast; however, little is known about the mechanism of its genetic evolution. This study generated the first genome assembly of E. muticus at the chromosomal level using a combination of PacBio SMRT, Illumina Nova-Seq, and Hi-C technologies. The final assembled genome size was 709.27 Mb, with a contig N50 of 25.07 Mb, GC content of 40.81%, heterozygosity rate of 1.18%, and repetitive sequence rate of 35.43%. E. muticus genome contained 21,949 protein-coding genes (97.92% of the genes were functionally annotated) and 24 chromosomes. There were 143 expansion gene families, 708 contraction gene families, and 4888 positively selected genes in the genome. Based on the comparative genomic analyses, we screened several candidate genes and pathways related to whip-like tail formation, innate immunity, and DNA repair in E. muticus. These findings preliminarily reveal some molecular evolutionary mechanisms of E. muticus at the genomic level and provide important reference genomic data for the genetic studies of other trichiurids.

Genome-wide support for incipient Tula hantavirus species within a single rodent host lineage

Evolutionary divergence of viruses is most commonly driven by co-divergence with their hosts or through isolation of transmission after host shifts. It remains mostly unknown, however, whether divergent phylogenetic clades within named virus species represent functionally equivalent byproducts of high evolutionary rates or rather incipient virus species. Here, we test these alternatives with genomic data from two widespread phylogenetic clades in Tula orthohantavirus (TULV) within a single evolutionary lineage of their natural rodent host, the common vole Microtus arvalis. We examined voles from forty-two locations in the contact region between clades for TULV infection by reverse transcription (RT)-PCR. Sequencing yielded twenty-three TULV Central North and twenty-one TULV Central South genomes, which differed by 14.9-18.5 per cent at the nucleotide and 2.2-3.7 per cent at the amino acid (AA) level without evidence of recombination or reassortment between clades. Geographic cline analyses demonstrated an abrupt (<1 km wide) transition between the parapatric TULV clades in continuous landscape. This transition was located within the Central mitochondrial lineage of M. arvalis, and genomic single nucleotide polymorphisms showed gradual mixing of host populations across it. Genomic differentiation of hosts was much weaker across the TULV Central North to South transition than across the nearby hybrid zone between two evolutionary lineages in the host. We suggest that these parapatric TULV clades represent functionally distinct, incipient species, which are likely differently affected by genetic polymorphisms in the host. This highlights the potential of natural viral contact zones as systems for investigating the genetic and evolutionary factors enabling or restricting the transmission of RNA viruses.

Effects of plant tissue permeability on invasion and population bottlenecks of a phytopathogen

Pathogen genetic diversity varies in response to environmental changes. However, it remains unclear whether plant barriers to invasion could be considered a genetic bottleneck for phytopathogen populations. Here, we implement a barcoding approach to generate a pool of 90 isogenic and individually barcoded Ralstonia solanacearum strains. We used 90 of these strains to inoculate tomato plants with different degrees of physical permeability to invasion (intact roots, wounded roots and xylem inoculation) and quantify the phytopathogen population dynamics during invasion. Our results reveal that the permeability of plant roots impacts the degree of population bottleneck, genetic diversity, and composition of Ralstonia populations. We also find that selection is the main driver structuring pathogen populations when barriers to infection are less permeable, i.e., intact roots, the removal of root physical and immune barriers results in the predominance of stochasticity in population assembly. Taken together, our study suggests that plant root permeability constitutes a bottleneck for phytopathogen invasion and genetic diversity.

Genetic difference between two Schistosoma japonicum isolates with contrasting cercarial shedding patterns revealed by whole genome sequencing

Schistosoma japonicum is one of the major infectious agents of human schistosomiasis, mainly endemic in China and the Philippines. We have previously reported the finding of two schistosome isolates, each with a different cercarial emergence pattern adapted to their different hosts. However, there are currently no whole-genome sequencing studies to investigate the underlining genetics of the adaptive traits. We sampled schistosomes in 2013 and 2020 from a hilly area Shitai (ST) and a marshland area Hexian (HX) of Anhui, China. Ten to 15 male or female adult worms from each site/year were sent for whole genome sequencing. Genetics were analyzed, and selection signals along genomes were detected. Gene enrichment analysis was performed for the genome regions under selection. The results revealed considerable genetic differentiation between the two isolates. The genome "windows" affected by natural selection were fewer in ST (64 windows containing 78 genes) than in HX (318 windows containing 276 genes). Twelve significantly enriched genes were identified in ST, but none in HX. These genes were mainly related to specific DNA binding and intercellular signaling transduction. Some functional region changes identified along the genome of the hilly schistosome may be related to its unique late afternoon cercarial emergence.

Chromosome-scale genome of the human blood fluke Schistosoma mekongi and its implications for public health

BACKGROUND

Schistosoma mekongi is a human blood fluke causing schistosomiasis that threatens approximately 1.5 million humans in the world. Nonetheless, the limited available S. mekongi genomic resources have hindered understanding of its biology and parasite-host interactions for disease management and pathogen control. The aim of our study was to integrate multiple technologies to construct a high-quality chromosome-level assembly of the S. mekongi genome.

METHODS

The reference genome for S. mekongi was generated through integrating Illumina, PacBio sequencing, 10 × Genomics linked-read sequencing, and high-throughput chromosome conformation capture (Hi-C) methods. In this study, we conducted de novo assembly, alignment, and gene prediction to assemble and annotate the genome. Comparative genomics allowed us to compare genomes across different species, shedding light on conserved regions and evolutionary relationships. Additionally, our transcriptomic analysis focused on genes associated with parasite-snail interactions in S. mekongi infection. We employed gene ontology (GO) enrichment analysis for functional annotation of these genes.

RESULTS

In the present study, the S. mekongi genome was both assembled into 8 pseudochromosomes with a length of 404 Mb, with contig N50 and scaffold N50 lengths of 1168 kb and 46,759 kb, respectively. We detected that 43% of the genome consists of repeat sequences and predicted 9103 protein-coding genes. We also focused on proteases, particularly leishmanolysin-like metalloproteases (M8), which are crucial in the invasion of hosts by 12 flatworm species. Through phylogenetic analysis, it was discovered that the M8 gene exhibits lineage-specific amplification among the genus Schistosoma. Lineage-specific expansion of M8 was observed in blood flukes. Additionally, the results of the RNA-seq revealed that a mass of genes related to metabolic and biosynthetic processes were up-regulated, which might be beneficial for cercaria production.

CONCLUSIONS

This study delivers a high-quality, chromosome-scale reference genome of S. mekongi, enhancing our understanding of the divergence and evolution of Schistosoma. The molecular research conducted here also plays a pivotal role in drug discovery and vaccine development. Furthermore, our work greatly advances the understanding of host-parasite interactions, providing crucial insights for schistosomiasis intervention strategies.

Ongoing shuffling of protein fragments diversifies core viral functions linked to interactions with bacterial hosts

Biological modularity enhances evolutionary adaptability. This principle is vividly exemplified by bacterial viruses (phages), which display extensive genomic modularity. Phage genomes are composed of independent functional modules that evolve separately and recombine in various configurations. While genomic modularity in phages has been extensively studied, less attention has been paid to protein modularity-proteins consisting of distinct building blocks that can evolve and recombine, enhancing functional and genetic diversity. Here, we use a set of 133,574 representative phage proteins and highly sensitive homology detection to capture instances of domain mosaicism, defined as fragment sharing between two otherwise unrelated proteins, and to understand its relationship with functional diversity in phage genomes. We discover that unrelated proteins from diverse functional classes frequently share homologous domains. This phenomenon is particularly pronounced within receptor-binding proteins, endolysins, and DNA polymerases. We also identify multiple instances of recent diversification via domain shuffling in receptor-binding proteins, neck passage structures, endolysins and some members of the core replication machinery, often transcending distant taxonomic and ecological boundaries. Our findings suggest that ongoing diversification via domain shuffling is reflective of a co-evolutionary arms race, driven by the need to overcome various bacterial resistance mechanisms against phages.

CRISPR-Cas9 genome editing in Steinernema entomopathogenic nematodes

Molecular tool development in traditionally non-tractable animals opens new avenues to study gene functions in the relevant ecological context. Entomopathogenic nematodes (EPN) Steinernema and their symbiotic bacteria of Xenorhabdus spp are a valuable experimental system in the laboratory and are applicable in the field to promote agricultural productivity. The infective juvenile (IJ) stage of the nematode packages mutualistic symbiotic bacteria in the intestinal pocket and invades insects that are agricultural pests. The lack of consistent and heritable genetics tools in EPN targeted mutagenesis severely restricted the study of molecular mechanisms underlying both parasitic and mutualistic interactions. Here, I report a protocol for CRISPR-Cas9 based genome-editing that is successful in two EPN species, S. carpocapsae and S. hermaphroditum . I adapted a gonadal microinjection technique in S. carpocapsae , which created on-target modifications of a homologue Sc-dpy-10 (cuticular collagen) by homology-directed repair. A similar delivery approach was used to introduce various alleles in S. hermaphroditum including Sh-dpy-10 and Sh-unc-22 (a muscle gene), resulting in visible and heritable phenotypes of dumpy and twitching, respectively. Using conditionally dominant alleles of Sh-unc-22 as a co-CRISPR marker, I successfully modified a second locus encoding Sh-Daf-22 (a homologue of human sterol carrier protein SCPx), predicted to function as a core enzyme in the biosynthesis of nematode pheromone that is required for IJ development. As a proof of concept, Sh-daf-22 null mutant showed IJ developmental defects in vivo ( in insecta) . This research demonstrates that Steinernema spp are highly tractable for targeted mutagenesis and has great potential in the study of gene functions under controlled laboratory conditions within the relevant context of its ecological niche.

Signatures of adaptive decreased virulence of deformed wing virus in an isolated population of wild honeybees (Apis mellifera)

Understanding the ecological and evolutionary processes that drive host-pathogen interactions is critical for combating epidemics and conserving species. The Varroa destructor mite and deformed wing virus (DWV) are two synergistic threats to Western honeybee (Apis mellifera) populations across the globe. Distinct honeybee populations have been found to self-sustain despite Varroa infestations, including colonies within the Arnot Forest outside Ithaca, NY, USA. We hypothesized that in these bee populations, DWV has been selected to produce an avirulent infection phenotype, allowing for the persistence of both host and disease-causing agents. To investigate this, we assessed the titre of viruses in bees from the Arnot Forest and managed apiaries, and assessed genomic variation and virulence differences between DWV isolates. Across groups, we found viral abundance was similar, but DWV genotypes were distinct. We also found that infections with isolates from the Arnot Forest resulted in higher survival and lower rates of symptomatic deformed wings, compared to analogous isolates from managed colonies, providing preliminary evidence to support the hypothesis of adaptive decreased viral virulence. Overall, this multi-level investigation of virus genotype and phenotype indicates that host ecological context can be a significant driver of viral evolution and host-pathogen interactions in honeybees.

Adaptive evolution in virulence effectors of the rice blast fungus Pyricularia oryzae

Plant pathogens secrete proteins called effectors that target host cellular processes to promote disease. Recently, structural genomics has identified several families of fungal effectors that share a similar three-dimensional structure despite remarkably variable amino-acid sequences and surface properties. To explore the selective forces that underlie the sequence variability of structurally-analogous effectors, we focused on MAX effectors, a structural family of effectors that are major determinants of virulence in the rice blast fungus Pyricularia oryzae. Using structure-informed gene annotation, we identified 58 to 78 MAX effector genes per genome in a set of 120 isolates representing seven host-associated lineages. The expression of MAX effector genes was primarily restricted to the early biotrophic phase of infection and strongly influenced by the host plant. Pangenome analyses of MAX effectors demonstrated extensive presence/absence polymorphism and identified gene loss events possibly involved in host range adaptation. However, gene knock-in experiments did not reveal a strong effect on virulence phenotypes suggesting that other evolutionary mechanisms are the main drivers of MAX effector losses. MAX effectors displayed high levels of standing variation and high rates of non-synonymous substitutions, pointing to widespread positive selection shaping the molecular diversity of MAX effectors. The combination of these analyses with structural data revealed that positive selection acts mostly on residues located in particular structural elements and at specific positions. By providing a comprehensive catalog of amino acid polymorphism, and by identifying the structural determinants of the sequence diversity, our work will inform future studies aimed at elucidating the function and mode of action of MAX effectors.

Redescription and new record of Paracapillaria (Ophidiocapillaria) najae (Nematoda: Trichuroidea) in the monocled cobra Naja kaouthia from central Thailand: morphological and molecular insights

The parasitic nematode Paracapillaria (Ophidiocapillaria) najae De, 1998, found in the Indian cobra Naja naja is redescribed and re-illustrated in the present study. The monocled cobra Naja kaouthia was discovered to be a new host for this parasite in central Thailand. A comprehensive description extending the morphological and molecular characteristics of the parasites is provided to aid species recognition in future studies. The morphometric characters of 41 parasites collected from 5 cobra specimens are compared with those described in the original studies. Phylogenetic analyses using mitochondrial cytochrome c oxidase subunit 1 and nuclear 18S ribosomal RNA genes were performed to provide novel information on the systematics of P. najae. Similar characteristics were observed in the examined nematode samples, despite being found in different hosts, confirming their identity as P. najae. The molecular genetic results support the species status of P. najae, indicating P. najae is well defined and separated from other related nematode species in the family Capillariidae. Morphological descriptions, genetic sequences, evolutionary relationships among capillariids and new host and distribution records of P. najae are discussed. Paracapillaria najae specimens found in the Thai cobra had some morphological variation, and sexual size dimorphism was also indicated. Paracapillaria najae was found to infect various cobra host species and appeared to be common throughout the Oriental regions, consistent with its hosts' distribution.

High functional allelic diversity and copy number in both MHC classes in the common buzzard

BACKGROUND

The major histocompatibility complex (MHC), which encodes molecules that recognize various pathogens and parasites and initiates the adaptive immune response in vertebrates, is renowned for its exceptional polymorphism and is a model of adaptive gene evolution. In birds, the number of MHC genes and sequence diversity varies greatly among taxa, believed due to evolutionary history and differential selection pressures. Earlier characterization studies and recent comparative studies suggest that non-passerine species have relatively few MHC gene copies compared to passerines. Additionally, comparative studies that have looked at partial MHC sequences have speculated that non-passerines have opposite patterns of selection on MHC class I (MHC-I) and class II (MHC-II) loci than passerines: namely, greater sequence diversity and signals of selection on MHC-II than MHC-I. However, new sequencing technology is revealing much greater MHC variation than previously expected while also facilitating full sequence variant detection directly from genomic data. Our study aims to take advantage of high-throughput sequencing methods to fully characterize both classes and domains of MHC of a non-passerine bird of prey, the common buzzard (Buteo buteo), to test predictions of MHC variation and differential selection on MHC classes.

RESULTS

Using genetic, genomic, and transcriptomic high-throughput sequencing data, we established common buzzards have at least three loci that produce functional alleles at both MHC classes. In total, we characterize 91 alleles from 113 common buzzard chicks for MHC-I exon 3 and 41 alleles from 125 chicks for MHC-IIB exon 2. Among these alleles, we found greater sequence polymorphism and stronger diversifying selection at MHC-IIB exon 2 than MHC-I exon 3, suggesting differential selection pressures on MHC classes. However, upon further investigation of the entire peptide-binding groove by including genomic data from MHC-I exon 2 and MHC-IIA exon 2, this turned out to be false. MHC-I exon 2 was as polymorphic as MHC-IIB exon 2 and MHC-IIA exon 2 was essentially invariant. Thus, comparisons between MHC-I and MHC-II that included both domains of the peptide-binding groove showed no differences in polymorphism nor diversifying selection between the classes. Nevertheless, selection analysis indicates balancing selection has been acting on common buzzard MHC and phylogenetic inference revealed that trans-species polymorphism is present between common buzzards and species separated for over 33 million years for class I and class II.

CONCLUSIONS

We characterize and confirm the functionality of unexpectedly high copy number and allelic diversity in both MHC classes of a bird of prey. While balancing selection is acting on both classes, there is no evidence of differential selection pressure on MHC classes in common buzzards and this result may hold more generally once more data for understudied MHC exons becomes available.

Evolution of MHC IIB Diversity Across Cichlid Fish Radiations

The genes of the major histocompatibility complex (MHC) are among the most polymorphic genes in vertebrates and crucial for their adaptive immune response. These genes frequently show inconsistencies between allelic genealogies and species phylogenies. This phenomenon is thought to be the result of parasite-mediated balancing selection maintaining ancient alleles through speciation events (trans-species polymorphism [TSP]). However, allele similarities may also arise from postspeciation mechanisms, such as convergence or introgression. Here, we investigated the evolution of MHC class IIB diversity in the cichlid fish radiations across Africa and the Neotropics by a comprehensive review of available MHC IIB DNA sequence information. We explored what mechanism explains the MHC allele similarities found among cichlid radiations. Our results showed extensive allele similarity among cichlid fish across continents, likely due to TSP. Functionality at MHC was also shared among species of the different continents. The maintenance of MHC alleles for long evolutionary times and their shared functionality may imply that certain MHC variants are essential in immune adaptation, even in species that diverged millions of years ago and occupy different environments.

Seascape genomics reveals limited dispersal and suggests spatially varying selection among European populations of sea lamprey (Petromyzon marinus)

Sea lamprey Petromyzon marinus is an anadromous and semelparous fish without homing behaviors. Despite being a freshwater, free-living organism for a large part of their life cycle, its adulthood is spent as a parasite of marine vertebrates. In their native European range, while it is well-established that sea lampreys comprise a single nearly-panmictic population, few studies have further explored the evolutionary history of natural populations. Here, we performed the first genome-wide characterization of sea lamprey's genetic diversity in their European natural range. The objectives were to investigate the connectivity among river basins and explore evolutionary processes mediating dispersal during the marine phase, with the sequencing of 186 individuals from 8 locations spanning the North Eastern Atlantic coast and the North Sea with double-digest RAD-sequencing, obtaining a total of 30,910 bi-allelic SNPs. Population genetic analyses reinforced the existence of a single metapopulation encompassing freshwater spawning sites within the North Eastern Atlantic and the North Sea, though the prevalence of private alleles at northern latitudes suggested some limits to the species' dispersal. Seascape genomics suggested a scenario where oxygen concentration and river runoffs impose spatially varying selection across their distribution range. Exploring associations with the abundance of potential hosts further suggested that hake and cod could also impose selective pressures, although the nature of such putative biotic interactions was unresolved. Overall, the identification of adaptive seascapes in a panmictic anadromous species could contribute to conservation practices by providing information for restoration activities to mitigate local extinctions on freshwater sites.

Understanding the evolution of immune genes in jawed vertebrates

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

Inspecting the interaction between human immunodeficiency virus and the immune system through genetic turnover

Chronic infections of the human immunodeficiency virus (HIV) create a very complex coevolutionary process, where the virus tries to escape the continuously adapting host immune system. Quantitative details of this process are largely unknown and could help in disease treatment and vaccine development. Here we study a longitudinal dataset of ten HIV-infected people, where both the B-cell receptors and the virus are deeply sequenced. We focus on simple measures of turnover, which quantify how much the composition of the viral strains and the immune repertoire change between time points. At the single-patient level, the viral-host turnover rates do not show any statistically significant correlation, however, they correlate if one increases the amount of statistics by aggregating the information across patients. We identify an anti-correlation: large changes in the viral pool composition come with small changes in the B-cell receptor repertoire. This result seems to contradict the naïve expectation that when the virus mutates quickly, the immune repertoire needs to change to keep up. However, a simple model of antagonistically evolving populations can explain this signal. If it is sampled at intervals comparable with the sweep time, one population has had time to sweep while the second cannot start a counter-sweep, leading to the observed anti-correlation. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Cryptosporidium uses multiple distinct secretory organelles to interact with and modify its host cell

Cryptosporidium is a leading cause of diarrheal disease in children and an important contributor to early childhood mortality. The parasite invades and extensively remodels intestinal epithelial cells, building an elaborate interface structure. How this occurs at the molecular level and the contributing parasite factors are largely unknown. Here, we generated a whole-cell spatial proteome of the Cryptosporidium sporozoite and used genetic and cell biological experimentation to discover the Cryptosporidium-secreted effector proteome. These findings reveal multiple organelles, including an original secretory organelle, and generate numerous compartment markers by tagging native gene loci. We show that secreted proteins are delivered to the parasite-host interface, where they assemble into different structures including a ring that anchors the parasite into its unique epicellular niche. Cryptosporidium thus uses a complex set of secretion systems during and following invasion that act in concert to subjugate its host cell.

Population genetic analysis of the microsporidium Ordospora colligata reveals the role of natural selection and phylogeography on its extremely compact and reduced genome

The determinants of variation in a species' genome-wide nucleotide diversity include historical, environmental, and stochastic aspects. This diversity can inform us about the species' past and present evolutionary dynamics. In parasites, the mode of transmission and the interactions with the host might supersede the effects of these aspects in shaping parasite genomic diversity. We used genomic samples from 10 populations of the microsporidian parasite Ordospora colligata to investigate present genomic diversity and how it was shaped by evolutionary processes, specifically, the role of phylogeography, co-phylogeography (with the host), natural selection, and transmission mode. Although very closely related microsporidia cause diseases in humans, O. colligata is specific to the freshwater crustacean Daphnia magna and has one of the smallest known eukaryotic genomes. We found an overlapping phylogeography between O. colligata and its host highlighting the long-term, intimate relationship between them. The observed geographic distribution reflects previous findings that O. colligata exhibits adaptations to colder habitats, which differentiates it from other microsporidian gut parasites of D. magna predominantly found in warmer areas. The co-phylogeography allowed us to calibrate the O. colligata phylogeny and thus estimate its mutation rate. We identified several genetic regions under potential selection. Our whole-genome study provides insights into the evolution of one of the most reduced eukaryotic genomes and shows how different processes shape genomic diversity of an obligate parasite.

The enemy within: lipid asymmetry in intracellular parasite-host interactions

Eukaryotic pathogens with an intracellular parasitic lifestyle are shielded from extracellular threats during replication and growth. In addition to many nutrients, parasites scavenge host cell lipids to establish complex membrane structures inside their host cells. To counteract the disturbance of the host cell plasma membrane they have evolved strategies to regulate phospholipid asymmetry. In this review, the function and importance of lipid asymmetry in the interactions of intracellular protozoan parasites with the target and immune cells of the host are highlighted. The malaria parasite Plasmodium infects red blood cells and extensively refurbishes these terminally differentiated cells. Cholesterol depletion and an altered intracellular calcium ion homeostasis can lead to disruption in erythrocyte membrane asymmetry and increased exposure of phosphatidylserine (PS). Binding to the PS receptor on monocytes and macrophages results in phagocytosis and destruction of infected erythrocytes. Leishmania parasites display apoptotic mimicry by actively enhancing PS exposure on their surface to trigger increased infection of macrophages. In extracellular Toxoplasma gondii a P4-type ATPase/CDC50 co-chaperone pair functions as a flippase important for exocytosis of specialised secretory organelles. Identification and functional analysis of parasite lipid-translocating proteins, i.e. flippases, floppases, and scramblases, will be central for the recognition of the molecular mechanisms of parasite/host interactions. Ultimately, a better understanding of parasitic diseases, host immunity, and immune escape by parasites require more research on the dynamics of phospholipid bilayers of parasites and the infected host cell.

Energy-harnessing problem solving of primordial life: Modeling the emergence of catalytic host-nested parasite life cycles

All life forms on earth ultimately descended from a primordial population dubbed the last universal common ancestor or LUCA via Darwinian evolution. Extant living systems share two salient functional features, a metabolism extracting and transforming energy required for survival, and an evolvable, informational polymer-the genome-conferring heredity. Genome replication invariably generates essential and ubiquitous genetic parasites. Here we model the energetic, replicative conditions of LUCA-like organisms and their parasites, as well as adaptive problem solving of host-parasite pairs. We show using an adapted Lotka-Volterra frame-work that three host-parasite pairs-individually a unit of a host and a parasite that is itself parasitized, therefore a nested parasite pair-are sufficient for robust and stable homeostasis, forming a life cycle. This nested parasitism model includes competition and habitat restriction. Its catalytic life cycle efficiently captures, channels and transforms energy, enabling dynamic host survival and adaptation. We propose a Malthusian fitness model for a quasispecies evolving through a host-nested parasite life cycle with two core features, rapid replacement of degenerate parasites and increasing evolutionary stability of host-nested parasite units from one to three pairs.

Gut microbiota of endangered Australian sea lion pups is unchanged by topical ivermectin treatment for endemic hookworm infection

The gut microbiota is essential for the development and maintenance of the hosts' immune system. Disturbances to the gut microbiota in early life stages can result in long-lasting impacts on host health. This study aimed to determine if topical ivermectin treatment for endemic hookworm (Uncinaria sanguinis) infection in endangered Australian sea lion (Neophoca cinerea) pups resulted in gut microbial changes. The gut microbiota was characterised for untreated (control) (n = 23) and treated (n = 23) Australian sea lion pups sampled during the 2019 and 2020/21 breeding seasons at Seal Bay, Kangaroo Island. Samples were collected pre- and post-treatment on up to four occasions over a four-to-five-month period. The gut microbiota of untreated (control) and treated pups in both seasons was dominated by five bacterial phyla, Fusobacteria, Firmicutes, Proteobacteria, Actinobacteria and Bacteroidetes. A significant difference in alpha diversity between treatment groups was seen in pups sampled during the 2020/21 breeding season (p = 0.008), with higher richness and diversity in treated pups. Modelling the impact of individual pup identification (ID), capture, pup weight (kg), standard length (cm), age and sex on beta diversity revealed that pup ID accounted for most of the variation (35% in 2019 and 42% in 2020/21), with pup ID, capture, and age being the only significant contributors to microbial variation (p < 0.05). There were no statistically significant differences in the composition of the microbiota between treatment groups in both the 2019 and 2020/21 breeding seasons, indicating that topical ivermectin treatment did not alter the composition of the gut microbiota. To our knowledge, this is the first study to characterise the gut microbiota of free-ranging Australian pinniped pups, compare the composition across multiple time points, and to consider the impact of parasitic treatment on overall diversity and microbial composition of the gut microbiota. Importantly, the lack of compositional changes in the gut microbiota with treatment support the utility of topical ivermectin as a safe and minimally invasive management strategy to enhance pup survival in this endangered species.

Could aging evolve as a pathogen control strategy?

Aging is often attributed to the detrimental side effects of beneficial traits but not a programmed adaptive process. Alternatively, the pathogen control hypothesis posits that defense against infectious diseases may provide a strong selection force for restriction of lifespan. Aging might have evolved to remove older individuals who carry chronic diseases that may transmit to their younger kin. Thus, selection for shorter lifespans may benefit kin's fitness. The pathogen control hypothesis addresses arguments typically raised against adaptive aging concepts: it explains the benefit of shorter lifespan and the absence of mutant variants that do not age. We discuss the consistency and explanatory power of this hypothesis and compare it with classic hypotheses of aging.

Integrating plant and fungal quantitative genetics to improve the ecological and agricultural applications of mycorrhizal symbioses

Finding and targeting genes that quantitatively contribute to agricultural and ecological processes progresses food production and conservation efforts. Typically, quantitative genetic approaches link variants in a single organism's genome with a trait of interest. Recently, genome-to-genome mapping has found genome variants interacting between species to produce the result of a multiorganism (including multikingdom) interaction. These were plant and bacterial pathogen genome interactions; plant-fungal coquantitative genetics have not yet been applied. Plant-mycorrhizae symbioses exist across most biomes, for a majority of land plants, including crop plants, and manipulate many traits from single organisms to ecosystems for which knowing the genetic basis would be useful. The availability of Rhizophagus irregularis mycorrhizal isolates, with genomic information, makes dual-genome methods with beneficial mutualists accessible and imminent.

Parasitic Infection and Immunity-A Special Biomedicines Issue

Infectious parasitic diseases that affect humans and animals remain a central health concern worldwide [...].

Spillover and genome selection of the gut parasite Nosema ceranae between honey bee species

Nosema ceranae is a honey bee gut parasite that has recently spilled to another honey bee host through trading. The impact of infection on the native host is minor, which is substantial in the novel host. In this study, artificial inoculation simulated the parasite transmission from the native to the novel host. We found that the parasite initiated proliferation earlier in the novel host than in the native host. Additionally, parasite gene expression was significantly higher when infecting the novel host compared with the native host, leading to a significantly higher number of spores. Allele frequencies were similar for spores of parasites infecting both native and novel hosts. This suggests that the high number of spores found in the novel host was not caused by a subset of more fit spores from native hosts. Native hosts also showed a higher number of up-regulated genes in response to infection when compared with novel hosts. Our data further showed that native hosts suppressed parasite gene expression and arguably sacrificed cells to limit the parasite. The results provide novel insights into host defenses and gene selection during a parasite spillover event.

Evolutionary aspects of the parasitoid life strategy, with a particular emphasis on fly–spider interactions

Parasitoids are significant ecological elements of terrestrial food webs and have evolved within seven insect orders. Interestingly, however, associations with spiders as hosts have evolved only in two insect orders, Diptera and Hymenoptera. Here, we summarize various aspects of host utilization by dipteran flies with an emphasis on associations with spiders. Our synthesis reveals that spider flies (family Acroceridae) have evolved a unique life strategy among all the parasitoid taxa associated with spiders, in which koinobiont small-headed flies utilize an indirect oviposition strategy. This indirect oviposition in spider flies is inherited from Nemestrinimorpha ancestors which appeared in the Late Triassic and is characterized by the evolution of planidial larvae. Further, we discuss the advantages and disadvantages of indirect oviposition in spider flies. On the one hand, indirect oviposition allows the fly to avoid contact/wrestling with spider hosts. On the other hand, larval survival is low because the planidium must actively seek out and infect a suitable host individually. The risk of failure to find a suitable spider host is offset by the fly’s extremely high fecundity.

Evolutionary gain and loss of a pathological immune response to parasitism

Parasites impose fitness costs on their hosts. Biologists often assume that natural selection favors infection-resistant hosts. Yet, when the immune response itself is costly, theory suggests that selection may sometimes favor loss of resistance, which may result in alternative stable states where some populations are resistant and others are tolerant. Intraspecific variation in immune costs is rarely surveyed in a manner that tests evolutionary patterns, and there are few examples of adaptive loss of resistance. Here, we show that when marine threespine stickleback colonized freshwater lakes, they gained resistance to the freshwater-associated cestode Schistocephalus solidus. Extensive peritoneal fibrosis and inflammation are a commonly observed phenotype that contributes to suppression of cestode growth and viability but also imposes a substantial cost on fecundity. Combining genetic mapping and population genomics, we find that opposing selection generates immune system differences between tolerant and resistant populations, consistent with divergent optimization.

MHC class II genotype-by-pathogen genotype interaction for infection prevalence in a natural rodent-Borrelia system

MHC genes are extraordinarily polymorphic in most taxa. Host-pathogen coevolution driven by negative frequency-dependent selection (NFDS) is one of the main hypotheses for the maintenance of such immunogenetic variation. Here, we test a critical but rarely tested assumption of this hypothesis-that MHC alleles affect resistance/susceptibility to a pathogen in a strain-specific way, that is, there is a host genotype-by-pathogen genotype interaction. In a field study of bank voles naturally infected with the tick-transmitted bacterium Borrelia afzelii, we tested for MHC class II (DQB) genotype-by-B. afzelii strain interactions for infection prevalence between 10 DQB alleles and seven strains. One allele (DQB*37) showed an interaction, such that voles carrying DQB*37 had higher prevalence of two strains and lower prevalence of one strain than individuals without the allele. These findings were corroborated by analyses of strain composition of infections, which revealed an effect of DQB*37 in the form of lower β diversity among infections in voles carrying the allele. Taken together, these results provide rare support at the molecular genetic level for a key assumption of models of antagonistic coevolution through NFDS.

Host-pathogen coevolution drives innate immune response to Aphanomyces astaci infection in freshwater crayfish: transcriptomic evidence

BACKGROUND

For over a century, scientists have studied host-pathogen interactions between the crayfish plague disease agent Aphanomyces astaci and freshwater crayfish. It has been hypothesised that North American crayfish hosts are disease-resistant due to the long-lasting coevolution with the pathogen. Similarly, the increasing number of latent infections reported in the historically sensitive European crayfish hosts seems to indicate that similar coevolutionary processes are occurring between European crayfish and A. astaci. Our current understanding of these host-pathogen interactions is largely focused on the innate immunity processes in the crayfish haemolymph and cuticle, but the molecular basis of the observed disease-resistance and susceptibility remain unclear. To understand how coevolution is shaping the host's molecular response to the pathogen, susceptible native European noble crayfish and invasive disease-resistant marbled crayfish were challenged with two A. astaci strains of different origin: a haplogroup A strain (introduced to Europe at least 50 years ago, low virulence) and a haplogroup B strain (signal crayfish in lake Tahoe, USA, high virulence). Here, we compare the gene expression profiles of the hepatopancreas, an integrated organ of crayfish immunity and metabolism.

RESULTS

We characterised several novel innate immune-related gene groups in both crayfish species. Across all challenge groups, we detected 412 differentially expressed genes (DEGs) in the noble crayfish, and 257 DEGs in the marbled crayfish. In the noble crayfish, a clear immune response was detected to the haplogroup B strain, but not to the haplogroup A strain. In contrast, in the marbled crayfish we detected an immune response to the haplogroup A strain, but not to the haplogroup B strain.

CONCLUSIONS

We highlight the hepatopancreas as an important hub for the synthesis of immune molecules in the response to A. astaci. A clear distinction between the innate immune response in the marbled crayfish and the noble crayfish is the capability of the marbled crayfish to mobilise a higher variety of innate immune response effectors. With this study we outline that the type and strength of the host immune response to the pathogen is strongly influenced by the coevolutionary history of the crayfish with specific A. astaci strains.

Combining GWAS and population genomic analyses to characterize coevolution in a legume-rhizobia symbiosis

The mutualism between legumes and rhizobia is clearly the product of past coevolution. However, the nature of ongoing evolution between these partners is less clear. To characterize the nature of recent coevolution between legumes and rhizobia, we used population genomic analysis to characterize selection on functionally annotated symbiosis genes as well as on symbiosis gene candidates identified through a two-species association analysis. For the association analysis, we inoculated each of 202 accessions of the legume host Medicago truncatula with a community of 88 Sinorhizobia (Ensifer) meliloti strains. Multistrain inoculation, which better reflects the ecological reality of rhizobial selection in nature than single-strain inoculation, allows strains to compete for nodulation opportunities and host resources and for hosts to preferentially form nodules and provide resources to some strains. We found extensive host by symbiont, that is, genotype-by-genotype, effects on rhizobial fitness and some annotated rhizobial genes bear signatures of recent positive selection. However, neither genes responsible for this variation nor annotated host symbiosis genes are enriched for signatures of either positive or balancing selection. This result suggests that stabilizing selection dominates selection acting on symbiotic traits and that variation in these traits is under mutation-selection balance. Consistent with the lack of positive selection acting on host genes, we found that among-host variation in growth was similar whether plants were grown with rhizobia or N-fertilizer, suggesting that the symbiosis may not be a major driver of variation in plant growth in multistrain contexts.

Immunogenomics of Killer Cell Immunoglobulin-Like Receptor (KIR) and HLA Class I: Coevolution and Consequences for Human Health

Interactions of killer cell immunoglobin-like receptors (KIR) with human leukocyte antigens (HLA) class I regulate effector functions of key cytotoxic cells of innate and adaptive immunity. The extreme diversity of this interaction is genetically determined, having evolved in the ever-changing environment of pathogen exposure. Diversity of KIR and HLA genes is further facilitated by their independent segregation on separate chromosomes. That fetal implantation relies on many of the same types of immune cells as infection control places certain constraints on the evolution of KIR interactions with HLA. Consequently, specific inherited combinations of receptors and ligands may predispose to specific immune-mediated diseases, including autoimmunity. Combinatorial diversity of KIR and HLA class I can also differentiate success rates of immunotherapy directed to these diseases. Progress toward both etiopathology and predicting response to therapy is being achieved through detailed characterization of the extent and consequences of the combinatorial diversity of KIR and HLA. Achieving these goals is more tractable with the development of integrated analyses of molecular evolution, function, and pathology that will establish guidelines for understanding and managing risks. Here, we present what is known about the coevolution of KIR with HLA class I and the impact of their complexity on immune function and homeostasis.

Unraveling coevolutionary dynamics using ecological genomics

Coevolutionary interactions, from the delicate co-dependency in mutualistic interactions to the antagonistic relationship of hosts and parasites, are a ubiquitous driver of adaptation. Surprisingly, little is known about the genomic processes underlying coevolution in an ecological context. However, species comprise genetically differentiated populations that interact with temporally variable abiotic and biotic environments. We discuss the recent advances in coevolutionary theory and genomics as well as shortcomings, to identify coevolving genes that take into account this spatial and temporal variability of coevolution, and propose a practical guide to understand the dynamic of coevolution using an ecological genomics lens.

How does genetic architecture affect eco-evolutionary dynamics? A theoretical perspective

Recent studies have revealed the importance of feedbacks between contemporary rapid evolution (i.e. evolution that occurs through changes in allele frequencies) and ecological dynamics. Despite its inherent interdisciplinary nature, however, studies on eco-evolutionary feedbacks have been mostly ecological and tended to focus on adaptation at the phenotypic level without considering the genetic architecture of evolutionary processes. In empirical studies, researchers have often compared ecological dynamics when the focal species under selection has a single genotype with dynamics when it has multiple genotypes. In theoretical studies, common approaches are models of quantitative traits where mean trait values change adaptively along the fitness gradient and Mendelian traits with two alleles at a single locus. On the other hand, it is well known that genetic architecture can affect short-term evolutionary dynamics in population genetics. Indeed, recent theoretical studies have demonstrated that genetic architecture (e.g. the number of loci, linkage disequilibrium and ploidy) matters in eco-evolutionary dynamics (e.g. evolutionary rescue where rapid evolution prevents extinction and population cycles driven by (co)evolution). I propose that theoretical approaches will promote the synthesis of functional genomics and eco-evolutionary dynamics through models that combine population genetics and ecology as well as nonlinear time-series analyses using emerging big data.

This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.

Genetic characterization of potential venom resistance proteins in California ground squirrels ( Otospermophilus beecheyi ) using transcriptome analyses

Understanding the molecular basis of adaptations in coevolving species requires identifying the genes that underlie reciprocally selected phenotypes, such as those involved in venom in snakes and resistance to the venom in their prey. In this regard, California ground squirrels (CGS; Otospermophilus beecheyi) are eaten by northern Pacific rattlesnakes (Crotalus oreganus oreganus), but individual squirrels may still show substantial resistance to venom and survive bites. A recent study using proteomics identified venom interactive proteins (VIPs) in the blood serum of CGS. These VIPs represent possible resistance proteins, but the sequences of genes encoding them are unknown despite the value of such data to molecular studies of coevolution. To address this issue, we analyzed a de novo assembled transcriptome from CGS liver tissue-where many plasma proteins are synthesized-and other tissues from this species. We then examined VIP sequences in terms of three characteristics that identify them as possible resistance proteins: evidence for positive selection, high liver expression, and nonsynonymous variation across CGS populations. Based on these characteristics, we identified five VIPs (i.e., α-2-macroglobulin, α-1-antitrypsin-like protein GS55-LT, apolipoprotein A-II, hibernation-associated plasma protein HP-20, and hibernation-associated plasma protein HP-27) as the most likely candidates for resistance proteins among VIPs identified to date. Four of these proteins have been previously implicated in conferring resistance to the venom in mammals, validating our approach. When combined with the detailed information available for rattlesnake venom proteins, these results set the stage for future work focused on understanding coevolutionary interactions at the molecular level between these species.

Evolution: Mitochondrial lodgers can take over in hermaphroditic snails

Mitochondria - the cell's power stations - are inherited uniparentally via eggs, not sperm. In hermaphroditic plants, they sometimes prevent their hosts from making pollen (and sperm), causing cytoplasmic male sterility. New evidence from a hermaphroditic freshwater snail now documents cytoplasmic male sterility in animals.

Bacterial origins of human cell-autonomous innate immune mechanisms

The cell-autonomous innate immune system enables animal cells to resist viral infection. This system comprises an array of sensors that, after detecting viral molecules, activate the expression of antiviral proteins and the interferon response. The repertoire of immune sensors and antiviral proteins has long been considered to be derived from extensive evolutionary innovation in vertebrates, but new data challenge this dogma. Recent studies show that central components of the cell-autonomous innate immune system have ancient evolutionary roots in prokaryotic genes that protect bacteria from phages. These include the cyclic GMP-AMP synthase (cGAS)-stimulator of interferon genes (STING) pathway, Toll/IL-1 receptor (TIR) domain-containing pathogen receptors, the viperin family of antiviral proteins, SAMHD1-like nucleotide-depletion enzymes, gasdermin proteins and key components of the RNA interference pathway. This Perspective details current knowledge of the elements of antiviral immunity that are conserved from bacteria to humans, and presents possible evolutionary scenarios to explain the observed conservation.

Association of Mannose-Binding Lectin 2 Gene Polymorphism with Tuberculosis Based on Mycobacterium tuberculosis Lineages

Purpose

Polymorphisms in MBL2 may contribute to the susceptibility to tuberculosis. The aim of the present study was to determine the associations of the polymorphisms of five loci (rs1800450, rs1800451, rs7096206, rs7095891, and rs11003125) in the MBL2 gene with susceptibility to tuberculosis and specific lineages of Mycobacterium tuberculosis causing tuberculosis in the Uyghur population of Xinjiang, China.

Methods

From January 2019 to January 2020, we enrolled 170 Uyghur tuberculosis patients as the case group and 147 Uyghur staff with no clinical symptoms as the control group from four designated tuberculosis hospitals in southern Xinjiang, China. The polymorphisms of five loci in MBL2 of human were detected by sequencing. Whole-genome sequencing was applied in 68 M. tuberculosis isolates from the case group and the data were used to perform genealogy analysis.

Results

The distributions of allele and genotype frequencies of five loci in MBL2 varied little between the case and control groups and varied little among the groups, including those infected with different lineages of M. tuberculosis and the control (except those of rs11003125), the P values were all >0.05. The distribution of alleles of rs11003125 was statistically different between patients infected with lineages 3 and 4 M. tuberculosis (χ²=7.037, P=0.008). The C allele and CC genotype of rs11003125 were found to be protective factors against lineage 4 infection when compared to lineage 3 (ORs were 0.190 and 0.158, respectively; 95% confidence intervals were 0.053~0.690 and 0.025~0.999, respectively).

Conclusion

Our results suggested that human’s susceptibility to tuberculosis is affected both by the host genetic polymorphisms and the lineage of the M. tuberculosis that people were exposed to. However, due to the limitation of the sample size in the present study, larger sample size and more rigorous design should be guaranteed in future studies.

Expansion of Cyclophyllidea Biodiversity in Rodents of Qinghai-Tibet Plateau and the "Out of Qinghai-Tibet Plateau" Hypothesis of Cyclophyllideans

The Cyclophyllidea comprises the most species-rich order of tapeworms (Platyhelminthes, Cestoda) and includes species with some of the most severe health impact on wildlife, livestock, and humans. We collected seven Cyclophyllidea specimens from rodents in Qinghai-Tibet Plateau (QTP) and its surrounding mountain systems, of which four specimens in QTP were unsequenced, representing “putative new species.” Their complete mitochondrial (mt) genomes were sequenced and annotated. Phylogenetic reconstruction of partial 28S rDNA, cox1 and nad1 datasets provided high bootstrap frequency support for the categorization of three “putative new species,” assigning each, respectively, to the genera Mesocestoides, Paranoplocephala, and Mosgovoyia, and revealing that some species and families in these three datasets, which contain 291 species from nine families, may require taxonomic revision. The partial 18S rDNA phylogeny of 29 species from Taeniidae provided high bootstrap frequency support for the categorization of the “putative new species” in the genus Hydatigera. Combined with the current investigation, the other three known Taeniidae species found in this study were Taenia caixuepengi, T. crassiceps, and Versteria mustelae and may be widely distributed in western China. Estimates of divergence time based on cox1 + nad1 fragment and mt protein-coding genes (PCGs) showed that the differentiation rate of Cyclophyllidea species was strongly associated with the rate of change in the biogeographic scenarios, likely caused by the uplift of the QTP; i.e., species differentiation of Cyclophyllidea might be driven by host-parasite co-evolution caused by the uplift of QTP. We propose an “out of QTP” hypothesis for the radiation of these cyclophyllidean tapeworms.

A myxozoan genome reveals mosaic evolution in a parasitic cnidarian

Background

Parasite evolution has been conceptualized as a process of genetic loss and simplification. Contrary to this model, there is evidence of expansion and conservation of gene families related to essential functions of parasitism in some parasite genomes, reminiscent of widespread mosaic evolution—where subregions of a genome have different rates of evolutionary change. We found evidence of mosaic genome evolution in the cnidarian Myxobolus honghuensis, a myxozoan parasite of fish, with extremely simple morphology.

Results

We compared M. honghuensis with other myxozoans and free-living cnidarians, and determined that it has a relatively larger myxozoan genome (206 Mb), which is less reduced and less compact due to gene retention, large introns, transposon insertion, but not polyploidy. Relative to other metazoans, the M. honghuensis genome is depleted of neural genes and has only the simplest animal immune components. Conversely, it has relatively more genes involved in stress resistance, tissue invasion, energy metabolism, and cellular processes compared to other myxozoans and free-living cnidarians. We postulate that the expansion of these gene families is the result of evolutionary adaptations to endoparasitism. M. honghuensis retains genes found in free-living Cnidaria, including a reduced nervous system, myogenic components, ANTP class Homeobox genes, and components of the Wnt and Hedgehog pathways.

Conclusions

Our analyses suggest that the M. honghuensis genome evolved as a mosaic of conservative, divergent, depleted, and enhanced genes and pathways. These findings illustrate that myxozoans are not as genetically simple as previously regarded, and the evolution of some myxozoans is driven by both genomic streamlining and expansion.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01249-8.

Epigenetics and Helicobacter pylori

Epigenetics regulates gene expression, cell type development during differentiation, and the cell response to environmental stimuli. To survive, bacteria need to evade the host immune response. Bacteria, including Helicobacter pylori (Hp), reach this target epigenetically, altering the chromatin of the host cells, in addition to several more approaches, such as DNA mutation and recombination. This review shows that Hp prevalently silences the genes of the human gastric mucosa by DNA methylation. Epigenetics includes different mechanisms. However, DNA methylation persists after DNA replication and therefore is frequently associated with the inheritance of repressed genes. Chromatin modification can be transmitted to daughter cells leading to heritable changes in gene expression. Aberrant epigenetic alteration of the gastric mucosa DNA remains the principal cause of gastric cancer. Numerous methylated genes have been found in cancer as well as in precancerous lesions of Hp-infected patients. These methylated genes inactivate tumor-suppressor genes. It is time for us to complain about our genetic and epigenetic makeups for our diseases.

Leapfrog dynamics in phage‐bacteria coevolution revealed by joint analysis of cross‐infection phenotypes and whole genome sequencing

Viruses and their hosts can undergo coevolutionary arms races where hosts evolve increased resistance and viruses evolve counter-resistance. Given these arms race dynamics (ARD), both players are predicted to evolve along a single trajectory as more recently evolved genotypes replace their predecessors. By coupling phenotypic and genomic analyses of coevolving populations of bacteriophage λ and Escherichia coli, we find conflicting evidence for ARD. Virus-host infection phenotypes fit the ARD model, yet genomic analyses revealed fluctuating selection dynamics. Rather than coevolution unfolding along a single trajectory, cryptic genetic variation emerges and is maintained at low frequency for generations until it eventually supplants dominant lineages. These observations suggest a hybrid 'leapfrog' dynamic, revealing weaknesses in the predictive power of standard coevolutionary models. The findings shed light on the mechanisms that structure coevolving ecological networks and reveal the limits of using phenotypic or genomic data alone to differentiate coevolutionary dynamics.