Emergence and clonal expansion of Vibrio aestuarianus lineages pathogenic for oysters in Europe.

Crassostrea gigas oysters represent a significant global food source with 4.7 million tons harvested per year. In 2001, the bacterium Vibrio aestuarianus subsp. francensis emerged as a pathogen that causes adult oyster mortality in France and Ireland. Its impact on oyster aquaculture has increased in Europe since its re-emergence in 2012. To better understand the evolutionary mechanisms leading to the emergence and persistence over time of this pathogen, we conducted a survey of mollusc diseases through national reference laboratories across Europe. We analysed 54 new genomes of Vibrio aestuarianus (Va) isolated from multiple environmental compartments since 2001, in areas with and without bivalve mortalities. We used a combination of comparative genomics and population genetics approaches and show that Va has a classical epidemic population structure from which the pathogenic Va francensis subspecies emerged and clonally expanded. Furthermore, we identified a specific cus-cop-containing island conferring copper resistance to Va francensis whose acquisition may have favoured the emergence of pathogenic lineages adapted and specialized to oysters.

Geographic Range Overlap Rather than Phylogenetic Distance Explains Rabies Virus Transmission among Closely Related Bat Species.

The cross-species transmission (CST) of pathogens can have dramatic consequences, as highlighted by recent disease emergence events affecting human, animal and plant health. Understanding the ecological and evolutionary factors that increase the likelihood of disease agents infecting and establishing in a novel host is therefore an important research area. Previous work across different pathogens, including rabies virus (RABV), found that increased evolutionary distance between hosts reduces the frequency of cross-species transmission and of permanent host shifts. However, whether this effect of host relatedness still holds for transmission among recently diverged hosts is not well understood. We aimed to ask if high host relatedness can still increase the probability of a host shift between more recently diverged hosts, and the importance of this effect relative to ecological predictors. We first addressed this question by quantifying the CST frequency of RABV between North American bat species within the genus Myotis, using a multi-decade data set containing 128 nucleoprotein (N) RABV sequences from ten host species. We compared RABV CST frequency within Myotis to the rates of CST between nine genera of North American bat species. We then examined whether host relatedness or host range overlap better explains the frequency of CST seen between Myotis species. We found that at the within genus scale, host range overlap, rather than host relatedness best explains the frequency of CST events. Moreover, we found evidence of CST occurring among a higher proportion of species, and CST more frequently resulting in sustained transmission in the novel host in the Myotis dataset compared to the multi-genus dataset. Our results suggest that among recently diverged species, the ability to infect a novel host is no longer restricted by physiological barriers but instead is limited by physical contact. Our results improve predictions of where future CST events for RABV might occur and clarify the relationship between host divergence and pathogen emergence.

Evaluation of tangential flow filtration coupled to long-read sequencing for ostreid herpesvirus type 1 genome assembly.

Whole-genome sequencing is widely used to better understand the transmission dynamics, the evolution and the emergence of new variants of viral pathogens. This can bring crucial information to stakeholders for disease management. Unfortunately, aquatic virus genomes are usually difficult to characterize because most of these viruses cannot be easily propagated in vitro. Developing methodologies for routine genome sequencing of aquatic viruses is timely given the ongoing threat of disease emergence. This is particularly true for pathogenic viruses infecting species of commercial interest that are widely exchanged between production basins or countries. For example, the ostreid herpesvirus type 1 (OsHV-1) is a Herpesvirus widely associated with mass mortality events of juvenile Pacific oyster Crassostrea gigas. Genomes of Herpesviruses are large and complex with long direct and inverted terminal repeats. In addition, OsHV-1 is unculturable. It therefore accumulates several features that make its genome sequencing and assembly challenging. To overcome these difficulties, we developed a tangential flow filtration (TFF) method to enrich OsHV-1 infective particles from infected host tissues. This virus purification allowed us to extract high molecular weight and high-quality viral DNA that was subjected to Illumina short-read and Nanopore long-read sequencing. Dedicated bioinformatic pipelines were developed to assemble complete OsHV-1 genomes with reads from both sequencing technologies. Nanopore sequencing allowed characterization of new structural variations and major viral isomers while having 99,98 % of nucleotide identity with the Illumina assembled genome. Our study shows that TFF-based purification method, coupled with Nanopore sequencing, is a promising approach to enable in field sequencing of unculturable aquatic DNA virus.

Meteorological and climatic variables predict the phenology of Ixodes ricinus nymph activity in France, accounting for habitat heterogeneity.

Ixodes ricinus ticks (Acari: Ixodidae) are the most important vector for Lyme borreliosis in Europe. As climate change might affect their distributions and activities, this study aimed to determine the effects of environmental factors, i.e., meteorological, bioclimatic, and habitat characteristics on host-seeking (questing) activity of I. ricinus nymphs, an important stage in disease transmissions, across diverse climatic types in France over 8 years. Questing activity was observed using a repeated removal sampling with a cloth-dragging technique in 11 sampling sites from 7 tick observatories from 2014 to 2021 at approximately 1-month intervals, involving 631 sampling campaigns. Three phenological patterns were observed, potentially following a climatic gradient. The mixed-effects negative binomial regression revealed that observed nymph counts were driven by different interval-average meteorological variables, including 1-month moving average temperature, previous 3-to-6-month moving average temperature, and 6-month moving average minimum relative humidity. The interaction effects indicated that the phenology in colder climates peaked differently from that of warmer climates. Also, land cover characteristics that support the highest baseline abundance were moderate forest fragmentation with transition borders with agricultural areas. Finally, our model could potentially be used to predict seasonal human-tick exposure risks in France that could contribute to mitigating Lyme borreliosis risk.

Genetic diversity and connectivity of the Ostreid herpesvirus 1 populations in France: A first attempt to phylogeographic inference for a marine mollusc disease.

The genetic diversity of viral populations is a key driver of the spatial and temporal diffusion of viruses; yet, studying the diversity of whole genomes from natural populations still remains a challenge. Phylodynamic approaches are commonly used for RNA viruses harboring small genomes but have only rarely been applied to DNA viruses with larger genomes. Here, we used the Pacific oyster mortality syndrome (a disease that affects oyster farms around the world) as a model to study the genetic diversity of its causative agent, the Ostreid herpesvirus 1 (OsHV-1) in the three main French oyster-farming areas. Using ultra-deep sequencing on individual moribund oysters and an innovative combination of bioinformatics tools, we de novo assembled twenty-one OsHV-1 new genomes. Combining quantification of major and minor genetic variations, phylogenetic analysis, and ancestral state reconstruction of discrete traits approaches, we assessed the connectivity of OsHV-1 viral populations between the three oyster-farming areas. Our results suggest that the Marennes-Oléron Bay represents the main source of OsHV-1 diversity, from where the virus has dispersed to other farming areas, a scenario consistent with current practices of oyster transfers in France. We demonstrate that phylodynamic approaches can be applied to aquatic DNA viruses to determine how epidemiological, immunological, and evolutionary processes act and potentially interact to shape their diversity patterns.

Genomic Diversity of the Ostreid Herpesvirus Type 1 Across Time and Location and Among Host Species.

The mechanisms underlying virus emergence are rarely well understood, making the appearance of outbreaks largely unpredictable. This is particularly true for pathogens with low per-site mutation rates, such as DNA viruses, that do not exhibit a large amount of evolutionary change among genetic sequences sampled at different time points. However, whole-genome sequencing can reveal the accumulation of novel genetic variation between samples, promising to render most, if not all, microbial pathogens measurably evolving and suitable for analytical techniques derived from population genetic theory. Here, we aim to assess the measurability of evolution on epidemiological time scales of the Ostreid herpesvirus 1 (OsHV-1), a double stranded DNA virus of which a new variant, OsHV-1 µVar, emerged in France in 2008, spreading across Europe and causing dramatic economic and ecological damage. We performed phylogenetic analyses of heterochronous (n = 21) OsHV-1 genomes sampled worldwide. Results show sufficient temporal signal in the viral sequences to proceed with phylogenetic molecular clock analyses and they indicate that the genetic diversity seen in these OsHV-1 isolates has arisen within the past three decades. OsHV-1 samples from France and New Zealand did not cluster together suggesting a spatial structuration of the viral populations. The genome-wide study of simple and complex polymorphisms shows that specific genomic regions are deleted in several isolates or accumulate a high number of substitutions. These contrasting and non-random patterns of polymorphism suggest that some genomic regions are affected by strong selective pressures. Interestingly, we also found variant genotypes within all infected individuals. Altogether, these results provide baseline evidence that whole genome sequencing could be used to study population dynamic processes of OsHV-1, and more broadly herpesviruses.

High dispersal capacity of Culicoides obsoletus (Diptera: Ceratopogonidae), vector of bluetongue and Schmallenberg viruses, revealed by landscape genetic analyses.

BACKGROUND: In the last two decades, recurrent epizootics of bluetongue virus and Schmallenberg virus have been reported in the western Palearctic region. These viruses affect domestic cattle, sheep, goats and wild ruminants and are transmitted by native hematophagous midges of the genus Culicoides (Diptera: Ceratopogonidae). Culicoides dispersal is known to be stratified, i.e. due to a combination of dispersal processes occurring actively at short distances and passively or semi-actively at long distances, allowing individuals to jump hundreds of kilometers. METHODS: Here, we aim to identify the environmental factors that promote or limit gene flow of Culicoides obsoletus, an abundant and widespread vector species in Europe, using an innovative framework integrating spatial, population genetics and statistical approaches. A total of 348 individuals were sampled in 46 sites in France and were genotyped using 13 newly designed microsatellite markers. RESULTS: We found low genetic differentiation and a weak population structure for C. obsoletus across the country. Using three complementary inter-individual genetic distances, we did not detect any significant isolation by distance, but did detect significant anisotropic isolation by distance on a north-south axis. We employed a multiple regression on distance matrices approach to investigate the correlation between genetic and environmental distances. Among all the environmental factors that were tested, only cattle density seems to have an impact on C. obsoletus gene flow. CONCLUSIONS: The high dispersal capacity of C. obsoletus over land found in the present study calls for a re-evaluation of the impact of Culicoides on virus dispersal, and highlights the urgent need to better integrate molecular, spatial and statistical information to guide vector-borne disease control.

"Frozen evolution" of an RNA virus suggests accidental release as a potential cause of arbovirus re-emergence.

The mechanisms underlying virus emergence are rarely well understood, making the appearance of outbreaks largely unpredictable. Bluetongue virus serotype 8 (BTV-8), an arthropod-borne virus of ruminants, emerged in livestock in northern Europe in 2006, spreading to most European countries by 2009 and causing losses of billions of euros. Although the outbreak was successfully controlled through vaccination by early 2010, puzzlingly, a closely related BTV-8 strain re-emerged in France in 2015, triggering a second outbreak that is still ongoing. The origin of this virus and the mechanisms underlying its re-emergence are unknown. Here, we performed phylogenetic analyses of 164 whole BTV-8 genomes sampled throughout the two outbreaks. We demonstrate consistent clock-like virus evolution during both epizootics but found negligible evolutionary change between them. We estimate that the ancestor of the second outbreak dates from the height of the first outbreak in 2008. This implies that the virus had not been replicating for multiple years prior to its re-emergence in 2015. Given the absence of any known natural mechanism that could explain BTV-8 persistence over this long period without replication, we hypothesise that the second outbreak could have been initiated by accidental exposure of livestock to frozen material contaminated with virus from approximately 2008. Our work highlights new targets for pathogen surveillance programmes in livestock and illustrates the power of genomic epidemiology to identify pathways of infectious disease emergence.

Contrasting selective patterns across the segmented genome of bluetongue virus in a global reassortment hotspot.

For segmented viruses, rapid genomic and phenotypic changes can occur through the process of reassortment, whereby co-infecting strains exchange entire segments creating novel progeny virus genotypes. However, for many viruses with segmented genomes, this process and its effect on transmission dynamics remain poorly understood. Here, we assessed the consequences of reassortment for selection on viral diversity through time using bluetongue virus (BTV), a segmented arbovirus that is the causative agent of a major disease of ruminants. We analysed ninety-two BTV genomes isolated across four decades from India, where BTV diversity, and thus opportunities for reassortment, are among the highest in the world. Our results point to frequent reassortment and segment turnover, some of which appear to be driven by selective sweeps and serial hitchhiking. Particularly, we found evidence for a recent selective sweep affecting segment 5 and its encoded NS1 protein that has allowed a single variant to essentially invade the full range of BTV genomic backgrounds and serotypes currently circulating in India. In contrast, diversifying selection was found to play an important role in maintaining genetic diversity in genes encoding outer surface proteins involved in virus interactions (VP2 and VP5, encoded by segments 2 and 6, respectively). Our results support the role of reassortment in driving rapid phenotypic change in segmented viruses and generate testable hypotheses for in vitro experiments aiming at understanding the specific mechanisms underlying differences in fitness and selection across viral genomes.

Endemic infection can shape exposure to novel pathogens: Pathogen co-occurrence networks in the Serengeti lions.

Pathogens are embedded in a complex network of microparasites that can collectively or individually alter disease dynamics and outcomes. Endemic pathogens that infect an individual in the first years of life, for example, can either facilitate or compete with subsequent pathogens thereby exacerbating or ameliorating morbidity and mortality. Pathogen associations are ubiquitous but poorly understood, particularly in wild populations. We report here on 10 years of serological and molecular data in African lions, leveraging comprehensive demographic and behavioural data to test if endemic pathogens shape subsequent infection by epidemic pathogens. We combine network and community ecology approaches to assess broad network structure and characterise associations between pathogens across spatial and temporal scales. We found significant non-random structure in the lion-pathogen co-occurrence network and identified both positive and negative associations between endemic and epidemic pathogens. Our results provide novel insights on the complex associations underlying pathogen co-occurrence networks.

Bluetongue virus spread in Europe is a consequence of climatic, landscape and vertebrate host factors as revealed by phylogeographic inference.

Spatio-temporal patterns of the spread of infectious diseases are commonly driven by environmental and ecological factors. This is particularly true for vector-borne diseases because vector populations can be strongly affected by host distribution as well as by climatic and landscape variables. Here, we aim to identify environmental drivers for bluetongue virus (BTV), the causative agent of a major vector-borne disease of ruminants that has emerged multiple times in Europe in recent decades. In order to determine the importance of climatic, landscape and host-related factors affecting BTV diffusion across Europe, we fitted different phylogeographic models to a dataset of 113 time-stamped and geo-referenced BTV genomes, representing multiple strains and serotypes. Diffusion models using continuous space revealed that terrestrial habitat below 300 m altitude, wind direction and higher livestock densities were associated with faster BTV movement. Results of discrete phylogeographic analysis involving generalized linear models broadly supported these findings, but varied considerably with the level of spatial partitioning. Contrary to common perception, we found no evidence for average temperature having a positive effect on BTV diffusion, though both methodological and biological reasons could be responsible for this result. Our study provides important insights into the drivers of BTV transmission at the landscape scale that could inform predictive models of viral spread and have implications for designing control strategies.

Linking social and spatial networks to viral community phylogenetics reveals subtype-specific transmission dynamics in African lions.

Heterogeneity within pathogen species can have important consequences for how pathogens transmit across landscapes; however, discerning different transmission routes is challenging. Here, we apply both phylodynamic and phylogenetic community ecology techniques to examine the consequences of pathogen heterogeneity on transmission by assessing subtype-specific transmission pathways in a social carnivore. We use comprehensive social and spatial network data to examine transmission pathways for three subtypes of feline immunodeficiency virus (FIVPle ) in African lions (Panthera leo) at multiple scales in the Serengeti National Park, Tanzania. We used FIVPle molecular data to examine the role of social organization and lion density in shaping transmission pathways and tested to what extent vertical (i.e., father- and/or mother-offspring relationships) or horizontal (between unrelated individuals) transmission underpinned these patterns for each subtype. Using the same data, we constructed subtype-specific FIVPle co-occurrence networks and assessed what combination of social networks, spatial networks or co-infection best structured the FIVPle network. While social organization (i.e., pride) was an important component of FIVPle transmission pathways at all scales, we find that FIVPle subtypes exhibited different transmission pathways at within- and between-pride scales. A combination of social and spatial networks, coupled with consideration of subtype co-infection, was likely to be important for FIVPle transmission for the two major subtypes, but the relative contribution of each factor was strongly subtype-specific. Our study provides evidence that pathogen heterogeneity is important in understanding pathogen transmission, which could have consequences for how endemic pathogens are managed. Furthermore, we demonstrate that community phylogenetic ecology coupled with phylodynamic techniques can reveal insights into the differential evolutionary pressures acting on virus subtypes, which can manifest into landscape-level effects.

Multiple independent transmission cycles of a tick-borne pathogen within a local host community.

Many pathogens are maintained by multiple host species and involve multiple strains with potentially different phenotypic characteristics. Disentangling transmission patterns in such systems is often challenging, yet investigating how different host species contribute to transmission is crucial to properly assess and manage disease risk. We aim to reveal transmission cycles of bacteria within the Borrelia burgdorferi species complex, which include Lyme disease agents. We characterized Borrelia genotypes found in 488 infected Ixodes ricinus nymphs collected in the Sénart Forest located near Paris (France). These genotypes were compared to those observed in three sympatric species of small mammals and network analyses reveal four independent transmission cycles. Statistical modelling shows that two cycles involving chipmunks, an introduced species, and non-sampled host species such as birds, are responsible for the majority of tick infections. In contrast, the cycle involving native bank voles only accounts for a small proportion of infected ticks. Genotypes associated with the two primary transmission cycles were isolated from Lyme disease patients, confirming the epidemiological threat posed by these strains. Our work demonstrates that combining high-throughput sequence typing with networks tools and statistical modeling is a promising approach for characterizing transmission cycles of multi-host pathogens in complex ecological settings.

Comparative population genomics of the Borrelia burgdorferi species complex reveals high degree of genetic isolation among species and underscores benefits and constraints to studying intra-specific epidemiological processes.

Lyme borreliosis, one of the most frequently contracted zoonotic diseases in the Northern Hemisphere, is caused by bacteria belonging to different genetic groups within the Borrelia burgdorferi species complex, which are transmitted by ticks among various wildlife reservoirs, such as small mammals and birds. These features make the Borrelia burgdorferi species complex an attractive biological model that can be used to study the diversification and the epidemiology of endemic bacterial pathogens. We investigated the potential of population genomic approaches to study these processes. Sixty-three strains belonging to three species within the Borrelia burgdorferi complex were isolated from questing ticks in Alsace (France), a region where Lyme disease is highly endemic. We first aimed to characterize the degree of genetic isolation among the species sampled. Phylogenetic and coalescent-based analyses revealed clear delineations: there was a ∼50 fold difference between intra-specific and inter-specific recombination rates. We then investigated whether the population genomic data contained information of epidemiological relevance. In phylogenies inferred using most of the genome, conspecific strains did not cluster in clades. These results raise questions about the relevance of different strategies when investigating pathogen epidemiology. For instance, here, both classical analytic approaches and phylodynamic simulations suggested that population sizes and migration rates were higher in B. garinii populations, which are normally associated with birds, than in B. burgdorferi s.s. populations. The phylogenetic analyses of the infection-related ospC gene and its flanking region provided additional support for this finding. Traces of recombination among the B. burgdorferi s.s. lineages and lineages associated with small mammals were found, suggesting that they shared the same hosts. Altogether, these results provide baseline evidence that can be used to formulate hypotheses regarding the host range of B. burgdorferi lineages based on population genomic data.

High-throughput sequence typing reveals genetic differentiation and host specialization among populations of the Borrelia burgdorferi species complex that infect rodents.

Lyme disease is a zoonosis caused by various species belonging to the Borrelia burgdorferi bacterial species complex. These pathogens are transmitted by ticks and infect multiple, taxonomically distinct, host species. From an epidemiological perspective, it is important to determine whether genetic variants within the species complex are able to spread freely through the whole host community or, instead, if certain variants are restricted to particular hosts. To this end, we characterized the genotypes of members of the B. burgdorferi species complex; the bacteria were isolated from more than two hundred individuals captured in the wild and belonging to three different rodent host species. For each individual, we used a high-throughput approach to amplify and sequence rplB, a housekeeping gene, and ospC, which is involved in infection. This approach allowed us to evaluate the genetic diversity both within and among species in the B. burgdorferi species complex. Strong evidence of genetic differentiation among host species was revealed by both genes, even though they are, a priori, not constrained by the same selective pressures. These data are discussed in the context of the advancements made possible by multi-locus high-throughput sequencing and current knowledge of Lyme disease epidemiology.