Episodes of gene flow and selection during the evolutionary history of domesticated barley.

BACKGROUND: Barley is one of the founder crops of Neolithic agriculture and is among the most-grown cereals today. The only trait that universally differentiates the cultivated and wild subspecies is 'non-brittleness' of the rachis (the stem of the inflorescence), which facilitates harvesting of the crop. Other phenotypic differences appear to result from facultative or regional selective pressures. The population structure resulting from these regional events has been interpreted as evidence for multiple domestications or a mosaic ancestry involving genetic interaction between multiple wild or proto-domesticated lineages. However, each of the three mutations that confer non-brittleness originated in the western Fertile Crescent, arguing against multiregional origins for the crop. RESULTS: We examined exome data for 310 wild, cultivated and hybrid/feral barley accessions and showed that cultivated barley is structured into six genetically-defined groups that display admixture, resulting at least in part from two or more significant passages of gene flow with distinct wild populations. The six groups are descended from a single founding population that emerged in the western Fertile Crescent. Only a few loci were universally targeted by selection, the identity of these suggesting that changes in seedling emergence and pathogen resistance could represent crucial domestication switches. Subsequent selection operated on a regional basis and strongly contributed to differentiation of the genetic groups. CONCLUSIONS: Identification of genetically-defined groups provides clarity to our understanding of the population history of cultivated barley. Inference of population splits and mixtures together with analysis of selection sweeps indicate descent from a single founding population, which emerged in the western Fertile Crescent. This founding population underwent relatively little genetic selection, those changes that did occur affecting traits involved in seedling emergence and pathogen resistance, indicating that these phenotypes should be considered as 'domestication traits'. During its expansion out of the western Fertile Crescent, the crop underwent regional episodes of gene flow and selection, giving rise to a modern genetic signature that has been interpreted as evidence for multiple domestications, but which we show can be rationalized with a single origin.

Multiregional origins of the domesticated tetraploid wheats.

We used genotyping-by-sequencing (GBS) to investigate the evolutionary history of domesticated tetraploid wheats. With a panel of 189 wild and domesticated wheats, we identified 1,172,469 single nucleotide polymorphisms (SNPs) with a read depth ≥3. Principal component analyses (PCAs) separated the Triticum turgidum and Triticum timopheevii accessions, as well as wild T. turgidum from the domesticated emmers and the naked wheats, showing that SNP typing by GBS is capable of providing robust information on the genetic relationships between wheat species and subspecies. The PCAs and a neighbour-joining analysis suggested that domesticated tetraploid wheats have closest affinity with wild emmers from the northern Fertile Crescent, consistent with the results of previous genetic studies on the origins of domesticated wheat. However, a more detailed examination of admixture and allele sharing between domesticates and different wild populations, along with genome-wide association studies (GWAS), showed that the domesticated tetraploid wheats have also received a substantial genetic input from wild emmers from the southern Levant. Taking account of archaeological evidence that tetraploid wheats were first cultivated in the southern Levant, we suggest that a pre-domesticated crop spread from this region to southeast Turkey and became mixed with a wild emmer population from the northern Fertile Crescent. Fixation of the domestication traits in this mixed population would account for the allele sharing and GWAS results that we report. We also propose that feralization of the component of the pre-domesticated population that did not acquire domestication traits has resulted in the modern wild population from southeast Turkey displaying features of both the domesticates and wild emmer from the southern Levant, and hence appearing to be the sole progenitor of domesticated tetraploids when the phylogenetic relationships are studied by methods that assume a treelike pattern of evolution.

Diversity of a cytokinin dehydrogenase gene in wild and cultivated barley.

The cytokinin dehydrogenase gene HvCKX2.1 is the regulatory target for the most abundant heterochromatic small RNAs in drought-stressed barley caryopses. We investigated the diversity of HvCKX2.1 in 228 barley landraces and 216 wild accessions and identified 14 haplotypes, five of these with ten or more members, coding for four different protein variants. The third largest haplotype was abundant in wild accessions (51 members), but absent from the landrace collection. Protein structure predictions indicated that the amino acid substitution specific to haplotype 3 could result in a change in the functional properties of the HvCKX2.1 protein. Haplotypes 1-3 have overlapping geographical distributions in the wild population, but the average rainfall amounts at the collection sites for haplotype 3 plants are significantly higher during November to February compared to the equivalent data for plants of haplotypes 1 and 2. We argue that the likelihood that haplotype 3 plants were excluded from landraces by sampling bias that occurred when the first wild barley plants were taken into cultivation is low, and that it is reasonable to suggest that plants with haplotype 3 are absent from the crop because these plants were less suited to the artificial conditions associated with cultivation. Although the cytokinin signalling pathway influences many aspects of plant development, the identified role of HvCKX2.1 in the drought response raises the possibility that the particular aspect of cultivation that mitigated against haplotype 3 relates in some way to water utilization. Our results therefore highlight the possibility that water utilization properties should be looked on as a possible component of the suite of physiological adaptations accompanying the domestication and subsequent evolution of cultivated barley.

Diversity of a wall-associated kinase gene in wild and cultivated barley.

Domestication of barley and other cereals was accompanied by an increase in seed size which has been ascribed to human selection, large seeds being preferred by early farmers or favoured by cultivation practices such as deep sowing. An alternative suggestion is that the increase in seed size was an indirect consequence of selection for plants with more vigorous growth. To begin to address the latter hypothesis we studied the diversity of HvWAK1, a wall-associated kinase gene involved in root proliferation, in 220 wild barley accessions and 200 domesticated landraces. A 3655-bp sequence comprising the gene and upstream region contained 69 single nucleotide polymorphisms (SNPs), one indel and four short tandem repeats. A network of 50 haplotypes revealed a complex evolutionary relationship, but with landraces largely restricted to two parts of the topology. SNPs in the HvWAK1 coding region resulted in nonsynonymous substitutions at nine positions in the translation product, but none of these changes were predicted to have a significant effect on the protein structure. In contrast, the region upstream of the coding sequence contained five SNPs that were invariant in the domesticated population, fixation of these SNPs decreasing the likelihood that the upstream of a pair of TATA boxes and transcription start sites would be used to promote transcription of HvWAK1. The sequence diversity therefore suggests that the cis-regulatory region of HvWAK1 might have been subject to selection during barley domestication. The extent of root proliferation has been linked with traits such as above-ground biomass, so selection for particular cis-regulatory variants of HvWAK1 would be consistent with the hypothesis that seed size increases during domestication were the indirect consequence of selection for plants with increased growth vigour.

A discriminatory test for the wheat B and G genomes reveals misclassified accessions of Triticum timopheevii and Triticum turgidum.

The tetraploid wheat species Triticum turgidum and Triticum timopheevii are morphologically similar, and misidentification of material collected from the wild is possible. We compared published sequences for the Ppd-A1, Ppd-B1 and Ppd-G1 genes from multiple accessions of T. turgidum and T. timopheevii and devised a set of four polymerase chain reactions (PCRs), two specific for Ppd-B1 and two for Ppd-G1. We used these PCRs with 51 accessions of T. timopheevii and 20 of T. turgidum. Sixty of these accessions gave PCR products consistent with their taxon identifications, but the other eleven accessions gave anomalous results: ten accessions that were classified as T. turgidum were identified as T. timopheevii by the PCRs, and one T. timopheevii accession was typed as T. turgidum. We believe that these anomalies are not due to errors in the PCR tests because the results agree with a more comprehensive analysis of genome-wide single nucleotide polymorphisms, which similarly suggest that these eleven accessions have been misclassified. Our results therefore show that the accepted morphological tests for discrimination between T. turgidum and T. timopheevii might not be entirely robust, but that species identification can be made cheaply and quickly by PCRs directed at the Ppd-1 gene.

Origin of the Aromatic Group of Cultivated Rice (Oryza sativa L.) Traced to the Indian Subcontinent.

The aromatic group of Asian cultivated rice is a distinct population with considerable genetic diversity on the Indian subcontinent and includes the popular Basmati types characterized by pleasant fragrance. Genetic and phenotypic associations with other cultivated groups are ambiguous, obscuring the origin of the aromatic population. From analysis of genome-wide diversity among over 1,000 wild and cultivated rice accessions, we show that aromatic rice originated in the Indian subcontinent from hybridization between a local wild population and examples of domesticated japonica that had spread to the region from their own center of origin in East Asia. Most present-day aromatic accessions have inherited their cytoplasm along with 29-47% of their nuclear genome from the local Indian rice. We infer that the admixture occurred 4,000-2,400 years ago, soon after japonica rice reached the region. We identify aus as the original crop of the Indian subcontinent, indica and japonica as later arrivals, and aromatic a specific product of local agriculture. These results prompt a reappraisal of our understanding of the emergence and development of rice agriculture in the Indian subcontinent.

Role of genetic introgression during the evolution of cultivated rice (Oryza sativa L.).

BACKGROUND: Models for the origins of cultivated rice currently fall into two groups: ones that identify independent domestications of the indica, japonica and possibly also the aus types, and others that propose that the domestication phenotype was initially acquired by japonica, the underlying alleles then transferred by introgression to other pre-domesticated populations, giving the indica and aus varieties. Identifying the impact of past gene flow on cultivated rice genomes is therefore crucial to distinguishing between these models and understanding the domestication history of rice. To this end, we used population-scale polymorphism data to identify the progenitor gene pools of indica, japonica and aus. Variation shared among the cultivated groups but absent from at least one progenitor population was identified, and genomic blocks putatively transferred by gene flow among cultivated groups mapped. RESULTS: Introgression signals were absent at the major domestication loci (Prog1, Rc, qSH1, qSH3, Sh4) of indica and aus, indicating that these loci were unaffected by gene flow from japonica. Other domestication-related loci (Ghd7, LABA1, Kala4, LG1) show signals of introgression from japonica or indica to aus. There is a strong signal for LABA1 in japonica, possibly indicating introgression from indica. The indica genome is the least affected by gene flow, with just a few short regions with allelic frequencies slightly altered by introgression from japonica. CONCLUSION: Introgression has occurred during the evolution of cultivated rice, but was not responsible for transfer of the key domestication alleles between the cultivated groups. The results are therefore consistent with models in which japonica, indica and aus were domesticated independently, with each of these cultivated groups acquiring the domestication alleles from standing variation in wild rice, without a significant contribution from inter-group gene flow.

Origin of rice (Oryza sativa L.) domestication genes.

A number of genes that contribute to the domestication traits of cultivated rice have been identified. These include Sh4, Rc, PROG1 and LABA1, which are associated with non-shattering rachis, white pericarp, erect growth and barbless awns, respectively. The mutations giving rise to the "domestication alleles" of these genes are either invariable in cultivated rice, or have variability that is strictly associated with the phenotypic trait. This observation forms the basis to those current rice domestication models that envisage a single origin for the domesticated phenotype. Such models assume that the domestication alleles are absent or rare in wild rice, emerged under cultivation and spread across all rice groups by introgressive hybridization. We examined whole-genome sequencing datasets for wild and cultivated rice to test the former two assumptions. We found that the rc and laba1 alleles occur in wild rice with broad geographical distribution, and reach frequencies as high as 13 and 15%, respectively. These results are in agreement with previous observations of the prog1 and sh4 domestication alleles in wild populations. We also show that the diversity of the genomic regions surrounding the rc, laba1, prog1 and sh4 alleles in wild accessions is greater than that in cultivated rice, suggesting that these alleles emerged prior to domestication. Our findings indicate that the possibility that independent rice groups obtained identical domestication alleles directly from the wild population needs to be considered.

A novel mutation conferring the nonbrittle phenotype of cultivated barley.

The nonbrittle rachis, resulting in a seed head which does not shatter at maturity, is one of the key phenotypes that distinguishes domesticated barley from its wild relatives. The phenotype is associated with two loci, Btr1 and Btr2, with all domesticated barleys thought to have either a 1 bp deletion in Btr1 or an 11 bp deletion in Btr2. We used a PCR genotyping method with 380 domesticated barley landraces to identify those with the Btr1 deletion and those with the Btr2 deletion. We discovered two landraces, from Serbia and Greece, that had neither deletion. Instead these landraces possess a novel point mutation in Btr1, changing a leucine to a proline in the protein product. We confirmed that plants carrying this mutation have the nonbrittle phenotype and identified wild haplotypes from the Gaziantep region of southeast Turkey as the closest wild relatives of these two landraces. The presence of a third mutation conferring the nonbrittle phenotype of domesticated barley shows that the origin of this trait is more complex than previously thought, and is consistent with recent models that view the transition to agriculture in southwest Asia as a protracted and multiregional process.

Inability of 'Whole Genome Amplification' to Improve Success Rates for the Biomolecular Detection of Tuberculosis in Archaeological Samples.

We assessed the ability of whole genome amplification (WGA) to improve the efficiency of downstream polymerase chain reactions (PCRs) directed at ancient DNA (aDNA) of members of the Mycobacterium tuberculosis complex (MTBC). Using extracts from a variety of bones and a tooth from human skeletons with or without lesions indicative of tuberculosis, from multiple time periods, we obtained inconsistent results. We conclude that WGA does not provide any advantage in studies of MTBC aDNA. The sporadic nature of our results are probably due to the fact that WGA is itself a PCR-based procedure which, although designed to deal with fragmented DNA, might be inefficient with the low concentration of templates in an aDNA extract. As such, WGA is subject to similar, if not the same, restrictions as PCR when applied to aDNA.

Sex identification of ancient DNA samples using a microfluidic device.

Ancient DNA is the name given to the degraded, fragmented, and chemically damaged biomolecules that can be recovered from archaeological remains of plants, animals, and humans. Where ancient human DNA has survived at archaeological sites, it can give valuable information and is especially useful for its potential to identify kinship, population affinities, pathogens, and biological sex. Here, we describe the operation of a microfluidic device for the sex identification of ancient DNA samples using an efficient sample handling process. DNA is extracted from powdered bone samples and abasic sites labeled with biotin. Streptavidin-coated superparamagnetic particles are used to isolate the labeled DNA prior to amplification of the Amelogenin sex marker.

Three geographically separate domestications of Asian rice.

Domesticated rice (Oryza sativa L.) accompanied the dawn of Asian civilization(1) and has become one of world's staple crops. From archaeological and genetic evidence various contradictory scenarios for the origin of different varieties of cultivated rice have been proposed, the most recent based on a single domestication(2,3). By examining the footprints of selection in the genomes of different cultivated rice types, we show that there were three independent domestications in different parts of Asia. We identify wild populations in southern China and the Yangtze valley as the source of the japonica gene pool, and populations in Indochina and the Brahmaputra valley as the source of the indica gene pool. We reveal a hitherto unrecognized origin for the aus variety in central India or Bangladesh. We also conclude that aromatic rice is a result of a hybridization between japonica and aus, and that the tropical and temperate versions of japonica are later adaptations of one crop. Our conclusions are in accord with archaeological evidence that suggests widespread origins of rice cultivation(1,4). We therefore anticipate that our results will stimulate a more productive collaboration between genetic and archaeological studies of rice domestication, and guide utilization of genetic resources in breeding programmes aimed at crop improvement.

Genotyping of ancient Mycobacterium tuberculosis strains reveals historic genetic diversity.

The evolutionary history of the Mycobacterium tuberculosis complex (MTBC) has previously been studied by analysis of sequence diversity in extant strains, but not addressed by direct examination of strain genotypes in archaeological remains. Here, we use ancient DNA sequencing to type 11 single nucleotide polymorphisms and two large sequence polymorphisms in the MTBC strains present in 10 archaeological samples from skeletons from Britain and Europe dating to the second-nineteenth centuries AD. The results enable us to assign the strains to groupings and lineages recognized in the extant MTBC. We show that at least during the eighteenth-nineteenth centuries AD, strains of M. tuberculosis belonging to different genetic groups were present in Britain at the same time, possibly even at a single location, and we present evidence for a mixed infection in at least one individual. Our study shows that ancient DNA typing applied to multiple samples can provide sufficiently detailed information to contribute to both archaeological and evolutionary knowledge of the history of tuberculosis.

Biomolecular identification of ancient Mycobacterium tuberculosis complex DNA in human remains from Britain and continental Europe.

Tuberculosis is known to have afflicted humans throughout history and re-emerged towards the end of the 20th century, to an extent that it was declared a global emergency in 1993. The aim of this study was to apply a rigorous analytical regime to the detection of Mycobacterium tuberculosis complex (MTBC) DNA in 77 bone and tooth samples from 70 individuals from Britain and continental Europe, spanning the 1st-19th centuries AD. We performed the work in dedicated ancient DNA facilities designed to prevent all types of modern contamination, we checked the authenticity of all products obtained by the polymerase chain reaction, and we based our conclusions on up to four replicate experiments for each sample, some carried out in an independent laboratory. We identified 12 samples that, according to our strict criteria, gave definite evidence for the presence of MTBC DNA, and another 22 that we classified as "probable" or "possible." None of the definite samples came from vertebrae displaying lesions associated with TB. Instead, eight were from ribs displaying visceral new bone formation, one was a tooth from a skeleton with rib lesions, one was taken from a skeleton with endocranial lesions, one from an individual with lesions to the sacrum and sacroiliac joint and the last was from an individual with no lesions indicative of TB or possible TB. Our results add to information on the past temporal and geographical distribution of TB and affirm the suitability of ribs for studying ancient TB.

Reticulated origin of domesticated emmer wheat supports a dynamic model for the emergence of agriculture in the fertile crescent.

We used supernetworks with datasets of nuclear gene sequences and novel markers detecting retrotransposon insertions in ribosomal DNA loci to reassess the evolutionary relationships among tetraploid wheats. We show that domesticated emmer has a reticulated genetic ancestry, sharing phylogenetic signals with wild populations from all parts of the wild range. The extent of the genetic reticulation cannot be explained by post-domestication gene flow between cultivated emmer and wild plants, and the phylogenetic relationships among tetraploid wheats are incompatible with simple linear descent of the domesticates from a single wild population. A more parsimonious explanation of the data is that domesticated emmer originates from a hybridized population of different wild lineages. The observed diversity and reticulation patterns indicate that wild emmer evolved in the southern Levant, and that the wild emmer populations in south-eastern Turkey and the Zagros Mountains are relatively recent reticulate descendants of a subset of the Levantine wild populations. Based on our results we propose a new model for the emergence of domesticated emmer. During a pre-domestication period, diverse wild populations were collected from a large area west of the Euphrates and cultivated in mixed stands. Within these cultivated stands, hybridization gave rise to lineages displaying reticulated genealogical relationships with their ancestral populations. Gradual movement of early farmers out of the Levant introduced the pre-domesticated reticulated lineages to the northern and eastern parts of the Fertile Crescent, giving rise to the local wild populations but also facilitating fixation of domestication traits. Our model is consistent with the protracted and dispersed transition to agriculture indicated by the archaeobotanical evidence, and also with previous genetic data affiliating domesticated emmer with the wild populations in southeast Turkey. Unlike other protracted models, we assume that humans played an intuitive role throughout the process.