The neutral effective migration rate in a mainland-island context

Complete mitochondrial exploration of Echinococcus multilocularis from French alveolar echinococcosis patients

Alveolar echinococcosis (AE) is a parasitosis that is expanding worldwide, including in Europe. The development of genotypic markers is essential to follow its spatiotemporal evolution. Sequencing of the commonly used mitochondrial genes cob, cox1, and nad2 shows low discriminatory power, and analysis of the microsatellite marker EmsB does not allow nucleotide sequence analysis. We aimed to develop a new method for the genotyping of Echinococcus multilocularis based on whole mitochondrial genome (mitogenome) sequencing, to determine the genetic diversity among 30 human visceral samples from French patients, and compare this method with those currently in use. Sequencing of the whole mitochondrial genome was carried out after amplification by PCR, using one uniplex and two multiplex reactions to cover the 13,738 bp of the mitogenome, combined with Illumina technology. Thirty complete mitogenome sequences were obtained from AE lesions. One showed strong identity with Asian genotypes (99.98% identity) in a patient who had travelled to China. The other 29 mitogenomes could be differentiated into 13 haplotypes, showing higher haplotype and nucleotide diversity than when using the cob, cox1, and nad2 gene sequences alone. The mitochondrial genotyping data and EmsB profiles did not overlap, probably because one method uses the mitochondrial genome and the other the nuclear genome. The pairwise fixation index (F_st) value between individuals living inside and those living outside the endemic area was high (F_st = 0.222, P = 0.002). This is consistent with the hypothesis of an expansion from historical endemic areas to peripheral regions.

Reproductive isolation via polygenic local adaptation in sub-divided populations: Effect of linkage disequilibria and drift

This paper considers how polygenic local adaptation and reproductive isolation between hybridizing populations is influenced by linkage disequilibria (LD) between loci, in scenarios where both gene flow and genetic drift counteract selection. It shows that the combined effects of multi-locus LD and genetic drift on allele frequencies at selected loci and on heterozygosity at neutral loci are predicted accurately by incorporating (deterministic) effective migration rates into the diffusion approximation (for selected loci) and into the structured coalescent (for neutral loci). Theoretical approximations are tested against individual-based simulations and used to investigate conditions for the maintenance of local adaptation on an island subject to one-way migration from a differently adapted mainland, and in an infinite-island population with two habitats under divergent selection. The analysis clarifies the conditions under which LD between sets of locally deleterious alleles allows these to be collectively eliminated despite drift, causing sharper and (under certain conditions) shifted migration thresholds for loss of adaptation. Local adaptation also has counter-intuitive effects on neutral (relative) divergence: F_ST is highest for a pair of subpopulations belonging to the same (rare) habitat, despite the lack of reproductive isolation between them.

Environmental adaptation often involves spatially heterogeneous selection at many genetic loci. Thus, the evolutionary consequences of hybridisation between populations adapted to different environments depend on the coupled dynamics of multiple loci under selection, migration and genetic drift, making them challenging to predict. Here, I introduce theoretical approximations that accurately capture the effect of such coupling on allele frequencies at individual loci, while also accounting for the stochastic effects of genetic drift. I then use these approximations to study hybridisation in a metapopulation consisting of many interconnected subpopulations, where each subpopulation belongs to one of two habitats under divergent selection. The analysis clarifies how subpopulations belonging to a rare habitat can maintain local adaptation despite high levels of migration if net selection against multi-locus genotypes is stronger than a threshold which depends on the relative abundances of the two habitats. Further, local adaptation in a metapopulation can significantly elevate F_ST between subpopulations belonging to the same habitat, even though these are not reproductively isolated. These findings highlight the importance of carefully considering the genetic architecture and spatial context of divergence when interpreting patterns of genomic differentiation between speciating populations.

[Spread and Development of Alveolar Echinococcosis in Germany, 1992-2018]

GOAL OF THE STUDY

Alveolar echinococcosis (AE) is a rare parasitosis caused by the pathogen Echinococcus multilocularis. There is an increase in the number of cases in many countries. The aim of the study was to investigate the current prevalence and the change in the geographical distribution pattern.

METHODOLOGY

Data were collected retrospectively for the period 1992-2018 using the registered cases in the national disease register for AE in Germany. Statistical analysis was performed using the SAS statistical analysis system version 9.4 (SAS Institute, Cary, N.C., USA).

RESULTS

The study population of n=569 patients included n=322 (56,59%) women and n=247 (43,40%) men. The mean average age of patients with alveolar echinococcosis at first presentation was 53,90±17,54 years (median: 56,00 years). The Moran's I test statistic showed a positive spatial autocorrelation for the period 1992-2018 corresponding to a heterogeneous distribution of disease cases in Germany (I=0,4165; Z=10,9591, p=0,001). An increase in age- and sex-specific prevalence could be determined for the entire study period (1992-2018). The overall prevalence in the period 1992-2018 was 0,71 cases per 100,000 population. The determination of the prevalence for the period 1992-2018 resulted in 0,31 cases for men and 0,40 cases for women per 100,000 population. In the period 1992-1996, no AE cases had been registered in 11/16 (68,8%) federal states (Berlin, Brandenburg, Bremen, Hamburg, Mecklenburg-Vorpommern, Rheinland-Pfalz, Saarland, Sachsen, Sachsen-Anhalt, Schleswig-Holstein und Thüringen). The evaluation recently shows an increased occurrence of cases in the federal states of Hessen, Rheinland-Pfalz and Nordrhein-Westfalen.

CONCLUSIONS

The analysis shows a rise in prevalence and an increasing number of cases outside the classic endemic areas of Baden-Württemberg and Bavaria.

The effectiveness of pseudomagic traits in promoting divergence and enhancing local adaptation

"Magic traits," in which the same trait is both under divergent ecological selection and forms the basis of assortative mating, have been sought after due to their supposed unique ability to promote divergence with gene flow. Here, we ask how unique magic traits are, by exploring whether a tightly linked complex of a locus under divergent selection and a locus that acts as a mating cue can mimic a magic trait in its divergence. We find that these "pseudomagic traits" can be very effective in promoting divergence; with tight linkage they are essentially as effective as a magic trait and with loose linkage, and even no linkage, divergence can still be enhanced. Distinguishing between magic and pseudomagic traits in empirical cases may thus not be important when inferring their role in divergence. The ability of divergence in the mating trait to drive divergence in the ecological trait by lowering the effective migration rate, which occurs somewhat even without linkage, is particularly striking; magic traits are typically considered to have the other direction of causality. Our results thus suggest that divergence in a mating trait can at least modestly increase local adaption by allowing more ecological divergence, particularly with tighter linkage.

Genetic Incompatibilities Between Mitochondria and Nuclear Genes: Effect on Gene Flow and Speciation

The process of speciation is, according to the biological species concept, the reduction in gene flow between genetically diverging populations. Most of the previous theoretical studies analyzed the effect of nuclear genetic incompatibilities on gene flow. There is, however, an increasing number of empirical examples suggesting that cytoplasmically inherited genetic elements play an important role in speciation. Here, we present a theoretical analysis of mitochondrial driven speciation, in which genetic incompatibilities occur between mitochondrial haplotypes and nuclear alleles. Four population genetic models with mainland-island structure were analyzed that differ with respect to the type of incompatibility and the underlying genetics. Gene flow reduction was measured on selectively neutral alleles of an unlinked locus and quantified by the effective migration rate. Analytical formulae for the different scenarios were derived using the fitness graph method. For the models with haploid genetics, we found that mito-nuclear incompatibilities (MtNI) are as strong as nuclear-nuclear incompatibilities (NNI) in reducing gene flow at the unlinked locus, but only if males and females migrate in equal number. For models with diploid genetics, we found that MtNI reduce gene flow stronger than NNI when incompatibilities are recessive, but weaker when they are dominant. For both haploid and diploid MtNI, we found that gene flow reduction is stronger if females are the migrating sex, but weaker than NNI when males are the migrating sex. These results encourage further examination on the role of mitochondria on genetic divergence and speciation and point toward specific factors (e.g., migrating sex) that could be the focus of an empirical test.

EWET: Data collection and interface for the genetic analysis of Echinococcus multilocularis based on EmsB microsatellite

Evolution and dispersion history on Earth of organisms can best be studied through biological markers in molecular epidemiological studies. The biological diversity of the cestode Echinococcus multilocularis was investigated in different cladistic approaches. First the morphological aspects were explored in connection with its ecology. More recently, molecular aspects were investigated to better understand the nature of the variations observed among isolates. The study of the tandemly repeated multilocus microsatellite EmsB allowed us to attain a high genetic diversity level where other classic markers have failed. Since 2006, EmsB data have been collected on specimens from various endemic foci of the parasite in Europe (in historic and newly endemic areas), Asia (China, Japan and Kyrgyzstan), and North America (Canada and Alaska). Biological data on the isolates and metadata were also recorded (e.g. host, geographical location, EmsB analysis, citation in the literature). In order to make available the data set of 1,166 isolates from classic and aberrant domestic and wild animal hosts (larval lesions and adult worms) and from human origin, an open web access interface, developed in PHP, and connected to a PostgreSQL database, was developed in the EmsB Website for the Echinococcus Typing (EWET) project. It allows researchers to access data collection, perform genetic analyses online (e.g. defining the genetic distance between their own samples and the samples in the database), consult distribution maps of EmsB profiles, and record and share their new EmsB genotyping data. In order to standardize the EmsB analyses performed in the different laboratories throughout the world, a calibrator was developed. The final aim of this project was to gather and arrange available data to permit to better understand the dispersion and transmission patterns of the parasite among definitive and intermediate hosts, in order to organize control strategies on the ground.

The effect of linkage on establishment and survival of locally beneficial mutations

We study invasion and survival of weakly beneficial mutations arising in linkage to an established migration-selection polymorphism. Our focus is on a continent-island model of migration, with selection at two biallelic loci for adaptation to the island environment. Combining branching and diffusion processes, we provide the theoretical basis for understanding the evolution of islands of divergence, the genetic architecture of locally adaptive traits, and the importance of so-called "divergence hitchhiking" relative to other mechanisms, such as "genomic hitchhiking", chromosomal inversions, or translocations. We derive approximations to the invasion probability and the extinction time of a de novo mutation. Interestingly, the invasion probability is maximized at a nonzero recombination rate if the focal mutation is sufficiently beneficial. If a proportion of migrants carries a beneficial background allele, the mutation is less likely to become established. Linked selection may increase the survival time by several orders of magnitude. By altering the timescale of stochastic loss, it can therefore affect the dynamics at the focal site to an extent that is of evolutionary importance, especially in small populations. We derive an effective migration rate experienced by the weakly beneficial mutation, which accounts for the reduction in gene flow imposed by linked selection. Using the concept of the effective migration rate, we also quantify the long-term effects on neutral variation embedded in a genome with arbitrarily many sites under selection. Patterns of neutral diversity change qualitatively and quantitatively as the position of the neutral locus is moved along the chromosome. This will be useful for population-genomic inference. Our results strengthen the emerging view that physically linked selection is biologically relevant if linkage is tight or if selection at the background locus is strong.

The consequences of dominance and gene flow for local adaptation and differentiation at two linked loci

For a subdivided population the consequences of dominance and gene flow for the maintenance of multilocus polymorphism, local adaptation, and differentiation are investigated. The dispersing population inhabits two demes in which selection acts in opposite direction. Fitness is determined additively by two linked diallelic loci with arbitrary intermediate dominance (no over- or underdominance). For weak as well as strong migration, the equilibrium structure is derived. As a special case, a continuous-time continent–island model (CI model) is analyzed, with one-way migration from the continent to the island. For this CI model, the equilibrium and stability configuration is obtained explicitly for weak migration, for strong migration, for independent loci, and for complete linkage. For independent loci, the possible bifurcation patterns are derived as functions of the migration rate. These patterns depend strongly on the degree of dominance. The effects of dominance, linkage, and migration on the amount of linkage disequilibrium (LD) and the degree of local adaptation are explored. Explicit formulas are obtained for D   (=x1x4−x2x3) and r2 (the squared correlation in allelic state). They demonstrate that dominant island alleles increase D and decrease r2. Local adaptation is elevated by dominance of the locally adaptive alleles. The effective migration rate at a linked neutral locus is calculated. If advantageous alleles are dominant, it is decreased only slightly below the actual migration rate. For a quantitative trait that is determined by two additive loci, the influence of dominance on measures of differentiation is studied. Explicit expressions for QST and two types of FST at equilibrium are deduced and their relation is discussed.

The consequences of gene flow for local adaptation and differentiation: a two-locus two-deme model

We consider a population subdivided into two demes connected by migration in which selection acts in opposite direction. We explore the effects of recombination and migration on the maintenance of multilocus polymorphism, on local adaptation, and on differentiation by employing a deterministic model with genic selection on two linked diallelic loci (i.e., no dominance or epistasis). For the following cases, we characterize explicitly the possible equilibrium configurations: weak, strong, highly asymmetric, and super-symmetric migration, no or weak recombination, and independent or strongly recombining loci. For independent loci (linkage equilibrium) and for completely linked loci, we derive the possible bifurcation patterns as functions of the total migration rate, assuming all other parameters are fixed but arbitrary. For these and other cases, we determine analytically the maximum migration rate below which a stable fully polymorphic equilibrium exists. In this case, differentiation and local adaptation are maintained. Their degree is quantified by a new multilocus version of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_\mathrm{ST}$$\end{document} and by the migration load, respectively. In addition, we investigate the invasion conditions of locally beneficial mutants and show that linkage to a locus that is already in migration-selection balance facilitates invasion. Hence, loci of much smaller effect can invade than predicted by one-locus theory if linkage is sufficiently tight. We study how this minimum amount of linkage admitting invasion depends on the migration pattern. This suggests the emergence of clusters of locally beneficial mutations, which may form ‘genomic islands of divergence’. Finally, the influence of linkage and two-way migration on the effective migration rate at a linked neutral locus is explored. Numerical work complements our analytical results.

The effects of linkage and gene flow on local adaptation: a two-locus continent-island model

Population subdivision and migration are generally considered to be important causes of linkage disequilibrium (LD). We explore the combined effects of recombination and gene flow on the amount of LD, the maintenance of polymorphism, and the degree of local adaptation in a subdivided population by analyzing a diploid, deterministic continent-island model with genic selection on two linked loci (i.e., no dominance or epistasis). For this simple model, we characterize explicitly all possible equilibrium configurations. Simple and intuitive approximations for many quantities of interest are obtained in limiting cases, such as weak migration, weak selection, weak or strong recombination. For instance, we derive explicit expressions for the measures D(=p(AB)-p(A)p(B)) and r(2) (the squared correlation in allelic state) of LD. They depend in qualitatively different ways on the migration rate. Remarkably high values of r(2) are maintained between weakly linked loci, especially if gene flow is low. We determine how the maximum amount of gene flow that admits preservation of the locally adapted haplotype, hence of polymorphism at both loci, depends on recombination rate and selection coefficients. We also investigate the evolution of differentiation by examining the invasion of beneficial mutants of small effect that are linked to an already present, locally adapted allele. Mutants of much smaller effect can invade successfully than predicted by naive single-locus theory provided they are at least weakly linked. Finally, the influence of linkage on the degree of local adaptation, the migration load, and the effective migration rate at a neutral locus is explored. We discuss possible consequences for the evolution of genetic architecture, in particular, for the emergence of clusters of tightly linked, slightly beneficial mutations and the evolution of recombination and chromosome inversions.

Sex-specific incompatibility generates locus-specific rates of introgression between species

Disruption of interactions among ensembles of epistatic loci has been shown to contribute to reproductive isolation among various animal and plant species. Under the Bateson-Dobzhansky-Muller model, such interspecific incompatibility arises as a by-product of genetic divergence in each species, and the Orr-Turelli model indicates that the number of loci involved in incompatible interactions may "snowball" over time. We address the combined effect of multiple incompatibility loci on the rate of introgression at neutral marker loci across the genome. Our analysis extends previous work by accommodating sex specificity: differences between the sexes in the expression of incompatibility, in rates of crossing over between neutral markers and incompatibility loci, and in transmission of markers or incompatibility factors. We show that the evolutionary process at neutral markers in a genome subject to incompatibility selection is well approximated by a purely neutral process with migration rates appropriately scaled to reflect the influence of selection targeted to incompatibility factors. We confirm that in the absence of sex specificity and functional epistasis among incompatibility factors, the barrier to introgression induced by multiple incompatibility factors corresponds to the product of the barriers induced by the factors individually. A new finding is that barriers to introgression due to sex-specific incompatibility depart in general from multiplicativity. Our partitioning of variation in relative reproductive rate suggests that such departures derive from associations between sex and incompatibility and between sex and neutral markers. Concordant sex-specific incompatibility (for example, greater impairment of male hybrids or longer map lengths in females) induces lower barriers (higher rates of introgression) than expected under multiplicativity, and discordant sex-specific incompatibility induces higher barriers.

Effects of polymorphism for locally adapted genes on rates of neutral introgression in structured populations

Adaptation to local conditions within demes balanced by migration can maintain polymorphisms for variants that reduce fitness in certain ecological contexts. Here, we address the effects of such polymorphisms on the rate of introgression of neutral marker genes, possibly genetically linked to targets of selection. Barriers to neutral gene flow are expected to increase with linkage to targets of local selection and with differences between demes in the frequencies of locally adapted alleles. This expectation is borne out under purifying and disruptive selection, regimes that promote monomorphism within demes. In contrast, overdominance within demes induces minimal barriers to neutral introgression even in the face of very large differences between demes in the frequencies of locally adapted alleles. Further, segregation distortion, a phenomenon observed in a number of interspecific hybrids, can in fact promote transmission by migrants to future generations at rates exceeding those of residents.

Effect of parasitic sex-ratio distorters on host gene frequencies in a mainland-island context

It was previously argued that infection by parasitic sex-ratio distorters can enhance both random genetic drift and genetic influx from outside the population. However, these two enhancement effects have been studied independently. Here, we study the equilibrium frequencies of alleles (neutral and selected) in a mainland-island scenario where both genetic drift and genetic influx are enhanced due to infection by a cytoplasmic feminizing element. Interestingly, our model reveals that at neutral loci, the two effects almost exactly cancel each other out, such that infection has only a very minor effect on the equilibrium frequency distributions of alleles. At selected loci, in contrast, the two effects are unbalanced and infection has conspicuous effects. Despite the cryptic effects of infection at neutral loci, we demonstrate that temporally spaced data can be used to evaluate the effect of infection on genetic drift and that on gene flow separately.

The concept of effective recombination rate and its application in speciation theory

The goal of this study is to develop a unifying theoretical framework to quantify the strength of reproductive isolation. We propose the use of the "effective recombination rate," which measures how fast associations of genes are broken by interlocus recombination. Applying the well-established theory of the effective migration rate, we derive two techniques to investigate the effective recombination rate in models of speciation: the weak migration approximation for parapatric scenarios and the weak recombination approximation for sympatric scenarios. We illustrate the use of these two methods by two examples each: (1) single-locus genetic incompatibility and (2) two-locus genetic incompatibility for the first method, and (3) assortative mating and (4) assortative mating combined with disruptive selection for the second method. An advantage of the effective recombination rate over previous approaches is that it integrates gene flow in both directions into a single index measuring the strength of isolation. This enables straightforward comparisons of speciation scenarios with the same or different geographic histories. The method also allows us to evaluate the relative contributions of F2 hybrid deficiency or linkage between multiple barriers in reproductive isolation.

Genetic diversity of the cestode Echinococcus multilocularis in red foxes at a continental scale in Europe

BACKGROUND

Alveolar echinococcosis (AE) is a severe helminth disease affecting humans, which is caused by the fox tapeworm Echinococcus multilocularis. AE represents a serious public health issue in larger regions of China, Siberia, and other regions in Asia. In Europe, a significant increase in prevalence since the 1990s is not only affecting the historically documented endemic area north of the Alps but more recently also neighbouring regions previously not known to be endemic. The genetic diversity of the parasite population and respective distribution in Europe have now been investigated in view of generating a fine-tuned map of parasite variants occurring in Europe. This approach may serve as a model to study the parasite at a worldwide level.

METHODOLOGY/PRINCIPAL FINDINGS

The genetic diversity of E. multilocularis was assessed based upon the tandemly repeated microsatellite marker EmsB in association with matching fox host geographical positions. Our study demonstrated a higher genetic diversity in the endemic areas north of the Alps when compared to other areas.

CONCLUSIONS/SIGNIFICANCE

The study of the spatial distribution of E. multilocularis in Europe, based on 32 genetic clusters, suggests that Europe can be considered as a unique global focus of E. multilocularis, which can be schematically drawn as a central core located in Switzerland and Jura Swabe flanked by neighbouring regions where the parasite exhibits a lower genetic diversity. The transmission of the parasite into peripheral regions is governed by a "mainland-island" system. Moreover, the presence of similar genetic profiles in both zones indicated a founder event.