Distinct geographic patterns of genetic diversity are maintained in wild barley (Hordeum vulgare ssp. spontaneum) despite migration

UDP-glucosyltransferase HvUGT13248 confers type II resistance to Fusarium graminearum in barley

Fusarium head blight (FHB) of barley (Hordeum vulgare) causes yield losses and accumulation of trichothecene mycotoxins (e.g. deoxynivalenol [DON]) in grains. Glucosylation of DON to the nontoxic DON-3-O-glucoside (D3G) is catalyzed by UDP-glucosyltransferases (UGTs), such as barley UGT13248. We explored the natural diversity of UGT13248 in 496 barley accessions and showed that all carried potential functional alleles of UGT13248, as no genotypes showed strongly increased seedling sensitivity to DON. From a TILLING population, we identified 2 mutant alleles (T368I and H369Y) that, based on protein modeling, likely affect the UDP-glucose binding of UGT13248. In DON feeding experiments, DON-to-D3G conversion was strongly reduced in spikes of these mutants compared to controls, and plants overexpressing UGT13248 showed increased resistance to DON and increased DON-to-D3G conversion. Moreover, field-grown plants carrying the T368I or H369Y mutations inoculated with Fusarium graminearum showed increased FHB disease severity and reduced D3G production. Barley is generally considered to have type II resistance that limits the spread of F. graminearum from the infected spikelet to adjacent spikelets. Point inoculation experiments with F. graminearum showed increased infection spread in T368I and H369Y across the spike compared to wild type, while overexpression plants showed decreased spread of FHB symptoms. Confocal microscopy revealed that F. graminearum spread to distant rachis nodes in T368I and H369Y mutants but was arrested at the rachis node of the inoculated spikelet in wild-type plants. Taken together, our data reveal that UGT13248 confers type II resistance to FHB in barley via conjugation of DON to D3G.

Genome resequencing and transcriptome profiling reveal molecular evidence of tolerance to water deficit in barley

INTRODUCTION

Frequent climate change-induced drought events are detrimental environmental stresses affecting global crop production and ecosystem health. Several efforts have facilitated crop breeding for resilient varieties to counteract stress. However, progress is hampered due to the complexity of drought tolerance; a greater variety of novel genes are required across varying environments. Tibetan annual wild barley is a unique and precious germplasm that is well adapted to abiotic stress and can provide elite genes for crop improvement in drought tolerance.

OBJECTIVES

To identify the genetic basis and unique mechanisms for drought tolerance in Tibetan wild barley.

METHODS

Whole genome resequencing and comparative RNA-seq approaches were performed to identify candidate genes associated with drought tolerance via investigating the genetic diversity and transcriptional variation between cultivated and Tibetan wild barley. Bioinformatics, population genetics, and gene silencing were conducted to obtain insights into ecological adaptation in barley and functions of key genes.

RESULTS

Over 20 million genetic variants and a total of 15,361 significantly affected genes were identified in our dataset. Combined genomic, transcriptomic, evolutionary, and experimental analyses revealed 26 water deficit resilience-associated genes in the drought-tolerant wild barley XZ5 with unique genetic variants and expression patterns. Functional prediction revealed Tibetan wild barley employs effective regulators to activate various responsive pathways with novel genes, such as Zinc-Induced Facilitator-Like 2 (HvZIFL2) and Peroxidase 11 (HvPOD11), to adapt to water deficit conditions. Gene silencing and drought tolerance evaluation in a natural barley population demonstrated that HvZIFL2 and HvPOD11 positively regulate drought tolerance in barley.

CONCLUSION

Our findings reveal functional genes that have been selected across barley's complex history of domestication to thrive in water deficit environments. This will be useful for molecular breeding and provide new insights into drought-tolerance mechanisms in wild relatives of major cereal crops.

Transcriptome and Metabolite Insights into Domestication Process of Cultivated Barley in China

The domestication process of cultivated barley in China remains under debate because of the controversial origins of barley. Here, we analyzed transcriptomic and non-targeted metabolic data from 29 accessions together with public resequencing data from 124 accessions to explore the domestication process of cultivated barley in China (Cb-C). These analyses revealed that both Cb-C and Tibetan wild barley (Wb-T) were the descendants of wild barley from the Near East Fertile Crescent (Wb-NE), yielding little support for a local origin of Wb-T. Wb-T was more likely an intermediate in the domestication process from Wb-NE to Cb-C. Wb-T contributed more genetically to Cb-C than Wb-NE, and was domesticated into Cb-C about 3300 years ago. These results together seem to support that Wb-T may be a feralized or hybrid form of cultivated barley from the Near East Fertile Crescent or central Asia. Additionally, the metabolite analysis revealed divergent metabolites of alkaloids and phenylpropanoids and these metabolites were specifically targeted for selection in the evolutionary stages from Wb-NE to Wb-T and from Wb-T to Cb-C. The key missense SNPs in the genes HORVU6Hr1G027650 and HORVU4Hr1G072150 might be responsible for the divergence of metabolites of alkaloids and phenylpropanoids during domestication. Our findings allow for a better understanding of the domestication process of cultivated barley in China.

Domestication shapes recombination patterns in tomato

Meiotic recombination is a biological process of key importance in breeding, to generate genetic diversity and develop novel or agronomically relevant haplotypes. In crop tomato, recombination is curtailed as manifested by linkage disequilibrium decay over a longer distance and reduced diversity compared with wild relatives. Here, we compared domesticated and wild populations of tomato and found an overall conserved recombination landscape, with local changes in effective recombination rate in specific genomic regions. We also studied the dynamics of recombination hotspots resulting from domestication and found that loss of such hotspots is associated with selective sweeps, most notably in the pericentromeric heterochromatin. We detected footprints of genetic changes and structural variants, among them associated with transposable elements, linked with hotspot divergence during domestication, likely causing fine-scale alterations to recombination patterns and resulting in linkage drag.

Natural variation in meiotic recombination rate shapes introgression patterns in intraspecific hybrids between wild and domesticated barley

Meiotic recombination rates vary considerably between species, populations and individuals. The genetic exchange between homologous chromosomes plays a major role in evolution by breaking linkage between advantageous and deleterious alleles in the case of introgressions. Identifying recombination rate modifiers is thus of both fundamental and practical interest to understand and utilize variation in meiotic recombination rates. We investigated recombination rate variation in a large intraspecific hybrid population (named HEB-25) derived from a cross between domesticated barley and 25 wild barley accessions. We observed quantitative variation in total crossover number with a maximum of a 1.4-fold difference between subpopulations and increased recombination rates across pericentromeric regions. The meiosis-specific α-kleisin cohesin subunit REC8 was identified as a candidate gene influencing crossover number and patterning. Furthermore, we quantified wild barley introgression patterns and revealed how local and genome-wide recombination rate variation shapes patterns of introgression. The identification of allelic variation in REC8 in combination with the observed changes in crossover patterning suggest a difference in how chromatin loops are tethered to the chromosome axis, resulting in reduced crossover suppression across pericentromeric regions. Local and genome-wide recombination rate variation is shaping patterns of introgressions and thereby directly influences the consequences of linkage drag.

Adaptive signals of flowering time pathways in wild barley from Israel over 28 generations

Flowering time is one of the most critical traits for plants' life cycles, which is influenced by various environment changes, such as global warming. Previous studies have suggested that to guarantee reproductive success, plants have shifted flowering times to adapt to global warming. Although many studies focused on the molecular mechanisms of early flowering, little was supported by the repeated sampling at different time points through the changing climate. To fully dissect the temporal and spatial evolutionary genetics of flowering time, we investigated nucleotide variation in ten flowering time candidate genes and nine reference genes for the same ten wild-barley populations sampled 28 years apart (1980-2008). The overall genetic differentiation was significantly greater in the descendant populations (2008) compared with the ancestral populations (1980); however, local adaptation tests failed to detect any single-nucleotide polymorphism (SNP)/indel under spatial-diversifying selection at either time point. By contrast, the WFABC (Wright-Fisher ABC-based approach) that detected 54 SNPs/indels was under strong selection during the past 28 generations. Moreover, all these 54 alleles were segregated in the ancestral populations, but fixed in the descendent populations. Among the top ten SNPs/indels, seven were located in genes of FT1 (FLOWERING TIME LOCUS T 1), CO1 (CONSTANS-LIKE PROTEIN 1), and VRN-H2 (VERNALIZATION-H2), which have been documented to be associated with flowering time regulation in barley cultivars. This study might suggest that all ten populations have undergone parallel evolution over the past few decades in response to global warming, and even an overwhelming local adaptation and ecological differentiation.

Distinct patterns of natural selection in Na+/H+ antiporter genes in Populus euphratica and Populus pruinosa

Salt tolerance genes constitute an important class of loci in plant genomes. Little is known about the extent to which natural selection in saline environments has acted upon these loci, and what types of nucleotide diversity such selection has given rise to. Here, we surveyed genetic diversity in three types of Na⁺/H⁺ antiporter gene (SOS, NhaD, and NHX, belonging to the cation/proton antiporter 1 family), which have well‐characterized essential roles in plant salt tolerance. Ten Na⁺/H⁺ antiporter genes and 16 neutral loci randomly selected as controls were sequenced from 17 accessions of two closely related members of the genus Populus, Populus euphratica and Populus pruinosa, section Turanga, which are native to northwest China. The results show that salt tolerance genes are common targets of natural selection in P. euphratica and P. pruinosa. Moreover, the patterns of nucleotide variation across the three types of Na⁺/H⁺ antiporter gene are distinctly different in these two closely related Populus species, and gene flow from P. pruinosa to P. euphratica is highly restricted. Our results suggest that natural selection played an important role in shaping the current distinct patterns of Na⁺/H⁺ antiporter genes, resulting in adaptive evolution in P. euphratica and P. pruinosa.

Molecular evidence of RNA polymerase II gene reveals the origin of worldwide cultivated barley

The origin and domestication of cultivated barley have long been under debate. A population-based resequencing and phylogenetic analysis of the single copy of RPB2 gene was used to address barley domestication, to explore genetic differentiation of barley populations on the worldwide scale, and to understand gene-pool exchanges during the spread and subsequent development of barley cultivation. Our results revealed significant genetic differentiation among three geographically distinct wild barley populations. Differences in haplotype composition among populations from different geographical regions revealed that modern cultivated barley originated from two major wild barley populations: one from the Near East Fertile Crescent and the other from the Tibetan Plateau, supporting polyphyletic origin of cultivated barley. The results of haplotype frequencies supported multiple domestications coupled with widespread introgression events that generated genetic admixture between divergent barley gene pools. Our results not only provide important insight into the domestication and evolution of cultivated barley, but also enhance our understanding of introgression and distinct selection pressures in different environments on shaping the genetic diversity of worldwide barley populations, thus further facilitating the effective use of the wild barley germplasm.

The Role of Deleterious Substitutions in Crop Genomes

Populations continually incur new mutations with fitness effects ranging from lethal to adaptive. While the distribution of fitness effects of new mutations is not directly observable, many mutations likely either have no effect on organismal fitness or are deleterious. Historically, it has been hypothesized that a population may carry many mildly deleterious variants as segregating variation, which reduces the mean absolute fitness of the population. Recent advances in sequencing technology and sequence conservation-based metrics for inferring the functional effect of a variant permit examination of the persistence of deleterious variants in populations. The issue of segregating deleterious variation is particularly important for crop improvement, because the demographic history of domestication and breeding allows deleterious variants to persist and reach moderate frequency, potentially reducing crop productivity. In this study, we use exome resequencing of 15 barley accessions and genome resequencing of 8 soybean accessions to investigate the prevalence of deleterious single nucleotide polymorphisms (SNPs) in the protein-coding regions of the genomes of two crops. We conclude that individual cultivars carry hundreds of deleterious SNPs on average, and that nonsense variants make up a minority of deleterious SNPs. Our approach identifies known phenotype-altering variants as deleterious more frequently than the genome-wide average, suggesting that putatively deleterious variants are likely to affect phenotypic variation. We also report the implementation of a SNP annotation tool BAD_Mutations that makes use of a likelihood ratio test based on alignment of all currently publicly available Angiosperm genomes.

Exome sequencing of geographically diverse barley landraces and wild relatives gives insights into environmental adaptation

After domestication, during a process of widespread range extension, barley adapted to a broad spectrum of agricultural environments. To explore how the barley genome responded to the environmental challenges it encountered, we sequenced the exomes of a collection of 267 georeferenced landraces and wild accessions. A combination of genome-wide analyses showed that patterns of variation have been strongly shaped by geography and that variant-by-environment associations for individual genes are prominent in our data set. We observed significant correlations of days to heading (flowering) and height with seasonal temperature and dryness variables in common garden experiments, suggesting that these traits were major drivers of environmental adaptation in the sampled germplasm. A detailed analysis of known flowering-associated genes showed that many contain extensive sequence variation and that patterns of single- and multiple-gene haplotypes exhibit strong geographical structuring. This variation appears to have substantially contributed to range-wide ecogeographical adaptation, but many factors key to regional success remain unidentified.

Evolutionary history of wild barley (Hordeum vulgare subsp. spontaneum) analyzed using multilocus sequence data and paleodistribution modeling

Studies of Hordeum vulgare subsp. spontaneum, the wild progenitor of cultivated barley, have mostly relied on materials collected decades ago and maintained since then ex situ in germplasm repositories. We analyzed spatial genetic variation in wild barley populations collected rather recently, exploring sequence variations at seven single-copy nuclear loci, and inferred the relationships among these populations and toward the genepool of the crop. The wild barley collection covers the whole natural distribution area from the Mediterranean to Middle Asia. In contrast to earlier studies, Bayesian assignment analyses revealed three population clusters, in the Levant, Turkey, and east of Turkey, respectively. Genetic diversity was exceptionally high in the Levant, while eastern populations were depleted of private alleles. Species distribution modeling based on climate parameters and extant occurrence points of the taxon inferred suitable habitat conditions during the ice-age, particularly in the Levant and Turkey. Together with the ecologically wide range of habitats, they might contribute to structured but long-term stable populations in this region and their high genetic diversity. For recently collected individuals, Bayesian assignment to geographic clusters was generally unambiguous, but materials from genebanks often showed accessions that were not placed according to their assumed geographic origin or showed traces of introgression from cultivated barley. We assign this to gene flow among accessions during ex situ maintenance. Evolutionary studies based on such materials might therefore result in wrong conclusions regarding the history of the species or the origin and mode of domestication of the crop, depending on the accessions included.

Transcriptome profiling reveals mosaic genomic origins of modern cultivated barley

The domestication of cultivated barley has been used as a model system for studying the origins and early spread of agrarian culture. Our previous results indicated that the Tibetan Plateau and its vicinity is one of the centers of domestication of cultivated barley. Here we reveal multiple origins of domesticated barley using transcriptome profiling of cultivated and wild-barley genotypes. Approximately 48-Gb of clean transcript sequences in 12 Hordeum spontaneum and 9 Hordeum vulgare accessions were generated. We reported 12,530 de novo assembled transcripts in all of the 21 samples. Population structure analysis showed that Tibetan hulless barley (qingke) might have existed in the early stage of domestication. Based on the large number of unique genomic regions showing the similarity between cultivated and wild-barley groups, we propose that the genomic origin of modern cultivated barley is derived from wild-barley genotypes in the Fertile Crescent (mainly in chromosomes 1H, 2H, and 3H) and Tibet (mainly in chromosomes 4H, 5H, 6H, and 7H). This study indicates that the domestication of barley may have occurred over time in geographically distinct regions.

Two genomic regions contribute disproportionately to geographic differentiation in wild barley

Genetic differentiation in natural populations is driven by geographic distance and by ecological or physical features within and between natural habitats that reduce migration. The primary population structure in wild barley differentiates populations east and west of the Zagros Mountains. Genetic differentiation between eastern and western populations is uneven across the genome and is greatest on linkage groups 2H and 5H. Genetic markers in these two regions demonstrate the largest difference in frequency between the primary populations and have the highest informativeness for assignment to each population. Previous cytological and genetic studies suggest there are chromosomal structural rearrangements (inversions or translocations) in these genomic regions. Environmental association analyses identified an association with both temperature and precipitation variables on 2H and with precipitation variables on 5H.

Two genomic regions contribute disproportionately to geographic differentiation in wild barley

Genetic differentiation in natural populations is driven by geographic distance and by ecological or physical features within and between natural habitats that reduce migration. The primary population structure in wild barley differentiates populations east and west of the Zagros Mountains. Genetic differentiation between eastern and western populations is uneven across the genome and is greatest on linkage groups 2H and 5H. Genetic markers in these two regions demonstrate the largest difference in frequency between the primary populations and have the highest informativeness for assignment to each population. Previous cytological and genetic studies suggest there are chromosomal structural rearrangements (inversions or translocations) in these genomic regions. Environmental association analyses identified an association with both temperature and precipitation variables on 2H and with precipitation variables on 5H.

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.

Resequencing data indicate a modest effect of domestication on diversity in barley: a cultigen with multiple origins

The levels of diversity and extent of linkage disequilibrium in cultivated species are largely determined by diversity in their wild progenitors. We report a comparison of nucleotide sequence diversity in wild and cultivated barley (Hordeum vulgare ssp. spontaneum and ssp. vulgare) at 7 nuclear loci totaling 9296bp, using sequence from Hordeum bulbosum to infer the ancestral state of mutations. The sample includes 36 accessions of cultivated barley, including 23 landraces (cultivated forms not subject to modern breeding) and 13 cultivated lines and genetic stocks compared to either 25 or 45 accessions of wild barley for the same loci. Estimates of nucleotide sequence diversity indicate that landraces retain >80% of the diversity in wild barley. The primary population structure in wild barley, which divides the species into eastern and western populations, is reflected in significant differentiation at all loci in wild accessions and at 3 of 7 loci in landraces. "Oriental" landraces have slightly higher diversity than "Occidental" landraces. Genetic assignment suggests more admixture from Occidental landraces into Oriental landraces than the converse, which may explain this difference. Based on θπ for silent sites, modern western cultivars have ~73% of the diversity found in landraces and ~71% of the diversity in wild barley.

Demographic factors shaped diversity in the two gene pools of wild common bean Phaseolus vulgaris L

Wild common bean (Phaseolus vulgaris L.) is distributed throughout the Americas from Mexico to northern Argentina. Within this range, the species is divided into two gene pools (Andean and Middle American) along a latitudinal gradient. The diversity of 24 wild common bean genotypes from throughout the geographic range of the species was described by using sequence data from 13 loci. An isolation-migration model was evaluated using a coalescent analysis to estimate multiple demographic parameters. Using a Bayesian approach, Andean and Middle American subpopulations with high percentage of parentages were observed. Over all loci, the Middle American gene pool was more diverse than the Andean gene pool (π(sil)=0.0089 vs 0.0068). The two subpopulations were strongly genetically differentiated over all loci (F(st)=0.29). It is estimated that the two current wild gene pools diverged from a common ancestor ∼111 000 years ago. Subsequently, each gene pool underwent a bottleneck immediately after divergence and lasted ∼40 000 years. The Middle American bottleneck population size was ∼46% of the ancestral population size, whereas the Andean was 26%. Continuous asymmetric gene flow was detected between the two gene pools with a larger number of migrants entering Middle American gene pool from the Andean gene pool. These results suggest that because of the complex population structure associated with the ancestral divergence, subsequent bottlenecks in each gene pool, gene pool-specific domestication and intense selection within each gene pool by breeders; association mapping would best be practised within each common bean gene pool.

A critical review of the protracted domestication model for Near-Eastern founder crops: linear regression, long-distance gene flow, archaeological, and archaeobotanical evidence

The recent review by Fuller et al. (2012a) in this journal is part of a series of papers maintaining that plant domestication in the Near East was a slow process lasting circa 4000 years and occurring independently in different locations across the Fertile Crescent. Their protracted domestication scenario is based entirely on linear regression derived from the percentage of domesticated plant remains at specific archaeological sites and the age of these sites themselves. This paper discusses why estimates like haldanes and darwins cannot be applied to the seven founder crops in the Near East (einkorn and emmer wheat, barley, peas, chickpeas, lentils, and bitter vetch). All of these crops are self-fertilizing plants and for this reason they do not fulfil the requirements for performing calculations of this kind. In addition, the percentage of domesticates at any site may be the result of factors other than those that affect the selection for domesticates growing in the surrounding area. These factors are unlikely to have been similar across prehistoric sites of habitation, societies, and millennia. The conclusion here is that single crop analyses are necessary rather than general reviews drawing on regression analyses based on erroneous assumptions. The fact that all seven of these founder crops are self-fertilizers should be incorporated into a comprehensive domestication scenario for the Near East, as self-fertilization naturally isolates domesticates from their wild progenitors.

Adaptive genetic differentiation in a predominantly self-pollinating species analyzed by transplanting into natural environment, crossbreeding and Q(ST)-F(ST) test

• Both genetic drift and natural selection result in genetic/phenotypic differentiation over space. I analyzed the role of local adaptation in the genetic differentiation of populations of the annual grass Hordeum spontaneum sampled along an aridity gradient. • The study included the introduction of plants having desert vs nondesert origin into natural (desert) environment, analysis of population differentiation in allozymes and random amplified polymorphic DNA (RAPD) markers vs phenotypic traits (Q(ST) -F(ST) comparison), and planting interpopulation hybrids under simulated desert conditions in a glasshouse. • The results of the home advantage test, Q(ST) -F(ST) comparison and crossbreeding were consistent with local adaptation; that is, that differentiation of the desert plants from plants of nondesert origin in phenotypic traits was adaptive, giving them home advantage. Each method used provided additional, otherwise unavailable, information, meaning that they should be viewed as complementary rather than alternative approaches. • Gene flow from adjacent populations (i.e. populations experiencing the desert environment) via seeds (but not pollen) had a positive effect on fitness by enhancing natural selection and counteracting drift. At the same time, the effect of genes from the species distributional core (nondesert plants) by either seed or pollen had a negative fitness effect despite its enriching effect on neutral diversity. The pattern of outbreeding depression observed in interpopulation hybrids (F(1) ) and their segregating progeny (F(2) ) was inconsistent with underdominance, but indicated the presence of additive, dominance and epistatic effects.

On the genome constitution and evolution of intermediate wheatgrass (Thinopyrum intermedium: Poaceae, Triticeae)

Background

The wheat tribe Triticeae (Poaceae) is a diverse group of grasses representing a textbook example of reticulate evolution. Apart from globally important grain crops, there are also wild grasses which are of great practical value. Allohexaploid intermediate wheatgrass, Thinopyrum intermedium (2n = 6x = 42), possesses many desirable agronomic traits that make it an invaluable source of genetic material useful in wheat improvement. Although the identification of its genomic components has been the object of considerable investigation, the complete genomic constitution and its potential variability are still being unravelled. To identify the genomic constitution of this allohexaploid, four accessions of intermediate wheatgrass from its native area were analysed by sequencing of chloroplast trnL-F and partial nuclear GBSSI, and genomic in situ hybridization.

Results

The results confirmed the allopolyploid origin of Thinopyrum intermedium and revealed new aspects in its genomic composition. Genomic heterogeneity suggests a more complex origin of the species than would be expected if it originated through allohexaploidy alone. While Pseudoroegneria is the most probable maternal parent of the accessions analysed, nuclear GBSSI sequences suggested the contribution of distinct lineages corresponding to the following present-day genera: Pseudoroegneria, Dasypyrum, Taeniatherum, Aegilops and Thinopyrum. Two subgenomes of the hexaploid have most probably been contributed by Pseudoroegneria and Dasypyrum, but the identity of the third subgenome remains unresolved satisfactorily. Possibly it is of hybridogenous origin, with contributions from Thinopyrum and Aegilops. Surprising diversity of GBSSI copies corresponding to a Dasypyrum-like progenitor indicates either multiple contributions from different sources close to Dasypyrum and maintenance of divergent copies or the presence of divergent paralogs, or a combination of both. Taeniatherum-like GBSSI copies are most probably pseudogenic, and the mode of their acquisition by Th. intermedium remains unclear.

Conclusions

Hybridization has played a key role in the evolution of the Triticeae. Transfer of genetic material via extensive interspecific hybridization and/or introgression could have enriched the species' gene pools significantly. We have shown that the genomic heterogeneity of intermediate wheatgrass is higher than has been previously assumed, which is of particular concern to wheat breeders, who frequently use it as a source of desirable traits in wheat improvement.

Population structure in sorghum accessions from West Africa differing in race and maturity class

Accounting for population structure to minimize spurious associations in association analyses is of crucial importance. With sorghum genomic sequence information being available, there is a growing interest in performing such association studies for a number of important agronomic traits using a candidate gene approach. The aims of our study were to conduct a systematic survey of molecular genetic diversity and analyze the population structure in cultivated sorghum [Sorghum bicolor (L.) Moench] accessions from West Africa. Our analysis included 219 West African cultivated sorghum accessions with differing maturity intended for a marker-trait association study. A total of 27 SSRs were used, which resulted in detection of 513 alleles. Genetic diversity estimates for the accessions were found to be high. The accessions were divided into two subgroups using a model-based approach. Our findings partly agree with previous studies in that the guinea race accessions could be distinguished clearly from other accessions included in the analysis. Race and geographical origin of the accessions may be responsible for the structure we observed in our material. The extent of linkage disequilibrium for all combinations of SSRs was in agreement with expectations based on the mating system.

Population genetics of genomics-based crop improvement methods

Many genome-wide association studies (GWAS) in humans are concluding that, even with very large sample sizes and high marker densities, most of the genetic basis of complex traits may remain unexplained. At the same time, recent research in plant GWAS is showing much greater success with fewer resources. Both GWAS and genomic selection (GS), a method for predicting phenotypes by the use of genome-wide marker data, are receiving considerable attention among plant breeders. In this review we explore how differences in population genetic histories, as well as past selection for traits of interest, have produced trait architectures and patterns of linkage disequilibrium (LD) that frequently differ dramatically between domesticated plants and humans, making detection of quantitative trait loci (QTL) effects in crops more rewarding and less costly than in humans.

Genetic analysis and ecological association of Hina genes based on single nucleotide polymorphisms (SNPs) in wild barley, Hordeum spontaneum

Specific primers were designed to amplify the sequences of Hina genes from 121 wild barley (Hordeum spontaneum) accessions belonging to 18 populations from Iran, Israel and Turkey. Forty-nine single nucleotide polymorphisms (SNPs), nine indels, and 26 haplotypes were determined by sequence analysis. The genetic polymorphism (P), genetic diversity (He), and Shannon's information index (I) in the 18 populations were 0.486, 0.181 and 0.269, respectively. Approximately 2/3 genetic variations of Hina genes were presented within populations, while approximately 1/3 genetic variations were observed between populations. Broad gene flow (Nm= 3.31) and low genetic variation (Gst= 0.0702) were detected. However, the genetic differentiation between populations was independent of geographical distances according to the Mantel test (p = 0.478). The result of Spearman rank correlations (r(s)) showed that the genetic indices (P, He and I) of Hina were not significantly correlated with ecological factors. Only eight SNP positions correlated significantly with ecological factors. Of the eight SNP positions that positively correlated with ecological factors, only one SNP (769, T-C) was located in the coding region; however, it was not responsible for the amino acid change.

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.

Development of genome-specific primers for homoeologous genes in allopolyploid species: the waxy and starch synthase II genes in allohexaploid wheat (Triticum aestivum L.) as examples

BACKGROUND

In allopolypoid crops, homoeologous genes in different genomes exhibit a very high sequence similarity, especially in the coding regions of genes. This makes it difficult to design genome-specific primers to amplify individual genes from different genomes. Development of genome-specific primers for agronomically important genes in allopolypoid crops is very important and useful not only for the study of sequence diversity and association mapping of genes in natural populations, but also for the development of gene-based functional markers for marker-assisted breeding. Here we report on a useful approach for the development of genome-specific primers in allohexaploid wheat.

FINDINGS

In the present study, three genome-specific primer sets for the waxy (Wx) genes and four genome-specific primer sets for the starch synthase II (SSII) genes were developed mainly from single nucleotide polymorphisms (SNPs) and/or insertions or deletions (Indels) in introns and intron-exon junctions. The size of a single PCR product ranged from 750 bp to 1657 bp. The total length of amplified PCR products by these genome-specific primer sets accounted for 72.6%-87.0% of the Wx genes and 59.5%-61.6% of the SSII genes. Five genome-specific primer sets for the Wx genes (one for Wx-7A, three for Wx-4A and one for Wx-7D) could distinguish the wild type wheat and partial waxy wheat lines. These genome-specific primer sets for the Wx and SSII genes produced amplifications in hexaploid wheat, cultivated durum wheat, and Aegilops tauschii accessions, but failed to generate amplification in the majority of wild diploid and tetraploid accessions.

CONCLUSIONS

For the first time, we report on the development of genome-specific primers from three homoeologous Wx and SSII genes covering the majority of the genes in allohexaploid wheat. These genome-specific primers are being used for the study of sequence diversity and association mapping of the three homoeologous Wx and SSII genes in natural populations of both hexaploid wheat and cultivated tetraploid wheat. The strategies used in this paper can be used to develop genome-specific primers for homoeologous genes in any allopolypoid species. They may be also suitable for (i) the development of gene-specific primers for duplicated paralogous genes in any diploid species, and (ii) the development of allele-specific primers at the same gene locus.

Fine-scale spatial genetic structure in a predominantly selfing plant: role of seed and pollen dispersal

We present a study of fine-scale spatial genetic structure (SGS) and assess the impact of seed and pollen dispersal on the pattern of genetic diversity in the predominantly selfing Hordeum spontaneum. The study included (1) direct measurement of dispersal in a controlled environment, and (2) analyses of SGS and estimation of the ratio of pollen to seed flow in three natural populations sampled in linear transects at fixed increasing inter-plant distances. Analysis of SGS with 10 nuclear SSRs showed in all three populations a significant autocorrelation for the distance classes of 1 or 2 m and a negative linear relationship between kinship coefficients, calculated for pairs of individuals, and logarithm of geographical distance between members of the pairs. Major seed dispersal (95%) was found to be within 1.2 m from the mother plant. Pollen flow, estimated from the comparison of nuclear and chloroplast variation, was spatially limited as much as was seed dispersal, and tended to be overestimated when measured at spatial scales exceeding that of SGS. We conclude that combined effects of selfing, occasional outcrossing, localized seed dispersal and high plant density create an equilibrium between drift and gene flow in this species resulting in SGS at a very fine spatial scale.

Patterns of genetic diversity and linkage disequilibrium in a highly structured Hordeum vulgare association-mapping population for the Mediterranean basin

Population structure and genome-wide linkage disequilibrium (LD) were investigated in 192 Hordeum vulgare accessions providing a comprehensive coverage of past and present barley breeding in the Mediterranean basin, using 50 nuclear microsatellite and 1,130 DArT((R)) markers. Both clustering and principal coordinate analyses clearly sub-divided the sample into five distinct groups centred on key ancestors and regions of origin of the germplasm. For given genetic distances, large variation in LD values was observed, ranging from closely linked markers completely at equilibrium to marker pairs at 50 cM separation still showing significant LD. Mean LD values across the whole population sample decayed below r (2) of 0.15 after 3.2 cM. By assaying 1,130 genome-wide DArT((R)) markers, we demonstrated that, after accounting for population substructure, current genome coverage of 1 marker per 1.5 cM except for chromosome 4H with 1 marker per 3.62 cM is sufficient for whole genome association scans. We show, by identifying associations with powdery mildew that map in genomic regions known to have resistance loci, that associations can be detected in strongly stratified samples provided population structure is effectively controlled in the analysis. The population we describe is, therefore, shown to be a valuable resource, which can be used in basic and applied research in barley.

Molecular population genetics and agronomic alleles in seed banks: searching for a needle in a haystack?

Seed banking has been the single most significant reaction of the research community to the alarming rates of plant genetic erosion occurring in the wild. One enduring challenge for a wiser utilization of the resources enclosed in seed banks, however, has been the estimation of their genetic potentials for agriculture's benefit. Key to detecting in landraces and/or wild relatives of modern crops any allelic variant lost during domestication and crop improvement is the use of molecular information to determine structure, evolution, and function of the genes harbouring these alleles. This paper reviews some of the theoretical and statistical issues surrounding the use of molecular population genetics tools for the detection of agronomical valuable alleles in seed banks. Emphasis is made on the technical limitations imposed by seed banking that may lessen the success of integrated and multi-disciplinary molecular approaches. The influence that population stratification and linkage disequilibrium exert on specific experimental designs for a better understanding of the evolutionary history of potential agronomic-related genes is also examined.

Tracing the geographic origins of major avocado cultivars

It has been difficult to infer the genetic history of avocado breeding, owing to the role of hybridization in the origin of contemporary avocado cultivars. To address this difficulty, we used the model-based clustering program, STRUCTURE, and nucleotide polymorphism in 5960 bp of sequence from 4 nuclear loci to examine population structure in 21 wild avocado accessions. The origins of 33 cultivars were inferred relative to the wild sample. Nucleotide sequence diversity in domesticated avocados ranged between 80% and 90% of that observed for the same loci in wild avocado, depending on the diversity statistic used for comparison. Substantial genetic differentiation among 3 geographic groups of wild germplasm corresponded to the classically defined horticultural races of avocado. Previously undetected genetic differentiation was revealed in wild populations from Central Mexico, where 2 subpopulations were distinguished based on elevation and latitude.

Reducing cloning artifacts for recovery of allelic sequences by T7 endonuclease I cleavage and single re-extension of PCR products--a benchmark

Occurrence of chimeric sequences and related artifacts in PCR cloning procedures gives us risks of over-estimation of haplotypes or alleles. Recombination among haplotypes occurs through template switching during PCR cycles or through random repair of mismatch sites on heteroduplex DNA by the host cell. To eliminate the chimeric cloning artifacts, we tested two alternative protocols using T7 endonuclease I cleavage of mismatch sites and re-extension of nascent strands. Though T7 endonuclease I effectively eliminated chimeric clones in some cases, it produced many short fragments. Protocol with single re-extension of PCR products successfully recovered non-recombinant clones with fewer short fragments. In spite of the improvement of allelic recovery through these two protocols, there were still a few recombinants that remained in both reaction mixtures, and thus interpretation of the results for haplotype diversity in a PCR-amplified DNA population should be cautionary. Because re-extension in a diluted reaction mixture is quick, inexpensive and effective, it is advisable to use this procedure for recovery of chromosomal alleles with PCR cloning.

Retroelement insertional polymorphisms, diversity and phylogeography within diploid, D-genome Aegilops tauschii (Triticeae, Poaceae) sub-taxa in Iran

BACKGROUND AND AIMS

The diploid goat grass Aegilops tauschii (2n = 2x = 14) is native to the Middle East and is the D-genome donor to hexaploid bread wheat. The aim of this study was to measure the diversity of different subspecies and varieties of wild Ae. tauschii collected across the major areas where it grows in Iran and to examine patterns of diversity related to the taxa and geography.

METHODS

Inter-retroelement amplified polymorphism (IRAP) markers were used to analyse the biodiversity of DNA from 57 accessions of Ae. tauschii from northern and central Iran, and two hexaploid wheats. Key Results Eight IRAP primer combinations amplified a total of 171 distinct DNA fragments between 180 and 3200 bp long from the accessions, of which 169 were polymorphic. On average, about eight fragments were amplified with each primer combination, with more bands being amplified from accessions from the north-west of the country than from other accessions.

CONCLUSIONS

The IRAP markers showed high levels of genetic diversity. Analysis of all accessions together did not allow the allocation of individuals to taxa based on morphology, but showed a tendency to put accessions from the north-west apart from others regions. It is speculated that this could be due to different activity of retroelements in the different regions. Within the two taxa with most accessions, there was a range of IRAP genotypes that could be correlated closely with geographical origin. This supports suggestions that the centre of origin of the species is towards the south-east of the Caspian Sea. IRAP is an appropriate marker system to evaluate genetic diversity and evolutionary relationships within the taxa, but it is too variable to define the taxa themselves, where more slowly evolving morphological, DNA sequence or chromosomal makers may be more appropriate.

Molecular phylogeography of domesticated barley traces expansion of agriculture in the Old World

Barley (Hordeum vulgare ssp. vulgare) was first cultivated 10,500 years ago in the Fertile Crescent and is one of the founder crops of Eurasian agriculture. Phylogeographic analysis of five nuclear loci and morphological assessment of two traits in >250 domesticated barley accessions reveal that landraces found in South and East Asia are genetically distinct from those in Europe and North Africa. A Bayesian population structure assessment method indicates that barley accessions are subdivided into six clusters and that barley landraces from 10 different geographical regions of Eurasia and North Africa show distinct patterns of distribution across these clusters. Using haplotype frequency data, it appears that the Europe/North Africa landraces are most similar to the Near East population (F ST = 0.15) as well as to wild barley (F ST = 0.11) and are strongly differentiated from all other Asian populations (F ST = 0.34-0.74). A neighbor-joining analysis using these F ST estimates also supports a division between European, North African, and Near East barley types from more easterly Asian accessions. There is also differentiation in the presence of a naked caryopsis and spikelet row number between eastern and western barley accessions. The data support the differential migration of barley from two domestication events that led to the origin of barley--one in the Fertile Crescent and another farther east, possibly at the eastern edge of the Iranian Plateau--with European and North African barley largely originating from the former and much of Asian barley arising from the latter. This suggests that cultural diffusion or independent innovation is responsible for the expansion of agriculture to areas of South and East Asia during the Neolithic revolution.

Plant domestication, a unique opportunity to identify the genetic basis of adaptation

Despite the fundamental role of plant domestication in human history and the critical importance of a relatively small number of crop plants to modern societies, we still know little about adaptation under domestication. Here we focus on efforts to identify the genes responsible for adaptation to domestication. We start from a historical perspective, arguing that Darwin's conceptualization of domestication and unconscious selection provides valuable insight into the evolutionary history of crops and also provides a framework to evaluate modern methods used to decipher the genetic mechanisms underlying phenotypic change. We then review these methods, framing the discussion in terms of the phenotype-genotype hierarchy. Top-down approaches, such as quantitative trait locus and linkage disequilibrium mapping, start with a phenotype of interest and use genetic analysis to identify candidate genes. Bottom-up approaches, alternatively, use population genetic analyses to identify potentially adaptive genes and then rely on standard bioinformatics and reverse genetic tools to connect selected genes to a phenotype. We discuss the successes, advantages, and challenges of each, but we conclude that bottom-up approaches to understanding domestication as an adaptive process hold greater promise both for the study of adaptation and as a means to identify genes that contribute to agronomically important traits.

Maize Sh2 gene is constrained by natural selection but escaped domestication

In Zea mays L., we studied the molecular evolution of Shrunken2 (Sh2), a gene that encodes the large subunits of a major enzyme in endosperm starch biosynthesis, ADP-glucose pyrophosphorylase. We compared 4669 bp of the Sh2 coding region on 50 accessions of maize and teosinte. Very few nucleotide polymorphisms were found when compared with other genes in Z. mays, revealing an effect of purifying selection in the whole species that predates domestication. Additionally, the comparison of Sh2 sequences in all Z. mays subspecies and outgroups Z. diploperennis and Tripsacum dactyloides suggests the occurrence of an ancient selective sweep in the Sh2 3' region. The amount and nature of nucleotide diversity are similar in both maize and teosinte, confirming previous results that suggested that Sh2 has not been involved in maize domestication. The very low level of nucleotide diversity as well as the highly conserved protein sequence suggest that natural selection retained effective Sh2 allele(s) long before agriculture started, making human selection inefficient on this gene.

Genetic evidence for a second domestication of barley (Hordeum vulgare) east of the Fertile Crescent

Cereal agriculture originated with the domestication of barley and early forms of wheat in the Fertile Crescent. There has long been speculation that barley was domesticated more than once. We use differences in haplotype frequency among geographic regions at multiple loci to infer at least two domestications of barley; one within the Fertile Crescent and a second 1,500-3,000 km farther east. The Fertile Crescent domestication contributed the majority of diversity in European and American cultivars, whereas the second domestication contributed most of the diversity in barley from Central Asia to the Far East.

Adaptive climatic molecular evolution in wild barley at the Isa defense locus

Wild barley (Hordeum spontaneum) represents a significant genetic resource for crop improvement in barley (Hordeum vulgare) and for the study of the evolution and domestication of plant populations. The Isa gene from barley has a putative role in plant defense. This gene encodes a bifunctional alpha-amylase/subtilisin inhibitor that inhibits the bacterial serine protease subtilisin, fungal xylanase, and the plant's own alpha-amylase. The inhibition of plant alpha-amylases suggests this protein may also be important for grain quality from a human perspective. We identified 16 SNPs in the coding region of the Isa locus of 178 wild barley accessions from eight climatically divergent sites across Israel. The pattern of SNPs suggested a large number of recombination events within this gene, indicating that the low-outcrossing rate of wild barley is not a barrier to recombinant haplotypes becoming established in the population. Seven amino acid substitutions were present in the coding region. Genetic diversity for each population was calculated by using Nei's diversity index, and a Spearman rank correlation was carried out to test the association between gene diversity and 16 ecogeographical factors. Highly significant correlations were found between diversity at the Isa locus and key water variables, evaporation, rainfall, humidity, and latitude. The pattern of association suggests selective sweeps in the wetter climates, with resulting low diversity and weaker selection or diversifying selection in the dryer climates resulting in much higher diversity.

Effects of introgression and recombination on haplotype structure and linkage disequilibrium surrounding a locus encoding Bymovirus resistance in barley

We present a detailed analysis of linkage disequilibrium (LD) in the physical and genetic context of the barley gene Hv-eIF4E, which confers resistance to the barley yellow mosaic virus (BYMV) complex. Eighty-three SNPs distributed over 132 kb of Hv-eIF4E and six additional fragments genetically mapped to its flanking region were used to derive haplotypes from 131 accessions. Three haplogroups were recognized, discriminating between the alleles rym4 and rym5, which each encode for a spectrum of resistance to BYMV. With increasing map distance, haplotypes of susceptible genotypes displayed diverse patterns driven mainly by recombination, whereas haplotype diversity within the subgroups of resistant genotypes was limited. We conclude that the breakdown of LD within 1 cM of the resistance gene was generated mainly by susceptible genotypes. Despite the LD decay, a significant association between haplotype and resistance to BYMV was detected up to a distance of 5.5 cM from the resistance gene. The LD pattern and the haplotype structure of the target chromosomal region are the result of interplay between low recombination and recent breeding history.

Capturing diversity in the cereals: many options but little promiscuity

It is generally recognized by geneticists and plant breeders alike that there is a need to further improve the ability to capture and manipulate genetic diversity. The effective harnessing of diversity in traditional breeding programmes is limited and, therefore, it is vital that meiotic recombination can be manipulated given that it plays a pivotal role in generating diversity. With the advent of a wider range of genomics technologies, our understanding of meiotic processes should increase rapidly. Although comparative genetics has been useful, particularly in the broader grass family, the development of physical maps, long-range sequencing and transcript profiles promises to unravel the complexities of genomes as large or larger than wheat. Highlighting the most significant findings to date, this review pools the knowledge on these tools and reproductive processes.

Multilocus analysis of nucleotide variation of Oryza sativa and its wild relatives: severe bottleneck during domestication of rice

Varying degrees of reduction of genetic diversity in crops relative to their wild progenitors occurred during the process of domestication. Such information, however, has not been available for the Asian cultivated rice (Oryza sativa) despite its importance as a staple food and a model organism. To reveal levels and patterns of nucleotide diversity and to elucidate the genetic relationship and demographic history of O. sativa and its close relatives (Oryza rufipogon and Oryza nivara), we investigated nucleotide diversity data from 10 unlinked nuclear loci in species-wide samples of these species. The results indicated that O. rufipogon and O. nivara possessed comparable levels of nucleotide variation ((sil) = 0.0077 approximately 0.0095) compared with the relatives of other crops. In contrast, nucleotide diversity of O. sativa was as low as (sil) = 0.0024 and even lower ((sil) = 0.0021 for indica and 0.0011 for japonica), if we consider the 2 subspecies separately. Overall, only 20-10% of the diversity in the wild species was retained in 2 subspecies of the cultivated rice (indica and japonica), respectively. Because statistic tests did not reject the assumption of neutrality for all 10 loci, we further used coalescent to simulate bottlenecks under various lengths and population sizes to better understand the domestication process. Consistent with the dramatic reduction in nucleotide diversity, we detected a severe domestication bottleneck and demonstrated that the sequence diversity currently found in the rice genome could be explained by a founding population of 1,500 individuals if the initial domestication event occurred over a 3,000-year period. Phylogenetic analyses revealed close genetic relationships and ambiguous species boundary of O. rufipogon and O. nivara, providing additional evidence to treat them as 2 ecotypes of a single species. Lowest linkage disequilibrium (LD) was found in the perennial O. rufipogon where the r(2) value dropped to a negligible level within 400 bp, and the highest in the japonica rice where LD extended to the entirely sequenced region ( approximately 900 bp), implying that LD mapping by genome scans may not be feasible in wild rice due to the high density of markers needed.