Transcriptomic insights into the molecular mechanism for response of wild emmer wheat to stripe rust fungus

Introduction

Continuous identification and application of novel resistance genes against stripe rust are of great importance for wheat breeding. Wild emmer wheat, Triticum dicoccoides, has adapted to a broad range of environments and is a valuable genetic resource that harbors important beneficial traits, including resistance to stripe rust caused by Puccinia striiformis f. sp. tritici (Pst). However, there has been a lack of systematic exploration of genes against Pst races in wild emmer wheat.

Methods

Genome-wide transcriptome profiles were conducted on two wild emmer wheat genotypes with different levels of resistance to (Pst (DR3 exhibiting moderate (Pst resistance, and D7 displaying high (Pst resistance). qRT-PCR was performed to verify findings by RNA-seq.

Results

A higher number of DEGs were identified in the moderately (Pst-resistant genotype, while the highly (Pst-resistant genotype exhibited a greater enrichment of pathways. Nonetheless, there were consistent patterns in the enrichment of pathways between the two genotypes at the same time of inoculation. At 24 hpi, a majority of pathways such as the biosynthesis of secondary metabolites, phenylpropanoid biosynthesis, phenylalanine metabolism, and alpha-Linolenic acid metabolism exhibited significant enrichment in both genotypes. At 72 hpi, the biosynthesis of secondary metabolites and circadian rhythm-plant pathways were notably and consistently enriched in both genotypes. The majority of (WRKY, MADs , and AP2-ERF families were found to be involved in the initial stage of response to Pst invasion (24 hpi), while the MYB, NAC, TCP, and b-ZIP families played a role in defense during the later stage of Pst infection (72 hpi).

Discussion

In this present study, we identified numerous crucial genes, transcription factors, and pathways associated with the response and regulation of wild emmer wheat to Pst infection. Our findings offer valuable information for understanding the function of crucial Pst-responsive genes, and will deepen the understanding of the complex resistance mechanisms against Pst in wheat.

Wild emmer wheat, the progenitor of modern bread wheat, exhibits great diversity in the VERNALIZATION1 gene

Wild emmer wheat is an excellent reservoir of genetic variability that can be utilized to improve cultivated wheat to address the challenges of the expanding world population and climate change. Bearing this in mind, we have collected a panel of 263 wild emmer wheat (WEW) genotypes across the Fertile Crescent. The genotypes were grown in different locations and phenotyped for heading date. Genome-wide association mapping (GWAS) was carried out, and 16 SNPs were associated with the heading date. As the flowering time is controlled by photoperiod and vernalization, we sequenced the VRN1 gene, the most important of the vernalization response genes, to discover new alleles. Unlike most earlier attempts, which characterized known VRN1 alleles according to a partial promoter or intron sequences, we obtained full-length sequences of VRN-A1 and VRN-B1 genes in a panel of 95 wild emmer wheat from the Fertile Crescent and uncovered a significant sequence variation. Phylogenetic analysis of VRN-A1 and VRN-B1 haplotypes revealed their evolutionary relationships and geographic distribution in the Fertile Crescent region. The newly described alleles represent an attractive resource for durum and bread wheat improvement programs.

dicoccoides) introduced into elite wheat varieties

Drought events or the combination of drought and heat conditions are expected to become more frequent due to global warming, and wheat yields may fall below their long-term average. One way to increase climate-resilience of modern high-yielding varieties is by their genetic improvement with beneficial alleles from crop wild relatives. In the present study, the effect of two beneficial QTLs introgressed from wild emmer wheat and incorporated in the three wheat varieties BarNir, Zahir and Uzan was studied under well-watered conditions and under drought stress using non-destructive High-throughput Phenotyping (HTP) throughout the life cycle in a single pot-experiment. Plants were daily imaged with RGB top and side view cameras and watered automatically. Further, at two time points, the quantum yield of photosystem II was measured with a top view FluorCam. The QTL carrying near isogenic lines (NILs) were compared with their corresponding parents by t-test for all non-invasively obtained traits and for the manually determined agronomic and yield parameters. Data quality of phenotypic traits (repeatability) in the controlled HTP experiment was above 85% throughout the life cycle and at maturity. Drought stress had a strong effect on growth in all wheat genotypes causing biomass reduction from 2% up to 70% at early and late points in the drought period, respectively. At maturity, the drought caused 47-55% decreases in yield-related traits grain weight, straw weight and total biomass and reduced TKW by 10%, while water use efficiency (WUE) increased under drought by 29%. The yield-enhancing effect of the introgressed QTLs under drought conditions that were previously demonstrated under field/screenhouse conditions in Israel, could be mostly confirmed in a greenhouse pot experiment using HTP. Daily precision phenotyping enabled to decipher the mode of action of the QTLs in the different genetic backgrounds throughout the entire wheat life cycle. Daily phenotyping allowed a precise determination of the timing and size of the QTLs effect (s) and further yielded information about which image-derived traits are informative at which developmental stage of wheat during the entire life cycle. Maximum height and estimated biovolume were reached about a week after heading, so experiments that only aim at exploring these traits would not need a longer observation period. To obtain information on different onset and progress of senescence, the CVa curves represented best the ongoing senescence of plants. The QTL on 7A in the BarNir background was found to improve yield under drought by increased biomass growth, a higher photosynthetic performance, a higher WUE and a "stay green effect."

Genome-Wide Survey and Functional Verification of the NAC Transcription Factor Family in Wild Emmer Wheat

The NAC transcription factor (TF) family is one of the largest TF families in plants, which has been widely reported in rice, maize and common wheat. However, the significance of the NAC TF family in wild emmer wheat (Triticum turgidum ssp. dicoccoides) is not yet well understood. In this study, a genome-wide investigation of NAC genes was conducted in the wild emmer genome and 249 NAC family members (TdNACs) were identified. The results showed that all of these genes contained NAM/NAC-conserved domains and most of them were predicted to be located on the nucleus. Phylogenetic analysis showed that these 249 TdNACs can be classified into seven clades, which are likely to be involved in the regulation of grain protein content, starch synthesis and response to biotic and abiotic stresses. Expression pattern analysis revealed that TdNACs were highly expressed in different wheat tissues such as grain, root, leaves and shoots. We found that TdNAC8470 was phylogenetically close to NAC genes that regulate either grain protein or starch accumulation. Overexpression of TdNAC8470 in rice showed increased grain starch concentration but decreased grain Fe, Zn and Mn contents compared with wild-type plants. Protein interaction analysis indicated that TdNAC8470 might interact with granule-bound starch synthase 1 (TdGBSS1) to regulate grain starch accumulation. Our work provides a comprehensive understanding of the NAC TFs family in wild emmer wheat and establishes the way for future functional analysis and genetic improvement of increasing grain starch content in wheat.

Exploiting Wild Emmer Wheat Diversity to Improve Wheat A and B Genomes in Breeding for Heat Stress Adaptation

Wheat is highly sensitive to temperature beyond the optimum. To improve wheat adaptation to heat stress, the best option is to exploit the diversity of wild wheat progenitors. This study aimed to identify germplasm and quantitative trait loci associated with heat stress tolerance from wild emmer wheat diversity. We evaluated a diverse set of multiple derivative lines harboring chromosome segments from nine wild emmer wheat parents under four environments: two optimum environments at Tottori, Japan and Dongola, Sudan, one moderate heat stress environment, and one severe heat stress environment at Wad Medani, Sudan. Genome-wide association analysis was conducted with 13,312 SNP markers. Strong marker-trait associations (MTAs) were identified for chlorophyll content at maturity on chromosomes 1A and 5B: these MTAs explained 28.8 and 26.8% of the variation, respectively. A region on chromosome 3A (473.7-638.4 Mbp) contained MTAs controlling grain yield, under optimum and severe heat stress. Under severe heat stress, regions on chromosomes 3A (590.4-713.3 Mbp) controlled grain yield, biomass, days to maturity and thousand kernel weight, and on 3B (744.0-795.2 Mbp) grain yield and biomass. Heat tolerance efficiency (HTE) was controlled by three MTAs, one each on chromosomes 2A, 2B, and 5A under moderate heat stress and one MTA on chromosome 3A under severe heat stress. Some of the MTAs found here were previously reported, but the new ones originated from the wild emmer wheat genomes. The favorable alleles identified from wild emmer wheat were absent or rare in the elite durum wheat germplasm being bred for heat stress tolerance. This study provides potential genetic materials, alleles, MTAs, and quantitative trait loci for enhancing wheat adaptation to heat stress. The derivative lines studied here could be investigated to enhance other stress tolerance such as drought and salinity.

Genetic variation and genetic control of intraspikelet differences in grain weight and seed dormancy in wild and domesticated emmer wheats

Seed dormancy, a vital strategy for wild plant species to adapt to an unpredictable environment in their natural habitats, was eliminated from cereals during the domestication process. Intraspikelet differences in grain size and seed dormancy have been observed in wild emmer wheat. To elucidate the genetic variation of these intraspikelet differences and to determine their genetic control, grain weight ratio (first florets/second florets) (GWR), germination rate, and germination index (GI) were analyzed in 67 wild and 82 domesticated emmer wheat accessions, as well as F₁ hybrids, F₂ populations, and F₃-F₆ populations derived from reciprocal crosses between wild and domesticated lines. Only the grains on the first florets of two-grained spikelets in wild accessions had varying degrees of dormancy with GI ranging from 0 to 1, which positively correlated with their GWR. This implies that wild emmer populations comprised genotypes with varying degrees of dormancy, including nondormant genotypes. According to segregations observed in F₂ populations, the intraspikelet grain weight difference was controlled by two independently inherited loci. Furthermore, low-GWR populations with low or high GI values could be selected in F₅ and F₆ generations, implying that the major loci associated with dormancy might be independent of intraspikelet grain weight difference.

Modifying root-to-shoot ratio improves root water influxes in wheat under drought stress

Drought intensity as experienced by plants depends upon soil moisture status and atmospheric variables such as temperature, radiation, and air vapour pressure deficit. Although the role of shoot architecture with these edaphic and atmospheric factors is well characterized, the extent to which shoot and root dynamic interactions as a continuum are controlled by genotypic variation is less well known. Here, we targeted these interactions using a wild emmer wheat introgression line (IL20) with a distinct drought-induced shift in the shoot-to-root ratio and its drought-sensitive recurrent parent Svevo. Using a gravimetric platform, we show that IL20 maintained higher root water influx and gas exchange under drought stress, which supported a greater growth. Interestingly, the advantage of IL20 in root water influx and transpiration was expressed earlier during the daily diurnal cycle under lower vapour pressure deficit and therefore supported higher transpiration efficiency. Application of a structural equation model indicates that under drought, vapour pressure deficit and radiation are antagonistic to transpiration rate, whereas the root water influx operates as a feedback for the higher atmospheric responsiveness of leaves. Collectively, our results suggest that a drought-induced shift in root-to-shoot ratio can improve plant water uptake potential in a short preferable time window during early morning when vapour pressure deficit is low and the light intensity is not a limiting factor for assimilation.

Deciphering the genetic basis of wheat seminal root anatomy uncovers ancestral axial conductance alleles

Root axial conductance, which describes the ability of water to move through the xylem, contributes to the rate of water uptake from the soil throughout the whole plant lifecycle. Under the rainfed wheat agro-system, grain-filling is typically occurring during declining water availability (i.e., terminal drought). Therefore, preserving soil water moisture during grain filling could serve as a key adaptive trait. We hypothesized that lower wheat root axial conductance can promote higher yields under terminal drought. A segregating population derived from a cross between durum wheat and its direct progenitor wild emmer wheat was used to underpin the genetic basis of seminal root architectural and functional traits. We detected 75 QTL associated with seminal roots morphological, anatomical and physiological traits, with several hotspots harbouring co-localized QTL. We further validated the axial conductance and central metaxylem QTL using wild introgression lines. Field-based characterization of genotypes with contrasting axial conductance suggested the contribution of low axial conductance as a mechanism for water conservation during grain filling and consequent increase in grain size and yield. Our findings underscore the potential of harnessing wild alleles to reshape the wheat root system architecture and associated hydraulic properties for greater adaptability under changing climate.

Genome-Wide Mapping of Loci for Adult-Plant Resistance to Stripe Rust in Durum Wheat Svevo Using the 90K SNP Array

Stripe rust is a foliar disease in wheat caused by Puccinia striiformis f. tritici. The best way to protect wheat from this disease is by growing resistant cultivars. Tetraploid wheat can serve as a good source of valuable genetic diversity for various traits. Here, we report the mapping of nine stripe rust resistance quantitative trait loci (QTL) effective against P. striiformis f. tritici in China and Israel. We used recombinant inbred lines (RILs) developed from a cross between the durum wheat cultivar Svevo and Triticum dicoccoides accession Zavitan. By genotyping the RIL population of 137 lines using the wheat 90K single-nucleotide polymorphism array, we mapped an adult-plant resistance locus QYrsv.swust-1BL.1, the most effective QTL, within a 0.75-centimorgan region in T. turgidum subsp. durum 'Svevo' on chromosome arm 1BL, corresponding to the region of 670.7 to 671.5 Mb on the Chinese Spring chromosome arm 1BL. Of the other eight minor-effect stripe rust QTL, seven were from Svevo and mapped on chromosomes 1A, 1B, 2B, 3A, 4A, and 5A, and one was from Zavitan and mapped on chromosome 2A. Several QTL with epistatic effects were identified as well. The markers linked to the resistance QTL can be useful in marker-assisted selection for incorporation of these resistance QTL into both durum and common wheat cultivars.

The Independent Domestication of Timopheev's Wheat: Insights from Haplotype Analysis of the Brittle rachis 1 (BTR1-A) Gene

Triticum turgidum and T. timopheevii are two tetraploid wheat species sharing T. urartu as a common ancestor, and domesticated accessions from both of these allopolyploids exhibit nonbrittle rachis (i.e., nonshattering spikes). We previously described the loss-of-function mutations in the Brittle Rachis 1 genes BTR1-A and BTR1-B in the A and B subgenomes, respectively, that are responsible for this most visible domestication trait in T. turgidum. Resequencing of a large panel of wild and domesticated T. turgidum accessions subsequently led to the identification of the two progenitor haplotypes of the btr1-A and btr1-B domesticated alleles. Here, we extended the haplotype analysis to other T. turgidum subspecies and to the BTR1 homologues in the related T. timopheevii species. Our results showed that all the domesticated wheat subspecies within T. turgidum share common BTR1-A and BTR1-B haplotypes, confirming their common origin. In T. timopheevii, however, we identified a novel loss-of-function btr1-A allele underlying a partially brittle spike phenotype. This novel recessive allele appeared fixed within the pool of domesticated Timopheev’s wheat but was also carried by one wild timopheevii accession exhibiting partial brittleness. The promoter region for BTR1-B could not be amplified in any T. timopheevii accessions with any T. turgidum primer combination, exemplifying the gene-level distance between the two species. Altogether, our results support the concept of independent domestication processes for the two polyploid, wheat-related species.

Genomic Architecture of Phenotypic Plasticity in Response to Water Stress in Tetraploid Wheat

Phenotypic plasticity is one of the main mechanisms of adaptation to abiotic stresses via changes in critical developmental stages. Altering flowering phenology is a key evolutionary strategy of plant adaptation to abiotic stresses, to achieve the maximum possible reproduction. The current study is the first to apply the linear regression residuals as drought plasticity scores while considering the variation in flowering phenology and traits under non-stress conditions. We characterized the genomic architecture of 17 complex traits and their drought plasticity scores for quantitative trait loci (QTL) mapping, using a mapping population derived from a cross between durum wheat (Triticum turgidum ssp. durum) and wild emmer wheat (T. turgidum ssp. dicoccoides). We identified 79 QTLs affected observed traits and their plasticity scores, of which 33 reflected plasticity in response to water stress and exhibited epistatic interactions and/or pleiotropy between the observed and plasticity traits. Vrn-B3 (TaTF1) residing within an interval of a major drought-escape QTL was proposed as a candidate gene. The favorable alleles for most of the plasticity QTLs were contributed by wild emmer wheat, demonstrating its high potential for wheat improvement. Our study presents a new approach for the quantification of plant adaptation to various stresses and provides new insights into the genetic basis of wheat complex traits under water-deficit stress.

Current Progress in Understanding and Recovering the Wheat Genes Lost in Evolution and Domestication

The modern cultivated wheat has passed a long evolution involving origin of wild emmer (WEM), development of cultivated emmer, formation of spelt wheat and finally establishment of modern bread wheat and durum wheat. During this evolutionary process, rapid alterations and sporadic changes in wheat genome took place, due to hybridization, polyploidization, domestication, and mutation. This has resulted in some modifications and a high level of gene loss. As a result, the modern cultivated wheat does not contain all genes of their progenitors. These lost genes are novel for modern wheat improvement. Exploring wild progenitor for genetic variation of important traits is directly beneficial for wheat breeding. WEM wheat (Triticum dicoccoides) is a great genetic resource with huge diversity for traits. Few genes and quantitative trait loci (QTL) for agronomic, quantitative, biotic and abiotic stress-related traits have already been mapped from WEM. This resource can be utilized for modern wheat improvement by integrating identified genes or QTLs through breeding.

A Post-Haustorial Defense Mechanism is Mediated by the Powdery Mildew Resistance Gene, PmG3M, Derived from Wild Emmer Wheat

The destructive wheat powdery mildew disease is caused by the fungal pathogen Blumeria graminis f. sp. tritici (Bgt). PmG3M, derived from wild emmer wheat Triticum dicoccoides accession G305-3M, is a major gene providing a wide-spectrum resistance against Bgt. PmG3M was previously mapped to wheat chromosome 6B using an F₆ recombinant inbred line (RIL) mapping population generated by crossing G305-3M with the susceptible T. durum wheat cultivar Langdon (LDN). In the current study, we aimed to explore the defense mechanisms conferred by PmG3M against Bgt. Histopathology of fungal development was characterized in artificially inoculated leaves of G305-3M, LDN, and homozygous RILs using fluorescence and light microscopy. G305-3M exhibited H₂O₂ accumulation typical of a hypersensitive response, which resulted in programmed cell death (PCD) in Bgt-penetrated epidermal cells, while LDN showed well-developed colonies without PCD. In addition, we observed a post-haustorial resistance mechanism that arrested the development of fungal feeding structures and pathogen growth in both G305-3M and resistant RIL, while LDN and a susceptible RIL displayed fully developed digitated haustoria and massive accumulation of fungal biomass. In contrast, both G305-3M and LDN exhibited callose deposition in attempt to prevent fungal invasion, supporting this as a mechanism of a basal defense response not associated with PmG3M resistance mechanism per se. The presented results shed light on the resistance mechanisms conferred by PmG3M against wheat powdery mildew.

Three previously characterized resistances to yellow rust are encoded by a single locus Wtk1

The wild emmer wheat (Triticum turgidum ssp. dicoccoides; WEW) yellow (stripe) rust resistance genes Yr15, YrG303, and YrH52 were discovered in natural populations from different geographic locations. They all localize to chromosome 1B but were thought to be non-allelic based on differences in resistance response. We recently cloned Yr15 as a Wheat Tandem Kinase 1 (WTK1) and show here that these three resistance loci co-segregate in fine-mapping populations and share an identical full-length genomic sequence of functional Wtk1. Independent ethyl methanesulfonate (EMS)-mutagenized susceptible yrG303 and yrH52 lines carried single nucleotide mutations in Wtk1 that disrupted function. A comparison of the mutations for yr15, yrG303, and yrH52 mutants showed that while key conserved residues were intact, other conserved regions in critical kinase subdomains were frequently affected. Thus, we concluded that Yr15-, YrG303-, and YrH52-mediated resistances to yellow rust are encoded by a single locus, Wtk1. Introgression of Wtk1 into multiple genetic backgrounds resulted in variable phenotypic responses, confirming that Wtk1-mediated resistance is part of a complex immune response network. WEW natural populations subjected to natural selection and adaptation have potential to serve as a good source for evolutionary studies of different traits and multifaceted gene networks.

Sympatric speciation of wild emmer wheat driven by ecology and chromosomal rearrangements

In plants, the mechanism for ecological sympatric speciation (SS) is little known. Here, after ruling out the possibility of secondary contact, we show that wild emmer wheat, at the microclimatically divergent microsite of "Evolution Canyon" (EC), Mt. Carmel, Israel, underwent triple SS. Initially, it split following a bottleneck of an ancestral population, and further diversified to three isolated populations driven by disruptive ecological selection. Remarkably, two postzygotically isolated populations (SFS1 and SFS2) sympatrically branched within an area less than 30 m at the tropical hot and dry savannoid south-facing slope (SFS). A series of homozygous chromosomal rearrangements in the SFS1 population caused hybrid sterility with the SFS2 population. We demonstrate that these two populations developed divergent adaptive mechanisms against severe abiotic stresses on the tropical SFS. The SFS2 population evolved very early flowering, while the SFS1 population alternatively evolved a direct tolerance to irradiance by improved ROS scavenging activity that potentially accounts for its evolutionary fate with unstable chromosome status. Moreover, a third prezygotically isolated sympatric population adapted on the abutting temperate, humid, cool, and forested north-facing slope (NFS), separated by 250 m from the SFS wild emmer wheat populations. The NFS population evolved multiple resistant loci to fungal diseases, including powdery mildew and stripe rust. Our study illustrates how plants sympatrically adapt and speciate under disruptive ecological selection of abiotic and biotic stresses.

Variation in phosphorus and sulfur content shapes the genetic architecture and phenotypic associations within the wheat grain ionome

Dissection of the genetic basis of wheat ionome is crucial for understanding the physiological and biochemical processes underlying mineral accumulation in seeds, as well as for efficient crop breeding. Most of the elements essential for plants are metals stored in seeds as chelate complexes with phytic acid or sulfur-containing compounds. We assume that the involvement of phosphorus and sulfur in metal chelation is the reason for strong phenotypic correlations within ionome. Adjustment of element concentrations for the effect of variation in phosphorus and sulfur seed content resulted in drastic change of phenotypic correlations between the elements. The genetic architecture of wheat grain ionome was characterized by quantitative trait loci (QTL) analysis using a cross between durum and wild emmer wheat. QTL analysis of the adjusted traits and two-trait analysis of the initial traits paired with either P or S considerably improved QTL detection power and accuracy, resulting in the identification of 105 QTLs and 617 QTL effects for 11 elements. Candidate gene search revealed some potential functional associations between QTLs and corresponding genes within their intervals. Thus, we have shown that accounting for variation in P and S is crucial for understanding of the physiological and genetic regulation of mineral composition of wheat grain ionome and can be implemented for other plants.

Variation in phosphorus and sulfur content shapes the genetic architecture and phenotypic associations within the wheat grain ionome

Dissection of the genetic basis of wheat ionome is crucial for understanding the physiological and biochemical processes underlying mineral accumulation in seeds, as well as for efficient crop breeding. Most of the elements essential for plants are metals stored in seeds as chelate complexes with phytic acid or sulfur-containing compounds. We assume that the involvement of phosphorus and sulfur in metal chelation is the reason for strong phenotypic correlations within ionome. Adjustment of element concentrations for the effect of variation in phosphorus and sulfur seed content resulted in drastic change of phenotypic correlations between the elements. The genetic architecture of wheat grain ionome was characterized by quantitative trait loci (QTL) analysis using a cross between durum and wild emmer wheat. QTL analysis of the adjusted traits and two-trait analysis of the initial traits paired with either P or S considerably improved QTL detection power and accuracy, resulting in the identification of 105 QTLs and 617 QTL effects for 11 elements. Candidate gene search revealed some potential functional associations between QTLs and corresponding genes within their intervals. Thus, we have shown that accounting for variation in P and S is crucial for understanding of the physiological and genetic regulation of mineral composition of wheat grain ionome and can be implemented for other plants.

Elevated mutation and selection in wild emmer wheat in response to 28 years of global warming

The realized threats of global warming to biodiversity have catalyzed the search for a solution to protect and conserve extant plant genetic resources. Part of the solution, however, is dependent on the knowledge of how plant populations respond genetically to these threats, which is largely lacking. We conducted a unique genomic characterization of genetic responses in 10 wild emmer wheat populations in Israel that were sampled twice in 1980 and 2008. After the 28 y of global warming, these populations displayed elevated selection, reduced diversity and temporal divergence, and carried increased mutational burdens forward. However, some populations still showed the ability to acquire beneficial alleles for future adaptation. The patterns of genetic response to rainfall and temperature were complex.

Global warming has been documented to threaten wild plants with strong selection pressures, but how plant populations respond genetically to the threats remains poorly understood. We characterized the genetic responses of 10 wild emmer wheat (Triticum dicoccoides Koern.; WEW) populations in Israel, sampling them in 1980 and again in 2008, through an exome capture analysis. It was found that these WEW populations were under elevated selection, displayed reduced diversity and temporal divergence, and carried increased mutational burdens forward. However, some populations still showed the ability to acquire beneficial alleles via selection or de novo mutation for future adaptation. Grouping populations with mean annual rainfall and temperature revealed significant differences in most of the 14 genetic estimates in either sampling year or over the 28 y. The patterns of genetic response to rainfall and temperature varied and were complex. In general, temperature groups displayed more temporal differences in genetic response than rainfall groups. The highest temperature group had more deleterious single nucleotide polymorphisms (dSNPs), higher nucleotide diversity, fewer selective sweeps, lower differentiation, and lower mutational burden. The least rainfall group had more dSNPs, higher nucleotide diversity, lower differentiation and higher mutational burden. These characterized genetic responses are significant, allowing not only for better understanding of evolutionary changes in the threatened populations, but also for realistic modeling of plant population adaptability and vulnerability to global warming.

Genome-Wide Association Study Reveals Novel Genomic Regions Associated With High Grain Protein Content in Wheat Lines Derived From Wild Emmer Wheat

Grain protein content (GPC) and yield are of two important traits in wheat, but their negative correlation has hampered their simultaneous improvement in conventional breeding. Wild emmer wheat (Triticum turgidum ssp. dicoccoides) is an important genetic resource for wheat quality improvement. In this study, we report a genome-wide association study (GWAS) using 13116 DArT-seq markers to characterize GPC in 161 wheat lines derived from wild emmer. Using a general linear model, we identified 141 markers that were significantly associated with GPC, and grouped into 48 QTL regions. Using both general linear model and mixed linear model, we identified four significant markers that were grouped into two novel QTL regions on chromosomes 2BS (QGpc.cd1-2B.1) and 7BL (QGpc.cd1-7B.2). The two QTLs have no negative effects on thousand kernel weight (TKW) and should be useful for simultaneous improvement of GPC and TKW in wheat breeding. Searches of public databases revealed 61 putative candidate/flanking genes related to GPC. The putative proteins of interest were grouped in four main categories: enzymes, kinase proteins, metal transport-related proteins, and disease resistance proteins. The linked markers and associated candidate genes provide essential information for cloning genes related to high GPC and performing marker-assisted breeding in wheat.

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.

Genome Based Meta-QTL Analysis of Grain Weight in Tetraploid Wheat Identifies Rare Alleles of GRF4 Associated with Larger Grains

The domestication and subsequent genetic improvement of wheat led to the development of large-seeded cultivated wheat species relative to their smaller-seeded wild progenitors. While increased grain weight (GW) continues to be an important goal of many wheat breeding programs, few genes underlying this trait have been identified despite an abundance of studies reporting quantitative trait loci (QTL) for GW. Here we perform a QTL analysis for GW using a population of recombinant inbred lines (RILs) derived from the cross between wild emmer wheat accession ‘Zavitan’ and durum wheat variety ‘Svevo’. Identified QTLs in this population were anchored to the recent Zavitan reference genome, along with previously published QTLs for GW in tetraploid wheat. This genome-based, meta-QTL analysis enabled the identification of a locus on chromosome 6A whose introgression from wild wheat positively affects GW. The locus was validated using an introgression line carrying the 6A GW QTL region from Zavitan in a Svevo background, resulting in >8% increase in GW compared to Svevo. Using the reference sequence for the 6A QTL region, we identified a wheat ortholog to OsGRF4, a rice gene previously associated with GW. The coding sequence of this gene (TtGRF4-A) contains four single nucleotide polymorphisms (SNPs) between Zavitan and Svevo, one of which reveals the Zavitan allele to be rare in a core collection of wild emmer and completely absent from the domesticated emmer genepool. Similarly, another wild emmer accession (G18-16) was found to carry a rare allele of TtGRF4-A that also positively affects GW and is characterized by a unique SNP absent from the entire core collection. These results exemplify the rich genetic diversity of wild wheat, posit TtGRF4-A as a candidate gene underlying the 6A GW QTL, and suggest that the natural Zavitan and G18-16 alleles of TtGRF4-A have potential to increase wheat yields in breeding programs.

Chromosome-based survey sequencing reveals the genome organization of wild wheat progenitor Triticum dicoccoides

Wild emmer wheat (Triticum turgidum ssp. dicoccoides) is the progenitor of wheat. We performed chromosome-based survey sequencing of the 14 chromosomes, examining repetitive sequences, protein-coding genes, miRNA/target pairs and tRNA genes, as well as syntenic relationships with related grasses. We found considerable differences in the content and distribution of repetitive sequences between the A and B subgenomes. The gene contents of individual chromosomes varied widely, not necessarily correlating with chromosome size. We catalogued candidate agronomically important loci, along with new alleles and flanking sequences that can be used to design exome sequencing. Syntenic relationships and virtual gene orders revealed several small-scale evolutionary rearrangements, in addition to providing evidence for the 4AL-5AL-7BS translocation in wild emmer wheat. Chromosome-based sequence assemblies contained five novel miRNA families, among 59 families putatively encoded in the entire genome which provide insight into the domestication of wheat and an overview of the genome content and organization.

Activation of seminal root primordia during wheat domestication reveals underlying mechanisms of plant resilience

Seminal roots constitute the initial wheat root system and provide the main route for water absorption during early stages of development. Seminal root number (SRN) varies among species. However, the mechanisms through which SRN is controlled and in turn contribute to environmental adaptation are poorly understood. Here, we show that SRN increased upon wheat domestication from 3 to 5 due to the activation of 2 root primordia that are suppressed in wild wheat, a trait controlled by loci expressed in the germinating embryo. Suppression of root primordia did not limit water uptake, indicating that 3 seminal roots is adequate to maintain growth during seedling development. The persistence of roots at their primordial state promoted seedling recovery from water stress through reactivation of suppressed primordia upon rehydration. Our findings suggest that under well-watered conditions, SRN is not a limiting factor, and excessive number of roots may be costly and maladaptive. Following water stress, lack of substantial root system suppresses growth and rapid recovery of the root system is essential for seedling recovery. This study underscores SRN as key adaptive trait that was reshaped upon domestication. The maintenance of roots at their primordial state during seedling development may be regarded as seedling protective mechanism against water stress.

Reciprocal Hosts' Responses to Powdery Mildew Isolates Originating from Domesticated Wheats and Their Wild Progenitor

The biotroph wheat powdery mildew, Blumeria graminis (DC.) E.O. Speer, f. sp. tritici Em. Marchal (Bgt), has undergone long and dynamic co-evolution with its hosts. In the last 10,000 years, processes involved in plant evolution under domestication, altered host-population structure. Recently both virulence and genomic profiling separated Bgt into two groups based on their origin from domestic host and from wild emmer wheat. While most studies focused on the Bgt pathogen, there is significant knowledge gaps in the role of wheat host diversity in this specification. This study aimed to fill this gap by exploring qualitatively and also quantitatively the disease response of diverse host panel to powdery mildew [105 domesticated wheat genotypes (Triticum turgidum ssp. dicoccum, T. turgidum ssp. durum, and T. aestivum) and 241 accessions of its direct progenitor, wild emmer wheat (T. turgidum ssp. dicoccoides)]. A set of eight Bgt isolates, originally collected from domesticated and wild wheat was used for screening this wheat collection. The isolates from domesticated wheat elicited susceptible to moderate plant responses on domesticated wheat lines and high resistance on wild genotypes (51.7% of the tested lines were resistant). Isolates from wild emmer elicited reciprocal disease responses: high resistance of domesticated germplasm and high susceptibility of the wild material (their original host). Analysis of variance of the quantitative phenotypic responses showed a significant Isolates × Host species interaction [P(F) < 0.0001] and further supported these findings. Furthermore, analysis of the range of disease severity values showed that when the group of host genotypes was inoculated with Bgt isolate from the reciprocal host, coefficient of variation was significantly higher than when inoculated with its own isolates. This trend was attributed to the role of major resistance genes in the latter scenario (high proportion of complete resistance). By testing the association between disease severity and geographical distance from the source of inoculum, we have found higher susceptibility in wild emmer close to the source. Both qualitative and quantitative assays showed a reciprocal resistance pattern in the wheat host and are well aligned with the recent findings of significant differentiation into wild-emmer and domesticated-wheat populations in the pathogen.

Solar Radiation-Associated Adaptive SNP Genetic Differentiation in Wild Emmer Wheat, Triticum dicoccoides

Whole-genome scans with large number of genetic markers provide the opportunity to investigate local adaptation in natural populations and identify candidate genes under positive selection. In the present study, adaptation genetic differentiation associated with solar radiation was investigated using 695 polymorphic SNP markers in wild emmer wheat originated in a micro-site at Yehudiyya, Israel. The test involved two solar radiation niches: (1) sun, in-between trees; and (2) shade, under tree canopy, separated apart by a distance of 2-4 m. Analysis of molecular variance showed a small (0.53%) but significant portion of overall variation between the sun and shade micro-niches, indicating a non-ignorable genetic differentiation between sun and shade habitats. Fifty SNP markers showed a medium (0.05 ≤ FST ≤ 0.15) or high genetic differentiation (FST > 0.15). A total of 21 outlier loci under positive selection were identified by using four different FST -outlier testing algorithms. The markers and genome locations under positive selection are consistent with the known patterns of selection. These results suggested that genetic differentiation between sun and shade habitats is substantial, radiation-associated, and therefore ecologically determined. Hence, the results of this study reflected effects of natural selection through solar radiation on EST-related SNP genetic diversity, resulting presumably in different adaptive complexes at a micro-scale divergence. The present work highlights the evolutionary theory and application significance of solar radiation-driven natural selection in wheat improvement.