Search

Disentangling the effects of selection and drift is a long-standing problem in population genetics. Simulations show that pervasive selection may bias the inference of demography. Ideally, models for the inference of demography and selection should account for the interaction between these two forces. With simulation-based likelihood-free methods such as Approximate Bayesian Computation (ABC), demography and selection parameters can be jointly estimated. We propose to use the ABC-Random Forests framework to jointly infer demographic and selection parameters from temporal population genomic data (e.g. experimental evolution, monitored populations, ancient DNA). Our framework allowed the separation of demography (census size, N) from the genetic drift (effective population size, Ne) and the estimation of genome-wide parameters of selection. Selection parameters informed us about the adaptive potential of a population (the scaled mutation rate of beneficial mutations, $\theta_{\mathrm{b}}$), the realized adaptation (the number of mutations under strong selection), and population fitness (genetic load). We applied this approach to a dataset of feral populations of honey bees (Apis mellifera) collected in California, and we estimated parameters consistent with the biology and the recent history of this species.

The breeding of resistant hosts based on the gene-for-gene interaction is crucial to address epidemics of plant pathogens in agroecosystems. Resistant host deployment strategies are developed and studied worldwide to decrease the probability of resistance breakdown and increase the resistance durability in various pathosystems. A major component of deployment strategies is the proportion of resistant hosts in the landscape. However, the impact of this proportion on resistance durability remains unclear for diploid pathogens with complex life cycles. In this study, we modelled pathogen population dynamics and genetic evolution at the virulence locus to assess the impact of the ploidy (haploid or diploid) and the pathogen’s life cycle (with or without host alternation) on resistance durability. Ploidy has a strong impact on evolutionary trajectories, with much greater stochasticity and delayed times of resistance breakdown for diploids. This result emphasises the importance of genetic drift in this system: as the virulent allele is recessive, positive selection on resistant hosts only applies to homozygous (virulent) individuals, which may lead to population collapse at low frequencies of the virulent allele. We also observed differences in the effect of host deployment depending on the pathogen’s life cycle. With host alternation, the probability that the pathogen population collapses strongly increases with the proportion of resistant hosts in the landscape. Therefore, resistance breakdown events occurring at high proportions of resistant hosts frequently amount to evolutionary rescue. Last, life cycles correspond to two selection regimes: without host alternation (soft selection) the resistance breakdown is mainly driven by the migration rate. Conversely, host alternation (hard selection) resembles an all-or-nothing game, with stochastic trajectories caused by the recurrent allele redistributions on the alternate host.

The deployment of plant varieties carrying resistance genes (R) exerts strong selection pressure on pathogen populations. Rapidly evolving avirulence genes (Avr) allow pathogens to escape R‐mediated plant immunity through a variety of mechanisms, leading to virulence. The poplar rust fungus Melampsora larici‐populina is a damaging pathogen of poplars in Europe. It underwent a major adaptive event in 1994, with the breakdown of the poplar RMlp7 resistance gene. Population genomics studies identified a locus in the genome of M. larici‐populina that probably corresponds to the candidate avirulence gene AvrMlp7. Here, to further characterize this effector, we used a population genetics approach on a comprehensive set of 281 individuals recovered throughout a 28‐year period encompassing the resistance breakdown event. Using two dedicated molecular tools, genotyping at the candidate locus highlighted two different alterations of a predominant allele found mainly before the resistance breakdown: a nonsynonymous mutation and a complete deletion of this locus. This results in six diploid genotypes: three genotypes related to the avirulent phenotype and three related to the virulent phenotype. The temporal survey of the candidate locus revealed that both alterations were found in association during the resistance breakdown event. They pre‐existed before the breakdown in a heterozygous state with the predominant allele cited above. Altogether, these results suggest that the association of both alterations at the candidate locus AvrMlp7 drove the poplar rust adaptation to RMlp7‐mediated immunity. This study demonstrates for the first time a case of adaptation from standing genetic variation in rust fungi during a qualitative resistance breakdown. [ABSTRACT FROM AUTHOR]

Rapid and repeatable polymorphism analyses have become a necessity with the current amount of genomic data that can be collected in many organisms. Traditionally, such analyses are conducted using a variety of tools in combination, often requiring numerous format translation and manipulation. Here, we present a massively updated version of our previous software package egglib, intended to alleviate such costly and error‐prone tinkering with the data. egglib has been streamlined into a python package and thoroughly updated and optimized to accommodate modern‐day sized dataset. We show the main characteristics of the package making it a tool of choice to perform population genetics analyses. Once the data are imported (whatever their encoding), they can be filtered, edited, analysed and compared to coalescent simulations very easily and efficiently. Furthermore, the list of diversity and polymorphism statistics that can now be calculated has been greatly expanded. The software and its full documentation are available at https://egglib.org/. [ABSTRACT FROM AUTHOR]

Abstract Background With the considerable growth of available nucleotide sequence data over the last decade, integrated and flexible analytical tools have become a necessity. In particular, in the field of population genetics, there is a strong need for automated and reliable procedures to conduct repeatable and rapid polymorphism analyses, coalescent simulations, data manipulation and estimation of demographic parameters under a variety of scenarios. Results In this context, we present EggLib (Evolutionary Genetics and Genomics Library), a flexible and powerful C++/Python software package providing efficient and easy to use computational tools for sequence data management and extensive population genetic analyses on nucleotide sequence data. EggLib is a multifaceted project involving several integrated modules: an underlying computationally efficient C++ library (which can be used independently in pure C++ applications); two C++ programs; a Python package providing, among other features, a high level Python interface to the C++ library; and the egglib script which provides direct access to pre-programmed Python applications. Conclusions EggLib has been designed aiming to be both efficient and easy to use. A wide array of methods are implemented, including file format conversion, sequence alignment edition, coalescent simulations, neutrality tests and estimation of demographic parameters by Approximate Bayesian Computation (ABC). Classes implementing different demographic scenarios for ABC analyses can easily be developed by the user and included to the package. EggLib source code is distributed freely under the GNU General Public License (GPL) from its website http://egglib.sourceforge.net/ where a full documentation and a manual can also be found and downloaded.

Crop varieties carrying qualitative resistance to targeted pathogens lead to strong selection pressure on parasites, often resulting in resistance breakdown. It is well known that qualitative resistance breakdowns modify pathogen population structure but few studies have analyzed the consequences on their quantitative aggressiveness‐related traits. The aim of this study was to characterize the evolution of these traits following a resistance breakdown in the poplar rust fungus, Melampsora larici‐populina. We based our experiment on three temporal populations sampled just before the breakdown event, immediately after and four years later. First, we quantified phenotypic differences among populations for a set of aggressiveness traits on a universally susceptible cultivar (infection efficiency, latent period, lesion size, mycelium quantity, and sporulation rate) and one morphological trait (mean spore volume). Then, we estimated heritability to establish which traits could be subjected to adaptive evolution and tested for evidence of selection. Our results revealed significant changes in the morphological trait but no variation in aggressiveness traits. By contrast, recent works have demonstrated that quantitative resistance (initially assumed more durable) could be eroded and lead to increased aggressiveness. Hence, this study is one example suggesting that the use of qualitative resistance may be revealed to be less detrimental to long‐term sustainable crop production. [ABSTRACT FROM AUTHOR]

The genetic consequences of range expansions have generally been investigated at wide geographical and temporal scales, long after the colonization event. A unique ecological system enabled us to both monitor the colonization dynamics and decipher the genetic footprints of expansion over a very short time period. Each year an epidemic of the poplar rust ( Melampsora larici-populina) expands clonally and linearly along the Durance River, in the Alps. The colonization dynamics observed in 2004 showed two phases with different genetic outcomes. Upstream, fast colonization maintained high genetic diversity. Downstream, the colonization wave progressively faltered, diversity eroded, and differentiation increased, as expected under recurrent founder events. In line with the high dispersal abilities of rust pathogens, we provide evidence for leapfrog dispersal of clones. Our results thus emphasize the importance of colonization dynamics in shaping spatial genetic structure in the face of high gene flow. [ABSTRACT FROM AUTHOR]

The Périgord black truffle ( Tuber melanosporum Vittad.), considered a gastronomic delicacy worldwide, is an ectomycorrhizal filamentous fungus that is ecologically important in Mediterranean French, Italian and Spanish woodlands. In this study, we developed a novel resource of single nucleotide polymorphisms ( SNPs) for T. melanosporum using Illumina high-throughput resequencing. The genome from six T. melanosporum geographical accessions was sequenced to a depth of approximately 20×. These geographical accessions were selected from different populations within the northern and southern regions of the geographical species distribution. Approximately 80% of the reads for each of the six resequenced geographical accessions mapped against the reference T. melanosporum genome assembly, estimating the core genome size of this organism to be approximately 110 Mbp. A total of 442 326 SNPs corresponding to 3540 SNPs/Mbps were identified as being included in all seven genomes. The SNPs occurred more frequently in repeated sequences (85%), although 4501 SNPs were also identified in the coding regions of 2587 genes. Using the ratio of nonsynonymous mutations per nonsynonymous site ( pN) to synonymous mutations per synonymous site ( pS) and Tajima's D index scanning the whole genome, we were able to identify genomic regions and genes potentially subjected to positive or purifying selection. The SNPs identified represent a valuable resource for future population genetics and genomics studies. [ABSTRACT FROM AUTHOR]

The poplar rust fungus Melampsora larici-populina causes significant yield reduction and severe economic losses in commercial poplar plantations. After several decades of breeding for qualitative resistance and subsequent breakdown of the released resistance genes, breeders now focus on quantitative resistance, perceived to be more durable. But quantitative resistance also can be challenged by an increase of aggressiveness in the pathogen. Thus, it is of primary importance to better understand the genetic architecture of aggressiveness traits. To this aim, our goal is to build a genetic linkage map for M. larici-populina in order to map quantitative trait loci related to aggressiveness. First, a large progeny of M. larici-populina was generated through selfing of the reference strain 98AG31 (which genome sequence is available) on larch plants, the alternate host of the poplar rust fungus. The progeny's meiotic origin was validated through a segregation analysis of 115 offspring with 14 polymorphic microsatellite markers, of which 12 segregated in the expected 1:2:1 Mendelian ratio. A microsatellite-based linkage disequilibrium analysis allowed us to identify one potential linkage group comprising two scaffolds. The whole genome of a subset of 47 offspring was resequenced using the Illumina HiSeq 2000 technology at a mean sequencing depth of 6X. The reads were mapped onto the reference genome of the parental strain and 144,566 SNPs were identified across the genome. Analysis of distribution and polymorphism of the SNPs along the genome led to the identification of 2580 recombination blocks. A second linkage disequilibrium analysis, using the recombination blocks as markers, allowed us to group 81 scaffolds into 23 potential linkage groups. These preliminary results showed that a high-density linkage map could be constructed by using high-quality SNPs based on low-coverage resequencing of a larger number of M. larici-populina offspring. [ABSTRACT FROM AUTHOR]

Approximate Bayesian Computation (ABC) has become a popular technique in evolutionary genetics for elucidating population structure and history due to its flexibility. The statistical inference framework has benefited from significant progress in recent years. In population genetics, however, its outcome depends heavily on the amount of information in the dataset, whether that be the level of genetic variation or the number of samples and loci. Here we look at the power to reject a simple constant population size coalescent model in favor of a bottleneck model in datasets of varying quality. Not only is this power dependent on the number of samples and loci, but it also depends strongly on the level of nucleotide diversity in the observed dataset. Whilst overall model choice in an ABC setting is fairly powerful and quite conservative with regard to false positives, detecting weaker bottlenecks is problematic in smaller or less genetically diverse datasets and limits the inferences possible in non-model organism where the amount of information regarding the two models is often limited. Our results show it is important to consider these limitations when performing an ABC analysis and that studies should perform simulations based on the size and nature of the dataset in order to fully assess the power of the study. [ABSTRACT FROM AUTHOR]

The ability of plants to track seasonal changes is largely dependent on genes assigned to the photoperiod pathway, and variation in those genes is thereby important for adaptation to local day length conditions. Extensive physiological data in several temperate conifer species suggest that populations are adapted to local light conditions, but data on the genes underlying this adaptation are more limited. Here we present nucleotide diversity data from 19 genes putatively involved in photoperiodic response in Norway spruce (Picea abies). Based on similarity to model plants the genes were grouped into three categories according to their presumed position in the photoperiod pathway: photoreceptors, circadian clock genes, and downstream targets. An HKA (Hudson, Kreitman and Aquade) test showed a significant excess of diversity at photoreceptor genes, but no departure from neutrality at circadian genes and downstream targets. Departures from neutrality were also tested with Tajima's D and Fay and Wu's H statistics under three demographic scenarios: the standard neutral model, a population expansion model, and a more complex population split model. Only one gene, the circadian clock gene PaPRR3 with a highly positive Tajima's D value, deviates significantly from all tested demographic scenarios. As the PaPRR3 gene harbours multiple non-synonymous variants it appears as an excellent candidate gene for control of photoperiod response in Norway spruce. [ABSTRACT FROM AUTHOR]

The symbiosis between legumes and nitrogen-fixing rhizobia co-opted pre-existing endomycorrhizal features. In particular, both symbionts release lipo-chitooligosaccharides ( LCOs) that are recognized by Lys M-type receptor kinases. We investigated the evolutionary history of rhizobial LCO receptor genes Mt LYK3- Lj NFR1 to gain insight into the evolutionary origin of the rhizobial symbiosis., We performed a phylogenetic analysis integrating gene copies from nonlegumes and legumes, including the non-nodulating, phylogenetically basal legume Cercis chinensis. Signatures of differentiation between copies were investigated through patterns of molecular evolution., We show that two rounds of duplication preceded the evolution of the rhizobial symbiosis in legumes. Molecular evolution patterns indicate that the resulting three paralogous gene copies experienced different selective constraints. In particular, one copy maintained the ancestral function, and another specialized into perception of rhizobial LCOs. It has been suggested that legume LCO receptors evolved from a putative ancestral defense-related chitin receptor through the acquisition of two kinase motifs. However, the phylogenetic analysis shows that these domains are actually ancestral, suggesting that this scenario is unlikely., Our study underlines the evolutionary significance of gene duplication and subsequent neofunctionalization in Mt LYK3- Lj NFR1 genes. We hypothesize that their ancestor was more likely a mycorrhizal LCO receptor, than a defense-related receptor kinase. [ABSTRACT FROM AUTHOR]

Transcription factors belonging to the CCAAT-box binding factor family (also known as the Nuclear Factor Y) are present in all higher eukaryotes. Studies in plants have revealed that each subunit of this heterotrimeric transcription factor is encoded by a gene belonging to a multigene family allowing a considerable modularity. In this review, we focus on recent findings concerning the expression patterns and potential functions of different members of these NF-Y protein families using a phylogenetic approach. During the course of evolution plant CCAAT-box binding factors seem to have diversified into at least two main groups. The first group has more general expression patterns and/or functions whereas the second group has acquired more specific expression patterns and/or functions and could play key roles in specific pathways. [Copyright &y& Elsevier]

While genomic erosion is common among intracellular symbionts, patterns of genome evolution in heritable extracellular endosymbionts remain elusive. We study vertically transmitted extracellular endosymbionts (Verminephrobacter, Betaproteobacteria) that form a beneficial, species-specific, and evolutionarily old (60–130 Myr) association with earthworms. We assembled a draft genome of Verminephrobacter aporrectodeae and compared it with the genomes of Verminephrobacter eiseniae and two nonsymbiotic close relatives (Acidovorax). Similar to V. eiseniae, the V. aporrectodeae genome was not markedly reduced in size and showed no A–T bias. We characterized the strength of purifying selection (ω = dN/dS) and codon usage bias in 876 orthologous genes. Symbiont genomes exhibited strong purifying selection (ω = 0.09 ± 0.07), although transition to symbiosis entailed relaxation of purifying selection as evidenced by 50% higher ω values and less codon usage bias in symbiont compared with reference genomes. Relaxation was not evenly distributed among functional gene categories but was overrepresented in genes involved in signal transduction and cell envelope biogenesis. The same gene categories also harbored instances of positive selection in the Verminephrobacter clade. In total, positive selection was detected in 89 genes, including also genes involved in DNA metabolism, tRNA modification, and TonB-dependent iron uptake, potentially highlighting functions important in symbiosis. Our results suggest that the transition to symbiosis was accompanied by molecular adaptation, while purifying selection was only moderately relaxed, despite the evolutionary age and stability of the host association. We hypothesize that biparental transmission of symbionts and rare genetic mixing during transmission can prevent genome erosion in heritable symbionts. [ABSTRACT FROM AUTHOR]

Legumes host their Rhizobium spp. symbiont in novel root organs called nodules. Nodules originate from differentiated root cortical cells that dedifferentiate and subsequently form nodule primordia, a process controlled by cytokinin. A whole-genome duplication has occurred at the root of the legume Papilionoideae subfamily. We hypothesize that gene pairs originating from this duplication event and are conserved in distinct Papilionoideae lineages have evolved symbiotic functions. A phylogenetic strategy was applied to search for such gene pairs to identify novel regulators of nodulation, using the cytokinin phosphorelay pathway as a test case. In this way, two paralogous type-A cytokinin response regulators were identified that are involved in root nodule symbiosis. Response Regulator9 (MtRR9) and MtRR11 in medicago (Medicago truncatula) and an ortholog in lotus (Lotus japonicus) are rapidly induced upon Rhizobium spp. Nod factor signaling. Constitutive expression of MtRR9 results in arrested primordia that have emerged from cortical, endodermal, and pericycle cells. In legumes, lateral root primordia are not exclusively formed from pericycle cells but also require the involvement of the root cortical cell layer. Therefore, the MtRR9-induced foci of cell divisions show a strong resemblance to lateral root primordia, suggesting an ancestral function of MtRR9 in this process. Together, these findings provide a proof of principle for the applied phylogenetic strategy to identify genes with a symbiotic function in legumes. [ABSTRACT FROM AUTHOR]

Climate changes will have an impact on food production and will require costly adaptive responses. Adapting to a changing environment will be particularly challenging in sub-Saharan Africa where climate change is expected to have a major impact. However, one important phenomenon that is often overlooked and is poorly documented is the ability of agrosystems to rapidly adapt to environmental variations. Such an adaptation could proceed by the adoption of new varieties or by the adaptation of varieties to a changing environment. In this study, we analyzed these two processes in one of the driest agro-ecosystems in Africa, the Sahel. We performed a detailed study in Niger where pearl millet is the main crop and covers 65% of the cultivated area. To assess how the agro-system is responding to recent recurrent drought, we analyzed samples of pearl millet landraces collected in the same villages in 1976 and 2003 throughout the entire cultivated area of Niger. We studied phenological and morphological differences in the 1976 and 2003 collections by comparing them over three cropping seasons in a common garden experiment. We found no major changes in the main cultivated varieties or in their genetic diversity. However, we observed a significant shift in adaptive traits. Compared to the 1976 samples, samples collected in 2003 displayed a shorter lifecycle, and a reduction in plant and spike size. We also found that an early flowering allele at the PHYC locus increased in frequency between 1976 and 2003. The increase exceeded the effect of drift and sampling, suggesting a direct effect of selection for earliness on this gene. We conclude that recurrent drought can lead to selection for earlier flowering in a major Sahelian crop. Surprisingly, these results suggest that diffusion of crop varieties is not the main driver of short term adaptation to climatic variation. [ABSTRACT FROM AUTHOR]

Symbiotic nitrogen-fixing rhizobia are able to trigger root deformation in their Fabaceae host plants, allowing their intracellular accommodation. They do so by delivering molecules called Nod factors. We analyzed the patterns of nucleotide polymorphism of five genes controlling early Nod factor perception and signaling in the Fabaceae Medicago truncatula to understand the selective forces shaping the evolution of these genes. We used 30 M. truncatula genotypes sampled in a genetically homogeneous region of the species distribution range. We first sequenced 24 independent loci and detected a genomewide departure from the hypothesis of neutrality and demographic equilibrium that suggests a population expansion. These data were used to estimate parameters of a simple demographic model incorporating population expansion. The selective neutrality of genes controlling Nod factor perception was then examined using a combination of two complementary neutrality tests, Tajima's D and Fay and Wu's standardized H. The joint distribution of D and Hexpected under neutrality was obtained under the fitted population expansion model. Only the gene DMI1, which is expected to regulate the downstream signal, shows a pattern consistent with a putative selective event. In contrast, the receptor-encoding genes NFP and NORK show no significant signatures of selection. Among the genes that we analyzed, only DMI1 should be viewed as a candidate for adaptation in the recent history of M. truncatula. [ABSTRACT FROM AUTHOR]

Background: The NODULATION RECEPTOR KINASE (NORK) gene encodes a Leucine-Rich Repeat (LRR)-containing receptor-like protein and controls the infection by symbiotic rhizobia and endomycorrhizal fungi in Legumes. The occurrence of numerous amino acid changes driven by directional selection has been reported in this gene, using a limited number of messenger RNA sequences, but the functional reason of these changes remains obscure. The Medicago genus, where changes in rhizobial associations have been previously examined, is a good model to test whether the evolution of NORK is influenced by rhizobial interactions. Results: We sequenced a region of 3610 nucleotides (encoding a 392 amino acid-long region of the NORK protein) in 32 Medicago species. We confirm that positive selection in NORK has occurred within the Medicago genus and find that the amino acid positions targeted by selection occur in sites outside of solvent-exposed regions in LRRs, and other sites in the N-terminal region of the protein. We tested if branches of the Medicago phylogeny where changes of rhizobial symbionts occurred displayed accelerated rates of amino acid substitutions. Only one branch out of five tested, leading to M. noeana, displays such a pattern. Among other branches, the most likely for having undergone positive selection is not associated with documented shift of rhizobial specificity. Conclusion: Adaptive changes in the sequence of the NORK receptor have involved the LRRs, but targeted different sites than in most previous studies of LRR proteins evolution. The fact that positive selection in NORK tends not to be associated to changes in rhizobial specificity indicates that this gene was probably not involved in evolving rhizobial preferences. Other explanations (e.g. coevolutionary arms race) must be tested to explain the adaptive evolution of NORK. [ABSTRACT FROM AUTHOR]

We investigate the genetic structure and molecular selection pattern of a sympatric population of Sinorhizobium meliloti and Sinorhizobium medicae. These bacteria fix nitrogen in association with plants of the genus Medicago. A set of 116 isolates were obtained from a soil sample, from root nodules of three groups of plants representing among-species, within-species and intraline diversity in the Medicago genus. Bacteria were characterized by sequencing at seven loci evenly distributed along the genome of both Sinorhizobium species, covering the chromosome and the two megaplasmids . We first test whether the diversity of host plants influence the bacterial diversity recovered. Using the same data set, we then analyse the selective pattern at each locus. There was no relationship between the diversity of Medicago plants that were used for sampling and the diversity of their symbionts. However, we found evidence of selection within each of the two main symbiotic regions, located on the two different megaplasmids. Purifying selection or a selective sweep was found to occur in the nod genomic region, which includes genes involved in nodulation specificity, whereas balancing selection was detected in the exo region, close to genes involved in exopolysaccharide production. Such pattern likely reflects the interaction between host plants and bacterial symbionts, with a possible conflict of interest between plants and cheater bacterial genotypes. Recombination appears to occur preferentially within and among loci located on megaplasmids, rather than within the chromosome. Thus, recombination may play an important role in resolving this conflict by allowing different selection patterns at different loci. [ABSTRACT FROM AUTHOR]

The author discusses the study conducted by M. L. Friesen on the symbiosis between nitrogen-fixing bacteria rhizobia and leguminous plants. The author states that the genetic variance of the effectiveness and competitiveness of rhizobia depends on the improved mutualism in both the plants and the bacteria. He also mentions that fitness alignment by host sanctions depends on the symbiosis as a whole in benefit of both the legume and rhizobium.

Rhizobium-root nodule symbiosis is generally considered to be unique for legumes. However, there is one exception, and that is Parasponia. In this nonlegume, the rhizobial nodule symbiosis evolved independently and is, as in legumes, induced by rhizobium Nod factors. We used Parasponia andersonii to identify genetic constraints underlying evolution of Nod factor signaling. Part of the signaling cascade, downstream of Nod factor perception, has been recruited from the more-ancient arbuscular endomycorrhizal symbiosis. However, legume Nod factor receptors that activate this common signaling pathway are not essential for arbuscular endomycorrhizae. Here, we show that in Parasponia a single Nod factor-like receptor is indispensable for both symbiotic interactions. Therefore, we conclude that the Nod factor perception mechanism also is recruited from the widespread endomycorrhizal symbiosis. [ABSTRACT FROM AUTHOR]