Connecting genomic patterns of local adaptation and niche suitability in teosintes.

The central abundance hypothesis predicts that local adaptation is a function of the distance to the centre of a species' geographic range. To test this hypothesis, we gathered genomic diversity data from 49 populations, 646 individuals and 33,464 SNPs of two wild relatives of maize, the teosintes Zea mays ssp. parviglumis and Zea. mays. ssp. mexicana. We examined the association between the distance to their climatic and geographic centroids and the enrichment of SNPs bearing signals of adaptation. We identified candidate adaptive SNPs in each population by combining neutrality tests and cline analyses. By applying linear regression models, we found that the number of candidate SNPs is positively associated with niche suitability, while genetic diversity is reduced at the limits of the geographic distribution. Our results suggest that overall, populations located at the limit of the species' niches are adapting locally. We argue that local adaptation to this limit could initiate ecological speciation processes and facilitate adaptation to global change.

Signatures of local adaptation in lowland and highland teosintes from whole-genome sequencing of pooled samples.

Spatially varying selection triggers differential adaptation of local populations. Here, we mined the determinants of local adaptation at the genomewide scale in the two closest maize wild relatives, the teosintes Zea mays ssp parviglumis and ssp. mexicana. We sequenced 120 individuals from six populations: two lowland, two intermediate and two highland populations sampled along two altitudinal gradients. We detected 8 479 581 single nucleotide polymorphisms (SNPs) covered in the six populations with an average sequencing depth per site per population ranging from 17.0× to 32.2×. Population diversity varied from 0.10 to 0.15, and linkage disequilibrium decayed very rapidly. We combined two differentiation-based methods, and correlation of allele frequencies with environmental variables to detect outlier SNPs. Outlier SNPs displayed significant clustering. From clusters, we identified 47 candidate regions. We further modified a haplotype-based method to incorporate genotype uncertainties in haplotype calling, and applied it to candidate regions. We retrieved evidence for selection at the haplotype level in 53% of our candidate regions, and in 70% of the cases the same haplotype was selected in the two lowland or the two highland populations. We recovered a candidate region located within a previously characterized inversion on chromosome 1. We found evidence of a soft sweep at a locus involved in leaf macrohair variation. Finally, our results revealed frequent colocalization between our candidate regions and loci involved in the variation of traits associated with plant-soil interactions such as root morphology, aluminium and low phosphorus tolerance. Soil therefore appears to be a major driver of local adaptation in teosintes.

Contrasted patterns of selection since maize domestication on duplicated genes encoding a starch pathway enzyme.

Maize domestication from teosinte (Zea mays ssp. parviglumis) was accompanied by an increase of kernel size in landraces. Subsequent breeding has led to a diversification of kernel size and starch content among major groups of inbred lines. We aim at investigating the effect of domestication on duplicated genes encoding a key enzyme of the starch pathway, the ADP-glucose pyrophosphorylase (AGPase). Three pairs of paralogs encode the AGPase small (SSU) and large (LSU) subunits mainly expressed in the endosperm, the embryo and the leaf. We first validated the putative sequence of LSU(leaf) through a comparative expression assay of the six genes. Second, we investigated the patterns of molecular evolution on a 2 kb coding region homologous among the six genes in three panels: teosintes, landraces, and inbred lines. We corrected for demographic effects by relying on empirical distributions built from 580 previously sequenced ESTs. We found contrasted patterns of selection among duplicates: three genes exhibit patterns of directional selection during domestication (SSU(end), LSU(emb)) or breeding (LSU(leaf)), two exhibit patterns consistent with diversifying (SSU(leaf)) and balancing selection (SSU(emb)) accompanying maize breeding. While patterns of linkage disequilibrium did not reveal sign of coevolution between genes expressed in the same organ, we detected an excess of non-synonymous substitutions in the small subunit functional domains highlighting their role in AGPase evolution. Our results offer a different picture on AGPase evolution than the one depicted at the Angiosperm level and reveal how genetic redundancy can provide flexibility in the response to selection.

Apparent mutational hotspots and long distance linkage disequilibrium resulting from a bottleneck.

Genome wide patterns of nucleotide diversity and recombination reveal considerable variation including hotspots. Some studies suggest that these patterns are primarily dictated by individual locus history related at a broader scale to the population demographic history. Because bottlenecks have occurred in the history of numerous species, we undertook a simulation approach to investigate their impact on the patterns of aggregation of polymorphic sites and linkage disequilibrium (LD). We developed a new index (Polymorphism Aggregation Index) to characterize this aggregation and showed that variation in the density of polymorphic sites results from an interplay between the bottleneck scenario and the recombination rate. Under particular conditions, aggregation is maximized and apparent mutation hotspots resulting in a 50-fold increase in polymorphic sites density can occur. In similar conditions, long distance LD can be detected.

Patterns of DNA sequence polymorphism along chromosome 1 of maize (Zea mays ssp. mays L.).

We measured sequence diversity in 21 loci distributed along chromosome 1 of maize (Zea mays ssp. mays L.). For each locus, we sequenced a common sample of 25 individuals representing 16 exotic landraces and nine U.S. inbred lines. The data indicated that maize has an average of one single nucleotide polymorphism (SNP) every 104 bp between two randomly sampled sequences, a level of diversity higher than that of either humans or Drosophila melanogaster. A comparison of genetic diversity between the landrace and inbred samples showed that inbreds retained 77% of the level of diversity of landraces, on average. In addition, Tajima's D values suggest that the frequency distribution of polymorphisms in inbreds was skewed toward fewer rare variants. Tests for selection were applied to all loci, and deviations from neutrality were detected in three loci. Sequence diversity was heterogeneous among loci, but there was no pattern of diversity along the genetic map of chromosome 1. Nonetheless, diversity was correlated (r = 0.65) with sequence-based estimates of the recombination rate. Recombination in our sample was sufficient to break down linkage disequilibrium among SNPs. Intragenic linkage disequilibrium declines within 100-200 bp on average, suggesting that genome-wide surveys for association analyses require SNPs every 100-200 bp.