- 1Nottingham BBSRC Wheat Research Centre, School of Biosciences, University of Nottingham, Loughborough, United Kingdom
- 2Earlham Institute, Norwich, United Kingdom
Triticum timopheevii (2n = 28, AtAtGG) is a tetraploid wild relative species with great potential to increase the genetic diversity of hexaploid wheat Triticum aestivum (2n = 42, AABBDD) for various important agronomic traits. A breeding scheme that propagated advanced backcrossed populations of wheat-T. timopheevii introgression lines through further backcrossing and self-fertilisation resulted in the generation of 99 introgression lines (ILs) that carried 309 homozygous segments from the At and G subgenomes of T. timopheevii. These introgressions contained 89 and 74 unique segments from the At and G subgenomes, respectively. These overlapping segments covered 98.9% of the T. timopheevii genome that has now been introgressed into bread wheat cv. Paragon including the entirety of all T. timopheevii chromosomes via varying sized segments except for chromosomes 3At, 4G, and 6G. Homozygous ILs contained between one and eight of these introgressions with an average of three per introgression line. These homozygous introgressions were detected through the development of a set of 480 chromosome-specific Kompetitive allele specific PCR (KASP) markers that are well-distributed across the wheat genome. Of these, 149 were developed in this study based on single nucleotide polymorphisms (SNPs) discovered through whole genome sequencing of T. timopheevii. A majority of these KASP markers were also found to be T. timopheevii subgenome specific with 182 detecting At subgenome and 275 detecting G subgenome segments. These markers showed that 98% of the At segments had recombined with the A genome of wheat and 74% of the G genome segments had recombined with the B genome of wheat with the rest recombining with the D genome of wheat. These results were validated through multi-colour in situ hybridisation analysis. Together these homozygous wheat-T. timopheevii ILs and chromosome-specific KASP markers provide an invaluable resource to wheat breeders for trait discovery to combat biotic and abiotic stress factors affecting wheat production due to climate change.
Introduction
The wild relatives of wheat provide genetic variation for most, if not all, traits of agronomic importance including resistances to abiotic stresses, e.g., Thinopyrum bessarabicum (King et al., 1997), biotic stresses, e.g., Aegilops speltoides (Niranjana et al., 2017), and Aegilops ventricosa (Cruz et al., 2016), novel variation for protein quality, e.g., Thinopyrum intermedium (Li et al., 2013) and mineral content, e.g., Aegilops sp. (Sharma et al., 2018) and increases in physiological traits such as photosynthetic capacity, e.g., Thinopyrum elongatum (Reynolds et al., 2001).
A wild relative that has proved to be of considerable interest due to the high levels of genetic variation it has been shown to contain is Triticum timopheevii subsp. timopheevii (Zhuk.) (2n = 4x = 28; AtAtGG). Traits of interest include resistance to leaf rust including the genes Lr10, Lr18, and Lr50 (McIntosh, 1983; Brown-Guedira et al., 2003; Leonova et al., 2010; Singh et al., 2017), stem rust including the genes Sr36, Sr37, and Sr40 (Allard and Shands, 1954; McIntosh and Gyarfas, 1971; Dyck, 1992), powdery mildew including the genes Pm2, Pm6, Pm27, and Pm37 (Allard and Shands, 1954; Jørgensen and Jensen, 1972; Peusha et al., 1995; Järve et al., 2000; Perugini et al., 2008; Qin et al., 2011), Fusarium head blight (Malihipour et al., 2016, 2017), black-point (Lehmensiek et al., 2004), Hessian fly, Septoria blotch, wheat curl mite and tan spot (Brown-Guedira et al., 1996). Triticum timopheevii has also been shown to be a source of resistance to abiotic stresses such as salt tolerance (Badridze et al., 2009; Yudina et al., 2016) and traits affecting grain quality such as milling yield and grain protein (Lehmensiek et al., 2008) and mineral content (Hu et al., 2017).
Triticum timopheevii belongs to the secondary gene pool of wheat as both genomes are very closely related to wheat although they show modifications. The A-genomes of durum wheat (Triticum durum Desf., 2n = 4x = 28; AABB), bread wheat (Triticum aestivum L., 2n = 6x = 42; AABBDD) and T. timopheevii are derived from Triticum urartu Thum. ex Gandil (Dvořák et al., 1993) while the B-genomes of durum and bread wheat and the G-genome of T. timopheevii are derived from an Ae. speltoides-like species (Dvorak and Zhang, 1990). A recent study showed that Ae. speltoides was not the direct donor of the G subgenome and diverged from the progenitor of the latter ca. 2.93 mya, although it was more closely related to it than to the B genome of wheat (Li et al., 2022). As in T. durum and T. aestivum, T. timopheevii has been shown to carry the 4AtL/5AtL translocation, but in addition, also contains a 6At/1G/4G cyclic translocation (Jiang and Gill, 1994; Maestra and Naranjo, 1999; Salina et al., 2006).
Chromosome pairing studies have shown that the A and At genomes are more closely related to each other than are the B and G genomes with A and At homoeologous pairing approximately 70% of the time compared to B–G homoeologues which pair approximately only 30% of the time (Feldman, 1966). The ease with which the genetic variation in a wild species can be accessed is partly dependent on how closely related the genome(s) of the wild relative is to that of wheat, as the closer the relationship the higher the levels of homoeologous recombination possible. Hence, the pairing studies above would suggest that when crossing T. timopheevii to wheat, introgression of the At genome into the A genome is more likely to occur as compared to introgression of the G genome into the B genome.
A comprehensive set of molecular markers, polymorphic between wheat and T. timopheevii, is essential both during the generation of introgressions (for identification and characterisation) and also for the use of the introgressions in a breeding programme (in order to track the introgressions through the different generations). We earlier reported the generation of T. timopheevii introgression lines (Devi et al., 2019) with introgression identification and characterisation based on the Axiom® Wheat Relative Genotyping Array (King et al., 2017). The array was extremely valuable for the initial screening of a large number of putative introgression lines but did not provide the flexibility required for screening different generations. A set of wheat chromosome-specific KASP markers was therefore generated for use on wheat/wild relative introgressions produced from 10 different wild relatives including T. timopheevii (Grewal et al., 2020a). These KASP markers were used to characterise 91 lines of wheat carrying homozygous introgressions from T. timopheevii. However, with the close relationship between the genomes of T. timopheevii and wheat, a relatively high level of recombination was expected, with the possibility of numerous smaller introgressions. Further KASP markers, polymorphic between T. timopheevii and bread wheat, have therefore been generated to provide a greater density of markers and to increase the detection power of the marker system.
Materials and Methods
Plant Material
Advanced back-crossed populations of wheat-T. timopheevii introgression lines, previously generated by Devi et al. (2019), were used as the starting material to generate homozygous ILs. Devi et al. (2019) crossed hexaploid wheat cv. Paragon ph1/ph1 mutant (2n = 2x = 14) [obtained from the Germplasm Resource Unit (GRU) at the John Innes Centre] with T. timopheevii, (accession PI 94760, obtained from the United States Department of Agriculture, USDA; 2n = 4x = 28), to produce F1 hybrids which were backcrossed with Paragon to generate BC1 populations. The BC1 individuals were recurrently backcrossed to produce BC2, BC3, and BC4 populations and self-fertilised at each generation.
Seeds were germinated in soil in seedling trays. Leaf tissues (1.5-inch leaf segment cut into pieces) were harvested from 3-week-old plants, immediately frozen on liquid nitrogen and stored at −80°C until nucleic acid extraction. Seedlings were vernalised for 4 weeks at 4°C after which the plants were transferred to 2L pots containing John Innes No. 2 soil, irrigated and fertilized using a drip irrigation system and maintained at 18–25°C under 16 h light and 8 h dark conditions. Plants were either self-fertilised or hand pollinated and crossed with the recurrent parent, wheat cv. Paragon, as the male parent.
Breeding Scheme
The breeding scheme used to develop the homozygous wheat-T. timopheevii ILs is shown in Figure 1. A set of advanced backcrossed families were obtained through a previous study (Devi et al., 2019) and genotyped with KASP markers. Plants from different backcross families, carrying overlapping T. timopheevii segments from each of the 14 linkage groups were selected for further backcrossing and self-fertilisation. In the next generation, six to eight progenies from each selected plant were genotyped and those containing heterozygous segments of interest but still carrying two or more other (heterozygous or homozygous) T. timopheevii segments were selected for further backcrossing and self-fertilisation. Plants with potentially one or two homozygous T. timopheevii introgressions were self-fertilised and selected as homozygous introgression line (ILs). This process of plant selection through genotyping of segregating backcross and self-fertilised families was repeated for another generation of plants. At this point any plants showing one or more homozygous T. timopheevii segments were selected as homozygous ILs.
Figure 1. Breeding scheme used to generate homozygous wheat-Triticum timopheevii introgression lines (ILs) from advanced backcrossed lines generated by Devi et al. (2019). For each year (one generation per year), the number of families, Kompetitive allele specific PCR (KASP) marker sets assessed, selected plants, and type of cross are shown. The backcrosses with the recurrent parents are represented by an × and a wheat plant, and the self-fertilisation rounds by an encircled ×. In each year a set of homozygous ILs were selected for additional rounds of self-fertilisation to create a panel of plants carrying overlapping T. timopheevii segments of varying sizes.
DNA Extraction
Leaf tissue samples were collected and freeze-dried in a deep-well plate and ground in a TissueLyser II (QIAGEN) for 4–6 mins at a frequency of 25 Hz. Genomic DNA for sequencing and genotyping was extracted as described in Grewal et al. (2022).
Genomic DNA extraction for generation of probes for genomic in situ hybridisation analysis, was carried using the same protocol with an additional step of purification with phenol/chloroform at the end.
Development of Kompetitive Allele Specific PCR Markers
Chromosome-Specific Single Nucleotide Polymorphism Discovery
DNA was isolated from T. timopheevii accession P95-99.1-1, as described above, and a PCR-free library was prepared and sequenced on Illumina HiSeq 2500 to produce 250 bp paired-end reads at a coverage of ∼5x. SNP discovery followed a bioinformatics pipeline that produced a set of SNPs suitable for KASP assay design. In detail, the reads were mapped to the wheat cv. Chinese Spring reference genome assembly RefSeq v1.0 (International Wheat Genome Sequencing Consortium [IWGSC], et al., 2018) using BWA MEM version 0.7.13 (Li, 2013) with the -M flag. The resulting alignments were filtered using SAMtools v1.4 (Li et al., 2009) to remove unmapped reads, supplementary alignments, improperly paired reads, and non-uniquely mapping reads (q < 10). PCR duplicates were removed using Picard’s MarkDuplicates (DePristo et al., 2011). SNPs were called using BCFtools (Li, 2011) using the multiallelic model (-m). To remove SNPs unsuitable for KASP assays, we retained homozygous SNPs with a read depth greater than 10 and 100% mapped reads supporting the alternative allele and each position 50 bp up and downstream of the SNP site having no INDELs, containing greater than five reads and ≥80% of reads supporting the reference allele. This aims to remove sites with nearby variants that could disrupt the KASP assay. To prevent amplification of off-target (not the target wheat chromosome) regions in the genome, the remaining SNP sites, along with 50 bp up and downstream were aligned to RefSeq v1.0 using BLASTn (Camacho et al., 2009). SNP sites were removed if in addition to the expected self-hit they produced an off-target hit with ≥97% identity.
Kompetitive Allele Specific PCR Assay Design
To design a KASP™ assay, the flanking sequence of a SNP was fed through the PolyMarker application (Ramirez-Gonzalez et al., 2015) which aligned the query sequence to wheat cv. Chinese Spring RefSeq v1.0 and provided two allele-specific primers and one common primer for each assay. A value of 1 in the “total_contigs” column of the output Primers file validated the query SNP to be in a unique sequence region in the wheat genome RefSeq v1.0 assembly. The position of the KASP markers on wheat chromosomes in shown in Figure 2 and was generated through RIdeogram (Hao et al., 2020).
Figure 2. A graphical representation of physical positions of the four sets of chromosome-specific KASP markers, across all wheat chromosomes, used in this study to characterise wheat- Triticum timopheevii introgression lines. Sets 1, 2 and 3 were developed in previous works (Grewal et al., 2020a,2021, 2022) and Set 4 was developed in this study to fill in the gaps between markers from the previous sets. Circles depict KASP markers that specifically detect At subgenome segments, squares depict markers that detect G subgenome segments and triangles represent markers that detect both subgenomes.
Genotyping With Kompetitive Allele Specific PCR Markers
For genotyping during the breeding pipeline, four sets of KASP markers were used. The first three were developed during previous studies (Grewal et al., 2020a,2021, 2022) and the fourth was developed during this work using chromosome-specific SNPs between wheat and T. timopheevii, discovered as described above.
Set 1 consisted of 250 chromosome-specific KASP markers previously reported to be polymorphic between wheat and T. timopheevii (codes between WRC0001-1000; Grewal et al., 2020a). KASP assays in Set 2 and 3 were originally designed to be polymorphic between wheat and other species and needed to be tested on T. timopheevii first to find those that were polymorphic between wheat and T. timopheevii. A subset of the KASP assays originally tested in those studies were selected for Set 2 and 3: those that did not fail at the PCR stage. Set 2 consisted of KASP assays designed to be tested on doubled haploid wheat-T. urartu ILs (codes between WRC1080-1308 and WRC1317-1393; Grewal et al., 2021) and Set 3 consisted of KASP assays designed to be tested on doubled haploid wheat-Amblyopyrum muticum ILs (codes between WRC1309-1316, WRC1394-1713, WRC1723-1872, WRC1894-1913, WRC1954-2113, and WRC2130-2169; Grewal et al., 2022). Set 4 consisted of the new KASP markers that were designed in this study (codes between WRC2170-2341 and WRC2366-2401).
The genotyping procedure was as described by Grewal et al. (2020b). Briefly, the genotyping reactions were set up using the automated PIPETMAX® 268 (Gilson, United Kingdom) and performed in a ProFlex PCR system (Applied Biosystems by Life Technology) in a final volume of 5 μl with 1 ng genomic DNA, 2.5 μl KASP reaction mix (ROX), 0.068 μl primer mix and 2.43 μl nuclease free water. PCR conditions were set as 15 min at 94°C; 10 touchdown cycles of 10 s at 94°C, 1 min at 65–57°C (dropping 0.8°C per cycle); and 35 cycles of 15 s at 94°C and 1 min at 57°C. Fluorescence detection of the reactions was performed using a QuantStudio 5 (Applied Biosystems, Waltham, MA, United States) and the data analysed using the QuantStudio™ Design and Analysis Software V1.5.0 (Applied Biosystems, Waltham, MA, United States).
ChromoMap (Anand and Rodriguez Lopez, 2022) was used to annotate each unique T. timopheevii segment introgressed into wheat and RIdeogram (Hao et al., 2020) was used to annotate the KASP marker positions on the wheat genome.
in situ Hybridisation
Preparation of the root-tip metaphase chromosome spreads for in situ hybridisation analysis was as described in Grewal et al. (2017) and King et al. (2017) with a few modifications. Briefly, two roots (2 cm) from a 2-day old seedling were excised and treated with nitrous oxide gas at 10 bar for 2 h. Treated roots were fixed in 90% acetic acid for 15 min and then washed three times with water, on ice. The root tip was dissected and digested in 20 μL of 1% pectolyase Y23 and 2% cellulase Onozuka R-10 (Yakult Pharmaceutical, Tokyo, Japan) solution for 50 min at 37°C and then washed three times with 70% ethanol. The root tips were crushed in 70% ethanol, and the cells collected by centrifugation at 3,000 × g for 1 min, briefly dried and then re-suspended in 20–30 μL of 100% acetic acid before being placed on ice. The cell suspension was dropped onto glass slides (6–7 μL per slide) in a moist box and dried slowly under cover. Slide denaturation was performed at 80°C for 5 min and probes were hybridised at 55°C for 16 h. Preparation of probes for genomic and fluorescent in situ hybridisation was as described below.
Post-hybridisation, all slides were washed with 2 × saline-sodium citrate (SSC) buffer for 1 min and counter-stained with DAPI (4′,6-diamidino-2-phenylindole). Metaphases were detected using a high-throughput, fully automated Zeiss Axio ImagerZ2 upright epifluorescence microscope (Carl Zeiss Ltd., Germany). Image capture was performed using a MetaSystems Coolcube 1m CCD camera and image analysis was carried out using Metafer4 (automated metaphase image capture) and ISIS (image processing) software (Metasystems GmbH, Germany).
Multi-Colour Genomic in situ Hybridisation
Genomic DNA from T. urartu (to detect the A-genome), Ae. speltoides (to detect the B-genome), and Aegilops tauschii (to detect the D-genome) were isolated as described above and labelled to be used as probes. Two micrograms of genomic DNA was labelled by nick translation (Rigby et al., 1977) as follows: (1) T. urartu with ChromaTide™ Alexa Fluor™ 488-5-dUTP (Invitrogen; C11397; coloured green), (2) Ae. speltoides with DEAC-dUTP (Jena Bioscience; NU-803-DEAC; coloured blueish purple), and (3) Ae. tauschii with ChromaTide™ Alexa Fluor™ 594-5-dUTP (Invitrogen; C11400; coloured red). Slides were probed using 150 ng of T. urartu, 150 ng of Ae. speltoides and 300 ng of Ae. tauschii, in the ratio 1:1:2 (green: blue: red). No blocking DNA was used.
Multi-Colour Fluorescent in situ Hybridisation
Multi-colour fluorescent in situ hybridisation (Mc-FISH) was carried out as described by Grewal et al. (2018). The pSc119.2 clone (McIntyre et al., 1990), a 120-bp repetitive DNA sequence derived from rye (Secale cereale L.) and inserted into the plasmid pBR322, and the pAs.1 clone (Rayburn and Gill, 1986), a 1 kb repetitive DNA sequence isolated from Ae. tauschii and inserted into the pUC8 plasmid, were nick-labelled (Rigby et al., 1977) with Alexa Fluor 488-5-dUTP (green) and Alexa Fluor 594-5-dUTP (red), respectively. Classification of pSc119.2- and pAs.1-labelled chromosomes followed the nomenclature of Badaeva et al. (2016).
Results
Development of Homozygous Introgression Lines
Advanced backcrossed populations of wheat-T. timopheevii introgression lines developed by Devi et al. (2019), segregating for various T. timopheevii segments, were propagated through backcrossing and self-fertilisation and genotyped at every generation with chromosome-specific KASP markers (Grewal et al., 2020a,2021, 2022) as shown in Figure 1. Plants carrying segments of various sizes and from different linkage groups were selected to create a panel of overlapping T. timopheevii introgressions from each of its linkage groups in the wheat background.
In 2018, a previously developed set of chromosome-specific KASP markers (Grewal et al., 2020a) consisting of 250 markers polymorphic between wheat and T. timopheevii, and designated as Set 1, were used to genotype 1,685 advanced backcrossed lines (Devi et al., 2019) belonging to 494 BCnFn families (Figure 1). The genotyping identified 101 plants that carried unique T. timopheevii segments but since these lines also had multiple segments or whole chromosomes from T. timopheevii, they were taken forward to another generation either through backcrossing with wheat cv. Paragon or through self-pollination.
The next generation, segregating for the desired T. timopheevii segments and consisting of 772 plants from 101 BCnFn families, were genotyped with the Set 1 KASP markers in 2019. As shown in Figure 2, this set of markers had significant gaps on many chromosomes and thus, additional chromosome-specific KASP markers were sought to fill in those gaps. Simultaneously to this work, KASP markers were being developed to characterise doubled haploid wheat ILs from T. urartu and Am. muticum (Grewal et al., 2021, 2022). These two groups of chromosome-specific KASP assays consisting of 224 markers developed for T. urartu and 425 markers developed for Am. muticum were tested on T. timopheevii and its derived segregating BCnFn families. From the first group, 30 were found to be diagnostic for T. timopheevii, designated as Set 2, and from the second group 51 markers were validated for T. timopheevii and designated as Set 3 (Figure 2 and Table 1). After genotyping with the three sets of KASP markers, 35 plants were selected as carrying one or two homozygous T. timopheevii segments and self-fertilised to create stable introgression lines (ILs). An additional 60 plants were selected as carrying segments of interest and were further backcrossed and self-fertilised either because the segment(s) were still heterozygous or to segregate multiple segments through a family to create ILs with fewer T. timopheevii segments.
Table 1. Number of chromosome-specific markers across the different sets of markers tested, average and maximum distance between markers and Triticum timopheevii subgenome-specific markers per wheat chromosome.
In 2020, the next generation consisting of 356 plants from 60 BCnFn families (Figure 1) were genotyped with Set 1, 2, and 3 KASP markers. The aim was to have a KASP marker at every 60–70 Mbp across the A and B genomes of wheat for genotyping of these lines. However, as shown in Figure 2, there were still significant gaps between markers from Sets 1–3 on some A and B genome chromosomes.
To fill in these gaps, new KASP markers were developed through discovery of chromosome-specific SNPs between T. timopheevii and wheat. After aligning T. timopheevii sequence reads against the wheat Chinese Spring RefSeqv1.0 assembly (International Wheat Genome Sequencing Consortium [IWGSC], et al., 2018), variant calling and filtering for KASP assay suitable flanking regions, 36,001 SNPs were generated between wheat and T. timopheevii. SNPs with significant off-target hits were removed to give 6957 chromosome-specific SNPs (Supplementary Data 1). Of these 208 SNPs were selected at positions that filled the most prominent gaps between markers in Sets 1–3, converted to KASP assays and tested on the segregating families. The results of genotyping showed that 39 assays failed at the PCR stage, 20 assays were not robust in separating heterozygotes and homozygotes and thus, 149 assays were able to successfully detect T. timopheevii segments in a wheat background. The latter were designated as Set 4 (Figure 2 and Table 1).
Based on the genotyping results of markers from Sets 1–4, a further 39 ILs were selected as carrying only homozygous T. timopheevii segments (Figure 1). Along with the 35 ILs from the previous generation, all 74 ILs were self-fertilised for two further generations and genotyped with all 480 KASP markers from Sets 1–4. These lines were code named Tim1 to Tim74.
A further 26 plants were selected, from this last generation, that had a combination of heterozygous and homozygous desirable T. timopheevii segments. Their corresponding 26 BCnFn families were taken forward through the germination of 89 plants. After genotyping with all 480 KASP markers, 25 were selected as homozygous for all introgressions and are currently undergoing the first of two rounds of self-fertilisation. These lines were code named Tim75 to Tim99 and will be undergo a further two rounds of multiplication.
Distribution and Specificity of Kompetitive Allele Specific PCR Markers
In total 480 chromosome-specific KASP markers were used to genotype the wheat-T. timopheevii ILs (Table 1, Figure 2, and Supplementary Data 2). These markers were designed to be well-distributed across the wheat chromosomes. The average distance between markers ranged from 21.51 to 29.68 Mbp on the A genome chromosomes, 18.85–27.80 Mbp on the B genome chromosomes and 27.86–56.65 Mbp on the D genome chromosomes (Table 1). The latter had the largest distances between two markers but there were also some significant gaps between markers on chromosomes 4A, 5B, and 7A (Table 1 and Figure 2).
In addition to being able to only detect one wheat chromosome (Grewal et al., 2020a), the KASP markers developed for T. timopheevii were, in the majority of cases, also able to specifically detect either the At or the G subgenome chromosomes of T. timopheevii. Of the 480 KASP markers used to genotype the ILs, 182 were found to detect only At subgenome segments, 275 detected G subgenome segments and 23 were non-specific for a T. timopheevii subgenome (Table 1, Figure 2, and Supplementary Data 2).
Introgressions From Triticum timopheevii in Wheat
A range of different sized segments from each T. timopheevii chromosome of the At and G subgenomes were found present in the ILs and are annotated in Figures 3, 4, respectively. In total, there were 89 unique segments from the At subgenome of T. timopheevii and a whole chromosome 4At (Figure 3). It was estimated that all At chromosomes, except chromosome 3At, were introgressed entirely into wheat through overlapping segments with approximately 93.4% of chromosome 3At also being transferred. Figure 4 shows the 74 unique segments from the G subgenome of T. timopheevii that were introgressed into the ILs. It was estimated that 92.4 and 98.3% of chromosomes 4G and 6G, respectively, were introgressed in overlapping segments in addition to the entirety of the rest of the G subgenome chromosomes. The most segments introgressed into wheat were from chromosome 5At (25) and chromosome 2G (20) of T. timopheevii. Together, approximately 98.9% of the T. timopheevii genome had been introgressed into bread wheat mainly via small and large overlapping segments. The full set of chromosome-specific KASP markers detected 183 introgressions from the At subgenome and 126 from the G subgenome of T. timopheevii in the 99 ILs developed in this work.
Figure 3. A graphical representation of the size and location of unique introgressions from the At genome of Triticum timopheevii (purple) in wheat A and D genome chromosomes (shades of blue) for each of the seven linkage groups. On the left of each ideogram is the label of each unique introgression and on the right is the number (in purple) of introgression lines (ILs) the segment is found in. A table on the bottom right shows the percentage of each At chromosome that has been introgressed into the wheat genome and the number of ILs that contains a homozygous segment from each of the At chromosomes. A scale at the bottom of each linkage group shows the size of the ideograms.
Figure 4. A graphical representation of the size and location of unique introgressions from the G genome of Triticum timopheevii (red) in wheat B and D genome chromosomes (shades of blue) for each of the seven linkage groups. On the left of each ideogram is the label of each unique introgression and on the right is the number (in red) of introgression lines (ILs) the segment is found in. A table on the bottom right shows the percentage of each G chromosome that has been introgressed into the wheat genome and the number of ILs that contains a homozygous segment from each of the G chromosomes. A scale at the bottom of each linkage group shows the size of the ideograms.
Table 2 details which of these segments are carried by each of the ILs Tim1-99 and a diagnostic KASP marker for each segment is provided in Supplementary Data 3. Genotyping of ILs Tim1 to 35 with Set 4 KASP markers identified additional homozygous T. timopheevii segments that had not been previously identified by markers in Sets 1–3. It was also found that ILs Tim14 and Tim29 had the same genotype according to the markers used. Each IL contained between one and eight homozygous introgressions from T. timopheevii with an average of three segments per IL (Table 2).
Table 2. Number and description of Triticum timopheevii introgressions present in each wheat-T. timopheevii introgression line (IL), corresponding to the segments depicted in Figures 3, 4).
Since the introgressions are homozygous, the absence of the wheat allele allows the chromosome-specific KASP markers to identify the wheat chromosomes involved in each recombination event. This showed that the majority of the At subgenome segments were introgressed into the A genome chromosomes of wheat (At.A; Figure 3) and the G subgenome segments into the B genome of wheat (G.B; Figure 4). However, only two At segments were found to have recombined with the D genome of wheat (At.D; Figure 3) compared to 19 segments from the G subgenome (G.D; Figure 4). Multi-colour GISH (mcGISH) was used to validate these D genome introgressions.
Figure 5A shows mcGISH images of introgressions 4At.D1 and 2At.D1 where small At subgenome segments from T. timopheevii chromosomes 4At and 2At (shown in green) had recombined with chromosomes 4D and 2D of wheat (shown in red), respectively. Ideograms of the recombinant chromosomes representing the KASP marker genotyping (of chromosomes 4D and 2D) indicating the introgressed segments in orange (Figure 5B) matched the corresponding recombination breakpoints as indicated with white arrows in Figure 5A.
Figure 5. Multi-colour in situ hybridisation analysis of root metaphase spreads of chromosomes from wheat- Triticum timopheevii introgression lines (ILs). (A,C) Multi-colour genomic in situ hybridisation (mcGISH) of ILs with At.D−G.B introgressions, respectively. Green probe detects the At subgenome and blue detects the G subgenome of T. timopheevii. Red detects the D genome of wheat. (E–G) Multi-colour fluorescence in situ hybridisation (mcFISH) of (E) wheat cv. Chinese Spring (CS), (F) T. timopheevii and (G) ILs where the green and red signals show pSc119.2 and pAs.1 binding sites, respectively. Images for panels (E,F) were taken from Devi et al. (2019). Panels (B,D,H) a graphical representation of the introgressions present in these lines obtained through chromosome-specific KASP marker analysis depicting ideograms of the corresponding ILs shown either above with mcGISH (A,C) or on the left with mcFISH (E–G). Size of introgressions from At and G subgenomes is indicated in orange in a wheat chromosome shown in black. White arrows indicate the corresponding recombination breakpoints in the mcGISH and mcFISH images.
Of the 19 G.D introgressions, mcGISH images of 6 introgressions, 3G.D1-D3 and 3G.D5-D7, each containing a different sized segments from chromosome 3G of T. timopheevii (shown in blue) recombined with chromosome 3D of wheat (shown in red) are shown in Figure 5C. Ideograms of the recombinant chromosomes representing the KASP marker genotyping (of chromosome 3D) indicating each of the introgressed segments in orange are shown below in Figure 5D which correspond to the recombination breakpoints of the introgressions as indicated by white arrows in Figure 5C.
It is not possible to detect At.A and G.B introgressions with mcGISH and thus, mc-FISH was used to validate the KASP marker genotyping by analysing three introgressions, 1At.A4, 2G.B3, and 3G.B3, which were large enough to span a diagnostic FISH signal for T. timopheevii. Images in Figures 5E–G show the FISH karyotypes for the indicated chromosomes in wheat cv. ‘Chinese Spring’, T. timopheevii (accession P95-99.1-1) and the wheat-T. timopheevii ILs, respectively. Figure 5H represents the KASP genotyping of those introgressions in each IL. The gain of the T. timopheevii signal (as indicated by the white arrow) in Tim14 validates the introgression of chromosome 1At in chromosome 1A of wheat which otherwise lacks this signal. The loss of the wheat signal on chromosome 2BS (as indicated by the white arrow) in Tim73 indicates the presence of the 2G segment from T. timopheevii which otherwise lacks this signal on chromosome 2GS. In Tim73, the gain of the T. timopheevii signal from chromosome 3G validates the 3G.D3 introgression in chromosome 3D which otherwise lacks this signal on 3DL.
Discussion
The main aims of this work were (i) the development of a set of chromosome-specific KASP markers, well spread across the A and B genomes of wheat, for the genotyping of the homozygous lines and (ii) to generate homozygous ILs from the advanced backcrossed lines generated by Devi et al. (2019), thus enabling the multiplication and distribution of the lines for trait phenotyping.
The selection of plants with segments from T. timopheevii of various sizes and from different linkage groups was enabled through the use of chromosome-specific KASP markers. However, the markers needed to be of a sufficient density across all linkage groups in order to detect all introgressions present (including very small segments) and to accurately assess the size of a segment. The detection of small segments is particularly important in wheat-T. timopheevii ILs as the two subgenomes of T. timopheevii share considerable similarities to two of the genomes of wheat having originated from the same progenitors, T. urartu (A and At genomes) (Dvořák et al., 1993) and an Ae. speltoides-like species (B and G genomes) (Dvorak and Zhang, 1990). The more closely related two genomes were, the higher the level of recombination is possible between the two.
Development of Chromosome-Specific Kompetitive Allele Specific PCR Markers
Prior to this work other studies have developed molecular marker resources for the detection of T. timopheevii segments in a wheat background (Ji et al., 2007; Uhrin et al., 2012; Timonova et al., 2013; Devi et al., 2019) but none have reported the density and distribution across the wheat genome that has been achieved in this work. The aim in this study was to have a KASP markers at every 60–70 Mbp across the A and B genomes of wheat in order to carry out accurate genotyping of the ILs. Sets 1, 2, and 3 were designed either for multi-species use (Set 1) or were designed with the target of filling gaps between wheat and other wild relatives using SNPs discovered in those species (Set 2 – T. urartu and Set 3 – Am. muticum). As can be seen in Figure 2, most of the chromosomes still had significant gaps between the markers of these three sets. The KASP markers of Set 4, however, were designed to target the gaps based on SNPs between wheat and T. timopheevii. The strategy of custom filtering to remove SNPs with significant off-target hits and only retain those suitable for KASP assay design resulted in 72% of the chromosome-specific KASP assays designed successfully able to detect the presence of T. timopheevii segments introgressed into wheat chromosomes. It is difficult to compare the success of this strategy with some of those used previously as the early panels of KASP assays were mainly designed to be polymorphic between wheat and more than one wild relative. However, the same strategy was recently followed to increase the number of SNPs obtained that were polymorphic between wheat and Am. muticum (Grewal et al., 2022). In that work only 48% of the KASP assays designed functioned to detect Am. muticum DNA in a wheat background. The difference in the success rates might have been due to Am. muticum being an outbreeder and thus highly heterozygous while T. timopheevii is an inbreeder with a much lower level of heterozygosity resulting in robust KASP assay primers with fewer SNPs in the flanking regions of the SNP.
Combined, the 4 sets of markers now give excellent genome coverage with only three relatively significant gaps still needing to be filled (1 each on 4A, 5B, and 7A) with the largest of the gaps on 7A (118.7 Mbp; Table 1). The new, higher density of markers identified new homozygous segments in ILs Tim1 to 35 which had previously been genotyped with KASP marker sets 1, 2, and 3, highlighting the importance of increased marker density. The extra KASP markers will also enable these newly identified introgressions to be tracked in breeding programmes. A high density of markers at the end of chromosome 3B short arm can also be seen in Figure 2. This region was specifically targetted as recent phenotyping has shown this region of chromosome 3G to be associated with resistance to Fusarium Head Blight. The higher density of SNPs in this region demonstrates that the strategy of SNP discovery and KASP development employed will be successful for filling the remaining gaps. It is possible that even with the new, higher density of KASP markers, small introgressions between two KASP markers might not have been detected.
Homozygous Triticum timopheevii Introgressions in the Wheat Background
Homozygous introgressions from all linkage groups of T. timopheevii were obtained (Figures 3, 4) covering 98.9% of the T. timopheevii genome. These introgressions covered the entire linkage groups of most At and G subgenomes except chromosomes 3At, 4G, and 6G. However, more than 92% of these three chromosomes was introgressed into wheat. In total, 163 unique introgressions were observed, 89 (55%) from the At subgenome and 74 (45%) from the G subgenome. This was expected as previous reports have shown At and A genome chromosomes pair more readily than G and B genome chromosomes (Feldman, 1966).
The highest number of unique homozygous introgressions were recovered from 5At (25) to 2G (20) which were present in 47 and 33 ILs, respectively, more than any other chromosome from the respective subgenomes. The high number of introgressions was expected for 2G as this chromosome shows a higher level of transmission through the gametes than expected (Nyquist, 1962; Timonova et al., 2012; Devi et al., 2019). Preferential transmission of chromosome 2S of Ae. speltoides, proposed to be more closely related to the T. timopheevii G genome than to the B genome of wheat (Rodriguez et al., 2000a; Li et al., 2022), is well documented (Tsujimoto and Tsunewaki, 1988; King et al., 2018) although the gene causing the preferential transmission in T. timopheevii appears to be weaker than the gene for preferential transmission in Ae. speltoides. A high incidence of introgressions from T. timopheevii linkage groups 1At, 2At, 2G, 5At, 5G, and 6G was also found in lines screened by Leonova et al. (2011) which when taken forward by Timonova et al. (2013) also showed a high transmission of the 2G segment. Analysis of the recent wheat pangenome also showed the presence of a large 2G introgression in chromosome 2B of wheat cultivar LongReach Lancer conferring stem rust resistance (Walkowiak et al., 2020).
The smallest number of introgressions was obtained for 4At (2 introgressions and a whole chromosome substitution). Recombination between chromosomes 4A of wheat and 4At of T. timopheevii is likely to have been reduced due to translocation differences between the two chromosomes. Although both carry the 4AL/5AL translocation, inherited from the A genome progenitor T. urartu, 4A of wheat has undergone additional paracentric and pericentric inversions (Devos et al., 1995; Dvorak et al., 2018) which are not present in T. timopheevii (Maestra and Naranjo, 1999). These structural differences would result in reduced levels of recombination between the two chromosomes resulting in fewer introgressions.
A reduction in the number of introgressions was also expected from T. timopheevii chromosomal regions of 3At, 6At, 1G, and 4G as four species-specific translocations involving these chromosomes occur within the T. timopheevii groups; 6AtS/1GS, 1GS/4GS, 4GS/4AtL, and 4AtL/3AtL, plus a paracentric inversion in 6AtL (Jiang and Gill, 1994; Maestra and Naranjo, 1999; Rodriguez et al., 2000b; Salina et al., 2006). Chromosomes 3A and 3At have also been reported to have diverged through mutations in their primary DNA structure (Dobrovol’skaia et al., 2009). From the regions involved in species-specific translocations two introgressions each were obtained for 3AL (3At.A9 and 3At.A10), 1GS (1G.B1 and 1G.B2), and 6AL (6At.A3 and 6At.A4) and one each for 6AtS (6At.A2) and 4GS (4G.B1; Figures 3, 4).
A number of authors have also reported different accessions of T. timopheevii as carrying karyotype rearrangements (Kawahara and Tanaka, 1977; Badaeva et al., 1990, 1994, 2022; Taketa and Kawahara, 1996). Thus, in addition to the species-specific translocations described above, the accession used in this work may have carried additional rearrangements affecting recombination and hence introgression. This could also explain the differences in levels of recombination not just between chromosomes, e.g., the small number of introgressions from 7G, but between chromosome arms as seen in Figure 3. For example, recombination has occurred along the whole length of chromosomes 2At/2A, 7At/7A, 2B/2G, and 3B/3D/3G (although lower in the pericentric regions), while recombination between 1At/1A, 5At/5A, and 5G/5B has occurred more frequently in the long arm.
Due to high levels of synteny between At and A, and G and B genomes (Dvořák et al., 1993; Tsunewaki et al., 1996), the KASP markers validated for A and B genome chromosomes of wheat also, in majority of the cases, detected At and G chromosome segments, respectively (Table 1 and Figure 2). The closeness of the two groups of genomes also resulted in the majority of At and G subgenome introgressions recombining with the A and B genomes of wheat, respectively. Although, 2% (2 out of 98) of the introgressions from the At subgenome and 26% (19 out of 74) of the introgressions obtained for the G subgenome of T. timopheevii had recombined with the D genome of wheat (Figures 3, 4). It is possible that synteny between regions of the G subgenome and regions of the wheat D genome may be greater than between the G and B genomes especially for those chromosomes where either species-specific translocations or accession-specific translocations have occurred in T. timopheevii.
The specific wheat chromosomes involved in these recombinants were identified through genotyping with the chromosome-specific KASP markers that lie within the introgressed regions. The At.D and G.D introgressions were validated through mcGISH (Figures 5A–D) since the D genome chromosomes can be easily distinguished against the At and G subgenomes of T. timopheevii unlike the A and B genomes of wheat. Due to the high sequence similarity between the At and A and G and B subgenomes, the mcGISH probes cannot distinguish between these groups of chromosomes. This meant that the At.A and G.B introgressions produced in this study could not be visualised using mcGISH. Instead, a few of these introgressions were validated using mcFISH (Figures 5E–H). Previous reports have done extensive karyotyping of all T. timopheevii chromosomes using mcFISH (Badaeva et al., 2016; Devi et al., 2019), however, characterising interspecific introgression using mcFISH is very laborious and time-consuming. This highlights the utility of the chromosome-specific KASP markers developed in this work which are capable of detecting the presence and site of T. timopheevii introgressions in wheat.
In addition to the 99 homozygous wheat-T. timopheevii ILs generated in this work (Figure 1 and Table 2), further ILs have been generated within the Wheat Research Centre at Nottingham. With the availability of a high density of KASP markers, these ILs will be screened for unique introgressions from T. timopheevii with particular emphasis on those linkage groups from which relatively few introgressions have so far been obtained. Although the accession of T. timopheevii used in this work was not originally selected based on phenotyping data, the homozygous lines described are being screened for a variety of traits including Fusarium head blight resistance (Steed et al., 2022). Seed of lines Tim1-15, Tim17-19, and Tim21-35 have been deposited at the Germplasm Resource Unit at the John Innes Centre,1 while lines Tim36-99 will be distributed after seed multiplication.
Data Availability Statement
Raw reads data for T. timopheevii has been made available through the Grassroots data repository hosted by the Earlham Institute and funded by DFW programme (https://opendata.earlham.ac.uk/wheat). The genotyping data is available from the corresponding authors on reasonable request.
Author Contributions
C-yY, SA, DS, S-HE, SG, IK, and JK: germplasm generation. C-yY, DS, and SA: DNA extraction. SG, SA, DS, MO, and JW: KASP marker development and genotyping. C-yY, MO, and NW: in situ hybridisation analysis. BC and AH: sequencing and SNP discovery. SG: data analysis and curation. IK, JK, and SG: conceptualisation and manuscript writing. IK and JK: funding acquisition. All authors have read and agreed to the published version of the manuscript.
Funding
This work was supported by the Biotechnology and Biological Sciences Research Council (grant numbers BB/J004596/1 and BB/P016855/1) as part of the Wheat Improvement Strategic Programme (WISP) and Developing Future Wheat (DFW) programme. MO and JW were supported by the Biotechnology and Biological Sciences Research Council (grant number BB/S012796/1) as part of the International Wheat Yield Partnership (IWYP) programme. BC was supported by the BBSRC funded Norwich Research Park Biosciences Doctoral Training Partnership grant BB/M011216/1 and NW was supported by BBSRC funded Nottingham Doctoral Training Partnership grant BB/T008369/1.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.919519/full#supplementary-material
Supplementary Table 1 | Lists of all the chromosome-specific SNPs, between wheat and Triticum timopheevii, discovered in this work.
Supplementary Table 2 | Full set of validated chromosome-specific KASP markers that are detect subgenomes of Triticum timopheevii in a wheat background.
Supplementary Table 3 | KASP markers diagnostic for every introgression in each of the introgression lines.
Footnotes
References
Allard, R., and Shands, R. (1954). Inheritance of resistance to stem rust and powdery mildew in cytologically stable spring wheats derived from Triticum timopheevii. Phytopathology 44, 266–274.
Anand, L., and Rodriguez Lopez, C. M. (2022). ChromoMap: an R package for interactive visualization of multi-omics data and annotation of chromosomes. BMC Bioinform. 23:33. doi: 10.1186/s12859-021-04556-z
Badaeva, E. D., Konovalov, F. A., Knüpffer, H., Fricano, A., Ruban, A. S., Kehel, Z., et al. (2022). Genetic diversity, distribution and domestication history of the neglected GGAtAt genepool of wheat. Theor. Appl. Genet. 135, 755–776. doi: 10.1007/s00122-021-03912-0
Badaeva, E. D., Ruban, A. S., Zoshchuk, S. A., Surzhikov, S. A., Knüpffer, H., and Kilian, B. (2016). Molecular cytogenetic characterization of Triticum timopheevii chromosomes provides new insight on genome evolution of T. zhukovskyi. Plant Syst. Evol. 302, 943–956. doi: 10.1007/s00606-016-1309-3
Badaeva, E., Badaev, N., Gill, B., and Filatenko, A. (1994). Intraspecific karyotype divergence in Triticum araraticum (Poaceae). Plant Syst. Evol. 192, 117–145. doi: 10.1007/BF00985912
Badaeva, E., Boguslavsky, R., Badaev, N., and Zelenin, A. (1990). Intraspecific chromosomal polymorphism of Triticum araraticum (Poaceae) detected by C-banding technique. Plant Syst. Evol. 169, 13–24.
Badridze, G., Weidner, A., Asch, F., and Börner, A. (2009). Variation in salt tolerance within a Georgian wheat germplasm collection. Genet. Resour. Crop Evol. 56, 1125–1130. doi: 10.1007/s10722-009-9436-0
Brown-Guedira, G. L., Singh, S., and Fritz, A. K. (2003). Performance and mapping of leaf rust resistance transferred to wheat from Triticum timopheevii subsp. armeniacum. Phytopathology 93, 784–789. doi: 10.1094/phyto.2003.93.7.784
Brown-Guedira, G., Gill, B., Bockus, W., Cox, T., Hatchett, J., Leath, S., et al. (1996). Evaluation of a collection of wild timopheevi wheat for resistance to disease and arthropod pests. Plant Dis. 80, 928–933. doi: 10.1094/PD-80-0928
Camacho, C., Coulouris, G., Avagyan, V., Ma, N., Papadopoulos, J., Bealer, K., et al. (2009). BLAST+: architecture and applications. BMC Bioinform. 10:421. doi: 10.1186/1471-2105-10-421
Cruz, C. D., Peterson, G. L., Bockus, W. W., Kankanala, P., Dubcovsky, J., Jordan, K. W., et al. (2016). The 2NS translocation from Aegilops ventricosa confers resistance to the Triticum Pathotype of Magnaporthe oryzae. Crop Sci. 56, 990–1000. doi: 10.2135/cropsci2015.07.0410
DePristo, M. A., Banks, E., Poplin, R., Garimella, K. V., Maguire, J. R., Hartl, C., et al. (2011). A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat. Genet. 43, 491–498. doi: 10.1038/ng.806
Devi, U., Grewal, S., Yang, C. Y., Hubbart-Edwards, S., Scholefield, D., Ashling, S., et al. (2019). Development and characterisation of interspecific hybrid lines with genome-wide introgressions from Triticum timopheevii in a hexaploid wheat background. BMC Plant Biol. 19:183. doi: 10.1186/s12870-019-1785-z
Devos, K. M., Dubcovsky, J., Dvorak, J., Chinoy, C. N., and Gale, M. D. (1995). Structural evolution of wheat chromosomes 4A, 5A, and 7B and its impact on recombination. Theor. Appl. Genet. 91, 282–288. doi: 10.1007/bf00220890
Dobrovol’skaia, O. B., Sourdille, P., Bernard, M., and Salina, E. (2009). Chromosome synteny of the A genome of two evolutionary wheat lines. Genetika 45, 1548–1555.
Dvorak, J., and Zhang, H.-B. (1990). Variation in repeated nucleotide sequences sheds light on the phylogeny of the wheat B and G genomes. Proc. Natl. Acad. Sci. U.S.A. 87, 9640–9644. doi: 10.1073/pnas.87.24.9640
Dvořák, J., Terlizzi, P. D., Zhang, H.-B., and Resta, P. (1993). The evolution of polyploid wheats: identification of the A genome donor species. Genome 36, 21–31. doi: 10.1139/g93-004
Dvorak, J., Wang, L., Zhu, T., Jorgensen, C. M., Luo, M. C., Deal, K. R., et al. (2018). Reassessment of the evolution of wheat chromosomes 4A, 5A, and 7B. Theor. Appl. Genet. 131, 2451–2462. doi: 10.1007/s00122-018-3165-8
Dyck, P. (1992). Transfer of a gene for stem rust resistance from Triticum araraticum to hexaploid wheat. Genome 35, 788–792. doi: 10.1139/g92-120
Feldman, M. (1966). Identification of unpaired chromosomes in F1 hybrids involving Triticum aestivum and T. timopheevii. Can. J. Genet. Cytol. 8, 144–151.
Grewal, S., Coombes, B., Joynson, R., Hall, A., Fellers, J., Yang, C.-Y., et al. (2022). Chromosome-specific KASP markers for detecting Amblyopyrum muticum segments in wheat introgression lines. Plant Genome 15:e20193. doi: 10.1002/tpg2.20193
Grewal, S., Gardiner, L. J., Ndreca, B., Knight, E., Moore, G., King, I. P., et al. (2017). Comparative mapping and targeted-capture sequencing of the gametocidal loci in Aegilops sharonensis. Plant Genome 10, lantgenome2016.2009.0090. doi: 10.3835/plantgenome2016.09.0090
Grewal, S., Guwela, V., Newell, C., Yang, C.-Y., Ashling, S., Scholefield, D., et al. (2021). Generation of doubled haploid wheat-Triticum urartu introgression lines and their characterisation using chromosome-specific KASP markers. Front. Plant Sci. 12:643636. doi: 10.3389/fpls.2021.643636
Grewal, S., Hubbart-Edwards, S., Yang, C., Devi, U., Baker, L., Heath, J., et al. (2020a). Rapid identification of homozygosity and site of wild relative introgressions in wheat through chromosome-specific KASP genotyping assays. Plant Biotechnol. J. 18, 743–755. doi: 10.1111/pbi.13241
Grewal, S., Othmeni, M., Walker, J., Hubbart-Edwards, S., Yang, C.-Y., Scholefield, D., et al. (2020b). Development of Wheat-Aegilops caudata introgression lines and their characterisation using genome-specific KASP markers. Front. Plant Sci. 11:606. doi: 10.3389/fpls.2020.00606
Grewal, S., Yang, C., Edwards, S. H., Scholefield, D., Ashling, S., Burridge, A. J., et al. (2018). Characterisation of Thinopyrum bessarabicum chromosomes through genome-wide introgressions into wheat. Theor. Appl. Genet. 131, 389–406. doi: 10.1007/s00122-017-3009-y
Hao, Z., Lv, D., Ge, Y., Shi, J., Weijers, D., Yu, G., et al. (2020). RIdeogram: drawing SVG graphics to visualize and map genome-wide data on the idiograms. PeerJ Comput. Sci. 6:e251. doi: 10.7717/peerj-cs.251
Hu, X., Liu, J., Zhang, L., Wu, B., Hu, J., Liu, D., et al. (2017). Zn and Fe concentration variations of grain and flag leaf and the relationship with NAM-G1 gene in Triticum timopheevii (Zhuk.) Zhuk. ssp. timopheevii. Cereal Res. Commun. 45, 421–431. doi: 10.1556/0806.45.2017.022
International Wheat Genome Sequencing Consortium [IWGSC],, R. Eversole, K., Feuillet, C., Keller, B., Rogers, J., et al. (2018). Shifting the limits in wheat research and breeding using a fully annotated reference genome. Science 361:eaar7191. doi: 10.1126/science.aar7191
Järve, K., Peusha, H. O., Tsymbalova, J., Tamm, S., Devos, K. M., and Enno, T. M. (2000). Chromosomal location of a Triticum timopheevii–derived powdery mildew resistance gene transferred to common wheat. Genome 43, 377–381. doi: 10.1139/g99-141
Ji, J. H., Cao, A. Z., Wang, H. Y., Qin, B., Wang, S. L., Kong, F., et al. (2007). Discrimination of the Triticum aestivum-T. timopheevii introgression lines using PCR-based molecular markers. Yi Chuan 29, 1256–1262. doi: 10.1360/yc-007-1256
Jiang, J., and Gill, B. S. (1994). Different species-specific chromosome translocations in Triticum timopheevii and T. turgidum support the diphyletic origin of polyploid wheats. Chromosome Res. 2, 59–64. doi: 10.1007/BF01539455
Jørgensen, J. H., and Jensen, C. (1972). Genes for resistance to wheat powdery mildew in derivatives of Triticum timopheevi and T. carthlicum. Euphytica 21, 121–128. doi: 10.1007/BF00040557
Kawahara, T., and Tanaka, M. (1977). Six chromosome types in Triticum araraticum Jakubz. differing with reciprocal translocations. Jpn. J. Genet. 52, 261–267. doi: 10.1266/jjg.52.261
King, I. P., Forster, B. P., Law, C. C., Cant, K. A., Orford, S. E., Gorham, J., et al. (1997). Introgression of salt-tolerance genes from Thinopyrum bessarabicum into wheat. New Phytol. 137, 75–81. doi: 10.1046/j.1469-8137.1997.00828.x
King, J., Grewal, S., Yang, C. Y., Hubbart, S., Scholefield, D., Ashling, S., et al. (2017). A step change in the transfer of interspecific variation into wheat from Ambylopyrum muticum. Plant Biotechnol. J. 15, 217–226. doi: 10.1111/pbi.12606
King, J., Grewal, S., Yang, C.-Y., Hubbart Edwards, S., Scholefield, D., Ashling, S., et al. (2018). Introgression of Aegilops speltoides segments in Triticum aestivum and the effect of the gametocidal genes. Ann. Bot. 121, 229–240. doi: 10.1093/aob/mcx149
Lehmensiek, A., Bovill, W., Banks, P., and Sutherland, M. (2008). “Molecular characterization of a Triticum timopheevii introgression in a Wentworth/Lang population,” in Proceedings of the 11th International Wheat Genetics Symposium, Vol. 3, eds R. Appels, R. Eastwood, E. Lagudah, P. Langridge, M. Mackay, L. McIntyre, et al. (Brisbane, QLD: Sydney University Press), 768.
Lehmensiek, A., Campbell, A., Sutherland, M., Williamson, P., Michalowitz, M., and Daggard, G. (2004). QTLs for black-point resistance in wheat and the identification of potential markers for use in breeding programmes. Plant Breed. 123, 410–416. doi: 10.1111/j.1439-0523.2004.01013.x
Leonova, I. N., Budashkina, E. B., Flath, K., Weidner, A., Börner, A., and Röder, M. S. (2010). Microsatellite mapping of a leaf rust resistance gene transferred to common wheat from Triticum timopheevii. Cereal Res. Commun. 38, 211–219. doi: 10.1556/CRC.38.2010.2.7
Leonova, I. N., Budashkina, E. B., Kalinina, N. P., Röder, M. S., Börner, A., and Salina, E. A. (2011). Triticum aestivum× Triticum timopheevii introgression lines as a source of pathogen resistance genes. Czech J. Genet. Plant Breed. 47:S49. doi: 10.17221/3254-CJGPB
Li, G.-R., Liu, C., Li, C.-H., Zhao, J.-M., Zhou, L., Dai, G., et al. (2013). Introgression of a novel Thinopyrum intermedium St-chromosome-specific HMW-GS gene into wheat. Mol. Breed. 31, 843–853. doi: 10.1007/s11032-013-9838-8
Li, H. (2011). A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics 27, 2987–2993. doi: 10.1093/bioinformatics/btr509
Li, H. (2013). Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv [Preprint]. doi: 10.48550/arXiv.1303.3997
Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., et al. (2009). The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078–2079. doi: 10.1093/bioinformatics/btp352
Li, L.-F., Zhang, Z.-B., Wang, Z.-H., Li, N., Sha, Y., Wang, X.-F., et al. (2022). Genome sequences of five Sitopsis species of Aegilops and the origin of polyploid wheat B subgenome. Mol. Plant 15, 488–503. doi: 10.1016/j.molp.2021.12.019
Maestra, B., and Naranjo, T. (1999). Structural chromosome differentiation between Triticum timopheevii and T. turgidum and T. aestivum. Theor. Appl. Genet. 98, 744–750. doi: 10.1007/s001220051130
Malihipour, A., Gilbert, J., Fedak, G., Brûlé-Babel, A., and Cao, W. (2016). Characterization of agronomic traits in a population of wheat derived from Triticum timopheevii and their association with Fusarium head blight. Eur. J. Plant Pathol. 144, 31–43. doi: 10.1007/s10658-015-0744-2
Malihipour, A., Gilbert, J., Fedak, G., Brûlé-Babel, A., and Cao, W. (2017). Mapping the A genome for QTL conditioning resistance to Fusarium head blight in a wheat population with Triticum timopheevii background. Plant Dis. 101, 11–19. doi: 10.1094/PDIS-02-16-0144-RE
McIntosh, R. (1983). “Genetic and cytogenetic studies involving Lr18 for resistance to Puccinia recondita,” in Proceedings of the Sixth International Wheat Genetics Symposium/Edited by Sadao Sakamoto, (Kyoto: Kyoto University).
McIntosh, R., and Gyarfas, J. (1971). Triticum timopheevii as a source of resistance to wheat stem rust. Z. Pflanzenzucht. 66, 240–248.
McIntyre, C. L., Pereira, S., Moran, L. B., and Appels, R. (1990). New Secale cereale (rye) DNA derivatives for the detection of rye chromosome segments in wheat. Genome 33, 635–640. doi: 10.1139/g90-094
Niranjana, M., Vinod, Sharma, J. B., Mallick, N., Tomar, S. M. S., and Jha, S. K. (2017). Cytogenetic analysis and mapping of leaf rust resistance in Aegilops speltoides Tausch derived bread wheat line Selection2427 carrying putative gametocidal gene(s). Genome 60, 1076–1085. doi: 10.1139/gen-2017-0107
Nyquist, W. (1962). Differential fertilization in the inheritance of stem rust resistance in hybrids involving a common wheat strain derived from Triticum timopheevii. Genetics 47:1109. doi: 10.1093/genetics/47.8.1109
Perugini, L. D., Murphy, J. P., Marshall, D., and Brown-Guedira, G. (2008). Pm37, a new broadly effective powdery mildew resistance gene from Triticum timopheevii. Theor. Appl. Genet. 116, 417–425. doi: 10.1007/s00122-007-0679-x
Peusha, K. O., Stephan, U., Hsam, S. L. K., Felsenstein, F. G., Enno, T. M., and Zeller, F. J. (1995). Identification of genes for resistance to powdery mildew in common wheat (Triticum aestivum L.). IV: breeding lines derived from wide crosses of russian cultivars with species T. timopheevii Zhuk., T. militinae Zhuk. et Migus H., T. dicoccum (Schrank) Schuebl., and Aegilops speltoides Tausch. Russian J. Genet. 31, 181–187.
Qin, B., Cao, A., Wang, H., Chen, T., You, F. M., Liu, Y., et al. (2011). Collinearity-based marker mining for the fine mapping of Pm6, a powdery mildew resistance gene in wheat. Theor. Appl. Genet. 123, 207–218. doi: 10.1007/s00122-011-1577-9
Ramirez-Gonzalez, R. H., Uauy, C., and Caccamo, M. (2015). PolyMarker: a fast polyploid primer design pipeline. Bioinformatics 31, 2038–2039. doi: 10.1093/bioinformatics/btv069
Rayburn, A. L., and Gill, B. S. (1986). Molecular identification of the D-genome chromosomes of wheat. J. Hered. 77, 253–255. doi: 10.1093/oxfordjournals.jhered.a110231
Reynolds, M. P., Calderini, D. F., Condon, A. G., and Rajaram, S. (2001). Physiological basis of yield gains in wheat associated with the LR19 translocation from Agropyron elongatum. Euphytica 119, 139–144. doi: 10.1023/A:1017521800795
Rigby, P. W., Dieckmann, M., Rhodes, C., and Berg, P. (1977). Labeling deoxyribonucleic acid to high specific activity in vitro by nick translation with DNA Polymerase I. J. Mol. Biol. 113, 237–251. doi: 10.1016/0022-2836(77)90052-3
Rodriguez, S., Maestra, B., Perera, E., Diez, M., and Naranjo, T. (2000a). Pairing affinities of the B-and G-genome chromosomes of polyploid wheats with those of Aegilops speltoides. Genome 43, 814–819. doi: 10.1139/g00-055
Rodriguez, S., Perera, E., Maestra, B., Díez, M., and Naranjo, T. (2000b). Chromosome structure of Triticum timopheevii relative to T. turgidum. Genome 43, 923–930.
Salina, E. A., Leonova, I. N., Efremova, T. T., and Röder, M. S. (2006). Wheat genome structure: translocations during the course of polyploidization. Funct. Integr. Genomics 6, 71–80. doi: 10.1007/s10142-005-0001-4
Sharma, P., Sheikh, I., Kumar, S., Verma, S. K., Kumar, R., Vyas, P., et al. (2018). Precise transfers of genes for high grain iron and zinc from wheat-Aegilops substitution lines into wheat through pollen irradiation. Mol. Breed. 38:81. doi: 10.1007/s11032-018-0836-8
Singh, A. K., Sharma, J. B., Vinod, Singh, P. K., Singh, A., and Mallick, N. (2017). Genetics and mapping of a new leaf rust resistance gene in Triticum aestivum L. × Triticum timopheevii Zhuk. derivative ‘Selection G12’. J. Genet. 96, 291–297. doi: 10.1007/s12041-017-0760-4
Steed, A., King, J., Grewal, S., Yang, C.-Y., Clarke, M., King, I., et al. (2022). Identification of Fusarium head blight resistance in Triticum timopheevii accessions and characterisation of wheat-T. timopheevii introgression lines for enhanced resistance. bioRxiv [Preprint]. doi: 10.1101/2022.05.13.491811
Taketa, S., and Kawahara, T. (1996). C-banding analysis on wild Emmer (Triticum dicoccoides Körn) strains with and without spontaneous reciprocal translocations. Theor. Appl. Genet. 92, 173–178. doi: 10.1007/BF00223373
Timonova, E. M., Leonova, I. N., Belan, I. A., Rosseeva, L. P., and Salina, E. A. (2012). The influence of particular chromosome regions of Triticum timopheevii on the formation of resistance to diseases and quantitative traits in common wheat. Russian J. Genet. Appl. Res. 2, 330–343. doi: 10.1134/S2079059712040119
Timonova, E. M., Leonova, I. N., Röder, M. S., and Salina, E. A. (2013). Marker-assisted development and characterization of a set of Triticum aestivum lines carrying different introgressions from the T. timopheevii genome. Mol. Breed. 31, 123–136. doi: 10.1007/s11032-012-9776-x
Tsujimoto, H., and Tsunewaki, K. (1988). Gametocidal genes in wheat and its relatives. III. Chromosome location and effects of two Aegilops speltoides-derived gametocidal genes in common wheat. Genome 30, 239–244. doi: 10.1139/g88-041
Tsunewaki, K., Wang, G.-Z., and Matsuoka, Y. (1996). Plasmon analysis of Triticum (wheat) and Aegilops. 1. Production of alloplasmic common wheats and their fertilities. Genes Genet. Syst. 71, 293–311. doi: 10.1266/ggs.71.293
Uhrin, A., Szakács, É, Láng, L., Bedő, Z., and Molnár-Láng, M. (2012). Molecular cytogenetic characterization and SSR marker analysis of a leaf rust resistant wheat line carrying a 6G(6B) substitution from Triticum timopheevii (Zhuk.). Euphytica 186, 45–55. doi: 10.1007/s10681-011-0483-1
Walkowiak, S., Gao, L., Monat, C., Haberer, G., Kassa, M. T., Brinton, J., et al. (2020). Multiple wheat genomes reveal global variation in modern breeding. Nature 588, 277–283. doi: 10.1038/s41586-020-2961-x
Keywords: wheat, wild relative introgression, Triticum timopheevii, KASP markers, GISH
Citation: King J, Grewal S, Othmeni M, Coombes B, Yang C-y, Walter N, Ashling S, Scholefield D, Walker J, Hubbart-Edwards S, Hall A and King IP (2022) Introgression of the Triticum timopheevii Genome Into Wheat Detected by Chromosome-Specific Kompetitive Allele Specific PCR Markers. Front. Plant Sci. 13:919519. doi: 10.3389/fpls.2022.919519
Received: 13 April 2022; Accepted: 12 May 2022;
Published: 01 June 2022.
Edited by:
Laura Toppino, Council for Research in Agriculture and the Analysis of Agrarian Economy (CREA), ItalyReviewed by:
Cheng Liu, Shandong Academy of Agricultural Sciences, ChinaElena A. Salina, Institute of Cytology and Genetics (RAS), Russia
Copyright © 2022 King, Grewal, Othmeni, Coombes, Yang, Walter, Ashling, Scholefield, Walker, Hubbart-Edwards, Hall and King. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Julie King, anVsaWUua2luZ0Bub3R0aW5naGFtLmFjLnVr; Surbhi Grewal, c3VyYmhpLmdyZXdhbEBub3R0aW5naGFtLmFjLnVr
†These authors have contributed equally to this work and share first authorship