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.

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.

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.

Genotype of a historic strain of Mycobacterium tuberculosis.

The use of ancient DNA in paleopathological studies of tuberculosis has largely been restricted to confirmation of disease identifications made by skeletal analysis; few attempts at obtaining genotype data from archaeological samples have been made because of the need to perform different PCRs for each genetic locus being studied in an ancient DNA extract. We used a next generation sequencing approach involving hybridization capture directed at specific polymorphic regions of the Mycobacterium tuberculosis genome to identify a detailed genotype for a historic strain of M. tuberculosis from an individual buried in the 19th century St. George's Crypt, Leeds, West Yorkshire, England. We obtained 664,500 sequencing by oligonucleotide ligation and detection (SOLiD) reads that mapped to the targeted regions of the M. tuberculosis genome; the coverage included 218 of 247 SNPs, 10 of 11 insertion/deletion regions, and the repeat elements IS1081 and IS6110. The accuracy of the SOLiD data was checked by conventional PCRs directed at 11 SNPs and two insertion/deletions. The data placed the historic strain of M. tuberculosis in a group that is uncommon today, but it is known to have been present in North America in the early 20th century. Our results show the use of hybridization capture followed by next generation sequencing as a means of obtaining detailed genotypes of ancient varieties of M. tuberculosis, potentially enabling meaningful comparisons between strains from different geographic locations and different periods in the past.

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.

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

BACKGROUND: Understanding the evolution of cultivated barley is important for two reasons. First, the evolutionary relationships between different landraces might provide information on the spread and subsequent development of barley cultivation, including the adaptation of the crop to new environments and its response to human selection. Second, evolutionary information would enable landraces with similar traits but different genetic backgrounds to be identified, providing alternative strategies for the introduction of these traits into modern germplasm. RESULTS: The evolutionary relationships between 651 barley landraces were inferred from the genotypes for 24 microsatellites. The landraces could be divided into nine populations, each with a different geographical distribution. Comparisons with ear row number, caryopsis structure, seasonal growth habit and flowering time revealed a degree of association between population structure and phenotype, and analysis of climate variables indicated that the landraces are adapted, at least to some extent, to their environment. Human selection and/or environmental adaptation may therefore have played a role in the origin and/or maintenance of one or more of the barley landrace populations. There was also evidence that at least some of the population structure derived from geographical partitioning set up during the initial spread of barley cultivation into Europe, or reflected the later introduction of novel varieties. In particular, three closely-related populations were made up almost entirely of plants with the daylength nonresponsive version of the photoperiod response gene PPD-H1, conferring adaptation to the long annual growth season of northern Europe. These three populations probably originated in the eastern Fertile Crescent and entered Europe after the initial spread of agriculture. CONCLUSIONS: The discovery of population structure, combined with knowledge of associated phenotypes and environmental adaptations, enables a rational approach to identification of landraces that might be used as sources of germplasm for breeding programs. The population structure also enables hypotheses concerning the prehistoric spread and development of agriculture to be addressed.

The genetic expectations of a protracted model for the origins of domesticated crops.

Until recently, domestication has been interpreted as a rapid process with little predomestication cultivation and a relatively rapid rise of the domestication syndrome. This interpretation has had a profound effect on the biological framework within which investigations into crop origins have been carried out. A major underlying assumption has been that artificial selection pressures were substantially stronger than natural selection pressures, resulting in genetic patterns of diversity that reflect genetic independence of geographic localities. Recent archaeobotanical evidence has overturned the notion of a rapid transition, resulting in a protracted model that undermines these assumptions. Conclusions of genome-wide multilocus studies remain problematic in their support of a rapid-transition model by indicating that domesticated crops appear to be associated by monophyly with only a single geographic locality. Simulations presented here resolve this conflict, indicating that the results observed in such studies are inevitable over time at a rate that is largely influenced by the long-term population size. Counterintuitively, multiple origin crops are shown to be more likely to produce monophyletic clades than crops of a single origin. Under the protracted transition, the importance of the rise of the domestication syndrome becomes paramount in producing the patterns of genetic diversity from which crop origins may be deduced. We identify four different interacting levels of organization that now need to be considered to track crop origins from modern genetic diversity, making crop origins a problem that could be addressed through system-based approaches.

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.

Population-based resequencing reveals that the flowering time adaptation of cultivated barley originated east of the Fertile Crescent.

Gene resequencing and association analysis present new opportunities to study the evolution of adaptive traits in crop plants. Here we apply these tools to an extensive set of barley accessions to identify a component of the molecular basis of the flowering time adaptation, a trait critical to plant survival. Using an association-based study to relate variation in flowering time to sequence-based polymorphisms in the Ppd-H1 gene, we identify a causative polymorphism (SNP48) that accounts for the observed variation in barley flowering time. This polymorphism also shows latitude-dependent geographical distribution, consistent with the expected clinal variation in phenotype with the nonresponsive form predominating in the north. Networks, genealogies, and phylogenetic trees drawn for the Ppd-H1 haplotypes reveal population structure both in wild barley and in domesticated barley landraces. The spatial distribution of these population groups indicates that phylogeographical analysis of European landraces can provide information relevant to the Neolithic spread of barley cultivation and also has implications for the origins of domesticated barley, including those with the nonresponsive ppd-H1 phenotype. Haplotypes containing the nonresponsive version of SNP48 are present in wild barley accessions, indicating that the nonresponsive phenotype of European landraces originated in wild barley. The wild accessions whose nonresponsive haplotypes are most closely similar to those of landraces are found in Iran, within a region suggested as an area for domestication of barley east of the Fertile Crescent but which has previously been thought to have contributed relatively little to the diversity of European cultivars.

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.

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.

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.

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.

Novel methodology for construction and pruning of quasi-median networks.

BACKGROUND: Visualising the evolutionary history of a set of sequences is a challenge for molecular phylogenetics. One approach is to use undirected graphs, such as median networks, to visualise phylogenies where reticulate relationships such as recombination or homoplasy are displayed as cycles. Median networks contain binary representations of sequences as nodes, with edges connecting those sequences differing at one character; hypothetical ancestral nodes are invoked to generate a connected network which contains all most parsimonious trees. Quasi-median networks are a generalisation of median networks which are not restricted to binary data, although phylogenetic information contained within the multistate positions can be lost during the preprocessing of data. Where the history of a set of samples contain frequent homoplasies or recombination events quasi-median networks will have a complex topology. Graph reduction or pruning methods have been used to reduce network complexity but some of these methods are inapplicable to datasets in which recombination has occurred and others are procedurally complex and/or result in disconnected networks. RESULTS: We address the problems inherent in construction and reduction of quasi-median networks. We describe a novel method of generating quasi-median networks that uses all characters, both binary and multistate, without imposing an arbitrary ordering of the multistate partitions. We also describe a pruning mechanism which maintains at least one shortest path between observed sequences, displaying the underlying relations between all pairs of sequences while maintaining a connected graph. CONCLUSION: Application of this approach to 5S rDNA sequence data from sea beet produced a pruned network within which genetic isolation between populations by distance was evident, demonstrating the value of this approach for exploration of evolutionary relationships.

Microsatellite typing of ancient maize: insights into the history of agriculture in southern South America.

Archaeological maize specimens from Andean sites of southern South America, dating from 400 to 1400 years before present, were tested for the presence of ancient DNA and three microsatellite loci were typed in the specimens that gave positive results. Genotypes were also obtained for 146 individuals corresponding to modern landraces currently cultivated in the same areas and for 21 plants from Argentinian lowland races. Sequence analysis of cloned ancient DNA products revealed a high incidence of substitutions appearing in only one clone, with transitions prevalent. In the archaeological specimens, there was no evidence of polymorphism at any one of the three microsatellite loci: each exhibited a single allelic variant, identical to the most frequent allele found in contemporary populations belonging to races Amarillo Chico, Amarillo Grande, Blanco and Altiplano. Affiliation between ancient specimens and a set of races from the Andean complex was further supported by assignment tests. The striking genetic uniformity displayed by the ancient specimens and their close relationship with the Andean complex suggest that the latter gene pool has predominated in the western regions of southern South America for at least the past 1400 years. The results support hypotheses suggesting that maize cultivation initially spread into South America via a highland route, rather than through the lowlands.

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.

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.