Genetic association of OPR genes with resistance to Hessian fly in hexaploid wheat
© Tan et al.; licensee BioMed Central Ltd. 2013
Received: 16 March 2013
Accepted: 17 May 2013
Published: 1 June 2013
Skip to main content
© Tan et al.; licensee BioMed Central Ltd. 2013
Received: 16 March 2013
Accepted: 17 May 2013
Published: 1 June 2013
Hessian fly (Mayetiola destructor) is one of the most destructive pests of wheat. The genes encoding 12-oxo-phytodienoic acid reductase (OPR) and lipoxygenase (LOX) play critical roles in insect resistance pathways in higher plants, but little is known about genes controlling resistance to Hessian fly in wheat.
In this study, 154 F6:8 recombinant inbred lines (RILs) generated from a cross between two cultivars, ‘Jagger’ and ‘2174’ of hexaploid wheat (2n = 6 × =42; AABBDD), were used to map genes associated with resistance to Hessian fly. Two QTLs were identified. The first one was a major QTL on chromosome 1A (QHf.osu-1A), which explained 70% of the total phenotypic variation. The resistant allele at this locus in cultivar 2174 could be orthologous to one or more of the previously mapped resistance genes (H9, H10, H11, H16, and H17) in tetraploid wheat. The second QTL was a minor QTL on chromosome 2A (QHf.osu-2A), which accounted for 18% of the total phenotypic variation. The resistant allele at this locus in 2174 is collinear to an Yr17-containing-fragment translocated from chromosome 2N of Triticum ventricosum (2n = 4 × =28; DDNN) in Jagger. Genetic mapping results showed that two OPR genes, TaOPR1-A and TaOPR2-A, were tightly associated with QHf.osu-1A and QHf.osu-2A, respectively. Another OPR gene and three LOX genes were mapped but not associated with Hessian fly resistance in the segregating population.
This study has located two major QTLs/genes in bread wheat that can be directly used in wheat breeding programs and has also provided insights for the genetic association and disassociation of Hessian fly resistance with OPR and LOX genes in wheat.
Hessian fly [Hf, Mayetiola destructor (Say)] is one of the most destructive pests of hexaploid wheat (T. aestivum L., AABBDD, 2n = 6 × =42) in the United States and worldwide . Hessian fly larvae live between leaf-sheaths at seedling stage in fall, inhibit wheat growth irreversibly, and the infested plants lose vigor and die after larvae become pupae. Hessian fly larvae also attack stalks in spring, and the attacked plants may break easily before harvest. Deployment of natural and genetic resistance in locally adapted wheat cultivars is the most effective, economical, and environmentally safe method to control this devastating insect.
At least 33 Hessian fly resistance genes have been identified, designated as H1 to H32 and Hdic. Among these Hessian fly resistance genes, however, only 8 genes (H1-H5, H7, H8, and H12) were identified in hexaploid wheat [3–8]. The remaining 25 genes were identified in distant and close relatives of hexaploid wheat. There are 15 genes (H6, H9–H11, H14–H20, H28, H29, H31, and Hdic) identified in tetraploid wheat species Triticum turgidum subsp. durum (AABB, 2n = 4 × =28) [9–19], and 6 of them (H6, H9-H11, H31, and Hdic) have been introgressed from tetraploid to hexaploid wheat [7, 18–20]. Other genes originating in wild diploid wheat or other relatives and introgressed to hexaploid wheat include six genes (H13, H22, H23, H24, H26, and H32) from Aegilops tauschii (DD, 2n = 2 × =14) [21–25], two genes (H21 and H25) from Secale cerale L. (RR, 2n = 2 × =14) [26, 27], one gene (H27) from Aegilops ventricosa (DvDvMvMv, 2n = 4 × =28, ), and one gene (H30) from Aegilops triuncialis (CCUU, 2n = 4 × =28) . These genes provide important resistance sources but are problematic in variety development programs when they are associated with alien linkage drag.
Simple sequence repeat (SSR) or microsatellite markers were used to map Hessian fly resistance genes in previous studies. Seven genes (H5, H9, H10, H11, H16, H17, and Hdic) were mapped on chromosome 1A, and these genes may comprise a cluster (or family) of Hessian fly resistance genes in the distal gene-rich region of wheat chromosome 1AS [19, 20, 30, 31]. The remaining mapped genes include eight genes (H3, H6, H12, H14, H15, H19, H28, and H29) on chromosome 5A [10, 17, 32, 33], three genes (H24, H26, and H32) on chromosome 3D [23–25], three genes (H18, H20 and H21) on chromosome 2B [16, 26, 29], two genes (H13 and H23) on chromosome 6D [23, 34], H7 on chromosome 5D , H22 on chromosome 1D , H25 on chromosome 6B , H31 on chromosome 5B , and H27 on chromosome 4D . Current understanding of SSR markers is that they can be effectively applied to marker-assisted selection (MAS) due to their relative simplicity. However, the ‘repeat’ feature of an SSR marker results in multiple and inconsistent locations of the same marker among divergent wheat cultivars . Only a gene marker that is developed for the specific functional polymorphism of a gene responsible for a given trait can provide the ultimate resolution needed for selection of the trait .
The best strategy to find the regulatory sites of a gene is to clone the gene/QTL. The cloned gene is then used to develop perfect gene markers for use in conventional breeding programs or to manipulate transgenic wheat. To date, 14 genes have been cloned from wheat using the positional cloning strategy [38, 39], but no gene has been cloned for resistance to Hessian fly. To clone a gene from hexaploid wheat using the map-based cloning approach remains a challenge because of the complexity imparted by three homoeologous genomes, the large genome size (17 Gb), high content of repetitive sequences (>80%), and the low polymorphism rate [40–43]. In recent studies, genetic association is also use to identify functional genes for important traits in wheat .
Recent progress in the application of high-throughput sequencing technologies and development of sophisticated genomic mapping tools has accelerated identification of agriculturally important genes in wheat. Numerous wheat ESTs assigned to wheat deletion bins have facilitated the physical mapping of a gene [45, 46]. The availability of genomic sequences from the model species rice  and Brachypodium distachyon, and synteny conservation in certain genomic regions between these species and wheat [48, 49] have allowed more precise determination of the physical distance between two markers flanking the target gene by a comparative genomics analysis. Newly available wheat genomic sequences, though at low coverage (http://www.cerealsdb.uk.net/), are useful for designing primers specific to a homoeologous gene for genomic mapping. With the integration of these techniques and tools, mutants of several genes including VRN-A1, VRN-D3, PPD-D1, Lr34-D, Yr17, and Pm3 have been identified in cultivated bread wheat. These genes each were initially mapped under the peak of a QTL for a relevant trait and their mutants were eventually confirmed. In this study, we aimed to determine if any genes known to confer insect resistance in plants are associated with a QTL for resistance in our mapping populations or in a genomic region where a QTL has been reported in published mapping populations.
Plants possess multiple defense mechanisms in response to mechanical damage due to insect attack. Jasmonic acid (JA) and its conjugates, jasmonates, play a central role in regulating defense responses of plants to insect herbivores . In higher plants, JA is synthesized via the octadecanoid pathway consisting of several enzymatic steps. The early steps of this pathway occur in chloroplasts where linolenic acid is converted to 12-oxo-phytodienoic acid (OPDA) by means of three enzymes, lipoxygenase (LOX), allene oxide synthase (AOS), and allene oxide cylase [54–56]. OPDA is subsequently reduced in a cyclopentenone ring by a peroxisome-localized enzyme, 12-oxo-phytodienoic acid reductase 3 (OPR3). The reaction product then undergoes three cycles of oxidation in the peroxisome, generating JA [57–59].
In this study, we mapped 3 OPR genes and 3 LOX genes in a population of recombinant inbred lines (RILs) that was generated from a cross between two locally adapted winter cultivars, ‘Jagger’ and ‘2174’, and showed demonstrable segregation for resistance to Hessian fly. The implication in association and disassociation of the QTLs for resistance to Hessian fly with these OPR and LOX genes is discussed.
The second QTL was a minor QTL on chromosome 2A (QHf.osu-2A) (Figure 2B). The peak of this QTL was associated with a group of several tightly linked SSR markers and the marker for Yr17 that was translocated from chromosome 2N of Triticum ventricosum (2n = 4 × =28; DDNN) in Jagger [52, 60]. The LOD value for QHf.osu-2A was 3.5 and this QTL accounted for 18% of the total phenotypic variation. Phenotypic data showed that those RILs which carried the Jagger allele had resistance at the level of 23.1%, whereas those RILs that carried the 2174 allele had resistance at the level of 51.9%. The resistant allele at this locus in 2174 is collinear to an Yr17-containing-fragment translocated from chromosome 2N in Jagger. Jagger attained the translocated chromosome 2N fragment conferring resistance to stripe rust . If the translocated fragment in Jagger was introgressed into a breeding line or a cultivar such as 2174, the novel breeding line could lose approximately 60% of Hessian fly resistance conferred by the OPR2-A gene on chromosome 2A in 2174.
The second OPR gene, TaOPR2-A, which showed 79% identity (Additional file 1) to TaOPR1-A within a range of 584 bp (E = 2e-112), showed complete linkage with a group of markers for Yr17 translocated from chromosome 2N and under the peak of QHf.osu-2A. The wheat EST BE403717 that was mapped in 1AS-3 FL 0.86 deletion bin was used to blast against GenBank Nucleotide Collection (nr/nt) databases, and it showed 82% identity within a range of 318 bp (E = 1e-83) to an orthologous OPR gene in rice BAC from chromosome 6 (GenBank: AP004741). Specific primers were designed based on the grouped sequences retrieved from the wheat genomic DNA databases. The primers OPR22-C1-F3 and OPR-R2 amplified a 634 bp fragment that was digested with HincII into 20 bp, 321 bp and 292 bp for the Jagger allele (GenBank: KF035084) and into 342 bp, 215 bp, and 77 bp for the 2174 allele (GenBank: KF035085) (Figure 4B).
Three LOX genes were mapped in this study. Initially, a wheat EST (GenBank: BF482663) that was mapped in the 1AS-3 FL 0.86 deletion bin was used to blast against wheat gDNA databases (http://www.cerealsdb.uk.net) and GenBank EST databases for the homoeologous and paralogous LOX genes. This wheat EST was expected to reside in a region collinear with rice chromosome 5. It showed 86% identity within a range of only 43 bp (E = 7e-05) to an orthologous gene in rice BAC from chromosome 5 (GenBank: AC136525), but it was more likely orthologous to other LOX genes in rice, with 87% identity in a range of only 294 bp (E = 5e-96) to one in chromosome 2 (GenBank: AP004184), with 72% identity (E = 2e-38) to two genes on chromosome 8 (GenBank: AP005816). Homoeologous chromosomes of hexaploid wheat could have three genes for each orthologous gene in rice. Primers were designed specifically to each of the groups retrieved from the sequences by BF482663.
The first LOX gene was mapped using primers LOX-C5-F5 and LOX-C5-R6 specific to chromosome 1A (Figure 3C). The primers amplified a fragment of approximately 2,500 bp that was digested with ScrFI into several fragments, distinguishable by an extra fragment of approximately 250 bp for the 2174 allele (GenBank: KF035089) compared with the Jagger allele (GenBank: KF035088) (Figure 4C). This LOX gene was mapped on chromosome 1A (TaLOX1-A) but outside of the QHf.osu-1A region (Figure 2A).
The second LOX gene was mapped using primers LOX0-F4 and LOX0-R6 that amplified a common band in both Jagger and 2174 and an additional lower band in Jagger (Figure 4F). PCR products from Jagger were purified from the agarose gel, and direct sequence analysis showed that it was part of a LOX gene (Additional file 2). The high quality sequence result indicated that it was from a pseudo gene in Jagger (GenBank: KF035090). The dominant marker for the LOX gene in Jagger was mapped with a linkage group of SSR markers on chromosome 2B (Figure 5B), and this gene was thus designed TaLOX2-B.
TaOPR1-A was used to screen the T. durum BAC library, and one positive clone (73J24) was sequenced with low coverage to find genes that could be used to determine collinear regions of rice and Brachypodium. Four OPR genes (GenBank: KF035074, KF035077, KF035080, and KF035081), and three genes encoding stress-induced receptor-like kinase (Srk) (GenBank: KF035093), catalytic domain of protein kinases (Pkc) (GenBank: KF035099), and calcineurin-like phosphoesterase-like (Clp) (GenBank: KF035096) respectively, were present in the wheat BAC clone (Figure 6C). Allelic variation was observed in the Srk gene for the Jagger allele (GenBank: KF035094) and the 2174 allele (GenBank: KF035095), the Clp gene for the Jagger allele (GenBank: KF035097) and the 2174 allele (GenBank: KF035098), and the Pkc gene for the Jagger allele (GenBank: KF035100) and the 2174 allele (GenBank: KF035101). Howerver, no crossover was found between the OPR genes and these linked genes in the the QHf.osu-1A region.
TaOPR1-A displayed the highest identity to orthologous OPR genes in a rice BAC from chromosome 6 (GenBank: AP004741), which contains 7 OPR genes arranged in tandem (Figure 6A). Except for the wheat OPR genes, however, the other wheat genes (Srk, Clp, and Pkc) were not present in the flanking regions of the OPR genes in the two BAC clones shown in Figure 6A or in other BAC clones containing OPRs in rice. TaOPR1-A displayed the highest identity to three orthologous OPR genes in Brachypodium (GenBank: NC016131, Figure 6B), but again, except for the wheat OPR genes and Srk gene, the other wheat genes (Clp and Pkc) were not observed in the collinear region of Brachypodium. The orthologous genes in Brachypodium were not present in the flanking regions of the OPR genes in rice either.
A resistance gene against Hessian fly has been repeatedly mapped to the end of the short arm of chromosome 1A in previous studies, and it was suggested that this genomic region contained a cluster of major dominant resistance genes against multiple Hessian fly biotypes, including H5, H9, H10, H11, H16, H17, and Hdic[19, 20, 30, 31]. However, all of the previous studies were performed in the tetraploid wheat T. durum. Our study is the first report that a resistance gene against Hessian fly exists in the short arm of chromosome 1A in bread wheat. The resistance gene at the QHf.osu-1A locus observed in hexaploid wheat cv. 2174 could be orthologous to one or more of the previously mapped resistance genes in tetraploid wheat. This study provides an effective resistance source to manage biotype GP that frequently damage bread wheat cultivars in the southern Great Plains. Further studies need to test if the QHf.osu-1A gene is allelic to genes H16 and H17 that confer resistance against Hessian fly biotype L, the most virulent and prevalent biotype in the eastern USA . The resistance gene at QHf.osu-1A in 2174 and its derived cultivars can be immediately utilized to control Hessian fly biotype GP in winter wheat improvement programs in the southern Great Plains in USA. The gene may also be useful for pyramiding to prolong the effective periods of other resistance genes.
Mapping of QHf.osu-1A showed that it is closely linked with Pm3a. Previous studies indicated that Pm3 was linked to H9 at a genetic distance of 4.5 cM . Several alleles have been characterized for the Pm3 gene . Leaf rust resistance gene Lr10, which has three paralogous genes and is effective against Puccinia triticina Eriks, was also mapped in the same chromosomal region on 1AS [64, 65]. The presence of a resistance gene cluster in the gene-rich distal region of 1AS in bread wheat has made it difficult to determine which gene is the candidate for QHf.osu-1A. An effort was made to construct a fine physical map for QHf.osu-1A. However, low collinearity of the gene order in this region among wheat, rice, or Brachypodium indicates that the fine physical map for QHf.osu-1A cannot be established by using genome information from rice or Brachypodium only. The recent sequencing of wheat genomes may provide a powerful tool in cloning QHf.osu-1A.
Two OPR genes were identified in association with two QTLs for resistance against Hessian fly, but we cannot yet decisively conclude if either of the candidate OPR genes is responsible for QHf.osu-1A. The two QTLs were discovered in the same population, and resistance to Hessian fly is inseparable between the two QTLs. Two specific RIL lines (#23 and #26) have been selected to backcross with the 2174 parental line to generate progeny that will segregate for a single Hessian fly gene. In these backcross populations, QHf.osu-1A will occur in the heterozygous state, whereas QHf.osu-2A will be fixed in a homozygous genetic background. We expect that the two alleles of QHf.osu-1A will produce clear segregation for resistance to Hessian fly. BC1F2-3 populations will be used also to determine degree of dominance for the Hessian fly gene.
TaLOX6-B was mapped between Xwmc398 and Xbarc198 toward the centromere of wheat chromosome 6B. It was reported that H25 was translocated from rye to the long arm of chromosome 6B of wheat , and thus, TaLOX6-B is not allelic to H25. TaLOX2-B is located between Xwms388 and Xwms120 on the long arm of chromosome 2B. H20 was transferred from T. durum to chromosome 2B, and H21 was translocated from rye to chromosome 2B in wheat. However, more information is needed to further clarify any allelic relationship of TaLOX6-B with either H20 or H21. A recent study indicates that a novel wheat gene encoding a lectin-like protein (Hfr-3) was associated with response to Hessian fly . The homoeologous Hfr-3 genes are located on group 7 chromosomes of the hexaploid wheat. Yet, no Hessian fly resistance gene was reported on chromosome 7B where TaOPR7-B was mapped.
In summary, a frequent objective of wheat breeding is to pyramid multiple genes for resistance to diseases and insects into a single cultivar. Compared with 2174, Jagger attained a translocated fragment that contains Lr37, Yr17 and Sr38 conferring resistance to leaf rust, stripe rust, and stem rust, but Jagger lost the allelic region containing the resistance gene on QHf.osu-2A against Hessian fly. The molecular marker can greatly facilitate the identification of such a candidate gene for the QHf.osu-1A locus, enabling opportunities for functional analysis which can then be used to better understand the defense mechanisms of wheat plants in response to this major insect.
A Kansas population of Hessian fly consisting of predominantly biotype GP  was used in this study. The population was maintained on Hessian fly susceptible wheat seedlings (‘Karl 92’) in the greenhouse. Karl 92 was used as a susceptible control and WGRC42 was used as a resistant control to test the Hessian fly population .
Based on observations on seedling plants, Jagger is susceptible to Hessian fly biotype GP infestation whereas 2174 is resistant; therefore, the Jagger × 2174 population of 154 F6:8 recombinant inbred lines (RILs) were used to map genes associated with resistance to Hessian fly. Wheat genomic DNA was extracted from leaf tissue of each RIL plant as described previously .
Parental lines and 154 F6:8 RILs were evaluated for phenotypic reaction to Hessian fly infestation in growth chambers at 18 ± 1°C with a 14 h:10 h (light: dark) photoperiod as described previously [20, 34, 69]. Briefly, seedling plants were infested with mated Hessian fly females. Three weeks after infestation, seedling plants were examined to identify resistant and susceptible phenotypes. Susceptible plants were stunted, dark green, and harbored live larvae, whereas resistant plants grew normally with light green color and dead larvae . Hessian fly reactions were recorded and a random subset of the RILs was confirmed with additional replications (three totally).
Resistance screening was carried out in greenhouse as described previously . Specifically, approximately 20 seeds of each wheat line were planted in a flat. A flat with soil was divided into 22 rows with a line divider, and each row was further divided into two half rows. Therefore in each flat, a total of 44 wheat lines could be planted. Along with wheat lines from a mapping population, two susceptible control wheat lines, Karl 92 and Jagger (the susceptible parent), and two resistant control wheat lines, 2174 (putative resistant donor) and WGRC42 (containing the R gene Hdic), were also planted in each flat. Wheat plants were infested at one leaf-stage (the second leaf just emerged) with Hessian fly adults under a cheese cloth tent. Females deposit eggs on the adaxial surface of the first leaf. Infestation was stopped when the egg density reaches 8 per leaf, which usually results in ~5 larvae per plant. Resistant and susceptible plants were phenotyped three weeks after infestation. Phenotypes of resistant and susceptible plants were distinct and easy to score. Resistant plants grew normally with light green color and susceptible plants were stunt with dark-green color. Plants with a resistant phenotype were further dissected to check for dead larvae. Plants without dead larvae were taken as escapes and were excluded from the data set.
A total of 404 SSR markers were assembled in linkage groups using MapMaker 3.0 program, and a QTL was identified using the WinQTLCart 2.5 program as previously described [37, 70]. Other unlinked markers were tested using correlation analysis to examine which markers might be related to Hessian fly resistance using SAS software (SAS 9.1, SAS Institute Inc. Cary, NC, USA).
Primers used for detecting allelic variation in TaOPR1-A, TaOPR2-A, TaOPR7-B, TaLOX1-A, TaLOX2-B, TaSrk, TaClp , and TaPkc in hexaploid wheat
PCR size (bp)
T m (°C)
All PCR reactions were performed in a 25 μl reaction under the following conditions: 1 denaturation cycle at 94°C for 5 min; 40 cycles of 94°C for 30 s, 52-62°C (based on primer annealing temperature) for 30 s, and 72°C for 30 s – 2.5 min (based on fragment length), then a final extension step at 72°C for 10 min before cooling to 4°C. PCR amplified fragments were separated on 2% agarose gels through electrophoresis in 1× TAE buffer or 6% acrylamide gel in 0.5× TBE buffer. DNA banding patterns were visualized under UV light with ethidium bromide staining.
12-oxo-phytodienoic acid reductase
Recombinant inbred lines
Quantitative trait loci
Simple sequence repeat
Allene oxide synthase
12-oxo-phytodienoic acid reductase 3
Expressed sequence tag
Bacterial artificial chromosome
Induced receptor-like kinase
Catalytic domain of protein kinases
This work was supported by USDA-NIFA T-CAP grant No. 2011-68002-30029, the Oklahoma Center of Advanced Science and Technology (OCAST, PAS07-002), and the Oklahoma Wheat Research Foundation.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.