A Summary of Two Decades of QTL and Candidate Genes That Control Seed Tocopherol Contents in Maize (Zea mays L.)

Tocopherols are secondary metabolites synthesized through the shikimate biosynthetic pathway in the plastids of most plants. It is well known that α-Tocopherol (vitamin E) has many health benefits for humans and animals; therefore, it is highly used in human and animal diets. Tocopherols vary considerably in most crop (and plant) species and within cultivars of the same species depending on environmental and growth conditions; tocopherol content is a polygenic, complex traits, and its inheritance is poorly understood. The objective of this review paper was to summarize all identified quantitative trait loci (QTL) that control seed tocopherols and related contents identified in maize (Zea mays) during the past two decades (2002-2022). Candidate genes identified within these QTL regions are also discussed. The QTL described here, and candidate genes identified within these genomic regions could be used in breeding programs to develop maize cultivars with high, beneficial levels of seed tocopherol contents.

Deep Learning Model for Classifying and Evaluating Soybean Leaf Disease Damage

Soybean (Glycine max (L.) Merr.) is a major source of oil and protein for human food and animal feed; however, soybean crops face diverse factors causing damage, including pathogen infections, environmental shifts, poor fertilization, and incorrect pesticide use, leading to reduced yields. Identifying the level of leaf damage aids yield projections, pesticide, and fertilizer decisions. Deep learning models (DLMs) and neural networks mastering tasks from abundant data have been used for binary healthy/unhealthy leaf classification. However, no DLM predicts and categorizes soybean leaf damage severity (five levels) for tailored pesticide use and yield forecasts. This paper introduces a novel DLM for accurate damage prediction and classification, trained on 2930 near-field soybean leaf images. The model quantifies damage severity, distinguishing healthy/unhealthy leaves and offering a comprehensive solution. Performance metrics include accuracy, precision, recall, and F1-score. This research presents a robust DLM for soybean damage assessment, supporting informed agricultural decisions based on specific damage levels and enhancing crop management and productivity.

No Pairwise Interactions of GmSNAP18, GmSHMT08 and AtPR1 with Suppressed AtPR1 Expression Enhance the Susceptibility of Arabidopsis to Beet Cyst Nematode

GmSNAP18 and GmSHMT08 are two major genes conferring soybean cyst nematode (SCN) resistance in soybean. Overexpression of either of these two soybean genes would enhance the susceptibility of Arabidopsis to beet cyst nematode (BCN), while overexpression of either of their corresponding orthologs in Arabidopsis, AtSNAP2 and AtSHMT4, would suppress it. However, the mechanism by which these two pairs of orthologous genes boost or inhibit BCN susceptibility of Arabidopsis still remains elusive. In this study, Arabidopsis with simultaneously overexpressed GmSNAP18 and GmSHMT0 suppressed the growth of underground as well as above-ground parts of plants. Furthermore, Arabidopsis that simultaneously overexpressed GmSNAP18 and GmSHMT08 substantially stimulated BCN susceptibility and remarkably suppressed expression of AtPR1 in the salicylic acid signaling pathway. However, simultaneous overexpression of GmSNAP18 and GmSHMT08 did not impact the expression of AtJAR1 and AtHEL1 in the jasmonic acid and ethylene signaling pathways. GmSNAP18, GmSHMT08, and a pathogenesis-related (PR) protein, GmPR08-Bet VI, in soybean, and AtSNAP2, AtSHMT4, and AtPR1 in Arabidopsis could interact pair-wisely for mediating SCN and BCN resistance in soybean and Arabidopsis, respectively. Both AtSNAP2 and AtPR1 were localized on the plasma membrane, and AtSHMT4 was localized both on the plasma membrane and in the nucleus of cells. Nevertheless, after interactions, AtSNAP2 and AtPR1 could partially translocate into the cell nucleus. GmSNAP18 interacted with AtSHMT4, and GmSHMT4 interacted with AtSNAP2. However, neither GmSNAP18 nor GmSHMT08 interacted with AtPR1. Thus, no pairwise interactions among α-SNAPs, SHMTs, and AtPR1 occurred in Arabidopsis overexpressing either GmSNAP18 or GmSHMT08, or both of them. Transgenic Arabidopsis overexpressing either GmSNAP18 or GmSHMT08 substantially suppressed AtPR1 expression, while transgenic Arabidopsis overexpressing either AtSNAP2 or AtSHMT4 remarkably enhanced it. Taken together, no pairwise interactions of GmSNAP18, GmSHMT08, and AtPR1 with suppressed expression of AtPR1 enhanced BCN susceptibility in Arabidopsis. This study may provide a clue that nematode-resistant or -susceptible functions of plant genes likely depend on both hosts and nematode species.

Genetic Mapping for QTL Associated with Seed Nickel and Molybdenum Accumulation in the Soybean 'Forrest' by 'Williams 82' RIL Population

Understanding the genetic basis of seed Ni and Mo is essential. Since soybean is a major crop in the world and a major source for nutrients, including Ni and Mo, the objective of the current research was to map genetic regions (quantitative trait loci, QTL) linked to Ni and Mo concentrations in soybean seed. A recombinant inbred line (RIL) population was derived from a cross between 'Forrest' and 'Williams 82' (F × W82). A total of 306 lines was used for genotyping using 5405 single nucleotides polymorphism (SNP) markers using Infinium SNP6K BeadChips. A two-year experiment was conducted and included the parents and the RIL population. One experiment was conducted in 2018 in North Carolina (NC), and the second experiment was conducted in Illinois in 2020 (IL). Logarithm of the odds (LOD) of ≥2.5 was set as a threshold to report identified QTL using the composite interval mapping (CIM) method. A wide range of Ni and Mo concentrations among RILs was observed. A total of four QTL (qNi-01, qNi-02, and qNi-03 on Chr 2, 8, and 9, respectively, in 2018, and qNi-01 on Chr 20 in 2020) was identified for seed Ni. All these QTL were significantly (LOD threshold > 2.5) associated with seed Ni, with LOD scores ranging between 2.71-3.44, and with phenotypic variance ranging from 4.48-6.97%. A total of three QTL for Mo (qMo-01, qMo-02, and qMo-03 on Chr 1, 3, 17, respectively) was identified in 2018, and four QTL (qMo-01, qMo-02, qMo-03, and qMo-04, on Chr 5, 11, 14, and 16, respectively) were identified in 2020. Some of the current QTL had high LOD and significantly contributed to the phenotypic variance for the trait. For example, in 2018, Mo QTL qMo-01 on Chr 1 had LOD of 7.8, explaining a phenotypic variance of 41.17%, and qMo-03 on Chr 17 had LOD of 5.33, with phenotypic variance explained of 41.49%. In addition, one Mo QTL (qMo-03 on Chr 14) had LOD of 9.77, explaining 51.57% of phenotypic variance related to the trait, and another Mo QTL (qMo-04 on Chr 16) had LOD of 7.62 and explained 49.95% of phenotypic variance. None of the QTL identified here were identified twice across locations/years. Based on a search of the available literature and of SoyBase, the four QTL for Ni, identified on Chr 2, 8, 9, and 20, and the five QTL associated with Mo, identified on Chr 1, 17, 11, 14, and 16, are novel and not previously reported. This research contributes new insights into the genetic mapping of Ni and Mo, and provides valuable QTL and molecular markers that can potentially assist in selecting Ni and Mo levels in soybean seeds.

Quantitative Trait Loci and Candidate Genes That Control Seed Sugars Contents in the Soybean 'Forrest' by 'Williams 82' Recombinant Inbred Line Population

Soybean seed sugars are among the most abundant beneficial compounds for human and animal consumption in soybean seeds. Higher seed sugars such as sucrose are desirable as they contribute to taste and flavor in soy-based food. Therefore, the objectives of this study were to use the 'Forrest' by 'Williams 82' (F × W82) recombinant inbred line (RIL) soybean population (n = 309) to identify quantitative trait loci (QTLs) and candidate genes that control seed sugar (sucrose, stachyose, and raffinose) contents in two environments (North Carolina and Illinois) over two years (2018 and 2020). A total of 26 QTLs that control seed sugar contents were identified and mapped on 16 soybean chromosomes (chrs.). Interestingly, five QTL regions were identified in both locations, Illinois and North Carolina, in this study on chrs. 2, 5, 13, 17, and 20. Amongst 57 candidate genes identified in this study, 16 were located within 10 Megabase (MB) of the identified QTLs. Amongst them, a cluster of four genes involved in the sugars' pathway was collocated within 6 MB of two QTLs that were detected in this study on chr. 17. Further functional validation of the identified genes could be beneficial in breeding programs to produce soybean lines with high beneficial sucrose and low raffinose family oligosaccharides.

Soybean gene co-expression network analysis identifies two co-regulated gene modules associated with nodule formation and development

Gene co-expression network analysis is an efficient systems biology approach for the discovery of novel gene functions and trait-associated gene modules. To identify clusters of functionally related genes involved in soybean nodule formation and development, we performed a weighted gene co-expression network analysis. Two nodule-specific modules (NSM-1 and NSM-2, containing 304 and 203 genes, respectively) were identified. The NSM-1 gene promoters were significantly enriched in cis-binding elements for ERF, MYB, and C2H2-type zinc transcription factors, whereas NSM-2 gene promoters were enriched in cis-binding elements for TCP, bZIP, and bHLH transcription factors, suggesting a role of these regulatory factors in the transcriptional activation of nodule co-expressed genes. The co-expressed gene modules included genes with potential novel roles in nodulation, including those involved in xylem development, transmembrane transport, the ethylene signalling pathway, cytoskeleton organization, cytokinesis and regulation of the cell cycle, regulation of meristem initiation and growth, transcriptional regulation, DNA methylation, and histone modifications. Functional analysis of two co-expressed genes using TILLING mutants provided novel insight into the involvement of unsaturated fatty acid biosynthesis and folate metabolism in nodule formation and development. The identified gene co-expression modules provide valuable resources for further functional genomics studies to dissect the genetic basis of nodule formation and development in soybean.

Soybean ZINC FINGER PROTEIN03 targets two SUPEROXIDE DISMUTASE1s and confers resistance to Phytophthora sojae

Phytophthora sojae causes Phytophthora root and stem rot disease of soybean (Glycine max), leading to huge annual yield loss worldwide, but resistance to Phytophthora sojae (Rps) genes remains elusive. Soybean cultivar "Yudou 29" is resistant to P. sojae strain PsMC1, and this study aimed to clone, identify, and characterize the Rps gene in Yudou 29 (RpsYD29) and clarify its functional mechanism. We map-based cloned RpsYD29 (ZINC FINGER PROTEIN03, GmZFP03) using the families of a cross between Yudou 29 and a P. sojae-susceptible soybean cultivar "Jikedou 2". P. sojae resistance of GmZFP03 was functionally validated by stable soybean genetic transformation and allele-phenotype association analysis. GmZFP03 was identified as a C2H2-type zinc finger protein transcription factor, showing 4 amino acid residue polymorphisms (V79F, G122-, G123-, and D125V) and remarkably different expression patterns between resistant and susceptible soybeans. Notably boosted activity and gene expression of superoxide dismutase (SOD) in resistant-type GmZFP03-expressed transgenic soybean, substantial enhancement of P. sojae resistance of wild-type soybean by exogenous SOD treatment, and GmZFP03 binding to and activation of 2 SOD1 (Glyma.03g242900 and Glyma.19g240400) promoters demonstrated the involvement of SOD1s in GmZFP03-mediated resistance to P. sojae strain PsMC1. Thus, this study cloned the soybean P. sojae-resistant GmZFP03, the product of which specifically targets 2 SOD1 promoters. GmZFP03 can be directly used for precise P. sojae-resistance soybean breeding.

Genome-wide identification and analysis of soybean acyl-ACP thioesterase gene family reveals the role of GmFAT to improve fatty acid composition in soybean seed

KEY MESSAGE

Soybean acyl-ACP thioesterase gene family have been characterized; GmFATA1A mutants were discovered to confer high oleic acid, while GmFATB mutants presented low palmitic and high oleic acid seed content. Soybean oil stability and quality are primarily determined by the relative proportions of saturated versus unsaturated fatty acids. Commodity soybean typically contains 11% palmitic acid, as the primary saturated fatty acids. Reducing palmitic acid content is the principal approach to minimize the levels of saturated fatty acids in soybean. Though high palmitic acid enhances oxidative stability of soybean oil, it is negatively correlated with oil and oleic acid content and can cause coronary heart diseases for humans. For plants, acyl-acyl carrier protein (ACP) thioesterases (TEs) are a group of enzymes to hydrolyze acyl group and release free fatty acid from plastid. Among them, GmFATB1A has become the main target to genetically reduce the palmitic acid content in soybean. However, the role of members in soybean acyl-ACP thioesterase gene family is largely unknown. In this study, we characterized two classes of TEs, GmFATA, and GmFATB in soybean. We also denominated two GmFATA members and discovered six additional members that belong to GmFATB gene family through phylogenetic, syntenic, and in silico analysis. Using TILLING-by-Sequencing+, we identified an allelic series of mutations in five soybean acyl-ACP thioesterase genes, including GmFATA1A, GmFATB1A, GmFATB1B, GmFATB2A, and GmFATB2B. Additionally, we discovered mutations at GmFATA1A to confer high oleic acid (up to 34.5%) content, while mutations at GmFATB presented low palmitic acid (as low as 5.6%) and high oleic acid (up to 36.5%) phenotypes. The obtained soybean mutants with altered fatty acid content can be used in soybean breeding program for improving soybean oil composition traits.

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.

Genome Wide MeDIP-Seq Profiling of Wild and Cultivated Olives Trees Suggests DNA Methylation Fingerprint on the Sensory Quality of Olive Oil

Secondary metabolites are particularly important to humans due to their pharmaceutical properties. Moreover, secondary metabolites are key compounds in climate change adaptation in long-living trees. Recently, it has been described that the domestication of Olea subspecies had no major selection signature on coding variants and was mainly related to changes in gene expression. In addition, the phenotypic plasticity in Olea subspecies was linked to the activation of transposable elements in the genes neighboring. Here, we investigated the imprint of DNA methylation in the unassigned fraction of the phenotypic plasticity of the Olea subspecies, using methylated DNA immuno-precipitation sequencing (MeDIP-seq) for a high-resolution genome-wide DNA methylation profiling of leaves and fruits during fruit development in wild and cultivated olives from Turkey. Notably, the methylation profiling showed a differential DNA methylation in secondary metabolism responsible for the sensory quality of olive oil. Here, we highlight for the first time the imprint of DNA methylation in modulating the activity of the Linoleate 9S lipoxygenase in the biosynthesis of volatile aromatic compounds. Unprecedently, the current study reveals the methylation status of the olive genome during fruit ripening.

Dissecting nematode resistance regions in soybean revealed pleiotropic effect of soybean cyst and reniform nematode resistance genes

Reniform nematode (RN, Rotylenchulus reniformis Linford & Oliveira) has emerged as one of the most important plant parasitic nematodes of soybean [Glycine max (L.) Merr.]. Planting resistant varieties is the most effective strategy for nematode management. The objective of this study was to identify quantitative trait loci (QTL) for RN resistance in an exotic soybean line, PI 438489B, using two linkage maps constructed from the Universal Soybean Linkage Panel (USLP 1.0) and next-generation whole-genome resequencing (WGRS) technology. Two QTL controlling RN resistance were identified-the soybean cyst nematode (SCN, Heterodera glycines) resistance gene GmSNAP18 at the rhg1 locus and its paralog GmSNAP11. Strong association between resistant phenotype and haplotypes of the GmSNAP11 and GmSNAP18 was observed. The results indicated that GmSNAP11 possibly could have epistatic effect on GmSNAP18, or vice versa, with the presence of a significant correlation in RN resistance of rhg1-a GmSNAP18 vs. rhg1-b GmSNAP18. Most importantly, our preliminary data suggested that GmSNAP18 and GmSNAP11 proteins physically interact in planta, suggesting that they belong to the same pathway for resistance. Unlike GmSNAP18, no indication of GmSNAP11 copy number variation was found. Moreover, gene-based single nucleotide polymorphism (SNP) markers were developed for rapid detection of RN or SCN resistance at these loci. Our analysis substantiates synergic interaction between GmSNAP11 and GmSNAP18 genes and confirms their roles in RN as well as SCN resistance. These results could contribute to a better understanding of evolution and subfunctionalization of genes conferring resistance to multiple nematode species and provide a framework for further investigations.

TILLING-by-Sequencing+ Reveals the Role of Novel Fatty Acid Desaturases (GmFAD2-2s) in Increasing Soybean Seed Oleic Acid Content

Soybean is the second largest source of oil worldwide. Developing soybean varieties with high levels of oleic acid is a primary goal of the soybean breeders and industry. Edible oils containing high level of oleic acid and low level of linoleic acid are considered with higher oxidative stability and can be used as a natural antioxidant in food stability. All developed high oleic acid soybeans carry two alleles; GmFAD2-1A and GmFAD2-1B. However, when planted in cold soil, a possible reduction in seed germination was reported when high seed oleic acid derived from GmFAD2-1 alleles were used. Besides the soybean fatty acid desaturase (GmFAD2-1) subfamily, the GmFAD2-2 subfamily is composed of five members, including GmFAD2-2A, GmFAD2-2B, GmFAD2-2C, GmFAD2-2D, and GmFAD2-2E. Segmental duplication of GmFAD2-1A/GmFAD2-1B, GmFAD2-2A/GmFAD2-2C, GmFAD2-2A/GmFAD2-2D, and GmFAD2-2D/GmFAD2-2C have occurred about 10.65, 27.04, 100.81, and 106.55 Mya, respectively. Using TILLING-by-Sequencing+ technology, we successfully identified 12, 8, 10, 9, and 19 EMS mutants at the GmFAD2-2A, GmFAD2-2B, GmFAD2-2C, GmFAD2-2D, and GmFAD2-2E genes, respectively. Functional analyses of newly identified mutants revealed unprecedented role of the five GmFAD2-2A, GmFAD2-2B, GmFAD2-2C, GmFAD2-2D, and GmFAD2-2E members in controlling the seed oleic acid content. Most importantly, unlike GmFAD2-1 members, subcellular localization revealed that members of the GmFAD2-2 subfamily showed a cytoplasmic localization, which may suggest the presence of an alternative fatty acid desaturase pathway in soybean for converting oleic acid content without substantially altering the traditional plastidial/ER fatty acid production.

TILLING-by-Sequencing+ to Decipher Oil Biosynthesis Pathway in Soybeans: A New and Effective Platform for High-Throughput Gene Functional Analysis

Reverse genetic approaches have been widely applied to study gene function in crop species; however, these techniques, including gel-based TILLING, present low efficiency to characterize genes in soybeans due to genome complexity, gene duplication, and the presence of multiple gene family members that share high homology in their DNA sequence. Chemical mutagenesis emerges as a genetically modified-free strategy to produce large-scale soybean mutants for economically important traits improvement. The current study uses an optimized high-throughput TILLING by target capture sequencing technology, or TILLING-by-Sequencing+ (TbyS+), coupled with universal bioinformatic tools to identify population-wide mutations in soybeans. Four ethyl methanesulfonate mutagenized populations (4032 mutant families) have been screened for the presence of induced mutations in targeted genes. The mutation types and effects have been characterized for a total of 138 soybean genes involved in soybean seed composition, disease resistance, and many other quality traits. To test the efficiency of TbyS+ in complex genomes, we used soybeans as a model with a focus on three desaturase gene families, GmSACPD, GmFAD2, and GmFAD3, that are involved in the soybean fatty acid biosynthesis pathway. We successfully isolated mutants from all the six gene family members. Unsurprisingly, most of the characterized mutants showed significant changes either in their stearic, oleic, or linolenic acids. By using TbyS+, we discovered novel sources of soybean oil traits, including high saturated and monosaturated fatty acids in addition to low polyunsaturated fatty acid contents. This technology provides an unprecedented platform for highly effective screening of polyploid mutant populations and functional gene analysis. The obtained soybean mutants from this study can be used in subsequent soybean breeding programs for improved oil composition traits.

Soybean TILLING-by-Sequencing+ reveals the role of novel GmSACPD members in unsaturated fatty acid biosynthesis while maintaining healthy nodules

Developing soybean lines with high levels of stearic acid is a primary goal of the soybean industry. Most high-stearic-acid soybeans carry different GmSACPD-C mutated alleles. However, due to the dual role of GmSACPD-C in seeds and nodule development, all derived deleterious GmSACPD-C mutant alleles are of extremely poor agronomic value because of defective nodulation. The soybean stearoyl-acyl carrier protein desaturase (GmSACPD) gene family is composed of five members. Comparative genomics analysis indicated that SACPD genes were duplicated and derived from a common ancestor that is still present in chlorophytic algae. Synteny analysis showed the presence of segment duplications between GmSACPD-A/GmSACPD-B, and GmSACPD-C/GmSACPD-D. GmSACPD-E was not contained in any duplicated segment and may be the result of tandem duplication. We developed a TILLING by Target Capture Sequencing (Tilling-by-Sequencing+) technology, a versatile extension of the conventional TILLING by sequencing, and successfully identified 12, 14, and 18 ethyl methanesulfonate mutants at the GmSACPD-A, GmSACPD-B, and GmSACPD-D genes, respectively. Functional analysis of all identified mutants revealed an unprecedented role of GmSACPD-A, GmSACPD-B, and GmSACPD-D in unsaturated fatty acid biosynthesis without affecting nodule development and structure. This discovery will positively impact the development of high-stearic-acid lines to enhance soybean nutritional value without potential developmental tradeoffs.

EMS-Induced Mutagenesis of Clostridium carboxidivorans for Increased Atmospheric CO2 Reduction Efficiency and Solvent Production

Clostridium carboxidivorans (P7) is one of the most important solvent-producing bacteria capable of fermenting syngas (CO, CO₂, and H₂) to produce chemical commodities when grown as an autotroph. This study aimed to develop ethyl methanesulfonate (EMS)-induced P7 mutants that were capable of growing in the presence of CO₂ as a unique source of carbon with increased solvent formation and atmospheric CO₂ reduction to limit global warming. Phenotypic analysis including growth and end product characterization of the P7 wild type (WT) demonstrated that this strain grew better at 25 °C than 37 °C when CO₂ served as the only source of carbon. In the current study, 55 mutagenized P7-EMS mutants were developed by using 100 mM and 120 mM EMS. Interestingly, using a forward genetic approach, three out of the 55 P7-EMS mutants showed a significant increase in ethanol, butyrate, and butanol production. The three P7-EMS mutants presented on average a 4.68-fold increase in concentrations of ethanol when compared to the P7-WT. Butyric acid production from 3 P7-EMS mutants contained an average of a 3.85 fold increase over the levels observed in the P7-WT cultures under the same conditions (CO₂ only). In addition, one P7-EMS mutant presented butanol production (0.23 ± 0.02 g/L), which was absent from the P7-WT under CO₂ conditions. Most of the P7-EMS mutants showed stability of the obtained end product traits after three transfers. Most importantly, the amount of reduced atmospheric CO₂ increased up to 8.72 times (0.21 g/Abs) for ethanol production and up to 8.73 times higher (0.16 g/Abs) for butyrate than the levels contained in the P7-WT. Additionally, to produce butanol, the P7-EMS_III-J mutant presented 0.082 g/Abs of CO₂ reduction. This study demonstrated the feasibility and effectiveness of employing EMS mutagenesis in generating solvent-producing anaerobic bacteria mutants with improved and novel product formation and increased atmospheric CO₂ reduction efficiency.

A pathogenesis-related protein GmPR08-Bet VI promotes a molecular interaction between the GmSHMT08 and GmSNAP18 in resistance to Heterodera glycines

Soybean cyst nematode (SCN, Heterodera glycines) is the most devastating pest affecting soybean production worldwide. SCN resistance requires both the GmSHMT08 and the GmSNAP18 in 'Peking'-type resistance. Here, we describe the molecular interaction between GmSHMT08 and GmSNAP18, which is potentiated by a pathogenesis-related protein GmPR08-Bet VI. Like GmSNAP18 and GmSHMT08, GmPR08-Bet VI expression was induced in response to SCN and its overexpression decreased SCN cysts by 65% in infected transgenic soybean roots. Overexpression of GmPR08-Bet VI did not have an effect on SCN resistance when the two cytokinin-binding sites in GmPR08-Bet VI were mutated, indicating a new role of GmPR08-Bet VI in SCN resistance. GmPR08-Bet VI was mapped to a QTL for resistance to SCN using different mapping populations. GmSHMT08, GmSNAP18 and GmPR08-Bet VI localize to the cytosol and plasma membrane. GmSNAP18 expression and localization hyper-accumulated at the plasma membrane and was specific to the root cells surrounding the nematode in SCN-resistant soybeans. Genes encoding key components of the salicylic acid signalling pathway were induced under SCN infection. GmSNAP18 and GmPR08-Bet VI were also induced under salicylic acid and cytokinin exogenous treatments, while GmSHMT08 was induced only when the resistant GmSNAP18 was present, pointing to the presence of a molecular crosstalk between SCN-resistant genes and defence genes. Expression analysis of GmSHMT08 and GmSNAP18 identified the need of a minimum expression requirement to trigger the SCN resistance reaction. These results provide insight into a new response mechanism towards plant nematode resistance involving haplotype compatibility, gene dosage and hormone signalling.

Identification of introduced and stably inherited DNA methylation variants in soybean associated with soybean cyst nematode parasitism

DNA methylation is a widespread epigenetic mark that contributes to transcriptome reprogramming during plant-pathogen interactions. However, the distinct role of DNA methylation in establishing resistant and susceptible responses remains largely unexplored. Here, we developed and used a pair of near-isogenic lines (NILs) to characterize DNA methylome landscapes of soybean roots during the susceptible and resistant interactions with soybean cyst nematode (SCN; Heterodera glycines). We also compared the methylomes of the NILs and their parents to identify introduced and stably inherited methylation variants. The genomes of the NILs were substantially differentially methylated under uninfected conditions. This difference was associated with differential gene expression that may prime the NIL responses to SCN infection. In response to SCN infection, the susceptible line exhibited reduced global methylation levels in both protein-coding genes and transposable elements, whereas the resistant line showed the opposite response, increased global methylation levels. Heritable and novel nonparental differentially methylated regions overlapping with genes associated with soybean response to SCN infection were identified and validated using transgenic hairy root system. Our analyses indicate that DNA methylation patterns associated with the susceptible and resistant interactions are highly specific and that novel and stably inherited methylation variants are of biological significance.

Mutations at the Serine Hydroxymethyltransferase Impact its Interaction with a Soluble NSF Attachment Protein and a Pathogenesis-Related Protein in Soybean

Resistance to soybean cyst nematodes (SCN) in “Peking-type” resistance is bigenic, requiring Rhg4-a and rhg1-a. Rhg4-a encodes a serine hydroxymethyltransferase (GmSHMT08) and rhg1-a encodes a soluble NSF attachment protein (GmSNAP18). Recently, it has been shown that a pathogenesis-related protein, GmPR08-Bet VI, potentiates the interaction between GmSHMT08 and GmSNAP18. Mutational analysis using spontaneously occurring and ethyl methanesulfonate (EMS)-induced mutations was carried out to increase our knowledge of the interacting GmSHMT08/GmSNAP18/GmPR08-Bet VI multi-protein complex. Mutations affecting the GmSHMT08 protein structure (dimerization and tetramerization) and interaction sites with GmSNAP18 and GmPR08-Bet VI proteins were found to impact the multi-protein complex. Interestingly, mutations affecting the PLP/THF substrate binding and catalysis did not affect the multi-protein complex, although they resulted in increased susceptibility to SCN. Most importantly, GmSHMT08 and GmSNAP18 from PI88788 were shown to interact within the cell, being potentiated in the presence of GmPR08-Bet VI. In addition, we have shown the presence of incompatibility between the GmSNAP18 (rhg1-b) of PI88788 and GmSHMT08 (Rhg4-a) from Peking. Components of the reactive oxygen species (ROS) pathway were shown to be induced in the SCN incompatible reaction and were mapped to QTLs for resistance to SCN using different mapping populations.

Assessment of Phenotypic Variations and Correlation among Seed Composition Traits in Mutagenized Soybean Populations

Soybean [Glycine max (L.) Merr.] seed is a valuable source of protein and oil worldwide. Traditionally, the natural variations were heavily used in conventional soybean breeding programs to select desired traits. However, traditional plant breeding is encumbered with low frequencies of spontaneous mutations. In mutation breeding, genetic variations from induced mutations provide abundant sources of alterations in important soybean traits; this facilitated the development of soybean germplasm with modified seed composition traits to meet the different needs of end users. In this study, a total of 2366 'Forrest'-derived M2 families were developed for both forward and reverse genetic studies. A subset of 881 M3 families was forward genetically screened to measure the contents of protein, oil, carbohydrates, and fatty acids. A total of 14 mutants were identified to have stable seed composition phenotypes observed in both M3 and M4 generations. Correlation analyses have been conducted among ten seed composition traits and compared to a collection of 103 soybean germplasms. Mainly, ethyl methanesulfonate (EMS) mutagenesis had a strong impact on the seed-composition correlation that was observed among the 103 soybean germplasms, which offers multiple benefits for the soybean farmers and industry to breed for desired multiple seed phenotypes.

Whole-genome re-sequencing reveals the impact of the interaction of copy number variants of the rhg1 and Rhg4 genes on broad-based resistance to soybean cyst nematode

Soybean cyst nematode (SCN) is the most devastating plant-parasitic nematode. Most commercial soybean varieties with SCN resistance are derived from PI88788. Resistance derived from PI88788 is breaking down due to narrow genetic background and SCN population shift. PI88788 requires mainly the rhg1-b locus, while 'Peking' requires rhg1-a and Rhg4 for SCN resistance. In the present study, whole genome re-sequencing of 106 soybean lines was used to define the Rhg haplotypes and investigate their responses to the SCN HG-Types. The analysis showed a comprehensive profile of SNPs and copy number variations (CNV) at these loci. CNV of rhg1 (GmSNAP18) only contributed towards resistance in lines derived from PI88788 and 'Cloud'. At least 5.6 copies of the PI88788-type rhg1 were required to confer SCN resistance, regardless of the Rhg4 (GmSHMT08) haplotype. However, when the GmSNAP18 copies dropped below 5.6, a 'Peking'-type GmSHMT08 haplotype was required to ensure SCN resistance. This points to a novel mechanism of epistasis between GmSNAP18 and GmSHMT08 involving minimum requirements for copy number. The presence of more Rhg4 copies confers resistance to multiple SCN races. Moreover, transcript abundance of the GmSHMT08 in root tissue correlates with more copies of the Rhg4 locus, reinforcing SCN resistance. Finally, haplotype analysis of the GmSHMT08 and GmSNAP18 promoters inferred additional levels of the resistance mechanism. This is the first report revealing the genetic basis of broad-based resistance to SCN and providing new insight into epistasis, haplotype-compatibility, CNV, promoter variation and its impact on broad-based disease resistance in plants.

Genome reorganization of the GmSHMT gene family in soybean showed a lack of functional redundancy in resistance to soybean cyst nematode

In soybeans, eighteen members constitute the serine hydroxymethyltransferase (GmSHMT) gene family, of which the cytosolic-targeted GmSHMT08c member has been reported to mediate resistance to soybean cyst nematode (SCN). This work presents a comprehensive study of the SHMT gene family members, including synteny, phylogeny, subcellular localizations, haplotypes, protein homology modeling, mutational, and expression analyses. Phylogenetic analysis showed that SHMT genes are divided into four classes reflecting their subcellular distribution (cytosol, nucleus, mitochondrion, and chloroplast). Subcellular localization of selected GmSHMT members supports their in-silico predictions and phylogenetic distribution. Expression and functional analyses showed that GmSHMT genes display many overlapping, but some divergent responses during SCN infection. Furthermore, mutational analysis reveals that all isolated EMS mutants that lose their resistance to SCN carry missense and nonsense mutations at the GmSHMT08c, but none of the Gmshmt08c mutants carried mutations in the other GmSHMT genes. Haplotype clustering analysis using the whole genome resequencing data from a collection of 106 diverse soybean germplams (15X) was performed to identify allelic variants and haplotypes within the GmSHMT gene family. Interestingly, only the cytosolic-localized GmSHMT08c presented SNP clusters that were associated with SCN resistance, supporting our mutational analysis. Although eight GmSHMT members respond to the nematode infestation, functional and mutational analysis has shown the absence of functional redundancy in resistance to SCN. Structural analysis and protein homology modeling showed the presence of spontaneous mutations at important residues within the GmSHMT proteins, suggesting the presence of altered enzyme activities based on substrate affinities. Due to the accumulation of mutations during the evolution of the soybean genome, the other GmSHMT members have undergone neofunctionalization and subfunctionalization events.

Systematic Mutagenesis of Serine Hydroxymethyltransferase Reveals an Essential Role in Nematode Resistance

Rhg4 is a major genetic locus that contributes to soybean cyst nematode (SCN) resistance in the Peking-type resistance of soybean (Glycine max), which also requires the rhg1 gene. By map-based cloning and functional genomic approaches, we previously showed that the Rhg4 gene encodes a predicted cytosolic serine hydroxymethyltransferase (GmSHMT08); however, the novel gain of function of GmSHMT08 in SCN resistance remains to be characterized. Using a forward genetic screen, we identified an allelic series of GmSHMT08 mutants that shed new light on the mechanistic aspects of GmSHMT08-mediated resistance. The new mutants provide compelling genetic evidence that Peking-type rhg1 resistance in cv Forrest is fully dependent on the GmSHMT08 gene and demonstrates that this resistance is mechanistically different from the PI 88788-type of resistance that only requires rhg1 We also demonstrated that rhg1-a from cv Forrest, although required, does not exert selection pressure on the nematode to shift from HG type 7, which further validates the bigenic nature of this resistance. Mapping of the identified mutations onto the SHMT structural model uncovered key residues for structural stability, ligand binding, enzyme activity, and protein interactions, suggesting that GmSHMT08 has additional functions aside from its main enzymatic role in SCN resistance. Lastly, we demonstrate the functionality of the GmSHMT08 SCN resistance gene in a transgenic soybean plant.

Neodiversification of homeologous CLAVATA1-like receptor kinase genes in soybean leads to distinct developmental outcomes

The CLAVATA pathway that regulates stem cell numbers of the shoot apical meristem has exclusively been studied in Arabidopsis; as such insight into other species is warranted. In this study, a GmCLV1A mutant (F-S562L) with altered lateral organ development, and two mutants of GmNARK, isolated from a Forrest M2 population (EMS-mutated soybean) were studied. GmCLV1A and GmNARK encode for LRR receptor kinases, and share 92% of protein sequence. While GmNARK is critical for systemic regulation of nodulation (new organ made on the root through symbiosis), we show that GmCLV1A functions locally and has no apparent function in nodulation or root development. However, a recessive, loss-of-function mutation (S562L) in a putative S-glycosylation site of GmCLV1A causes stem nodal identity alterations as well as flower and pod abnormalities (deformed flower and pod). The mutant also exhibits a homeotic phenotype, displaying abnormal leaf development/number, vein-derived leaf emergence, and a thick, faciated stem. The mutant phenotype is also temperature-sensitive. Interestingly, a novel truncated version of GmCLV1A was identified upstream of GmCLV1A that is absent from GmNARK, but is present upstream of the GmNARK orthologues, MtSUNN and PvNARK. Taken together, our findings indicate that GmCLV1A acts on shoot architecture, whereas GmNARK, functions in controlling nodule numbers.

Stearoyl-Acyl Carrier Protein Desaturase Mutations Uncover an Impact of Stearic Acid in Leaf and Nodule Structure

Stearoyl-acyl carrier protein desaturase (SACPD-C) has been reported to control the accumulation of seed stearic acid; however, no study has previously reported its involvement in leaf stearic acid content and impact on leaf structure and morphology. A subset of an ethyl methanesulfonate mutagenized population of soybean (Glycine max) 'Forrest' was screened to identify mutants within the GmSACPD-C gene. Using a forward genetics approach, one nonsense and four missense Gmsacpd-c mutants were identified to have high levels of seed, nodule, and leaf stearic acid content. Homology modeling and in silico analysis of the GmSACPD-C enzyme revealed that most of these mutations were localized near or at conserved residues essential for diiron ion coordination. Soybeans carrying Gmsacpd-c mutations at conserved residues showed the highest stearic acid content, and these mutations were found to have deleterious effects on nodule development and function. Interestingly, mutations at nonconserved residues show an increase in stearic acid content yet retain healthy nodules. Thus, random mutagenesis and mutational analysis allows for the achievement of high seed stearic acid content with no associated negative agronomic characteristics. Additionally, expression analysis demonstrates that nodule leghemoglobin transcripts were significantly more abundant in soybeans with deleterious mutations at conserved residues of GmSACPD-C. Finally, we report that Gmsacpd-c mutations cause an increase in leaf stearic acid content and an alteration of leaf structure and morphology in addition to differences in nitrogen-fixing nodule structure.

Genetic Analysis of Root and Shoot Traits in the ‘Essex’ By ‘Forrest’ Recombinant Inbred Line (RIL) Population of Soybean [Glycine max (L.) Merr.]

Crop productivity is severely reduced by water deficit and drought in many plant species including soybean. Improved root and shoot traits can contribute to drought tolerance ability of the plant. This research was conducted to identify QTL that underlie several root and shoot traits in the ‘Essex’ by ‘Forrest’ (ExF RILs, n=94) recombinant inbred line (RIL) soybean population. Field collected samples were used for gathering phenotypic data of basal root thickness (BRT), lateral root number (LRN), maximum root length (MRL), root fresh weight (RFW), root dry weight (RDW), shoot fresh weight (SFW), shoot dry weight (SDW), and calculating RFW/SFW, and RDW/SDW ratios. All traits and ratios were compared against DNA markers using the composite interval mapping (CIM). A total of 12 QTL: 3 for MRL, 1 QTL for LRN, 1 QTL for BRT, 2 QTL for RFW, 2 QTL for RDW, 4 QTL for SFW, 3 QTL for SDW, and 3 QTL for SFW/SDW were identified and mapped on different linkage groups (LGs) A2, B2, C2, D1a, F, G, and N. The LOD scores of these QTL ranged from 2.5 to 6.0. No QTL were associated with RFW/RDW. The root and shoot trait QTL of this study may benefit breeding programs for producing cultivars tolerant to water deficit and high yield. Preliminary analyses of genes the QTL regions using GO annotation gave insight into genes that may underlie some of these QTLs.

The soybean GmSNAP18 gene underlies two types of resistance to soybean cyst nematode

Two types of resistant soybean (Glycine max (L.) Merr.) sources are widely used against soybean cyst nematode (SCN, Heterodera glycines Ichinohe). These include Peking-type soybean, whose resistance requires both the rhg1-a and Rhg4 alleles, and PI 88788-type soybean, whose resistance requires only the rhg1-b allele. Multiple copy number of PI 88788-type GmSNAP18, GmAAT, and GmWI12 in one genomic segment simultaneously contribute to rhg1-b resistance. Using an integrated set of genetic and genomic approaches, we demonstrate that the rhg1-a Peking-type GmSNAP18 is sufficient for resistance to SCN in combination with Rhg4. The two SNAPs (soluble NSF attachment proteins) differ by only five amino acids. Our findings suggest that Peking-type GmSNAP18 is performing a different role in SCN resistance than PI 88788-type GmSNAP18. As such, this is an example of a pathogen resistance gene that has evolved to underlie two types of resistance, yet ensure the same function within a single plant species.

Characterization of the Soluble NSF Attachment Protein gene family identifies two members involved in additive resistance to a plant pathogen

Proteins with Tetratricopeptide-repeat (TPR) domains are encoded by large gene families and distributed in all plant lineages. In this study, the Soluble NSF-Attachment Protein (SNAP) subfamily of TPR containing proteins is characterized. In soybean, five members constitute the SNAP gene family: GmSNAP18, GmSNAP11, GmSNAP14, GmSNAP02, and GmSNAP09. Recently, GmSNAP18 has been reported to mediate resistance to soybean cyst nematode (SCN). Using a population of recombinant inbred lines from resistant and susceptible parents, the divergence of the SNAP gene family is analysed over time. Phylogenetic analysis of SNAP genes from 22 diverse plant species showed that SNAPs were distributed in six monophyletic clades corresponding to the major plant lineages. Conservation of the four TPR motifs in all species, including ancestral lineages, supports the hypothesis that SNAPs were duplicated and derived from a common ancestor and unique gene still present in chlorophytic algae. Syntenic analysis of regions harbouring GmSNAP genes in soybean reveals that this family expanded from segmental and tandem duplications following a tetraploidization event. qRT-PCR analysis of GmSNAPs indicates a co-regulation following SCN infection. Finally, genetic analysis demonstrates that GmSNAP11 contributes to an additive resistance to SCN. Thus, GmSNAP11 is identified as a novel minor gene conferring resistance to SCN.

Characterization of the FAD2 Gene Family in Soybean Reveals the Limitations of Gel-Based TILLING in Genes with High Copy Number

Soybean seed oil typically contains 18-20% oleic acid. Increasing the content of oleic acid is beneficial for health and biodiesel production. Mutations in FAD2-1 genes have been reported to increase seed oleic acid content. A subset of 1,037 mutant families from a mutagenized soybean cultivar (cv.) Forrest population was screened using reverse genetics (TILLING) to identify mutations within FAD2 genes. Although no fad2 mutants were identified using gel-based TILLING, four fad2-1A and one fad2-1B mutants were identified to have high seed oleic acid content using forward genetic screening and subsequent target sequencing. TILLING has been successfully used as a non-transgenic reverse genetic approach to identify mutations in genes controlling important agronomic traits. However, this technique presents limitations in traits such as oil composition due to gene copy number and similarities within the soybean genome. In soybean, FAD2 are present as two copies, FAD2-1 and FAD2-2. Two FAD2-1 members: FAD2-1A and FAD2-1B; and three FAD2-2 members: FAD2-2A, FAD2-2B, and FAD2-2C have been reported. Syntenic, phylogenetic, and in silico analysis revealed two additional members constituting the FAD2 gene family: GmFAD2-2D and GmFAD2-2E, located on chromosomes 09 and 15, respectively. They are presumed to have diverged from other FAD2-2 members localized on chromosomes 19 (GmFAD2-2A and GmFAD2-2B) and 03 (GmFAD2-2C). This work discusses alternative solutions to the limitations of gel-based TILLING in functional genomics due to high copy number and multiple paralogs of the FAD2 gene family in soybean.

The Cardamine hirsuta genome offers insight into the evolution of morphological diversity

Finding causal relationships between genotypic and phenotypic variation is a key focus of evolutionary biology, human genetics and plant breeding. To identify genome-wide patterns underlying trait diversity, we assembled a high-quality reference genome of Cardamine hirsuta, a close relative of the model plant Arabidopsis thaliana. We combined comparative genome and transcriptome analyses with the experimental tools available in C. hirsuta to investigate gene function and phenotypic diversification. Our findings highlight the prevalent role of transcription factors and tandem gene duplications in morphological evolution. We identified a specific role for the transcriptional regulators PLETHORA5/7 in shaping leaf diversity and link tandem gene duplication with differential gene expression in the explosive seed pod of C. hirsuta. Our work highlights the value of comparative approaches in genetically tractable species to understand the genetic basis for evolutionary change.

A SNARE-Like Protein and Biotin Are Implicated in Soybean Cyst Nematode Virulence

Phytoparasitic nematodes that are able to infect and reproduce on plants that are considered resistant are referred to as virulent. The mechanism(s) that virulent nematodes employ to evade or suppress host plant defenses are not well understood. Here we report the use of a genetic strategy (allelic imbalance analysis) to associate single nucleotide polymorphisms (SNPs) with nematode virulence genes in Heterodera glycines, the soybean cyst nematode (SCN). To accomplish this analysis, a custom SCN SNP array was developed and used to genotype SCN F3-derived populations grown on resistant and susceptible soybean plants. Three SNPs reproducibly showed allele imbalances between nematodes grown on resistant and susceptible plants. Two candidate SCN virulence genes that were tightly linked to the SNPs were identified. One SCN gene encoded biotin synthase (HgBioB), and the other encoded a bacterial-like protein containing a putative SNARE domain (HgSLP-1). The two genes mapped to two different linkage groups. HgBioB contained sequence polymorphisms between avirulent and virulent nematodes. However, the gene encoding HgSLP-1 had reduced copy number in virulent nematode populations and appears to produce multiple forms of the protein via intron retention and alternative splicing. We show that HgSLP-1 is an esophageal-gland protein that is secreted by the nematode during plant parasitism. Furthermore, in bacterial co-expression experiments, HgSLP-1 co-purified with the SCN resistance protein Rhg1 α-SNAP, suggesting that these two proteins physically interact. Collectively our data suggest that multiple SCN genes are involved in SCN virulence, and that HgSLP-1 may function as an avirulence protein and when absent it helps SCN evade host defenses.

SNP identification and marker assay development for high-throughput selection of soybean cyst nematode resistance

BACKGROUND

Soybean cyst nematode (SCN) is the most economically devastating pathogen of soybean. Two resistance loci, Rhg1 and Rhg4 primarily contribute resistance to SCN race 3 in soybean. Peking and PI 88788 are the two major sources of SCN resistance with Peking requiring both Rhg1 and Rhg4 alleles and PI 88788 only the Rhg1 allele. Although simple sequence repeat (SSR) markers have been reported for both loci, they are linked markers and limited to be applied in breeding programs due to accuracy, throughput and cost of detection methods. The objectives of this study were to develop robust functional marker assays for high-throughput selection of SCN resistance and to differentiate the sources of resistance.

RESULTS

Based on the genomic DNA sequences of 27 soybean lines with known SCN phenotypes, we have developed Kompetitive Allele Specific PCR (KASP) assays for two Single nucleotide polymorphisms (SNPs) from Glyma08g11490 for the selection of the Rhg4 resistance allele. Moreover, the genomic DNA of Glyma18g02590 at the Rhg1 locus from 11 soybean lines and cDNA of Forrest, Essex, Williams 82 and PI 88788 were fully sequenced. Pairwise sequence alignment revealed seven SNPs/insertion/deletions (InDels), five in the 6th exon and two in the last exon. Using the same 27 soybean lines, we identified one SNP that can be used to select the Rhg1 resistance allele and another SNP that can be employed to differentiate Peking and PI 88788-type resistance. These SNP markers have been validated and a strong correlation was observed between the SNP genotypes and reactions to SCN race 3 using a panel of 153 soybean lines, as well as a bi-parental population, F5-derived recombinant inbred lines (RILs) from G00-3213xLG04-6000.

CONCLUSIONS

Three functional SNP markers (two for Rhg1 locus and one for Rhg4 locus) were identified that could provide genotype information for the selection of SCN resistance and differentiate Peking from PI 88788 source for most germplasm lines. The robust KASP SNP marker assays were developed. In most contexts, use of one or two of these markers is sufficient for high-throughput marker-assisted selection of plants that will exhibit SCN resistance.

Quantitative trait loci underlying resistance to sudden death syndrome (SDS) in MD96-5722 by 'Spencer' recombinant inbred line population of soybean

The best way to protect yield loss of soybean [Glycine max (L.) Merr.] due to sudden death syndrome (SDS), caused by Fusarium virguliforme (Aoki, O'Donnel, Homma & Lattanzi), is the development and use of resistant lines. Mapping quantitative trait loci (QTL) linked to SDS help developing resistant soybean germplasm through molecular marker-assisted selection strategy. QTL for SDS presented herein are from a high-density SNP-based genetic linkage map of MD 96-5722 (a.k.a 'Monocacy') by 'Spencer' recombinant inbred line using SoySNP6K Illumina Infinium BeadChip genotyping array. Ninety-four F5:7 lines were evaluated for 2 years (2010 and 2011) at two locations (Carbondale and Valmeyer) in southern Illinois, USA to identify QTL controlling SDS resistance using disease index (DX). Composite interval mapping identified 19 SDS controlling QTL which were mapped on 11 separate linkage group (LG) or chromosomes (Chr) out of 20 LG or Chr of soybean genome. Many of these significant QTL identified in one environment/year were confirmed in another year or environment, which suggests a common genetic effects and modes of the pathogen. These new QTL are useful sources for SDS resistance studies in soybean breeding, complementing previously reported loci.

Quantitative trait loci for seed isoflavone contents in 'MD96-5722' by 'Spencer' recombinant inbred lines of soybean

Isoflavones from soybeans [ Glycine max (L.) Merr.] have a significant impact on human health to reduce the risk of several major diseases. Breeding soybean for high isoflavone content in the seed is possible through marker-assisted selection (MAS) which can be based on quantitative trait loci (QTL). The objective of this study was to identify QTL controlling isoflavone content in a set of 'MD96-5722' by 'Spencer' recombinant inbred line (RIL) populations of soybean. Wide variations were found for seed concentrations of daidzein, glycitein, genistein, and total isoflavones among the RIL populations. Three QTL were identified on three different linkage groups (LG) represented by three different chromosomes (Chr). One QTL that controlled daidzein content was identified on LG A1 (Chr 5), and two QTL that underlay glycitein content were identified on LG K (Chr 9) and LG B2 (Chr 14). Identified QTL could be functional in developing soybean with preferable isoflavone concentrations in the seeds through MAS.

Quantitative trait loci for seed isoflavone contents in 'MD96-5722' by 'Spencer' recombinant inbred lines of soybean

Isoflavones from soybeans [ Glycine max (L.) Merr.] have a significant impact on human health to reduce the risk of several major diseases. Breeding soybean for high isoflavone content in the seed is possible through marker-assisted selection (MAS) which can be based on quantitative trait loci (QTL). The objective of this study was to identify QTL controlling isoflavone content in a set of 'MD96-5722' by 'Spencer' recombinant inbred line (RIL) populations of soybean. Wide variations were found for seed concentrations of daidzein, glycitein, genistein, and total isoflavones among the RIL populations. Three QTL were identified on three different linkage groups (LG) represented by three different chromosomes (Chr). One QTL that controlled daidzein content was identified on LG A1 (Chr 5), and two QTL that underlay glycitein content were identified on LG K (Chr 9) and LG B2 (Chr 14). Identified QTL could be functional in developing soybean with preferable isoflavone concentrations in the seeds through MAS.

A soybean cyst nematode resistance gene points to a new mechanism of plant resistance to pathogens

Soybean (Glycine max (L.) Merr.) is an important crop that provides a sustainable source of protein and oil worldwide. Soybean cyst nematode (Heterodera glycines Ichinohe) is a microscopic roundworm that feeds on the roots of soybean and is a major constraint to soybean production. This nematode causes more than US$1 billion in yield losses annually in the United States alone, making it the most economically important pathogen on soybean. Although planting of resistant cultivars forms the core management strategy for this pathogen, nothing is known about the nature of resistance. Moreover, the increase in virulent populations of this parasite on most known resistance sources necessitates the development of novel approaches for control. Here we report the map-based cloning of a gene at the Rhg4 (for resistance to Heterodera glycines 4) locus, a major quantitative trait locus contributing to resistance to this pathogen. Mutation analysis, gene silencing and transgenic complementation confirm that the gene confers resistance. The gene encodes a serine hydroxymethyltransferase, an enzyme that is ubiquitous in nature and structurally conserved across kingdoms. The enzyme is responsible for interconversion of serine and glycine and is essential for cellular one-carbon metabolism. Alleles of Rhg4 conferring resistance or susceptibility differ by two genetic polymorphisms that alter a key regulatory property of the enzyme. Our discovery reveals an unprecedented plant resistance mechanism against a pathogen. The mechanistic knowledge of the resistance gene can be readily exploited to improve nematode resistance of soybean, an increasingly important global crop.

Soybean cyst nematode resistance in soybean is independent of the Rhg4 locus LRR-RLK gene

To test the function of candidate genes in soybean for resistance to the soybean cyst nematode (SCN), a large collection of EMS-mutants from the SCN-resistant soybean cultivar "Forrest" was developed for Targeting Induced Local Lesions IN Genomes (TILLING). Additionally, due to the complexity of the soybean genome, an integrated set of genomic and genetic analysis tools was employed to complement the TILLING approach. The efficiency of this integrated set of tools was tested using a candidate soybean gene for resistance to SCN, encoding a leucine-rich repeat receptor-like kinase (LRR-RLK) that was identified by map-based cloning at the Rhg4 locus. The Rhg4 locus is one of the major quantitative trait loci controlling soybean resistance against SCN race 3 (HG type 0) in cv. Forrest, but the gene(s) sequence for resistance remains to be determined. Using TILLING, a Forrest mutant containing a nonsense mutation in the LRR domain of the candidate resistance protein was identified and confirmed; however, the SCN-resistant phenotype of the mutant was not altered. Haplotyping and EcoTILLING of recombinant inbred lines along with complementation analysis corroborated the TILLING result and ruled out the possibility of functional redundancy by a second copy of the LRR-RLK gene identified in the soybean genome. This study validates the use of TILLING, in combination with an integrated set of genomic tools, as an efficient means of testing candidate genes for SCN resistance in soybean.

Characterization of a plant-transformation-ready large-insert BIBAC library of Arabidopsis and bombardment transformation of a large-insert BIBAC of the library into tobacco

Plant-transformation-ready, large-insert binary bacterial artificial chromosome (BIBAC) libraries are of significance for functional and network analysis of large genomic regions, gene clusters, large-spanning genes, and complex loci in the post-genome era. Here, we report the characterization of a plant-transformation-ready BIBAC library of the sequenced Arabidopsis genome for which such a library is not available to the public, the transformation of a large-insert BIBAC of the library into tobacco by biolistic bombardment, and the expression analysis of its containing genes in transgenic plants. The BIBAC library was constructed from nuclear DNA partially digested with BamHI in the BIBAC vector pCLD04541. It contains 6144 clones and has a mean insert size of 108 kb, representing 5.2× equivalents of the Arabidopsis genome or a probability of greater than 99% of obtaining at least one positive clone from the library using a single-copy sequence as a probe. The transformation of the large-insert BIBAC and analyses of the transgenic plants showed that not only did transgenic plants have intact BIBAC DNA, but also could the BIBAC be transmitted stably into progenies and its containing genes be expressed actively. These results suggest that the large-insert BIBAC library, combined with the biolistic bombardment transformation method, could provide a useful tool for large-scale functional analysis of the Arabidopsis genome sequence and applications in plant-molecular breeding.

Multigeneic QTL: the laccase encoded within the soybean Rfs2/rhg1 locus inferred to underlie part of the dual resistance to cyst nematode and sudden death syndrome

Multigeneic QTL present significant problems to analysis. Resistance to soybean (Glycine max (L) Merr.) sudden death syndrome (SDS) caused by Fusarium virguliforme was partly underlain by QRfs2 that was clustered with, or pleiotropic to, the multigeneic rhg1 locus providing resistance to soybean cyst nematode (SCN; Heterodera glycines). A group of five genes were found between the two markers that delimited the Rfs2/rhg1 locus. One of the five genes was predicted to encode an unusual diphenol oxidase (laccase; EC 1.10.3.2). The aim of this study was to characterize this member of the soybean laccase gene-family and explore its involvement in SDS resistance. A genomic clone and a full length cDNA was isolated from resistant cultivar 'Forrest' that were different among susceptible cultivars 'Asgrow 3244' and 'Williams 82' at four residues R/H168, I/M271, R/H330, E/K470. Additional differences were found in six of the seven introns and the promoter region. Transcript abundance (TA) among genotypes that varied for resistance to SDS or SCN did not differ significantly. Therefore the protein activity was inferred to underlie resistance. Protein expressed in yeast pYES2/NTB had weak enzyme activity with common substrates but good activity with root phenolics. The Forrest isoform may underlie both QRfs2 and rhg1.

Fine mapping and DNA fiber FISH analysis locates the tobamovirus resistance gene L3 of Capsicum chinense in a 400-kb region of R-like genes cluster embedded in highly repetitive sequences

The tobamovirus resistance gene L(3) of Capsicum chinense was mapped using an intra-specific F2 population (2,016 individuals) of Capsicum annuum cultivars, into one of which had been introduced the C. chinense L(3) gene, and an inter-specific F2 population (3,391 individuals) between C. chinense and Capsicum frutescence. Analysis of a BAC library with an AFLP marker closely linked to L(3)-resistance revealed the presence of homologs of the tomato disease resistance gene I2. Partial or full-length coding sequences were cloned by degenerate PCR from 35 different pepper I2 homologs and 17 genetic markers were generated in the inter-specific combination. The L(3) gene was mapped between I2 homolog marker IH1-04 and BAC-end marker 189D23M, and located within a region encompassing two different BAC contigs consisting of four and one clones, respectively. DNA fiber FISH analysis revealed that these two contigs are separated from each other by about 30 kb. DNA fiber FISH results and Southern blotting of the BAC clones suggested that the L(3) locus-containing region is rich in highly repetitive sequences. Southern blot analysis indicated that the two BAC contigs contain more than ten copies of the I2 homologs. In contrast to the inter-specific F2 population, no recombinant progeny were identified to have a crossover point within two BAC contigs consisting of seven and two clones in the intra-specific F2 population. Moreover, distribution of the crossover points differed between the two populations, suggesting linkage disequilibrium in the region containing the L locus.

In silico comparison of transcript abundances during Arabidopsis thaliana and Glycine max resistance to Fusarium virguliforme

BACKGROUND

Sudden death syndrome (SDS) of soybean (Glycine max L. Merr.) is an economically important disease, caused by the semi-biotrophic fungus Fusarium solani f. sp. glycines, recently renamed Fusarium virguliforme (Fv). Due to the complexity and length of the soybean-Fusarium interaction, the molecular mechanisms underlying plant resistance and susceptibility to the pathogen are not fully understood. F. virguliforme has a very wide host range for the ability to cause root rot and a very narrow host range for the ability to cause a leaf scorch. Arabidopsis thaliana is a host for many types of phytopathogens including bacteria, fungi, viruses and nematodes. Deciphering the variations among transcript abundances (TAs) of functional orthologous genes of soybean and A. thaliana involved in the interaction will provide insights into plant resistance to F. viguliforme.

RESULTS

In this study, we reported the analyses of microarrays measuring TA in whole plants after A. thaliana cv 'Columbia' was challenged with fungal pathogen F. virguliforme. Infection caused significant variations in TAs. The total number of increased transcripts was nearly four times more than that of decreased transcripts in abundance. A putative resistance pathway involved in responding to the pathogen infection in A. thaliana was identified and compared to that reported in soybean.

CONCLUSION

Microarray experiments allow the interrogation of tens of thousands of transcripts simultaneously and thus, the identification of plant pathways is likely to be involved in plant resistance to Fusarial pathogens. Dissection of the set functional orthologous genes between soybean and A. thaliana enabled a broad view of the functional relationships and molecular interactions among plant genes involved in F. virguliforme resistance.

TILLING to detect induced mutations in soybean

BACKGROUND

Soybean (Glycine max L. Merr.) is an important nitrogen-fixing crop that provides much of the world's protein and oil. However, the available tools for investigation of soybean gene function are limited. Nevertheless, chemical mutagenesis can be applied to soybean followed by screening for mutations in a target of interest using a strategy known as Targeting Induced Local Lesions IN Genomes (TILLING). We have applied TILLING to four mutagenized soybean populations, three of which were treated with ethyl methanesulfonate (EMS) and one with N-nitroso-N-methylurea (NMU).

RESULTS

We screened seven targets in each population and discovered a total of 116 induced mutations. The NMU-treated population and one EMS mutagenized population had similar mutation density (approximately 1/140 kb), while another EMS population had a mutation density of approximately 1/250 kb. The remaining population had a mutation density of approximately 1/550 kb. Because of soybean's polyploid history, PCR amplification of multiple targets could impede mutation discovery. Indeed, one set of primers tested in this study amplified more than a single target and produced low quality data. To address this problem, we removed an extraneous target by pretreating genomic DNA with a restriction enzyme. Digestion of the template eliminated amplification of the extraneous target and allowed the identification of four additional mutant alleles compared to untreated template.

CONCLUSION

The development of four independent populations with considerable mutation density, together with an additional method for screening closely related targets, indicates that soybean is a suitable organism for high-throughput mutation discovery even with its extensively duplicated genome.

Construction and characterization of two bacterial artificial chromosome libraries of pea (Pisum sativum L.) for the isolation of economically important genes

Pea (Pisum sativum L.) has a genome of about 4 Gb that appears to share conserved synteny with model legumes having genomes of 0.2-0.4 Gb despite extensive intergenic expansion. Pea plant inventory (PI) accession 269818 has been used to introgress genetic diversity into the cultivated germplasm pool. The aim here was to develop pea bacterial artificial chromosome (BAC) libraries that would enable the isolation of genes involved in plant disease resistance or control of economically important traits. The BAC libraries encompassed about 3.2 haploid genome equivalents consisting of partially HindIII-digested DNA fragments with a mean size of 105 kb that were inserted in 1 of 2 vectors. The low-copy oriT-based T-DNA vector (pCLD04541) library contained 55 680 clones. The single-copy oriS-based vector (pIndigoBAC-5) library contained 65 280 clones. Colony hybridization of a universal chloroplast probe indicated that about 1% of clones in the libraries were of chloroplast origin. The presence of about 0.1% empty vectors was inferred by white/blue colony plate counts. The usefulness of the libraries was tested by 2 replicated methods. First, high-density filters were probed with low copy number sequences. Second, BAC plate-pool DNA was used successfully to PCR amplify 7 of 9 published pea resistance gene analogs (RGAs) and several other low copy number pea sequences. Individual BAC clones encoding specific sequences were identified. Therefore, the HindIII BAC libraries of pea, based on germplasm accession PI 269818, will be useful for the isolation of genes underlying disease resistance and other economically important traits.

Comparative sequence analysis of Solanum and Arabidopsis in a hot spot for pathogen resistance on potato chromosome V reveals a patchwork of conserved and rapidly evolving genome segments

Background

Quantitative phenotypic variation of agronomic characters in crop plants is controlled by environmental and genetic factors (quantitative trait loci = QTL). To understand the molecular basis of such QTL, the identification of the underlying genes is of primary interest and DNA sequence analysis of the genomic regions harboring QTL is a prerequisite for that. QTL mapping in potato (Solanum tuberosum) has identified a region on chromosome V tagged by DNA markers GP21 and GP179, which contains a number of important QTL, among others QTL for resistance to late blight caused by the oomycete Phytophthora infestans and to root cyst nematodes.

Results

To obtain genomic sequence for the targeted region on chromosome V, two local BAC (bacterial artificial chromosome) contigs were constructed and sequenced, which corresponded to parts of the homologous chromosomes of the diploid, heterozygous genotype P6/210. Two contiguous sequences of 417,445 and 202,781 base pairs were assembled and annotated. Gene-by-gene co-linearity was disrupted by non-allelic insertions of retrotransposon elements, stretches of diverged intergenic sequences, differences in gene content and gene order. The latter was caused by inversion of a 70 kbp genomic fragment. These features were also found in comparison to orthologous sequence contigs from three homeologous chromosomes of Solanum demissum, a wild tuber bearing species. Functional annotation of the sequence identified 48 putative open reading frames (ORF) in one contig and 22 in the other, with an average of one ORF every 9 kbp. Ten ORFs were classified as resistance-gene-like, 11 as F-box-containing genes, 13 as transposable elements and three as transcription factors. Comparing potato to Arabidopsis thaliana annotated proteins revealed five micro-syntenic blocks of three to seven ORFs with A. thaliana chromosomes 1, 3 and 5.

Conclusion

Comparative sequence analysis revealed highly conserved collinear regions that flank regions showing high variability and tandem duplicated genes. Sequence annotation revealed that the majority of the ORFs were members of multiple gene families. Comparing potato to Arabidopsis thaliana annotated proteins suggested fragmented structural conservation between these distantly related plant species.