The receptor like kinase at Rhg1-a/Rfs2 caused pleiotropic resistance to sudden death syndrome and soybean cyst nematode as a transgene by altering signaling responses

Deciphering the genetic architecture of resistance to Corynespora cassiicola in soybean (Glycine max L.) by integrating genome-wide association mapping and RNA-Seq analysis

Target spot caused by Corynespora cassiicola is a problematic disease in tropical and subtropical soybean (Glycine max) growing regions. Although resistant soybean genotypes have been identified, the genetic mechanisms underlying target spot resistance has not yet been studied. To address this knowledge gap, this is the first genome-wide association study (GWAS) conducted using the SoySNP50K array on a panel of 246 soybean accessions, aiming to unravel the genetic architecture of resistance. The results revealed significant associations of 14 and 33 loci with resistance to LIM01 and SSTA C. cassiicola isolates, respectively, with six loci demonstrating consistent associations across both isolates. To identify potential candidate genes within GWAS-identified loci, dynamic transcriptome profiling was conducted through RNA-Seq analysis. The analysis involved comparing gene expression patterns between resistant and susceptible genotypes, utilizing leaf tissue collected at different time points after inoculation. Integrating results of GWAS and RNA-Seq analyses identified 238 differentially expressed genes within a 200 kb region encompassing significant quantitative trait loci (QTLs) for disease severity ratings. These genes were involved in defense response to pathogen, innate immune response, chitinase activity, histone H3-K9 methylation, salicylic acid mediated signaling pathway, kinase activity, and biosynthesis of flavonoid, jasmonic acid, phenylpropanoid, and wax. In addition, when combining results from this study with previous GWAS research, 11 colocalized regions associated with disease resistance were identified for biotic and abiotic stress. This finding provides valuable insight into the genetic resources that can be harnessed for future breeding programs aiming to enhance soybean resistance against target spot and other diseases simultaneously.

Breeding for disease resistance in soybean: a global perspective

KEY MESSAGE

This review provides a comprehensive atlas of QTLs, genes, and alleles conferring resistance to 28 important diseases in all major soybean production regions in the world. Breeding disease-resistant soybean [Glycine max (L.) Merr.] varieties is a common goal for soybean breeding programs to ensure the sustainability and growth of soybean production worldwide. However, due to global climate change, soybean breeders are facing strong challenges to defeat diseases. Marker-assisted selection and genomic selection have been demonstrated to be successful methods in quickly integrating vertical resistance or horizontal resistance into improved soybean varieties, where vertical resistance refers to R genes and major effect QTLs, and horizontal resistance is a combination of major and minor effect genes or QTLs. This review summarized more than 800 resistant loci/alleles and their tightly linked markers for 28 soybean diseases worldwide, caused by nematodes, oomycetes, fungi, bacteria, and viruses. The major breakthroughs in the discovery of disease resistance gene atlas of soybean were also emphasized which include: (1) identification and characterization of vertical resistance genes reside rhg1 and Rhg4 for soybean cyst nematode, and exploration of the underlying regulation mechanisms through copy number variation and (2) map-based cloning and characterization of Rps11 conferring resistance to 80% isolates of Phytophthora sojae across the USA. In this review, we also highlight the validated QTLs in overlapping genomic regions from at least two studies and applied a consistent naming nomenclature for these QTLs. Our review provides a comprehensive summary of important resistant genes/QTLs and can be used as a toolbox for soybean improvement. Finally, the summarized genetic knowledge sheds light on future directions of accelerated soybean breeding and translational genomics studies.

Enhanced resistance to soybean cyst nematode in transgenic soybean via host-induced silencing of vital Heterodera glycines genes

Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is the most economically damaging pathogen affecting soybean production worldwide. Host-induced gene silencing provides a promising approach to confer resistance to plant parasitic nematodes. In the present study, we produced stable transgenic soybean plants individually harboring the inverted repeats of three essential H. glycines genes, Hg-rps23, Hg-snb1, and Hg-cpn1, and evaluated their resistance to SCN infection. Molecular characterization confirmed the stable integration of the hairpin double stranded (ds) RNA in host plants. Inoculation assays with SCN race 3 showed significant reduction of female index (FI, 11.84 ~ 17.47%) on the roots of T4 transgenic plants, with 73.29 ~ 81.90% reduction for the three RNA interference (RNAi) constructs, compared to non-transformed plants (NT, 65.43%). Enhanced resistance to SCN race 3 was further confirmed in subsequent generations (T5) of transgenic soybean. Moreover, when inoculated with SCN race 4 which was considered highly virulent to most of soybean germplasms and varieties, transgenic soybean plants also exhibited reduced FIs (9.96 ~ 23.67%) and increased resistance, relative to the NT plants (46.46%). Consistently, significant down-regulation in transcript levels of the Hg-rps23, Hg-snb1, Hg-cpn1 genes were observed in the nematodes feeding on the transgenic roots, suggesting a broad-spectrum resistance mediated by the host-mediated silencing of vital H. glycines genes. There were no significant differences in morphological traits between transgenic and NT soybean plants under conditions with negligible SCN infection. In summary, our results demonstrate the effectiveness of host-induced silencing of essential H. glycines genes to enhance broad-spectrum SCN resistance in stable transgenic soybean plants, without negative consequences on the agronomic performance.

Genetic Diversity and Selection Footprints in the Genome of Brazilian Soybean Cultivars

Although Brazil is currently the largest soybean producer in the world, only a small number of studies have analyzed the genetic diversity of Brazilian soybean. These studies have shown the existence of a narrow genetic base. The objectives of this work were to analyze the population structure and genetic diversity, and to identify selection signatures in the genome of soybean germplasms from different companies in Brazil. A panel consisting of 343 soybean lines from Brazil, North America, and Asia was genotyped using genotyping by sequencing (GBS). Population structure was assessed by Bayesian and multivariate approaches. Genetic diversity was analyzed using metrics such as the fixation index, nucleotide diversity, genetic dissimilarity, and linkage disequilibrium. The software BayeScan was used to detect selection signatures between Brazilian and Asian accessions as well as among Brazilian germplasms. Region of origin, company of origin, and relative maturity group (RMG) all had a significant influence on population structure. Varieties belonging to the same company and especially to the same RMG exhibited a high level of genetic similarity. This result was exacerbated among early maturing accessions. Brazilian soybean showed significantly lower genetic diversity when compared to Asian accessions. This was expected, because the crop's region of origin is its main genetic diversity reserve. We identified 7 genomic regions under selection between the Brazilian and Asian accessions, and 27 among Brazilian varieties developed by different companies. Associated with these genomic regions, we found 96 quantitative trait loci (QTLs) for important soybean breeding traits such as flowering, maturity, plant architecture, productivity components, pathogen resistance, and seed composition. Some of the QTLs associated with the markers under selection have genes of great importance to soybean's regional adaptation. The results reported herein allowed to expand the knowledge about the organization of the genetic variability of the Brazilian soybean germplasm. Furthermore, it was possible to identify genomic regions under selection possibly associated with the adaptation of soybean to Brazilian environments.

Molecular mapping of quantitative disease resistance loci for soybean partial resistance to Phytophthora sansomeana

KEY MESSAGE

Two soybean QDRL were identified with additive interaction to P. sansomeana isolate MPS17-22. Further analyses uncovered four interaction patterns between the two QDRL and seven additional P. sansomeana isolates. Phytophthora sansomeana is a recently recognized species that contributes to root rot in soybean. Previous studies indicated that P. sansomeana is widely distributed among soybean growing regions and has a much wider host range than P. sojae, a well-known pathogen of soybean. Unlike P. sojae, no known disease resistance genes have been documented that can effectively control P. sansomeana. Therefore, it is important to identify resistance that can be quickly integrated into future soybean varieties. E13901 is an improved soybean line that confers partial resistance to P. sansomeana. A mapping population of 228 F4:5 families was developed from a cross between E13901 and a susceptible improved soybean variety E13390. Using a composite interval mapping method, two quantitative disease resistance loci (QDRL) were identified on Chromosomes 5 (designated qPsan5.1) and 16 (designated qPsan16.1), respectively. qPsan5.1 was mapped at 54.71 cM between Gm05_32565157_T_C and Gm05_32327497_T_C. qPsan5.1 was contributed by E13390 and explained about 6% of the disease resistance variation. qPsan16.1 was located at 39.01 cM between Gm16_35700223_G_T and Gm16_35933600/ Gm16_35816475. qPsan16.1 was from E13901 and could explain 5.5% of partial disease resistance. Further analysis indicated an additive interaction of qPsan5.1 and qPsan16.1 against P. sansomeana isolate MPS17-22. Marker assisted resistance spectrum analysis and progeny tests verified the two QDRL and their interaction patterns with other P. sansomeana isolates. Both QDRL can be quickly integrated into soybean varieties using marker assisted selection.

Comparative genomic analyses of two segregating mutants reveal seven genes likely involved in resistance to Fusarium equiseti in soybean via whole genome re-sequencing

KEY MESSAGE

The candidate genes involved in resistance to Fusarium equiseti in soybean PI 437654 were identified through comparative genomic analyses of mutants via whole genome re-sequencing. The fungus Fusarium infects each stage of the growth and development of soybean and causes soybean (Glycine max (L.)) seed and root rot and seedling damping-off and wilt with a large quantity of annual yield loss worldwide. It is very important to identify the resistant genes in soybean to prevent and control this pathogen. One Fusarium equiseti isolate was previously identified to be incompatible with 'PI 437654' but compatible with a Chinese soybean cultivar 'Zhonghuang 13'. In this study, with the infection of this isolate on the seedling roots of developed PI 437654 mutants, 6 mutants were identified from 500 mutants to significantly alter their phenotypes to F. equiseti compared to wild-type PI 437654. Then, two identified segregating mutants were selected to directly perform whole genome re-sequencing. Finally, through comparative genomic analyses 7 genes including one cluster of 4 nucleotide binding site-leucine-rich repeat genes on one genomic region of chromosome 7, a 60S ribosomal protein L12 gene and 2 uncharacterized genes were identified to be likely involved in the resistance to F. equiseti. These genes will facilitate the breeding of resistant germplasm resources and the identification of resistance of soybean to Fusarium spp.

Predicting Soybean Yield and Sudden Death Syndrome Development Using At-Planting Risk Factors

In the United States, sudden death syndrome (SDS) of soybean is caused by the fungal pathogen Fusarium virguliforme and is responsible for important yield losses each year. Understanding the risk of SDS development and subsequent yield loss could provide growers with valuable information for management of this challenging disease. Current management strategies for F. virguliforme use partially resistant cultivars, fungicide seed treatments, and extended crop rotations with diverse crops. The aim of this study was to develop models to predict SDS severity and soybean yield loss using at-planting risk factors to integrate with current SDS management strategies. In 2014 and 2015, field studies were conducted in adjacent fields in Decatur, MI, which were intensively monitored for F. virguliforme and nematode quantities at-planting, plant health throughout the growing season, end-of-season SDS severity, and yield using an unbiased grid sampling scheme. In both years, F. virguliforme and soybean cyst nematode (SCN) quantities were unevenly distributed throughout the field. The distribution of F. virguliforme at-planting had a significant correlation with end-of-season SDS severity in 2015, and a significant correlation to yield in 2014 (P < 0.05). SCN distributions at-planting were significantly correlated with end-of-season SDS severity and yield in 2015 (P < 0.05). Prediction models developed through multiple linear regression showed that F. virguliforme abundance (P < 0.001), SCN egg quantity (P < 0.001), and year (P < 0.01) explained the most variation in end-of-season SDS (R2 = 0.32), whereas end-of-season SDS (P < 0.001) and end-of-season root dry weight (P < 0.001) explained the most variation in soybean yield (R2 = 0.53). Further, multivariate analyses support a synergistic relationship between F. virguliforme and SCN, enhancing the severity of foliar SDS. These models indicate that it is possible to predict patches of SDS severity using at-planting risk factors. Verifying these models and incorporating additional data types may help improve SDS management and forecast soybean markets in response to SDS threats.

Molecular Basis of Soybean Resistance to Soybean Aphids and Soybean Cyst Nematodes

Soybean aphid (SBA; Aphis glycines Matsumura) and soybean cyst nematode (SCN; Heterodera glycines Ichninohe) are major pests of the soybean (Glycine max [L.] Merr.). Substantial progress has been made in identifying the genetic basis of limiting these pests in both model and non-model plant systems. Classical linkage mapping and genome-wide association studies (GWAS) have identified major and minor quantitative trait loci (QTLs) in soybean. Studies on interactions of SBA and SCN effectors with host proteins have identified molecular cues in various signaling pathways, including those involved in plant disease resistance and phytohormone regulations. In this paper, we review the molecular basis of soybean resistance to SBA and SCN, and we provide a synthesis of recent studies of soybean QTLs/genes that could mitigate the effects of virulent SBA and SCN populations. We also review relevant studies of aphid–nematode interactions, particularly in the soybean–SBA–SCN system.

Effective identification of soybean candidate genes involved in resistance to soybean cyst nematode via direct whole genome re-sequencing of two segregating mutants

KEY MESSAGE

Three soybean candidate genes involved in resistance to soybean cyst nematode race 4 were identified via direct whole genome re-sequencing of two segregating mutants. The genes conferring resistance to soybean cyst nematode (SCN) race 4 (Hg type 1.2.3.5.7) in soybean (Glycine max L. Merr.) remains unknown. Next generation sequencing-based methods identify a wide range of targets, it is difficult to identify genes underlying traits. Use of the MutMap and QTL-seq methods to identify trait candidate genes needs backcrossing and is very time-consuming. Here we report a simple method to effectively identify candidate genes involved in resistance to SCN race 4. Two ethane methylsulfonate mutagenized mutants of soybean 'PI 437654', whose SCN race 4-infection phenotype altered, were selected. Six relevant whole genomes were re-sequenced, and then calling of genomic variants (SNPs and InDels) was conducted and compared to 'Williams 82'. The comparison eliminated many genomic variants from the mutant lines that overlapped two non-phenotypic but mutant progeny plants, wild-type PI 437654 and 'Zhonghuang 13'. Finally, only 27 mutations were found among 10 genes. Of these 10 genes, 3 genes, Glyma.09g054000, Glyma.16g065700 and Glyma.18g192200 were overlapped between two phenotypic mutant progeny plants. Therefore, the three genes may be the candidate genes involved in resistance of PI 437654 to soybean cyst nematode race 4. This method simplifies the effective identification of candidate genes.

Temporal Dynamics of Fusarium virguliforme Colonization of Soybean Roots

Soybean sudden death syndrome (SDS) caused by Fusarium virguliforme is one of the most yield limiting soybean diseases in the United States. SDS disease symptoms include root rot and foliar symptoms induced by fungal toxins. Soybean cultivar resistance is one of the most effective SDS disease management options, but no cultivar displays complete resistance. Soybean SDS foliar symptoms are the primary phenotype used to screen and breed for SDS resistance. Root rot or root colonization measures are seldom utilized, partly due to the lack of convenient and accurate methods for quantification of F. virguliforme. In this study, greenhouse and field experiments were conducted to determine the temporal dynamics of F. virguliforme colonization of soybean roots using quantitative real-time PCR (qPCR). The infection coefficient (IC), or ratio of F. virguliforme DNA to soybean DNA, was determined in soybean cultivars with different SDS foliar resistance ratings. In greenhouse experiments, F. virguliforme was detected in all cultivars 7 days after planting (DAP), with a peak in IC at 14 DAP. All soybean cultivars developed SDS foliar symptoms, but F. virguliforme soybean root colonization levels did not significantly correlate with SDS foliar symptom severity. In field experiments, SDS foliar symptoms developed among soybean cultivars in alignment with provided foliar resistance ratings; however, the F. virguliforme IC were not significantly different between SDS foliar symptomatic and asymptomatic cultivars. F. virguliforme was detected in all cultivars at the first sample collection point 25 DAP (V3 vegetative growth stage), and the IC increased throughout the season, peaking at the last sample collection point 153 DAP (postharvest). Collectively, appearance and disease severity ratings of SDS foliar symptoms were not associated with F. virguliforme quantity in roots, suggesting a need to include F. virguliforme root colonization in breeding efforts to screen soybean germplasm for F. virguliforme root infection resistance. The findings also demonstrates root colonization of the pathogen on nonsymptomatic soybean cultivars leading to persistence of the pathogen in the field, and possible hidden yield loss.

QTL mapping and epistatic interaction analysis of field resistance to sudden death syndrome (Fusarium virguliforme) in soybean

KEY MESSAGE

Two interactive quantitative trait loci (QTLs) controlled the field resistance to sudden death syndrome (SDS) in soybean. The interaction between them was confirmed. Sudden death syndrome (SDS), caused by Fusarium virguliforme, is a major disease of soybean [Glycine max (L.) Merr.] in the United States. Breeding for soybean resistance to SDS is the most cost-effective method to manage the disease. The objective of this study was to identify and characterize quantitative trait loci (QTLs) underlying field resistance to SDS in a recombinant inbred line population from the cross GD2422 × LD01-5907. This population was genotyped with 1786 polymorphic single nucleotide polymorphisms (SNPs) using SoySNP6 K iSelect BeadChip and evaluated for SDS resistance in a naturally infested field. Four SDS resistance QTLs were mapped on Chromosomes 4, 8, 12 and 18. The resistant parent, LD01-5907, contributed the resistance alleles for the QTLs on Chromosomes 8 and 18 (qSDS-8 and qSDS-18), while the other parent, GD2422, provided the resistance alleles for the QTLs on Chromosomes 4 and 12 (qSDS-4 and qSDS-12). The minor QTL on Chromosome 12 (qSDS-12) is novel. The QTL on Chromosomes 8 and 18 (qSDS-8 and qSDS-18) overlapped with two soybean cyst nematode resistance-related loci, Rhg4 and Rhg1, respectively. A significant interaction between qSDS-8 and qSDS-18 was detected by disease incidence. Individual effects together with the interaction effect explained around 70% of the phenotypic variance. The epistatic interaction of qSDS-8 and qSDS-18 was confirmed by the field performance across multiple years. Furthermore, the resistance alleles at qSDS-8 and qSDS-18 were demonstrated to be recessive. The SNP markers linked to these QTLs will be useful for marker-assisted breeding to enhance the SDS resistance.

Field evaluation of three sources of genetic resistance to sudden death syndrome of soybean

KEY MESSAGE

Despite numerous challenges, field testing of three sources of genetic resistance to sudden death syndrome of soybean provides information to more effectively improve resistance to this disease in cultivars. Sudden death syndrome (SDS) of soybean [Glycine max (L.) Merrill] is a disease that causes yield loss in soybean growing regions across the USA and worldwide. While several quantitative trait loci (QTL) for SDS resistance have been mapped, studies to further evaluate these QTL are limited. The objective of our research was to map SDS resistance QTL and to test the effect of mapped resistance QTL on foliar symptoms when incorporated into elite soybean backgrounds. We mapped a QTL from Ripley to chromosome 10 (CHR10) and a QTL from PI507531 to chromosomes 1 and 18 (CHR1 and 18). Six populations were then developed to test the following QTL: cqSDS-001, with resistance originating from PI567374, CHR10, CHR1, and CHR18. The populations which segregated for resistant and susceptible QTL alleles were field tested in multiple environments and evaluated for SDS foliar symptoms. While foliar disease development was variable across environments and populations, a significant effect of each QTL on disease was detected within at least one environment. This includes the detection of cqSDS-001 in three genetic backgrounds. The QTL allele from the resistant parents was associated with greater resistance than the susceptible alleles for all QTL and backgrounds with the exception of the allele for CHR18, where the opposite occurred. This study highlights the importance and difficulties of evaluating QTL and the need for multi-year SDS field testing. The information presented in this study can aid breeders in making decisions to improve resistance to SDS.

Mapping of new quantitative trait loci for sudden death syndrome and soybean cyst nematode resistance in two soybean populations

KEY MESSAGE

Novel QTL conferring resistance to both the SDS and SCN was detected in two RIL populations. Dual resistant RILs could be used in breeding programs for developing resistant soybean cultivars. Soybean cultivars, susceptible to the fungus Fusarium virguliforme, which causes sudden death syndrome (SDS), and to the soybean cyst nematode (SCN) (Heterodera glycines), suffer yield losses valued over a billion dollars annually. Both pathogens may occur in the same production fields. Planting of cultivars genetically resistant to both pathogens is considered one of the most effective means to control the two pathogens. The objective of the study was to map quantitative trait loci (QTL) underlying SDS and SCN resistances. Two recombinant inbred line (RIL) populations were developed by crossing 'A95-684043', a high-yielding maturity group (MG) II line resistant to SCN, with 'LS94-3207' and 'LS98-0582' of MG IV, resistant to both F. virguliforme and SCN. Two hundred F7 derived recombinant inbred lines from each population AX19286 (A95-684043 × LS94-3207) and AX19287 (A95-684043 × LS98-0582) were screened for resistance to each pathogen under greenhouse conditions. Five hundred and eighty and 371 SNP markers were used for mapping resistance QTL in each population. In AX19286, one novel SCN resistance QTL was mapped to chromosome 8. In AX19287, one novel SDS resistance QTL was mapped to chromosome 17 and one novel SCN resistance QTL was mapped to chromosome 11. Previously identified additional SDS and SCN resistance QTL were also detected in the study. Lines possessing superior resistance to both pathogens were also identified and could be used as germplasm sources for breeding SDS- and SCN-resistant soybean cultivars.

Integration of sudden death syndrome resistance loci in the soybean genome

KEY MESSAGE

Complexity and inconsistencies in resistance mapping publications of soybean sudden death syndrome (SDS) result in interpretation difficulty. This review integrates SDS mapping literature and proposes a new nomenclature system for reproducible SDS resistance loci. Soybean resistance to sudden death syndrome (SDS) is composed of foliar resistance to phytotoxins and root resistance to pathogen invasion. There are more than 80 quantitative trait loci (QTL) and dozens of single nucleotide polymorphisms (SNPs) associated with soybean resistance to SDS. The validity of these QTL and SNPs is questionable because of the complexity in phenotyping methodologies, the disease synergism between SDS and soybean cyst nematode (SCN), the variability from the interactions between soybean genotypes and environments, and the inconsistencies in the QTL nomenclature. This review organizes SDS mapping results and proposes the Rfv (resistance to Fusarium virguliforme) nomenclature based on supporting criteria described in the text. Among ten reproducible loci receiving our Rfv nomenclature, Rfv18-01 is mostly supported by field studies and it co-localizes to the SCN resistance locus rhg1. The possibility that Rfv18-01 is a pleiotropic resistance locus and the concern about Rfv18-01 being confounded with Rhg1 is discussed. On the other hand, Rfv06-01, Rfv06-02, Rfv09-01, Rfv13-01, and Rfv16-01 were identified both by screening soybean leaves against phytotoxic culture filtrates and by evaluating SDS severity in fields. Future phenotyping using leaf- and root-specific resistance screening methodologies may improve the precision of SDS resistance, and advanced genetic studies may further clarify the interactions among soybean genotypes, F. virguliforme, SCN, and environments. The review provides a summary of the SDS resistance literature and proposes a framework for communicating SDS resistance loci for future research considering molecular interactions and genetic breeding for soybean SDS resistance.

'MN1606SP' by 'Spencer' filial soybean population reveals novel quantitative trait loci and interactions among loci conditioning SDS resistance

KEY MESSAGE

Four novel QTL and interactions among QTL were identified in this research, using as a parent line the most SDS-resistant genotype within soybean cultivars of the US early maturity groups. Soybean sudden death syndrome (SDS) reduces soybean yield in most of the growing areas of the world. The causal agent of SDS, soilborne fungus Fusarium virguliforme (Fv), releases phytotoxins taken up by the plant to produce chlorosis and necrosis in the leaves. Planting resistant cultivars is the most successful management practice to control the disease. The objective of this study was to identify quantitative trait loci (QTL) associated with the resistance response of MN1606SP to SDS. A mapping population of F 2:3 lines created by crossing the highly resistant cultivar 'MN1606SP' and the susceptible cultivar 'Spencer' was phenotyped in the greenhouse at three different planting times, each with three replications. Plants were artificially inoculated using SDS infested sorghum homogeneously mixed with the soil. Data were collected on three disease criteria, foliar disease incidence (DI), foliar leaf scorch disease severity (DS), and root rot severity. Disease index (DX) was calculated as DI × DS. Ten QTL were identified for the different disease assessment criteria, three for DI, four for DX, and three for root rot severity. Three QTL identified for root rot severity and one QTL for disease incidence are considered novel, since no previous reports related to these QTL are available. Among QTL, two interactions were detected between four different QTL. The interactions suggest that resistance to SDS is not only dependent on additive gene effects. The novel QTL and the interactions observed in this study will be useful to soybean breeders for improvement of SDS resistance in soybean germplasm.

QTLomics in Soybean: A Way Forward for Translational Genomics and Breeding

Food legumes play an important role in attaining both food and nutritional security along with sustainable agricultural production for the well-being of humans globally. The various traits of economic importance in legume crops are complex and quantitative in nature, which are governed by quantitative trait loci (QTLs). Mapping of quantitative traits is a tedious and costly process, however, a large number of QTLs has been mapped in soybean for various traits albeit their utilization in breeding programmes is poorly reported. For their effective use in breeding programme it is imperative to narrow down the confidence interval of QTLs, to identify the underlying genes, and most importantly allelic characterization of these genes for identifying superior variants. In the field of functional genomics, especially in the identification and characterization of gene responsible for quantitative traits, soybean is far ahead from other legume crops. The availability of genic information about quantitative traits is more significant because it is easy and effective to identify homologs than identifying shared syntenic regions in other crop species. In soybean, genes underlying QTLs have been identified and functionally characterized for phosphorous efficiency, flowering and maturity, pod dehiscence, hard-seededness, α-Tocopherol content, soybean cyst nematode, sudden death syndrome, and salt tolerance. Candidate genes have also been identified for many other quantitative traits for which functional validation is required. Using the sequence information of identified genes from soybean, comparative genomic analysis of homologs in other legume crops could discover novel structural variants and useful alleles for functional marker development. The functional markers may be very useful for molecular breeding in soybean and harnessing benefit of translational research from soybean to other leguminous crops. Thus, soybean crop can act as a model crop for translational genomics and breeding of quantitative traits in legume crops. In this review, we summarize current status of identification and characterization of genes underlying QTLs for various quantitative traits in soybean and their significance in translational genomics and breeding of other legume crops.

Proteomic Profiling and the Predicted Interactome of Host Proteins in Compatible and Incompatible Interactions Between Soybean and Fusarium virguliforme

Sudden death syndrome (SDS) is a complex of two diseases of soybean (Glycine max), caused by the soil borne pathogenic fungus Fusarium virguliforme. The root rot and leaf scorch diseases both result in significant yield losses worldwide. Partial SDS resistance has been demonstrated in multiple soybean cultivars. This study aimed to highlight proteomic changes in soybean roots by identifying proteins which are differentially expressed in near isogenic lines (NILs) contrasting at the Rhg1/Rfs2 locus for partial resistance or susceptibility to SDS. Two-dimensional gel electrophoresis resolved approximately 1000 spots on each gel; 12 spots with a significant (P < 0.05) difference in abundance of 1.5-fold or more were picked, trypsin-digested, and analyzed using quadruple time-of-flight tandem mass spectrometry. Several spots contained more than one protein, so that 18 distinct proteins were identified overall. A functional analysis performed to categorize the proteins depicted that the major pathways altered by fungal infection include disease resistance, stress tolerance, and metabolism. This is the first report which identifies proteins whose abundances are altered in response to fungal infection leading to SDS. The results provide valuable information about SDS resistance in soybean plants, and plant partial resistance responses in general. More importantly, several of the identified proteins could be good candidates for the development of SDS-resistant soybean plants.

Four receptor-like cytoplasmic kinases regulate development and immunity in rice

Receptor-like cytoplasmic kinases (RLCKs) represent a large family of proteins in plants. However, few RLCKs have been well characterized. Here, we report the functional characterization of four rice RLCKs - OsRLCK57, OsRLCK107, OsRLCK118 and OsRLCK176 from subfamily VII. These OsRLCKs interact with the rice brassinosteroid receptor, OsBRI1 in yeast cell, but not the XA21 immune receptor. Transgenic lines silenced for each of these genes have enlarged leaf angles and are hypersensitive to brassinolide treatment compared to wild type rice. Transgenic plants silenced for OsRLCK57 had significantly fewer tillers and reduced panicle secondary branching, and lines silenced for OsRLCK107 and OsRLCK118 produce fewer seeds. Silencing of these genes decreased Xa21 gene expression and compromised XA21-mediated immunity to Xanthomonas oryzae pv. oryzae. Our study demonstrates that these OsRLCKs negatively regulate BR signalling, while positively regulating immune responses by contributing to the expression of the immune receptor XA21.

Genome-wide association and epistasis studies unravel the genetic architecture of sudden death syndrome resistance in soybean

Soybean [Glycine max (L.) Merr.] is an economically important crop that is grown worldwide. Sudden death syndrome (SDS), caused by Fusarium virguliforme, is one of the top yield-limiting diseases in soybean. However, the genetic basis of SDS resistance, especially with respect to epistatic interactions, is still unclear. To better understand the genetic architecture of soybean SDS resistance, genome-wide association and epistasis studies were performed using a population of 214 germplasm accessions and 31,914 SNPs from the SoySNP50K Illumina Infinium BeadChip. Twelve loci and 12 SNP-SNP interactions associated with SDS resistance were identified at various time points after inoculation. These additive and epistatic loci together explained 24-52% of the phenotypic variance. Disease-resistant, pathogenesis-related and chitin- and wound-responsive genes were identified in the proximity of peak SNPs, including stress-induced receptor-like kinase gene 1 (SIK1), which is pinpointed by a trait-associated SNP and encodes a leucine-rich repeat-containing protein. We report that the proportion of phenotypic variance explained by identified loci may be considerably improved by taking epistatic effects into account. This study shows the necessity of considering epistatic effects in soybean SDS resistance breeding using marker-assisted and genomic selection approaches. Based on our findings, we propose a model for soybean root defense against the SDS pathogen. Our results facilitate identification of the molecular mechanism underlying SDS resistance in soybean, and provide a genetic basis for improvement of soybean SDS resistance through breeding strategies based on additive and epistatic effects.

Genetic characteristics of soybean resistance to HG type 0 and HG type 1.2.3.5.7 of the cyst nematode analyzed by genome-wide association mapping

BACKGROUND

Soybean cyst nematode (SCN, Heterodera glycines Ichinohe) is one of the most fatal pests of soybean (Glycine max (L.) Merr.) worldwide and causes huge loss of soybean yield each year. Multiple sources of resistance are urgently needed for effective management of SCN via the development of resistant cultivars. The aim of the present study was to investigate the genetic architecture of resistance to SCN HG Type 0 (race 3) and HG Type 1.2.3.5.7 (race 4) in landraces and released elite soybean cultivars mostly from China.

RESULTS

A total of 440 diverse soybean landraces and elite cultivars were screened for resistance to SCN HG Type 0 and HG Type 1.2.3.5.7. Exactly 131 new sources of SCN resistance were identified. Lines were genotyped by SNP markers detected by the Specific Locus Amplified Fragment Sequencing (SLAF-seq) approach. A total of 36,976 SNPs were identified with minor allele frequencies (MAF) > 4% that were present in 97% of all the genotypes. Genome-wide association mapping showed that a total of 19 association signals were significantly related to the resistance for the two HG Types. Of the 19 association signals, eight signals overlapped with reported QTL including Rhg1 and Rhg4 genes. Another eight were located in the linked regions encompassing known QTL. Three QTL were found that were not previously reported. The average value of female index (FI) of soybean accessions with resistant alleles was significantly lower than those with susceptible alleles for each peak SNP. Disease resistance proteins with leucine rich regions, cytochrome P450s, protein kinases, zinc finger domain proteins, RING domain proteins, MYB and WRKY transcription activation families were identified. Such proteins may participate in the resistant reaction to SCN and were frequently found in the tightly linked genomic regions of the peak SNPs.

CONCLUSIONS

GWAS extended understanding of the genetic architecture of SCN resistance in multiple genetic backgrounds. Nineteen association signals were obtained for the resistance to the two Hg Types of SCN. The multiple beneficial alleles from resistant germplasm sources will be useful for the breeding of cultivars with improved resistance to SCN. Analysis of genes near association signals may facilitate the recognition of the causal gene(s) underlying SCN resistances.

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.

Soybean kinome: functional classification and gene expression patterns

The protein kinase (PK) gene family is one of the largest and most highly conserved gene families in plants and plays a role in nearly all biological functions. While a large number of genes have been predicted to encode PKs in soybean, a comprehensive functional classification and global analysis of expression patterns of this large gene family is lacking. In this study, we identified the entire soybean PK repertoire or kinome, which comprised 2166 putative PK genes, representing 4.67% of all soybean protein-coding genes. The soybean kinome was classified into 19 groups, 81 families, and 122 subfamilies. The receptor-like kinase (RLK) group was remarkably large, containing 1418 genes. Collinearity analysis indicated that whole-genome segmental duplication events may have played a key role in the expansion of the soybean kinome, whereas tandem duplications might have contributed to the expansion of specific subfamilies. Gene structure, subcellular localization prediction, and gene expression patterns indicated extensive functional divergence of PK subfamilies. Global gene expression analysis of soybean PK subfamilies revealed tissue- and stress-specific expression patterns, implying regulatory functions over a wide range of developmental and physiological processes. In addition, tissue and stress co-expression network analysis uncovered specific subfamilies with narrow or wide interconnected relationships, indicative of their association with particular or broad signalling pathways, respectively. Taken together, our analyses provide a foundation for further functional studies to reveal the biological and molecular functions of PKs in soybean.

Orthologous plant microRNAs: microregulators with great potential for improving stress tolerance in plants

Small RNAs that are highly conserved across many plant species are involved in stress responses. Plants are exposed to many types of unfavorable conditions during their life cycle that result in some degree of stress. Recent studies on microRNAs (miRNAs) have highlighted their great potential as regulators of stress tolerance in plants. One of the possible ways in which plants counter environmental stresses is by altering their gene expression by the action of miRNAs. miRNAs regulate the expression of target genes by hybridizing to their nascent reverse complementary sequences marking them for cleavage in the nucleus or translational repression in the cytoplasm. Some miRNAs have been reported to be key regulators in biotic as well as abiotic stress responses across many species. The present review highlights some of the regulatory roles of orthologous plant miRNAs in response to various types of stress conditions.

Genome-wide association mapping of quantitative resistance to sudden death syndrome in soybean

BACKGROUND

Sudden death syndrome (SDS) is a serious threat to soybean production that can be managed with host plant resistance. To dissect the genetic architecture of quantitative resistance to the disease in soybean, two independent association panels of elite soybean cultivars, consisting of 392 and 300 unique accessions, respectively, were evaluated for SDS resistance in multiple environments and years. The two association panels were genotyped with 52,041 and 5,361 single nucleotide polymorphisms (SNPs), respectively. Genome-wide association mapping was carried out using a mixed linear model that accounted for population structure and cryptic relatedness.

RESULT

A total of 20 loci underlying SDS resistance were identified in the two independent studies, including 7 loci localized in previously mapped QTL intervals and 13 novel loci. One strong peak of association on chromosome 18, associated with all disease assessment criteria across the two panels, spanned a physical region of 1.2 Mb around a previously cloned SDS resistance gene (GmRLK18-1) in locus Rfs2. An additional variant independently associated with SDS resistance was also found in this genomic region. Other peaks were within, or close to, sequences annotated as homologous to genes previously shown to be involved in plant disease resistance. The identified loci explained an average of 54.5% of the phenotypic variance measured by different disease assessment criteria.

CONCLUSIONS

This study identified multiple novel loci and refined the map locations of known loci related to SDS resistance. These insights into the genetic basis of SDS resistance can now be used to further enhance durable resistance to SDS in soybean. Additionally, the associations identified here provide a basis for further efforts to pinpoint causal variants and to clarify how the implicated genes affect SDS resistance in soybean.

Usefulness of 10 genomic regions in soybean associated with sudden death syndrome resistance

Sudden death syndrome (SDS) is an important soybean [Glycine max (L) Merrill] disease caused by the soilborne fungus Fusarium virguliforme. Currently, 14 quantitative trait loci (QTL) had been confirmed associated with resistance or tolerance to SDS. The objective of the study was to evaluate usefulness of 10 of these QTL in controlling disease expression. Six populations were developed providing a total of 321 F2-derived lines for the study. Recombinant inbred lines (RIL) used as parents were obtained from populations of 'Essex' × 'Forrest' (EF), 'Flyer' × 'Hartwig' (FH), and 'Pyramid' × 'Douglas' (PD). Disease resistance was evaluated in the greenhouse at three different planting times, each with four replications, using sorghum infested with F. virguliforme homogeneously mixed in the soil (Luckew et al., Crop Sci 52:2215-2223, 2012). Four disease assessment criteria-foliar disease incidence (DI), foliar leaf scorch disease severity (DS), area under the disease progress curve (AUDPC), and root rot severity-were used. QTL were identified in more than one of the disease assessment criteria, mainly associated with lines in the most resistant categories. Five QTL (qRfs4, qRfs5, qRfs7, qRfs12, and Rfs16) were associated with at least one of the disease assessments across multiple populations. Of the five, qRfs4 was associated with DI, AUDPC, and root rot severity, and Rfs16 with AUDPC and root rot severity. The findings suggest it may be possible for plant breeders to focus on stacking a subset of the previously identified QTL to improve resistance to SDS in soybean.

War of the worms: how plants fight underground attacks

Sedentary plant-parasitic nematodes (PPNs) establish specialized feeding cells within roots to maintain long-term relationships with their hosts. However, feeding cells degenerate prematurely in plants that harbor resistance (R) genes against these parasites reducing their life span and ability to reproduce. Recognition of the nematode, mediated directly or indirectly by plant R proteins, occurs via nematode secreted effectors and evokes a resistance response, which is referred to as effector-triggered immunity (ETI). Recent breakthroughs in nematode effector biology shed new light on key players mediating ETI and have identified those involved in plant defense suppression as novel targets for engineering resistance in transgenic plants. Advances in plant genetics and genomics has facilitated the discovery of R genes to nematodes. Nevertheless, underlying resistance mechanisms remain poorly understood and are confounded by recently identified R genes that do not fit previously proposed paradigms. Thus, there is still much to be learned about how plants fight off underground attacks from PPNs. In coming years, we can expect breakthroughs in our understanding of the nature and mechanisms of plant resistance and nematode virulence as we explore these novel R genes.

Homo-dimerization and ligand binding by the leucine-rich repeat domain at RHG1/RFS2 underlying resistance to two soybean pathogens

BACKGROUND

The protein encoded by GmRLK18-1 (Glyma_18_02680 on chromosome 18) was a receptor like kinase (RLK) encoded within the soybean (Glycine max L. Merr.) Rhg1/Rfs2 locus. The locus underlies resistance to the soybean cyst nematode (SCN) Heterodera glycines (I.) and causal agent of sudden death syndrome (SDS) Fusarium virguliforme (Aoki). Previously the leucine rich repeat (LRR) domain was expressed in Escherichia coli.

RESULTS

The aims here were to evaluate the LRRs ability to; homo-dimerize; bind larger proteins; and bind to small peptides. Western analysis suggested homo-dimers could form after protein extraction from roots. The purified LRR domain, from residue 131-485, was seen to form a mixture of monomers and homo-dimers in vitro. Cross-linking experiments in vitro showed the H274N region was close (<11.1 A) to the highly conserved cysteine residue C196 on the second homo-dimer subunit. Binding constants of 20-142 nM for peptides found in plant and nematode secretions were found. Effects on plant phenotypes including wilting, stem bending and resistance to infection by SCN were observed when roots were treated with 50 pM of the peptides. Far-Western analyses followed by MS showed methionine synthase and cyclophilin bound strongly to the LRR domain. A second LRR from GmRLK08-1 (Glyma_08_g11350) did not show these strong interactions.

CONCLUSIONS

The LRR domain of the GmRLK18-1 protein formed both a monomer and a homo-dimer. The LRR domain bound avidly to 4 different CLE peptides, a cyclophilin and a methionine synthase. The CLE peptides GmTGIF, GmCLE34, GmCLE3 and HgCLE were previously reported to be involved in root growth inhibition but here GmTGIF and HgCLE were shown to alter stem morphology and resistance to SCN. One of several models from homology and ab-initio modeling was partially validated by cross-linking. The effect of the 3 amino acid replacements present among RLK allotypes, A87V, Q115K and H274N were predicted to alter domain stability and function. Therefore, the LRR domain of GmRLK18-1 might underlie both root development and disease resistance in soybean and provide an avenue to develop new variants and ligands that might promote reduced losses to SCN.

Molecular mapping and breeding with microsatellite markers

In genetics databases for crop plant species across the world, there are thousands of mapped loci that underlie quantitative traits, oligogenic traits, and simple traits recognized by association mapping in populations. The number of loci will increase as new phenotypes are measured in more diverse genotypes and genetic maps based on saturating numbers of markers are developed. A period of locus reevaluation will decrease the number of important loci as those underlying mega-environmental effects are recognized. A second wave of reevaluation of loci will follow from developmental series analysis, especially for harvest traits like seed yield and composition. Breeding methods to properly use the accurate maps of QTL are being developed. New methods to map, fine map, and isolate the genes underlying the loci will be critical to future advances in crop biotechnology. Microsatellite markers are the most useful tool for breeders. They are codominant, abundant in all genomes, highly polymorphic so useful in many populations, and both economical and technically easy to use. The selective genotyping approaches, including genotype ranking (indexing) based on partial phenotype data combined with favorable allele data and bulked segregation event (segregant) analysis (BSA), will be increasingly important uses for microsatellites. Examples of the methods for developing and using microsatellites derived from genomic sequences are presented for monogenic, oligogenic, and polygenic traits. Examples of successful mapping, fine mapping, and gene isolation are given. When combined with high-throughput methods for genotyping and a genome sequence, the use of association mapping with microsatellite markers will provide critical advances in the analysis of crop traits.

Recombination suppression at the dominant Rhg1/Rfs2 locus underlying soybean resistance to the cyst nematode

Host resistance to "yellow dwarf" or "moonlight" disease cause by any population (Hg type) of Heterodera glycines I., the soybean cyst nematode (SCN), requires a functional allele at rhg1. The host resistance encoded appears to mimic an apoptotic response in the giant cells formed at the nematode feeding site about 24-48 h after nematode feeding commences. Little is known about how the host response to infection is mediated but a linked set of 3 genes has been identified within the rhg1 locus. This study aimed to identify the role of the genes within the locus that includes a receptor-like kinase (RLK), a laccase and an ion antiporter. Used were near isogeneic lines (NILs) that contrasted at their rhg1 alleles, gene-based markers, and a new Hg type 0 and new recombination events. A syntenic gene cluster on Lg B1 was found. The effectiveness of SNP probes from the RLK for distinguishing homolog sequence variants on LgB1 from alleles at the rhg1 locus on LgG was shown. The resistant allele of the rhg1 locus was shown to be dominant in NILs. None of the recombination events were within the cluster of the three candidate genes. Finally, rhg1 was shown to reduce the plant root development. A model for rhg1 as a dominant multi-gene resistance locus based on the developmental control was inferred.