A point mutation and large deletion at the candidate avirulence locus AvrMlp7 in the poplar rust fungus correlate with poplar RMlp7 resistance breakdown.

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.

EggLib 3: A python package for population genetics and genomics.

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/.

Genomic Signatures of a Major Adaptive Event in the Pathogenic Fungus Melampsora larici-populina.

The recent availability of genome-wide sequencing techniques has allowed systematic screening for molecular signatures of adaptation, including in nonmodel organisms. Host-pathogen interactions constitute good models due to the strong selective pressures that they entail. We focused on an adaptive event which affected the poplar rust fungus Melampsora larici-populina when it overcame a resistance gene borne by its host, cultivated poplar. Based on 76 virulent and avirulent isolates framing narrowly the estimated date of the adaptive event, we examined the molecular signatures of selection. Using an array of genome scan methods based on different features of nucleotide diversity, we detected a single locus exhibiting a consistent pattern suggestive of a selective sweep in virulent individuals (excess of differentiation between virulent and avirulent samples, linkage disequilibrium, genotype-phenotype statistical association, and long-range haplotypes). Our study pinpoints a single gene and further a single amino acid replacement which may have allowed the adaptive event. Although our samples are nearly contemporary to the selective sweep, it does not seem to have affected genome diversity further than the immediate vicinity of the causal locus, which can be explained by a soft selective sweep (where selection acts on standing variation) and by the impact of recombination in mitigating the impact of selection. Therefore, it seems that properties of the life cycle of M. larici-populina, which entails both high genetic diversity and outbreeding, has facilitated its adaptation.

Evolution of morphological but not aggressiveness-related traits following a major resistance breakdown in the poplar rust fungus, Melampsora larici-populina.

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.

The escalatory Red Queen: Population extinction and replacement following arms race dynamics in poplar rust.

Host-parasite systems provide convincing examples of Red Queen co-evolutionary dynamics. Yet, a key process underscored in Van Valen's theory - that arms race dynamics can result in extinction - has never been documented. One reason for this may be that most sampling designs lack the breadth needed to illuminate the rapid pace of adaptation by pathogen populations. In this study, we used a 25-year temporal sampling to decipher the demographic history of a plant pathogen: the poplar rust fungus, Melampsora larici-populina. A major adaptive event occurred in 1994 with the breakdown of R7 resistance carried by several poplar cultivars widely planted in Western Europe since 1982. The corresponding virulence rapidly spread in M. larici-populina populations and nearly reached fixation in northern France, even on susceptible hosts. Using both temporal records of virulence profiles and temporal population genetic data, our analyses revealed that (i) R7 resistance breakdown resulted in the emergence of a unique and homogeneous genetic group, the so-called cultivated population, which predominated in northern France for about 20 years, (ii) selection for Vir7 individuals brought with it multiple other virulence types via hitchhiking, resulting in an overall increase in the population-wide number of virulence types and (iii) - above all - the emergence of the cultivated population superseded the initial population which predominated at the same place before R7 resistance breakdown. Our temporal analysis illustrates how antagonistic co-evolution can lead to population extinction and replacement, hence providing direct evidence for the escalation process which is at the core of Red Queen dynamics.

Intrapopulation genomics in a model mutualist: Population structure and candidate symbiosis genes under selection in Medicago truncatula.

Bottom-up evolutionary approaches, including geographically explicit population genomic analyses, have the power to reveal the mechanistic basis of adaptation. Here, we conduct a population genomic analysis in the model legume, Medicago truncatula, to characterize population genetic structure and identify symbiosis-related genes showing evidence of spatially variable selection. Using RAD-seq, we generated over 26,000 SNPs from 191 accessions from within three regions of the native range in Europe. Results from STRUCTURE analysis identify five distinct genetic clusters with divisions that separate east and west regions in the Mediterranean basin. Much of the genetic variation is maintained within sampling sites, and there is evidence for isolation by distance. Extensive linkage disequilibrium was identified, particularly within populations. We conducted genetic outlier analysis with FST -based genome scans and a Bayesian modeling approach (PCAdapt). There were 70 core outlier loci shared between these distinct methods with one clear candidate symbiosis related gene, DMI1. This work sets that stage for functional experiments to determine the important phenotypes that selection has acted upon and complementary efforts in rhizobium populations.

Genetic signatures of a range expansion in natura: when clones play leapfrog.

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.

A survey of genome-wide single nucleotide polymorphisms through genome resequencing in the Périgord black truffle (Tuber melanosporum Vittad.).

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.

Patterns of genomic variation in the poplar rust fungus Melampsora larici-populina identify pathogenesis-related factors.

Melampsora larici-populina is a fungal pathogen responsible for foliar rust disease on poplar trees, which causes damage to forest plantations worldwide, particularly in Northern Europe. The reference genome of the isolate 98AG31 was previously sequenced using a whole genome shotgun strategy, revealing a large genome of 101 megabases containing 16,399 predicted genes, which included secreted protein genes representing poplar rust candidate effectors. In the present study, the genomes of 15 isolates collected over the past 20 years throughout the French territory, representing distinct virulence profiles, were characterized by massively parallel sequencing to assess genetic variation in the poplar rust fungus. Comparison to the reference genome revealed striking structural variations. Analysis of coverage and sequencing depth identified large missing regions between isolates related to the mating type loci. More than 611,824 single-nucleotide polymorphism (SNP) positions were uncovered overall, indicating a remarkable level of polymorphism. Based on the accumulation of non-synonymous substitutions in coding sequences and the relative frequencies of synonymous and non-synonymous polymorphisms (i.e., PN/PS ), we identify candidate genes that may be involved in fungal pathogenesis. Correlation between non-synonymous SNPs in genes encoding secreted proteins (SPs) and pathotypes of the studied isolates revealed candidate genes potentially related to virulences 1, 6, and 8 of the poplar rust fungus.

Genome-wide patterns of segregation and linkage disequilibrium: the construction of a linkage genetic map of the poplar rust fungus Melampsora larici-populina.

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.

Amount of information needed for model choice in Approximate Bayesian Computation.

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.

Patterns of nucleotide diversity at photoperiod related genes in Norway spruce [Picea abies (L.) Karst].

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.

Evolution of a symbiotic receptor through gene duplications in the legume-rhizobium mutualism.

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 LysM-type receptor kinases. We investigated the evolutionary history of rhizobial LCO receptor genes MtLYK3-LjNFR1 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 MtLYK3-LjNFR1 genes. We hypothesize that their ancestor was more likely a mycorrhizal LCO receptor, than a defense-related receptor kinase.

Detecting selection along environmental gradients: analysis of eight methods and their effectiveness for outbreeding and selfing populations.

Thanks to genome-scale diversity data, present-day studies can provide a detailed view of how natural and cultivated species adapt to their environment and particularly to environmental gradients. However, due to their sensitivity, up-to-date studies might be more sensitive to undocumented demographic effects such as the pattern of migration and the reproduction regime. In this study, we provide guidelines for the use of popular or recently developed statistical methods to detect footprints of selection. We simulated 100 populations along a selective gradient and explored different migration models, sampling schemes and rates of self-fertilization. We investigated the power and robustness of eight methods to detect loci potentially under selection: three designed to detect genotype-environment correlations and five designed to detect adaptive differentiation (based on F(ST) or similar measures). We show that genotype-environment correlation methods have substantially more power to detect selection than differentiation-based methods but that they generally suffer from high rates of false positives. This effect is exacerbated whenever allele frequencies are correlated, either between populations or within populations. Our results suggest that, when the underlying genetic structure of the data is unknown, a number of robust methods are preferable. Moreover, in the simulated scenario we used, sampling many populations led to better results than sampling many individuals per population. Finally, care should be taken when using methods to identify genotype-environment correlations without correcting for allele frequency autocorrelation because of the risk of spurious signals due to allele frequency correlations between populations.

CCAAT-box binding transcription factors in plants: Y so many?

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.

Contrasted patterns of selective pressure in three recent paralogous gene pairs in the Medicago genus (L.).

BACKGROUND: Gene duplications are a molecular mechanism potentially mediating generation of functional novelty. However, the probabilities of maintenance and functional divergence of duplicated genes are shaped by selective pressures acting on gene copies immediately after the duplication event. The ratio of non-synonymous to synonymous substitution rates in protein-coding sequences provides a means to investigate selective pressures based on genic sequences. Three molecular signatures can reveal early stages of functional divergence between gene copies: change in the level of purifying selection between paralogous genes, occurrence of positive selection, and transient relaxed purifying selection following gene duplication. We studied three pairs of genes that are known to be involved in an interaction with symbiotic bacteria and were recently duplicated in the history of the Medicago genus (Fabaceae). We sequenced two pairs of polygalacturonase genes (Pg11-Pg3 and Pg11a-Pg11c) and one pair of auxine transporter-like genes (Lax2-Lax4) in 17 species belonging to the Medicago genus, and sought for molecular signatures of differentiation between copies. RESULTS: Selective histories revealed by these three signatures of molecular differentiation were found to be markedly different between each pair of paralogs. We found sites under positive selection in the Pg11 paralogs while Pg3 has mainly evolved under purifying selection. The most recent paralogs examined Pg11a and Pg11c, are both undergoing positive selection and might be acquiring new functions. Lax2 and Lax4 paralogs are both under strong purifying selection, but still underwent a temporary relaxation of purifying selection immediately after duplication. CONCLUSIONS: This study illustrates the variety of selective pressures undergone by duplicated genes and the effect of age of the duplication. We found that relaxation of selective constraints immediately after duplication might promote adaptive divergence.

Rhizobium-legume symbiosis shares an exocytotic pathway required for arbuscule formation.

Endosymbiotic interactions are characterized by the formation of specialized membrane compartments, by the host in which the microbes are hosted, in an intracellular manner. Two well-studied examples, which are of major agricultural and ecological importance, are the widespread arbuscular mycorrhizal symbiosis and the Rhizobium-legume symbiosis. In both symbioses, the specialized host membrane that surrounds the microbes forms a symbiotic interface, which facilitates the exchange of, for example, nutrients in a controlled manner and, therefore, forms the heart of endosymbiosis. Despite their key importance, the molecular and cellular mechanisms underlying the formation of these membrane interfaces are largely unknown. Recent studies strongly suggest that the Rhizobium-legume symbiosis coopted a signaling pathway, including receptor, from the more ancient arbuscular mycorrhizal symbiosis to form a symbiotic interface. Here, we show that two highly homologous exocytotic vesicle-associated membrane proteins (VAMPs) are required for formation of the symbiotic membrane interface in both interactions. Silencing of these Medicago VAMP72 genes has a minor effect on nonsymbiotic plant development and nodule formation. However, it blocks symbiosome as well as arbuscule formation, whereas root colonization by the microbes is not affected. Identification of these VAMP72s as common symbiotic regulators in exocytotic vesicle trafficking suggests that the ancient exocytotic pathway forming the periarbuscular membrane compartment has also been coopted in the Rhizobium-legume symbiosis.

EggLib: processing, analysis and simulation tools for population genetics and genomics.

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.

Purifying selection and molecular adaptation in the genome of Verminephrobacter, the heritable symbiotic bacteria of earthworms.

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.