Aegilops umbellulata introgression carrying leaf rust and stripe rust resistance genes Lr76 and Yr70 located to 9.47-Mb region on 5DS telomeric end through a combination of chromosome sorting and sequencing.

KEY MESSAGE: Lr76 and Yr70 have been fine mapped using the sequence of flow-sorted recombinant 5D chromosome from wheat-Ae. umbellulata introgression line. The alien introgression has been delineated to 9.47-Mb region on short arm of wheat chromosome 5D. Leaf rust and stripe rust are among the most damaging diseases of wheat worldwide. Wheat cultivation based on limited number of rust resistance genes deployed over vast areas expedites the emergence of new pathotypes warranting a continuous deployment of new resistance genes. In this paper, fine mapping of Aegilops umbellulata-derived leaf rust and stripe rust resistance genes Lr76 and Yr70 is being reported. We flow sorted and paired-end sequenced 5U chromosome of Ae. umbellulata, recombinant chromosome 5D (5DIL) from wheat-Ae. umbellulata introgression line pau16057 and 5DRP of recurrent parent WL711. Chromosome 5U reads were mapped against the reference Chinese Spring chromosome 5D sequence, and alien-specific SNPs were identified. Chromosome 5DIL and 5DRP sequences were de novo assembled, and alien introgression-specific markers were designed by selecting 5U- and 5D-specific SNPs. Overall, 27 KASP markers were mapped in high-resolution population consisting of 1404 F5 RILs. The mapping population segregated for single gene each for leaf rust and stripe rust resistance. The physical order of the SNPs in pau16057 was defined by projecting the 27 SNPs against the IWGSC RefSeq v1.0 sequence. Based on this physical map, the size of Ae. umbellulata introgression was determined to be 9.47 Mb on the distal most end of the short arm of chromosome 5D. This non-recombining alien segment carries six NB-LRR encoding genes based on NLR annotation of assembled chromosome 5DIL sequence and IWGSC RefSeq v1.1 gene models. The presence of SNPs and other sequence variations in these genes between pau16057 and WL711 suggested that they are candidates for Lr76 and Yr70.

Divergence between bread wheat and Triticum militinae in the powdery mildew resistance QPm.tut-4A locus and its implications for cloning of the resistance gene.

A segment of Triticum militinae chromosome 7G harbors a gene(s) conferring powdery mildew resistance which is effective at both the seedling and the adult plant stages when transferred into bread wheat (T. aestivum). The introgressed segment replaces a piece of wheat chromosome arm 4AL. An analysis of segregating materials generated to positionally clone the gene highlighted that in a plant heterozygous for the introgression segment, only limited recombination occurs between the introgressed region and bread wheat 4A. Nevertheless, 75 genetic markers were successfully placed within the region, thereby confining the gene to a 0.012 cM window along the 4AL arm. In a background lacking the Ph1 locus, the localized rate of recombination was raised 33-fold, enabling the reduction in the length of the region containing the resistance gene to a 480 kbp stretch harboring 12 predicted genes. The substituted segment in the reference sequence of bread wheat cv. Chinese Spring is longer (640 kbp) and harbors 16 genes. A comparison of the segments' sequences revealed a high degree of divergence with respect to both their gene content and nucleotide sequence. Of the 12 T. militinae genes, only four have a homolog in cv. Chinese Spring. Possible candidate genes for the resistance have been identified based on function predicted from their sequence.

The in silico identification and characterization of a bread wheat/Triticum militinae introgression line.

The capacity of the bread wheat (Triticum aestivum) genome to tolerate introgression from related genomes can be exploited for wheat improvement. A resistance to powdery mildew expressed by a derivative of the cross-bread wheat cv. Tähti × T. militinae (Tm) is known to be due to the incorporation of a Tm segment into the long arm of chromosome 4A. Here, a newly developed in silico method termed rearrangement identification and characterization (RICh) has been applied to characterize the introgression. A virtual gene order, assembled using the GenomeZipper approach, was obtained for the native copy of chromosome 4A; it incorporated 570 4A DArTseq markers to produce a zipper comprising 2132 loci. A comparison between the native and introgressed forms of the 4AL chromosome arm showed that the introgressed region is located at the distal part of the arm. The Tm segment, derived from chromosome 7G, harbours 131 homoeologs of the 357 genes present on the corresponding region of Chinese Spring 4AL. The estimated number of Tm genes transferred along with the disease resistance gene was 169. Characterizing the introgression's position, gene content and internal gene order should not only facilitate gene isolation, but may also be informative with respect to chromatin structure and behaviour studies.

The wheat Phs-A1 pre-harvest sprouting resistance locus delays the rate of seed dormancy loss and maps 0.3 cM distal to the PM19 genes in UK germplasm.

The precocious germination of cereal grains before harvest, also known as pre-harvest sprouting, is an important source of yield and quality loss in cereal production. Pre-harvest sprouting is a complex grain defect and is becoming an increasing challenge due to changing climate patterns. Resistance to sprouting is multi-genic, although a significant proportion of the sprouting variation in modern wheat cultivars is controlled by a few major quantitative trait loci, including Phs-A1 in chromosome arm 4AL. Despite its importance, little is known about the physiological basis and the gene(s) underlying this important locus. In this study, we characterized Phs-A1 and show that it confers resistance to sprouting damage by affecting the rate of dormancy loss during dry seed after-ripening. We show Phs-A1 to be effective even when seeds develop at low temperature (13 °C). Comparative analysis of syntenic Phs-A1 intervals in wheat and Brachypodium uncovered ten orthologous genes, including the Plasma Membrane 19 genes (PM19-A1 and PM19-A2) previously proposed as the main candidates for this locus. However, high-resolution fine-mapping in two bi-parental UK mapping populations delimited Phs-A1 to an interval 0.3 cM distal to the PM19 genes. This study suggests the possibility that more than one causal gene underlies this major pre-harvest sprouting locus. The information and resources reported in this study will help test this hypothesis across a wider set of germplasm and will be of importance for breeding more sprouting resilient wheat varieties.

Characterization of repetitive DNA landscape in wheat homeologous group 4 chromosomes.

BACKGROUND: The number and complexity of repetitive elements varies between species, being in general most represented in those with larger genomes. Combining the flow-sorted chromosome arms approach to genome analysis with second generation DNA sequencing technologies provides a unique opportunity to study the repetitive portion of each chromosome, enabling comparisons among them. Additionally, different sequencing approaches may produce different depth of insight to repeatome content and structure. In this work we analyze and characterize the repetitive sequences of Triticum aestivum cv. Chinese Spring homeologous group 4 chromosome arms, obtained through Roche 454 and Illumina sequencing technologies, hereinafter marked by subscripts 454 and I, respectively. Repetitive sequences were identified with the RepeatMasker software using the interspersed repeat database mips-REdat_v9.0p. The input sequences consisted of our 4DS454 and 4DL454 scaffolds and 4ASI, 4ALI, 4BSI, 4BLI, 4DSI and 4DLI contigs, downloaded from the International Wheat Genome Sequencing Consortium (IWGSC). RESULTS: Repetitive sequences content varied from 55% to 63% for all chromosome arm assemblies except for 4DLI, in which the repeat content was 38%. Transposable elements, small RNA, satellites, simple repeats and low complexity sequences were analyzed. SSR frequency was found one per 24 to 27 kb for all chromosome assemblies except 4DLI, where it was three times higher. Dinucleotides and trinucleotides were the most abundant SSR repeat units. (GA)n/(TC)n was the most abundant SSR except for 4DLI where the most frequently identified SSR was (CCG/CGG)n. Retrotransposons followed by DNA transposons were the most highly represented sequence repeats, mainly composed of CACTA/En-Spm and Gypsy superfamilies, respectively. This whole chromosome sequence analysis allowed identification of three new LTR retrotransposon families belonging to the Copia superfamily, one belonging to the Gypsy superfamily and two TRIM retrotransposon families. Their physical distribution in wheat genome was analyzed by fluorescent in situ hybridization (FISH) and one of them, the Carmen retrotransposon, was found specific for centromeric regions of all wheat chromosomes. CONCLUSION: The presented work is the first deep report of wheat repetitive sequences analyzed at the chromosome arm level, revealing the first insight into the repeatome of T. aestivum chromosomes of homeologous group 4.

Chromosomal genomics facilitates fine mapping of a Russian wheat aphid resistance gene.

KEY MESSAGE: Making use of wheat chromosomal resources, we developed 11 gene-associated markers for the region of interest, which allowed reducing gene interval and spanning it by four BAC clones. Positional gene cloning and targeted marker development in bread wheat are hampered by high complexity and polyploidy of its nuclear genome. Aiming to clone a Russian wheat aphid resistance gene Dn2401 located on wheat chromosome arm 7DS, we have developed a strategy overcoming problems due to polyploidy and enabling efficient development of gene-associated markers from the region of interest. We employed information gathered by GenomeZipper, a synteny-based tool combining sequence data of rice, Brachypodium, sorghum and barley, and took advantage of a high-density linkage map of Aegilops tauschii. To ensure genome- and locus-specificity of markers, we made use of survey sequence assemblies of isolated wheat chromosomes 7A, 7B and 7D. Despite the low level of polymorphism of the wheat D subgenome, our approach allowed us to add in an efficient and cost-effective manner 11 new gene-associated markers in the Dn2401 region and narrow down the target interval to 0.83 cM. Screening 7DS-specific BAC library with the flanking markers revealed a contig of four BAC clones that span the Dn2401 region in wheat cultivar 'Chinese Spring'. With the availability of sequence assemblies and GenomeZippers for each of the wheat chromosome arms, the proposed strategy can be applied for focused marker development in any region of the wheat genome.

Molecular mapping of stripe rust resistance gene Yr51 in chromosome 4AL of wheat.

KEY MESSAGE: This manuscript describes the chromosomal location of a new source of stripe rust resistance in wheat. DNA markers closely linked with the resistance locus were identified and validated. A wheat landrace, AUS27858, from the Watkins collection showed high levels of resistance against Australian pathotypes of Puccinia striiformis f. sp. tritici. It was reported to carry two genes for stripe rust resistance, tentatively named YrAW1 and YrAW2. One hundred seeds of an F3 line (HSB#5515; YrAW1yrAW1) that showed monogenic segregation for stripe rust response were sown and harvested individually to generate monogenically segregating population (MSP) #5515. Stripe rust response variation in MSP#5515 conformed to segregation at a single locus. Bulked segregant analysis using high-throughput DArT markers placed YrAW1 in chromosome 4AL. MSP#5515 was advanced to F6 and phenotyped for detailed mapping. Novel wheat genomic resources including chromosome-specific sequence and genome zipper were employed to develop markers specific for the long arm of chromosome 4A. These markers were used for further saturation of the YrAW1 carrying region. YrAW1 was delimited by 3.7 cM between markers owm45F3R3 and sun104. Since there was no other stripe rust resistance gene located in chromosome 4AL, YrAW1 was formally named Yr51. Reference stock for Yr51 was lodged at the Australian Winter Cereal Collection, Tamworth, Australia and it was accessioned as AUS91456. Marker sun104 was genotyped on a set of Australian and Indian wheat cultivars and was shown to lack the resistance-linked sun104-225 bp allele. Marker sun104 is currently being used for marker-assisted backcrossing of Yr51 in Australian and Indian wheat backgrounds.

Next-generation sequencing and syntenic integration of flow-sorted arms of wheat chromosome 4A exposes the chromosome structure and gene content.

Wheat is the third most important crop for human nutrition in the world. The availability of high-resolution genetic and physical maps and ultimately a complete genome sequence holds great promise for breeding improved varieties to cope with increasing food demand under the conditions of changing global climate. However, the large size of the bread wheat (Triticum aestivum) genome (approximately 17 Gb/1C) and the triplication of genic sequence resulting from its hexaploid status have impeded genome sequencing of this important crop species. Here we describe the use of mitotic chromosome flow sorting to separately purify and then shotgun-sequence a pair of telocentric chromosomes that together form chromosome 4A (856 Mb/1C) of wheat. The isolation of this much reduced template and the consequent avoidance of the problem of sequence duplication, in conjunction with synteny-based comparisons with other grass genomes, have facilitated construction of an ordered gene map of chromosome 4A, embracing ≥85% of its total gene content, and have enabled precise localization of the various translocation and inversion breakpoints on chromosome 4A that differentiate it from its progenitor chromosome in the A genome diploid donor. The gene map of chromosome 4A, together with the emerging sequences of homoeologous wheat chromosome groups 4, 5 and 7, represent unique resources that will allow us to obtain new insights into the evolutionary dynamics between homoeologous chromosomes and syntenic chromosomal regions.

The chromatin determinants and Ph1 gene effect at wheat sites with contrasting recombination frequency.

INTRODUCTION: Meiotic recombination is one of the most important processes of evolution and adaptation to environmental conditions. Even though there is substantial knowledge about proteins involved in the process, targeting specific DNA loci by the recombination machinery is not well understood. OBJECTIVES: This study aims to investigate a wheat recombination hotspot (H1) in comparison with a "regular" recombination site (Rec7) on the sequence and epigenetic level in conditions with functional and non-functional Ph1 locus. METHODS: The DNA sequence, methylation pattern, and recombination frequency were analyzed for the H1 and Rec7 in three mapping populations derived by crossing introgressive wheat line 8.1 with cv. Chinese Spring (with Ph1 and ph1 alleles) and cv. Tähti. RESULTS: The H1 and Rec7 loci are 1.586 kb and 2.538 kb long, respectively. High-density mapping allowed to delimit the Rec7 and H1 to 19 and 574 bp and 593 and 571 bp CO sites, respectively. A new method (ddPing) allowed screening recombination frequency in almost 66 thousand gametes. The screening revealed a 5.94-fold higher recombination frequency at the H1 compared to the Rec7. The H1 was also found out of the Ph1 control, similarly as gamete distortion. The recombination was strongly affected by larger genomic rearrangements but not by the SNP proximity. Moreover, chromatin markers for open chromatin and DNA hypomethylation were found associated with crossover occurrence except for the CHH methylation. CONCLUSION: Our results, for the first time, allowed study of wheat recombination directly on sequence, shed new light on chromatin landmarks associated with particular recombination sites, and deepened knowledge about role of the Ph1 locus in control of wheat recombination processes. The results are suggesting more than one recombination control pathway. Understanding this phenomenon may become a base for more efficient wheat genome manipulation, gene pool enrichment, breeding, and study processes of recombination itself.

Subgenomic analysis of microRNAs in polyploid wheat.

In this study, a survey of miRNAs using the next-generation sequencing data was performed at subgenomic level. After analyzing shotgun sequences from chromosome 4A of bread wheat (Triticum aestivum L.), a total of 68 different miRNAs were predicted in silico, of which 37 were identified in wheat for the first time. The long arm of the chromosome was found to harbor a higher variety (51) and representation (3,928) of miRNAs compared with the short arm (49; 2,226). Out of the 68 miRNAs, 32 were detected to be common to both arms, revealing the presence of separate miRNA clusters in the two chromosome arms. The differences in degree of representation of the different miRNAs were found to be highly variable, ranging 592-fold, which may have an effect on target regulation. Targets were retrieved for 62 (out of 68) of wheat-specific, newly identified miRNAs indicated that fundamental aspects of plant morphology such as height and flowering were predicted to be affected. In silico expression blast analysis indicated 24 (out of 68) were found to give hits to expressed sequences. This is the first report of species- and chromosome-specific miRNAs.

Fine mapping, phenotypic characterization and validation of non-race-specific resistance to powdery mildew in a wheat-Triticum militinae introgression line.

Introgression of several genomic loci from tetraploid Triticum militinae into bread wheat cv. Tähti has increased resistance of introgression line 8.1 to powdery mildew in seedlings and adult plants. In our previous work, only a major quantitative trait locus (QTL) on chromosome 4AL of the line 8.1 contributed significantly to resistance, whereas QTL on chromosomes 1A, 1B, 2A, 5A and 5B were detected merely on a suggestive level. To verify and characterize all QTLs in the line 8.1, a mapping population of double haploid lines was established. Testing for seedling resistance to 16 different races/mixtures of Blumeria graminis f. sp. tritici revealed four highly significant non-race-specific resistance QTL including the main QTL on chromosome 4AL, and a race-specific QTL on chromosome 5B. The major QTL on chromosome 4AL (QPm.tut-4A) as well as QTL on chromosome 5AL and a newly detected QTL on 7AL were highly effective at the adult stage. The QPm.tut-4A QTL accounts on average for 33-49 % of the variation in resistance in the double haploid population. Interactions between the main QTL QPm.tut-4A and the minor QTL were evaluated and discussed. A population of 98 F(2) plants from a cross of susceptible cv. Chinese Spring and the line 8.1 was created that allowed mapping the QPm.tut-4A locus to the proximal 2.5-cM region of the introgressed segment on chromosome 4AL. The results obtained in this work make it feasible to use QPm.tut-4A in resistance breeding and provide a solid basis for positional cloning of the major QTL.

Identification, High-Density Mapping, and Characterization of New Major Powdery Mildew Resistance Loci From the Emmer Wheat Landrace GZ1.

Powdery mildew is one of the most devastating diseases of wheat which significantly decreases yield and quality. Identification of new sources of resistance and their implementation in breeding programs is the most effective way of disease control. Two major powdery mildew resistance loci conferring resistance to all races in seedling and adult plant stages were identified in the emmer wheat landrace GZ1. Their positions, effects, and transferability were verified using two linkage maps (1,510 codominant SNP markers) constructed from two mapping populations (276 lines in total) based on the resistant GZ1 line. The dominant resistance locus QPm.GZ1-7A was located in a 90 cM interval of chromosome 7AL and explains up to 20% of the trait variation. The recessive locus QPm.GZ1-2A, which provides total resistance, explains up to 40% of the trait variation and was located in the distal part of chromosome 2AL. The locus was saturated with 14 PCR-based markers and delimited to a 0.99 cM region which corresponds to 4.3 Mb of the cv. Zavitan reference genome and comprises 55 predicted genes with no apparent candidate for the QPm.GZ1-2A resistance gene. No recessive resistance gene or allele was located at the locus before, suggesting the presence of a new powdery mildew resistance gene in the GZ1. The mapping data and markers could be used for the implementation of the locus in breeding. Moreover, they are an ideal base for cloning and study of host-pathogen interaction pathways determined by the resistance genes.

A multi gene sequence-based phylogeny of the Musaceae (banana) family.

BACKGROUND: The classification of the Musaceae (banana) family species and their phylogenetic inter-relationships remain controversial, in part due to limited nucleotide information to complement the morphological and physiological characters. In this work the evolutionary relationships within the Musaceae family were studied using 13 species and DNA sequences obtained from a set of 19 unlinked nuclear genes. RESULTS: The 19 gene sequences represented a sample of ~16 kb of genome sequence (~73% intronic). The sequence data were also used to obtain estimates for the divergence times of the Musaceae genera and Musa sections. Nucleotide variation within the sample confirmed the close relationship of Australimusa and Callimusa sections and showed that Eumusa and Rhodochlamys sections are not reciprocally monophyletic, which supports the previous claims for the merger between the two latter sections. Divergence time analysis supported the previous dating of the Musaceae crown age to the Cretaceous/Tertiary boundary (~ 69 Mya), and the evolution of Musa to ~50 Mya. The first estimates for the divergence times of the four Musa sections were also obtained. CONCLUSIONS: The gene sequence-based phylogeny presented here provides a substantial insight into the course of speciation within the Musaceae. An understanding of the main phylogenetic relationships between banana species will help to fine-tune the taxonomy of Musaceae.

Control of flowering time and spike development in cereals: the earliness per se Eps-1 region in wheat, rice, and Brachypodium.

The earliness per se gene Eps-A ( m ) 1 from diploid wheat Triticum monococcum affects heading time, spike development, and spikelet number. In this study, the Eps1 orthologous regions from rice, Aegilops tauschii, and Brachypodium distachyon were compared as part of current efforts to clone this gene. A single Brachypodium BAC clone spanned the Eps-A ( m ) 1 region, but a gap was detected in the A. tauschii physical map. Sequencing of the Brachypodium and A. tauschii BAC clones revealed three genes shared by the three species, which showed higher identity between wheat and Brachypodium than between them and rice. However, most of the structural changes were detected in the wheat lineage. These included an inversion encompassing the wg241-VatpC region and the presence of six unique genes. In contrast, only one unique gene (and one pseudogene) was found in Brachypodium and none in rice. Three genes were present in both Brachypodium and wheat but were absent in rice. Two of these genes, Mot1 and FtsH4, were completely linked to the earliness per se phenotype in the T. monococcum high-density genetic map and are candidates for Eps-A ( m ) 1. Both genes were expressed in apices and developing spikes, as expected for Eps-A ( m ) 1 candidates. The predicted MOT1 protein showed amino acid differences between the parental T. monococcum lines, but its effect is difficult to predict. Future steps to clone the Eps-A ( m ) 1 gene include the generation of mot1 and ftsh4 mutants and the completion of the T. monococcum physical map to test for the presence of additional candidate genes.

Fine Mapping of Lr49 Using 90K SNP Chip Array and Flow-Sorted Chromosome Sequencing in Wheat.

Leaf rust, caused by Puccinia triticina, threatens global wheat production due to the constant evolution of virulent pathotypes that defeat commercially deployed all stage-resistance (ASR) genes in modern cultivars. Hence, the deployment of combinations of adult plant resistance (APR) and ASR genes in new wheat cultivars is desirable. Adult plant resistance gene Lr49 was previously mapped on the long arm of chromosome 4B of cultivar VL404 and flanked by microsatellite markers barc163 (8.1 cM) and wmc349 (10.1 cM), neither of which was sufficiently closely linked for efficient marker assisted selection. This study used high-density SNP genotyping and flow sorted chromosome sequencing to fine-map the Lr49 locus as a starting point to develop a diagnostic marker for use in breeding and to clone this gene. Marker sunKASP_21 was mapped 0.4 cM proximal to Lr49, whereas a group of markers including sunKASP_24 were placed 0.6 cM distal to this gene. Testing of the linked markers on 75 Australian and 90 European cultivars with diverse genetic backgrounds showed that sunKASP_21 was most strongly associated with Lr49. Our results also show that the Lr49 genomic region contains structural variation relative to the reference stock Chinese Spring, possibly an inverted genomic duplication, which introduces a new set of challenges for the Lr49 cloning.

Identification of COS markers specific for Thinopyrum elongatum chromosomes preliminary revealed high level of macrosyntenic relationship between the wheat and Th. elongatum genomes.

Thinopyrum elongatum (Host) D.R. Dewey has served as an important gene source for wheat breeding improvement for many years. The exact characterization of its chromosomes is important for the detailed analysis of prebreeding materials produced with this species. The major aim of this study was to identify and characterize new molecular markers to be used for the rapid analysis of E genome chromatin in wheat background. Sixty of the 169 conserved orthologous set (COS) markers tested on diverse wheat-Th. elongatum disomic/ditelosomic addition lines were assigned to various Th. elongatum chromosomes and will be used for marker-assisted selection. The macrosyntenic relationship between the wheat and Th. elongatum genomes was investigated using EST sequences. Several rearrangements were revealed in homoeologous chromosome groups 2, 5, 6 and 7, while chromosomes 1 and 4 were conserved. Molecular cytogenetic and marker analysis showed the presence of rearranged chromosome involved in 6ES and 2EL arms in the 6E disomic addition line. The selected chromosome arm-specific COS markers will make it possible to identify gene introgressions in breeding programmes and will also be useful in the development of new chromosome-specific markers, evolutionary analysis and gene mapping.

Correction: Heritable heading time variation in wheat lines with the same number of Ppd-B1 gene copies.

[This corrects the article DOI: 10.1371/journal.pone.0183745.].

Heritable heading time variation in wheat lines with the same number of Ppd-B1 gene copies.

The ability of plants to identify an optimal flowering time is critical for ensuring the production of viable seeds. The main environmental factors that influence the flowering time include the ambient temperature and day length. In wheat, the ability to assess the day length is controlled by photoperiod (Ppd) genes. Due to its allohexaploid nature, bread wheat carries the following three Ppd-1 genes: Ppd-A1, Ppd-B1 and Ppd-D1. While photoperiod (in)sensitivity controlled by Ppd-A1 and Ppd-D1 is mainly determined by sequence changes in the promoter region, the impact of the Ppd-B1 alleles on the heading time has been linked to changes in the copy numbers (and possibly their methylation status) and sequence changes in the promoter region. Here, we report that plants with the same number of Ppd-B1 copies may have different heading times. Differences were observed among F7 lines derived from crossing two spring hexaploid wheat varieties. Several lines carrying three copies of Ppd-B1 headed 16 days later than other plants in the population with the same number of gene copies. This effect was associated with changes in the gene expression level and methylation of the Ppd-B1 gene.

Haplotype Analysis of the Pre-harvest Sprouting Resistance Locus Phs-A1 Reveals a Causal Role of TaMKK3-A in Global Germplasm.

Pre-harvest sprouting (PHS) is an important cause of quality loss in many cereal crops and is particularly prevalent and damaging in wheat. Resistance to PHS is therefore a valuable target trait in many breeding programs. The Phs-A1 locus on wheat chromosome arm 4AL has been consistently shown to account for a significant proportion of natural variation to PHS in diverse mapping populations. However, the deployment of sprouting resistance is confounded by the fact that different candidate genes, including the tandem duplicated Plasma Membrane 19 (PM19) genes and the mitogen-activated protein kinase kinase 3 (TaMKK3-A) gene, have been proposed to underlie Phs-A1. To further define the Phs-A1 locus, we constructed a physical map across this interval in hexaploid and tetraploid wheat. We established close proximity of the proposed candidate genes which are located within a 1.2 Mb interval. Genetic characterization of diverse germplasm used in previous genetic mapping studies suggests that TaMKK3-A, and not PM19, is the major gene underlying the Phs-A1 effect in European, North American, Australian and Asian germplasm. We identified the non-dormant TaMKK3-A allele at low frequencies within the A-genome diploid progenitor Triticum urartu genepool, and show an increase in the allele frequency in modern varieties. In United Kingdom varieties, the frequency of the dormant TaMKK3-A allele was significantly higher in bread-making quality varieties compared to feed and biscuit-making cultivars. Analysis of exome capture data from 58 diverse hexaploid wheat accessions identified fourteen haplotypes across the extended Phs-A1 locus and four haplotypes for TaMKK3-A. Analysis of these haplotypes in a collection of United Kingdom and Australian cultivars revealed distinct major dormant and non-dormant Phs-A1 haplotypes in each country, which were either rare or absent in the opposing germplasm set. The diagnostic markers and haplotype information reported in the study will help inform the choice of germplasm and breeding strategies for the deployment of Phs-A1 resistance into breeding germplasm.