QTL underlying the resistance to soybean aphid (Aphis glycines Matsumura) through isoflavone-mediated antibiosis in soybean cultivar 'Zhongdou 27'

Fine-Mapping and Functional Analyses of a Candidate Gene Controlling Isoflavone Content in Soybeans Seed

Isoflavones, one of the most important secondary metabolites produced by soybeans (Glycine max (L.) Merr.), are important for a variety of biological processes, and are beneficial for human health. To identify genetic loci underlying soybean isoflavone content, a mapping population containing 119 F_5:18 recombinant inbred lines, derived by crossing soybean cultivar "Zhongdou27" with "Dongong8004," was used. We identified 15 QTLs associated with isoflavone contents. A novel loci, qISO19-1, was mapped onto soybean chromosome 19 and was fine-mapped to a 62.8 kb region using a BC₂F₂ population. We considered GmMT1 as a candidate gene for the qISO19-1 locus due to the significant positive correlation recovered between its expression level and isoflavone content in the seeds of 43 soybean germplasms. Overexpression of GmMT1 in Arabidopsis and soybean cultivars increased isoflavone contents. Transgenic soybeans overexpressing GmMT1 also exhibited improved resistance to pathogenic infection, while transgenic Arabidopsis resisted salt and drought stress.

Investigation of Metabolic Resistance to Soybean Aphid (Aphis glycines Matsumura) Feeding in Soybean Cultivars

Simple Summary

This project examined the interaction between soybean aphids and Ontario-grown soybean cultivars to determine which leaf metabolites were most associated with aphid resistance. Tolerance and resistance were determined by measuring the growth and reproduction of aphids and leaf feeding damage over 10-day and 4-week infestation periods. Chromatographic techniques were used for the analysis of legume-specific plant natural products and isoflavonoids, and high-resolution mass spectrometry was used for the identification of free amino acids in aphid-resistant, tolerant, and susceptible soybean cultivars. There was a low correlation between isoflavonoid leaf concentrations and aphid resistance in the soybean varieties studied; however, the aphid-resistant cultivars were determined to have lower free amino acid concentrations, indicating that lower nutrient quality could be responsible for the resistance observed. Identifying these cultivars is important for managing aphid populations and provides an additional tool for soybean integrated pest management.

Abstract

Soybean aphid (Aphis glycines) is a major soybean (Glycine max) herbivore pest in many soybean growing regions. High numbers of aphids on soybean can cause severe reductions in yield. The management of soybean aphids includes monitoring, insecticide applications when required, and the use of resistant cultivars. Soybean aphid-resistant soybean varieties are associated with genes that confer one or more categories of resistance to soybean aphids, including antibiosis (affects survival, growth, and fecundity), antixenosis (affects behaviour such as feeding), and tolerance (plant can withstand greater damage without economic loss). The genetic resistance of soybean to several herbivores has been associated with isoflavonoid phytoalexins; however, this correlation has not been observed in soybean varieties commonly grown in southern Ontario, Canada. Isoflavonoids in the leaves of 18 cultivars in the early growth stage were analyzed by HPLC and the concentration by fresh weight was used to rate the potential resistance to aphids. Greenhouse and growth cabinet trials determined that the cultivars with greater resistance to aphids were Harosoy 63 and OAC Avatar. The most susceptible cultivar was Maple Arrow, whereas Pagoda and Conrad were more tolerant to aphid feeding damage. Overall, there was a low correlation between the number of aphids per leaf, feeding damage, and leaf isoflavonoid levels. Metabolite profiling by high-resolution LC-MS determined that the most resistant cultivar had on average lower levels of certain free amino acids (Met, Tyr, and His) relative to the most susceptible cultivar. This suggests that within the tested cultivars, nutritional quality stimulates aphid feeding more than isoflavonoids negatively affect aphid feeding or growth. These findings provide a better understanding of soybean host plant resistance and suggest ways to improve soybean resistance to aphid feeding through the breeding or metabolic engineering of leaf metabolites.

Exploring the Dynamics of Virulent and Avirulent Aphids: A Case for a 'Within Plant' Refuge

The soybean aphid, Aphis glycines (Hemiptera: Aphididae), is an invasive pest that can cause severe yield loss to soybeans in the North Central United States. A tactic to counter this pest is the use of aphid-resistant soybean varieties. However, the frequency of virulent biotypes that can survive on resistant varieties is expected to increase as more farmers use these varieties. Soybean aphids can alter soybean physiology primarily by two mechanisms, feeding facilitation, and the obviation of resistance, favoring subsequent colonization by additional conspecifics. We developed a nonlocal, differential equation population model to explore the dynamics of these biological mechanisms on soybean plants coinfested with virulent and avirulent aphids. We then use demographic parameters from laboratory experiments to perform numerical simulations via the model. We used this model to determine that initial conditions are an important factor in the season-long cooccurrence of both biotypes. The initial population of both biotypes above the resistance threshold or avirulent aphid close to resistance threshold and high virulent aphid population results in coexistence of the aphids throughout the season. These simulations successfully mimicked aphid dynamics observed in the field- and laboratory-based microcosms. The model showed an increase in colonization of virulent aphids increases the likelihood that aphid resistance is suppressed, subsequently increasing the survival of avirulent aphids. This interaction produced an indirect, positive interaction between the biotypes. These results suggest the potential for a 'within plant' refuge that could contribute to the sustainable use of aphid-resistant soybeans.

Interaction between Rag genes results in a unique synergistic transcriptional response that enhances soybean resistance to soybean aphids

BACKGROUND

Pyramiding different resistance genes into one plant genotype confers enhanced resistance at the phenotypic level, but the molecular mechanisms underlying this effect are not well-understood. In soybean, aphid resistance is conferred by Rag genes. We compared the transcriptional response of four soybean genotypes to aphid feeding to assess how the combination of Rag genes enhanced the soybean resistance to aphid infestation.

RESULTS

A strong synergistic interaction between Rag1 and Rag2, defined as genes differentially expressed only in the pyramid genotype, was identified. This synergistic effect in the Rag1/2 phenotype was very evident early (6 h after infestation) and involved unique biological processes. However, the response of susceptible and resistant genotypes had a large overlap 12 h after aphid infestation. Transcription factor (TF) analyses identified a network of interacting TF that potentially integrates signaling from Rag1 and Rag2 to produce the unique Rag1/2 response. Pyramiding resulted in rapid induction of phytochemicals production and deposition of lignin to strengthen the secondary cell wall, while repressing photosynthesis. We also identified Glyma.07G063700 as a novel, strong candidate for the Rag1 gene.

CONCLUSIONS

The synergistic interaction between Rag1 and Rag2 in the Rag1/2 genotype can explain its enhanced resistance phenotype. Understanding molecular mechanisms that support enhanced resistance in pyramid genotypes could facilitate more directed approaches for crop improvement.

Soybean (Glycine max L Merr) host-plant defenses and resistance to the two-spotted spider mite (Tetranychus urticae Koch)

In southern Ontario, Canada, the two-spotted spider mite (Tetranychus urticae) is an emerging pest of soybean (Glycine max) due to the increasing incidence of warmer, drier weather conditions. One key strategy to manage soybean pests is breeding resistant cultivars. Resistance to pathogens and herbivores in soybean has been associated with isoflavonoid phytoalexins, a group of specialized metabolites commonly associated with root, leaf and seed tissues. A survey of 18 Ontario soybean cultivars for spider mite resistance included evaluations of antibiosis and tolerance in relation to isoflavonoid and other metabolites detected in the leaves. Ten-day and 4-week trials beginning with early growth stage plants were used to compare survival, growth, fecundity as well as damage to leaves. Two-spotted spider mite (TSSM) counts were correlated with HPLC measurements of isoflavonoid concentration in the leaves and global metabolite profiling by high resolution LC-MS to identify other metabolites unique to the most resistant (R) and susceptible (S) cultivars. Within 10 days, no significant difference (P>0.05) in resistance to TSSM was determined between cultivars, but after 4 weeks, one cultivar, OAC Avatar, was revealed to have the lowest number of adult TSSMs and their eggs. Other cultivars showing partial resistance included OAC Wallace and OAC Lakeview, while Pagoda was the most tolerant to TSSM feeding. A low, positive correlation between isoflavonoid concentrations and TSSM counts and feeding damage indicated these compounds alone do not explain the range of resistance or tolerance observed. In contrast, other metabolite features were significantly different (P<0.05) in R versus S cultivars. In the presence of TSSM, the R cultivars had significantly greater (P<0.05) concentrations of the free amino acids Trp, Val, Thr, Glu, Asp and His relative to S cultivars. Furthermore, the R cultivar metabolites detected are viable targets for more in-depth analysis of their potential roles in TSSM defense.

The Soybean High Density 'Forrest' by 'Williams 82' SNP-Based Genetic Linkage Map Identifies QTL and Candidate Genes for Seed Isoflavone Content

Isoflavones are secondary metabolites that are abundant in soybean and other legume seeds providing health and nutrition benefits for both humans and animals. The objectives of this study were to construct a single nucleotide polymorphism (SNP)-based genetic linkage map using the ‘Forrest’ by ‘Williams 82’ (F×W82) recombinant inbred line (RIL) population (n = 306); map quantitative trait loci (QTL) for seed daidzein, genistein, glycitein, and total isoflavone contents in two environments over two years (NC-2018 and IL-2020); identify candidate genes for seed isoflavone. The FXW82 SNP-based map was composed of 2075 SNPs and covered 4029.9 cM. A total of 27 QTL that control various seed isoflavone traits have been identified and mapped on chromosomes (Chrs.) 2, 4, 5, 6, 10, 12, 15, 19, and 20 in both NC-2018 (13 QTL) and IL-2020 (14 QTL). The six QTL regions on Chrs. 2, 4, 5, 12, 15, and 19 are novel regions while the other 21 QTL have been identified by other studies using different biparental mapping populations or genome-wide association studies (GWAS). A total of 130 candidate genes involved in isoflavone biosynthetic pathways have been identified on all 20 Chrs. And among them 16 have been identified and located within or close to the QTL identified in this study. Moreover, transcripts from four genes (Glyma.10G058200, Glyma.06G143000, Glyma.06G137100, and Glyma.06G137300) were highly abundant in Forrest and Williams 82 seeds. The identified QTL and four candidate genes will be useful in breeding programs to develop soybean cultivars with high beneficial isoflavone contents.

Antixenosis in Glycine max (L.) Merr against Acyrthosiphon pisum (Harris)

To reveal the antixenosis potential against the pea aphid Acyrthosiphon pisum (Harris) (Hemiptera: Aphididae) we analyzed the pea aphid survival and probing behavior, and the quantitative and qualitative variation of flavonoids in the leaves of selected soybean Glycine max (L.) Merr (Fabaceae) cultivars 'Aldana', 'Annushka', 'Augusta', 'Madlen', 'Mavka', 'Simona', 'Violetta', and 'Viorica'. Aphid survival was drastically impeded on all cultivars. The electronic monitoring of aphid probing using the Electrical Penetration Graph (EPG) technique revealed that on all soybean cultivars, A. pisum readily probed into leaf tissues but the probes were usually terminated before reaching vascular tissues, which demonstrates the activity of antixenosis mechanisms in peripheral tissues epidermis and/or mesophyll in soybean leaves. The potency of antixenosis factors differed among soybean cultivars, which was reflected in differences in aphid survival and frequency and duration of phloem sap ingestion. Seven flavonoids were found: apigenin, daidzein, genistein, glycitein, isorhamnetin, kaempferol, and rutin, which occurred in different amount and proportion in individual cultivars. The content of apigenin and genistein in all soybean cultivars studied probably made them relatively unacceptable to A. pisum. Kaempferol in 'Aldana' might be responsible for the observed strong antixenosis resistance of this cultivar to A. pisum. The results of our survey provide the first detailed data that can be used for future studies.

The Source of Rag5-Mediated Resistance to Soybean Aphids Is Located in the Stem

The soybean aphid (Aphis glycines) continues to threaten soybean production in the United States. A suite of management strategies, such as planting aphid-resistant cultivars, has been successful in controlling soybean aphids. Several Rag genes (resistance against A. glycines) have been identified, and two are currently being deployed in commercial soybean cultivars. However, the mechanisms underlying Rag-mediated resistance are yet to be identified. In this study, we sought to determine the nature of resistance conferred by the Rag5 gene using behavioral, molecular biology, physiological, and biochemical approaches. We confirmed previous findings that plants carrying the Rag5 gene were resistant to soybean aphids in whole plant assays, and this resistance was absent in detached leaf assays. Analysis of aphid feeding behaviors using the electrical penetration graph technique on whole plants and detached leaves did not reveal differences between the Rag5 plants and Williams 82, a susceptible cultivar. In reciprocal grafting experiments, aphid populations were lower in the Rag5/rag5 (Scion/Root stock) chimera, suggesting that Rag5-mediated resistance is derived from the shoots. Further evidence for the role of stems comes from poor aphid performance in detached stem plus leaf assays. Gene expression analysis revealed that biosynthesis of the isoflavone kaempferol is upregulated in both leaves and stems in resistant Rag5 plants. Moreover, supplementing with kaempferol restored resistance in detached stems of plants carrying Rag5. This study demonstrates for the first time that Rag5-mediated resistance against soybean aphids is likely derived from stems.

Identification of a candidate gene associated with isoflavone content in soybean seeds using genome-wide association and linkage mapping

Isoflavone, a secondary metabolite produced by Glycine max (L.) Merr. (soybean), is valuable for human and plant health. The genetic architecture of soybean isoflavone content remains unclear, however, despite several mapping studies. We generated genomic data for 200 soybean cultivars and 150 recombinant inbred lines (RILs) to localize putative loci associated with soybean seed isoflavone content. Using a genome-wide association study (GWAS), we identified 87 single-nucleotide polymorphisms (SNPs) that were significantly associated with isoflavone concentration. Using linkage mapping, we identified 37 quantitative trait loci (QTLs) underlying the content of four isoflavones found in the RILs. A major locus on chromosome 8 (qISO8-1) was co-located by both the GWAS and linkage mapping. qISO8-1 was fine mapped to a 99.5-kb region, flanked by SSR_08_1651 and SSR_08_1656, in a BC2 F5 population. GmMPK1, encoding a mitogen-activated protein kinase, was identified as the causal gene in qISO8-1, and two natural GmMPK1 polymorphisms were significantly associated with isoflavone content. Overexpression of GmMPK1 in soybean hairy roots resulted in increased isoflavone concentrations. Overexpressing GmMPK1 in transgenic soybeans had greater resistance to Phytophthora root rot, suggesting that GmMPK1 might increase soybean resistance to biotic stress by influencing isoflavone content. Our results not only increase our understanding of the genetic architecture of soybean seed isoflavone content, but also provide a framework for the future marker-assisted breeding of high isoflavone content in soybean cultivars.

Role of reactive oxygen species and isoflavonoids in soybean resistance to the attack of the southern green stink bug

Southern green stink bugs (Nezara viridula L.) are one of the major pests in many soybean producing areas. They cause a decrease in yield and affect seed quality by reducing viability and vigor. Alterations have been reported in the oxidative response and in the secondary metabolites in different plant species due to insect damage. However, there is little information available on soybean-stink bug interactions. In this study we compare the response of undamaged and damaged seeds by Nezara viridula in two soybean cultivars, IAC-100 (resistant) and Davis (susceptible), grown under greenhouse conditions. Pod hardness, H₂O₂ generation, enzyme activities in guaiacol peroxidase (GPOX), catalase (CAT) and superoxide dismutase (SOD) as well as lipoxygenase expression and isoflavonoid production were quantified. Our results showed a greater resistance of IAC-100 to pod penetration, a decrease in peroxide content after stink bug attack, and higher GPOX, CAT and SOD activities in seeds due to the genotype and to the genotype-interaction with the herbivory treatment. Induction of LOX expression in both cultivars and higher production of isoflavonoids in IAC-100 were also detected. It was then concluded that the herbivory stink bug induces pathways related to oxidative stress and to the secondary metabolites in developing seeds of soybean and that differences between cultivars hold promise for a plant breeding program.

Genome scan detection of selective sweeps among biotypes of the soybean aphid, Aphis glycines, with differing virulence to resistance to A. glycines (Rag) traits in soybean, Glycine max

Multiple biotypes of soybean aphid, Aphis glycines, occur in North America adapted for survival (virulence) on soybean, Glycine max, with one or more different resistance to A. glycines (Rag) traits. The degree of genome-wide variance between biotypes and the basis of virulence remains unknown, but the latter is hypothesized to involve secreted effector proteins. Between 167,249 and 217,750 single nucleotide polymorphisms (SNPs) were predicted from whole genome re-sequencing of A. glycines avirulent biotype 1 (B1) and virulent B2, B3 and B4 colony-derived iso-female lines when compared to the draft B1 genome assembly, Ag_bt1_v6.0. Differences in nucleotide diversity indices (π) estimated within 1000 bp sliding windows demonstrated that 226 of 353 (64.0%) regions most differentiated between B1 and ≥ 2 virulent biotypes, representing < 0.1% of the 308 Mb assembled genome size, are located on 15 unordered scaffolds. Furthermore, these 226 intervals were coincident and show a significant association with 326 of 508 SNPs with significant locus-by-locus FST estimates between biotype populations (r = 0.6271; F1,70 = 45.36, P < 0.001) and genes showing evidence of directions selection (πN/πS > 2.0; r = 0.6233; F1,70 = 50.20, P < 0.001). A putative secreted effector glycoprotein is encoded in proximity to genome intervals of low estimated π (putative selective sweep) within avirulent B1 compared to all three virulent biotypes. Additionally, SNPs are clustered in or in proximity to genes putatively involved in intracellular protein cargo transport and the regulation of secretion. Results of this study indicate that factors on a small number of scaffolds of the A. glycines genome may contribute to variance in virulence towards Rag traits in G. max.

The Resistant Soybean-Aphis glycines Interaction: Current Knowledge and Prospects

Soybean aphids (Aphis glycines Matsumura) are invasive insect pests of soybean, and they cause significant yield losses. Resistance to soybean aphids is conferred by Resistance to Aphis glycines (Rag) genes. Since the first discovery of aphid-resistant soybean genotypes in 2004, several studies have attempted to characterize Rag genes from aphid-resistant soybean genotypes. To date, 12 Rag genes and four quantitative trait loci for aphid resistance have been reported on soybean chromosomes 07, 08, 13, 16, and 17. Although candidate genes have been proposed for several discovered Rag loci, additional studies are needed to pinpoint, validate, and further explain the potential mechanisms of Rag gene action. A major challenge to utilizing host plant resistance is the discovery of virulent aphid biotypes that can colonize aphid-resistant soybean. This occurrence suggests the need for additional studies to devise strategies to enhance the effectiveness of aphid-resistant soybean. In this mini review, we discuss current knowledge on the resistant soybean-Aphis glycines interaction, potential mechanisms of Rag gene action, opportunities to discover new Rag genes, and prospects for utilization of host plant resistance to manage soybean aphids. A clearer understanding of host plant resistance to soybean aphids will guide researchers on strategies for developing soybean varieties with more durable aphid resistance, reducing the present challenge of virulent aphid biotypes.

Characterization and application of a gall midge resistance gene (Gm6) from Oryza sativa 'Kangwenqingzhan'

The resistance gene Gm6 was mapped and characterized using near-isogenic and pyramided lines, followed by marker-assisted selection to develop lines with resistance to both gall midge and brown planthopper. The Asian rice gall midge (GM; Orseolia oryzae; Diptera: Cecidomyiidae) is a major destructive pest affecting rice cultivation regions. The characterization of GM-resistance genes and the breeding of resistant varieties are together considered the most efficient strategy for managing this insect. Here, the Gm6 resistance gene derived from the Kangwenqingzhan (KW) variety was found to be located on the long arm of chromosome 4 using the F₂ population of 9311/KW. The region was narrowed to a 90-kb segment flanked by the markers YW91 and YW3-4 using backcrossing populations. Based on no-choice feeding and host choice tests, GM development and growth in near-isogenic lines (NILs) were severely restricted compared to that in the 9311 control. On day 8, the average GM body length was 0.69 mm and 0.56 mm on NILs and 9311, respectively, and the differences were more significant at later time points. However, GM insects exhibited no host preference between NILs and 9311, and there was normal egg hatching on the resistant plants. We developed pyramided lines carrying BPH27, BPH36, and Gm6 by crossing and backcrossing with marker-assisted selection. These lines were similar to the KW parent in terms of agronomic traits while also exhibiting high resistance to brown planthopper (BPH) and GM. The present mapping and characterization of Gm6 will facilitate map-based cloning of this important resistance gene and its application in the breeding of insect-resistant rice varieties.

Plant phenolics mediated bottom-up effects of elevated CO2 on Acyrthosiphon pisum and its parasitoid Aphidius avenae

Elevated concentrations of atmospheric CO₂ can alter plant secondary metabolites, which play important roles in the interactions among plants, herbivorous insects and natural enemies. However, few studies have examined the cascading effects of host plant secondary metabolites on tri-trophic interactions under elevated CO₂ (eCO₂ ). In this study, we determined the effects of eCO₂ on the growth and foliar phenolics of Medicago truncatula and the cascading effects on two color genotypes of Acyrthosiphon pisum (pink vs. green) and their parasitoid Aphidius avenae in the field open-top chambers. Our results showed that eCO₂ increased photosynthetic rate, nodule number, yield and the total phenolic content of M. truncatula. eCO₂ had contrasting effects on two genotypes of A. pisum; the green genotype demonstrated increased population abundance, fecundity, growth and feeding efficiency, while the pink genotype showed decreased fitness and these were closely associated with the foliar genstein content. Furthermore, eCO₂ decreased the parasitic rate of A. avenae independent of aphid genotypes. eCO₂ prolonged the emergence time and reduced the emergence rate and percentage of females when associated with the green genotype, but little difference, except for increased percentage of females, was observed in A. avenae under eCO₂ when associated with the pink genotype, indicating that parasitoids can perceive and discriminate the qualities of aphid hosts. We concluded that eCO₂ altered plant phenolics and thus the performance of aphids and parasitoids. Our results indicate that plant phenolics vary by different abiotic and biotic stimuli and could potentially deliver the cascading effects of eCO₂ to the higher trophic levels. Our results also suggest that the green genotype is expected to perform better in future eCO₂ because of decreased plant resistance after its infestation and decreased parasitic rate.

BABA application improves soybean resistance to aphid through activation of phenylpropanoid metabolism and callose deposition

BACKGROUND

Beta-aminobutyric acid (BABA) confer plant resistance to a broad spectrum of biotic and abiotic stresses. The soybean aphid (SBA), is native to eastern Asia and is a predominant insect pest of soybean. Both isoflavone and lignin pathway are important branches of the general phenylpropanoid pathway, which would be likely associated with resistance against soybean aphid. However, little is known about the role of the phenylpropanoid pathway in defense response to SBA as induced by BABA application.

RESULTS

The application of BABA effectively enhanced soybean resistance against Aphis glycines, the soybean aphid. Consistent with significantly increased content of isoflavones, especially genistein, the related biosynthetic genes were upregulated by use of BABA. Lignin, another important defense component against arthropods, accumulated at a high level and four lignin biosynthesis related genes were also activated. Additionally, BABA application augmented the expression of callose synthase genes and increased callose deposition in SBA-infested seedlings. In non-caged and caged tests, SBA numbers were significantly reduced in BABA-treated seedlings.

CONCLUSION

These results demonstrate that application of BABA has an obvious positive effect on soybean resistance to aphids, and this defense response partly depends on the potentiation of isoflavone biosynthesis and callose deposition. © 2019 Society of Chemical Industry.

Increases in Genistein in Medicago sativa Confer Resistance against the Pisum Host Race of Acyrthosiphon pisum

Interspecific interaction with host plants have important consequences for the host race formation of herbivorous insects. Plant secondary metabolites, particularly those that are involved in host races specializing on plants, warrant the theory of host specialization. Acyrthosiphon pisum comprises various host races that adapt to different Fabaceae plants, which provides an ideal system for determining the behavioral and physiological mechanisms underlying host-adaptive diversification. The current study evaluated the effects of host transfer on population fitness, feeding behavior and the transcriptome-wide gene expression of the two host races of A. pisum, one of which was originally from Medicago sativa and the other from Pisum sativum. The results showed that the Pisum host race of A. pisum had a lower population abundance and feeding efficiency than the Medicago host race in terms of a longer penetration time and shorter duration times of phloem ingestion when fed on M. sativa. In contrast, few differences were found in the population abundance and feeding behavior of A. pisum between the two host races when fed on P. sativum. Meanwhile, of the nine candidate phenolic compounds, only genistein was significantly affected by aphid infestation; higher levels of genistein were detected in M. sativa after feeding by the Pisum host race, but these levels were reduced relative to uninfested controls after feeding by the Medicago host race, which suggested that genistein may be involved in the specialization of the aphid host race on M. sativa. Further exogenous application of genistein in artificial diets showed that the increase in genistein reduced the survival rate of the Pisum host race but had little effect on that of the Medicago host race. The transcriptomic profiles indicated that the transcripts of six genes with functions related to detoxification were up-regulated in the Pisum host race relative to the Medicago host race of A. pisum. These results suggested that the inducible plant phenolics and associated metabolic process in aphids resulted in their differential adaptations to their Fabaceae host.

Transcriptional and Chemical Changes in Soybean Leaves in Response to Long-Term Aphid Colonization

Soybean aphids (Aphis glycines Matsumura) are specialized insects that feed on soybean (Glycine max) phloem sap. Transcriptome analyses have shown that resistant soybean plants mount a fast response that limits aphid feeding and population growth. Conversely, defense responses in susceptible plants are slower and it is hypothesized that aphids block effective defenses in the compatible interaction. Unlike other pests, aphids can colonize plants for long periods of time; yet the effect on the plant transcriptome after long-term aphid feeding has not been analyzed for any plant-aphid interaction. We analyzed the susceptible and resistant (Rag1) transcriptome response to aphid feeding in soybean plants colonized by aphids (biotype 1) for 21 days. We found a reduced resistant response and a low level of aphid growth on Rag1 plants, while susceptible plants showed a strong response consistent with pattern-triggered immunity. GO-term analyses identified chitin regulation as one of the most overrepresented classes of genes, suggesting that chitin could be one of the hemipteran-associated molecular pattern that triggers this defense response. Transcriptome analyses also indicated the phenylpropanoid pathway, specifically isoflavonoid biosynthesis, was induced in susceptible plants in response to long-term aphid feeding. Metabolite analyses corroborated this finding. Aphid-treated susceptible plants accumulated daidzein, formononetin, and genistein, although glyceollins were present at low levels in these plants. Choice experiments indicated that daidzein may have a deterrent effect on aphid feeding. Mass spectrometry imaging showed these isoflavones accumulate likely in the mesophyll cells or epidermis and are absent from the vasculature, suggesting that isoflavones are part of a non-phloem defense response that can reduce aphid feeding. While it is likely that aphid can initially block defense responses in compatible interactions, it appears that susceptible soybean plants can eventually mount an effective defense in response to long-term soybean aphid colonization.

Elevated CO2 and O3 alter the feeding efficiency of Acyrthosiphon pisum and Aphis craccivora via changes in foliar secondary metabolites

Elevated CO2 and O3 can affect aphid performance via altering plant nutrients, however, little is known about the role of plant secondary metabolites in this process, especially for aphids feeding behaviors. We determined the effects of elevated CO2 and O3 on the growth and phenolics of alfalfa (Medicago sativa) and feeding behaviors of the pea aphids (Acyrthosiphon pisum) and cowpea aphids (Aphis craccivora). Elevated CO2 improved plant growth, but could not completely offset the negative effects of elevated O3. Elevated O3 increased foliar genistin content at the vegetative stage, increased ferulic acid at the reproductive stage, and elevated CO2 increased those at both stages. Simultaneously elevated CO2 and O3 increased foliar ferulic acid content at the reproductive stage and increased genistin content at both stages. For pea aphids, feeding efficiency was reduced under elevated CO2 at the reproductive stage and decreased under elevated O3 at the vegetative stage. For cowpea aphids, feeding efficiency was increased under elevated CO2 at the vegetative stage and decreased under elevated O3 at both stages. Simultaneously elevated CO2 and O3 decreased both aphids feeding efficiency. We concluded that CO2 and O3 independently or interactively had different effects on two aphids feeding behaviors through altering foliar ferulic acid and genistin contents.

Fine mapping of the soybean aphid-resistance genes Rag6 and Rag3c from Glycine soja 85-32

KEY MESSAGE

Rag6 and Rag3c were delimited to a 49-kb interval on chromosome 8 and a 150-kb interval on chromosome 16, respectively. Structural variants in the exons of candidate genes were identified. The soybean aphid, an invasive species, has significantly threatened soybean production in North America since 2000. Host-plant resistance is known as an ideal management strategy for aphids. Two novel aphid-resistance loci, Rag6 and Rag3c, from Glycine soja 85-32, were previously detected in a 10.5-cM interval on chromosome 8 and a 7.5-cM interval on chromosome 16, respectively. Defining the exact genomic position of these two genes is critical for improving the effectiveness of marker-assisted selection for aphid resistance and for identification of the functional genes. To pinpoint the locations of Rag6 and Rag3c, four populations segregating for Rag6 and Rag3c were used to fine map these two genes. The availability of the Illumina Infinium SoySNP50K/8K iSelect BeadChip, combined with single-nucleotide polymorphism (SNP) markers discovered through the whole-genome re-sequencing of E12901, facilitated the fine mapping process. Rag6 was refined to a 49-kb interval on chromosome 8 with four candidate genes, including three clustered nucleotide-binding site leucine-rich repeat (NBS-LRR) genes and an amine oxidase encoding gene. Rag3c was refined to a 150-kb interval on chromosome 16 with 11 candidate genes, two of which are a LRR gene and a lipase gene. Moreover, by sequencing the whole-genome exome-capture of the resistant source (E12901), structural variants were identified in the exons of the candidate genes of Rag6 and Rag3c. The closely linked SNP markers and the candidate gene information presented in this study will be significant resources for integrating Rag6 and Rag3c into elite cultivars and for future functional genetics studies.

Mapping novel aphid resistance QTL from wild soybean, Glycine soja 85-32

KEY MESSAGE

Two novel QTLs conferring aphid resistance were mapped and validated on soybean chromosomes 8 and 16, respectively. Closely linked markers were developed to assist breeding for aphid resistance. Soybean aphid, Aphis glycines Matsumura, is a highly destructive pest for soybean production. E08934, a soybean advanced breeding line derived from the wild soybean Glycine soja 85-32, has shown strong resistance to aphids. To dissect the genetic basis of aphid resistance in E08934, a mapping population (070020) consisting of 140 F3-derived lines was developed by crossing E08934 with an aphid-susceptible line E00003. This mapping population was evaluated for aphid resistance in a greenhouse trial in 2010 and three field trials in 2009, 2010, and 2011, respectively. The broad-sense heritability across the field trials was 0.84. In the mapping population 070020, two major quantitative trait loci (QTL) were detected as significantly associated with aphid resistance, and designated as Rag6 and Rag3c, respectively. Rag6 was mapped to a 10.5 centiMorgan (cM) interval between markers MSUSNP08-2 and Satt209 on chromosome 8, explaining 19.5-46.4% of the phenotypic variance in different trials. Rag3c was located at a 7.5 cM interval between markers MSUSNP16-10 and Sat_370 on chromosome 16, explaining 12.5-22.9% of the phenotypic variance in different trials. Rag3c had less resistance effect than Rag6 across all the trials. Furthermore, Rag6 and Rag3c were confirmed in two validation populations with different genetic backgrounds. No significant interaction was detected between Rag6 and Rag3c in either the mapping population or the validation populations. Both Rag6 and Rag3c were indicated as conferring antibiosis resistance to aphids by a no-choice test. The new aphid-resistance gene(s) derived from the wild germplasm G. soja 85-32 are valuable in improving soybeans for aphid resistance.

Exploiting Phenylpropanoid Derivatives to Enhance the Nutraceutical Values of Cereals and Legumes

Phenylpropanoids are a diverse chemical class with immense health benefits that are biosynthesized from the aromatic amino acid L-phenylalanine. This article reviews the progress for accessing variation in phenylpropanoids in germplasm collections, the genetic and molecular basis of phenylpropanoid biosynthesis, and the development of cultivars dense in seed-phenylpropanoids. Progress is also reviewed on high-throughput assays, factors that influence phenylpropanoids, the site of phenylpropanoids accumulation in seed, Genotype × Environment interactions, and on consumer attitudes for the acceptance of staple foods rich in phenylpropanoids. A paradigm shift was noted in barley, maize, rice, sorghum, soybean, and wheat, wherein cultivars rich in phenylpropanoids are grown in Europe and North and Central America. Studies have highlighted some biological constraints that need to be addressed for development of high-yielding cultivars that are rich in phenylpropanoids. Genomics-assisted breeding is expected to facilitate rapid introgression into improved genetic backgrounds by minimizing linkage drag. More research is needed to systematically characterize germplasm pools for assessing variation to support crop genetic enhancement, and assess consumer attitudes to foods rich in phenylpropanoids.

Investigation of the Chemical Interface in the Soybean-Aphid and Rice-Bacteria Interactions Using MALDI-Mass Spectrometry Imaging

Mass spectrometry imaging (MSI) is an emerging technology for high-resolution plant biology. It has been utilized to study plant-pest interactions, but limited to the surface interfaces. Here we expand the technology to explore the chemical interactions occurring inside the plant tissues. Two sample preparation methods, imprinting and fracturing, were developed and applied, for the first time, to visualize internal metabolites of leaves in matrix-assisted laser desorption ionization (MALDI)-MSI. This is also the first time nanoparticle-based ionization was implemented to ionize diterpenoid phytochemicals that were difficult to analyze with traditional organic matrices. The interactions between rice-bacterium and soybean-aphid were investigated as two model systems to demonstrate the capability of high-resolution MSI based on MALDI. Localized molecular information on various plant- or pest-derived chemicals provided valuable insight for the molecular processes occurring during the plant-pest interactions. Specifically, salicylic acid and isoflavone based resistance was visualized in the soybean-aphid system and antibiotic diterpenoids in rice-bacterium interactions.

Mapping isoflavone QTL with main, epistatic and QTL × environment effects in recombinant inbred lines of soybean

Soybean (Glycine max (L.) Merr.) isoflavone is important for human health and plant defense system. To identify novel quantitative trait loci (QTL) and epistatic QTL underlying isoflavone content in soybean, F5:6, F5:7 and F5:8 populations of 130 recombinant inbred (RI) lines, derived from the cross of soybean cultivar 'Zhong Dou 27' (high isoflavone) and 'Jiu Nong 20' (low isoflavone), were analyzed with 95 new SSR markers. A new linkage map including 194 SSR markers and covering 2,312 cM with mean distance of about 12 cM between markers was constructed. Thirty four QTL for both individual and total seed isoflavone contents of soybean were identified. Six, seven, ten and eleven QTL were associated with daidzein (DZ), glycitein (GC), genistein (GT) and total isoflavone (TI), respectively. Of them 23 QTL were newly identified. The qTIF_1 between Satt423 and Satt569 shared the same marker Satt569 with qDZF_2, qGTF_1 and qTIF_2. The qGTD2_1 between Satt186 and Satt226 was detected in four environments and explained 3.41%-10.98% of the phenotypic variation. The qGTA2_1, overlapped with qGCA2_1 and detected in four environments, was close to the previously identified major QTL for GT, which were responsible for large a effects. QTL (qDZF_2, qGTF_1 and qTIF_2) between Satt144-Satt569 were either clustered or pleiotropic. The qGCM_1, qGTM_1 and qTIM_1 between Satt540-Sat_244 explained 2.02%-9.12% of the phenotypic variation over six environments. Moreover, the qGCE_1 overlapped with qGTE_1 and qTIE_1, the qTIH_2 overlapped with qGTH_1, qGCI_1 overlapped with qDZI_1, qTIL_1 overlapped with qGTL_1, and qTIO_1 overlapped with qGTO_1. In this study, some of unstable QTL were detected in different environments, which were due to weak expression of QTL, QTL by environment interaction in the opposite direction to a effects, and/or epistasis. The markers identified in multi-environments in this study could be applied in the selection of soybean cultivars for higher isoflavone content and in the map-based gene cloning.

Expression quantitative trait loci infer the regulation of isoflavone accumulation in soybean (Glycine max L. Merr.) seed

BACKGROUND

Mapping expression quantitative trait loci (eQTL) of targeted genes represents a powerful and widely adopted approach to identify putative regulatory variants. Linking regulation differences to specific genes might assist in the identification of networks and interactions. The objective of this study is to identify eQTL underlying expression of four gene families encoding isoflavone synthetic enzymes involved in the phenylpropanoid pathway, which are phenylalanine ammonia-lyase (PAL; EC 4.3.1.5), chalcone synthase (CHS; EC 2.3.1.74), 2-hydroxyisoflavanone synthase (IFS; EC1.14.13.136) and flavanone 3-hydroxylase (F3H; EC 1.14.11.9). A population of 130 recombinant inbred lines (F5:11), derived from a cross between soybean cultivar 'Zhongdou 27' (high isoflavone) and 'Jiunong 20' (low isoflavone), and a total of 194 simple sequence repeat (SSR) markers were used in this study. Overlapped loci of eQTLs and phenotypic QTLs (pQTLs) were analyzed to identify the potential candidate genes underlying the accumulation of isoflavone in soybean seed.

RESULTS

Thirty three eQTLs (thirteen cis-eQTLs and twenty trans-eQTLs) underlying the transcript abundance of the four gene families were identified on fifteen chromosomes. The eQTLs between Satt278-Sat_134, Sat_134-Sct_010 and Satt149-Sat_234 underlie the expression of both IFS and CHS genes. Five eQTL intervals were overlapped with pQTLs. A total of eleven candidate genes within the overlapped eQTL and pQTL were identified.

CONCLUSIONS

These results will be useful for the development of marker-assisted selection to breed soybean cultivars with high or low isoflavone contents and for map-based cloning of new isoflavone related genes.

Quantitative trait loci analysis of individual and total isoflavone contents in soybean seeds

Soybean isoflavones play diverse roles in human health, including cancers, osteoporosis, heart disease, menopausal symptoms and pabulums. The objective of this study was to identify the quantitative trait loci (QTL) associated with the isoflavones daidzein (DC), genistein (GeC), glycitein (GlC) and total isoflavone contents (TIC) in soybean seeds. A population of 184 F2:10 recombinant inbred lines derived from a 'Xiaoheidou' x 'GR8836' cross was planted in pot and field conditions to evaluate soybean isoflavones. Twenty-one QTL were detected by composite interval mapping. Several QTL were associated with the traits for DC, GeC, GlC and TIC only. QDGeGlTIC4_1 and QDGlTIC12_1 are reported first in this study and were associated with the DC, GeC, GlC and TIC traits simultaneously. The QTL identified have potential value for marker-assisted selection to develop soybean varieties with desirable isoflavone content.

DL-β-aminobutyric acid-induced resistance in soybean against Aphis glycines Matsumura (Hemiptera: Aphididae)

Priming can improve plant innate capability to deal with the stresses caused by both biotic and abiotic factors. In this study, the effect of DL-β-amino-n-butyric acid (BABA) against Aphis glycines Matsumura, the soybean aphid (SA) was evaluated. We found that 25 mM BABA as a root drench had minimal adverse impact on plant growth and also efficiently protected soybean from SA infestation. In both choice and non-choice tests, SA number was significantly decreased to a low level in soybean seedlings drenched with 25 mM BABA compared to the control counterparts. BABA treatment resulted in a significant increase in the activities of several defense enzymes, such as phenylalanine ammonia-lyase (PAL), peroxidase (POX), polyphenol oxidase (PPO), chitinase (CHI), and β-1, 3-glucanase (GLU) in soybean seedlings attacked by aphid. Meanwhile, the induction of 15 defense-related genes by aphid, such as AOS, CHS, MMP2, NPR1-1, NPR1-2, and PR genes, were significantly augmented in BABA-treated soybean seedlings. Our study suggest that BABA application is a promising way to enhance soybean resistance against SA.

Characterization and genetic dissection of resistance to spotted alfalfa aphid (Therioaphis trifolii) in Medicago truncatula

Aphids cause significant yield losses in agricultural crops worldwide. Medicago truncatula, a model legume, cultivated pasture species in Australia and close relative of alfalfa (Medicago sativa), was used to study the defence response against Therioaphis trifolii f. maculate [spotted alfalfa aphid (SAA)]. Aphid performance and plant damage were compared among three accessions. A20 is highly susceptible, A17 has moderate resistance, and Jester is strongly resistant. Subsequent analyses using A17 and A20, reciprocal F1s and an A17×A20 recombinant inbred line (RIL) population revealed that this moderate resistance is phloem mediated and involves antibiosis and tolerance but not antixenosis. Electrical penetration graph analysis also identified a novel waveform termed extended potential drop, which occurred following SAA infestation of M. truncatula. Genetic dissection using the RIL population revealed three quantitative trait loci on chromosomes 3, 6, and 7 involved in distinct modes of aphid defence including antibiosis and tolerance. An antibiosis locus resides on linkage group 3 (LG3) and is derived from A17, whereas a plant tolerance and antibiosis locus resides on LG6 and is derived from A20, which exhibits strong temporary tolerance. The loci identified reside in regions harbouring classical resistance genes, and introgression of these loci in current medic cultivars may help provide durable resistance to SAA, while elucidation of their molecular mechanisms may provide valuable insight into other aphid-plant interactions.

Genetic mapping of a novel gene for soybean aphid resistance in soybean (Glycine max [L.] Merr.) line P203 from China

The soybean aphid (SA: Aphis glycines Matsumura) is a worldwide pest of soybean (Glycine max [L.] Merr.). The objectives of this study were to identify the type of aphid resistance and the resistance phenotype in soybean line 'P203', and to map the relative position of the gene involved. Compared with cultivars 'P746' and 'Dongnong 47', P203 was demonstrated to possess antixenosis resistance. P203 prevented aphids from reproducing in a choice test, but the resistance level decreased significantly in a no-choice test at 11 and 21 days after infestation. Analysis of 273 Dongnong 47/P203 F2 plants and confirmed using 260 F2:3 families revealed that a single dominant gene from P203 was positioned between marker loci Sat_377 and Satt409 on chromosome 8. The gene was further mapped to a 1.57 Mb interval flanked by marker loci BARCSOYSSR_08_1451 and BARCSOYSSR_08_1527. We developed five new SSR markers in the target interval and the resistance locus mapped between new markers SSR_08_75 and SSR_08_88 with the genetic distance of 1.1 and 1.0 cM corresponding to a physical distance of 192 kb on the Williams 82 8X draft genome assembly (Glyma1.01). A single serine/threonine protein kinase gene is present in this region, suggesting that the SA resistance mechanism in P203 may be different from those previously reported. Therefore, the resistance gene could very well be novel, and could be valuable in soybean aphid resistance breeding programs.