New insights into the evolution of wheat avenin-like proteins in wild emmer wheat (Triticum dicoccoides).

Fifteen full-length wheat grain avenin-like protein coding genes (TaALP) were identified on chromosome arms 7AS, 4AL, and 7DS of bread wheat with each containing five genes. Besides the a- and b-type ALPs, a c type was identified in the current paper. Both a and b types have two subunits, named x and y types. The five genes on each of the three chromosome arms consisted of two x-type genes, two y-type genes, and one c-type gene. The a-type genes were typically of 520 bp in length, whereas the b types were of 850 bp in length, and the c type was of 470 bp in length. The ALP gene transcript levels were significantly up-regulated in Blumeria graminis f. sp. tritici (Bgt)-infected wheat grain caryopsis at early grain filling. Wild emmer wheat [(WEW), Triticum dicoccoides] populations were focused on in our paper to identify allelic variations of ALP genes and to study the influence of natural selection on certain alleles. Consequently, 25 alleles were identified for TdALP-bx-7AS, 13 alleles were identified for TdALP-ax-7AS, 7 alleles were identified for TdALP-ay-7AS, and 4 alleles were identified for TdALP-ax-4AL Correlation studies on TdALP gene diversity and ecological stresses suggested that environmental factors contribute to the ALP polymorphism formation in WEW. Many allelic variants of ALPs in the endosperm of WEW are not present in bread wheat and therefore could be utilized in breeding bread wheat varieties for better quality and elite plant defense characteristics.

Prolamin superfamily proteins share a conserved pattern of cysteine residues, including the sulfur-rich prolamins of the Triticeae, the cereal α-amylase/trypsin inhibitors, 2S storage albumins, puroindolines, grain softness proteins, α-globulins, and a group of hydroxyproline-rich cell wall proteins, which might all have originated from a small number of ancestral genes. According to Shewry and Halford (1), the gliadins, members of the prolamin superfamily, include members with a large repetitive domain and a conserved set of cysteine residues (α- and γ-gliadins), members with a repetitive domain but no cysteine (ω-gliadins), and members with novel low molecular weight gliadins (LMWGs) also known as avenin-like proteins (ALPs) that contain a conserved cysteine pattern but with no repetitive domains (2).

LMWGs are proteo-lipid-like hydrophobic proteins, similar to albumin-like and globulin-like proteins (3, 4). Genes encoding LMWGs are located in bread wheat on chromosomes 7A, 4A, and 7D (3). This observation supports the 4A/7B chromosome interchange hypothesis because there is a similar chromosomal distribution of peroxidase genes (5–7). In 2001, Anderson and others (8) cloned five genes that shared complex relationships with the gliadins. One cloned gene, 11dc7, corresponded to one group of LMWGs described by Salcedo et al. (4). Rocher et al. (9) reported two similar LMWG proteins, rye-15 and rye-18 that showed weak immune reactivity with antibodies in serum from celiac patients. Clark et al. described the identification of a functional class of genes relevant to wheat grain end use belonging to a novel glutenin/gliadin seed storage protein (10). Kan et al. (11) identified two highly expressed transcripts encoding a- and b-type ALPs in wheat but with typically much higher expression in the Aegilops species. Overexpression of type-b ALPs in transgenic wheat improved dough mixing properties (12). ALP-coding genes were mapped to the short arms of chromosomes 7A and 7D and to the long arm of chromosome 4A in bread wheat (13). Importantly, alleles on 7A have been found with differential effects on dough quality, and its allele-specific markers have been developed to track the allelic effects (13). Recently, wheat ALP proteins were discovered with a significant Fusarium head blight resistant function, which highlighted the divergent functions of this gliadin domain containing protein family (14).

Wild emmer wheat (WEW), T. dicoccoides, is the progenitor of cultivated tetraploid and hexaploid wheats. It evolved in the northern ecogeographical region of the upper Jordan River in the eastern Upper Galilee Mountains and Golan Heights. Here, we studied 21 WEW populations from across their natural range in Israel. These were screened for allelic variation of ALP genes with the aim of identifying alleles useful for bread wheat improvement and determining the regional ecological influences on allele formation. The 21 Israeli populations used in this paper had previously been studied by Nevo and coworkers (15–22). They identified local and regional ecological differences, genetic differences, and allozymic polymorphisms. These early studies identified adaptive allozyme diversity induced by abiotic and biotic stresses, highlighting the influence of selection on the adaptive nature of allozymic variation, and thereby negating the neutral theory of evolution.

Results

Identification of ALP Homologous Genes from Wheat Genome Phylogeny.

We start the results section with an analysis of ALP homologs from wheat genomes to place the TdALP analysis from WEW populations () in a broader context of the Triticeae and other plants. In allohexaploid bread wheat, 15 unique full-length TaALP genes cDNA were mapped to chromosome groups 4 and 7 (). Besides ALP genes reported in published studies (13, 23, 24), all other genes were cloned and sequenced (). Alignment of the translated amino acid sequences encoded by the 15 full-length TaALP genes showed that ALPs vary in length from 150 to 285 amino acids (). Their signal peptides were predicted and listed in . According to the domain classification based on the pFam database, ALPs are characterized by possessing gliadin domains (PF13016) as well as α-amylase inhibitors and seed storage protein subfamily domains (PF00234) (25). Based on alignment analysis and comparison with the reported a- and b-type TaALP genes, a new type was found and named c type in this paper (Fig. 1 and ). The homogeneity of the TaALP genes within the same subgroup is based on their high sequence identities (>86.43%) (). To investigate the evolutionary relationships among ALP genes from the Triticeae, a ML phylogeny was constructed based on the deduced amino acid sequence alignment of 46 genes, including ALP-related sequences. As shown (Fig. 1), besides the outgroup of puroindolines and avenin-3 in the monophyletic group, three major ALP gene clades, type a (blue), type b (red), and type c (green) for the 15 genes in bread wheat as well as one gene copy of B. distachyon were classified. Type-a ALP and type-b ALPs can be further divided into x and y subgroups. For the type-c clade, three TaALP genes on chromosomes 4A, 7A, and 7D from bread wheat and one orthologous barley gene HvALP on chromosome 7H, one from T. dicoccoides chromosome 7A, and genes from T. monoccocum, T. urartu, A. speltoides, A. sharonesis, A. tauchii, and B. distachyon were all closely related. Three subunits from T. urartu are clustered, indicating three type-c ALP genes in that species. Within the type-b clade y-type subgroup, one ortholog for T. monococcum and two homeologous genes on chromosomes 4A and 7A for T. dicoccoides were found. However, other orthologs of type-b ALPs were not identified in databases. For the x-type subgroup, one ortholog from A. taucchii and three homeologous TaALP genes from bread wheat genomes were identified. For the type-a clade y-type subgroups, 11 genes were identified—three from bread wheat, two from T. dicoccoides (4A and 7A), one ortholog from barley (7H), and one each from T. monoccocum, T. urartu, A. speltoides, A. sharonesis, and A. tauchii. The same patterns were found for the type-a ALP clade x-type subunits. All of the gene sequences used in this phylogenetic analysis () contain a gliadin domain classified as PF13016. An unrooted ML phylogenetic analysis identified a wheat ALP clade (gray), avenin-3, a gliadin clade (purple), and a clade comprising millet, sorghum, and maize prolamins that have a gliadin domain (yellow). The TaALP clade was separated into three subgroups comprising type a, type b, and type c, whereas the rice prolamin, α-amylase inhibitor, grain softness protein, and puroindolines diverged from ALPs much earlier in the evolutionary history. This most recent common ancestor of the wheat and barley ALP clade (gray) and avenin-3 and gliadin clades (purple) can be traced to much earlier in the PF13016 domain evolutionary history (). Much earlier are the common ancestors of the members of the yellow clade for the millet, sorghum, and maize prolamins with a gliadin domain. As a result, Triticeae prolamins (ALPs, gliadins, and avenin-3) are more closely related to one another than to Panicoideae prolamins. A neighbor joining (NJ) analysis was performed on gliadin domains, AAI-LTSS domains, and LTP2 domains of monocots and lower plants ().

Fig. 1.

Phylogenetic analyses of the ALPs. Maximum likelihood (ML) phylogeny of ALPs of the bread wheat (Triricum aestivum), T. dicoccoides, Triticum urartu, Triricum monococcum, Aegilops speltoides, Aegilops. sharonesis, Aegilops tauschii, Brachypodium distachyon, and Hordeum vulgare based on amino acid sequence alignment.

Transcriptional Analyses of TaALP Genes in Bread Wheat Under Bgt Infection.

Gene expression dynamics of TaALP genes under biotic stress were studied to select gene loci for detailed evolutionary study. The relative expression of TaALP genes in the lemma and grain of Spitfire × Mace doubled haploid (DH) lines 130, 131, and 187 were studied at 2 and 10 d after pollination (DAP) under Bgt infection (). Gene-specific primers were designed for the 15 TaALP genes (). The parent wheat cvs. Spitfire and Mace displayed allelic variations at three TaALP loci, whereas the four DH lines were selected with consistent allelic compositions for the TaALP genes (). The healthy wheat lines (DH line 241) were chosen as the control, and their ALP expressions at 2, 7, and 10 DAP were shown in . The TaALP genes were up-regulated in Bgt affected DH lines 130, 131, and 187 compared with the DH line 241 (), and relative expression of diverse types of TaALP genes showed significant positive linear correlations under Bgt infection (). Based on the transcriptional study above, we selected four TaALP genes (bx/ay/ax-7AS and -ax-4AL) for allele screening and evolution study across the 21 WEW populations.

Gene Cloning and Sequencing Analyses of Four Selected TdALP Genes in WEW.

Cloning and sequencing of the four selected TdALP genes in WEW (TdALP-bx/ay/ax-7AS and -ax-4AL) revealed a surprisingly rich diversity. A total of 49 alleles were identified, including 25 bx-7AS genes, 13 ax-7AS genes, 7 ay-7AS genes, and 4 ax-4AL genes (). For the TdALP-bx-7AS gene, among the 25 haplotypes (), 14 bx-7AS genes were assumed to be pseudogenes. Other alleles, bx-7AS-a*, -d, -g, -k, -m, -r, -s, -t, -x, -y, and -z, were functional genes, and the amino acid translations indicated continuous reads from initiation to termination. Amino acid A/T replacement at position 12, a Q insertion at position 35, and I/S and M/W replacements at positions 58 and 60 occurred for genes bx-7AS-g, -s, -t, and -r. The Q/H replacement at position 205 also occurred in several alleles (bx-7AS-m, -g, -s, and -t). The Q insertion at position 35 occurred for alleles bx-7AS-d and -k. For bx-7AS-k encoded type-b ALP proteins, G/C replacement at the N-terminal region occurred. The coding sequences and deduced amino acid sequences of TdALP-ax/ay-7AS and ax-4AL were aligned, and the particular single-nucleotide polymorphisms (SNPs) and indels are shown in . The amino acid alignments of the 11 functional alleles of TdALP-bx-7AS are shown in , whereas the type-a alleles are shown in .

Population Genetics in Relation to Water and Edaphic Effects on TdALP Gene Diversity.

The genetic diversity among different populations of WEW in Israel was assessed by comparing the TdALP gene alleles identified from each population and the corresponding He and I indices. As shown in , 45 of 49 alleles were present in WEW populations in Israel. Overall, the four TdALP gene loci were polymorphic in most populations. The mean number of alleles per locus ranged from one to eight (). The genetic variation of TdALP genes displayed a clear region-specific pattern, corresponding to a He index ranging from 0 to 0.64 and an I index from 0 to 1.26 (Fig. 2).

Fig. 2.

Population genetics of WEW based on TdALP diversity. (A) Genetic indices of 21 WEW populations; (B) Spearman rank correlations of genetic indices at each WEW and climatic variables; (C) So and Rn are two environmental parameters which can best predict the TdALP diversity in WEW populations. Note: ** Correlation is significant at the 0.01 level (one tailed). * Correlation is significant at the 0.05 level (one tailed). He, Nei’s—1973 gene diversity; I, Shannon’s Information index (46); na, Observed number of alleles; ne, Effective number of alleles (45).

To investigate the association of genetic diversity of TdALP genes in WEW with environmental variables, first, one-tailed Spearman correlation analysis was performed to analyze the association of genetic diversity of TdALP (He and I) with various environmental variables (Fig. 2 and ). The results showed that eight variables, including altitude, mean annual temperature, mean August temperature, mean January temperature, mean annual rainfall (Rn), mean annual evaporation, interannual variability of rainfall and soil type (So), and mean interannual variability of rainfall, were significantly correlated with the TdALP genotype (Fig. 2). Second, the eight variables listed above plus two further geographical parameters were tested by backward multiple linear regression (MR) analysis for which 15 of the 21 WEW populations were included (). When the three-variable model was used (), a significant regression equation was found: F(3, 11) = 8.99, P < 0.01 with an R2 of 0.710, t(15) = 5.62, and P < 0.01. Accordingly, latitude (50%), Rn (99%), and So (56%) were identified as the variables that could explain the highest proportion of the TdALP genetic diversity (He) among different populations. Noteworthy is the standardized coefficient for the Rn variable calculated as 99%, which suggests that Rn dominates the other two variables in their relevance to TdALP diversity. As such, a two-variable model was applied to the dataset (). A significant regression equation was found: F(2, 12) = 10.69, P < 0.01 with an R2 of 0.640, t(15) = 5.90, and P < 0.001. The results showed that Rn and So could explain 65% and 35% of TdALP diversity (He), respectively, which suggest that the two-variable model fits the dataset much better than the three-variable model. To further validate that Rn and So are variables mostly associated with the TdALP genetic diversity, other backward MR analyses were performed. Similarly, a significant regression equation was found: F(2, 12) = 11.59, P < 0.01 with an R2 of 0.659, t(15) = 4.98, and P < 0.001 (). The results showed that Rn and So contribute 57% and 47%, respectively, to the genetic diversity (I) (Fig. 2). These results are comparable to the He index calculation (Fig. 2).

TdALP-bx-7AS Gene Clustering Analysis and Correlation with Environmental Factors.

A NJ analysis was performed based on 25 TdALP-bx-7AS gene sequence alignments (Fig. 3). Natural selection pressure on the bx-7AS gene in WEW was examined by measuring the ratio of nonsynonymous to synonymous substitutions (dN:dS = ω) (26). The silent alleles (−) and the functional alleles (+) are listed in . The branch ω value of Haplo 4, 9, 12, and 14 is 0.45, for Haplo 15 and 25 the branch ω value is 0.0001, for Haplo 1, 7, and 10 the branch ω value is 0.0001, for Haplo 2, 8, 17, and 23 the branch ω value is 2.37, for Haplo 18, 21, and 22 the branch ω value is 2.23, for Haplo 3, 5, 6, 10, 11, 13, 16, and 19 the branch ω value is 1.93, and for Haplo 24 the branch ω value is 1, displaying a neutral selection. The results indicated that all of the functional alleles (ω < 1 branches) are under purifying selection (except for Haplo 11, clustered with other silent alleles), whereas all of the silent alleles (ω > 1 branches) are under positive selection. The t test results (one tailed with equal variance) () showed that the environmental factors had a P value < 0.05, indicating significant correlations of the TdALP-bx-7AS functional allele/silent allele (±) with environmental factors. The microenvironments selecting functional alleles (+) were significantly different from the microenvironments favoring silent alleles (−). The P values of the mean number of rainy days, the mean number of dew nights in summer, and the mean annual evaporation were >0.05, indicating no significant correlations.

Fig. 3.

Phylogenetic analysis of the TdALP-bx-7AS gene in WEW populations. (A) NJ phylogenetic analysis and natural selection tests of 25 haplotypes of the TaALP-bx-7AS gene in WEW; (B) UPGMA phylogenetic analysis based on TaALP-bx-7AS gene variation in 21 WEW populations.

Unweighted Pair Group Method with Arithmetic Mean Phylogenetic Analysis of TdALP-bx/ay/ax-7AS in WEW Populations.

Genetic variation of TdALP-bx-7AS among different WEW populations was analyzed. An unweighted pair group method with arithmetic mean (UPGMA) tree was developed based on TdALP-bx/ax/ay-7AS sequence alignments (Fig. 3). In total, eight alleles (bx-7AS-f, -g, -h, -i, -j, -k, -l, and -m) were identified in the Yehudiyya population. These alleles were distributed among different branches, suggesting a significant degree of genetic variation. Some alleles are related to adaption to shade (-j and -h) under tree canopies and sun (-k and -l) between trees (17, 20–22). Some alleles collected from different WEW populations clustered together. For example, alleles bx-7AS-g (Yehudiyya, Tabigha terra rossa, and Mt. Gerizm), -s (Bet Oren, Daliyya, and Mt. Hermon), and -t (J’aba and Taiyiba) were grouped together. The bx-7AS-m allele was found in Yehudiyya (both sun and shade), Mt. Gerizm, and Amirim, whereas -n was present in Tabigha basalt and Ammiad-85. The bx-7AS-ps from Nesher, Tabigha terra rossa, Mt. Gilboa, Kokhav Hashahar, and Bet Meir (mostly xeric populations) were grouped together with -j from Yehudiyya shade, suggesting they were related populations from sites of similar ecogeographic backgrounds with respect to rainfall and terra rossa soil. The remaining 12 bx-7AS genes were from populations of diverse ecogeographic backgrounds, and these were clustered together in a separate group. Notably, bx-7AS-f dominates those populations from Yehudiyya (allele frequency of 28.22%, ) and ranks the highest in Tabigha (with both terra rossa and basalt soils). Specifically, taking So into consideration, alleles bx-7AS-g, -p, and -i were found in Tabigha terra rossa not in the abutting basalt soil. In contrast, -n was found in Tabigha basalt soil, demonstrating adaption to this So. The results of the UPGMA analysis for TdALP-ax-7AS and TdALP-ay-7AS genes are shown in , respectively.

Genetic Distance Analyses Among Different WEW Populations.

Pairwise genetic distances (p distance), based on the normalized gene sequence identity of TdALP-bx/ay/ax-7AS in WEW populations, were calculated using MEGA 7.0 software (). Overall, the target populations demonstrated a close distance with each other, ranging from 0.013 to 0.291 (). The genetic distance between WEW populations was found for J’aba, >0.195 with Bet Meir, Gamla, Gitit, Mt. Gilboa, Qazrin, Sanhedriyya, Tabigha 1979, Taiyiba, and Kokhav Hashahar. The maximum pairwise genetic distance (p distance = 0.291) among different populations was identified between populations J’aba and Kokhav Hashahar, indicating significant variations in closer but marginal steppic populations. However, one interesting phenomenon was observed, some physically distant populations displayed relatively lower p distances than some physically close populations. For example, populations Gamla and Yehudiyya sun, which were separated by only 9.5 km (), have a p distance of 0.212. In contrast, the p distance between Mt. Gerizim and Gamla (156 km, ) is 0.155, and, most significantly, Yehudiyya sun and J’aba (130 km, ) have the lowest p distance (0.013) (). Noteworthy, the p distance between Yehudiyya sun and Yehudiyya shade is 0.133, whereas that of Tabigha terra rossa and basalt is 0.163 (), suggesting that the So plays a more significant role than temperature.

Discussion

Origin, Mechanism, and Phylogeny of ALP Gene Evolution.

The ALPs' antifungal functions revealed by Zhang et al. (14) indicate a vital importance to identify genetic diversity of ALPs for potential exploration in wheat breeding. Similar to most grain storage protein genes, the current paper revealed that abundant ALP alleles have accumulated during evolution. Previous studies have shown that unequal crossover or gene slippage of insertions or deletions of blocks often happened during duplication events for HMW-GS genes (27, 28), which might also help to explain the emergence, expansion, and the allelic variations of the ALP genes. Domain replication observed for ALP genes (b type) serves a similar function to gene duplication, establishing gene variability driven by evolutionary forces. As for many protease inhibitor gene families, instead of complete gene duplication, including promoter and terminator sequences and possible reintegration at a distinct locus, there is duplication of the inhibitory domain sequence with the domains remaining fused (29–31). Fifteen full-length TaALP genes were clustered into three major subgroups (types a, b, and c) in our phylogenetic analyses (Fig. 1 and ). In addition, our results revealed the existence of intra- and interchromosomal ALP genes in WEW and bread wheat. There are three copies of type-c genes in T. urartu (Fig. 1), whereas only one copy in each chromosome of bread wheat and WEW, indicating that wheat-specific ALP gene duplication/elimination events most likely occurred in a diploid wheat ancestor, leading to the loss of these genes in WEW as well as in bread wheat. Similar wheat-specific gene duplication events and/or chromosomal translocations are also likely to be responsible for the origin of the multi-ALP genes.

The Importance of Natural Population in Highlighting Genetic Adaptations.

Experimental populations evolving under natural selection represent an important resource for studying the genetic basis of adaptation. Our analysis of TdALP gene variation and evolution in WEW was based on ecogenetic analyses of Israeli and Golan Heights populations as was demonstrated previously in a study of allozyme evolution in these populations (15–22). Our results demonstrated that polymorphisms in ALP genes in WEW correlated with the ecogeographic distribution of the genotypes. Observations were consistent with previous results on HMW-GS, LMW-GS, gliadins, and α-/β-amylase inhibitors (32–37). Some geographically close populations were very different in their TdALP structures at the considered loci (). An example can be found at Tabigha, now designated the evolution slope where two divergent Sos—the calcareous terra rossa soil and the volcanic basaltic soil (15, 38)—influenced the TdALP composition of WEW populations occupying these different soils. Alleles bx-7AS-g, ay-7AS-b, and ax-7AS-c occurred only in plants from terra rossa soil, whereas bx-7AS-n, ay-7AS-c, and ax-7AS-b occurred only in basaltic soil. The absence of a significant relationship between geographic separation and genetic distance attests to a sharp local ecological differentiation rather than a gradual change in allele frequencies across the range of WEW in Israel. Genetic diversity did not follow the simple isolation by the distance model of Wright et al. (39). Quite often, a greater genetic difference occurred between physically close populations than between distant populations. This was clearly demonstrated by the proximal populations located at Tabigha (two Sos) (39) and Yehudiyya (sun vs. shade) (17, 20–22). Thus, the genetic structure of WEW populations in Israel is mosaic. This patchy genetic distribution appears to reflect the underlying ecological, climatic, edaphic, and biotic heterogeneities on both micro- and macroscales (16, 18, 40). Microenvironmental variation coupled with a limited migration of T. dicoccoides may explain the dramatic genetic divergence of the two populations at the Tabigha site (15, 38). Specific SNP positions detected in TdALP genes were found to be highly effective in distinguishing genotypes and populations of WEW originating from diverse ecogeographic sites. These results suggest that genetic variations at these SNP positions in the TdALP were, at least, partly ecologically determined.

Natural Selection of TdALP-bx-7AS Genes in WEW.

Significant diversities at the TdALP-bx-7AS gene locus were detected both between and within WEW populations. The bx-7AS genes were naturally selected across populations supported by a different ratio of dN:dS (ω) (26). Environmental factors significantly correlate with the functional and silent alleles for the bx-7AS locus (). A sharp genetic divergence over short geographic distances compared with a small genetic divergence between long geographic distances also suggested that the SNPs were subjected to natural selection, and ecological factors played an important evolutionary role in gene polymorphism formation (15–22). Natural selection of orthologous genes can be assessed by comparing the ratio of ω in protein coding sequences (41). Ecological stresses have often been proposed as inducing active and rapid evolutionary changes. Compared with positive natural selection, purifying selection acts against mutations that have deleterious effects on protein structure. The ω value of Haplo 24 equals 1, displaying neutral selection (Fig. 3). The results indicated that all of the functional alleles are smaller than 1, suggesting that natural selection may have eliminated most of the deleterious effects caused by purifying selection (Fig. 3). On the other hand, all of the silent alleles and one functional allele Haplo 11 are greater than 1, suggesting that positive selection gives rise to dominant alleles in several WEW populations (Fig. 3). Altitude plays a significant role in population TdALP divergence as evidenced by the populations of Mt. Hermon, Rosh Pinna, Gamla, Bat Shelomo, and Tabigha, located at altitudes of 1300, 700, 200, 75, and 0 m, respectively (). The results showed that the populations located below 700 m (Tagbiha, 0 m) tend to have a higher level of genetic diversity with the He and I being 0.67 and 1.26, respectively (Fig. 2). In contrast, populations collected above 900 m, such as Mt. Hermon, were not polymorphic at all (Fig. 2). Along with altitude, several other environmental factors differed among these populations, such as abiotic climatic conditions, water availability, So, and biotic factors including parasites, pathogens, and competitors () (16, 42).

Genetic Distance and Evolution of TdALP in WEW.

The relationship between TdALP genetic distance and geographical distance indicated that the estimates of genetic distance (p distance) were geographically independent. Sharp genetic divergences (long p distance) over very short geographic distances against small genetic divergence (short p distance) between long physical distances were observed. For example, the genetic distance between populations of Tabigha (terra rossa and basalt) and Gamla located only about 9.5 km apart from p distance = 0.212, was 16 times higher than that between Yehudiyyasun and J’aba (130 km, ). Environmental stress can greatly influence plant susceptibility to herbivores and pathogens, and drought stress can promote outbreaks of fungal diseases and plant-eating insects (43). Different herbivore-related and pathogen-related selection pressures at these ecological locations may influence polymorphism of insect-resistant and pathogen-resistant loci in WEW (36). It can be concluded that the variation in ALP genetic diversity between populations is due to selective forces. The genetic structure of WEW populations in Israel is a mosaic (16, 17, 19, 33, 34, 44, 45). Thus, higher levels of polymorphisms and genetic variations of TdALP within and between populations can be explained as adaptive complexes generated by natural selection and coevolution with biotic or abiotic pressure.

Conclusions and Prospects.

Our molecular characterization of the TdALP gene family in WEW allows several conclusions to be made about the origin of ALP genes. Future challenges of crop improvement can be overcome by effectively utilizing the immense resources of genetic diversity unraveled by the evolution and allele analysis in natural populations of the wheat progenitors. The drivers of ALP allelic variations in WEW populations appear to be intimately linked to the environment in which the populations originated. These results suggest: (i) during the evolutionary history of WEW, diversifying natural selection through climatic (e.g., annual rain fall and temperature) and edaphic factors (So) was a major agent of genetic structure and differentiation at TdALP loci; and (ii) WEW populations harbor large amounts of genetic diversity exploitable for wheat improvement. Furthermore, at the transcriptional level, we found that most members of this multifunction large gene family are transcriptionally active at multiple stages of bread wheat development as well as under conditions of pathogen infection (powdery mildew). The allelic diversity associated with the germplasm-originating environmental conditions may provide a solution to fight the negative impact of the global warming complexities on modern wheat production. These genetic resources provide potential values for improving wheat cultivars under uncertain environmental conditions in the future.

Materials and Methods

The geographical sources of 21 WEW (T. dicoccoides) populations analyzed are shown (). WEW seeds were obtained from the GenBank of the Institute of Evolution, University of Haifa, Mount Carmel, Haifa, Israel. Mace × Spitfire DH population wheat lines were used for ALP gene transcriptional analysis. Details of experimental procedures, such as ALP genes' cloning, phylogenetic analysis, quantitative reverse transcription polymerase chain reaction, population genetic analysis, and statistical analysis, are described in .