Original Research ARTICLE
Origin of worldwide cultivated barley revealed by NAM-1 gene and grain protein content
- 1College of Plant Science and Technology, Huazhong Agricultural University, Wuhan, China
- 2Department of Biology, Saint Mary’s University, Halifax, NS, Canada
The origin, evolution, and distribution of cultivated barley provides powerful insights into the historic origin and early spread of agrarian culture. Here, population-based genetic diversity and phylogenetic analyses were performed to determine the evolution and origin of barley and how domestication and subsequent introgression have affected the genetic diversity and changes in cultivated barley on a worldwide scale. A set of worldwide cultivated and wild barleys from Asia and Tibet of China were analyzed using the sequences for NAM-1 gene and gene-associated traits-grain protein content (GPC). Our results showed Tibetan wild barley distinctly diverged from Near Eastern barley, and confirmed that Tibet is one of the origin and domestication centers for cultivated barley, and in turn supported a polyphyletic origin of domesticated barley. Comparison of haplotype composition among geographic regions revealed gene flow between Eastern and Western barley populations, suggesting that the Silk Road might have played a crucial role in the spread of genes. The GPC in the 118 cultivated and 93 wild barley accessions ranged from 6.73 to 12.35% with a mean of 9.43%. Overall, wild barley had higher averaged GPC (10.44%) than cultivated barley. Two unique haplotypes (Hap2 and Hap7) caused by a base mutations (at position 544) in the coding region of the NAM-1 gene might have a significant impact on the GPC. Single nucleotide polymorphisms and haplotypes of NAM-1 associated with GPC in barley could provide a useful method for screening GPC in barley germplasm. The Tibetan wild accessions with lower GPC could be useful for malt barley breeding.
Wild barley, Hordeum spontaneum C. Koch, is the progenitor of cultivated barley, Hordeum vulgare L. As one of the earliest domesticated crops, barley has been one of most important staple crops in old world Neolithic agriculture upon which early agriculture was built (Harlan and Zohary, 1966; Zohary and Hopf, 2000). The domestication of barley is fundamental to understanding the origin and early diffusion of agrarian culture (Morrell and Clegg, 2007).
The geographic range of H. spontaneum was clearly defined, and the Fertile Crescent has been suggested as the only location where barley was domesticated by a large number of researchers (Harlan and Zohary, 1966; Nevo et al., 1984, 1986; Pakniyat et al., 1997; Badr et al., 2000; Nevo, 2006; Zohary et al., 2012). However, the centre of origin of barley has not been fully resolved. H. spontaneum, the wild ancestor of cultivated barley, has been discovered in several geographically distinct locations other than the Fertile Crescent, such as Morocco, Algeria, Libya, Egypt, Crete, Ethiopia, and Tibet, leading to the proposal of a multicentric origin for this crop (Åberg, 1938; Bekele, 1983; Molina-Cano et al., 1987; Von Bothmer et al., 1995; Molina-Cano et al., 2002; Fuller et al., 2011; Von Bothmer and Komatsuda, 2011), and was supported by additional studies (Takahashi and Hayashi, 1964; Azhaguvel and Komatsuda, 2007; Fuller et al., 2011). For instance, archeological evidence suggests a diffuse “center” of origin for barley (Fuller et al., 2011), and non-brittle rachis in oriental and occidental lines is controlled by two distinct genetic loci, indicating independent origins of oriental and occidental barleys (Takahashi and Hayashi, 1964; Azhaguvel and Komatsuda, 2007). In addition, numerous studies have reported clear genetic differentiation among barley populations from Eastern and Southern Asia and those from Western Asia, Europe, and North Africa (Azhaguvel and Komatsuda, 2007; Saisho and Purugganan, 2007; Fu and Peterson, 2011; Morrell et al., 2013). Recent molecular evidence indicated that an additional center of wild barley domestication event may occur in Central Asia at the eastern edge of the Iranian Plateau. It is assumed that this area constitutes the center of origin for cultivated barley from South and East Asia (Morrell and Clegg, 2007; Jones et al., 2008). Most noticeably, the Qinghai–Tibetan Plateau has been proposed as one of the centers of origin of cultivated barley. Recent Diversity array technology (DArT) data and transcriptome profiling suggested that the Tibetan Plateau and its vicinity was one of the centers of barley domestication, and further confirmed the multiple origins of cultivated barley (Dai et al., 2012, 2014).
The resequencing of candidate gene loci within diverse populations has implication for understanding the origin of barley domestication (Morrell and Clegg, 2007; Jones et al., 2008). Domestication has genetically not only transformed the brittle rachis, tenacious glume, and non-free thresh-ability, but also modified yield and yield components in many crops (Peng et al., 2011; Yan et al., 2015). Grain protein content (GPC) is a very important quality determinant in many cereals. In barley, GPC is closely related to feeding quality as well as malting and brewing processes. High protein concentration is attributed to feeding quality, while low protein content is favorable for barley malt and beer production (See et al., 2002). However, it is difficult to improve simultaneously grain yield (GY) and grain protein concentration (Bogard et al., 2010). Previous studies have demonstrated a strong genetic negative correlation between GPC and yield in various cereals (Oury and Godin, 2007; Bogard et al., 2011; Blanco et al., 2012; Martre et al., 2015). The content of grain protein (GPC) is determined by numerous factors including environmental effects, the elements determined yield, and complex genetic system (Simmonds, 1995). Whole-plant senescence processes overlap with grain filling, and the synchronization of these two processes is important for affecting yield and quality characteristics such as grain protein content (GPC; Distelfeld et al., 2014).
The NAC transcription factor family plays a central role in regulating organ and organism senescence (Uauy et al., 2006; Balazadeh et al., 2010; Zhu et al., 2012; Jensen and Skriver, 2014). It has been shown that a NAC transcription factor (TtNAM-B1) was related to the GPC of wheat (Uauy et al., 2006; Peng et al., 2013). In barley, the relationships between plant senescence, NAC gene expression and GPC were also studied (Cai et al., 2013). Two orthologs genes of a NAC transcription factor, HvNAM-1 and HvNAM-2 (GenBank accession number DQ869678 and DQ869679) have been mapped on chromosomes 6H and 2H, respectively (Uauy et al., 2006; Distelfeld et al., 2008). Allelic variation of the NAM-1 gene seems to be related to protein content (Distelfeld et al., 2008; Jukanti and Fischer, 2008; Jamar et al., 2010). Recently, a genome-wide association study (GWAS) and a multi-platform candidate gene-based association analysis for cultivated and wild barley found that the two HvNAM genes might play a role in controlling GPC in barley (Cai et al., 2013).
The origin and domestication of cultivated barley have been widely discussed; however, the debate on these subjects still remains. To provide further evidence to determine the evolution and origin of barley and how barley domestication and subsequent introgression have affected the genetic diversity of cultivated barley, we characterized nucleotide diversity of the NAM-1 gene and gene-associated traits-GPC in wild barley from Southwest Asia, Central Asia, Tibet of China, and cultivated barley from different parts of the world. Our primary goals were (i) to examine genetic differentiation between cultivated barley and wild-barley populations; (ii) to elucidate the origin and spread of worldwide cultivated barley; and (iii) to investigate the GPC in barley populations and its relationship to NAM-1 gene.
Materials and Methods
Total of 218 barley accessions were used in this study, including 94 accessions of wild barley (H. spontaneum) and 124 lines of cultivated barley (H. vulgare). The wild barley populations included: 53 wild barley accessions from the Southwest Asia (Israel, Jordan, Ethiopia, Lebanon, Azerbaijan, Syria, Turkey, and Iraq); 21 wild barley from Central Asia (Iran, Afghanistan, Tajikistan, and Pakistan), and 20 wild barley from Tibet of China. One hundred and twenty-four cultivated barley (H. vulgare) accessions were from 18 countries (61 form Eastern Asia, 8 from South America, 18 form North America, 10 from Mediterranean coast, 5 form Australia, and 22 from Europe). Those materials were provided by USDA (United States Department of Agriculture) and the Huazhong Agricultural University barley germplasm collection. The accessions names and their geographical origin were given in the Supplementary Table S1.
DNA Extraction, NAM-1 Gene Amplification, and Sequencing
The seeds were planted in pots with nutrient soil, and grown in a growth chamber with 14 h of light at 22°C and 10 h of darkness at 18°C prior to DNA extraction. Young leaves were collected from 5 to 10 plants of each accession and ground to a fine powder with liquid nitrogen and stored at –80°C until DNA extraction. DNA was extracted according to the method of Stein et al. (2001).
The NAM-1 gene sequence was amplified using polymerase chain reaction (PCR). The forward and reverse primer sequences are Nam-1HF: 5′-TATCAAGCGCCGTAATTTCC-3′ and Nam-1HR: 5′-ATACTGCCGACGTTTCTGCT-3′, respectively (Ren et al., 2013). Amplification of DNA was carried out in 40 μl reaction mixture containing 60 ng template DNA, 0.2 μM of each primer, 1.5 mM MgCl2, 0.2 mM of each deoxynucleoside triphosphate (dNTP), and 1.5 unit of high-fidelity polymerase ExTaq (TaKaRa, Dalian, China). PCR amplification was performed with an initial denaturing of 4 min at 95°C, followed by 40 cycles of 1 min at 95°C, 1 min at 52°C, 2 min at 72°C, and ending with an 8 min extension at 72°C.
Amplified products were electrophoresed in 1.0% agarose gel, and purified using the QIAquickTM PCR purification kit (QIAGEN Inc) according to the manufacturer’s protocol. DNA was commercially sequenced at the Beijing TsingKe BioTech Co., Ltd (Beijing, China). The amplification and sequencing were repeated three times to exclude sequencing errors introduced by TaqDNA polymerase during PCR amplification. In addition, both forward and reverse strands were sequenced independently, and further checked for data quality using Chromas 2.32 (Technelysium Pty. Ltd.). Sequences of NAM-1 genes in 214 barley accessions were shown in Supplementary Data Sheet S1.
Grain Protein Content (GPC) Measurement
The GPC in 211 barley accessions was determined. All barley accessions were planted at the experimental field of Ezhou (Hubei, China, 114.41°E. 30.06°N) in the early winter of 2013, and were cultivated in the same conditions with identical agronomic managements until maturity and harvested. Mature grains were fully ground and passed through a 0.5-mm screen after dehydration until constant mass is reached. The GPC was measured using the method of Kjeldahl (1883), three measurements were done for each sample. Protein content is calculated with a factor of 6.25 for N content (Mariotti et al., 2008).
Multiple sequences were aligned using ClustalX (Thompson et al., 1997). Genetic diversity was estimated by Tajima’s (1989) π and Watterson’s (1975) statistics using DnaSP version 5.0 (Librado and Rozas, 2009), and the tests of neutral evolution were performed as described by Tajima (1989) and Fu and Li (1993). Phylogenetic analysis was performed with the computer program MEGA 6 (Tamura et al., 2013). The phylogenetic tree of the 214 accessions was constructed using neighbor-joining (Saitou and Nei, 1987) methods with Tajima–Nei model. The confidence of each clade was evaluated by the bootstrap values with 1,000 replicates. Statistical analyses were performed by using SAS 9.0 software (SAS Institute, Cary, North Carolina, USA). Significance between groups was evaluated by one-way analysis of variance (ANOVA) followed by a Newman–Keuls post hoc test, a P value of <0.05 was considered statistically significant.
Haplotype Frequency Analysis in Barley Populations
Eight Single nucleotide polymorphisms (SNPs) and 10 distinguishable haplotypes were identified in 214 barley accessions (Table 1). Haplotypes in three wild barley populations and six cultivated barley populations were compared and are shown in Table 2. Ten haplotypes were detected across 94 wild barley accessions. Of these, seven haplotypes were detected in the Southwest Asian wild population, five in the Tibetan wild barley population, and two in the Central Asian wild barley population. Moreover, seven haplotypes were population-specific, four specific to the Southwest Asian wild population, and three specific to the Tibetan wild barley population. In contrast, only three haplotypes were observed in cultivars, with all three present in the North American and European cultivated barley populations, and two of them in the other four cultivar populations. No population-specific haplotype was found in the cultivated barley populations. The three haplotypes identified in cultivars were also identified in both Southwest Asian and Tibetan wild barley populations, suggesting that two major wild-barley populations, the Near East Fertile Crescent and Tibetan Plateau populations, might have contributed to the origin of cultivated barley.
TABLE 2. Haplotype frequencies of NAM-1 gene in three wild barley populations and six landrace barley populations.
The haplotype frequencies present in all accessions ranged from 0.005 to 0.449. Among all the haplotypes, three haplotypes were detected in 193 barley accessions. The haplotype Hap2 appeared in 96 accessions (44.9%), Hap7 was observed in 80 accessions (37.4%), and Hap1 was present in 17 accessions (7.9%). Of the seven haplotypes present in <2% of the accessions sampled, five were unique to the specific wild populations, i.e., Hap3 and Hap8 was unique to the Southwest Asian wild barley population and the Tibetan wild barley populations, respectively. The frequencies of NAM-1 haplotypes differed markedly among different geographic populations. The haplotype Hap1 was the most frequent one in the Southwest Asian wild barley population (0.264), but rare in the North American and European landraces (0.063 and 0.10, respectively), and absent in the remaining six populations. The rare haplotypes either confined or occurred in specific geographic regions, i.e., the Hap2 was detected in a few landrace accessions from North America (0.063), and the Hap3 (0.075), Hap5 (0.057), and Hap6 (0.057) were population-specific to wild barley in Southwest Asia; Hap8 and Hap10 were unique to the wild population in Tibet (Table 2; Figure 1).
FIGURE 1. Geographic distribution of wild barley populations and landrace populations. NAM-1 haplotype frequencies among nine geographic regions were displayed in pie diagrams and the exact proportions of each are given in percent by the corresponding color code.
Genetic Diversity Analysis and Neutrality Test
Genetic variation analyses revealed that wild barley, except Central Asian populations, showed higher haplotype diversity than landraces. The highest diversity was detected in the Southwest Asian wild barley population (0.722), followed by the Tibetan wild barley population (0.679). The pre-site nucleotide diversity ranged from 0.00023 ± 0.00023 (East Asian landrace population) to 0.00141 ± 0.00067 (Southwest Asian wild population). Higher values were also discovered in the Southwest Asian and Tibetan wild barley populations. These estimates corresponded well to the number of NAM-1 haplotypes (7 and 5, respectively), as well as the diversity based on the number of segregating sites (0.00100 and 0.00103, respectively; Table 3).
Tajima (1989) and Fu and Li’s (1993) neutrality tests were performed to determine whether the observed genetic diversity follows an equilibrium neutral model. Both positive values were obtained from the East Asian, South American, Mediterranean Coast and European landrace populations as well as Tibetan wild barley populations. In contrast, negative values for both tests were obtained from the North American and Australian landrace populations. Moreover, Tajima’s D values were negative for the Central Asian and Southwest Asian wild barley populations (–0.61772 and –0.72915, respectively), while the Fu and Li test values were positive for these two populations. None of the values were statistically significant (P = 0.05; Table 3).
Phylogenetic Analysis of NAM-1 Gene
A Phylogenetic tree was constructed to depict genetic relationships among the 214 samples based on the NAM-1 gene (Figure 2). The neighbor-joining analysis placed these samples into two major clusters, one comprised of the majority of wild barley accessions (red bar in Figure 2) and another comprised of the majority of cultivated barley accessions (green bar in Figure 2). However, some Tibetan wild barleys were distinct from the Near Eastern and Central Asian wild barleys, and appeared in the cultivars-dominated cluster. The third cluster comprised of four accession of wild barley from Tibet, three accessions of wild barley from Southwest Asia, and one accession of wild barley from Central Asia.
FIGURE 2. Phylogenetic tree of 214 barley accessions based on the NAM-1 gene. Two major clusters, one comprised of a majority of wild barley accessions (represented in green bar) and another comprised of a majority of cultivated barley accessions (represented in red bar) are separated. The square stands for wild barley accessions: Tibet (Wb-T, red), Southwest Asia (Wb-S, purple), and Central Asia (Wb-C, orange), respectively; the triangle indicates landrace barleys: East Asia (Lb-EA, black), North America (Lb-NA, blue), South America (Lb-SA, pink), Mediterranean Coast Areas (Lb-MA, green), Europe (Lb-EU, yellow), and Australia (Lb-AU, orange).
The Variation of Protein Content
The mean value of GPC and the differences among populations are shown in Figure 3. The GPC in 118 cultivated and 93 wild barley accessions ranged from 6.73 to 12.35% with a mean of 9.43%. Overall, wild barley had higher averaged GPC (10.44%) than cultivated barley. A significantly statistical difference was found between the group of landraces and the group of Southwest Asian and Central Asian wild barleys; however, no significant difference was found within each of these.
FIGURE 3. The means of grain protein content (GPC) variation among nine populations. Different letters (from a–c) on top of the histogram bars correspond to classes of which the population belongs, based on the Newman–Keuls test. Error bars indicate standard error. The populations used are Wb-T (Wild barley of Tibet), Wb-C (Wild barley of Central Asia), and Wb-S (Wild barley of Southwest Asia); Lb-EA (Landrace barley of East Asia), Lb-NA (Landraces of North America), Lb-SA (Landraces of South America), Lb-MA (Landraces of the Mediterranean Coast Areas), Lb-EU (Landraces of Europe), and Lb-AU (Landraces of Australia).
Association between SNP and GPC
To determine association between SNP of NAM-1 and GPC among barley populations, a sequence from NCBI (accession number DQ869678) was used as a reference to identify SNP in our populations. The SNPs identified were summarized in Table 4. Three SNPs were recognized in the Tibetan wild barley population and are located at position 375, 544, and 1190, respectively; the SNPs at position 544 and 1190 are within the coding sequence. Seven SNPs were found in the Southwest Asian wild barley population, of which two SNPs (position 473 and 823) were located in a non-coding region, and five SNPs (position 507, 544, 616, 1190, and 1253) in coding; Only one SNP at position 544 within the coding sequence was obtained in the Central Asian wild barley population. Two SNPs were identified (position 507 and 544) from landraces, with one in a coding region. The haplotype 2 of NAM-1 was associated with the highest GPC, while the haplotype 7 of NAM-1 with lowest GPC (Figure 4).
TABLE 4. The single nucleotide polymorphism (SNP) and its positions relative to the reference sequence among distinct populations.
FIGURE 4. Boxplot of grain protein content (GPC) variation among 176 barley accessions grouped according to the haplotype of Hap2 and Hap7. Lines across the boxes depict the medians. Boxes indicate the interquartile range. Whiskers represent 95% confidence intervals.
Tibetan Wild Barley Diverged from Southwest Asian Barley
Previous studies have reported clear genetic differentiation between oriental and occidental barleys (Kahler and Allard, 1981; Zhang et al., 1992a,b; Ma, 2002; Komatsuda et al., 2004; Azhaguvel and Komatsuda, 2007). Ma (2002) showed morphological difference between two-rowed wild barley of the Qinghai–Tibet Plateau and the Middle Eastern barley. Genome-wide diversity analysis has shown significant genetic differences between wild barley from the Near East and Tibet (Dai et al., 2012, 2014). Our NAM-1 gene data revealed significant genetic differentiation among wild populations. Most dramatic differences in haplotype composition in wild barley occurred between the Tibetan and Southwest Asian barley, three haplotypes specific to the Tibetan wild barley population, and four haplotypes specific to the Southwest Asian wild barley population were detected (Table 2). No population-specific haplotype was found in wild barley from Central Asia. A higher genetic diversity and population-specific haplotype observed in the Tibetan wild barley, further supports that Tibetan wild barley is distinctly diverged from the Southwest Asian barley.
The Origin of Cultivated Barley
It has been well recognized that the Near East Fertile Crescent (Harlan and Zohary, 1966; Badr et al., 2000; Zohary et al., 2012) and Central Asia (Morrell and Clegg, 2007; Saisho and Purugganan, 2007) are the primary evolutionary centers of wild barley as well as a domestication center of cultivated forms. However, increasing evidence suggested that Tibet of China is an additional domestication center of cultivated barley (Dai et al., 2012, 2014; Nevo, 2013). The present study not only supports the status of the Fertile Crescent in domestication of cultivated barley, but also reveals Tibet as one of the centers of domestication of cultivated barley, thus further supporting the concept of polyphyletic domestication of barley. Hypotheses concerning the origin of barley have suggested that the varieties growing in the original center generally contain large amounts of dominant genetic genes (Wang et al., 2009). The region with the highest level of genetic diversity in wild barley is also most likely the center of origin for cultivated barley (Wang et al., 2009). In our study, high levels of nucleotide diversity, haplotype diversity and number of haplotypes were detected in the Tibetan and Southwest Asian wild barley populations (Table 3). Furthermore, the haplotype analysis showed that the worldwide cultivars shared the same haplotypes with the Southwest Asia wild barley, and specially, the wild barley from the Tibet (Table 2; Figure 1). A close relationship between worldwide domesticated barley and the Tibetan wild barley was revealed in our study, suggesting that Tibetan wild barley is one of the ancestors of domesticated barley. Our results corroborated previous finding that cultivated barley is not only derived from wild-barley genotypes in the Fertile Crescent, but also from those in Tibet of China (Dai et al., 2014).
Where Chinese cultivated barley originated from still remains to be addressed. At present, two hypotheses have been suggested. One suggested the Chinese cultivated barley was introduced from the Near East (Harlan, 1971, 1976; Badr et al., 2000), and another indicated that Chinese cultivated barley might have originated from two-rowed or six-rowed wild barley from Tibet (Zhang et al., 1994; Feng et al., 2006). Recent genome-wide diversity array data suggested that Chinese hulless and six-rowed barleys were domesticated in the Tibetan Plateau and its vicinity (Dai et al., 2012). Our results showed that 29 out of 61 landraces shared the same haplotype (Hap7) with the Tibetan wild barley, and phylogenetic analysis also revealed a close relationship between them. Haplotype Hap2 in the rest of the Chinese landraces was not only present in the Tibetan wild barley population, but also in the Central Asian and Southwest Asian wild barley populations (Table 2). Thus, our results not only supported that the Tibetan wild barley is the ancestor of Chinese domesticated barley, but also suggested that the Near East Fertile Crescent wild barley might have contributed to the origin of Chinese cultivars. This is in agreement with previous findings that landraces with majority western ancestry were relatively commonly encountered among Asian samples (Morrell et al., 2013) and the Oriental landraces have high proportion of admixed ancestry (Morrell and Clegg, 2007).
Implications for the Spread of Barley Cultivation in the World
Various hypotheses about the world spread of domesticated barley have been proposed. Badr et al. (2000) suggested that the border region between Israel and Jordan might be the region where barley was brought into cultivation and subsequently migrated to the area of the Himalayas. Some studies argued that barley was domesticated in this region and subsequently expanded westward into Europe and North Africa and eastward into Asia 8000 years ago (Von Bothmer et al., 2003). Morrell and Clegg (2007) proposed that the Fertile Crescent domestication contributed the majority of diversity in European and American cultivars, whereas the second domestication, 1500–3000 km farther east contributed most of the diversity in barley from Central Asia to the Far East. The trade of barley between the New World and Europe was supported by eIF4E gene data (Hofinger et al., 2011).
Our study provided interesting insights into historic global cultural / trade routes of barley. First, a haplotype that is private to the Southwest Asian wild barley population was also detected in the North American and European landrace barley populations (Table 2; Figure 1), corroborating assumptions made by Morrell and Clegg (2007) that Fertile Crescent domestication contributed the majority of diversity in European and American cultivars. In addition, we were surprised to see a haplotype that is exclusively found in Tibetan wild barley population is pervasive in all landrace populations (Table 2; Figure 1). It seems likely that the ancestral carrier(s) of this haplotype was initially introduced from the Tibet region to other geographic regions, which might explain the high levels of similarity between Eastern malting barley and European cultivars reported by Ordon et al. (1997). Our results suggested that the gene pool of Tibetan wild barley has been widely circulated, and has significantly contributed to the gene pool of global cultivated barley. Moreover, it may be assumed that Central Asia is the sole route for wild barley migration between the Near East and the Qinghai–Tibetan Plateau (Dai et al., 2012), as deduced in our study from haplotype Hap2 that was widespread and frequently found in Southwest Asia, Central Asia, and Tibetan wild barley populations, and worldwide cultivars (Table 2; Figure 1). Thus, our results supported the most likely scenario that the gene pool of the cultivated barley includes contributions of wild barleys from both the Near East and Tibet (Dai et al., 2012, 2014). Meanwhile, we suggest that the gene flow between Eastern and Western cultivars has occurred via the Silk Road, which started from China and moved westward, through the Eurasian civilization zones, Central Asia, and the Roman empire to Europe (Ma, 1998). The Silk Road might be an important barley transition route between the Orient and the Occident as previously proposed (Harold, 2007; Dai et al., 2012).
The spread of agriculture from domestication region involved the dispersal of crop plants well beyond their progenitors’ native range, and may have required adaptation to new environments (Jones et al., 2008). Ecologically, Tibetan wild barley is adapted to cold and dry environments, these characteristics may also be an important reason for its successfully spread all over the world (Dai et al., 2012).
Natural Variation in Barley Population
Crop domestication is the outcome of complex independent or combined processes of artificial and natural selection that lead to plants adapted to cultivation and to meet the requirements of human consumption (Dai et al., 2014). Gene pools undergoing domestication experienced dramatic changes in allele frequencies due to genetic drift or selection, and some allelic combinations may be lost (Wang et al., 2014). In our studies, a total of 10 distinct haplotypes were discovered, only 3 haplotypes were detected in the diverse set of 120 domesticated barleys from across the world, while more haplotypes occurred in wild barley accessions. This result agreed with previous reports (Kilian et al., 2006; Jakob et al., 2014) and further confirmed that most alleles in wild types have been lost in the domesticated forms. In addition, reduction in haplotype diversity, nucleotide diversity, and pre-site nucleotide diversity in domesticated lines was in accord with previous findings that H. spontaneum has a higher genetic diversity than H. vulgare landraces (Russell et al., 2004; Jin et al., 2011; Jakob et al., 2014), which might be caused by genetic bottlenecks acting on neutrally evolving loci either during the domestication process or during subsequent breeding, or both (Badr and El-Shazly, 2012). A shift toward more positive values of Tajima’s D in the domesticated relative to wild populations is indicative of reduced genetic diversity in the domesticated forms (Hufford et al., 2012; Morrell et al., 2013). A similar pattern was observed in our study, where Tajima’s D values are negative in wild barley populations of Southwest Asia, while they are more positive in domesticated barley populations of East Asia, South America, Mediterranean Coast Areas, and Europe. This similar shift from negative Tajima’s D at the majority of loci in the wild toward positive values in landraces, was also mentioned by Morrell et al. (2013). However, a positive Tajima’s D observed in Tibetan wild barley populations (Table 3) might be due to the fact that a rare-allele advantage resulted in an accumulating allelic frequency up to an intermediate level may have been involved in balancing selection, thus causing a positive value of Tajima’s D, as suggested by Chung et al. (2010). In general, deviation from neutrality with Tajima’s D was not significant (at P > 0.05) for any barley populations in our study. Domesticated barley populations of North America and Australia showed high negative Tajima’s D values (Table 3), suggesting that purifying selection might act on these populations. Insignificance may be attributed to the low number of SNPs (Table 4) observed, which weakens the neutrality test (Xia et al., 2013).
Association between GPC and NAM-1 Gene
GPC, as a key factor for quality in cereals, is influenced to a large extent by both genotype and environment (Smith, 1990; Jaradat, 1991). In our study, GPC was significantly different between domesticated and the wild barley population of Southwest Asia and Central Asia. Our findings support the studies that showed H. spontaneum with higher GPC values than cultivated barley (Jaradat, 1991; Jamar et al., 2010). However, wild barley population of Tibet had a lower GPC, which contradicts the study by Cai et al. (2013).
GPC in barley is influenced by both genetic and environmental factors (Bertholdsson, 1999; Distelfeld et al., 2008; Jamar et al., 2010; Cai et al., 2013). That allelic variation of the NAM-1 gene is an important genetic factor was demonstrated by Distelfeld et al. (2008), who found that two amino acid (aa) substitutions in HvNAM-1 might be associated with the GPC in barley. Jamar et al. (2010) suggested that allelic variation of NAM-1 gene might be associated with GPC variation in the genus Hordeum. Recent GWAS showed a significant correlation between haplotypes of HvNAM1, HvNAM2, and GPC in barley (Cai et al., 2013). In our study, two unique haplotypes (Hap2 and Hap7) might have a significant impact on the GPC (Table 4; Figure 4; Supplementary Data Sheet S2). The SNP at position 544 is within the coding region and causes a non-synonymous change with aa substitution between Alanine (A) and Proline (P). This substitution occurred in the C subdomain of NAC domain in the N-terminal, and may have an impact on protein folding (Jamar et al., 2010). The SNP at position 544 was also identified by Cai et al. (2013), in which the 102th aa proline (P) was replaced by alanine (A), resulting in low GPC. Therefore, we suggested that aa substitution in Hap2 and Hap7 might have an impact on DNA-binding ability of NAM-1 gene, and further affect GPC. In addition, the differences in GPC of a particular group with no polymorphisms at the NAM-1 gene might suggest that expression of the NAM-1 gene or other genes are also important (Jamar et al., 2010).
In summary, our results showed significant genetic differentiation among wild populations. Our data supported that Tibet is a center of origin and domestication centre for cultivated barleys, and suggested that the Silk Road might have played an important role in gene flow between Eastern and Western barley. Moreover, SNPs and haplotypes of NAM-1 associated with GPC in barley could provide a useful method for screening GPC in barley germplasm. The Tibetan wild accessions with lower GPC could be useful for malt barley breeding.
Conceived and designed the experiments: DS, GS. Performed the experiments: YW, XR. Analyzed the data: YW, XR, DS, and GS. Contributed reagents/materials/analysis tools: YW, XR, DS, and GS. Wrote the paper: YW, XR, DS, and GS.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This project was supported in part by the National Natural Science Foundation of China (Grant Nos. 31301310 and 31228017), an earmarked fund for China Agriculture Research System (CARS-5) and PhD Candidate Research Innovation Project of Huazhong Agricultural University (Program No. 2014bs07).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2015.00803
Azhaguvel, P., and Komatsuda, T. (2007). A phylogenetic analysis based on nucleotide sequence of a marker linked to the brittle rachis locus indicates a diphyletic origin of barley. Ann. Bot. 100, 1009–1015. doi: 10.1093/aob/mcm129
Badr, A., and El-Shazly, H. (2012). Molecular approaches to origin, ancestry and domestication history of crop plants: barley and clover as examples. J. Genet. Eng. Biotechnol. 10, 1–12. doi: 10.1016/j.jgeb.2011.08.002
Badr, A., Müller, K., Schäfer-Pregl, R., El Rabey, H., Effgen, S., Ibrahim, H., et al. (2000). On the origin and domestication history of barley (Hordeum vulgare). Mol. Biol. Evol. 17, 499–510. doi: 10.1093/oxfordjournals.molbev.a026330
Balazadeh, S., Siddiqui, H., Allu, A. D., Matallana-Ramirez, L. P., Caldana, C., Mehrnia, M., et al. (2010). A gene regulatory network controlled by the NAC transcription factor ANAC092/AtNAC2/ORE1 during salt-promoted senescence. Plant J. 62, 250–264. doi: 10.1111/j.1365-313X.2010.04151.x
Bekele, E. (1983). A differential rate of regional distribution of barley flavonoid patterns in Ethiopia, and a view on the center of origin of barley. Hereditas 98, 269–280. doi: 10.1111/j.1601-5223.1983.tb00605.x
Bertholdsson, N. (1999). Characterization of malting barley cultivars with more or less stable grain protein content under varying environmental conditions. Eur. J. Agron. 10, 1–8. doi: 10.1016/S1161-0301(98)00043-4
Blanco, A., Mangini, G., Giancaspro, A., Giove, S., Colasuonno, P., Simeone, R., et al. (2012). Relationships between grain protein content and grain yield components through quantitative trait locus analyses in a recombinant inbred line population derived from two elite durum wheat cultivars. Mol. Breed. 30, 79–92. doi: 10.1007/s11032-011-9600-z
Bogard, M., Allard, V., Brancourt-Hulmel, M., Heumez, E., Machet, J. M., Jeuffroy, M. H., et al. (2010). Deviation from the grain protein concentration–grain yield negative relationship is highly correlated to post-anthesis N uptake in winter wheat. J. Exp. Bot. 61, 4303–4312. doi: 10.1093/jxb/erq238
Bogard, M., Jourdan, M., Allard, V., Martre, P., Perretant, M. R., Ravel, C., et al. (2011). Anthesis date mainly explained correlations between post-anthesis leaf senescence, grain yield, and grain protein concentration in a winter wheat population segregating for flowering time QTLs. J. Exp. Bot. 62, 3621–3636. doi: 10.1093/jxb/err061
Chung, J. D., Lin, T. P., Cheng, Y. P., and Hwang, S. Y. (2010). Geographical variation and differential selection modes of paralogs of chloroplast small heat shock protein genes in Machilus kusanoi (Lauraceae). Bot. Stud. 51, 27–36.
Dai, F., Chen, Z. H., Wang, X., Li, Z., Jin, G., Wu, D., et al. (2014). Transcriptome profiling reveals mosaic genomic origins of modern cultivated barley. Proc. Natl. Acad. Sci. U.S.A. 111, 13403–13408. doi: 10.1073/pnas.1414335111
Dai, F., Nevo, E., Wu, D., Comadran, J., Zhou, M., Qiu, L., et al. (2012). Tibet is one of the centers of domestication of cultivated barley. Proc. Natl. Acad. Sci. U.S.A. 109, 16969–16973. doi: 10.1073/pnas.1215265109
Distelfeld, A., Korol, A., Dubcovsky, J., Uauy, C., Blake, T., and Fahima, T. (2008). Colinearity between the barley grain protein content (GPC) QTL on chromosome arm 6HS and the wheat Gpc-B1 region. Mol. Breed. 22, 25–38. doi: 10.1007/s11032-007-9153-3
Fu, Y. B., and Peterson, G. W. (2011). Genetic diversity analysis with 454 pyrosequencing and genomic reduction confirmed the eastern and western division in the cultivated barley gene pool. Plant Genome 4, 226–237. doi: 10.3835/plantgenome2011.08.0022
Fuller, D. Q., Willcox, G., and Allaby, R. G. (2011). Cultivation and domestication had multiple origins: arguments against the core area hypothesis for the origins of agriculture in the Near East. World Archaeol. 43, 628–652. doi: 10.1080/00438243.2011.624747
Hofinger, B. J., Russell, J. R., Bass, C. G., Baldwin, T., Dos Reis, M., Hedley, P. E., et al. (2011). An exceptionally high nucleotide and haplotype diversity and a signature of positive selection for the eIF4E resistance gene in barley are revealed by allele mining and phylogenetic analyses of natural populations. Mol. Ecol. 20, 3653–3668. doi: 10.1111/j.1365-294X.2011.05201.x
Hufford, M. B., Xu, X., Van Heerwaarden, J., Pyhäjärvi, T., Chia, J. M., Cartwright, R. A., et al. (2012). Comparative population genomics of maize domestication and improvement. Nat. Genet. 44, 808–811. doi: 10.1038/ng.2309
Jakob, S. S., Rödder, D., Engler, J. O., Shaaf, S., Özkan, H., Blattner, F. R., et al. (2014). Evolutionary history of wild barley (Hordeum vulgare subsp. spontaneum) analyzed using multilocus sequence data and paleodistribution modeling. Genome Biol. Evol. 6, 685–702. doi: 10.1093/gbe/evu047
Jamar, C., Loffet, F., Frettinger, P., Ramsay, L., Fauconnier, M. L., and Du Jardin, P. (2010). NAM-1gene polymorphism and grain protein content in Hordeum. J. Plant Physiol. 167, 497–501. doi: 10.1016/j.jplph.2009.10.014
Jensen, M. K., and Skriver, K. (2014). NAC transcription factor gene regulatory and protein–protein interaction networks in plant stress responses and senescence. IUBMB Life 66, 156–166. doi: 10.1002/iub.1256
Jin, X., Cai, S., Han, Y., Wang, J., Wei, K., and Zhang, G. (2011). Genetic variants of HvGlb1 in Tibetan annual wild barley and cultivated barley and their correlation with malt quality. J. Cereal Sci. 53, 59–64. doi: 10.1016/j.jcs.2010.09.006
Jones, H., Leigh, F. J., Mackay, I., Bower, M. A., Smith, L. M., Charles, M. P., et al. (2008). Population-based resequencing reveals that the flowering time adaptation of cultivated barley originated east of the Fertile Crescent. Mol. Biol. Evol. 25, 2211–2219. doi: 10.1093/molbev/msn167
Jukanti, A. K., and Fischer, A. M. (2008). A high-grain protein content locus on barley (Hordeum vulgare) chromosome 6 is associated with increased flag leaf proteolysis and nitrogen remobilization. Physiol. Plant. 132, 426–439. doi: 10.1111/j.1399-3054.2007.01044.x
Kilian, B., Özkan, H., Kohl, J., Von Haeseler, A., Barale, F., Deusch, O., et al. (2006). Haplotype structure at seven barley genes: relevance to gene pool bottlenecks, phylogeny of ear type and site of barley domestication. Mol. Genet. Genomics 276, 230–241. doi: 10.1007/s00438-006-0136-6
Komatsuda, T., Maxim, P., Senthil, N., and Mano, Y. (2004). High-density AFLP map of nonbrittle rachis 1 (btr1) and 2 (btr2) genes in barley (Hordeum vulgare L.). Theor. Appl. Genet. 109, 986–995. doi: 10.1007/s00122-004-1710-0
Ma, D. (1998). “The great silk exchange: how the world was connected and developed,” in Pacific Centuries: Pacific and Pacific Rim Economic History Since the Sixteenth Century, eds D. O. Flynn, L. Frost, and A. J. H. Latham (New York, NY: Routledge Press), 38–69.
Martre, P., He, J., Le Gouis, J., and Semenov, M. A. (2015). In silico system analysis of physiological traits determining grain yield and protein concentration for wheat as influenced by climate and crop management. J. Exp. Bot. 66, 3581–3598. doi: 10.1093/jxb/erv049
Molina-Cano, J., Fra-Mon, P., Salcedo, G., Aragoncillo, C., De Togores, F. R., and García-Olmedo, F. (1987). Morocco as a possible domestication center for barley: biochemical and agromorphological evidence. Theor. Appl. Genet. 73, 531–536. doi: 10.1007/BF00289190
Molina-Cano, J., Igartua, E., Casas, A., and Moralejo, M. (2002). “New views on the origin of cultivated barley,” in Barley Science: Recent Advances from Molecular Biology to Agronomy of Yield and Quality, eds G. A. Slafer, J. L. Molina-Cano, R. Savin, J. L. Araus, and I. Romagosa (Binghamton, NY: The Haworth Press), 15–29.
Morrell, P. L., and Clegg, M. T. (2007). Genetic evidence for a second domestication of barley (Hordeum vulgare) east of the Fertile Crescent. Proc. Natl. Acad. Sci. U.S.A. 104, 3289–3294. doi: 10.1073/pnas.0611377104
Morrell, P. L., Gonzales, A. M., Meyer, K. K., and Clegg, M. T. (2013). Resequencing data indicate a modest effect of domestication on diversity in barley: a cultigen with multiple origins. J. Hered. 105, 253–264. doi: 10.1093/jhered/est083
Nevo, E. (2013). “Evolution of wild barley and barley improvement,” in Advance in Barley Sciences, eds G. Zhang, C. Li, and X. Liu (Dordrecht: Zhejiang University Press and Springer Verlag), 1–23. doi: 10.1007/978-94-007-4682-4_1
Nevo, E., Beiles, A., Gutterman, Y., Storch, N., and Kaplan, D. (1984). Genetic resources of wild cereals in Israel and vicinity. II. Phenotypic variation within and between populations of wild barley, Hordeum spontaneum. Euphytica 33, 737–756. doi: 10.1007/BF00021901
Nevo, E., Beiles, A., and Zohary, D. (1986). Genetic resources of wild barley in the Near East: structure, evolution and application in breeding. Biol. J. Linn. Soc. 27, 355–380. doi: 10.1111/j.1095-8312.1986.tb01742.x
Ordon, F., Schiemann, A., and Friedt, W. (1997). Assessment of the genetic relatedness of barley accessions (Hordeum vulgare s.l) resistant to soil-borne mosaic-inducing viruses (BaMMV, BaYMV, BaYMV-2) using RAPDs. Theor. Appl. Genet. 94, 325–330. doi: 10.1007/s001220050419
Oury, F. X., and Godin, C. (2007). Yield and grain protein concentration in bread wheat: how to use the negative relationship between the two characters to identify favourable genotypes? Euphytica 157, 45–57. doi: 10.1007/s10681-007-9395-5
Pakniyat, H., Powell, W., Baird, E., Handley, L., Robinson, D., Scrimgeour, C., et al. (1997). AFLP variation in wild barley (Hordeum spontaneum C. Koch) with reference to salt tolerance and associated ecogeography. Genome 40, 332–341. doi: 10.1139/g97-046
Russell, J., Booth, A., Fuller, J., Harrower, B., Hedley, P., Machray, G., et al. (2004). A comparison of sequence-based polymorphism and haplotype content in transcribed and anonymous regions of the barley genome. Genome 47, 389–398. doi: 10.1139/g03-125
Stein, N., Herren, G., and Keller, B. (2001). A new DNA extraction method for high-throughput marker analysis in a large-genome species such as Triticum aestivum. Plant Breed. 120, 354–356. doi: 10.1046/j.1439-0523.2001.00615.x
Thompson, J. D., Gibson, T. J., Plewniak, F., Jeanmougin, F., and Higgins, D. G. (1997). The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 25, 4876–4882. doi: 10.1093/nar/25.24.4876
Uauy, C., Distelfeld, A., Fahima, T., Blechl, A., and Dubcovsky, J. (2006). A NAC gene regulating senescence improves grain protein, zinc, and iron content in wheat. Science 314, 1298–1301. doi: 10.1126/science.1133649
Von Bothmer, R., Jacobsen, N., Baden, C., Jørgensen, R., and Linde-Laursen, I. (1995). An Ecogeographical Study of the Genus Hordeum. Systematic and Ecogeographic Studies on Crop Genepools 7. Rome: International Plant Genetic Resources Institute.
Von Bothmer, R., and Komatsuda, T. (2011). “Barley origin and related species,” in Barley: Production, Improvement and Uses, ed. S. E. Ullrich (Chichester: Wiley-Blackwell Press), 14–62. doi: 10.1002/9780470958636.ch2
Von Bothmer, R., Sato, K., Kamatsuda, T., Yasuda, S., and Fischbeck, G. (2003). “The domestication of cultivated barley,” in Diversity in Barley (Hordeum vulgare), eds R. Von Bothmer, T. van Hintum, H. Knüpffer, and K. Sato (Amsterdam: Elsevier Press), 9–27. doi: 10.1016/S0168-7972(03)80004-X
Wang, A., Yu, Z., and Ding, Y. (2009). Genetic diversity analysis of wild close relatives of barley from tibet and the middle east by ISSR and SSR markers. C. R. Biol. 332, 393–403. doi: 10.1016/j.crvi.2008.11.007
Wang, L., Xu, J., Xia, T., Zhang, H., Liu, D., and Shen, Y. (2014). Population structure and linkage disequilibrium in six-rowed barley landraces from the Qinghai-Tibetan plateau. Crop Sci. 54, 2011–2022. doi: 10.2135/cropsci2013.08.0517
Xia, Y., Li, R., Ning, Z., Bai, G., Siddique, K. H., Yan, G., et al. (2013). Single nucleotide polymorphisms in HSP17. 8 and their association with agronomic traits in barley. PLoS ONE 8:e56816. doi: 10.1371/journal.pone.0056816
Zhang, Q., Dai, X., and Saghai Maroof, M. (1992a). Comparative assessment of genetic variation at 6 isozyme loci in barley from two centers of diversity: ethiopia and Tibet. Acta Genet. Sin. 19, 236–243.
Zhang, Q., Maroof, M. S., and Yang, P. (1992b). Ribosomal DNA polymorphisms and the oriental-occidental genetic differentiation in cultivated barley. Theor. Appl. Genet. 84, 682–687. doi: 10.1007/BF00224168
Zhang, Q., Yang, G., Dai, X., and Sun, J. (1994). A comparative analysis of genetic polymorphism in wild and cultivated barley from Tibet using isozyme and ribosomal DNA markers. Genome 37, 631–638. doi: 10.1139/g94-090
Zohary, D., Hopf, M., and Weiss, E. (2012). Domestication of Plants in the Old World: The Origin and Spread of Domesticated Plants in Southwest Asia, Europe, and the Mediterranean Basin. Oxford: Oxford University Press.
Keywords: origin, spread, barley, NAM-1 gene, grain protein content
Citation: Wang Y, Ren X, Sun D and Sun G (2015) Origin of worldwide cultivated barley revealed by NAM-1 gene and grain protein content. Front. Plant Sci. 6:803. doi: 10.3389/fpls.2015.00803
Received: 12 July 2015; Accepted: 15 September 2015;
Published: 30 September 2015.
Edited by:James J. Cai, Texas A&M University, USA
Reviewed by:Tina T. Hu, Princeton University, USA
Stefanie M. Ickert-Bond, University of Alaska Fairbanks, USA
Chuanchao Wang, Max Planck Institute for the Science of Human History, Germany
Lingxiang Wang, Fudan University, China
Junhua Peng, China National Seed Group Co. Ltd., China
Copyright © 2015 Wang, Ren, Sun and Sun. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dongfa Sun, College of Plant Science and Technology, Huazhong Agricultural University, Wuhan 430070, China, email@example.com; Genlou Sun, Department of Biology, Saint Mary’s University, 923 Robie Street, Halifax, NS B3H 3C3, Canada, firstname.lastname@example.org