Genomic microsatellite adaptive divergence of wild barley by microclimatic stress in ‘Evolution Canyon’, Israel

The effect of habitat and climatic on microsatellite diversity and allele length variation

Microsatellite loci have been shown to vary according to environment. We studied allelic length variation and diversity in eight microsatellite loci along a sharp climatic and habitat gradient in Israel, using the eastern spadefoot toad (Pelobates syriacus) as our model system. We found a gradual increase in allele lengths from north to south. We used a distance-based redundancy analysis (dbRDA) to associate between allele length and habitat and climatic measures and found that geography and annual climate explained the significant variation in allele length. We also used additional measurements pertaining to demography, heterozygosity and allelic diversity to explore four different hypotheses that might explain the variations in allele length. Our results suggest that the changes we observed in allele lengths may not be purely random but could be influenced by the differential mutation rate and/or local environmental conditions operating at the different locations.

Physical geography, isolation by distance and environmental variables shape genomic variation of wild barley (Hordeum vulgare L. ssp. spontaneum) in the Southern Levant

Determining the extent of genetic variation that reflects local adaptation in crop-wild relatives is of interest for the purpose of identifying useful genetic diversity for plant breeding. We investigated the association of genomic variation with geographical and environmental factors in wild barley (Hordeum vulgare L. ssp. spontaneum) populations of the Southern Levant using genotyping by sequencing (GBS) of 244 accessions in the Barley 1K+ collection. The inference of population structure resulted in four genetic clusters that corresponded to eco-geographical habitats and a significant association between lower gene flow rates and geographical barriers, e.g. the Judaean Mountains and the Sea of Galilee. Redundancy analysis (RDA) revealed that spatial autocorrelation explained 45% and environmental variables explained 15% of total genomic variation. Only 4.5% of genomic variation was solely attributed to environmental variation if the component confounded with spatial autocorrelation was excluded. A synthetic environmental variable combining latitude, solar radiation, and accumulated precipitation explained the highest proportion of genomic variation (3.9%). When conditioned on population structure, soil water capacity was the most important environmental variable explaining 1.18% of genomic variation. Genome scans with outlier analysis and genome-environment association studies were conducted to identify adaptation signatures. RDA and outlier methods jointly detected selection signatures in the pericentromeric regions, which have reduced recombination, of the chromosomes 3H, 4H, and 5H. However, selection signatures mostly disappeared after correction for population structure. In conclusion, adaptation to the highly diverse environments of the Southern Levant over short geographical ranges had a limited effect on the genomic diversity of wild barley. This highlighted the importance of nonselective forces in genetic differentiation.

Transcribed microsatellite allele lengths are often correlated with gene expression in natural sunflower populations

Microsatellites are common in genomes of most eukaryotic species. Due to their high mutability, an adaptive role for microsatellites has been considered. However, little is known concerning the contribution of microsatellites towards phenotypic variation. We used populations of the common sunflower (Helianthus annuus) at two latitudes to quantify the effect of microsatellite allele length on phenotype at the level of gene expression. We conducted a common garden experiment with seed collected from sunflower populations in Kansas and Oklahoma followed by an RNA-Seq experiment on 95 individuals. The effect of microsatellite allele length on gene expression was assessed across 3,325 microsatellites that could be consistently scored. Our study revealed 479 microsatellites at which allele length significantly correlates with gene expression (eSTRs). When irregular allele sizes not conforming to the motif length were removed, the number of eSTRs rose to 2,379. The percentage of variation in gene expression explained by eSTRs ranged from 1%-86% when controlling for population and allele-by-population interaction effects at the 479 eSTRs. Of these eSTRs, 70.4% are in untranslated regions (UTRs). A gene ontology (GO) analysis revealed that eSTRs are significantly enriched for GO terms associated with cis- and trans-regulatory processes. Our findings suggest that a substantial number of transcribed microsatellites can influence gene expression.

Sympatric speciation of wild emmer wheat driven by ecology and chromosomal rearrangements

In plants, the mechanism for ecological sympatric speciation (SS) is little known. Here, after ruling out the possibility of secondary contact, we show that wild emmer wheat, at the microclimatically divergent microsite of "Evolution Canyon" (EC), Mt. Carmel, Israel, underwent triple SS. Initially, it split following a bottleneck of an ancestral population, and further diversified to three isolated populations driven by disruptive ecological selection. Remarkably, two postzygotically isolated populations (SFS1 and SFS2) sympatrically branched within an area less than 30 m at the tropical hot and dry savannoid south-facing slope (SFS). A series of homozygous chromosomal rearrangements in the SFS1 population caused hybrid sterility with the SFS2 population. We demonstrate that these two populations developed divergent adaptive mechanisms against severe abiotic stresses on the tropical SFS. The SFS2 population evolved very early flowering, while the SFS1 population alternatively evolved a direct tolerance to irradiance by improved ROS scavenging activity that potentially accounts for its evolutionary fate with unstable chromosome status. Moreover, a third prezygotically isolated sympatric population adapted on the abutting temperate, humid, cool, and forested north-facing slope (NFS), separated by 250 m from the SFS wild emmer wheat populations. The NFS population evolved multiple resistant loci to fungal diseases, including powdery mildew and stripe rust. Our study illustrates how plants sympatrically adapt and speciate under disruptive ecological selection of abiotic and biotic stresses.

Genomic adaptation to drought in wild barley is driven by edaphic natural selection at the Tabigha Evolution Slope

Ecological divergence at a microsite suggests adaptive evolution, and this study examined two abutting wild barley populations, each 100 m across, differentially adapted to drought tolerance on two contrasting soil types, Terra Rossa and basalt at the Tabigha Evolution Slope, Israel. We resequenced the genomes of seven and six wild barley genotypes inhabiting the Terra Rossa and basalt soils, respectively, and identified a total of 69,192,653 single-nucleotide variants (SNVs) and insertions/deletions in comparison with a reference barley genome. Comparative genomic analysis between these abutting wild barley populations involved 19,615,087 high-quality SNVs. The results revealed dramatically different selection sweep regions relevant to drought tolerance driven by edaphic natural selection within 2,577 selected genes in these regions, including key drought-responsive genes associated with ABA synthesis and degradation (such as Cytochrome P450 protein) and ABA receptor complex (such as PYL2, SNF1-related kinase). The genetic diversity of the wild barley population inhabiting Terra Rossa soil is much higher than that from the basalt soil. Additionally, we identified different sets of genes for drought adaptation in the wild barley populations from Terra Rossa soil and from wild barley populations from Evolution Canyon I at Mount Carmel. These genes are associated with abscisic acid signaling, signaling and metabolism of reactive oxygen species, detoxification and antioxidative systems, rapid osmotic adjustment, and deep root morphology. The unique mechanisms for drought adaptation of the wild barley from the Tabigha Evolution Slope may be useful for crop improvement, particularly for breeding of barley cultivars with high drought tolerance.

Geographical and environmental determinants of the genetic structure of wild barley in southeastern Anatolia

Despite the global value of barley, compared to its wild progenitor, genetic variation in this crop has been drastically reduced due to the process of domestication, selection and improvement. In the medium term, this will negatively impact both the vulnerability and yield stability of barley against biotic and abiotic stresses under climate change. Returning to the crop wild relatives (CWR) as sources of new and beneficial alleles is a clear option for enhancing the resilience of diversity and adaptation to climate change. Southeastern Anatolia constitutes an important part of the natural distribution of wild barley in the Fertile Crescent where important crops were initially domesticated. In this study, we investigated genetic diversity in a comprehensive collection of 281 geo-referenced wild barley individuals from 92 collection sites with sample sizes ranging from 1 to 9 individuals per site, collected from southeastern Anatolia and 131 domesticated genotypes from 49 different countries using 40 EST-SSR markers. A total of 375 alleles were detected across entire collection, of which 283 were carried by domesticated genotypes and 316 alleles were present in the wild gene pool. The number of unique alleles in the wild and in the domesticated gene pool was 92 and 59, respectively. The population structure at K = 3 suggested two groups of wild barley namely G1-W consisting wild barley genotypes from the western part and G1-E comprising those mostly from the eastern part of the study area, with a sharp separation from the domesticated gene pool. The geographic and climatic factors jointly showed significant effects on the distribution of wild barley. Using a Latent Factor Mixed Model, we identified four candidate loci potentially involved in adaptation of wild barley to three environmental factors: temperature seasonality, mean temperature of driest quarter, and precipitation of coldest quarter. These loci are probably the targets of genomic regions, with potential roles against abiotic stresses.

Adaptive evolution of benzoxazinoids in wild emmer wheat, Triticum dicoccoides, at "Evolution Canyon", Mount Carmel, Israel

Background

"Evolution Canyon" (ECI) at Lower Nahal Oren, Mount Carmel, Israel, is an optimal natural microscale model for unraveling evolution-in-action, highlighting the evolutionary processes of biodiversity evolution, adaptation, and incipient sympatric speciation. A major model organism in ECI is the tetraploid wild emmer wheat, Triticum dicoccoides (TD), the progenitor of cultivated emmer and durum wheat. TD displays dramatic interslope adaptive evolutionary divergence on the tropical, savannoid-hot and dry south-facing, "African" slope (AS), and on the temperate, forested, cool and humid, north-facing, "European" slope (ES), separated on average by 250 m. From the perspective of chemical evolution and metabolomics, it is important to unravel interslope divergence in biologically relevant secondary metabolites between the abutting slope populations. Here, in TD we examined hydroxamic acid (Hx), which is a family of secondary cereal metabolites, and plays a major role in defending the plant against fungi, insects and weeds.

Results

Our examination revealed that higher concentrations of DIBOA and DIMBOA were found in seedlings growing in the same greenhouse from seeds collected from the cool and humid forested ES, whereas the seedlings of seeds collected from the savannoid AS (both in root and shoot tissues), showed no DIMBOA. Remarkably, only DIBOA appears in both shoots and roots of the AS seedlings. It rises to a peak and then decreases in both organs and in seedlings from both slopes. The DIMBOA, which appears only in the ES seedlings, rises to a peak and decreases in the shoot, but increased and remained in a plateau in the root, till the end of the experiment.

Conculsions/Significance

The results suggest stronger genetic resistance of defense compounds DIBOA and DIMBOA against biotic stresses (fungi and other pathogens) by ES seedlings. However, AS seedlings responded earlier but were to the same biotic stresses. The genetic difference found in AS seedlings was caused by the main adaptive selection in AS, which was against climatic, abiotic stresses, and was weaker, or not at all, against biotic stresses. The distinct genetic interslope differences appear important and is very significant and are elaborated in the discussion.

Sympatric speciation of spiny mice, Acomys, unfolded transcriptomically at Evolution Canyon, Israel

Spiny mice, Acomys cahirinus, colonized Israel 30,000 y ago from dry tropical Africa and inhabited rocky habitats across Israel. Earlier, we had shown by mtDNA that A. cahirinus incipiently sympatrically speciates at Evolution Canyon I (EC I) in Mount Carmel, Israel because of microclimatic interslope divergence. The EC I microsite consists of a dry and hot savannoid "African" slope (AS) and an abutting humid and cool-forested "European" slope (ES). Here, we substantiate incipient SS in A. cahirinus at EC I based on the entire transcriptome, showing that multiple slope-specific adaptive complexes across the transcriptome result in two divergent clusters. Tajima's D distribution of the abutting Acomys interslope populations shows that the ES population is under stronger positive selection, whereas the AS population is under balancing selection, harboring higher genetic polymorphisms. Considerable sites of the two populations were differentiated with a coefficient of FST = 0.25-0.75. Remarkably, 24 and 37 putatively adaptively selected genes were detected in the AS and ES populations, respectively. The AS genes involved DNA repair, growth arrest, neural cell differentiation, and heat-shock proteins adapting to the local AS stresses of high solar radiation, drought, and high temperature. In contrast, the ES genes involved high ATP associated with energetics stress. The sharp ecological interslope divergence led to strong slope-specific selection overruling the interslope gene flow. Earlier tests suggested slope-specific mate choice. Habitat interslope-adaptive selection across the transcriptome and mate choice substantiate sympatric speciation (SS), suggesting its prevalence at EC I and commonality in nature.

Developing informative microsatellite markers for non-model species using reference mapping against a model species' genome

Microsatellites have a wide range of applications from behavioral biology, evolution, to agriculture-based breeding programs. The recent progress in the next-generation sequencing technologies and the rapidly increasing number of published genomes may greatly enhance the current applications of microsatellites by turning them from anonymous to informative markers. Here we developed an approach to anchor microsatellite markers of any target species in a genome of a related model species, through which the genomic locations of the markers, along with any functional genes potentially linked to them, can be revealed. We mapped the shotgun sequence reads of a non-model rodent species Apodemus semotus against the genome of a model species, Mus musculus, and presented 24 polymorphic microsatellite markers with detailed background information for A. semotus in this study. The developed markers can be used in other rodent species, especially those that are closely related to A. semotus or M. musculus. Compared to the traditional approaches based on DNA cloning, our approach is likely to yield more loci for the same cost. This study is a timely demonstration of how a research team can efficiently generate informative (neutral or function-associated) microsatellite markers for their study species and unique biological questions.

Natural selection causes adaptive genetic resistance in wild emmer wheat against powdery mildew at "Evolution Canyon" microsite, Mt. Carmel, Israel

Background

“Evolution Canyon” (ECI) at Lower Nahal Oren, Mount Carmel, Israel, is an optimal natural microscale model for unraveling evolution in action highlighting the basic evolutionary processes of adaptation and speciation. A major model organism in ECI is wild emmer, Triticum dicoccoides, the progenitor of cultivated wheat, which displays dramatic interslope adaptive and speciational divergence on the tropical-xeric “African” slope (AS) and the temperate-mesic “European” slope (ES), separated on average by 250 m.

Methods

We examined 278 single sequence repeats (SSRs) and the phenotype diversity of the resistance to powdery mildew between the opposite slopes. Furthermore, 18 phenotypes on the AS and 20 phenotypes on the ES, were inoculated by both Bgt E09 and a mixture of powdery mildew races.

Results

In the experiment of genetic diversity, very little polymorphism was identified intra-slope in the accessions from both the AS or ES. By contrast, 148 pairs of SSR primers (53.23%) amplified polymorphic products between the phenotypes of AS and ES. There are some differences between the two wild emmer wheat genomes and the inter-slope SSR polymorphic products between genome A and B. Interestingly, all wild emmer types growing on the south-facing slope (SFS=AS) were susceptible to a composite of Blumeria graminis, while the ones growing on the north-facing slope (NFS=ES) were highly resistant to Blumeria graminis at both seedling and adult stages.

Conclusion/Significance

Remarkable inter-slope evolutionary divergent processes occur in wild emmer wheat, T. dicoccoides at EC I, despite the shot average distance of 250 meters. The AS, a dry and hot slope, did not develop resistance to powdery mildew, whereas the ES, a cool and humid slope, did develop resistance since the disease stress was strong there. This is a remarkable demonstration in host-pathogen interaction on how resistance develops when stress causes an adaptive result at a micro-scale distance.

The effects of microsatellite selection on linked sequence diversity

The genome-wide scan for selection is an important method for identifying loci involved in adaptive evolution. However, theory that underlies standard scans for selection assumes a simple mutation model. In particular, recurrent mutation of the selective target is not considered. Although this assumption is reasonable for single-nucleotide variants (SNVs), a microsatellite targeted by selection will reliably violate this assumption due to high mutation rate. Moreover, the mutation rate of microsatellites is generally high enough to ensure that recurrent mutation is pervasive rather than occasional. It is therefore unclear if positive selection targeting microsatellites can be detected using standard scanning statistics. Examples of functional variation at microsatellites underscore the significance of understanding the genomic effects of microsatellite selection. Here, we investigate the joint effects of selection and complex mutation on linked sequence diversity, comparing simulations of microsatellite selection and SNV-based selective sweeps. We find that selection on microsatellites is generally difficult to detect using popular summaries of the site frequency spectrum, and, under certain conditions, using popular methods such as the integrated haplotype statistic and SweepFinder. However, comparisons of the number of haplotypes (K) and segregating sites (S) often provide considerable power to detect selection on microsatellites. We apply this knowledge to a scan of autosomes in the human CEU population (CEPH population sampled from Utah). In addition to the most commonly reported targets of selection in European populations, we identify numerous novel genomic regions that bear highly anomalous haplotype configurations. Using one of these regions—intron 1 of MAGI2—as an example, we show that the anomalous configuration is coincident with a perfect CA repeat of length 22. We conclude that standard genome-wide scans will commonly fail to detect mutationally complex targets of selection but that comparisons of K and S will, in many cases, facilitate their identification.

Genetic divergence in domesticated and non-domesticated gene regions of barley chromosomes

Little is known about the genetic divergence in the chromosomal regions with domesticated and non-domesticated genes. The objective of our study is to examine the effect of natural selection on shaping genetic diversity of chromosome region with domesticated and non-domesticated genes in barley using 110 SSR markers. Comparison of the genetic diversity loss between wild and cultivated barley for each chromosome showed that chromosome 5H had the highest divergence of 35.29%, followed by 3H, 7H, 4H, 2H, 6H. Diversity ratio was calculated as (diversity of wild type – diversity of cultivated type)/diversity of wild type×100%. It was found that diversity ratios of the domesticated regions on 5H, 1H and 7H were higher than those of non-domesticated regions. Diversity ratio of the domesticated region on 2H and 4H is similar to that of non-domesticated region. However, diversity ratio of the domesticated region on 3H is lower than that of non-domesticated region. Averaged diversity among six chromosomes in domesticated region was 33.73% difference between wild and cultivated barley, and was 27.56% difference in the non-domesticated region. The outcome of this study advances our understanding of the evolution of crop chromosomes.

Genetic characterization of wild barley populations (Hordeum vulgare ssp. spontaneum) from Kazakhstan based on genome wide SNP analysis

The wild ancestral form of barley, Hordeum vulgare ssp. spontaneum, is a valuable source for gene enrichment of cultivated barley. The purpose of this work was to study the area of distribution as well as the extent and structure of genetic variation of wild barley populations grown in Kazakhstan. It was found that distribution of wild barley populations in Kazakhstan is restricted to the most southern province. A genome wide single nucleotide polymorphism (SNP) analysis was performed in order to study the level of the genetic diversity in 96 accessions representing 14 wild barley populations from Kazakhstan and 25 accessions from the Middle East which is the center of diversity of this subspecies. The oligonucleotide pooled assay was used to genotype 384 SNPs distributed throughout the genome. In total 233 polymorphic SNPs were selected for further statistical analysis. The level of genetic diversity of wild barley populations from Kazakhstan was predictably narrower (He = 0.19 ± 0.01) in comparison with wild barley samples from the Middle East (He = 0.29 ± 0.01). The results suggested that H. vulgare ssp. spontaneum populations in Kazakhstan probably represent a recent spread of a limited number of plants from the primary distribution area and might be well adapted to winter low temperature.

Evolution in action: adaptation and incipient sympatric speciation with gene flow across life at “Evolution Canyon”, Israel

Various major evolutionary problems are still open, controversial or unsettled. These include even the basic evolutionary processes of adaptation and speciation. The “Evolution Canyon” model is a microscale natural laboratory that can highlight some of the basic problems requiring clarification (Nevo list of “Evolution Canyon” publications at http://evolution.haifa.ac.il). This is especially true if an interdisciplinary approach is practiced including ecological functional genomics, transcriptomics, proteomics, metabolomics and phenomics. Here I overview and reanalyze the incipient sympatric adaptive ecological speciation of five model organisms at “Evolution Canyon”, across life: the soil bacterium, Bacillus simplex; wild barley, the progenitor of cultivated barley, Hordeum spontaneum; the tiny beetle Oryzaephilus surinamensis; the cosmopolitan fruit-fly, Drosophila melanogaster, and the Africa-originated spiny mouse, Acomys cahirinus. All five models of organisms display evolution in action of microclimatic adaptation and incipient sympatric adaptive ecological speciation on the tropical and temperate abutting slopes, separated on average by only 250 meters. Some distant species converge in their micro-climatic adaptations to the hot and dry “African”, south-facing slope (SFS or AS) and to the cool and humid “European”, north-facing slope (NSF or ES). Natural selection overrules ongoing inter-slope gene-flow between the free interbreeding populations within and between slopes, and leads to adaptive incipient sympatric ecological speciation on the dramatically opposite abutting xeric savannoid and mesic forested slopes.

DNA sequence variation of wild barley Hordeum spontaneum (L.) across environmental gradients in Israel

Wild barley Hordeum spontaneum (L.) shows a wide geographic distribution and ecological diversity. A key question concerns the spatial scale at which genetic differentiation occurs and to what extent it is driven by natural selection. The Levant region exhibits a strong ecological gradient along the North-South axis, with numerous small canyons in an East-West direction and with small-scale environmental gradients on the opposing North- and South-facing slopes. We sequenced 34 short genomic regions in 54 accessions of wild barley collected throughout Israel and from the opposing slopes of two canyons. The nucleotide diversity of the total sample is 0.0042, which is about two-thirds of a sample from the whole species range (0.0060). Thirty accessions collected at 'Evolution Canyon' (EC) at Nahal Oren, close to Haifa, have a nucleotide diversity of 0.0036, and therefore harbor a large proportion of the genetic diversity. There is a high level of genetic clustering throughout Israel and within EC, which roughly differentiates the slopes. Accessions from the hot and dry South-facing slope have significantly reduced genetic diversity and are genetically more distinct from accessions from the North-facing slope, which are more similar to accessions from other regions in Northern Israel. Statistical population models indicate that wild barley within the EC consist of three separate genetic clusters with substantial gene flow. The data indicate a high level of population structure at large and small geographic scales that shows isolation-by-distance, and is also consistent with ongoing natural selection contributing to genetic differentiation at a small geographic scale.

Genome diversity in wild grasses under environmental stress

Patterns of diversity distribution in the Isa defense locus in wild-barley populations suggest adaptive selection at this locus. The extent to which environmental selection may act at additional nuclear-encoded defense loci and within the whole chloroplast genome has now been examined by analyses in two grass species. Analysis of genetic diversity in wild barley (Hordeum spontaneum) defense genes revealed much greater variation in biotic stress-related genes than abiotic stress-related genes. Genetic diversity at the Isa defense locus in wild populations of weeping ricegrass [Microlaena stipoides (Labill.) R. Br.], a very distant wild-rice relative, was more diverse in samples from relatively hotter and drier environments, a phenomenon that reflects observations in wild barley populations. Whole-chloroplast genome sequences of bulked weeping ricegrass individuals sourced from contrasting environments showed higher levels of diversity in the drier environment in both coding and noncoding portions of the genome. Increased genetic diversity may be important in allowing plant populations to adapt to greater environmental variation in warmer and drier climatic conditions.

Evolutionary history of barley cultivation in Europe revealed by genetic analysis of extant landraces

Background

Understanding the evolution of cultivated barley is important for two reasons. First, the evolutionary relationships between different landraces might provide information on the spread and subsequent development of barley cultivation, including the adaptation of the crop to new environments and its response to human selection. Second, evolutionary information would enable landraces with similar traits but different genetic backgrounds to be identified, providing alternative strategies for the introduction of these traits into modern germplasm.

Results

The evolutionary relationships between 651 barley landraces were inferred from the genotypes for 24 microsatellites. The landraces could be divided into nine populations, each with a different geographical distribution. Comparisons with ear row number, caryopsis structure, seasonal growth habit and flowering time revealed a degree of association between population structure and phenotype, and analysis of climate variables indicated that the landraces are adapted, at least to some extent, to their environment. Human selection and/or environmental adaptation may therefore have played a role in the origin and/or maintenance of one or more of the barley landrace populations. There was also evidence that at least some of the population structure derived from geographical partitioning set up during the initial spread of barley cultivation into Europe, or reflected the later introduction of novel varieties. In particular, three closely-related populations were made up almost entirely of plants with the daylength nonresponsive version of the photoperiod response gene PPD-H1, conferring adaptation to the long annual growth season of northern Europe. These three populations probably originated in the eastern Fertile Crescent and entered Europe after the initial spread of agriculture.

Conclusions

The discovery of population structure, combined with knowledge of associated phenotypes and environmental adaptations, enables a rational approach to identification of landraces that might be used as sources of germplasm for breeding programs. The population structure also enables hypotheses concerning the prehistoric spread and development of agriculture to be addressed.

Characterization and comparison of gene-based simple sequence repeats across Brassica species

Simple sequence repeats (SSRs) are important components of eukaryotic genomes and may play an important role in regulating gene expression. However, the characteristics of genic SSRs and the effect of interspecific hybridization and polyploidization on genic SSRs seem not to have received desired attention in terms of scientific investigations. To determine the features of genic SSRs and elucidate their role in polyploidization process of the Brassica family, we identified SSRs in Plant Genome Database-assembled unique transcripts (PUTs) of Brassica species. A higher density of SSRs and a greater number of compound motif SSRs and mononucleotide motif types with large average number of repeats were detected in allotetraploid Brassica napus than in the diploid parental species (Brassica rapa and Brassica oleracea). In addition, a greater proportion of SSR-PUTs were found to be associated with the stress response and developmental processes in B. napus than in the parents. A negative correlation between the repeat number and the motif type and the total length, and a positive correlation between the repeat number and the total length of SSRs were observed. PUT-SSR might be generated from A/T-rich regions. The successful development of 123 pairs of SSR primers for Brassica PUTs showed that SSR-PUTs could be exploited as gene-based SSR functional markers for application in Brassica breeding. These results indicate that interspecific hybridization and polyploidization could trigger the amplification of SSRs, and long SSRs might become shorter to enable the plant to adapt to environmental and artificial selection.

Genetic structure and ecogeographical adaptation in wild barley (Hordeum chilense Roemer et Schultes) as revealed by microsatellite markers

BACKGROUND

Multi-allelic microsatellite markers have become the markers of choice for the determination of genetic structure in plants. Synteny across cereals has allowed the cross-species and cross-genera transferability of SSR markers, which constitute a valuable and cost-effective tool for the genetic analysis and marker-assisted introgression of wild related species. Hordeum chilense is one of the wild relatives with a high potential for cereal breeding, due to its high crossability (both interspecies and intergenera) and polymorphism for adaptation traits. In order to analyze the genetic structure and ecogeographical adaptation of this wild species, it is necessary to increase the number of polymorphic markers currently available for the species. In this work, the possibility of using syntenic wheat SSRs as a new source of markers for this purpose has been explored.

RESULTS

From the 98 wheat EST-SSR markers tested for transferability and polymorphism in the wild barley genome, 53 primer pairs (54.0%) gave cross-species transferability and 20 primer pairs (20.4%) showed polymorphism. The latter were used for further analysis in the H. chilense germplasm. The H. chilense-Triticum aestivum addition lines were used to test the chromosomal location of the new polymorphic microsatellite markers. The genetic structure and diversity was investigated in a collection of 94 H. chilense accessions, using a set of 49 SSR markers distributed across the seven chromosomes. Microsatellite markers showed a total of 351 alleles over all loci. The number of alleles per locus ranged from two to 27, with a mean of 7.2 alleles per locus and a mean Polymorphic Information Content (PIC) of 0.5.

CONCLUSIONS

According to the results, the germplasm can be divided into two groups, with morphological and ecophysiological characteristics being key determinants of the population structure. Geographic and ecological structuring was also revealed in the analyzed germplasm. A significant correlation between geographical and genetic distance was detected in the Central Chilean region for the first time in the species. In addition, significant ecological influence in genetic distance has been detected for one of the population structure groups (group II) in the Central Chilean region. Finally, the association of the SSR markers with ecogeographical variables was investigated and one marker was found significantly associated with precipitation. These findings have a potential application in cereal breeding.

Epistasis in natural populations of a predominantly selfing plant

Populations of predominantly selfing plant species often show spatial genetic structure but little is known whether epistatic gene interactions are spatially structured. To detect a possible epistatic effect and a spatial scale at which it operates, we created artificial crosses between plants spanning a range of fixed distances from 1 to 400 m in three populations of wild barley. The self-pollinated and crossed progeny (F(1)) and two generations of segregated progeny (F(2) and F(3)) were tested in experimentally simulated population environments for relative performance (RP). The measured fitness traits included number of seeds, total seed weight and seed germination. For any of these traits, there was no association between RP of F(1), F(2) and F(3) plants and either pairwise kinship coefficients or crossing distance. In contrast, in all three populations, we found lower seed viability of outcrossed as compared with self-pollinated genotypes in the first generation of segregation. However, in the F(3) generation this outbreeding effect disappeared in the two populations and greatly decreased in the third population. For seed production, heterosis in F(1) and outbreeding depression in F(2) were observed only in the population with unusually high number of heterozygotes. Our findings support the view that in selfing species a spatial mosaic of various locally abundant genotypes represents not randomly fixed combinations of alleles but the co-adapted gene complexes that were sieved by selection, while heterozygotes are characteristic for the transient phase of this process, when segregation and purging of maladaptive genotypes have not yet occurred.

Strong correlation of wild barley (Hordeum spontaneum) population structure with temperature and precipitation variation

In this study, we present the genetic analysis of a new collection of wild barley (Hordeum spontaneum) using 42 simple sequence repeat (SSR) markers that represent the seven chromosomes. The Barley1K (B1K) infrastructure consists of 1020 accessions collected in a hierarchical sampling mode (HSM) from 51 sites across Israel and represents the wide adaptive niche of the modern barley's ancestor. According to the genetic structure analysis, the sampled sites can be divided into seven groups, and sampled microsites located on opposing slopes or in different soil types did not show significant genetic differentiation. Although the genetic analysis indicates a simple isolation-by-distance model among the populations, examination of the genetic populations' structure with abiotic parameters in an ordination analysis revealed that the combination of elevation, mid-day temperature and rainfall explains a high proportion of the variance in the principal components analysis. Our findings demonstrate that the current populations have therefore been shaped and distinguished by non-selective forces such as migration; however, we suggest that aridity and temperature gradients played major roles as selective forces in the adaptation of wild barley in this part of the Fertile Crescent. This unique collection is a prelude for the investigation of the molecular basis underlying plant adaptation and responsiveness to harsh environments.

Adaptive microclimatic structural and expressional dehydrin 1 evolution in wild barley, Hordeum spontaneum, at 'Evolution Canyon', Mount Carmel, Israel

'Evolution Canyon' (ECI) at Lower Nahal Oren, Mount Carmel, Israel, is an optimal natural microscale model for unravelling evolution in action highlighting the twin evolutionary processes of adaptation and speciation. A major model organism in ECI is wild barley, Hordeum spontaneum, the progenitor of cultivated barley, which displays dramatic interslope adaptive and speciational divergence on the 'African' dry slope (AS) and the 'European' humid slope (ES), separated on average by 200 m. Here we examined interslope single nucleotide polymorphism (SNP) sequences and the expression diversity of the drought resistant dehydrin 1 gene (Dhn1) between the opposite slopes. We analysed 47 plants (genotypes), 4-10 individuals in each of seven stations (populations) in an area of 7000 m(2), for Dhn1 sequence diversity located in the 5' upstream flanking region of the gene. We found significant levels of Dhn1 genic diversity represented by 29 haplotypes, derived from 45 SNPs in a total of 708 bp sites. Most of the haplotypes, 25 out of 29 (= 86.2%), were represented by one genotype; hence, unique to one population. Only a single haplotype was common to both slopes. Genetic divergence of sequence and haplotype diversity was generally and significantly different among the populations and slopes. Nucleotide diversity was higher on the AS, whereas haplotype diversity was higher on the ES. Interslope divergence was significantly higher than intraslope divergence. The applied Tajima D rejected neutrality of the SNP diversity. The Dhn1 expression under dehydration indicated interslope divergent expression between AS and ES genotypes, reinforcing Dhn1 associated with drought resistance of wild barley at 'Evolution Canyon'. These results are inexplicable by mutation, gene flow, or chance effects, and support adaptive natural microclimatic selection as the major evolutionary divergent driving force.

Genetic diversity and stress of Ricotia lunaria in "Evolution Canyon," Israel

We examined the genetic diversity and divergence of Ricotia lunaria, a family relative species of Arabidopsis thaliana, sampled from 6 stations on 2 opposing slopes, the south-facing slope ("African" or AS) and north-facing slope ("European" or ES), separated on average by 200 m, at "Evolution Canyon," Lower Nahal Oren, Mount Carmel, Israel, along a transect presenting sharply differing microclimates. The density of R. lunaria populations was slope specific: a higher density and smaller plants were observed on the AS. In addition, the density was positively correlated with annual plant cover. The interslope and intraslope genetic diversities of R. lunaria populations were examined using the amplified fragment length polymorphism (AFLP) technique with 5 primer pairs. Ricotia lunaria populations inhabiting the ES and AS differed, and among the 468 scored loci, 304 (65%) were polymorphic (at P >or= 0.05 level). Polymorphism values obtained for AS and ES populations were similar (52% vs. 56%), but different loci were polymorphic in different populations; 40% of polymorphic loci were identical on both the ES and AS, 16% were polymorphic for the ES only, and 12% were polymorphic only for the AS. The AFLP results grouped the analyzed genotypes into 2 distinct clusters: one cluster included the plants belonging to the AS and the other included ES plants. The unbiased estimate of Nei genetic distances (D) indicated significantly higher interslope (D = 0.124 +/- 0.011) than intraslope (D = 0.076 +/- 0.015) differences (P < 0.001 in t-test). Correspondingly, mean intraslope gene flow was significantly higher than the interslope gene flow (2.9 +/- 0.6 vs. 1.9 +/- 0.2). Natural selection appears to adaptively diverge the plant ecotypes on the opposite slope, both phenotypically and genotypically. This includes significant divergence in flowering time likely to initiate incipient sympatric speciation.

Ecological genomics of natural plant populations: the Israeli perspective

The genomic era revolutionized evolutionary population biology. The ecological genomics of the wild progenitors of wheat and barley reviewed here was central in the research program of the Institute of Evolution, University of Haifa, since 1975 ( http://evolution.haifa.ac.il ). We explored the following questions: (1) How much of the genomic and phenomic diversity of wild progenitors of cultivars (wild emmer wheat, Triticum dicoccoides, the progenitor of most wheat, plus wild relatives of the Aegilops species; wild barley, Hordeum spontaneum, the progenitor of cultivated barley; wild oat, Avena sterilis, the progenitor of cultivated oats; and wild lettuce species, Lactuca, the progenitor and relatives of cultivated lettuce) are adaptive and processed by natural selection at both coding and noncoding genomic regions? (2) What is the origin and evolution of genomic adaptation and speciation processes and their regulation by mutation, recombination, and transposons under spatiotemporal variables and stressful macrogeographic and microgeographic environments? (3) How much genetic resources are harbored in the wild progenitors for crop improvement? We advanced ecological genetics into ecological genomics and analyzed (regionally across Israel and the entire Near East Fertile Crescent and locally at microsites, focusing on the "Evolution Canyon" model) hundreds of populations and thousands of genotypes for protein (allozyme) and deoxyribonucleic acid (DNA) (coding and noncoding) diversity, partly combined with phenotypic diversity. The environmental stresses analyzed included abiotic (climatic and microclimatic, edaphic) and biotic (pathogens, demographic) stresses. Recently, we introduced genetic maps, cloning, and transformation of candidate genes. Our results indicate abundant genotypic and phenotypic diversity in natural plant populations. The organization and evolution of molecular and organismal diversity in plant populations, at all genomic regions and geographical scales, are nonrandom and are positively correlated with, and partly predictable by, abiotic and biotic environmental heterogeneity and stress. Biodiversity evolution, even in small isolated populations, is primarily driven by natural selection including diversifying, balancing, cyclical, and purifying selection regimes interacting with, but, ultimately, overriding the effects of mutation, migration, and stochasticity. The progenitors of cultivated plants harbor rich genetic resources and are the best hope for crop improvement by both classical and modern biotechnological methods. Future studies should focus on the interplay between structural and functional genome organization focusing on gene regulation.

Patterns of genetic and eco-geographical diversity in Spanish barleys

The pool of Western Mediterranean landraces has been under-utilised for barley breeding so far. The objectives of this study were to assess genetic diversity in a core collection of inbred lines derived from Spanish barley landraces to establish its relationship to barleys from other origins, and to correlate the distribution of diversity with geographical and climatic factors. To this end, 64 SSR were used to evaluate the polymorphism among 225 barley (Hordeum vulgare ssp. vulgare) genotypes, comprising two-row and six-row types. These included 159 landraces from the Spanish barley core collection (SBCC) plus 66 cultivars, mainly from European countries, as a reference set. Out of the 669 alleles generated, a large proportion of them were unique to the six-row Spanish barleys. An analysis of molecular variance revealed a clear genetic divergence between the six-row Spanish barleys and the reference cultivars, whereas this was not evident for the two-row barleys. A model-based clustering analysis identified an underlying population structure, consisting of four main populations for the whole genotype set, and suggested further possible subdivision within two of these populations. Most of the six-row Spanish landraces clustered into two groups that corresponded to geographic regions with contrasting environmental conditions. The existence of wide genetic diversity in Spanish germplasm, possibly related to adaptation to a broad range of environmental conditions, and its divergence from current European cultivars confirm its potential as a new resource for barley breeders, and make the SBCC a valuable tool for the study of adaptation in barley.

Patterns of thermal adaptation of Bacillus simplex to the microclimatically contrasting slopes of 'Evolution Canyons' I and II, Israel

Identification of selective forces that drive evolution and speciation of bacteria in natural habitats is a central issue in bacterial ecology and evolution. Exploring the adaptive evolution of Bacillus simplex at 'Evolution Canyons' I and II, Israel, we report here on the impact of high heat stress on the speciation progress of individual evolutionary lineages. These canyons represent similar ecological replicates, separated by 40 km, in which the orientation of the sun yields a strong sun-exposed and hot 'African' south-facing slope (SFS) versus a rather cooler and mesic-lush 'European' north-facing slope (NFS) within a distance of only 50-100 m at the bottom and 400 m at the top. Among 131 strains studied, in Luria-Bertani broth, 'African' strains grow better than 'European' strains at a stressful high temperature (43.25 degrees C). The results suggest that adaptation to the hotter and more stressful SFS is continuously ongoing. The patterns of heat adaptation override the phylogenetic history of individual lineages. A positive correlation of growth rates at 43.25 degrees C and 20 degrees C, more markedly among 'African' strains, reflects probably the broader temperature range on the SFS. Summarizing, the hot temperature stress on the 'African' slope is a major environmental force driving the twin evolutionary processes of adaptation and speciation of B. simplex at 'Evolution Canyon'. Finally, we discuss the data in light of current controversies on species concepts.

The rise, fall and renaissance of microsatellites in eukaryotic genomes

Microsatellites are among the most versatile of genetic markers, being used in an impressive number of biological applications. However, the evolutionary dynamics of these markers remain a source of contention. Almost 20 years after the discovery of these ubiquitous simple sequences, new genomic data are clarifying our understanding of the structure, distribution and variability of microsatellites in genomes, especially for the eukaryotes. While these new data provide a great deal of descriptive information about the nature and abundance of microsatellite sequences within eukaryotic genomes, there have been few attempts to synthesise this information to develop a global concept of evolution. This review provides an up-to-date account of the mutational processes, biases and constraints believed to be involved in the evolution of microsatellites, particularly with respect to the creation and degeneration of microsatellites, which we assert may be broadly viewed as a life cycle. In addition, we identify areas of contention that require further research and propose some possible directions for future investigation.

On The Necessity to Study Natural Bacterial Populations-The Model of Bacillus Simplex From "Evolution Canyons" I and II, Israel

How do bacteria evolve and speciate in natural environments? How does bacterial evolution relate to bacterial systematics? Exploring these answers is essential because bacteria profoundly impact life in general and, in particular, that of humans. Much insight into bacterial microevolution has come from theoretical and computational studies and from multigenerational laboratory systems ("Experimental Evolution"). These studies, however, do not take into account the diversity of modes of how bacteria can evolve under the complexity of the real world, i.e., nature. We argue, therefore, that for a comprehensive understanding of bacterial microevolution, it is essential to study natural populations. We underline our argument by introducing theBacillus simplexmodel from "Evolution Canyon", Israel. This metapopulation splits into different evolutionary lineages that have adapted to the microclimatically different slopes of "Evolution Canyon". It was shown that temperature stress is a major environmental factor driving theB. simplexadaptation and speciation progress. Therefore, this model population has proven highly suitable to study bacterial microevolution in natural habitats. Finally, we discuss theB. simplexintrapopulation divergence of lineages in light of current controversies on bacterial species concepts and taxon identification.

"Evolution Canyon": A Microcosm of Life's Evolution Focusing on Adaptation and Speciation

Local microcosmic natural laboratories, dubbed "Evolution Canyon" (EC) models, reinforce studies of regional and global macrocosmic ecological theaters across life and unravelevolution in action.The EC model laboratories permit genomic, proteomic, and phenomic studies highlighting speciation and adaptation at a microscale. Critical transplant experiment tests can evaluate interslope differential fitness. Novel techniques of genetic mapping, sequence nucleotide polymorphisms (SNPs), and wide genome coding and noncoding expressions can unravel evolutionary dynamics. Finally, fundamental problems such as stress effects on nonrandom mutations, lateral transfers, splicing variations, sex, and social evolutions, and adaptive strategies of prokaryotes and eukaryotes are testable. We are studying four "Evolution Canyons" (EC I-IV) in the Carmel, Galilee, Negev, and Golan mountains. We've identified 2,500 species in EC I (Carmel) from bacteria to mammals in an area of 7,000 m. Higher terrestrial species richness was found on the more stressful tropical "African" slope (AS). Aquatic species richness was higher on the milder, temperate "European" slope (ES). In soil fungi we found interslope local and regional adaptive divergence in sex, melanism, and conidia. In nine out of 14 (64%) model organisms across life we identified largely higher genetic polymorphisms on the more stressful "African" slope. In some model species we found higher mutation rates, gene conversion, recombination, DNA repair, and larger genome size as well as interslope divergent micro-satellites, molecular polymorphisms, retrotransposons, and wide-genome gene expression on the more stressful AS. Remarkably, interslope incipient sympatric speciation was found across life. The "Evolution Canyon" model represents the Israeli ecological equivalent of the Galapagos Islands. Micro-climatic selection overrides drift and drives both interslope adaptive radiation and incipient sympatric speciation. The EC model could potentially highlight many mysteries of evolutionary biology.

Has Simple Sequence Repeat Mutability Been Selected to Facilitate Evolution?

While adaptation and speciation begin with heritable variation, the underlying processes of mutation remain poorly understood. One particularly interesting source for prolific and adaptively meaningful variation is presented by the exceptionally high mutability of simple sequence repeats (SSRs: microsatellites and minisatellites). Frequent mutations at SSR sites alter the number of tandem repeats and create extensive polymorphism. Although most SSR variants are commonly presumed to be neutral, SSR variation has been shown to influence many biochemical, morphological, physiological, and behavioral characters, with at least a few examples offering evidence of response to selection. The type and degree of phenotypic variation depend upon each SSR's motif and on its location in exon, intron, or regulatory region, but the generation of abundant repeat-number variation is intrinsic to all of these repetitive sequences. Given the widespread distribution of SSRs within most genomes and their potential to modify almost any aspect of gene function, we believe that SSR mutability can facilitate evolutionary adaptation. Furthermore, we argue that the properties of SSRs allow natural selection to favor, indirectly, the mutability of these sites, in contrast to a conventional expectation that selection normally minimizes mutation rates by balancing the cost of deleterious mutations against the cost of replication fidelity. We believe that SSR mutability is not an "accident" of DNA replication, but has been adjusted and selected for this role. SSRs thus have a true biological function as general-purpose "tuning knobs" whereby mutations provide reversible adjustment for many quantitative and qualitative traits.