Adaptive introgression from maize has facilitated the establishment of teosinte as a noxious weed in Europe.

Global trade has considerably accelerated biological invasions. The annual tropical teosintes, the closest wild relatives of maize, were recently reported as new agricultural weeds in two European countries, Spain and France. Their prompt settlement under climatic conditions differing drastically from that of their native range indicates rapid genetic evolution. We performed a phenotypic comparison of French and Mexican teosintes under European conditions and showed that only the former could complete their life cycle during maize cropping season. To test the hypothesis that crop-to-wild introgression triggered such rapid adaptation, we used single nucleotide polymorphisms to characterize patterns of genetic variation in French, Spanish, and Mexican teosintes as well as in maize germplasm. We showed that both Spanish and French teosintes originated from Zea mays ssp. mexicana race "Chalco," a weedy teosinte from the Mexican highlands. However, introduced teosintes differed markedly from their Mexican source by elevated levels of genetic introgression from the high latitude Dent maize grown in Europe. We identified a clear signature of divergent selection in a region of chromosome 8 introgressed from maize and encompassing ZCN8, a major flowering time gene associated with adaptation to high latitudes. Moreover, herbicide assays and sequencing revealed that French teosintes have acquired herbicide resistance via the introgression of a mutant herbicide-target gene (ACC1) present in herbicide-resistant maize cultivars. Altogether, our results demonstrate that adaptive crop-to-wild introgression has triggered both rapid adaptation to a new climatic niche and acquisition of herbicide resistance, thereby fostering the establishment of an emerging noxious weed.

Genome-wide association study identifies favorable SNP alleles and candidate genes for frost tolerance in pea.

BACKGROUND: Frost is a limiting abiotic stress for the winter pea crop (Pisum sativum L.) and identifying the genetic determinants of frost tolerance is a major issue to breed varieties for cold northern areas. Quantitative trait loci (QTLs) have previously been detected from bi-parental mapping populations, giving an overview of the genome regions governing this trait. The recent development of high-throughput genotyping tools for pea brings the opportunity to undertake genetic association studies in order to capture a higher allelic diversity within large collections of genetic resources as well as to refine the localization of the causal polymorphisms thanks to the high marker density. In this study, a genome-wide association study (GWAS) was performed using a set of 365 pea accessions. Phenotyping was carried out by scoring frost damages in the field and in controlled conditions. The association mapping collection was also genotyped using an Illumina Infinium® BeadChip, which allowed to collect data for 11,366 single nucleotide polymorphism (SNP) markers. RESULTS: GWAS identified 62 SNPs significantly associated with frost tolerance and distributed over six of the seven pea linkage groups (LGs). These results confirmed 3 QTLs that were already mapped in multiple environments on LG III, V and VI with bi-parental populations. They also allowed to identify one locus, on LG II, which has not been detected yet and two loci, on LGs I and VII, which have formerly been detected in only one environment. Fifty candidate genes corresponding to annotated significant SNPs, or SNPs in strong linkage disequilibrium with the formers, were found to underlie the frost damage (FD)-related loci detected by GWAS. Additionally, the analyses allowed to define favorable haplotypes of markers for the FD-related loci and their corresponding accessions within the association mapping collection. CONCLUSIONS: This study led to identify FD-related loci as well as corresponding favorable haplotypes of markers and representative pea accessions that might to be used in winter pea breeding programs. Among the candidate genes highlighted at the identified FD-related loci, the results also encourage further attention to the presence of C-repeat Binding Factors (CBF) as potential genetic determinants of the frost tolerance locus on LG VI.

Genetic diversity and trait genomic prediction in a pea diversity panel.

BACKGROUND: Pea (Pisum sativum L.), a major pulse crop grown for its protein-rich seeds, is an important component of agroecological cropping systems in diverse regions of the world. New breeding challenges imposed by global climate change and new regulations urge pea breeders to undertake more efficient methods of selection and better take advantage of the large genetic diversity present in the Pisum sativum genepool. Diversity studies conducted so far in pea used Simple Sequence Repeat (SSR) and Retrotransposon Based Insertion Polymorphism (RBIP) markers. Recently, SNP marker panels have been developed that will be useful for genetic diversity assessment and marker-assisted selection. RESULTS: A collection of diverse pea accessions, including landraces and cultivars of garden, field or fodder peas as well as wild peas was characterised at the molecular level using newly developed SNP markers, as well as SSR markers and RBIP markers. The three types of markers were used to describe the structure of the collection and revealed different pictures of the genetic diversity among the collection. SSR showed the fastest rate of evolution and RBIP the slowest rate of evolution, pointing to their contrasted mode of evolution. SNP markers were then used to predict phenotypes -the date of flowering (BegFlo), the number of seeds per plant (Nseed) and thousand seed weight (TSW)- that were recorded for the collection. Different statistical methods were tested including the LASSO (Least Absolute Shrinkage ans Selection Operator), PLS (Partial Least Squares), SPLS (Sparse Partial Least Squares), Bayes A, Bayes B and GBLUP (Genomic Best Linear Unbiased Prediction) methods and the structure of the collection was taken into account in the prediction. Despite a limited number of 331 markers used for prediction, TSW was reliably predicted. CONCLUSION: The development of marker assisted selection has not reached its full potential in pea until now. This paper shows that the high-throughput SNP arrays that are being developed will most probably allow for a more efficient selection in this species.

Engineering multi-degrading bacterial communities to bioremediate soils contaminated with pesticides residues.

Parallel to the important use of pesticides in conventional agriculture there is a growing interest for green technologies to clear contaminated soil from pesticides and their degradation products. Bioaugmentation i. e. the inoculation of degrading micro-organisms in polluted soil, is a promising method still in needs of further developments. Specifically, improvements in the understanding of how degrading microorganisms must overcome abiotic filters and interact with the autochthonous microbial communities are needed in order to efficiently design bioremediation strategies. Here we designed a protocol aiming at studying the degradation of two herbicides, glyphosate (GLY) and isoproturon (IPU), via experimental modifications of two source bacterial communities. We used statistical methods stemming from genomic prediction to link community composition to herbicides degradation potentials. Our approach proved to be efficient with correlation estimates over 0.8 - between model predictions and measured pesticide degradation values. Multi-degrading bacterial communities were obtained by coalescing bacterial communities with high GLY or IPU degradation ability based on their community-level properties. Finally, we evaluated the efficiency of constructed multi-degrading communities to remove pesticide contamination in a different soil. While results are less clear in the case of GLY, we showed an efficient transfer of degrading capacities towards the receiving soil even at relatively low inoculation levels in the case of IPU. Altogether, we developed an innovative protocol for building multi-degrading simplified bacterial communities with the help of genomic prediction tools and coalescence, and proved their efficiency in a contaminated soil.

Genomic consequences of transitions from cross- to self-fertilization on the efficacy of selection in three independently derived selfing plants.

BACKGROUND: Transitions from cross- to self-fertilization are associated with increased genetic drift rendering weakly selected mutations effectively neutral. The effect of drift is predicted to reduce selective constraints on amino acid sequences of proteins and relax biased codon usage. We investigated patterns of nucleotide variation to assess the effect of inbreeding on the accumulation of deleterious mutations in three independently evolved selfing plants. Using high-throughput sequencing, we assembled the floral transcriptomes of four individuals of Eichhornia (Pontederiaceae); these included one outcrosser and two independently derived selfers of E. paniculata, and E. paradoxa, a selfing outgroup. The dataset included ~8000 loci totalling ~3.5 Mb of coding DNA. RESULTS: Tests of selection were consistent with purifying selection constraining evolution of the transcriptome. However, we found an elevation in the proportion of non-synonymous sites that were potentially deleterious in the E. paniculata selfers relative to the outcrosser. Measurements of codon usage in high versus low expression genes demonstrated reduced bias in both E. paniculata selfers. CONCLUSIONS: Our findings are consistent with a small reduction in the efficacy of selection on protein sequences associated with transitions to selfing, and reduced selection in selfers on synonymous changes that influence codon usage.

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.

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

De novo sequence assembly and characterization of the floral transcriptome in cross- and self-fertilizing plants.

BACKGROUND: The shift from cross-fertilization to predominant self-fertilization is among the most common evolutionary transitions in the reproductive biology of flowering plants. Increased inbreeding has important consequences for floral morphology, population genetic structure and genome evolution. The transition to selfing is usually characterized by a marked reduction in flower size and the loss of traits involved in pollinator attraction and the avoidance of self-fertilization. Here, we use short-read sequencing to assemble, de novo, the floral transcriptomes of three genotypes of Eichhornia paniculata, including an outcrosser and two genotypes from independently derived selfers, and a single genotype of the sister species E. paradoxa. By sequencing mRNA from tissues sampled at various stages of flower development, our goal was to sequence and assemble the floral transcriptome and identify differential patterns of gene expression. RESULTS: Our 24 Mbp assembly resulted in ~27,000 contigs that averaged ~900 bp in length. All four genotypes had highly correlated gene expression, but the three E. paniculata genotypes were more correlated with one another than each was to E. paradoxa. Our analysis identified 269 genes associated with floral development, 22 of which were differentially expressed in selfing lineages relative to the outcrosser. Many of the differentially expressed genes affect floral traits commonly altered in selfing plants and these represent a set of potential candidate genes for investigating the evolution of the selfing syndrome. CONCLUSIONS: Our study is among the first to demonstrate the use of Illumina short read sequencing for de novo transcriptome assembly in non-model species, and the first to implement this technology for comparing floral transcriptomes in outcrossing and selfing plants.

The population genomics of plant adaptation.

There has been an enormous increase in the amount of data on DNA sequence polymorphism available for many organisms in the last decade. New sequencing technologies provide great potential for investigating natural selection in plants using population genomic approaches. However, plant populations frequently show significant departures from the assumptions of standard models used to detect selection and many forms of directional selection do not fit with classical population genetics theory. Here, we explore the extent to which plant populations show departures from standard model assumptions, and the implications this has for detecting selection on molecular variation. A growing number of multilocus studies of nucleotide variation suggest that changes in population size, particularly bottlenecks, and strong subdivision may be common in plants. This demographic variation presents important challenges for models used to infer selection. In addition, selection from standing genetic variation and multiple independent adaptive substitutions can further complicate efforts to understand the nature of selection. We discuss emerging patterns from plant studies and propose that, rather than treating population history as a nuisance variable when testing for selection, the interaction between demography and selection is of fundamental importance for evolutionary studies of plant populations using molecular data.

Patterns of Genetic Structure and Linkage Disequilibrium in a Large Collection of Pea Germplasm.

Pea (Pisum sativum, L.) is a major pulse crop used both for animal and human alimentation. Owing to its association with nitrogen-fixing bacteria, it is also a valuable component for low-input cropping systems. To evaluate the genetic diversity and the scale of linkage disequilibrium (LD) decay in pea, we genotyped a collection of 917 accessions, gathering elite cultivars, landraces, and wild relatives using an array of ∼13,000 single nucleotide polymorphisms (SNP). Genetic diversity is broadly distributed across three groups corresponding to wild/landraces peas, winter types, and spring types. At a finer subdivision level, genetic groups relate to local breeding programs and type usage. LD decreases steeply as genetic distance increases. When considering subsets of the data, LD values can be higher, even if the steep decay remains. We looked for genomic regions exhibiting high level of differentiation between wild/landraces, winter, and spring pea, respectively. Two regions on linkage groups 5 and 6 containing 33 SNPs exhibit stronger differentiation between winter and spring peas than would be expected under neutrality. Interestingly, QTL for resistance to cold acclimation and frost resistance have been identified previously in the same regions.

Co-inoculation of a Pea Core-Collection with Diverse Rhizobial Strains Shows Competitiveness for Nodulation and Efficiency of Nitrogen Fixation Are Distinct traits in the Interaction.

Pea forms symbiotic nodules with Rhizobium leguminosarum sv. viciae (Rlv). In the field, pea roots can be exposed to multiple compatible Rlv strains. Little is known about the mechanisms underlying the competitiveness for nodulation of Rlv strains and the ability of pea to choose between diverse compatible Rlv strains. The variability of pea-Rlv partner choice was investigated by co-inoculation with a mixture of five diverse Rlv strains of a 104-pea collection representative of the variability encountered in the genus Pisum. The nitrogen fixation efficiency conferred by each strain was determined in additional mono-inoculation experiments on a subset of 18 pea lines displaying contrasted Rlv choice. Differences in Rlv choice were observed within the pea collection according to their genetic or geographical diversities. The competitiveness for nodulation of a given pea-Rlv association evaluated in the multi-inoculated experiment was poorly correlated with its nitrogen fixation efficiency determined in mono-inoculation. Both plant and bacterial genetic determinants contribute to pea-Rlv partner choice. No evidence was found for co-selection of competitiveness for nodulation and nitrogen fixation efficiency. Plant and inoculant for an improved symbiotic association in the field must be selected not only on nitrogen fixation efficiency but also for competitiveness for nodulation.

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.

Pedigree-free estimates of heritability in the wild: promising prospects for selfing populations.

Estimating the genetic variance available for traits informs us about a population's ability to evolve in response to novel selective challenges. In selfing species, theory predicts a loss of genetic diversity that could lead to an evolutionary dead-end, but empirical support remains scarce. Genetic variability in a trait is estimated by correlating the phenotypic resemblance with the proportion of the genome that two relatives share identical by descent ('realized relatedness'). The latter is traditionally predicted from pedigrees (Φ A : expected value) but can also be estimated using molecular markers (average number of alleles shared). Nevertheless, evolutionary biologists, unlike animal breeders, remain cautious about using marker-based relatedness coefficients to study complex phenotypic traits in populations. In this paper, we review published results comparing five different pedigree-free methods and use simulations to test individual-based models (hereafter called animal models) using marker-based relatedness coefficients, with a special focus on the influence of mating systems. Our literature review confirms that Ritland's regression method is unreliable, but suggests that animal models with marker-based estimates of relatedness and genomic selection are promising and that more testing is required. Our simulations show that using molecular markers instead of pedigrees in animal models seriously worsens the estimation of heritability in outcrossing populations, unless a very large number of loci is available. In selfing populations the results are less biased. More generally, populations with high identity disequilibrium (consanguineous or bottlenecked populations) could be propitious for using marker-based animal models, but are also more likely to deviate from the standard assumptions of quantitative genetics models (non-additive variance).