Parallel adaptation of rabbit populations to myxoma virus

Copy number variation contributes to parallel local adaptation in an invasive plant

Adaptation is a critical determinant of the diversification, persistence, and geographic range limits of species. Yet the genetic basis of adaptation is often unknown and potentially underpinned by a wide range of mutational types-from single nucleotide changes to large-scale alterations of chromosome structure. Copy number variation (CNV) is thought to be an important source of adaptive genetic variation, as indicated by decades of candidate gene studies that point to CNVs underlying rapid adaptation to strong selective pressures. Nevertheless, population-genomic studies of CNVs face unique logistical challenges not encountered by other forms of genetic variation. Consequently, few studies have systematically investigated the contributions of CNVs to adaptation at a genome-wide scale. We present a genome-wide analysis of CNV contributing to the adaptation of an invasive weed, Ambrosia artemisiifolia. CNVs show clear signatures of parallel local adaptation between North American (native) and European (invaded) ranges, implying widespread reuse of CNVs during adaptation to shared heterogeneous patterns of selection. We used a local principal component analysis (PCA) to genotype CNV regions in whole-genome sequences of samples collected over the last two centuries. We identified 16 large CNV regions of up to 11.85 megabases in length, eight of which show signals of rapid evolutionary change, with pronounced frequency shifts between historic and modern populations. Our results provide compelling genome-wide evidence that CNV underlies rapid adaptation over contemporary timescales of natural populations.

An interdisciplinary approach to improving conservation outcomes for parasites

Parasites represent a significant proportion of Earth's biodiversity and play important roles in the ecology and biology of ecosystems and hosts, making them an important target for conservation. Despite increasing calls to prioritize protection for parasites in the academic literature, they remain undervalued and underrepresented in global biodiversity conservation efforts, not least due to the perception that the interests of parasite and host conservation are opposing and the common misconception that parasites are a threat, rather than a benefit, to conservation. We considered whether taking an interdisciplinary approach to parasite conservation research will generate novel insights and solutions concerning why and how parasite conservation should be practiced for the benefit of parasites, their hosts, ecosystems, and people. We argue that 2 of the main barriers to more widespread parasite conservation are the knowledge gap concerning the role of sociocultural factors affecting the willingness to enact parasite conservation and the lack of a consistent and cohesive philosophical basis for parasite conservation. Possible sociocultural barriers to parasite conservation include misconceptions of the risks posed by parasites, taxonomic bias, differences in conservation values, economic constraints, and technical challenges. The use of social science can generate insights into levels of awareness and support for parasite conservation and improve understanding of how human values and attitudes mediate conservation practices concerning parasites. Such knowledge will have a critical role in addressing sociocultural barriers and improving support for parasite conservation. Issues with the current philosophical basis for parasite conservation include contradictory accounts of which parasites merit conservation, insufficient explanation of how different conservation values apply to parasite biodiversity, and the existence of a false antagonism between host and parasite conservation. Greater engagement with philosophical work on environmental ethics and biological unitization will strengthen existing arguments for parasite conservation and will support conservation decision-making processes.

Towards a unified approach in managing resistance to vaccines, drugs, and pesticides

Everywhere, pests and pathogens evolve resistance to our control efforts, impairing human health and welfare. Developing sustainable solutions to this problem requires working with evolved immune and ecological systems, rather than against these evolutionary forces. We advocate a transdisciplinary approach to resistance based on an evolutionary foundation informed by the concepts of integrated pest management and One Health. Diverse, multimodal management approaches create a more challenging environment for the evolution of resistance. Given our permanent evolutionary and ecological relationships with pests and pathogens, responses to most biological threats to health and agriculture should seek sustainable harm reduction rather than eradication.

Estimating the impact of sarcoptic mange epidemic on the population size of wild raccoon dogs (Nyctereutes procyonoides) from wildlife rescue data

The impact of infectious diseases on host populations is often not quantified because it is difficult to observe the host population and infectious disease dynamics. To address this problem, we developed a state-space model to simultaneously estimate host population and disease dynamics using wildlife rescue data. Using this model, we aimed to quantify the impact of sarcoptic mange on a Japanese raccoon dog population by estimating the change in their relative population size. We classified the status of rescued raccoon dogs into four categories: i) rescued due to infection with mange, ii) rescued due to traffic accidents without mange, iii) rescued due to traffic accidents with mange, and iv) rescued due to causes other than traffic accidents or mange. We modelled the observation process for each category and fitted the model to the reported number of raccoon dogs rescued between 1990 and 2010 at three wildlife rescue facilities in Kanagawa Prefecture, Japan. The mortality rate induced by mange was estimated to be 1.09 (95% credible interval (CI): 0.47-1.72) per year. The estimated prevalence of sarcoptic mange ranged between 4 and 80% in the study period. When a substantial prevalence of mange was observed (1995-2002), the host population size decreased by 91.2% (95% credible intervals: 86.3-94.7). We show that the impact of infectious disease outbreak on the wildlife population can be estimated from the time-series data of wildlife rescue events due to multiple causes. Our estimates suggest that sarcoptic mange triggered a substantial decrease in the Japanese wild raccoon dog populations.

Linking the spatial and genomic structure of adaptive potential for conservation management: a review

We unified the recent literature with the goal to contribute to the discussion on how genetic diversity might best be conserved. We argue that this decision will be guided by how genomic variation is distributed among manageable populations (i.e., its spatial structure), the degree to which adaptive potential is best predicted by variation across the entire genome or the subset of that variation that is identified as putatively adaptive (i.e., its genomic structure), and whether we are managing species as single entities or as collections of diversifying lineages. The distribution of genetic variation and our ultimate goal will have practical implications for on-the-ground management. If adaptive variation is largely polygenic or responsive to change, its spatial structure might be broadly governed by the forces determining genome-wide variation (linked selection, drift, and gene flow), making measurement and prioritization straightforward. If we are managing species as single entities, then population-level prioritization schemes are possible so as to maximize future pooled genetic variation. We outline one such scheme based on the popular Shapley value from cooperative game theory that considers the relative genetic contribution of a population to an unknown future collection of populations.

Genetic Independence of Naturally Correlated Variation in Resistance to Endemic and Novel Pathogens

The emergence of new diseases is an urgent concern, but hosts can also vary in resistance to pathogens that are novel to them, facilitating evolutionary rescue. However, little is known about the genetic source for polymorphic resistance to novel pathogens or its relationship to defences against endemic diseases. With anther-smut disease from wild plant populations, we used selection experiments and genetic analyses to show that resistances to novel and endemic pathogens are genetically independent, despite being positively correlated in nature. Moreover, novel-pathogen resistance presented a much simpler genetic basis and more rapid response to selection. We demonstrate that polymorphic resistance to a newly introduced disease is genetically determined and not an extension of defences against the related endemic pathogen, challenging the conventional view of nonhost resistance.

Impact of putatively beneficial genomic loci on gene expression in little brown bats (Myotis lucifugus, Le Conte, 1831) affected by white-nose syndrome

Genome-wide scans for selection have become a popular tool for investigating evolutionary responses in wildlife to emerging diseases. However, genome scans are susceptible to false positives and do little to demonstrate specific mechanisms by which loci impact survival. Linking putatively resistant genotypes to observable phenotypes increases confidence in genome scan results and provides evidence of survival mechanisms that can guide conservation and management efforts. Here we used an expression quantitative trait loci (eQTL) analysis to uncover relationships between gene expression and alleles associated with the survival of little brown bats (Myotis lucifugus) despite infection with the causative agent of white-nose syndrome. We found that 25 of the 63 single-nucleotide polymorphisms (SNPs) associated with survival were related to gene expression in wing tissue. The differentially expressed genes have functional annotations associated with the innate immune system, metabolism, circadian rhythms, and the cellular response to stress. In addition, we observed differential expression of multiple genes with survival implications related to loci in linkage disequilibrium with focal SNPs. Together, these findings support the selective function of these loci and suggest that part of the mechanism driving survival may be the alteration of immune and other responses in epithelial tissue.

Practical Suggestions for Assessing Rabbit Haemorrhagic Disease Virus 2 Risk to Endangered Native Lagomorphs in North America and Southern Africa

A new form of the rabbit haemorrhagic disease virus, RHDV2, first observed in European rabbits, has spread widely among different species of hares in Europe, jackrabbits and cottontails in North America, and hares in southern Africa. However, only limited surveillance studies have been undertaken so far. It is suggested that methods developed for controlling the disease in farmed rabbits in Europe and studying the efficacy of RHDV as a biological control agent in Australia could facilitate epidemiological research on those recently affected lagomorph species. This would enable the assessment of the risk of RHDV2 to native lagomorphs, including endangered species, and the determination of the main host species of RHDV2. Because RHDV2 has not spread equally through all lagomorph species, epidemiological studies could give insights into factors important for determining host susceptibility.

Drosophila melanogaster pigmentation demonstrates adaptive phenotypic parallelism but genomic unpredictability over multiple timescales

Populations are capable of responding to environmental change over ecological timescales via adaptive tracking. However, the translation from patterns of allele frequency change to rapid adaptation of complex traits remains unresolved. We used abdominal pigmentation in Drosophila melanogaster as a model phenotype to address the nature, genetic architecture, and repeatability of rapid adaptation in the field. We show that D. melanogaster pigmentation evolves as a highly parallel and deterministic response to shared environmental gradients across latitude and season in natural North American populations. We then experimentally evolved replicate, genetically diverse fly populations in field mesocosms to remove any confounding effects of demography and/or cryptic structure that may drive patterns in wild populations; we show that pigmentation rapidly responds, in parallel, in fewer than ten generations. Thus, pigmentation evolves concordantly in response to spatial and temporal climatic gradients. We next examined whether phenotypic differentiation was associated with allele frequency change at loci with established links to genetic variance in pigmentation in natural populations. We found that across all spatial and temporal scales, phenotypic patterns were associated with variation at pigmentation-related loci, and the sets of genes we identified in each context were largely nonoverlapping. Therefore, our findings suggest that parallel phenotypic evolution is associated with an unpredictable genomic response, with distinct components of the polygenic architecture shifting across each environmental gradient to produce redundant adaptive patterns.

Significance Statement

Shifts in global climate conditions have heightened our need to understand the dynamics and pace of adaptation in natural populations. In order to anticipate the population-level response to rapidly changing environmental conditions, we need to understand whether trait evolution is predictable over short timescales, and whether the genetic basis of adaptation is shared or distinct across multiple timescales. Here, we explored parallelism in the adaptive response of a complex phenotype, D. melanogaster pigmentation, to shared conditions that varied over multiple spatiotemporal scales. Our results demonstrate that while phenotypic adaptation proceeds as a predictable response to environmental gradients, even over short timescales, the genetic basis of the adaptive response is variable and nuanced across spatial and temporal contexts.

Intergenomic signatures of coevolution between Tasmanian devils and an infectious cancer

Coevolution is common and frequently governs host-pathogen interaction outcomes. Phenotypes underlying these interactions often manifest as the combined products of the genomes of interacting species, yet traditional quantitative trait mapping approaches ignore these intergenomic interactions. Devil facial tumor disease (DFTD), an infectious cancer afflicting Tasmanian devils (Sarcophilus harrisii), has decimated devil populations due to universal host susceptibility and a fatality rate approaching 100%. Here, we used a recently developed joint genome-wide association study (i.e., co-GWAS) approach, 15 y of mark-recapture data, and 960 genomes to identify intergenomic signatures of coevolution between devils and DFTD. Using a traditional GWA approach, we found that both devil and DFTD genomes explained a substantial proportion of variance in how quickly susceptible devils became infected, although genomic architectures differed across devils and DFTD; the devil genome had fewer loci of large effect whereas the DFTD genome had a more polygenic architecture. Using a co-GWA approach, devil-DFTD intergenomic interactions explained ~3× more variation in how quickly susceptible devils became infected than either genome alone, and the top genotype-by-genotype interactions were significantly enriched for cancer genes and signatures of selection. A devil regulatory mutation was associated with differential expression of a candidate cancer gene and showed putative allele matching effects with two DFTD coding sequence variants. Our results highlight the need to account for intergenomic interactions when investigating host-pathogen (co)evolution and emphasize the importance of such interactions when considering devil management strategies.

First Detection of Benign Rabbit Caliciviruses in Chile

Pathogenic lagoviruses (Rabbit hemorrhagic disease virus, RHDV) are widely spread across the world and are used in Australia and New Zealand to control populations of feral European rabbits. The spread of the non-pathogenic lagoviruses, e.g., rabbit calicivirus (RCV), is less well studied as the infection results in no clinical signs. Nonetheless, RCV has important implications for the spread of RHDV and rabbit biocontrol as it can provide varying levels of cross-protection against fatal infection with pathogenic lagoviruses. In Chile, where European rabbits are also an introduced species, myxoma virus was used for localised biocontrol of rabbits in the 1950s. To date, there have been no studies investigating the presence of lagoviruses in the Chilean feral rabbit population. In this study, liver and duodenum rabbit samples from central Chile were tested for the presence of lagoviruses and positive samples were subject to whole RNA sequencing and subsequent data analysis. Phylogenetic analysis revealed a novel RCV variant in duodenal samples that likely originated from European RCVs. Sequencing analysis also detected the presence of a rabbit astrovirus in one of the lagovirus-positive samples.

Varying conjunctival immune response adaptations of house finch populations to a rapidly evolving bacterial pathogen

Pathogen adaptations during host-pathogen co-evolution can cause the host balance between immunity and immunopathology to rapidly shift. However, little is known in natural disease systems about the immunological pathways optimised through the trade-off between immunity and self-damage. The evolutionary interaction between the conjunctival bacterial infection Mycoplasma gallisepticum (MG) and its avian host, the house finch (Haemorhous mexicanus), can provide insights into such adaptations in immune regulation. Here we use experimental infections to reveal immune variation in conjunctival tissue for house finches captured from four distinct populations differing in the length of their co-evolutionary histories with MG and their disease tolerance (defined as disease severity per pathogen load) in controlled infection studies. To differentiate contributions of host versus pathogen evolution, we compared house finch responses to one of two MG isolates: the original VA1994 isolate and a more evolutionarily derived one, VA2013. To identify differential gene expression involved in initiation of the immune response to MG, we performed 3'-end transcriptomic sequencing (QuantSeq) of samples from the infection site, conjunctiva, collected 3-days post-infection. In response to MG, we observed an increase in general pro-inflammatory signalling, as well as T-cell activation and IL17 pathway differentiation, associated with a decrease in the IL12/IL23 pathway signalling. The immune response was stronger in response to the evolutionarily derived MG isolate compared to the original one, consistent with known increases in MG virulence over time. The host populations differed namely in pre-activation immune gene expression, suggesting population-specific adaptations. Compared to other populations, finches from Virginia, which have the longest co-evolutionary history with MG, showed significantly higher expression of anti-inflammatory genes and Th1 mediators. This may explain the evolution of disease tolerance to MG infection in VA birds. We also show a potential modulating role of BCL10, a positive B- and T-cell regulator activating the NFKB signalling. Our results illuminate potential mechanisms of house finch adaptation to MG-induced immunopathology, contributing to understanding of the host evolutionary responses to pathogen-driven shifts in immunity-immunopathology trade-offs.

Genomic Tools in Biological Invasions: Current State and Future Frontiers

Human activities are accelerating rates of biological invasions and climate-driven range expansions globally, yet we understand little of how genomic processes facilitate the invasion process. Although most of the literature has focused on underlying phenotypic correlates of invasiveness, advances in genomic technologies are showing a strong link between genomic variation and invasion success. Here, we consider the ability of genomic tools and technologies to (i) inform mechanistic understanding of biological invasions and (ii) solve real-world issues in predicting and managing biological invasions. For both, we examine the current state of the field and discuss how genomics can be leveraged in the future. In addition, we make recommendations pertinent to broader research issues, such as data sovereignty, metadata standards, collaboration, and science communication best practices that will require concerted efforts from the global invasion genomics community.

Independently evolved pollution resistance in four killifish populations is largely explained by few variants of large effect

The genetic architecture of phenotypic traits can affect the mode and tempo of trait evolution. Human-altered environments can impose strong natural selection, where successful evolutionary adaptation requires swift and large phenotypic shifts. In these scenarios, theory predicts that adaptation is due to a few adaptive variants of large effect, but empirical studies that have revealed the genetic architecture of rapidly evolved phenotypes are rare, especially for populations inhabiting polluted environments. Fundulus killifish have repeatedly evolved adaptive resistance to extreme pollution in urban estuaries. Prior studies, including genome scans for signatures of natural selection, have revealed some of the genes and pathways important for evolved pollution resistance, and provide context for the genotype-phenotype association studies reported here. We created multiple quantitative trait locus (QTL) mapping families using progenitors from four different resistant populations, and using RAD-seq genetically mapped variation in sensitivity (developmental perturbations) following embryonic exposure to a model toxicant PCB-126. We found that one to two large-effect QTL loci accounted for resistance to PCB-mediated developmental toxicity. QTLs harbored candidate genes that govern the regulation of aryl hydrocarbon receptor (AHR) signaling. One QTL locus was shared across all populations and another was shared across three populations. One QTL locus showed strong signatures of recent natural selection in the corresponding wild population but another QTL locus did not. Some candidate genes for PCB resistance inferred from genome scans in wild populations were identified as QTL, but some key candidate genes were not. We conclude that rapidly evolved resistance to the developmental defects normally caused by PCB-126 is governed by few genes of large effect. However, other aspects of resistance beyond developmental phenotypes may be governed by additional loci, such that comprehensive resistance to PCB-126, and to the mixtures of chemicals that distinguish urban estuaries more broadly, may be more genetically complex.

MHCtools 1.5: Analysis of MHC Sequencing Data in R

The genes of the major histocompatibility complex (MHC) play a vital role in the vertebrate immune system and have attracted considerable interest in evolutionary biology. While the MHC has been characterized in detail in humans (human leukocyte antigen, HLA) and in model organisms such as the mouse, studies in non-model organisms often lack prior knowledge about structure, genetic variability, and evolutionary properties of this locus. MHC genotyping in non-model species commonly relies on PCR-based amplicon sequencing, and while several published protocols facilitate generation of MHC sequence data, there is a lack of transparent and standardized tools for downstream data analysis.Here, I present the R package MHCtools version 1.5, which contains 15 tools that (i) assist accurate MHC genotyping from high-throughput amplicon sequencing data, and provide standardized methods to analyze (ii) MHC diversity, (iii) MHC supertypes, and (iv) MHC haplotypes.I hope that MHCtools will be helpful in future studies of the MHC in non-model species and that it may help to advance our understanding of the important roles of the MHC in ecology and evolution.

Temporal collections to study invasion biology

Biological invasions represent an extraordinary opportunity to study evolution. This is because accidental or deliberate species introductions have taken place for centuries across large geographical scales, frequently prompting rapid evolutionary transitions in invasive populations. Until recently, however, the utility of invasions as evolutionary experiments has been hampered by limited information on the makeup of populations that were part of earlier invasion stages. Now, developments in ancient and historical DNA technologies, as well as the quickening pace of digitization for millions of specimens that are housed in herbaria and museums globally, promise to help overcome this obstacle. In this review, we first introduce the types of temporal data that can be used to study invasions, highlighting the timescale captured by each approach and their respective limitations. We then discuss how ancient and historical specimens as well as data available from prior invasion studies can be used to answer questions on mechanisms of (mal)adaptation, rates of evolution, or community-level changes during invasions. By bridging the gap between contemporary and historical invasive populations, temporal data can help us connect pattern to process in invasion science. These data will become increasingly important if invasions are to achieve their full potential as experiments of evolution in nature.

Effect of epidemic diseases on wild animal conservation

Under the background of global species extinction, the impact of epidemic diseases on wild animal protection is increasingly prominent. Here, we review and synthesize the literature on this topic, and discuss the relationship between diseases and biodiversity. Diseases usually reduce species diversity by decreasing or extinction of species populations, but also accelerate species evolution and promote species diversity. At the same time, species diversity can regulate disease outbreaks through dilution or amplification effects. The synergistic effect of human activities and global change is emphasized, which further aggravates the complex relationship between biodiversity and diseases. Finally, we emphasize the importance of active surveillance of wild animal diseases, which can protect wild animals from potential diseases, maintain population size and genetic variation, and reduce the damage of diseases to the balance of the whole ecosystem and human health. Therefore, we suggest that a background survey of wild animal populations and their pathogens should be carried out to assess the impact of potential outbreaks on the population or species level. The mechanism of dilution and amplification effect between species diversity and diseases of wild animals should be further studied to provide a theoretical basis and technical support for human intervention measures to change biodiversity. Most importantly, we should closely combine the protection of wild animals with the establishment of an active surveillance, prevention, and control system for wild animal epidemics, in an effort to achieve a win-win situation between wild animal protection and disease control.

Revisiting the Role of Genetic Variation in Adaptation

AbstractEvolutionary biologists have thought about the role of genetic variation during adaptation for a very long time-before we understood the organization of the genetic code, the provenance of genetic variation, and how such variation influenced the phenotypes on which natural selection acts. Half a century after the discovery of the structure of DNA and the unraveling of the genetic code, we have a rich understanding of these problems and the means to both delve deeper and widen our perspective across organisms and natural populations. The 2022 Vice Presidential Symposium of the American Society of Naturalists highlighted examples of recent insights into the role of genetic variation in adaptive processes, which are compiled in this special section. The work was conducted in different parts of the world, included theoretical and empirical studies with diverse organisms, and addressed distinct aspects of how genetic variation influences adaptation. In our introductory article to the special section, we discuss some important recent insights about the generation and maintenance of genetic variation, its impacts on phenotype and fitness, its fate in natural populations, and its role in driving adaptation. By placing the special section articles in the broader context of recent developments, we hope that this overview will also serve as a useful introduction to the field.

Selection on standing genetic variation mediates convergent evolution in extremophile fish

Hydrogen sulfide is a toxic gas that disrupts numerous biological processes, including energy production in the mitochondria, yet fish in the Poecilia mexicana species complex have independently evolved sulfide tolerance several times. Despite clear evidence for convergence at the phenotypic level in these fishes, it is unclear if the repeated evolution of hydrogen sulfide tolerance is the result of similar genomic changes. To address this gap, we used a targeted capture approach to sequence genes associated with sulfide processes and toxicity from five sulfidic and five nonsulfidic populations in the species complex. By comparing sequence variation in candidate genes to a reference set, we identified similar population structure and differentiation, suggesting that patterns of variation in most genes associated with sulfide processes and toxicity are due to demographic history and not selection. But the presence of tree discordance for a subset of genes suggests that several loci are evolving divergently between ecotypes. We identified two differentiation outlier genes that are associated with sulfide detoxification in the mitochondria that have signatures of selection in all five sulfidic populations. Further investigation into these regions identified long, shared haplotypes among sulfidic populations. Together, these results reveal that selection on standing genetic variation in putatively adaptive genes may be driving phenotypic convergence in this species complex.

Detecting parallel polygenic adaptation to novel evolutionary pressure in wild populations: a case study in Atlantic cod (Gadus morhua)

Populations can adapt to novel selection pressures through dramatic frequency changes in a few genes of large effect or subtle shifts in many genes of small effect. The latter (polygenic adaptation) is expected to be the primary mode of evolution for many life-history traits but tends to be more difficult to detect than changes in genes of large effect. Atlantic cod (Gadus morhua) were subjected to intense fishing pressure over the twentieth century, leading to abundance crashes and a phenotypic shift toward earlier maturation across many populations. Here, we use spatially replicated temporal genomic data to test for a shared polygenic adaptive response to fishing using methods previously applied to evolve-and-resequence experiments. Cod populations on either side of the Atlantic show covariance in allele frequency change across the genome that are characteristic of recent polygenic adaptation. Using simulations, we demonstrate that the degree of covariance in allele frequency change observed in cod is unlikely to be explained by neutral processes or background selection. As human pressures on wild populations continue to increase, understanding and attributing modes of adaptation using methods similar to those demonstrated here will be important in identifying the capacity for adaptive responses and evolutionary rescue. This article is part of the theme issue 'Detecting and attributing the causes of biodiversity change: needs, gaps and solutions'.

Rapid adaptation to a novel pathogen through disease tolerance in a wild songbird

Animal hosts can adapt to emerging infectious disease through both disease resistance, which decreases pathogen numbers, and disease tolerance, which limits damage during infection without limiting pathogen replication. Both resistance and tolerance mechanisms can drive pathogen transmission dynamics. However, it is not well understood how quickly host tolerance evolves in response to novel pathogens or what physiological mechanisms underlie this defense. Using natural populations of house finches (Haemorhous mexicanus) across the temporal invasion gradient of a recently emerged bacterial pathogen (Mycoplasma gallisepticum), we find rapid evolution of tolerance (<25 years). In particular, populations with a longer history of MG endemism have less pathology but similar pathogen loads compared with populations with a shorter history of MG endemism. Further, gene expression data reveal that more-targeted immune responses early in infection are associated with tolerance. These results suggest an important role for tolerance in host adaptation to emerging infectious diseases, a phenomenon with broad implications for pathogen spread and evolution.

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.

How predictable is adaptation from standing genetic variation? Experimental evolution in Drosophila highlights the central role of redundancy and linkage disequilibrium

Experimental evolution is well-suited to test the predictability of evolution without the confounding effects of inaccurate forecasts about future environments. Most of the literature about parallel (and thus predictable) evolution has been carried out in asexual microorganisms, which adapt by de novo mutations. Nevertheless, parallel evolution has also been studied in sexual species at the genomic level. Here, I review the evidence for parallel evolution in Drosophila, the best-studied obligatory outcrossing model for adaptation from standing genetic variation in the laboratory. Similar to asexual microorganisms, evidence for parallel evolution varies between the focal hierarchical levels. Selected phenotypes consistently respond in a very predicable way, but the underlying allele frequency changes are much less predictable. The most important insight is that the predictability of the genomic selection response for polygenic traits depends highly on the founder population and to a much lesser extent on the selection regime. This implies that predicting adaptive genomic response is challenging and requires a good understanding of the adaptive architecture (including linkage disequilibrium) in the ancestral populations. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Polygenic routes lead to parallel altitudinal adaptation in Heliosperma pusillum (Caryophyllaceae)

Understanding how organisms adapt to the environment is a major goal of modern biology. Parallel evolution-the independent evolution of similar phenotypes in different populations-provides a powerful framework to investigate the evolutionary potential of populations, the constraints of evolution, its repeatability and therefore its predictability. Here, we quantified the degree of gene expression and functional parallelism across replicated ecotype formation in Heliosperma pusillum (Caryophyllaceae), and gained insights into the architecture of adaptive traits. Population structure analyses and demographic modelling support a previously formulated hypothesis of parallel polytopic divergence of montane and alpine ecotypes. We detect a large proportion of differentially expressed genes (DEGs) underlying divergence within each replicate ecotype pair, with a strikingly low number of shared DEGs across pairs. Functional enrichment of DEGs reveals that the traits affected by significant expression divergence are largely consistent across ecotype pairs, in strong contrast to the nonshared genetic basis. The remarkable redundancy of differential gene expression indicates a polygenic architecture for the diverged adaptive traits. We conclude that polygenic traits appear key to opening multiple routes for adaptation, widening the adaptive potential of organisms.

Increasing ecological heterogeneity can constrain biopesticide resistance evolution

Microbial biopesticides containing living parasites are valuable emerging crop protection technologies against insect pests, but they are vulnerable to resistance evolution. Fortunately, the fitness of alleles that provide resistance, including to parasites used in biopesticides, frequently depends on parasite identity and environmental conditions. This context-specificity suggests a sustainable approach to biopesticide resistance management through landscape diversification. To mitigate resistance risks, we advocate increasing the range of biopesticides available to farmers, whilst simultaneously encouraging other aspects of landscape-wide crop heterogeneity that can generate variable selection on resistance alleles. This approach requires agricultural stakeholders to prioritize diversity as well as efficiency, both within agricultural landscapes and the biocontrol marketplace.

Pathological changes and viral antigen distribution in tissues of Iberian hare (Lepus granatensis) naturally infected with the emerging recombinant myxoma virus (ha-MYXV)

BACKGROUND

A cross-species jump was confirmed in 2018, when a novel recombinant myxoma virus (MYXV) (ha-MYXV) caused high mortality in Iberian hare (Lepus granatensis) in the Iberian Peninsula.

METHOD

The aim of this study was to evaluate the main lesions, tissular distribution and target cells of ha-MYXV in Iberian hare. Gross postmortem examinations and histological and immunohistochemical studies to detect ha-MYXV were carried out in 28 animals that were confirmed as ha-MYXV positive by PCR.

RESULTS

The main macroscopic lesions were bilateral blepharoconjunctivitis, epistaxis, intense congestion and oedema in several organs and some internal haemorrhages. Visible myxomas were not found. Histopathological examination revealed hyperplastic epidermis with predominant hyperkeratosis and myxoid matrix in the dermis. ha-MYXV-positive keratinocytes showed hydropic degeneration and cytoplasmic inclusion bodies. Alveolar oedema, interstitial pneumonia, dramatic lymphoid depletion in the spleen and necrosis in the liver and testis were observed. ha-MYXV was mainly detected in epithelial and myxoma cells in the skin, and also in macrophages, lymphocytes, fibroblasts and endothelial cells in several organs, as well as in hepatocytes and Leydig cells.

LIMITATIONS

A non-homogeneous number of samples were included in all the animals. Future experimental studies with controlled variables are necessary.

CONCLUSION

These findings correspond to an unusual form of myxomatosis, characterised by an acute or hyperacute presentation.

Role of cytokines in poxvirus host tropism and adaptation

Poxviruses are a diverse family of double-stranded DNA viruses that cause mild-to-severe disease in selective hosts, including humans. Although most poxviruses are restricted to their hosts, some members can leap host species and cause zoonotic diseases and, therefore, are genuine threats to human and animal health. The recent global spread of monkeypox in humans suggests that zoonotic poxviruses can adapt to a new host, spread rapidly in the new host, and evolve to better evade host innate barriers. Unlike many other viruses, poxviruses express an extensive repertoire of self-defense proteins that play a vital role in the evasion of host innate and adaptive immune responses in their newest host species. The function of these viral immune modulators and host-specific cytokine responses can result in different host tropism and poxvirus disease progression. Here, we review the role of different cytokines that control poxvirus host tropism and adaptation.

Retrospective serological and molecular survey of myxoma or antigenically related virus in the Iberian hare, Lepus granatensis

The 2018 outbreak of myxomatosis in the Iberian hare (Lepus granatensis) has been hypothesized to originate from a species jump of the rabbit-associated myxoma virus (MYXV), after natural recombination with an unknown poxvirus. Iberian hares were long considered resistant to myxomatosis as no prior outbreaks were reported. To provide insights into the emergence of this recombinant virus (ha-MYXV), we investigated serum samples from 451 Iberian hares collected over two time periods almost two decades apart, 1994-1999 and 2017-2019 for the presence of antibodies and MYXV-DNA. First, we screened all serum samples using a rabbit commercial indirect ELISA (iELISA) and then tested a subset of these samples in parallel using indirect immunofluorescence test (IFT), competitive ELISA (cELISA) and qPCR targeting M000.5L/R gene conserved in MYXV and ha-MYXV. The cut-off of iELISA relative index 10 = 6.1 was selected from a semiparametric finite mixture analysis aiming to minimize the probability of false positive results. Overall, MYXV related-antibodies were detected in 57 hares (12.6%) including 38 apparently healthy hares (n = 10, sampled in 1994-1999, none MYXV-DNA positive, and n = 28 sampled in 2017-2019 of which four were also ha-MYXV-DNA positive) and 19 found-dead and ha-MYXV-DNA-positive sampled in 2018-2019. Interestingly, four seronegative hares sampled in 1997 were MYXV-DNA positive by qPCR, the result being confirmed by sequencing of three of them. For the Iberian hares hunted or live trapped (both apparently health), seroprevalence was significantly higher in 2017-2019 (13.0%, CI95% 9.2-18.2%) than in 1994-1999 (5.4%, CI95% 3.0-9.6%) (p = .009). Within the second period, seroprevalence was significantly higher in 2019 compared to 2017 (24.7 vs 1.7% considering all the sample, p = .007), and lower during the winter than the autumn (p < .001). While our molecular and serological results show that Iberian hares have been in contact with MYXV or an antigenically similar virus at least since 1996, they also show an increase in seroprevalence in 2018-2019. The remote contact with MYXV may have occurred with strains that circulated in rabbits, or with unnoticed strains already circulating in Iberian hare populations. This work strongly suggests the infection of Iberian hares with MYXV or an antigenically related virus, at least 20 years before the severe virus outbreaks were registered in 2018.

Spatiotemporal monitoring of myxomatosis in European wild rabbit (Oryctolagus cuniculus) in Spanish Mediterranean ecosystems

A long-term active epidemiological surveillance programme was conducted to determine seroprevalence to myxoma virus (MYXV), infection prevalence and spatiotemporal patterns and factors associated with MYXV circulation in wild rabbits (Oryctolagus cuniculus) in Spanish Mediterranean ecosystems. A total of 2376 animals were sampled over four study periods: 2009-2012 (P1), 2012-2015 (P2), 2015-2018 (P3) and 2018-2021 (P4). Antibodies against MYXV were detected by a commercial indirect ELISA in 59.9% (1424/2376; 95% CI: 58.0-61.9) of wild rabbits. At least one seropositive animal was detected on 131 (96.3%) of 136 game estates sampled. MYXV infection was confirmed by PCR in 94 of 1063 (8.8%; 95% CI: 7.3-10.7) wild rabbits. Circulation of the novel recombinant MYXV (ha-MYXV) was not found in wild rabbits analysed during P4. Five statistically significant spatiotemporal clusters of high MYXV seroprevalence were identified using a Bernoulli model: one in P2 and four in P3. A generalized linear mixed model (GLMM) analysis identified sampling season (autumn), age (adult and juvenile), outbreaks of myxomatosis in the month prior to sampling, mean annual temperature, humidity and seropositivity to rabbit haemorrhagic disease virus as factors potentially linked with MYXV seropositivity. GLMM analysis identified outbreaks of myxomatosis in the month prior to sampling, MYXV seropositivity and presence of lesions compatible with myxomatosis as factors associated with MYXV infection. The results indicate high exposure, widespread but non-homogeneous distribution, and endemic circulation of MYXV in wild rabbit populations in southern Spain during the last decade. Prevalence of antibodies against MYXV showed fluctuations both within the year and over the study periods, revealing variations in the immunity of wild rabbit populations in Mediterranean ecosystems that could increase the risk of MYXV re-emergence in immunologically naïve populations. The present study highlights the importance of long-term surveillance to better understand the epidemiology of MYXV in wild lagomorphs.

EGFL7 drives the evolution of resistance to EGFR inhibitors in lung cancer by activating NOTCH signaling

Accumulating evidence supports evolutionary trait of drug resistance. Like resilience in other systems, most tumor cells experience drug-tolerant state before full resistance acquired. However, the underlying mechanism is still poorly understood. Here, we identify that EGF like domain multiple 7 (EGFL7) is a responsive gene to epidermal growth factor receptor (EGFR) kinase inhibition during a period when tumors are decimated. Moreover, our data reveal that the adaptive increase of EGFL7 during this process is controlled by the depression of nonsense-mediated mRNA decay (NMD) pathway. Upregulation of EGFL7 activates NOTCH signaling in lung cancer cells, which slows down the decrease of c-Myc caused by EGFR inhibition, thereby helping the survival of cancer cells. Our data, taken together, demonstrate that EGFL7 is a driver gene for resistance to EGFR kinase inhibition, and suggest that targeting EGFL7/NOTCH signaling may improve the clinical benefits of EGFR inhibitors in patients with EGFR mutant tumors.

Divergent Evolutionary Pathways of Myxoma Virus in Australia: Virulence Phenotypes in Susceptible and Partially Resistant Rabbits Indicate Possible Selection for Transmissibility

To characterize the ongoing evolution of myxoma virus in Australian rabbits, we used experimental infections of laboratory rabbits to determine the virulence and disease phenotypes of recent virus isolates. The viruses, collected between 2012 and 2015, fell into three lineages, one of which, lineage c, experienced a punctuated increase in evolutionary rate. All viruses were capable of causing acute death with aspects of neutropenic septicemia, characterized by minimal signs of myxomatosis, the occurrence of pulmonary edema and bacteria invasions throughout internal organs, but with no inflammatory response. For the viruses of highest virulence all rabbits usually died at this point. In more attenuated viruses, some rabbits died acutely, while others developed an amyxomatous phenotype. Rabbits that survived for longer periods developed greatly swollen cutaneous tissues with very high virus titers. This was particularly true of lineage c viruses. Unexpectedly, we identified a line of laboratory rabbits with some innate resistance to myxomatosis and used these in direct comparisons with the fully susceptible rabbit line. Importantly, the same disease phenotype occurred in both susceptible and resistant rabbits, although virulence was shifted toward more attenuated grades in resistant animals. We propose that selection against inflammation at cutaneous sites prolongs virus replication and enhances transmission, leading to the amyxomatous phenotype. In some virus backgrounds this creates an immunosuppressive state that predisposes to high virulence and acute death. The alterations in disease pathogenesis, particularly the overwhelming bacterial invasions that characterize the modern viruses, suggest that their virulence grades are not directly comparable with earlier studies.

IMPORTANCE The evolution of the myxoma virus (MYXV) following its release as a biological control for European rabbits in Australia is the textbook example of the coevolution of virus virulence and host resistance. However, most of our knowledge of MYXV evolution only covers the first few decades of its spread in Australia and often with little direct connection between how changes in virus phenotype relate to those in the underlying virus genotype. By conducting detailed experimental infections of recent isolates of MYXV in different lines of laboratory rabbits, we examined the ongoing evolution of MYXV disease phenotypes. Our results reveal a wide range of phenotypes, including an amyxomatous type, as well as the impact of invasive bacteria, that in part depended on the level of rabbit host resistance. These results provide a unique insight into the complex virus and host factors that combine to shape disease phenotype and viral evolution.

Myxoma virus lacking the host range determinant M062 stimulates cGAS-dependent type 1 interferon response and unique transcriptomic changes in human monocytes/macrophages

The evolutionarily successful poxviruses possess effective and diverse strategies to circumvent or overcome host defense mechanisms. Poxviruses encode many immunoregulatory proteins to evade host immunity to establish a productive infection and have unique means of inhibiting DNA sensing-dependent type 1 interferon (IFN-I) responses, a necessity given their dsDNA genome and exclusively cytoplasmic life cycle. We found that the key DNA sensing inhibition by poxvirus infection was dominant during the early stage of poxvirus infection before DNA replication. In an effort to identify the poxvirus gene products which subdue the antiviral proinflammatory responses (e.g., IFN-I response), we investigated the function of one early gene that is the known host range determinant from the highly conserved poxvirus host range C7L superfamily, myxoma virus (MYXV) M062. Host range factors are unique features of poxviruses that determine the species and cell type tropism. Almost all sequenced mammalian poxviruses retain at least one homologue of the poxvirus host range C7L superfamily. In MYXV, a rabbit-specific poxvirus, the dominant and broad-spectrum host range determinant of the C7L superfamily is the M062R gene. The M062R gene product is essential for MYXV infection in almost all cells tested from different mammalian species and specifically inhibits the function of host Sterile αMotif Domain-containing 9 (SAMD9), as M062R-null (ΔM062R) MYXV causes abortive infection in a SAMD9-dependent manner. In this study we investigated the immunostimulatory property of the ΔM062R. We found that the replication-defective ΔM062R activated host DNA sensing pathway during infection in a cGAS-dependent fashion and that knocking down SAMD9 expression attenuated proinflammatory responses. Moreover, transcriptomic analyses showed a unique feature of the host gene expression landscape that is different from the dsDNA alone-stimulated inflammatory state. This study establishes a link between the anti-neoplastic function of SAMD9 and the regulation of innate immune responses.

Poxviruses encode a group of genes called host range determinants to maintain or expand their host tropism. The mechanism by which many viral host range factors function remains elusive. Some host range factors possess immunoregulatory functions responsible for evading or subduing host immune defense mechanisms. Most known immunoregulatory proteins encoded by poxviruses are dispensable for viral replication in vitro. The uniqueness of MYXV M062R is that it is essential for viral infection in vitro and belongs to one of the most conserved poxvirus host range families, the C7L superfamily. There is one known host target of the MYXV M062 protein, SAMD9. SAMD9 is constitutively expressed in mammalian cells and exclusively present in the cytoplasm with an anti-neoplastic function. Humans with deleterious mutations in SAMD9 present disease that ranges from lethality at a young age to a predisposition to myelodysplastic syndromes (MDS) that often require bone marrow transplantation. More importantly, SAMD9 serves as an important antiviral intrinsic molecule to many viruses. The cellular function of SAMD9 remains unclear mostly due to the difficulty of studying this protein, i.e., its large size, long half-life, and its constitutive expression in most cells. In this study we used M062R-null MYXV as a tool to study SAMD9 function and report a functional link between SAMD9 and the regulation of the proinflammatory responses triggered by cGAS-dependent DNA sensing.

A single introduction of wild rabbits triggered the biological invasion of Australia

Biological invasions are a major cause of environmental and economic disruption. While ecological factors are key determinants of their success, the role of genetics has been more challenging to demonstrate. The colonization of Australia by the European rabbit is one of the most iconic and devastating biological invasions in recorded history. Here, we show that despite numerous introductions over a 70-y period, this invasion was triggered by a single release of a few animals that spread thousands of kilometers across the continent. We found genetic support for historical accounts that these were English rabbits imported in 1859 by a settler named Thomas Austin and traced the origin of the invasive population back to his birthplace in England. We also find evidence of additional introductions that established local populations but have not spread geographically. Combining genomic and historical data we show that, contrary to the earlier introductions, which consisted mostly of domestic animals, the invasive rabbits had wild ancestry. In New Zealand and Tasmania, rabbits also became a pest several decades after being introduced. We argue that the common denominator of these invasions was the arrival of a new genotype that was better adapted to the natural environment. These findings demonstrate how the genetic composition of invasive individuals can determine the success of an introduction and provide a mechanism by which multiple introductions can be required for a biological invasion.

Standing Genetic Diversity and Transmission Bottleneck Size Drive Adaptation in Bacteriophage Qβ

A critical issue to understanding how populations adapt to new selective pressures is the relative contribution of the initial standing genetic diversity versus that generated de novo. RNA viruses are an excellent model to study this question, as they form highly heterogeneous populations whose genetic diversity can be modulated by factors such as the number of generations, the size of population bottlenecks, or exposure to new environment conditions. In this work, we propagated at nonoptimal temperature (43 °C) two bacteriophage Qβ populations differing in their degree of heterogeneity. Deep sequencing analysis showed that, prior to the temperature change, the most heterogeneous population contained some low-frequency mutations that had previously been detected in the consensus sequences of other Qβ populations adapted to 43 °C. Evolved populations with origin in this ancestor reached similar growth rates, but the adaptive pathways depended on the frequency of these standing mutations and the transmission bottleneck size. In contrast, the growth rate achieved by populations with origin in the less heterogeneous ancestor did depend on the transmission bottleneck size. The conclusion is that viral diversification in a particular environment may lead to the emergence of mutants capable of accelerating adaptation when the environment changes.

Sequence Capture From Historical Museum Specimens: Maximizing Value for Population and Phylogenomic Studies

The application of high-throughput, short-read sequencing to degraded DNA has greatly increased the feasibility of generating genomic data from historical museum specimens. While many published studies report successful sequencing results from historical specimens; in reality, success and quality of sequence data can be highly variable. To examine predictors of sequencing quality, and methodological approaches to improving data accuracy, we generated and analyzed genomic sequence data from 115 historically collected museum specimens up to 180 years old. Data span both population genomic and phylogenomic scales, including historically collected specimens from 34 specimens of four species of Australian rock-wallabies (genus Petrogale) and 92 samples from 79 specimens of Australo-Papuan murine rodents (subfamily Murinae). For historical rodent specimens, where the focus was sampling for phylogenomics, we found that regardless of specimen age, DNA sequence libraries prepared from toe pad or bone subsamples performed significantly better than those taken from the skin (in terms of proportion of reads on target, number of loci captured, and data accuracy). In total, 93% of DNA libraries from toe pad or bone subsamples resulted in reliable data for phylogenetic inference, compared to 63% of skin subsamples. For skin subsamples, proportion of reads on target weakly correlated with collection year. Then using population genomic data from rock-wallaby skins as a test case, we found substantial improvement in final data quality by mapping to a high-quality “closest sister” de novo assembly from fresh tissues, compared to mapping to a sample-specific historical de novo assembly. Choice of mapping approach also affected final estimates of the number of segregating sites and Watterson's θ, both important parameters for population genomic inference. The incorporation of accurate and reliable sequence data from historical specimens has important outcomes for evolutionary studies at both population and phylogenomic scales. By assessing the outcomes of different approaches to specimen subsampling, library preparation and bioinformatic processing, our results provide a framework for increasing sequencing success for irreplaceable historical specimens.

MHCtools - an R package for MHC high-throughput sequencing data: genotyping, haplotype and supertype inference, and downstream genetic analyses in non-model organisms

The major histocompatibility complex (MHC) plays a central role in the vertebrate adaptive immune system and has been of long-term interest in evolutionary biology. While several protocols have been developed for MHC genotyping, there is a lack of transparent and standardized tools for downstream analysis of MHC data. Here, we present the R package MHCtools and demonstrate the use of its functions to (i) assist accurate MHC genotyping from high-throughput amplicon-sequencing data, (ii) infer functional MHC supertypes using bootstrapped clustering analysis, (iii) identify segregating MHC haplotypes from family data, and (iv) analyse functional and genetic distances between MHC sequences. We employed MHCtools to analyse MHC class I (MHC-I) amplicon data of 559 great reed warblers (Acrocephalus arundinaceus). We identified 390 MHC-I alleles which clustered into 14 functional supertypes. A phylogenetic analysis and analyses of positive selection suggested that the MHC-I alleles belonged to several distinct functional groups. We furthermore identified 107 segregating haplotypes among 116 families, and found substantial variation in diversity with 4-21 MHC-I alleles and 3-13 MHC-I supertypes per haplotype. Finally, we show that the great reed warbler haplotypes harboured combinations of MHC-I supertypes with greater functional divergence than observed in simulated populations of possible haplotypes, a result that is in accordance with the divergent allele advantage hypothesis. Our study demonstrates the power of MHCtools to support genotyping and analysis of MHC in non-model species, which we hope will encourage broad implementation among researchers in MHC genetics and evolution.

Long term patterns of association between MHC and helminth burdens in the bank vole support Red Queen dynamics

Major histocompatibility complex (MHC) genes encode proteins crucial for adaptive immunity of vertebrates. Negative frequency-dependent selection (NFDS), resulting from adaptation of parasites to common MHC types, has been hypothesized to maintain high, functionally relevant polymorphism of MHC, but demonstration of this relationship has remained elusive. In particular, differentiation of NFDS from fluctuating selection, resulting from changes in parasite communities in time and space (FS), has proved difficult in short-term studies. Here, we used temporal data, accumulated through long-term monitoring of helminths infecting bank voles (Myodes glareolus), to test specific predictions of NFDS on MHC class II. Data were collected in three, moderately genetically differentiated subpopulations in Poland, which were characterized by some stable spatiotemporal helminth communities but also events indicating introduction of new species and loss of others. We found a complex association between individual MHC diversity and species richness, where intermediate numbers of DRB supertypes correlated with lowest species richness, but the opposite was true for DQB supertypes - arguing against universal selection for immunogenetic optimality. We also showed that particular MHC supertypes explain a portion of the variance in prevalence and abundance of helminths, but this effect was subpopulation-specific, which is consistent with both NFDS and FS. Finally, in line with NFDS, we found that certain helminths that have recently colonized or spread in a given subpopulation, more frequently or intensely infected voles with MHC supertypes that have been common in the recent past. Overall, our results highlight complex spatial and temporal patterns of MHC-parasite associations, the latter being consistent with Red Queen coevolutionary dynamics.

The genetics of adaptation in freshwater Eurasian shad (Alosa)

Studying the genetics of phenotypic convergence can yield important insights into adaptive evolution. Here, we conducted a comparative genomic study of four lineages (species and subspecies) of anadromous shad (Alosa) that have independently evolved life cycles entirely completed in freshwater. Three naturally diverged (A. fallax lacustris, A. f. killarnensis, and A. macedonica), and the fourth (A. alosa) was artificially landlocked during the last century. To conduct this analysis, we assembled and annotated a draft of the A. alosa genome and generated whole‐genome sequencing for 16 anadromous and freshwater populations of shad. Widespread evidence for parallel genetic changes in freshwater populations within lineages was found. In freshwater A. alosa, which have only been diverging for tens of generations, this shows that parallel adaptive evolution can rapidly occur. However, parallel genetic changes across lineages were comparatively rare. The degree of genetic parallelism was not strongly related to the number of shared polymorphisms between lineages, thus suggesting that other factors such as divergence among ancestral populations or environmental variation may influence genetic parallelism across these lineages. These overall patterns were exemplified by genetic differentiation involving a paralog of ATPase‐α1 that appears to be under selection in just two of the more distantly related lineages studied, A. f. lacustris and A. alosa. Our findings provide insights into the genetic architecture of adaptation and parallel evolution along a continuum of population divergence.

Identification of a Novel Myxoma Virus C7-Like Host Range Factor That Enabled a Species Leap from Rabbits to Hares

Myxoma virus (MYXV) is naturally found in rabbit Sylvilagus species and is known to cause lethal myxomatosis in European rabbits (Oryctolagus cuniculus). In 2019, an MYXV strain (MYXV strain Toledo [MYXV-Tol]) causing myxomatosis-like disease in Iberian hares (Lepus granatensis) was identified. MYXV-Tol acquired a recombinant region of ∼2.8 kb harboring several new genes, including a novel host range gene (M159) that we show to be an orthologous member of the vaccinia virus C7 host range family. Here, to test whether M159 alone has enabled MYXV to alter its host range to Iberian hares, several recombinant viruses were generated, including an MYXV-Tol ΔM159 (knockout) strain. While MYXV-Tol underwent fully productive infection in hare HN-R cells, neither the wild-type MYXV-Lau strain (lacking M159) nor vMyxTol-ΔM159 (deleted for M159) was able to infect and replicate, showing that the ability of MYXV-Tol to infect these cells and replicate depends on the presence of M159. Similar to other C7L family members, M159 was shown to be expressed as an early/late gene but was translocated into the nucleus at later time points, indicating that further studies are needed to elucidate its role in the nucleus. Finally, in rabbit cells, the M159 protein did not contribute to increased replication but was able to upregulate the replication levels of MYXV in nonpermissive and semipermissive human cancer cells, suggesting that the M159-targeted pathway is conserved across mammalian species. Altogether, these observations demonstrate that the M159 protein plays a critical role in determining the host specificity of MYXV-Tol in hare and human cells by imparting new host range functions. IMPORTANCE The coevolution of European rabbit populations and MYXV is a textbook example of an arms race between a pathogen and a host. Recently, a recombinant MYXV (MYXV-Tol) crossed the species barrier by jumping from leporid species to another species, causing lethal myxomatosis-like disease. Given the highly pathogenic nature of this new virus in hares and the incidences of other poxvirus cross-species spillovers into other animals, including humans, it is important to understand how and why MYXV-Tol was able to become virulent in a new host species. The results presented clearly demonstrate that M159 is the key factor allowing MYXV-Tol replication in hare cells by imparting new host range functions. These results have the potential to improve current knowledge about the virulence of poxviruses and provide a platform to better understand the new MYXV-Tol, rendering the virus capable of leaping into a new host species.

Genomic associations with poxvirus across divergent island populations in Berthelot's pipit

Understanding the mechanisms and genes that enable animal populations to adapt to pathogens is important from an evolutionary, health and conservation perspective. Berthelot's pipit (Anthus berthelotii) experiences extensive and consistent spatial heterogeneity in avian pox infection pressure across its range of island populations, thus providing an excellent system with which to examine how pathogen-mediated selection drives spatial variation in immunogenetic diversity. Here we test for evidence of genetic variation associated with avian pox at both an individual and population-level. At the individual level, we find no evidence that variation in MHC class I and TLR4 (both known to be important in recognising viral infection) was associated with pox infection within two separate populations. However, using genotype-environment association (Bayenv) in conjunction with genome-wide (ddRAD-seq) data, we detected strong associations between population-level avian pox prevalence and allele frequencies of single nucleotide polymorphisms (SNPs) at a number of sites across the genome. These sites were located within genes involved in cellular stress signalling and immune responses, many of which have previously been associated with responses to viral infection in humans and other animals. Consequently, our analyses indicates that pathogen-mediated selection may play a role in shaping genomic variation among relatively recently colonised island bird populations and highlights the utility of genotype-environment associations for identifying candidate genes potentially involved in host-pathogen interactions.

Cross-protection, infection, and case fatality rates in wild European rabbits experimentally challenged with different rabbit haemorrhagic disease viruses

Rabbit haemorrhagic disease virus 2 (RHDV2) is now the dominant calicivirus circulating in wild rabbit populations in Australia. This study compared the infection and case fatality rates of RHDV2 and two RHDVs in wild rabbits, as well as their ability to overcome immunity to the respective other strains. Wild rabbits were allocated to groups either blindly or based on prescreening for RHDV/RHDV2 antibodies at capture. Rabbits were monitored regularly until their death or humane killing at 7 days post infection. Liver and eyeball samples were collected for lagovirus testing and aging rabbits, respectively. At capture, rabbits showed high seroprevalence to RHDV2 but not to RHDV. In RHDV/RHDV2 seronegative rabbits at capture, infection rates were highest in those inoculated with RHDV2 (81.8%, 18/22), followed by K5 (53.8%, 7/13) and CZECH (40.0%, 2/5), but these differences were not statistically significant. In rabbits with previous exposure to RHDV2 at capture, infection rates were highest when inoculated with K5 (59.6%, 31/52) followed by CZECH (46.0%, 23/50), with infection rates higher in younger rabbits for both viruses. In RHDV/RHDV2 seronegative rabbits at capture, case fatality rates were highest for those inoculated with K5 (71.4%), followed by RHDV2 (50.0%) and CZECH (50.0%). In rabbits with previous exposure to RHDV2 at capture, case fatality rates were highest in rabbits inoculated with K5 (12.9%) followed by CZECH (8.7%), with no case fatalities following RHDV2 inoculation. Case fatality rates did not differ significantly between inoculums in either serostatus group at capture. Based on multivariable modelling, time to death post RHDV inoculation increased in rabbits with recent RHDV2 exposure compared to seronegative rabbits and with age. The results suggest that RHDV2 may cause higher mortalities than other variants in seronegative rabbit populations but that K5 may be more effective in reducing rabbit populations in an RHDV2-dominant landscape. This article is protected by copyright. All rights reserved.

On the genetic architecture of rapidly adapting and convergent life history traits in guppies

The genetic basis of traits shapes and constrains how adaptation proceeds in nature; rapid adaptation can proceed using stores of polygenic standing genetic variation or hard selective sweeps, and increasing polygenicity fuels genetic redundancy, reducing gene re-use (genetic convergence). Guppy life history traits evolve rapidly and convergently among natural high- and low-predation environments in northern Trinidad. This system has been studied extensively at the phenotypic level, but little is known about the underlying genetic architecture. Here, we use four independent F2 QTL crosses to examine the genetic basis of seven (five female, two male) guppy life history phenotypes and discuss how these genetic architectures may facilitate or constrain rapid adaptation and convergence. We use RAD-sequencing data (16,539 SNPs) from 370 male and 267 female F2 individuals. We perform linkage mapping, estimates of genome-wide and per-chromosome heritability (multi-locus associations), and QTL mapping (single-locus associations). Our results are consistent with architectures of many loci of small-effect for male age and size at maturity and female interbrood period. Male trait associations are clustered on specific chromosomes, but female interbrood period exhibits a weak genome-wide signal suggesting a potentially highly polygenic component. Offspring weight and female size at maturity are also associated with a single significant QTL each. These results suggest rapid, repeatable phenotypic evolution of guppies may be facilitated by polygenic trait architectures, but subsequent genetic redundancy may limit gene re-use across populations, in agreement with an absence of strong signatures of genetic convergence from recent analyses of wild guppies.

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.

Comparing genome scans among species of the stickleback order reveals three different patterns of genetic diversity

Comparing genome scans among species is a powerful approach for investigating the patterns left by evolutionary processes. In particular, this offers a way to detect candidate genes that drive convergent evolution. We compared genome scan results to investigate if patterns of genetic diversity and divergence are shared among divergent species within the stickleback order (Gasterosteiformes): the threespine stickleback (Gasterosteus aculeatus), ninespine stickleback (Pungitius pungitus), and tubesnout (Aulorhynchus flavidus). Populations were sampled from the southern and northern edges of each species' range, to identify patterns associated with latitudinal changes in genetic diversity. Weak correlations in genetic diversity (F ST and expected heterozygosity) and three different patterns in the genomic landscape were found among these species. Additionally, no candidate genes for convergent evolution were detected. This is a counterexample to the growing number of studies that have shown overlapping genetic patterns, demonstrating that genome scan comparisons can be noisy due to the effects of several interacting evolutionary forces.

Long-read sequencing reveals the evolutionary drivers of intra-host diversity across natural RNA mycovirus infections

Intra-host dynamics are a core component of virus evolution but most intra-host data come from a narrow range of hosts or experimental infections. Gaining broader information on the intra-host diversity and dynamics of naturally occurring virus infections is essential to our understanding of evolution across the virosphere. Here we used PacBio long-read HiFi sequencing to characterize the intra-host populations of natural infections of the RNA mycovirus Cryphonectria hypovirus 1 (CHV1). CHV1 is a biocontrol agent for the chestnut blight fungus (Cryphonectria parasitica), which co-invaded Europe alongside the fungus. We characterized the mutational and haplotypic intra-host virus diversity of thirty-eight natural CHV1 infections spread across four locations in Croatia and Switzerland. Intra-host CHV1 diversity values were shaped by purifying selection and accumulation of mutations over time as well as epistatic interactions within the host genome at defense loci. Geographical landscape features impacted CHV1 inter-host relationships through restricting dispersal and causing founder effects. Interestingly, a small number of intra-host viral haplotypes showed high sequence similarity across large geographical distances unlikely to be linked by dispersal.