Barley landraces are characterized by geographically heterogeneous genomic origins.

BACKGROUND: The genetic provenance of domesticated plants and the routes along which they were disseminated in prehistory have been a long-standing source of debate. Much of this debate has focused on identifying centers of origins for individual crops. However, many important crops show clear genetic signatures of multiple domestications, inconsistent with geographically circumscribed centers of origin. To better understand the genetic contributions of wild populations to domesticated barley, we compare single nucleotide polymorphism frequencies from 803 barley landraces to 277 accessions from wild populations. RESULTS: We find that the genetic contribution of individual wild populations differs across the genome. Despite extensive human movement and admixture of barley landraces since domestication, individual landrace genomes indicate a pattern of shared ancestry with geographically proximate wild barley populations. This results in landraces with a mosaic of ancestry from multiple source populations rather than discrete centers of origin. We rule out recent introgression, suggesting that these contributions are ancient. The over-representation in landraces of genomic segments from local wild populations suggests that wild populations contributed locally adaptive variation to primitive varieties. CONCLUSIONS: This study increases our understanding of the evolutionary process associated with the transition from wild to domesticated barley. Our findings indicate that cultivated barley is comprised of multiple source populations with unequal contributions traceable across the genome. We detect putative adaptive variants and identify the wild progenitor conferring those variants.

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.

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.

Tracing the geographic origins of weedy Ipomoea purpurea in the southeastern United States.

Ipomoea purpurea (common morning glory) is an annual vine native to Mexico that is well known for its large, showy flowers. Humans have spread morning glories worldwide, owing to the horticultural appeal of morning glory flowers. Ipomoea purpurea is an opportunistic colonizer of disturbed habitats including roadside and agricultural settings, and it is now regarded as a noxious weed in the Southeastern US. Naturalized populations in the Southeastern United States are highly polymorphic for a number of flower color morphs, unlike native Mexican populations that are typically monomorphic for the purple color morph. Although I. purpurea was introduced into the United States from Mexico, little is known about the specific geographic origins of US populations relative to the Mexican source. We use resequencing data from 11 loci and 30 I. purpurea accessions collected from the native range of the species in Central and Southern Mexico and 8 accessions from the Southeastern United States to infer likely geographic origins in Mexico. Based on genetic assignment analysis, haplotype composition, and the degree of shared polymorphism, I. purpurea samples from the Southeastern United States are genetically most similar to samples from the Valley of Mexico and Veracruz State. This supports earlier speculation that I. purpurea in the Southeastern United States was likely to have been introduced by European colonists from sources in Central Mexico.

Nucleotide sequence diversity of floral pigment genes in Mexican populations of Ipomoea purpurea (morning glory) accord with a neutral model of evolution.

The common morning glory (Ipomoea purpurea) is an annual vine native to Central and Southern Mexico. The genetics of flower color polymorphisms and interactions with the biotic environment have been extensively studied in I. purpurea and in its sister species I. nil. In this study, we examine nucleotide sequence polymorphism in 11 loci, 9 of which are known to participate in a pathway that produces floral pigments. A sample of 30 I. purpurea accessions from the native range of Central and Southern Mexico comprise the data, along with one accession from each of the two sister species I. alba and I. nil. We observe moderate levels of nucleotide sequence polymorphism of ~1%. The ratio of recombination to mutation parameter estimates (ρ/θ) of ~2.5 appears consistent with a mixed-mating system. Ipomoea resequencing data from these genic regions are noteworthy in providing a good fit to the standard neutral model of molecular evolution. The derived silent site frequency spectrum is very close to that predicted by coalescent simulations of a drift-mutation process, and Tajima's D values are not significantly different from expectations under neutrality.

Nucleotide diversity maps reveal variation in diversity among wheat genomes and chromosomes.

BACKGROUND: A genome-wide assessment of nucleotide diversity in a polyploid species must minimize the inclusion of homoeologous sequences into diversity estimates and reliably allocate individual haplotypes into their respective genomes. The same requirements complicate the development and deployment of single nucleotide polymorphism (SNP) markers in polyploid species. We report here a strategy that satisfies these requirements and deploy it in the sequencing of genes in cultivated hexaploid wheat (Triticum aestivum, genomes AABBDD) and wild tetraploid wheat (Triticum turgidum ssp. dicoccoides, genomes AABB) from the putative site of wheat domestication in Turkey. Data are used to assess the distribution of diversity among and within wheat genomes and to develop a panel of SNP markers for polyploid wheat. RESULTS: Nucleotide diversity was estimated in 2114 wheat genes and was similar between the A and B genomes and reduced in the D genome. Within a genome, diversity was diminished on some chromosomes. Low diversity was always accompanied by an excess of rare alleles. A total of 5,471 SNPs was discovered in 1791 wheat genes. Totals of 1,271, 1,218, and 2,203 SNPs were discovered in 488, 463, and 641 genes of wheat putative diploid ancestors, T. urartu, Aegilops speltoides, and Ae. tauschii, respectively. A public database containing genome-specific primers, SNPs, and other information was constructed. A total of 987 genes with nucleotide diversity estimated in one or more of the wheat genomes was placed on an Ae. tauschii genetic map, and the map was superimposed on wheat deletion-bin maps. The agreement between the maps was assessed. CONCLUSIONS: In a young polyploid, exemplified by T. aestivum, ancestral species are the primary source of genetic diversity. Low effective recombination due to self-pollination and a genetic mechanism precluding homoeologous chromosome pairing during polyploid meiosis can lead to the loss of diversity from large chromosomal regions. The net effect of these factors in T. aestivum is large variation in diversity among genomes and chromosomes, which impacts the development of SNP markers and their practical utility. Accumulation of new mutations in older polyploid species, such as wild emmer, results in increased diversity and its more uniform distribution across the genome.

Allele-specific PCR can improve the efficiency of experimental resolution of heterozygotes in resequencing studies.

Resolution of the two haplotypes present in an individual that is heterozygous at a locus has been a difficult problem for nucleotide sequence-based population genetic studies. Here, we demonstrate a method in which allele-specific polymerase chain reaction (AS-PCR) and computational phasing are combined for relatively high-throughput, efficient resolution of phase in resequencing studies. Using data from multiple loci that were fully experimentally phased, we demonstrate that the popular computational tool PHASE can accurately phase heterozygous individuals with common SNPs (single nucleotide polymorphisms) and/or common haplotypes. However, we also demonstrate that experimental phasing with AS-PCR can efficiently supplement computational phasing, providing a rapid means to phase individuals with rare SNPs or haplotypes and with heterozygous insertion/deletion polymorphisms. By following simple stepwise procedures, AS-PCR can result in much more efficient and accurate experimental phasing of haplotypes than is possible with traditional methods such as cloning.

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.

Persea americana (avocado): bringing ancient flowers to fruit in the genomics era.

The avocado (Persea americana) is a major crop commodity worldwide. Moreover, avocado, a paleopolyploid, is an evolutionary "outpost" among flowering plants, representing a basal lineage (the magnoliid clade) near the origin of the flowering plants themselves. Following centuries of selective breeding, avocado germplasm has been characterized at the level of microsatellite and RFLP markers. Nonetheless, little is known beyond these general diversity estimates, and much work remains to be done to develop avocado as a major subtropical-zone crop. Among the goals of avocado improvement are to develop varieties with fruit that will "store" better on the tree, show uniform ripening and have better post-harvest storage. Avocado transcriptome sequencing, genome mapping and partial genomic sequencing will represent a major step toward the goal of sequencing the entire avocado genome, which is expected to aid in improving avocado varieties and production, as well as understanding the evolution of flowers from non-flowering seed plants (gymnosperms). Additionally, continued evolutionary and other comparative studies of flower and fruit development in different avocado strains can be accomplished at the gene expression level, including in comparison with avocado relatives, and these should provide important insights into the genetic regulation of fruit development in basal angiosperms.

Nucleotide diversity and linkage disequilibrium in wild avocado (Persea americana Mill.).

Resequencing studies provide the ultimate resolution of genetic diversity because they identify all mutations in a gene that are present within the sampled individuals. We report a resequencing study of Persea americana, a subtropical tree species native to Meso- and Central America and the progenitor of cultivated avocado. The sample includes 21 wild accessions from Mexico, Costa Rica, Ecuador, and the Dominican Republic. Estimated levels of nucleotide polymorphism and linkage disequilibrium (LD) are obtained from fully resolved haplotype data from 4 nuclear loci that span 5960 nucleotide sites. Results show that, although avocado is a subtropical tree crop and a predominantly outcrossing plant, the overall level of genetic variation is not exceptionally high (nucleotide diversity at silent sites, pi(sil) = 0.0102) compared with available estimates from temperate plant species. Intralocus LD decays rapidly to half the initial value within about 1 kb. Estimates of recombination rate (based on the sequence data) show that the rate is not exceptionally high when compared with annual plants such as wild barley or maize. Interlocus LD is significant owing to substantial population structure induced by mixing of the 3 botanical races of avocado.

Error detection in SNP data by considering the likelihood of recombinational history implied by three-site combinations.

MOTIVATION: Errors in nucleotide sequence and SNP genotyping data are problematic when inferring haplotypes. Previously published methods for error detection in haplotype data make use of pedigree information; however, for many samples, individuals are not related by pedigree. This article describes a method for detecting errors in haplotypes by considering the recombinational history implied by the patterns of variation, three SNPs at a time. RESULTS: Coalescent simulations provide evidence that the method is robust to high levels of recombination as well as homologous gene conversion, indicating that patterns produced by both proximate and distant SNPs may be useful for detecting unlikely three-site haplotypes. AVAILABILITY: The perl script implementing the described method is called EDUT (Error Detection Using Triplets) and is available on request from the authors. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

Estimating the contribution of mutation, recombination and gene conversion in the generation of haplotypic diversity.

Recombination occurs through both homologous crossing over and homologous gene conversion during meiosis. The contribution of recombination relative to mutation is expected to be dramatically reduced in inbreeding organisms. We report coalescent-based estimates of the recombination parameter (rho) relative to estimates of the mutation parameter (theta) for 18 genes from the highly self-fertilizing grass, wild barley, Hordeum vulgare ssp. spontaneum. Estimates of rho/theta are much greater than expected, with a mean rho/theta approximately 1.5, similar to estimates from outcrossing species. We also estimate rho with and without the contribution of gene conversion. Genotyping errors can mimic the effect of gene conversion, upwardly biasing estimates of the role of conversion. Thus we report a novel method for identifying genotyping errors in nucleotide sequence data sets. We show that there is evidence for gene conversion in many large nucleotide sequence data sets including our data that have been purged of all detectable sequencing errors and in data sets from Drosophila melanogaster, D. simulans, and Zea mays. In total, 13 of 27 loci show evidence of gene conversion. For these loci, gene conversion is estimated to contribute an average of twice as much as crossing over to total recombination.

Low levels of linkage disequilibrium in wild barley (Hordeum vulgare ssp. spontaneum) despite high rates of self-fertilization.

High levels of inbreeding cause populations to become composed of homozygous, inbred lines. High levels of homozygosity limit the effectiveness of recombination, and therefore, retard the rate of decay of linkage (gametic phase) disequilibrium (LD) among mutations. Inbreeding and recombination interact to shape the expected pattern of LD. The actual extent of nucleotide sequence level LD within inbreeding species has only been studied in Arabidopsis, a weedy species whose global range has recently expanded. In the present study, we examine the levels of LD within and between 18 nuclear genes in 25 accessions from across the geographic range of wild barley, a species with a selfing rate of approximately 98%. In addition to examination of intralocus LD, we employ a resampling method to determine whether interlocus LD exceeds expectations. We demonstrate that, for the majority of wild barley loci, intralocus LD decays rapidly, i.e., at a rate similar to that observed in the outcrossing species, Zea mays (maize). Excess interlocus LD is observed at 15% of two-locus combinations; almost all interlocus LD involves loci with significant geographic structuring of mutational variation.

Evolutionary dynamics of the DNA-binding domains in putative R2R3-MYB genes identified from rice subspecies indica and japonica genomes.

The molecular evolution of the R2R3-MYB gene family is of great interest because it is one of the most important transcription factor gene families in the plant kingdom. Comparative analyses of a gene family may reveal important adaptive changes at the protein level and thereby provide insights that relate structure to function. We have performed a range of comparative and bioinformatics analyses on R2R3-MYB genes identified from the rice (Oryza sativa subsp. japonica and indica) and Arabidopsis genome sequences. The study provides an initial framework to investigate how different evolutionary lineages in a gene family evolve new functions. Our results reveal a remarkable excess of non-synonymous substitutions, an indication of adaptive selection on protein structure that occurred during the evolution of both helix1 and helix2 of rice R2R3-MYB DNA-binding domains. These flexible alpha-helix regions associated with high frequencies of excess non-synonymous substitutions may play critical roles in the characteristic packing of R2R3-MYB DNA-binding domains and thereby modify the protein-DNA interaction process resulting in the recognition of novel DNA-binding sites. Furthermore, a co-evolutionary pattern is found between the second alpha-helix of the R2 domain and the second alpha-helix of the R3 domain by examining all the possible alpha-helix pairings in both the R2 and R3 domains. This points to the functional importance of pairing interactions between related secondary structures.

Genes that determine flower color: the role of regulatory changes in the evolution of phenotypic adaptations.

A central goal of evolutionary genetics is to trace the causal pathway between mutations at particular genes and adaptation at the phenotypic level. The proximate objective is to identify adaptations through the analysis of molecular sequence data from specific candidate genes or their regulatory elements. In this paper, we consider the molecular evolution of floral color in the morning glory genus (Ipomoea) as a model for relating molecular and phenotypic evolution. To begin, flower color variation usually conforms to simple Mendelian transmission, thus facilitating genetic and molecular analyses. Population genetic studies of flower color polymorphisms in the common morning glory (Ipomoea purpurea) have shown that some morphs are subject to complex patterns of selection. Striking differences in floral color and morphology are also associated with speciation in the genus Ipomoea. The molecular bases for these adaptive shifts can be dissected because the biosynthetic pathways that determine floral pigmentation are well understood and many of the genes of flavonoid biosynthesis have been isolated and extensively studied. We present a comparative analysis of the level of gene expression in Ipomoea for several key genes in flavonoid biosynthesis. Specifically we ask: how frequently are adaptive shifts in flower color phenotypes associated with changes in regulation of gene expression versus mutations in structural genes? The results of this study show that most species differences in this crucial phenotype are associated with changes in the regulation of gene expression.

Excess non-synonymous substitutions suggest that positive selection episodes occurred during the evolution of DNA-binding domains in the Arabidopsis R2R3-MYB gene family.

It has been suggested that evolutionary changes in regulatory genes may be the predominant molecular mechanism governing both physiological and morphological evolution. R2R3-AtMYB is one of the largest transcription factor gene families in Arabidopsis. Using inferred ancestral sequences we show that several lineages in the R2R3-AtMYB phylogeny experienced excess non-synonymous nucleotide substitution upon gene duplication, indicating episodes of positive selection driving adaptive shifts early in the evolution of this gene family. A noise reduction technique was then used to determine individual sites in DNA-binding domains (R2 domain and R3 domain) of R2R3-AtMYB protein sequence that were favored by frequent non-synonymous substitutions. The analyses reveal that the first helix (helix1) and the second helix (helix2) in both R2 and R3 domains are characterized by more frequent non-synonymous substitutions, and thus experienced significantly higher positive selection pressure than the third helix (helix3) in both domains. Previous MYB protein structure studies have suggested that helix1 and helix2 in both R2 and R3 domains are involved in the characteristic packing of R2R3-AtMYB DNA-binding domains. This suggests that excess non-synonymous substitutions in these helices could have resulted in MYB recognition of novel gene target sites.

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

Mutations arise in a single individual and at a single point in time and space. The geographic distribution of mutations reflects both historical population size and frequency of migration. We employ coalescence-based methods to coestimate effective population size, frequency of migration, and level of recombination compatible with observed genealogical relationships in sequence data from nine nuclear genes in wild barley (Hordeum vulgare ssp. spontaneum), a highly self-fertilizing grass species. In self-fertilizing plants, gamete dispersal is severely limited; dissemination occurs primarily through seed dispersal. Also, heterozygosity is greatly reduced, which renders recombination less effective at randomizing genetic variation and causes larger portions of the genome to trace a similar history. Despite these predicted effects of this mating system, the majority of loci show evidence of recombination. Levels of nucleotide variation and the patterns of geographic distribution of mutations in wild barley are highly heterogeneous across loci. Two of the nine sampled loci maintain highly diverged, geographic region-specific suites of mutations. Two additional loci include region-specific haplotypes with a much shallower coalescence. Despite inbreeding, sessile growth habit, and the observation of geographic structure at almost half of sampled loci, parametric estimates of migration suggest that seed dispersal is sufficient for migration across the approximately 3,500-km range of the species. Recurrent migration is also evident based on the geographic distribution of mutational variation at some loci. At one locus a single haplotype has spread rapidly enough to occur, unmodified by mutation, across the range of the species.

Tracing floral adaptations from ecology to molecules.

Flowers have long fascinated humans. The scientific study of floral biology unifies many diverse areas of research, ranging from systematics to ecology, and from genetics to molecular biology. Despite this unity, few plant species offer the experimental versatility to encompass all levels of biological investigation in a single system. An exception is the morning glory genus Ipomoea, in which a broad picture of floral evolution, ranging from ecology to molecular biology, is emerging.