Revisiting the historical scenario of a disease dissemination using genetic data and Approximate Bayesian Computation methodology: The case of Pseudocercospora fijiensis invasion in Africa.

The reconstruction of geographic and demographic scenarios of dissemination for invasive pathogens of crops is a key step toward improving the management of emerging infectious diseases. Nowadays, the reconstruction of biological invasions typically uses the information of both genetic and historical information to test for different hypotheses of colonization. The Approximate Bayesian Computation framework and its recent Random Forest development (ABC-RF) have been successfully used in evolutionary biology to decipher multiple histories of biological invasions. Yet, for some organisms, typically plant pathogens, historical data may not be reliable notably because of the difficulty to identify the organism and the delay between the introduction and the first mention. We investigated the history of the invasion of Africa by the fungal pathogen of banana Pseudocercospora fijiensis, by testing the historical hypothesis against other plausible hypotheses. We analyzed the genetic structure of eight populations from six eastern and western African countries, using 20 microsatellite markers and tested competing scenarios of population foundation using the ABC-RF methodology. We do find evidence for an invasion front consistent with the historical hypothesis, but also for the existence of another front never mentioned in historical records. We question the historical introduction point of the disease on the continent. Crucially, our results illustrate that even if ABC-RF inferences may sometimes fail to infer a single, well-supported scenario of invasion, they can be helpful in rejecting unlikely scenarios, which can prove much useful to shed light on disease dissemination routes.

Highly pathogenic avian influenza A/H5N1 Clade 2.3.2.1c virus in poultry in Cameroon, 2016-2017.

In May 2016, highly pathogenic avian influenza virus of the subtype A/H5N1 was detected in Cameroon in an industrial poultry farm at Mvog-Betsi, Yaoundé (Centre region), with a recorded sudden increase of deaths among chickens, and an overall mortality rate of 75%. The virus spread further and caused new outbreaks in some parts of the country. In total, 21 outbreaks were confirmed from May 2016 to March 2017 (six in the Centre, six in the West, eight in the South and one in the Adamaoua regions). This resulted in an estimated total loss of 138,252 birds (44,451 deaths due to infection and 93,801 stamped out). Only domestic birds (chickens, ducks and geese) were affected in farms as well as in poultry markets. The outbreaks occurred in three waves, the first from May to June 2016, the second in September 2016 and the last wave in March 2017. The topology of the phylogeny based on the haemagglutinin gene segment indicated that the causative H5N1 viruses fall within the genetic clade 2.3.2.1c, within the same group as the A/H5N1 viruses collected in Niger in 2015 and 2016. More importantly, the gene constellation of four representative viruses showed evidence of H5N1/H9N2 intra-clade reassortment. Additional epidemiological and genetic data from affected countries in West Africa are needed to better trace the origin, spread and evolution of A/H5N1 in Cameroon. RESEARCH HIGHLIGHTS HPAI A/H5N1 was detected in May 2016 in domestic chickens in Yaoundé-Cameroon. Twenty-one outbreaks in total were confirmed from May 2016 to March 2017. The causative H5N1 viruses fall within the genetic clade 2.3.2.1c. The viral gene constellation showed evidence of H5N1/H9N2 intra-clade reassortment.

Deciphering the demographic history of allochronic differentiation in the pine processionary moth Thaumetopoea pityocampa.

Understanding the processes of adaptive divergence, which may ultimately lead to speciation, is a major question in evolutionary biology. Allochronic differentiation refers to a particular situation where gene flow is primarily impeded by temporal isolation between early and late reproducers. This process has been suggested to occur in a large array of organisms, even though it is still overlooked in the literature. We here focused on a well-documented case of incipient allochronic speciation in the winter pine processionary moth Thaumetopoea pityocampa. This species typically reproduces in summer and larval development occurs throughout autumn and winter. A unique, phenologically shifted population (SP) was discovered in 1997 in Portugal. It was proved to be strongly differentiated from the sympatric "winter population" (WP), but its evolutionary history could only now be explored. We took advantage of the recent assembly of a draft genome and of the development of pan-genomic RAD-seq markers to decipher the demographic history of the differentiating populations and develop genome scans of adaptive differentiation. We showed that the SP diverged relatively recently, that is, few hundred years ago, and went through two successive bottlenecks followed by population size expansions, while the sympatric WP is currently experiencing a population decline. We identified outlier SNPs that were mapped onto the genome, but none were associated with the phenological shift or with subsequent adaptations. The strong genetic drift that occurred along the SP lineage certainly challenged our capacity to reveal functionally important loci.

Genome scans reveal candidate regions involved in the adaptation to host plant in the pea aphid complex.

A major goal in evolutionary biology is to uncover the genetic basis of adaptation. Divergent selection exerted on ecological traits may result in adaptive population differentiation and reproductive isolation and affect differentially the level of genetic divergence along the genome. Genome-wide scan of large sets of individuals from multiple populations is a powerful approach to identify loci or genomic regions under ecologically divergent selection. Here, we focused on the pea aphid, a species complex of divergent host races, to explore the organization of the genomic divergence associated with host plant adaptation and ecological speciation. We analysed 390 microsatellite markers located at variable distances from predicted genes in replicate samples of sympatric populations of the pea aphid collected on alfalfa, red clover and pea, which correspond to three common host-adapted races reported in this species complex. Using a method that accounts for the hierarchical structure of our data set, we found a set of 11 outlier loci that show higher genetic differentiation between host races than expected under the null hypothesis of neutral evolution. Two of the outliers are close to olfactory receptor genes and three other nearby genes encoding salivary proteins. The remaining outliers are located in regions with genes of unknown functions, or which functions are unlikely to be involved in interactions with the host plant. This study reveals genetic signatures of divergent selection across the genome and provides an inventory of candidate genes responsible for plant specialization in the pea aphid, thereby setting the stage for future functional studies.

Population differentiation of sessile oak at the altitudinal front of migration in the French Pyrenees.

To assess the effects of altitude on the level and structure of genetic diversity, a genetic survey was conducted in 12 populations of sessile oak (Quercus petraea) located between 130 and 1660 m in two parallel valleys on the northern side of the Pyrenees Mountains. Genetic diversity was monitored at 16 nuclear microsatellite loci and 5 chloroplast DNA (cpDNA) markers. The cpDNA survey suggested that extant populations in both valleys shared the same source populations from the plain. There was no visible trend of nuclear genetic diversity along altitude, even if indirect estimates of effective population sizes revealed a consistent reduction at higher altitudes. Population differentiation, although low, was mostly present among populations of the same valleys and reached similar levels than differentiation across the range of distribution of sessile oak. Contribution to the overall differentiation in the valleys was mostly due to the genetic divergence of the highest populations and the altitudinal variation of allelic frequencies at a few loci. Bayesian inference of migration between groups of populations showed that gene flow is preferentially unidirectional from lower altitudes in one valley to other groups of populations. Finally, we found evidence of clonal reproduction in high altitude populations. The introgression of Quercus robur and Quercus pubescens was also more frequent at the altitudinal margin suggesting that this mechanism may have contributed to the present migration and adaptation of Q. petraea and may also facilitate its future upslope shift in the context of climate change.

Quantifying male-biased dispersal among social groups in the collared peccary (Pecari tajacu) using analyses based on mtDNA variation.

Recent advances in the statistical analysis of microsatellite data permit calculation of sex-specific dispersal rates through sex- and age-specific comparisons of genetic variation. This approach, developed for the analysis of data derived from co-dominant autosomal markers, should be applicable to a sex-specific marker such as mitochondrial DNA. To test this premise, we amplified a 449 bp control region DNA sequence from the mitochondrial genome of the collared peccary (Pecari tajacu), and estimated intra-class correlations among herds sampled from three Texas populations. Analyses on data partitioned by breeding group showed a clear signal of male-biased dispersal; sex-specific fixation indices associated with genetic variation among social groups within populations yielded values for females (F(GP)=0.91), which were significantly larger than values for males (F(GP)=0.24; P=0.0015). The same general pattern emerged when the analyses were conducted on age classes (albeit nonsignificantly), as well as categories of individuals that were predicted a posteriori to be dispersers (adult males) and philopatric (adult females and all immatures). By extending a previously published methodology based on biparentally inherited markers to matrilineally inherited haploid data, we calculated sex-specific rates of contemporary dispersal among social groups within populations (m(male symbol)=0.37). These results support the idea that mitochondrial DNA haplotype frequency data can be used to estimate sex-specific instantaneous dispersal rates in a social species.

DetSel 1.0: a computer program to detect markers responding to selection.

Estimating population parameters from polymorphism frequency data requires neutral genetic markers. Any departure from neutrality may invalidate the inferences drawn from such analyses. We recently discussed the possibility of identifying markers that show deviation from neutral expectations in pairwise comparisons of diverging populations. We are now releasing a user-friendly software package that implements all the necessary steps to identify the signature of selection among molecular markers in a set of polymorphism data. This software can be downloaded free of charge at http://www.univ-montp2.fr/~genetix/detsel/detsel.html.

Interpretation of variation across marker loci as evidence of selection.

Population structure and history have similar effects on the genetic diversity at all neutral loci. However, some marker loci may also have been strongly influenced by natural selection. Selection shapes genetic diversity in a locus-specific manner. If we could identify those loci that have responded to selection during the divergence of populations, then we may obtain better estimates of the parameters of population history by excluding these loci. Previous attempts were made to identify outlier loci from the distribution of sample statistics under neutral models of population structure and history. Unfortunately these methods depend on assumptions about population structure and history that usually cannot be verified. In this article, we define new population-specific parameters of population divergence and construct sample statistics that are estimators of these parameters. We then use the joint distribution of these estimators to identify outlier loci that may be subject to selection. We found that outlier loci are easier to recognize when this joint distribution is conditioned on the total number of allelic states represented in the pooled sample at each locus. This is so because the conditional distribution is less sensitive to the values of nuisance parameters.

Two-locus identity probabilities and identity disequilibrium in a partially selfing subdivided population.

Measures of association of genes at different loci (linkage disequilibrium) are widely used to determine whether the structure of natural populations is clonal or not, to map genes from population data, or to test for the homogeneity of response of molecular markers to background selection, for example. However, the usual definitions of parameters for gametic associations may not be suitable for all these purposes. In this paper, we derive the recursion equations for one- and two-locus identity probabilities in an infinite island model. We study the role of drift, gene flow, partial selfing and mutation model on the expected association of genes across loci. We define the 'within-subpopulation identity disequilibrium' as the difference between the joint two-locus probability of identity in state and the expected product of one-locus identity probabilities. We evaluate this parameter as a function of recombination rate, effective size, gene flow and selfing rate. Within-subpopulation identity disequilibrium attains maximum values for intermediate immigration rates, whatever the selfing rate. Moreover, identity disequilibrium may be very small, even for high selfing rates. We discuss the implications of these findings for the analysis of data from natural populations.

Estimation of effective population size and migration rate from one- and two-locus identity measures.

Standard methods for inferring demographic parameters from genetic data are based mainly on one-locus theory. However, the association of genes at different loci (e.g., two-locus identity disequilibrium) may also contain some information about demographic parameters of populations. In this article, we define one- and two-locus parameters of population structure as functions of one- and two-locus probabilities for the identity in state of genes. Since these parameters are known functions of demographic parameters in an infinite island model, we develop moment-based estimators of effective population size and immigration rate from one- and two-locus parameters. We evaluate this method through simulation. Although variance and bias may be quite large, increasing the number of loci on which the estimates are derived improves the method. We simulate an infinite allele model and a K allele model of mutation. Bias and variance are smaller with increasing numbers of alleles per locus. This is, to our knowledge, the first attempt of a joint estimation of local effective population size and immigration rate.

Relationships between transposable elements based upon the integrase-transposase domains: is there a common ancestor?

The integrase domain of RNA-mediated elements (class I) and the transposase domain of DNA-mediated transposable elements (class II) were compared. A number of elements contain the DDE signature, which plays an important role in their integration. The possible relationships between mariner-Tc1 and IS elements, retrotransposons, and retroviruses were analyzed from an alignment of this region. The mariner-Tc1 superfamily, and LTR retrotransposons and retroviruses were found to be monophyletic groups. However, the IS elements of bacteria were found in several groups. These results were used to propose an evolutionary history that suggests a common ancestor for some integrases and transposases.