A Genome-Wide Association Study of the Maize Hypersensitive Defense Response Identifies Genes That Cluster in Related Pathways

Much remains unknown of molecular events controlling the plant hypersensitive defense response (HR), a rapid localized cell death that limits pathogen spread and is mediated by resistance (R-) genes. Genetic control of the HR is hard to quantify due to its microscopic and rapid nature. Natural modifiers of the ectopic HR phenotype induced by an aberrant auto-active R-gene (Rp1-D21), were mapped in a population of 3,381 recombinant inbred lines from the maize nested association mapping population. Joint linkage analysis was conducted to identify 32 additive but no epistatic quantitative trait loci (QTL) using a linkage map based on more than 7000 single nucleotide polymorphisms (SNPs). Genome-wide association (GWA) analysis of 26.5 million SNPs was conducted after adjusting for background QTL. GWA identified associated SNPs that colocalized with 44 candidate genes. Thirty-six of these genes colocalized within 23 of the 32 QTL identified by joint linkage analysis. The candidate genes included genes predicted to be in involved programmed cell death, defense response, ubiquitination, redox homeostasis, autophagy, calcium signalling, lignin biosynthesis and cell wall modification. Twelve of the candidate genes showed significant differential expression between isogenic lines differing for the presence of Rp1-D21. Low but significant correlations between HR-related traits and several previously-measured disease resistance traits suggested that the genetic control of these traits was substantially, though not entirely, independent. This study provides the first system-wide analysis of natural variation that modulates the HR response in plants.     

Use of virus-induced gene silencing to characterize genes involved in modulating hypersensitive cell death in maize

Plant disease resistance proteins (R-proteins) detect specific pathogen-derived molecules, triggering a defence response often including a rapid localized cell death at the point of pathogen penetration called the hypersensitive response (HR). The maize Rp1-D21 gene encodes a protein that triggers a spontaneous HR causing spots on leaves in the absence of any pathogen. Previously, we used fine mapping and functional analysis in a Nicotiana benthamiana transient expression system to identify and characterize a number of genes associated with variation in Rp1-D21-induced HR. Here we describe a system for characterizing genes mediating HR, using virus-induced gene silencing (VIGS) in a maize line carrying Rp1-D21. We assess the roles of 12 candidate genes. Three of these genes, SGT1, RAR1, and HSP90, are required for HR induced by a number of R-proteins across several plant–pathogen systems. We confirmed that maize HSP90 was required for full Rp1-D21-induced HR. However, suppression of SGT1 expression unexpectedly increased the severity of Rp1-D21-induced HR while suppression of RAR1 expression had no measurable effect. We confirmed the effects on HR of two genes we had previously validated in the N. benthamiana system, hydroxycinnamoyltransferase and caffeoyl CoA O-methyltransferase. We further showed the suppression the expression of two previously uncharacterized, candidate genes, IQ calmodulin binding protein (IQM3) and vacuolar protein sorting protein 37, suppressed Rp1-D21-induced HR. This approach is an efficient way to characterize the roles of genes modulating the hypersensitive defence response and other dominant lesion phenotypes in maize.     

Homologues of the maize rust resistance gene Rp1-D are genetically associated with a major rust resistance QTL in sorghum

As part of a comparative mapping study between sugarcane and sorghum, a sugarcane cDNA clone with homology to the maize Rp1-D rust resistance gene was mapped in sorghum. The cDNA probe hybridised to multiple loci, including one on sorghum linkage group (LG) E in a region where a major rust resistance QTL had been previously mapped. Partial sorghum Rp1-D homologues were isolated from genomic DNA of rust-resistant and -susceptible progeny selected from a sorghum mapping population. Sequencing of the Rp1-D homologues revealed five discrete sequence classes: three from resistant progeny and two from susceptible progeny. PCR primers specific to each sequence class were used to amplify products from the progeny and confirmed that the five sequence classes mapped to the same locus on LG E. Cluster analysis of these sorghum sequences and available sugarcane, maize and sorghum Rp1-D homologue sequences showed that the maize Rp1-D sequence and the partial sugarcane Rp1-D homologue were clustered with one of the sorghum resistant progeny sequence classes, while previously published sorghum Rp1-D homologue sequences clustered with the susceptible progeny sequence classes. Full-length sequence information was obtained for one member of a resistant progeny sequence class ( Rp1-SO) and compared with the maize Rp1-D sequence and a previously identified sorghum Rp1 homologue ( Rph1-2). There was considerable similarity between the two sorghum sequences and less similarity between the sorghum and maize sequences. These results suggest a conservation of function and gene sequence homology at the Rp1 loci of maize and sorghum and provide a basis for convenient PCR-based screening tools for putative rust resistance alleles in sorghum.     

A gene for resistance to the maize streak virus in the African CIMMYT maize inbred line CML202

Resistance to maize streak virus (MSV) is an essential trait of improved maize varieties in sub-Saharan Africa. We mapped quantitative trait loci (QTL) for resistance to MSV in a population of 196 F_2:3 lines derived from a cross between the maize inbred lines CML202 (resistant) from CIMMYT-Zimbabwe and Lo951 (susceptible) from Italy. Field tests were planted at two locations in Zimbabwe, inoculated with viruliferous leaf hoppers (Cicadulina mbila), and scored twice (21 and 83 days after infesting, DAI) on a 1–5 scale. The mean final streak intensity (score 2) of the parent lines was 2.2 (CML202) and 4.8 (Lo951). Genotype × location interaction was large for score 1 but negligible for score 2. Consequently, the heritability was higher for score 2 (0.93) than for score 1 (0.62). By composite interval mapping across locations, using a linkage map with 110 RFLP loci, four significant (LOD 3.0) QTL were identified for score 1 on chromosomes (C) 1, 2, 3, and 4, respectively. All four were contributed by CML202. For score 2, only the QTL on C 1 was significant (LOD =37), explaining 59% of the phenotypic and 64% of the genotypic variance. The QTL's partially dominant gene action was consistent with the nearly intermediate resistance of the F₁ generation (relative heterosis for resistance 12%). The presence of one major QTL is consistent with the bimodal frequency distribution of the mapping population showing a clear 3:1 segregation. This gene seems to be allelic or identical to Msv1, a major resistance gene which was previously identified in the same genomic region in Tzi4, an inbred line from IITA. Inbred CML202 had lower final disease ratings than Tzi4. The greater resistance of CML202 may be due to allelic differences at the msv1 locus or due to the minor QTL on C 2, 3, and 4 which were not detected in Tzi4.z y Trigo (International Maize and Wheat Improvement Center); IITA, International Institute of Tropical Agriculture; IRAT, Institute de Recherches Agronomiques Tropicales et des Cultures Vivrières; KARI, Kenya Agricultural Research Institute; MSV, maize streak virus; QTL, quantitative trait locus/loci     

Disease Lesion Mimicry Caused by Mutations in the Rust Resistance Gene rp1

The rp1 locus of maize controls race-specific resistance to the common rust fungus Puccinia sorghi. Four mutant or recombinant Rp1 alleles (rp1-NC3, Rp1-D21, Rp1-MD19, and Rp1-Kr1N) were identified. They condition necrotic phenotypes in the absence of the rust pathogen. These Rp1 lesion mimics fall into three different phenotypic classes: (1) The rp1-NC3 and Rp1-D21 alleles require rust infection or other biotic stimulus to initiate necrotic lesions. These alleles react strongly to all maize rust biotypes tested and also to nonhost rusts. (2) The Rp1-MD19 allele, which has a similar phenotype, also requires a biotic stimulus to initiate lesions. However, Rp1-MD19 shows the race specificity of the Rp1-D gene. (3) The Rp1-Kr1N allele specifies a diffuse necrotic phenotype in the absence of any biotic stimulus and a race-specific reaction when inoculated with maize rust.     

On the use of mathematically-derived traits in QTL mapping

Mathematically-derived traits from two or more component traits, either by addition, subtraction, multiplication, or division, have been frequently used in genetics and breeding. When used in quantitative trait locus (QTL) mapping, derived traits sometimes show discrepancy with QTL identified for the component traits. We used three QTL distributions and three genetic effects models, and an actual maize mapping population, to investigate the efficiency of using derived traits in QTL mapping, and to understand the genetic and biological basis of derived-only QTL, i.e., QTL identified for a derived trait but not for any component trait. Results indicated that the detection power of the four putative QTL was consistently greater than 90% for component traits in simulated populations, each consisting of 200 recombinant inbred lines. Lower detection power and higher false discovery rate (FDR) were observed when derived traits were used. In an actual maize population, simulations were designed based on the observed QTL distributions and effects. When derived traits were used, QTL detected for both component and derived traits had comparable power, but those detected for component traits but not for derived traits had low detection power. The FDR from subtraction and division in the maize population were higher than the FDR from addition and multiplication. The use of derived traits increased the gene number, caused higher-order gene interactions than observed in component traits, and possibly complicated the linkage relationship between QTL as well. The increased complexity of the genetic architecture with derived traits may be responsible for the reduced detection power and the increased FDR. Derived-only QTL identified in practical genetic populations can be explained either as minor QTL that are not significant in QTL mapping of component traits, or as false positives.     

Mapping quantitative trait loci using naturally occurring genetic variance among commercial inbred lines of maize (Zea mays L.)

Many commercial inbred lines are available in crops. A large amount of genetic variation is preserved among these lines. The genealogical history of the inbred lines is usually well documented. However, quantitative trait loci (QTL) responsible for the genetic variances among the lines are largely unexplored due to lack of statistical methods. In this study, we show that the pedigree information of the lines along with the trait values and marker information can be used to map QTL without the need of further crossing experiments. We develop a Monte Carlo method to estimate locus-specific identity-by-descent (IBD) matrices. These IBD matrices are further incorporated into a mixed-model equation for variance component analysis. QTL variance is estimated and tested at every putative position of the genome. The actual QTL are detected by scanning the entire genome. Applying this new method to a well-documented pedigree of maize (Zea mays L.) that consists of 404 inbred lines, we mapped eight QTL for the maize male flowering trait, growing degree day heat units to pollen shedding (GDUSHD). These detected QTL contributed >80% of the variance observed among the inbred lines. The QTL were then used to evaluate all the inbred lines using the best linear unbiased prediction (BLUP) technique. Superior lines were selected according to the estimated QTL allelic values, a technique called marker-assisted selection (MAS). The MAS procedure implemented via BLUP may be routinely used by breeders to select superior lines and line combinations for development of new cultivars.     

Genomic Dissection of Leaf Angle in Maize (Zea mays L.) Using a Four-Way Cross Mapping Population

Increasing grain yield by the selection for optimal plant architecture has been the key focus in modern maize breeding. As a result, leaf angle, an important determinant of plant architecture, has been significantly improved to adapt to the ever-increasing plant density in maize production over the past several decades. To extend our understanding on the genetic mechanisms of leaf angle in maize, we developed the first four-way cross mapping population, consisting of 277 lines derived from four maize inbred lines with varied leaf angles. The four-way cross mapping population together with the four parental lines were evaluated for leaf angle in two environments. In this study, we reported linkage maps built in the population and quantitative trait loci (QTL) on leaf angle detected by inclusive composite interval mapping (ICIM). ICIM applies a two-step strategy to effectively separate the cofactor selection from the interval mapping, which controls the background additive and dominant effects at the same time. A total of 14 leaf angle QTL were identified, four of which were further validated in near-isogenic lines (NILs). Seven of the 14 leaf angle QTL were found to overlap with the published leaf angle QTL or genes, and the remaining QTL were unique to the four-way population. This study represents the first example of QTL mapping using a four-way cross population in maize, and demonstrates that the use of specially designed four-way cross is effective in uncovering the basis of complex and polygenetic trait like leaf angle in maize.     

Precise mapping of quantitative trait loci for resistance to southern leaf blight, caused by Cochliobolus heterostrophus race O, and flowering time using advanced intercross maize lines

The intermated B73 × Mo17 (IBM) population, an advanced intercross recombinant inbred line population derived from a cross between the maize lines B73 (susceptible) and Mo17 (resistant), was evaluated in four environments for resistance to southern leaf blight (SLB) disease caused by Cochliobolus heterostrophus race O. Two environments were artificially inoculated, while two were not inoculated and consequently had substantially lower disease pressure. Four common SLB resistance quantitative trait loci (QTL) were identified in all environments, two in bin 3.04 and one each in bins 1.10 and 8.02/3. There was no significant correlation between disease resistance and days to anthesis. A direct comparison was made between SLB QTL detected in two populations, independently derived from the same parental cross: the IBM advanced intercross population and a conventional recombinant inbred line population. Several QTL for SLB resistance were detected in both populations, with the IBM providing between 5 and, in one case, 50 times greater mapping resolution.     

Boron toxicity in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.). I. Quantitative trait locus (QTL) analysis of tolerance to boron toxicity

Boron toxicity tolerance of rice plants was studied. Modern japonica subspecies such as Koshihikari, Nipponbare, and Sasanishiki were tolerant, whereas indica subspecies such as Kasalath and IR36 were intolerant to excessive application of boron (B), even though their shoot B contents under B toxicity were not significantly different. Recombinant inbred lines (RILs) of japonica Nekken-1 and indica IR36 were used for quantitative trait locus (QTL) analysis to identify the gene responsible for B toxicity tolerance. A major QTL that could explain 45% of the phenotypic variation was detected in chromosome 4. The QTL was confirmed using a population derived from a recombinant inbred line which is heterogenic at the QTL region. The QTL was also confirmed in other chromosome segment substitution lines (CSSLs).     

Maize leaf temperature responses to drought: Thermal imaging and quantitative trait loci (QTL) mapping

Leaf temperature has been shown to vary when plants are subjected to water stress conditions. Recent advances in infrared thermography have increased the probability of recording drought tolerant responses more accurately. The aims of this study were to identify the effects of drought on leaf temperature using infrared thermography. Furthermore, the genomic regions responsible for the expression of leaf temperature variation in maize seedlings (Zea mays L.) were explored. The maize inbred lines Zong3 and 87-1 were evaluated using infrared thermography and exhibited notable differences in leaf temperature response to water stress. Correlation analysis indicated that leaf temperature response to water stress played an integral role in maize biomass accumulation. Additionally, a mapping population of 187 recombinant inbred lines (RILs) derived from a cross between Zong3 and 87-1 was constructed to identify quantitative trait loci (QTL) responsible for physiological traits associated with seedling water stress. Leaf temperature differences (LTD) and the drought tolerance index (DTI) of shoot fresh weight (SFW) and shoot dry weight (SDW) were the traits evaluated for QTL analysis in maize seedlings. A total of nine QTL were detected by composite interval mapping (CIM) for the three traits (LTD, RSFW and RSDW). Two co-locations responsible for both RSFW and RSDW were detected on chromosomes 1 and 2, respectively, which showed common signs with their trait correlations. Another co-location was detected on chromosome 9 between LTD and shoot biomass, which provided genetic evidence that leaf temperature affects biomass accumulation. Additionally, the utility of a thermography system for drought tolerance breeding in maize was discussed.     

Molecular and Functional Analyses of a Maize Autoactive NB-LRR Protein Identify Precise Structural Requirements for Activity

Plant disease resistance is often mediated by nucleotide binding-leucine rich repeat (NLR) proteins which remain auto-inhibited until recognition of specific pathogen-derived molecules causes their activation, triggering a rapid, localized cell death called a hypersensitive response (HR). Three domains are recognized in one of the major classes of NLR proteins: a coiled-coil (CC), a nucleotide binding (NB-ARC) and a leucine rich repeat (LRR) domains. The maize NLR gene Rp1-D21 derives from an intergenic recombination event between two NLR genes, Rp1-D and Rp1-dp2 and confers an autoactive HR. We report systematic structural and functional analyses of Rp1 proteins in maize and N. benthamiana to characterize the molecular mechanism of NLR activation/auto-inhibition. We derive a model comprising the following three main features: Rp1 proteins appear to self-associate to become competent for activity. The CC domain is signaling-competent and is sufficient to induce HR. This can be suppressed by the NB-ARC domain through direct interaction. In autoactive proteins, the interaction of the LRR domain with the NB-ARC domain causes de-repression and thus disrupts the inhibition of HR. Further, we identify specific amino acids and combinations thereof that are important for the auto-inhibition/activity of Rp1 proteins. We also provide evidence for the function of MHD2, a previously uncharacterized, though widely conserved NLR motif. This work reports several novel insights into the precise structural requirement for NLR function and informs efforts towards utilizing these proteins for engineering disease resistance.     

Detection and integration of quantitative trait loci for grain yield components and oil content in two connected recombinant inbred line populations of high-oil maize

Improvement in grain yield is an important objective in high-oil maize breeding. In this study, one high-oil maize inbred was crossed with two normal maize inbreds to produce two connected recombinant inbred line (RIL) populations with 282 and 263 F_7:8 families, respectively. The field experiments were conducted under four environments, and eight grain yield components and grain oil content were evaluated. Two genetic linkage maps were constructed using 216 and 208 polymorphic SSR markers. Quantitative trait loci (QTL) were detected for all traits under each environment and in combined analysis. Meta-analysis was used to integrate genetic maps and detected QTL in both populations. A total of 199 QTL were detected, 122 in population 1 and 87 in population 2. Seven, 11 and 19 QTL showed consistency across five environments, across two RIL populations and with respective F_2:3 generations, respectively. 183 QTL were integrated in 28 meta-QTL (mQTL). QTL with contributions over 15% were consistently detected in 3–4 cases and integrated in mQTL. Each mQTL included 3–19 QTL related to 1–4 traits, reflecting remarkable QTL co-location for grain yield components and oil content. Further research and marker-assisted selection (MAS) should be concentrated on 37 consistent QTL and four genetic regions of mQTL with more than 10 QTL at bins 3.04–3.05, 7.02, 8.04–8.05 and 9.04–9.05. Near-isogenic lines for 100-grain-weight QTL at bin 7.02–7.03, for ear-length QTL at bin 7.02–7.03 and for rows-per-ear QTL at bin 3.08 are now in construction using MAS. Co-located candidate genes could facilitate the identification of candidate genes for grain yield in maize.     

Genetic relationships between resistance to stalk-tunneling by the European corn borer and cell-wall components in maize population B73×B52

The objective of this study was to assess the relationships among quantitative trait loci (QTL) detected for European corn borer (ECB) tunneling and cell-wall components (CWC) neutral detergent fiber (NDF), acid detergent fiber (ADF), and acid detergent lignin (ADL) content in leaf-sheath and stalk tissues in a maize recombinant inbred line population derived from inbred lines B73 and B52. Most of the QTL for ECB resistance (10/13) were at QTL positions for one or more CWC. Of the 12 QTL for NDF and ADF in leaf-sheaths, five for each trait were at or near QTL for ECB tunneling. Four of these five QTL for NDF and ADF mapped to common locations. Four of the eight leaf-sheath ADL QTL were detected in the same genomic regions as ECB QTL. For stalk tissue, four regions contained common/overlapping QTL for ECB tunneling, NDF, and ADF. Six such regions were observed for stalk ADL and ECB tunneling. Seven of the ten QTL associated with both CWC and ECB tunneling contributed to the negative correlations observed between these traits, while relatively few QTL effects were positively correlated. This suggests that while CWC contribute to ECB resistance in this population, other mechanisms and other genes also are involved. Several QTL contributing to the negative correlations between ECB tunneling and CWC in the leaf-sheaths mapped to similar positions as QTL detected in tropical maize populations for resistance to leaf-feeding by Diatraea grandiosella Dyar and Diatraea saccharalis Fabricus. These regions may contain genes involved in the synthesis of cellulose, hemicellulose, and lignin in the leaf-blades and leaf-sheaths of maize.     

Identification of novel QTL contributing resistance to aflatoxin accumulation in maize

The toxic metabolic product aflatoxin produced by the opportunistic fungus Aspergillus flavus (Link:Fr) in maize (Zea mays L.) can cause disease and economic harm when levels exceed very minute quantities. The selection of resistant germplasm has great potential to reduce the problem, but the highly quantitative nature of the trait makes this a difficult endeavor. The identification of aflatoxin accumulation resistance quantitative trait loci (QTL) from resistant donor lines and the discovery of linked markers could speed this task. To identify marker–trait associations for marker-assisted breeding, a genetic mapping population of F_2:3 families was developed from Mp715, a maize inbred line resistant to aflatoxin accumulation, and T173, a susceptible, southern-adapted maize inbred line. QTL, some with large phenotypic effects, were identified in multiple years on chromosomes 1, 3, 5, and 10, and smaller QTL identified in only 1 year were found on chromosomes 4 and 9. The phenotypic effect of each QTL ranged from 2.7 to 18.5%, and models created with multiple QTL could explain up to 45.7% of the phenotypic variation across years, indicating that the variation associated with the trait can be manipulated using molecular markers.     

Allelic and haplotypic diversity at the rp1 rust resistance locus of maize

The maize Rp1 rust resistance locus is a complex consisting of a family of closely related resistance genes. The number of Rp1 paralogs in different maize lines (haplotypes) varied from a single gene in some stocks of the inbred A188 to >50 genes in haplotypes carrying the Rp1-A and Rp1-H specificities. The sequences of paralogs in unrelated haplotypes differ, indicating that the genetic diversity of Rp1-related genes is extremely broad in maize. Two unrelated haplotypes with five or nine paralogs had identical resistance phenotypes (Rp1-D) encoded in genes that differed by three nucleotides resulting in a single amino acid substitution. Genes in some haplotypes are more similar to each other than to any of the genes in other haplotypes indicating that they are evolving in a concerted fashion.