Evolution and genetic characterization of Seoul virus in wild rats Rattus norvegicus from an urban park in Lyon, France 2020-2022.

BACKGROUND: Seoul virus (SEOV) is an orthohantavirus primarily carried by rats. In humans, it may cause hemorrhagic fever with renal syndrome (HFRS). Its incidence is likely underestimated and given the expansion of urban areas, a better knowledge of SEOV circulation in rat populations is called for. Beyond the need to improve human case detection, we need to deepen our comprehension of the ecological, epidemiological, and evolutionary processes involved in the transmission of SEOV. METHODOLOGY / PRINCIPAL FINDINGS: We performed a comprehensive serological and molecular characterization of SEOV in Rattus norvegicus in a popular urban park within a large city (Lyon, France) to provide essential information to design surveillance strategies regarding SEOV. We sampled rats within the urban park of 'La Tête d'Or' in Lyon city from 2020 to 2022. We combined rat population genetics, immunofluorescence assays, SEOV high-throughput sequencing (S, M, and L segments), and phylogenetic analyses. We found low structuring of wild rat populations within Lyon city. Only one sampling site within the park (building created in 2021) showed high genetic differentiation and deserves further attention. We confirmed the circulation of SEOV in rats from the park with high seroprevalence (17.2%) and high genetic similarity with the strain previously described in 2011 in Lyon city. CONCLUSION/SIGNIFICANCE: This study confirms the continuous circulation of SEOV in a popular urban park where the risk for SEOV transmission to humans is present. Implementing a surveillance of this virus could provide an efficient early warning system and help prepare risk-based interventions. As we reveal high gene flow between rat populations from the park and the rest of the city, we advocate for SEOV surveillance to be conducted at the scale of the entire city.

Characterizing viral within-host diversity in fast and non-equilibrium demo-genetic dynamics.

High-throughput sequencing has opened the route for a deep assessment of within-host genetic diversity that can be used, e.g., to characterize microbial communities and to infer transmission links in infectious disease outbreaks. The performance of such characterizations and inferences cannot be analytically assessed in general and are often grounded on computer-intensive evaluations. Then, being able to simulate within-host genetic diversity across time under various demo-genetic assumptions is paramount to assess the performance of the approaches of interest. In this context, we built an original model that can be simulated to investigate the temporal evolution of genotypes and their frequencies under various demo-genetic assumptions. The model describes the growth and the mutation of genotypes at the nucleotide resolution conditional on an overall within-host viral kinetics, and can be tuned to generate fast non-equilibrium demo-genetic dynamics. We ran simulations of this model and computed classic diversity indices to characterize the temporal variation of within-host genetic diversity (from high-throughput amplicon sequences) of virus populations under three demographic kinetic models of viral infection. Our results highlight how demographic (viral load) and genetic (mutation, selection, or drift) factors drive variations in within-host diversity during the course of an infection. In particular, we observed a non-monotonic relationship between pathogen population size and genetic diversity, and a reduction of the impact of mutation on diversity when a non-specific host immune response is activated. The large variation in the diversity patterns generated in our simulations suggests that the underlying model provides a flexible basis to produce very diverse demo-genetic scenarios and test, for instance, methods for the inference of transmission links during outbreaks.

A hotspot of Toxoplasma gondii Africa 1 lineage in Benin: How new genotypes from West Africa contribute to understand the parasite genetic diversity worldwide.

Through international trades, Europe, Africa and South America share a long history of exchanges, potentially of pathogens. We used the worldwide parasite Toxoplasma gondii to test the hypothesis of a historical influence on pathogen genetic diversity in Benin, a West African country with a longstanding sea trade history. In Africa, T. gondii spatial structure is still non-uniformly studied and very few articles have reported strain genetic diversity in fauna and clinical forms of human toxoplasmosis so far, even in African diaspora. Sera from 758 domestic animals (mainly poultry) in two coastal areas (Cotonou and Ouidah) and two inland areas (Parakou and Natitingou) were tested for T. gondii antibodies using a Modified Agglutination Test (MAT). The hearts and brains of 69 seropositive animals were collected for parasite isolation in a mouse bioassay. Forty-five strains were obtained and 39 genotypes could be described via 15-microsatellite genotyping, with a predominance of the autochthonous African lineage Africa 1 (36/39). The remaining genotypes were Africa 4 variant TUB2 (1/39) and two identical isolates (clone) of Type III (2/39). No difference in terms of genotype distribution between inland and coastal sampling sites was found. In particular, contrarily to what has been described in Senegal, no type II (mostly present in Europe) was isolated in poultry from coastal cities. This result seems to refute a possible role of European maritime trade in Benin despite it was one of the most important hubs during the slave trade period. However, the presence of the Africa 1 genotype in Brazil, predominant in Benin, and genetic analyses suggest that the triangular trade was a route for the intercontinental dissemination of genetic strains from Africa to South America. This supports the possibility of contamination in humans and animals with potentially imported virulent strains.

Distribution and evolution of the major viruses infecting cucurbitaceous and solanaceous crops in the French Mediterranean area.

Plant viral diseases represent a significant burden to plant health, and their highest impact in Mediterranean agriculture is on vegetables grown under intensive horticultural practices. In order to understand better virus evolution and emergence, the most prevalent viruses were mapped in the main cucurbitaceous (melon, squashes) and solanaceous (tomato, pepper) crops and in some wild hosts in the French Mediterranean area, and virus diversity, evolution and population structure were studied through molecular epidemiology approaches. Surveys were performed in summer 2016 and 2017, representing a total of 1530 crop samples and 280 weed samples. The plant samples were analysed using serological and molecular approaches, including high-throughput sequencing (HTS). The viral species and their frequency in crops were quite similar to those of surveys conducted ten years before in the same areas. Contrary to other Mediterranean countries, aphid-transmitted viruses remain the most prevalent in France whereas whitefly-transmitted ones have not yet emerged. However, HTS analysis of viral evolution revealed the appearance of undescribed viral variants, especially for watermelon mosaic virus (WMV) in cucurbits, or variants not present in France before, as for cucumber mosaic virus (CMV) in solanaceous crops. Deep sequencing also revealed complex virus populations within individual plants with frequent recombination or reassortment. The spatial genetic structure of cucurbit aphid-borne yellows virus (CABYV) was related to the landscape structure, whereas in the case of WMV, the recurrence of introduction events and probable human exchanges of plant material resulted in complex spatial pattern of genetic variation.

Exploiting Genetic Information to Trace Plant Virus Dispersal in Landscapes.

During the past decade, knowledge of pathogen life history has greatly benefited from the advent and development of molecular epidemiology. This branch of epidemiology uses information on pathogen variation at the molecular level to gain insights into a pathogen's niche and evolution and to characterize pathogen dispersal within and between host populations. Here, we review molecular epidemiology approaches that have been developed to trace plant virus dispersal in landscapes. In particular, we highlight how virus molecular epidemiology, nourished with powerful sequencing technologies, can provide novel insights at the crossroads between the blooming fields of landscape genetics, phylogeography, and evolutionary epidemiology. We present existing approaches and their limitations and contributions to the understanding of plant virus epidemiology.

Modeling spatiotemporal dynamics of outbreaking species: influence of environment and migration in a locust.

Many pest species exhibit huge fluctuations in population abundance. Understanding their large-scale and long-term dynamics is necessary to develop effective control and management strategies. Occupancy models represent a promising approach to unravel interactions between environmental factors and spatiotemporal dynamics of outbreaking populations. Here, we investigated population dynamics of the Australian plague locust, Chortoicetes terminifera, using density data collected between 1988 and 2010 by the Australian Plague Locust Commission over more than 3 million km2 in eastern Australia. We applied multistate and autologistic multi-season occupancy models to test competing hypotheses about environmental and demographic processes affecting the large-scale dynamics of the Australian plague locust. We found that rainfall and land cover predictors best explained the spatial variability in outbreak probability across eastern Australia. Outbreaks are more likely to occur in temperate than tropical regions, with a faster and more continuous response to rainfall in desert than in agricultural areas. Our results also support the hypothesis that migration tends to propagate outbreaks only locally (over distances lower than 400 km) rather than across climatic regions. Our study suggests that locust outbreak forecasting and management systems could be improved by implementing key environmental factors and migration in hierarchical spatial models. Finally, our modeling framework can be seen as a step towards bridging the gap between mechanistic and more phenomenological models in the spatial analysis of fluctuating populations.

Short-term variations in gene flow related to cyclic density fluctuations in the common vole.

In highly fluctuating populations with complex social systems, genetic patterns are likely to vary in space and time due to demographic and behavioural processes. Cyclic rodents are extreme examples of demographically instable populations that often exhibit strong social organization. In such populations, kin structure and spacing behaviour may vary with density fluctuations and impact both the composition and spatial structure of genetic diversity. In this study, we analysed the multiannual genetic structure of a cyclic rodent, Microtus arvalis, using a sample of 875 individuals trapped over three complete cycles (from 1999 to 2007) and genotyped at 10 microsatellite loci. We tested the predictions that genetic diversity and gene flow intensity vary with density fluctuations. We found evidences for both spatial scale-dependant variations in genetic diversity and higher gene flow during high density. Moreover, investigation of sex-specific relatedness patterns revealed that, although dispersal is biased toward males in this species, distances moved by both sexes were lengthened during high density. Altogether, these results suggest that an increase in migration with density allows to restore the local loss of genetic diversity occurring during low density. We then postulate that this change in migration results from local competition, which enhances female colonization of empty spaces and male dispersal among colonies.

Dispersal, landscape and travelling waves in cyclic vole populations.

Travelling waves (TW) are among the most striking ecological phenomena emerging in oscillating populations. Despite much theory, understanding how real-world TW arise remains a challenge for ecology. Herein, we analyse 16-year time series of cyclic vole populations collected at 314 localities covering 2500 km² in France. We found evidence for a linear front TW spreading at a speed of 7.4 km year(-1) along a north-west/south-east direction and radiating away from a major landscape discontinuity as predicted by recent theory. The spatial signature of vole dispersal was assessed using genetic data collected at 14 localities. Both data sets were handled using similar autocorrelation approaches. Our results revealed a remarkable congruence of the spatial extent and direction of anisotropy of both demographic and genetic structures. Our results constitute the first empirical evidence that effective dispersal is limited in the direction of TW while most of the individual exchanges occur along the wave front.

Massive nest-box supplementation boosts fecundity, survival and even immigration without altering mating and reproductive behaviour in a rapidly recovered bird population.

Habitat restoration measures may result in artificially high breeding density, for instance when nest-boxes saturate the environment, which can negatively impact species' demography. Potential risks include changes in mating and reproductive behaviour such as increased extra-pair paternity, conspecific brood parasitism, and polygyny. Under particular cicumstances, these mechanisms may disrupt reproduction, with populations dragged into an extinction vortex. With the use of nuclear microsatellite markers, we investigated the occurrence of these potentially negative effects in a recovered population of a rare secondary cavity-nesting farmland bird of Central Europe, the hoopoe (Upupa epops). High intensity farming in the study area has resulted in a total eradication of cavity trees, depriving hoopoes from breeding sites. An intensive nest-box campaign rectified this problem, resulting in a spectacular population recovery within a few years only. There was some concern, however, that the new, high artificially-induced breeding density might alter hoopoe mating and reproductive behaviour. As the species underwent a serious demographic bottleneck in the 1970-1990s, we also used the microsatellite markers to reconstitute the demo-genetic history of the population, looking in particular for signs of genetic erosion. We found i) a low occurrence of extra-pair paternity, polygyny and conspecific brood parasitism, ii) a high level of neutral genetic diversity (mean number of alleles and expected heterozygosity per locus: 13.8 and 83%, respectively) and, iii) evidence for genetic connectivity through recent immigration of individuals from well differentiated populations. The recent increase in breeding density did thus not induce so far any noticeable detrimental changes in mating and reproductive behaviour. The demographic bottleneck undergone by the population in the 1970s-1990s was furthermore not accompanied by any significant drop in neutral genetic diversity. Finally, genetic data converged with a concomitant demographic study to evidence that immigration strongly contributed to local population recovery.

Challenges to assessing connectivity between massive populations of the Australian plague locust.

Linking demographic and genetic dispersal measures is of fundamental importance for movement ecology and evolution. However, such integration can be difficult, particularly for highly fecund species that are often the target of management decisions guided by an understanding of population movement. Here, we present an example of how the influence of large population sizes can preclude genetic approaches from assessing demographic population structuring, even at a continental scale. The Australian plague locust, Chortoicetes terminifera, is a significant pest, with populations on the eastern and western sides of Australia having been monitored and managed independently to date. We used microsatellites to assess genetic variation in 12 C. terminifera population samples separated by up to 3000 km. Traditional summary statistics indicated high levels of genetic diversity and a surprising lack of population structure across the entire range. An approximate Bayesian computation treatment indicated that levels of genetic diversity in C. terminifera corresponded to effective population sizes conservatively composed of tens of thousands to several million individuals. We used these estimates and computer simulations to estimate the minimum rate of dispersal, m, that could account for the observed range-wide genetic homogeneity. The rate of dispersal between both sides of the Australian continent could be several orders of magnitude lower than that typically considered as required for the demographic connectivity of populations.

Crosses prior to parthenogenesis explain the current genetic diversity of tropical plant-parasitic Meloidogyne species (Nematoda: Tylenchida).

The tropical and subtropical parthenogenetic plant-parasitic nematodes Meloidogyne are polyphagous major agricultural pests. Implementing proper pest management approaches requires a good understanding of mechanisms, population structure, evolutionary patterns and species identification. A comparative analysis of the mitochondrial vs nuclear diversity was conducted on a selected set of Meloidogyne lines from various geographic origins. Mitochondrial co2-16S sequences and AFLP markers of total DNA were applied because of their ability to evidence discrete genetic variation between closely related isolates. Several distinct maternal lineages were present, now associated with different genetic backgrounds. Relative discordances were found when comparing mitochondrial and nuclear diversity patterns. These patterns are most likely related to crosses within one ancestral genetic pool, followed by the establishment of parthenogenesis. In this case, they mirror the genetic backgrounds of the original individuals. Another aspect could be that species emergence was recent or on process from this original genetic pool and that the relatively short time elapsed since then and before parthenogenesis settlement did not allow for lineage sorting. This could also be compatible with the hypothesis of hybrids between closely related species. This genetic pool would correspond to a species as defined by the species interbreeding concept, but also including the grey area of species boundaries. This complex process has implications on the way genotypic and phenotypic diversity should be addressed. The phenotype of parthenogenetic lines is at least for part determined by the ancestral amphimictic genetic background. A direct consequence is, therefore, in terms of risk management, the limited confidence one can have on the direct association of an agronomic threat to a simple typing or species delineation. Risk management strategies and tools must thus consider this complexity when designing quarantine implementation, resistance breeding programmes or molecular diagnostic.

Stress and demographic decline: a potential effect mediated by impairment of reproduction and immune function in cyclic vole populations.

The stress response is initially adaptive, operating to maintain homeostasis. However, chronic long-term exposure to stressors may have detrimental effects. We proposed that chronic stress may be a major factor in demographic vole cycles, inducing decline in high-density populations. We monitored four populations of the fossorial water vole Arvicola scherman, which undergo pluriannual demographic cycles in the Jura Mountains (France). Sampling was conducted during the high densities and the decline. We measured fecal corticosterone metabolites (FCMs) to assess stress levels and injected phytohemagglutinin to estimate the cell-mediated immune response. We demonstrated that stress levels increase between the high densities and the decline in most of the vole populations. At the individual level, FCM concentrations varied with reproductive status and body condition. During the outbreak, we observed significantly higher levels of FCM concentrations in nulliparous females than in females that had previously reproduced. Moreover, a significant negative correlation was observed between concentrations of FCMs and both immunocompetence and body condition during population decline. These results might reflect an impairment of the female reproductive capability in high densities and accelerated senescence affecting immune function during decline, both arising from chronic stress.

Modeling survival and mark loss in molting animals: recapture, dead recoveries, and exuvia recoveries.

Capture-mark-recapture (CMR) analyses aim primarily at estimating relevant life history parameters, despite the fact that some individuals are not always recaptured, even if alive on the study site. Applying such approaches to species with a complex life cycle, such as insects, remains challenging because each change of stage tends to cause mark loss through molting. We developed a multistate model based on three exclusive events ("dead", "surviving and molting", and "surviving and staying in the same larval stage") to estimate probabilities of survival and mark loss. Estimates of biologically relevant parameters were derived from those of the probabilities of transition between these states. The model was applied to data from radio-tracking diodes glued on grasshoppers. The estimates of recapture probabilities decreased throughout the season for animals remaining alive, while the detection of dead animals and lost diodes was exhaustive. The survival probability was higher for larvae than for adults (0.98 vs. 0.96), and mark loss was stronger in larvae than in adults (0.09 vs. 0.06). We show that the survival rate of a species with a high rate of mark loss can be estimated using multistate models, provided that marks can be recovered after being lost. These models are flexible enough to test for several effects that potentially affect survival and mark loss probabilities.

Linking demography and host dispersal to Trichuris arvicolae distribution in a cyclic vole species.

Spatial structure in the distribution of pathogen infection can influence both epidemiology and host-parasite coevolutionary processes. It may result from the spatial heterogeneity of intrinsic and extrinsic factors, or from the local population dynamics of hosts and parasites. In this study, we investigated the effects of landscape, host dispersal and demography (population abundance and phase of the fluctuation) on the distribution of a gastro-intestinal nematode Trichuris arvicolae in the fossorial water vole Arvicola terrestris sherman. This rodent exhibits outbreaks occurring regularly in Franche-Comté (France). Thirteen out-of-phase populations were studied in autumn 2003. They exhibited highly different T. arvicolae prevalences. The heterogeneity in prevalences was not explained by population structure, landscape or vole abundance, but by the phase of the vole population fluctuations. Populations at the end of the high density phase showed null prevalence whereas populations in increase or outbreak phases exhibited higher prevalences. Population genetic analyses based on microsatellites revealed significant differentiation between vole populations, and higher dispersal rates of young voles compared with old ones. These younger individuals were also infected more frequently than older voles. This suggested a role of host dispersal in the distribution of T. arvicolae. However, there was a strong discrepancy between the spatial patterns of prevalence and of host genetics or demographic phase. Genetic differentiation and differences in demographic phase exhibited significant spatial autocorrelations whereas prevalence did not. We concluded that the distribution of T. arvicolae is influenced by vole dispersal, although this effect might be overwhelmed by local adaptation processes or environmental conditions.