Microbe-dependent heterosis in maize.

Hybrids account for nearly all commercially planted varieties of maize and many other crop plants because crosses between inbred lines of these species produce first-generation [F1] offspring that greatly outperform their parents. The mechanisms underlying this phenomenon, called heterosis or hybrid vigor, are not well understood despite over a century of intensive research. The leading hypotheses-which focus on quantitative genetic mechanisms (dominance, overdominance, and epistasis) and molecular mechanisms (gene dosage and transcriptional regulation)-have been able to explain some but not all of the observed patterns of heterosis. Abiotic stressors are known to impact the expression of heterosis; however, the potential role of microbes in heterosis has largely been ignored. Here, we show that heterosis of root biomass and other traits in maize is strongly dependent on the belowground microbial environment. We found that, in some cases, inbred lines perform as well by these criteria as their F1 offspring under sterile conditions but that heterosis can be restored by inoculation with a simple community of seven bacterial strains. We observed the same pattern for seedlings inoculated with autoclaved versus live soil slurries in a growth chamber and for plants grown in steamed or fumigated versus untreated soil in the field. In a different field site, however, soil steaming increased rather than decreased heterosis, indicating that the direction of the effect depends on community composition, environment, or both. Together, our results demonstrate an ecological phenomenon whereby soil microbes differentially impact the early growth of inbred and hybrid maize.

Characterization of a Host-Specific Toxic Activity Produced by Bipolaris cookei, Causal Agent of Target Leaf Spot of Sorghum.

Target leaf spot (TLS) of sorghum, caused by the necrotrophic fungus Bipolaris cookei, can cause severe yield loss in many parts of the world. We grew B. cookei in liquid culture and observed that the resulting culture filtrate (CF) was differentially toxic when infiltrated into the leaves of a population of 288 diverse sorghum lines. In this population, we found a significant correlation between high CF sensitivity and susceptibility to TLS. This suggests that the toxin produced in culture may play a role in the pathogenicity of B. cookei in the field. We demonstrated that the toxic activity is light sensitive and, surprisingly, insensitive to pronase, suggesting that it is not proteinaceous. We identified the two sorghum genetic loci most associated with the response to CF in this population. Screening seedlings with B. cookei CF could be a useful approach for prescreening germplasm for TLS resistance.

Maize metacaspases modulate the defense response mediated by the NLR protein Rp1-D21 likely by affecting its subcellular localization.

Plants usually employ resistance (R) genes to defend against the infection of pathogens, and most R genes encode intracellular nucleotide-binding, leucine-rich repeat (NLR) proteins. The recognition between R proteins and their cognate pathogens often triggers a rapid localized cell death at the pathogen infection sites, termed the hypersensitive response (HR). Metacaspases (MCs) belong to a cysteine protease family, structurally related to metazoan caspases. MCs play crucial roles in plant immunity. However, the underlying molecular mechanism and the link between MCs and NLR-mediated HR are not clear. In this study, we systematically investigated the MC gene family in maize and identified 11 ZmMCs belonging to two types. Further functional analysis showed that the type I ZmMC1 and ZmMC2, but not the type II ZmMC9, suppress the HR-inducing activity of the autoactive NLR protein Rp1-D21 and of its N-terminal coiled-coil (CCD21 ) signaling domain when transiently expressed in Nicotiana benthamiana. ZmMC1 and ZmMC2 physically associate with CCD21 in vivo. We further showed that ZmMC1 and ZmMC2, but not ZmMC9, are predominantly localized in a punctate distribution in both N. benthamiana and maize (Zea mays) protoplasts. Furthermore, the co-expression of ZmMC1 and ZmMC2 with Rp1-D21 and CCD21 causes their re-distribution from being uniformly distributed in the nucleocytoplasm to a punctate distribution co-localizing with ZmMC1 and ZmMC2. We reveal a novel role of plant MCs in modulating the NLR-mediated defense response and derive a model to explain it.

Multistatic fiber-based system for measuring the Mueller matrix bidirectional reflectance distribution function.

Bidirectionality effects can be a significant confounding factor when measuring hyperspectral reflectance data. The bidirectional reflectance distribution function (BRDF) can effectively characterize the reflectivity of surfaces to correct remote sensing measurements. However, measuring BRDFs can be time-consuming, especially when collecting Mueller matrix BRDF (mmBRDF) measurements of a surface via conventional goniometric techniques. In this paper, we present a system for collecting mmBRDF measurements using static optical fiber detectors that sample the hemisphere surrounding an object. The entrance to each fiber contains a polarization state analyzer configuration, allowing for the simultaneous acquisition of the Stokes vector intensity components at many altitudinal and azimuthal viewing positions. We describe the setup, calibration, and data processing used for this system and present its performance as applied to mmBRDF measurements of a ground glass diffuser.

Genotypic and phenotypic characterization of a large, diverse population of maize near-isogenic lines.

Genome-wide association (GWA) studies can identify quantitative trait loci (QTL) putatively underlying traits of interest, and nested association mapping (NAM) can further assess allelic series. Near-isogenic lines (NILs) can be used to characterize, dissect and validate QTL, but the development of NILs is costly. Previous studies have utilized limited numbers of NILs and introgression donors. We characterized a panel of 1270 maize NILs derived from crosses between 18 diverse inbred lines and the recurrent inbred parent B73, referred to as the nested NILs (nNILs). The nNILs were phenotyped for flowering time, height and resistance to three foliar diseases, and genotyped with genotyping-by-sequencing. Across traits, broad-sense heritability (0.4-0.8) was relatively high. The 896 genotyped nNILs contain 2638 introgressions, which span the entire genome with substantial overlap within and among allele donors. GWA with the whole panel identified 29 QTL for height and disease resistance with allelic variation across donors. To date, this is the largest and most diverse publicly available panel of maize NILs to be phenotypically and genotypically characterized. The nNILs are a valuable resource for the maize community, providing an extensive collection of introgressions from the founders of the maize NAM population in a B73 background combined with data on six agronomically important traits and from genotyping-by-sequencing. We demonstrate that the nNILs can be used for QTL mapping and allelic testing. The majority of nNILs had four or fewer introgressions, and could readily be used for future fine mapping studies.

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.

Correction: Molecular and functional analyses of a maize autoactive NB-LRR protein identify precise structural requirements for activity.

[This corrects the article DOI: 10.1371/journal.ppat.1004674.].

Maize Homologs of CCoAOMT and HCT, Two Key Enzymes in Lignin Biosynthesis, Form Complexes with the NLR Rp1 Protein to Modulate the Defense Response.

Disease resistance (R) genes encode nucleotide binding Leu-rich-repeat (NLR) proteins that confer resistance to specific pathogens. Upon pathogen recognition they trigger a defense response that usually includes a so-called hypersensitive response (HR), a rapid localized cell death at the site of pathogen infection. Intragenic recombination between two maize (Zea mays) NLRs, Rp1-D and Rp1-dp2, resulted in the formation of a hybrid NLR, Rp1-D21, which confers an autoactive HR in the absence of pathogen infection. From a previous quantitative trait loci and genome-wide association study, we identified genes encoding two key enzymes in lignin biosynthesis, hydroxycinnamoyltransferase (HCT) and caffeoyl CoA O-methyltransferase (CCoAOMT), adjacent to the nucleotide polymorphisms that were highly associated with variation in the severity of Rp1-D21-induced HR We have previously shown that the two maize HCT homologs suppress the HR conferred by Rp1-D21 in a heterologous system, very likely through physical interaction. Here, we show, similarly, that CCoAOMT2 suppresses the HR induced by either the full-length or by the N-terminal coiled-coil domain of Rp1-D21 also likely via physical interaction and that the metabolic activity of CCoAOMT2 is unlikely to be necessary for its role in suppressing HR. We also demonstrate that CCoAOMT2, HCTs, and Rp1 proteins can form in the same complexes. A model is derived to explain the roles of CCoAOMT and HCT in Rp1-mediated defense resistance.

The Genetics of Leaf Flecking in Maize and Its Relationship to Plant Defense and Disease Resistance.

Physiological leaf spotting, or flecking, is a mild-lesion phenotype observed on the leaves of several commonly used maize (Zea mays) inbred lines and has been anecdotally linked to enhanced broad-spectrum disease resistance. Flecking was assessed in the maize nested association mapping (NAM) population, comprising 4,998 recombinant inbred lines from 25 biparental families, and in an association population, comprising 279 diverse maize inbreds. Joint family linkage analysis was conducted with 7,386 markers in the NAM population. Genome-wide association tests were performed with 26.5 million single-nucleotide polymorphisms (SNPs) in the NAM population and with 246,497 SNPs in the association population, resulting in the identification of 18 and three loci associated with variation in flecking, respectively. Many of the candidate genes colocalizing with associated SNPs are similar to genes that function in plant defense response via cell wall modification, salicylic acid- and jasmonic acid-dependent pathways, redox homeostasis, stress response, and vesicle trafficking/remodeling. Significant positive correlations were found between increased flecking, stronger defense response, increased disease resistance, and increased pest resistance. A nonlinear relationship with total kernel weight also was observed whereby lines with relatively high levels of flecking had, on average, lower total kernel weight. We present evidence suggesting that mild flecking could be used as a selection criterion for breeding programs trying to incorporate broad-spectrum disease resistance.

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.

Maize Homologs of Hydroxycinnamoyltransferase, a Key Enzyme in Lignin Biosynthesis, Bind the Nucleotide Binding Leucine-Rich Repeat Rp1 Proteins to Modulate the Defense Response.

In plants, most disease resistance genes encode nucleotide binding Leu-rich repeat (NLR) proteins that trigger a rapid localized cell death called a hypersensitive response (HR) upon pathogen recognition. The maize (Zea mays) NLR protein Rp1-D21 derives from an intragenic recombination between two NLRs, Rp1-D and Rp1-dp2, and confers an autoactive HR in the absence of pathogen infection. From a previous quantitative trait loci and genome-wide association study, we identified a single-nucleotide polymorphism locus highly associated with variation in the severity of Rp1-D21-induced HR. Two maize genes encoding hydroxycinnamoyltransferase (HCT; a key enzyme involved in lignin biosynthesis) homologs, termed HCT1806 and HCT4918, were adjacent to this single-nucleotide polymorphism. Here, we show that both HCT1806 and HCT4918 physically interact with and suppress the HR conferred by Rp1-D21 but not other autoactive NLRs when transiently coexpressed in Nicotiana benthamiana. Other maize HCT homologs are unable to confer the same level of suppression on Rp1-D21-induced HR. The metabolic activity of HCT1806 and HCT4918 is unlikely to be necessary for their role in suppressing HR. We show that the lignin pathway is activated by Rp1-D21 at both the transcriptional and metabolic levels. We derive a model to explain the roles of HCT1806 and HCT4918 in Rp1-mediated disease resistance.

A Genome-Wide Association Study for Partial Resistance to Maize Common Rust.

Quantitative resistance to maize common rust (causal agent Puccinia sorghi) was assessed in an association mapping population of 274 diverse inbred lines. Resistance to common rust was found to be moderately correlated with resistance to three other diseases and with the severity of the hypersensitive defense response previously assessed in the same population. Using a mixed linear model accounting for the confounding effects of population structure and flowering time, genome-wide association tests were performed based at 246,497 single-nucleotide polymorphism loci. Three loci associated with maize common rust resistance were identified. Candidate genes at each locus had predicted roles, mainly in cell wall modification. Other candidate genes included a resistance gene and a gene with a predicted role in regulating accumulation of reactive oxygen species.

Cytoplasmic and Nuclear Localizations Are Important for the Hypersensitive Response Conferred by Maize Autoactive Rp1-D21 Protein.

Disease resistance (R) genes have been isolated from many plant species. Most encode nucleotide binding leucine-rich repeat (NLR) proteins that trigger a rapid localized programmed cell death called the hypersensitive response (HR) upon pathogen recognition. Despite their structural similarities, different NLR are distributed in a range of subcellular locations, and analogous domains play diverse functional roles. The autoactive maize NLR gene Rp1-D21 derives from an intragenic recombination between two NLR genes, Rp1-D and Rp1-dp2, and confers a HR independent of the presence of a pathogen. Rp1-D21 and its N-terminal coiled coil (CC) domain (CCD21) confer autoactive HR when transiently expressed in Nicotiana benthamiana. Rp1-D21 was predominantly localized in cytoplasm with a small amount in the nucleus, while CCD21 was localized in both nucleus and cytoplasm. Targeting of Rp1-D21 or CCD21 predominantly to either the nucleus or the cytoplasm abolished HR-inducing activity. Coexpression of Rp1-D21 or CCD21 constructs confined, respectively, to the nucleus and cytoplasm did not rescue full activity, suggesting nucleocytoplasmic movement was important for HR induction. This work emphasizes the diverse structural and subcellular localization requirements for activity found among plant NLR R genes.

Limits on the reproducibility of marker associations with southern leaf blight resistance in the maize nested association mapping population.

BACKGROUND: A previous study reported a comprehensive quantitative trait locus (QTL) and genome wide association study (GWAS) of southern leaf blight (SLB) resistance in the maize Nested Association Mapping (NAM) panel. Since that time, the genomic resources available for such analyses have improved substantially. An updated NAM genetic linkage map has a nearly six-fold greater marker density than the previous map and the combined SNPs and read-depth variants (RDVs) from maize HapMaps 1 and 2 provided 28.5 M genomic variants for association analysis, 17 fold more than HapMap 1. In addition, phenotypic values of the NAM RILs were re-estimated to account for environment-specific flowering time covariates and a small proportion of lines were dropped due to genotypic data quality problems. Comparisons of original and updated QTL and GWAS results confound the effects of linkage map density, GWAS marker density, population sample size, and phenotype estimates. Therefore, we evaluated the effects of changing each of these parameters individually and in combination to determine their relative impact on marker-trait associations in original and updated analyses. RESULTS: Of the four parameters varied, map density caused the largest changes in QTL and GWAS results. The updated QTL model had better cross-validation prediction accuracy than the previous model. Whereas joint linkage QTL positions were relatively stable to input changes, the residual values derived from those QTL models (used as inputs to GWAS) were more sensitive, resulting in substantial differences between GWAS results. The updated NAM GWAS identified several candidate genes consistent with previous QTL fine-mapping results. CONCLUSIONS: The highly polygenic nature of resistance to SLB complicates the identification of causal genes. Joint linkage QTL are relatively stable to perturbations of data inputs, but their resolution is generally on the order of tens or more Mbp. GWAS associations have higher resolution, but lower power due to stringent thresholds designed to minimize false positive associations, resulting in variability of detection across studies. The updated higher density linkage map improves QTL estimation and, along with a much denser SNP HapMap, greatly increases the likelihood of detecting SNPs in linkage with causal variants. We recommend use of the updated genetic resources and results but emphasize the limited repeatability of small-effect associations.

Multivariate analysis of maize disease resistances suggests a pleiotropic genetic basis and implicates a GST gene.

Plants are attacked by pathogens representing diverse taxonomic groups, such that genes providing multiple disease resistance (MDR) are expected to be under positive selection pressure. To address the hypothesis that naturally occurring allelic variation conditions MDR, we extended the framework of structured association mapping to allow for the analysis of correlated complex traits and the identification of pleiotropic genes. The multivariate analytical approach used here is directly applicable to any species and set of traits exhibiting correlation. From our analysis of a diverse panel of maize inbred lines, we discovered high positive genetic correlations between resistances to three globally threatening fungal diseases. The maize panel studied exhibits rapidly decaying linkage disequilibrium that generally occurs within 1 or 2 kb, which is less than the average length of a maize gene. The positive correlations therefore suggested that functional allelic variation at specific genes for MDR exists in maize. Using a multivariate test statistic, a glutathione S-transferase (GST) gene was found to be associated with modest levels of resistance to all three diseases. Resequencing analysis pinpointed the association to a histidine (basic amino acid) for aspartic acid (acidic amino acid) substitution in the encoded protein domain that defines GST substrate specificity and biochemical activity. The known functions of GSTs suggested that variability in detoxification pathways underlie natural variation in maize MDR.

The quantification of southern corn leaf blight disease using deep UV fluorescence spectroscopy and autoencoder anomaly detection techniques.

Southern leaf blight (SLB) is a foliar disease caused by the fungus Cochliobolus heterostrophus infecting maize plants in humid, warm weather conditions. SLB causes production losses to corn producers in different regions of the world such as Latin America, Europe, India, and Africa. In this paper, we demonstrate a non-destructive method to quantify the signs of fungal infection in SLB-infected corn plants using a deep UV (DUV) fluorescence spectrometer, with a 248.6 nm excitation wavelength, to acquire the emission spectra of healthy and SLB-infected corn leaves. Fluorescence emission spectra of healthy and diseased leaves were used to train an Autoencoder (AE) anomaly detection algorithm-an unsupervised machine learning model-to quantify the phenotype associated with SLB-infected leaves. For all samples, the signature of corn leaves consisted of two prominent peaks around 450 nm and 325 nm. However, SLB-infected leaves showed a higher response at 325 nm compared to healthy leaves, which was correlated to the presence of C. heterostrophus based on disease severity ratings from Visual Scores (VS). Specifically, we observed a linear inverse relationship between the AE error and the VS (R2 = 0.94 and RMSE = 0.935). With improved hardware, this method may enable improved quantification of SLB infection versus visual scoring based on e.g., fungal spore concentration per unit area and spatial localization.

The interaction between abiotic and biotic soil factors drive heterosis expression in maize.

Heterosis or hybrid vigor refers to the superior phenotypes of hybrids relative to their parental inbred lines. Recently, soil microbes were identified as an environmental driver of maize heterosis. While manipulation of the soil microbial community consistently altered heterosis, the direction of the effect appeared to be dependent on the microbiome composition, environment, or both. Abiotic factors are well-known modifiers of heterosis expression, however, how the interactive effects between the soil microbial community and abiotic factors contribute to heterosis are poorly understood. To disentangle the proposed mechanisms by which microbes influence heterosis, we characterize the variation in heterosis expression when maize was grown in soil inocula derived from active maize farms or prairies. While we did not observe consistent differences in heterosis among plants grown in these inocula, our observations reaffirm that microbial effects on heterosis are likely specific to the local microbial community. The introduction of a nutrient amendment resulted in greater heterosis expression in the presence of an agricultural inoculum but not a prairie inoculum. We also observed an effect of soil inocula and nutrient treatment on the composition of bacterial and fungal communities in the root endosphere. In addition, the interaction between soil and nutrient treatment significantly affected bacterial community composition, whereas fungal community composition was only marginally affected by this interaction. These results further suggest that the soil microbial community plays a role in maize heterosis expression but that the abiotic environment is likely a larger driver.

Analysis of quantitative disease resistance to southern leaf blight and of multiple disease resistance in maize, using near-isogenic lines.

Maize inbred lines NC292 and NC330 were derived by repeated backcrossing of an elite source of southern leaf blight (SLB) resistance (NC250P) to the SLB-susceptible line B73, with selection for SLB resistance among and within backcross families at each generation. Consequently, while B73 is very SLB susceptible, its sister lines NC292 and NC330 are both SLB resistant. Previously, we identified the 12 introgressions from NC250P that differentiate NC292 and NC330 from B73. The goals of this study were to determine the effects of each introgression on resistance to SLB and to two other foliar fungal diseases of maize, northern leaf blight and gray leaf spot. This was achieved by generating and testing a set of near isogenic lines carry single or combinations of just two or three introgressions in a B73 background. Introgressions 3B, 6A, and 9B (bins 3.03-3.04, 6.01, and 9.02-9.03) all conferred significant levels of SLB resistance in the field. Introgression 6A was the only introgression that had a significant effect on juvenile plant resistance to SLB. Introgressions 6A and 9B conferred resistance to multiple diseases.

Multivariate mixed linear model analysis of longitudinal data: an information-rich statistical technique for analyzing plant disease resistance.

ABSTRACT The mixed linear model (MLM) is an advanced statistical technique applicable to many fields of science. The multivariate MLM can be used to model longitudinal data, such as repeated ratings of disease resistance taken across time. In this study, using an example data set from a multi-environment trial of northern leaf blight disease on 290 maize lines with diverse levels of resistance, multivariate MLM analysis was performed and its utility was examined. In the population and environments tested, genotypic effects were highly correlated across disease ratings and followed an autoregressive pattern of correlation decay. Because longitudinal data are often converted to the univariate measure of area under the disease progress curve (AUDPC), comparisons between univariate MLM analysis of AUDPC and multivariate MLM analysis of longitudinal data were made. Univariate analysis had the advantage of simplicity and reduced computational demand, whereas multivariate analysis enabled a comprehensive perspective on disease development, providing the opportunity for unique insights into disease resistance. To aid in the application of multivariate MLM analysis of longitudinal data on disease resistance, annotated program syntax for model fitting is provided for the software ASReml.

Shades of gray: the world of quantitative disease resistance.

A thorough understanding of quantitative disease resistance (QDR) would contribute to the design and deployment of durably resistant crop cultivars. However, the molecular mechanisms that control QDR remain poorly understood, largely due to the incomplete and inconsistent nature of the resistance phenotype, which is usually conditioned by many loci of small effect. Here, we discuss recent advances in research on QDR. Based on inferences from analyses of the defense response and from the few isolated QDR genes, we suggest several plausible hypotheses for a range of mechanisms underlying QDR. We propose that a new generation of genetic resources, complemented by careful phenotypic analysis, will produce a deeper understanding of plant defense and more effective utilization of natural resistance alleles.