Genetic architecture of rust resistance in a wheat (Triticum turgidum) diversity panel

INTRODUCTION

Wheat rust diseases are widespread and affect all wheat growing areas around the globe. Breeding strategies focus on incorporating genetic disease resistance. However, pathogens can quickly evolve and overcome the resistance genes deployed in commercial cultivars, creating a constant need for identifying new sources of resistance.

METHODS

We have assembled a diverse tetraploid wheat panel comprised of 447 accessions of three Triticum turgidum subspecies and performed a genome-wide association study (GWAS) for resistance to wheat stem, stripe, and leaf rusts. The panel was genotyped with the 90K Wheat iSelect single nucleotide polymorphism (SNP) array and subsequent filtering resulted in a set of 6,410 non-redundant SNP markers with known physical positions.

RESULTS

Population structure and phylogenetic analyses revealed that the diversity panel could be divided into three subpopulations based on phylogenetic/geographic relatedness. Marker-trait associations (MTAs) were detected for two stem rust, two stripe rust and one leaf rust resistance loci. Of them, three MTAs coincide with the known rust resistance genes Sr13, Yr15 and Yr67, while the other two may harbor undescribed resistance genes.

DISCUSSION

The tetraploid wheat diversity panel, developed and characterized herein, captures wide geographic origins, genetic diversity, and evolutionary history since domestication making it a useful community resource for mapping of other agronomically important traits and for conducting evolutionary studies.

Discovery of stripe rust resistance with incomplete dominance in wild emmer wheat using bulked segregant analysis sequencing

Durable crop disease resistance is an essential component of global food security. Continuous pathogen evolution leads to a breakdown of resistance and there is a pressing need to characterize new resistance genes for use in plant breeding. Here we identified an accession of wild emmer wheat (Triticum turgidum ssp. dicoccoides), PI 487260, that is highly resistant to multiple stripe rust isolates. Genetic analysis revealed resistance was conferred by a single, incompletely dominant gene designated as Yr84. Through bulked segregant analysis sequencing (BSA-Seq) we identified a 52.7 Mb resistance-associated interval on chromosome 1BS. Detected variants were used to design genetic markers for recombinant screening, further refining the interval of Yr84 to a 2.3-3.3 Mb in tetraploid wheat genomes. This interval contains 34 candidate genes encoding for protein domains involved in disease resistance responses. Furthermore, KASP markers closely-linked to Yr84 were developed to facilitate marker-assisted selection for rust resistance breeding.

TdPm60 identified in wild emmer wheat is an ortholog of Pm60 and constitutes a strong candidate for PmG16 powdery mildew resistance

We identified TdPm60 alleles from wild emmer wheat (WEW), an ortholog of Pm60 from T. urartu, which constitutes a strong candidate for PmG16 mildew resistance. Deployment of PmG16 in Israeli modern bread wheat cultivar Ruta improved the resistance to several local Bgt isolates. Wild emmer wheat (WEW), the tetraploid progenitor of durum and bread wheat, is a valuable genetic resource for resistance to powdery mildew fungal disease caused by Blumeria graminis f. sp. tritici (Bgt). PmG16 gene, derived from WEW, confers high resistance to most tested Bgt isolates. We mapped PmG16 to a 1.4-cM interval between the flanking markers uhw386 and uhw390 on Chromosome 7AL. Based on gene annotation of WEW reference genome Zavitan_V1, 34 predicted genes were identified within the ~ 3.48-Mb target region. Six genes were annotated as associated with disease resistance, of which TRIDC7AG077150.1 was found to be highly similar to Pm60, previously cloned from Triticum urartu, and resides in the same syntenic region. The functional molecular marker (FMM) for Pm60 (M-Pm60-S1) co-segregated with PmG16, suggesting the Pm60 ortholog from WEW (designated here as TdPm60) as a strong candidate for PmG16. Sequence alignment identified only eight SNPs that differentiate between TdPm60 and TuPm60. Furthermore, TdPm60 was found to be present also in the WEW donor lines of the powdery mildew resistance genes MlIW172 and MlIW72, mapped to the same region of Chromosome 7AL as PmG16, suggesting that TdPm60 constitutes a candidate also for these genes. Furthermore, screening of additional 230 WEW accessions with Pm60 specific markers revealed 58 resistant accessions from the Southern Levant that harbored TdPm60, while none of the susceptible accessions showed the presence of this gene. Deployment of PmG16 in Israeli modern bread wheat cultivar Ruta conferred resistance against several local Bgt isolates.

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.

GenoTypeMapper: graphical genotyping on genetic and sequence-based maps

BACKGROUND

The rising availability of assemblies of large genomes (e.g. bread and durum wheat, barley) and their annotations deliver the basis to graphically present genome organization of parents and progenies on a physical scale. Genetic maps are a very important tool for breeders but often represent distorted models of the actual chromosomes, e.g., in centromeric and telomeric regions. This biased picture might lead to imprecise assumptions and estimations about the size and complexity of genetic regions and the selection of suitable molecular markers for the incorporation of traits in breeding populations or near-isogenic lines (NILs). Some software packages allow the graphical illustration of genotypic data, but to the best of our knowledge, suitable software packages that allow the comparison of genotypic data on the physical and genetic scale are currently unavailable.

RESULTS

We developed a simple Java-based-software called GenoTypeMapper (GTM) for comparing genotypic data on genetic and physical maps and tested it for effectiveness on data of two NILs that carry QTL-regions for drought stress tolerance from wild emmer on chromosome 2BS and 7AS. Both NILs were more tolerant to drought stress than their recurrent parents but exhibited additional undesirable traits such as delayed heading time.

CONCLUSIONS

In this article, we illustrate that the software easily allows users to display and identify additional chromosomal introgressions in both NILs originating from the wild emmer parent. The ability to detect and diminish linkage drag can be of particular interest for pre-breeding purposes and the developed software is a well-suited tool in this respect. The software is based on a simple allele-matching algorithm between the offspring and parents of a crossing scheme. Despite this simple approach, GTM seems to be the only software that allows us to analyse, illustrate and compare genotypic data of offspring of different crossing schemes with up to four parents in two different maps. So far, up to 500 individuals with a maximum number of 50,000 markers can be examined with the software. The main limitation that hampers the performance of the software is the number of markers that are examined in parallel. Since each individual must be analysed separately, a maximum of ten individuals can currently be displayed in a single run. On a computer with an Intel five processor of the 8th generation, GTM can reliably either analyse a single individual with up to 12,000 markers or ten individuals with up to 3,600 markers in less than five seconds. Future work aims to improve the performance of the software so that more complex crossing schemes with more parents and more markers can be analysed.

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.

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.

Grain protein content and thousand kernel weight QTLs identified in a durum × wild emmer wheat mapping population tested in five environments

Genetic dissection of GPC and TKW in tetraploid durum × WEW RIL population, based on high-density SNP genetic map, revealed 12 GPC QTLs and 11 TKW QTLs, with favorable alleles for 11 and 5 QTLs, respectively, derived from WEW. Wild emmer wheat (Triticum turgidum ssp. dicoccoides, WEW) was shown to exhibit high grain protein content (GPC) and therefore possess a great potential for improvement of cultivated wheat nutritional value. Genetic dissection of thousand kernel weight (TKW) and grain protein content (GPC) was performed using a high-density genetic map constructed based on a recombinant inbred line (RIL) population derived from a cross between T. durum var. Svevo and WEW acc. Y12-3. Genotyping of 208 F₆ RILs with a 15 K wheat single nucleotide polymorphism (SNP) array yielded 4166 polymorphic SNP markers, of which 1510 were designated as skeleton markers. A total map length of 2169 cM was obtained with an average distance of 1.5 cM between SNPs. A total of 12 GPC QTLs and 11 TKW QTLs were found under five different environments. No significant correlations were found between GPC and TKW across all environments. Four major GPC QTLs with favorable alleles from WEW were found on chromosomes 4BS, 5AS, 6BS and 7BL. The 6BS GPC QTL coincided with the physical position of the NAC transcription factor TtNAM-B1, underlying the cloned QTL, Gpc-B1. Comparisons of the physical intervals of the GPC QTLs described here with the results previously reported in other durum × WEW RIL population led to the discovery of seven novel GPC QTLs. Therefore, our research emphasizes the importance of GPC QTL dissection in diverse WEW accessions as a source of novel alleles for improvement of GPC in cultivated wheat.

Grain protein content and thousand kernel weight QTLs identified in a durum × wild emmer wheat mapping population tested in five environments

Genetic dissection of GPC and TKW in tetraploid durum × WEW RIL population, based on high-density SNP genetic map, revealed 12 GPC QTLs and 11 TKW QTLs, with favorable alleles for 11 and 5 QTLs, respectively, derived from WEW. Wild emmer wheat (Triticum turgidum ssp. dicoccoides, WEW) was shown to exhibit high grain protein content (GPC) and therefore possess a great potential for improvement of cultivated wheat nutritional value. Genetic dissection of thousand kernel weight (TKW) and grain protein content (GPC) was performed using a high-density genetic map constructed based on a recombinant inbred line (RIL) population derived from a cross between T. durum var. Svevo and WEW acc. Y12-3. Genotyping of 208 F₆ RILs with a 15 K wheat single nucleotide polymorphism (SNP) array yielded 4166 polymorphic SNP markers, of which 1510 were designated as skeleton markers. A total map length of 2169 cM was obtained with an average distance of 1.5 cM between SNPs. A total of 12 GPC QTLs and 11 TKW QTLs were found under five different environments. No significant correlations were found between GPC and TKW across all environments. Four major GPC QTLs with favorable alleles from WEW were found on chromosomes 4BS, 5AS, 6BS and 7BL. The 6BS GPC QTL coincided with the physical position of the NAC transcription factor TtNAM-B1, underlying the cloned QTL, Gpc-B1. Comparisons of the physical intervals of the GPC QTLs described here with the results previously reported in other durum × WEW RIL population led to the discovery of seven novel GPC QTLs. Therefore, our research emphasizes the importance of GPC QTL dissection in diverse WEW accessions as a source of novel alleles for improvement of GPC in cultivated wheat.

Cloning of the wheat Yr15 resistance gene sheds light on the plant tandem kinase-pseudokinase family

Yellow rust, caused by Puccinia striiformis f. sp. tritici (Pst), is a devastating fungal disease threatening much of global wheat production. Race-specific resistance (R)-genes are used to control rust diseases, but the rapid emergence of virulent Pst races has prompted the search for a more durable resistance. Here, we report the cloning of Yr15, a broad-spectrum R-gene derived from wild emmer wheat, which encodes a putative kinase-pseudokinase protein, designated as wheat tandem kinase 1, comprising a unique R-gene structure in wheat. The existence of a similar gene architecture in 92 putative proteins across the plant kingdom, including the barley RPG1 and a candidate for Ug8, suggests that they are members of a distinct family of plant proteins, termed here tandem kinase-pseudokinases (TKPs). The presence of kinase-pseudokinase structure in both plant TKPs and the animal Janus kinases sheds light on the molecular evolution of immune responses across these two kingdoms.

Algal Bio-Indication in Assessment of Hydrological Impact on Ecosystem in Wetlands of “Slavyansky Resort”

Algal bio-indication is commonly used in water quality assessment but can also help in assessing the impact of hydrology on freshwater wetland ecosystems.We identified 350 species and infraspecific taxa of algae from nine taxonomic divisions (Cyanoprokaryota, Chrysophyta, Euglenophyta,Dinophyta,Xanthophyta,Cryptophyta,Bacillariophyta,Chlorophyta,Charophyta) in 121 phytoplankton samples collected between 2007-2013 from seven lakes in the wetlands of the Regional Landscape Park “Slavyansky Resort”, Ukraine. The algal species richness and phytoplankton biomass decreased as water salinity increased. In turn the water salinity was influenced by the inflow of groundwater, karst fracture and by the alluvial water tributaries of a paleoriver that affects the formation processes of lake-spring sulphide mud from the resort, which is often used for therapeutic purposes.