Common loci underlie field resistance to soybean sudden death syndrome in Forrest, Pyramid, Essex, and Douglas

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.

Molecular Breeding to Overcome Biotic Stresses in Soybean: Update

Soybean (Glycine max (L.) Merr.) is an important leguminous crop and biotic stresses are a global concern for soybean growers. In recent decades, significant development has been carried outtowards identification of the diseases caused by pathogens, sources of resistance and determination of loci conferring resistance to different diseases on linkage maps of soybean. Host-plant resistance is generally accepted as the bestsolution because of its role in the management of environmental and economic conditions of farmers owing to low input in terms of chemicals. The main objectives of soybean crop improvement are based on the identification of sources of resistance or tolerance against various biotic as well as abiotic stresses and utilization of these sources for further hybridization and transgenic processes for development of new cultivars for stress management. The focus of the present review is to summarize genetic aspects of various diseases caused by pathogens in soybean and molecular breeding research work conducted to date.

Understanding Root Rot Disease in Agricultural Crops

Root rot diseases remain a major global threat to the productivity of agricultural crops. They are usually caused by more than one type of pathogen and are thus often referred to as a root rot complex. Fungal and oomycete species are the predominant participants in the complex, while bacteria and viruses are also known to cause root rot. Incorporating genetic resistance in cultivated crops is considered the most efficient and sustainable solution to counter root rot, however, resistance is often quantitative in nature. Several genetics studies in various crops have identified the quantitative trait loci associated with resistance. With access to whole genome sequences, the identity of the genes within the reported loci is becoming available. Several of the identified genes have been implicated in pathogen responses. However, it is becoming apparent that at the molecular level, each pathogen engages a unique set of proteins to either infest the host successfully or be defeated or contained in attempting so. In this review, a comprehensive summary of the genes and the potential mechanisms underlying resistance or susceptibility against the most investigated root rots of important agricultural crops is presented.

Fine mapping of a Phytophthora-resistance locus RpsGZ in soybean using genotyping-by-sequencing

BACKGROUND

Phytophthora root rot (PRR) caused by Phytophthora sojae (P. sojae) is one of the most serious limitations to soybean production worldwide. The identification of resistance gene(s) and their incorporation into elite varieties is an effective approach for breeding to prevent soybean from being harmed by this disease. A valuable mapping population of 228 F8:11 recombinant inbred lines (RILs) derived from a cross of the resistant cultivar Guizao1 and the susceptible cultivar BRSMG68 and a high-density genetic linkage map with an average distance of 0.81 centimorgans (cM) between adjacent bin markers in this population were used to map and explore candidate gene(s).

RESULTS

PRR resistance in Guizao1 was found to be controlled by a single Mendelian locus and was finely mapped to a 367.371-kb genomic region on chromosome 3 harbouring 19 genes, including 7 disease resistance (R)-like genes, in the reference Willliams 82 genome. Quantitative real-time PCR assays of possible candidate genes revealed that Glyma.03 g05300 was likely involved in PRR resistance.

CONCLUSIONS

These findings from the fine mapping of a novel Rps locus will serve as a basis for the cloning and transfer of resistance genes in soybean and the breeding of P. sojae-resistant soybean cultivars through marker-assisted selection.

Linkage Mapping for Foliar Necrosis of Soybean Sudden Death Syndrome

Sudden death syndrome (SDS) foliar symptoms consist of foliar chlorosis, foliar necrosis, leaf marginal curling, and premature defoliation, but resistance screening has been evaluated mostly based on the overall SDS foliar severity rather than on a specific foliar symptom. This study generated an F2 population derived from crossing the susceptible variety Sloan and the resistant germplasm line PI 243518, which exhibits resistance to both foliar chlorosis and necrosis. A total of 400 F2 lines were evaluated for foliar chlorosis, foliar necrosis, and overall SDS foliar symptoms, separately. Genotyping-by-sequencing was applied to obtain single nucleotide polymorphisms (SNPs) in the F2 population, and linkage mapping using 135 F2 lines with 969 high-quality SNPs identified a locus on chromosome 13 for foliar necrosis and SDS foliar symptoms. The locus partially overlaps with loci previously reported for SDS on chromosome 13, which is the third time the region from 15.98 to 21.00 Mbp has been reproduced independently and therefore qualifies this locus for a new nomenclature proposed as Rfv13-02. In summary, this study generated a new biparental population that enables not only the discovery of a locus for foliar necrosis and SDS foliar symptoms on chromosome 13 but also the potential for advanced exploration of SDS foliar resistance derived from the germplasm line PI 243518.

A high-density genetic map constructed using specific length amplified fragment (SLAF) sequencing and QTL mapping of seed-related traits in sesame (Sesamum indicum L.)

BACKGROUND

Sesame (Sesamum indicum L., 2n = 2x = 26) is an important oilseed crop with high oil content but small seed size. To reveal the genetic loci of the quantitative seed-related traits, we constructed a high-density single nucleotide polymorphism (SNP) linkage map of an F₂ population by using specific length amplified fragment (SLAF) technique and determined the quantitative trait loci (QTLs) of seed-related traits for sesame based on the phenotypes of F₃ progeny.

RESULTS

The genetic map comprised 2159 SNP markers distributed on 13 linkage groups (LGs) and was 2128.51 cM in length, with an average distance of 0.99 cM between adjacent markers. QTL mapping revealed 19 major-effect QTLs with the phenotypic effect (R²) more than 10%, i.e., eight QTLs for seed coat color, nine QTLs for seed size, and two QTLs for 1000-seed weight (TSW), using composite interval mapping method. Particularly, LG04 and LG11 contained collocated QTL regions for the seed coat color and seed size traits, respectively, based on their close or identical locations. In total, 155 candidate genes for seed coat color, 22 for seed size traits, and 54 for TSW were screened and analyzed.

CONCLUSIONS

This report presents the first QTL mapping of seed-related traits in sesame using an F₂ population. The results reveal the location of specific markers associated with seed-related traits in sesame and provide the basis for further seed quality traits research.

Legume genomics and transcriptomics: From classic breeding to modern technologies

Legumes are essential and play a significant role in maintaining food standards and augmenting physiochemical soil properties through the biological nitrogen fixation process. Biotic and abiotic factors are the main factors limiting legume production. Classical breeding methodologies have been explored extensively about the problem of truncated yield in legumes but have not succeeded at the desired rate. Conventional breeding improved legume genotypes but with more resources and time. Recently, the invention of next-generation sequencing (NGS) and high-throughput methods for genotyping have opened new avenues for research and developments in legume studies. During the last decade, genome sequencing for many legume crops documented. Sequencing and re-sequencing of important legume species have made structural variation and functional genomics conceivable. NGS and other molecular techniques such as the development of markers; genotyping; high density genetic linkage maps; quantitative trait loci (QTLs) identification, expressed sequence tags (ESTs), single nucleotide polymorphisms (SNPs); and transcription factors incorporated into existing breeding technologies have made possible the accurate and accelerated delivery of information for researchers. The application of genome sequencing, RNA sequencing (transcriptome sequencing), and DNA sequencing (re-sequencing) provide considerable insights for legume development and improvement programs.

Moreover, RNA-Seq helps to characterize genes, including differentially expressed genes, and can be applied for functional genomics studies, especially when there is limited information available for the studied genomes. Genome-based crop development studies and the availability of genomics data as well as decision-making gears look be specific for breeding programs. This review mainly presents an overview of the path from classical breeding to new emerging genomics tools, which will trigger and accelerate genomics-assisted breeding for recognition of novel genes for yield and quality characters for sustainable legume crop production.

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.

Genome wide association study identifies novel single nucleotide polymorphic loci and candidate genes involved in soybean sudden death syndrome resistance

Fusarium virguliforme is a soil borne root pathogen that causes sudden death syndrome (SDS) in soybean [Glycine max (L.) Merrill]. Once the fungus invades the root xylem tissues, the pathogen secretes toxins that cause chlorosis and necrosis in foliar tissues leading to defoliation, flower and pod drop and eventually death of plants. Resistance to F. virguliforme in soybean is partial and governed by over 80 quantitative trait loci (QTL). We have conducted genome-wide association study (GWAS) for a group of 254 plant introductions lines using a panel of approximately 30,000 SNPs and identified 19 single nucleotide polymorphic loci (SNPL) that are associated with 14 genomic regions encoding foliar SDS and eight SNPL associated with seven genomic regions for root rot resistance. Of the identified 27 SNPL, six SNPL for foliar SDS resistance and two SNPL for root rot resistance co-mapped to previously identified QTL for SDS resistance. This study identified 13 SNPL associated with eight novel genomic regions containing foliar SDS resistance genes and six SNPL with five novel regions for root-rot resistance. This study identified five genes carrying nonsynonymous mutations: (i) three of which mapped to previously identified QTL for foliar SDS resistance and (ii) two mapped to two novel regions containing root rot resistance genes. Of the three genes mapped to QTL for foliar SDS resistance genes, two encode LRR-receptors and third one encodes a novel protein with unknown function. Of the two genes governing root rot resistance, Glyma.01g222900.1 encodes a soybean-specific LEA protein and Glyma.10g058700.1 encodes a heparan-alpha-glucosaminide N-acetyltransferase. In the LEA protein, a conserved serine residue was substituted with asparagine; and in the heparan-alpha-glucosaminide N-acetyltransferase, a conserved histidine residue was substituted with an arginine residue. Such changes are expected to alter functions of these two proteins regulated through phosphorylation. The five genes with nonsynonymous mutations could be considered candidate SDS resistance genes and should be suitable molecular markers for breeding SDS resistance in soybean. The study also reports desirable plant introduction lines and novel genomic regions for enhancing SDS resistance in soybean.

Different loci associated with root and foliar resistance to sudden death syndrome (Fusarium virguliforme) in soybean

KEY MESSAGE

Different loci associated with root resistance to F. virguliforme colonization and foliar resistance to phytotoxin damage in soybean. Use of resistant cultivars is the most efficacious approach to manage soybean sudden death syndrome (SDS), caused by Fusarium virguliforme. The objectives of this study were to (1) map the loci associated with root and foliar resistance to F. virguliforme infection and (2) decipher the relationships between root infection, foliar damage, and plot yield. A mapping population consisting of 153 F4-derived recombinant inbred lines from the cross U01-390489 × E07080 was genotyped by SoySNP6 K BeadChip assay. Both foliar damage and F. virguliforme colonization in roots were investigated in the field, and a weak positive correlation was identified between them. Foliar damage had a stronger negative correlation with plot yield than F. virguliforme colonization. Twelve loci associated with foliar damage were identified, and four of them were associated with multiple traits across environments. In contrast, only one locus associated with root resistance to F. virguliforme colonization was identified and mapped on Chromosome 18. It colocalized with the locus associated with foliar damage in the same environment. The locus on Chromosome 6, qSDS6-2, and the locus on Chromosome 18, qSDS18-1, were associated with resistance to SDS phytotoxins and resistance to F. virguliforme colonization of roots, respectively. Both loci affected plot yield. Foliar damage-related traits, especially disease index, are valuable indicators for SDS resistance breeding because of consistency of the identified loci and their stronger correlation with plot yield. The information provided by this study will facilitate marker-assisted selection to improve SDS resistance in soybean.

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.

Tanscriptomic Study of the Soybean-Fusarium virguliforme Interaction Revealed a Novel Ankyrin-Repeat Containing Defense Gene, Expression of Whose during Infection Led to Enhanced Resistance to the Fungal Pathogen in Transgenic Soybean Plants

Fusarium virguliforme causes the serious disease sudden death syndrome (SDS) in soybean. Host resistance to this pathogen is partial and is encoded by a large number of quantitative trait loci, each conditioning small effects. Breeding SDS resistance is therefore challenging and identification of single-gene encoded novel resistance mechanisms is becoming a priority to fight this devastating this fungal pathogen. In this transcriptomic study we identified a few putative soybean defense genes, expression of which is suppressed during F. virguliforme infection. The F. virguliforme infection-suppressed genes were broadly classified into four major classes. The steady state transcript levels of many of these genes were suppressed to undetectable levels immediately following F. virguliforme infection. One of these classes contains two novel genes encoding ankyrin repeat-containing proteins. Expression of one of these genes, GmARP1, during F. virguliforme infection enhances SDS resistance among the transgenic soybean plants. Our data suggest that GmARP1 is a novel defense gene and the pathogen presumably suppress its expression to establish compatible interaction.

Genome-wide association mapping of partial resistance to Phytophthora sojae in soybean plant introductions from the Republic of Korea

BACKGROUND

Phytophthora root and stem rot is one of the most yield-limiting diseases of soybean [Glycine max (L.) Merr], caused by the oomycete Phytophthora sojae. Partial resistance is controlled by several genes and, compared to single gene (Rps gene) resistance to P. sojae, places less selection pressure on P. sojae populations. Thus, partial resistance provides a more durable resistance against the pathogen. In previous work, plant introductions (PIs) originating from the Republic of Korea (S. Korea) have shown to be excellent sources for high levels of partial resistance against P. sojae.

RESULTS

Resistance to two highly virulent P. sojae isolates was assessed in 1395 PIs from S. Korea via a greenhouse layer test. Lines exhibiting possible Rps gene immunity or rot due to other pathogens were removed and the remaining 800 lines were used to identify regions of quantitative resistance using genome-wide association mapping. Sixteen SNP markers on chromosomes 3, 13 and 19 were significantly associated with partial resistance to P. sojae and were grouped into seven quantitative trait loci (QTL) by linkage disequilibrium blocks. Two QTL on chromosome 3 and three QTL on chromosome 19 represent possible novel loci for partial resistance to P. sojae. While candidate genes at QTL varied in their predicted functions, the coincidence of QTLs 3-2 and 13-1 on chromosomes 3 and 13, respectively, with Rps genes and resistance gene analogs provided support for the hypothesized mechanism of partial resistance involving weak R-genes.

CONCLUSIONS

QTL contributing to partial resistance towards P. sojae in soybean germplasm originating from S. Korea were identified. The QTL identified in this study coincide with previously reported QTL, Rps genes, as well as novel loci for partial resistance. Molecular markers associated with these QTL can be used in the marker-assisted introgression of these alleles into elite cultivars. Annotations of genes within QTL allow hypotheses on the possible mechanisms of partial resistance to P. sojae.

Quantitative trait loci underlying host responses of soybean to Fusarium virguliforme toxins that cause foliar sudden death syndrome

KEY MESSAGE

Soybean deploys multiple genetic mechanisms to confer tolerance to Fusarium virguliforme toxins. This study revealed that F. virguliforme culture filtrates could be used in mapping QTL underlying foliar SDS resistance. Sudden death syndrome (SDS) is a major soybean disease throughout most of the soybean growing regions in the world including the United States. The disease is caused by the fungal pathogen, Fusarium virguliforme (Fv). The fungus produces several toxins that are responsible for development of interveinal leaf chlorosis and necrosis, which are typical foliar SDS symptoms. Growing of resistant cultivars has been the most effective method in controlling the disease. The objective of the present study was to identify quantitative trait loci (QTL) underlying host responses of soybean to Fv toxins present in culture filtrates. To accomplish this objective, two recombinant inbred line (RIL) populations, AX19286 (A95-684043 × LS94-3207) and AX19287 (A95-684043 × LS98-0582), segregating for SDS resistance were evaluated for foliar symptom development by applying two screening protocols, the stem cutting and the root feeding assays. The AX19286 population revealed two major and seven minor QTL for SDS resistance. In the AX19287 population, we identified five major QTL and three minor QTL. The two QTL mapped to Chromosome 7 [molecular linkage group (MLG) M] and Chromosome 20 (MLG I) are most likely novel, and were detected through screening of the AX19287 population with stem cutting and root feeding assays, respectively. This study established that Fv culture filtrates could be employed in mapping QTL underlying foliar SDS resistance. The outcomes of the research also suggest that multiple genetic mechanisms might be used by soybean to overcome the toxic effects of the toxins secreted by the pathogen into culture filtrates.

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.

Association mapping and genomic prediction for resistance to sudden death syndrome in early maturing soybean germplasm

Sudden death syndrome (SDS), caused by Fusarium virguliforme, has spread to northern soybean growing regions in the US causing significant yield losses. The objectives of this study were to identify loci underlying variation in plant responses to SDS through association mapping (AM) and to assess prediction accuracy of genomic selection (GS) in a panel of early maturing soybean germplasm. A set of 282 soybean breeding lines was selected from the University of Minnesota soybean breeding program and then genotyped using a genome-wide panel of 1536 single-nucleotide polymorphism markers. Four resistance traits, root lesion severity (RLS), foliar symptom severity (FSS), root retention (RR), and dry matter reduction (DMR), were evaluated using soil inoculation in the greenhouse. AM identified significant peaks in genomic regions of known SDS resistance quantitative trait loci cqSDS001, cqRfs4, and SDS11-2. Additionally, two novel loci, one on chromosome 3 and another on chromosome 18, were tentatively identified. A ninefold cross-validation scheme was used to assess the prediction accuracy of GS for SDS resistance. The prediction accuracy of single-trait GS (ST-GS) was 0.64 for RLS, but less than 0.30 for RR, DMR, and FSS. Compared to ST-GS, none of multi-trait GS (MT-GS) models significantly improved the prediction accuracy due to weak correlations between the four traits. This study suggests both AM and GS hold promise for implementation in genetic improvement of SDS resistance in existing soybean breeding programs.

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.

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.

QTLs for resistance to Phomopsis seed decay are associated with days to maturity in soybean (Glycine max)

Phomopsis seed decay (PSD), primarily caused by Phomopsis longicolla, is a major contributor to poor soybean seed quality and significant yield loss, particularly in early maturing soybean genotypes. However, it is not yet known whether PSD resistance is associated with early maturity. This study was conducted to identify quantitative trait loci (QTLs) for resistance to PSD and days to maturity using a recombinant inbred line (RIL) population derived from a cross between the PSD-resistant Taekwangkong and the PSD-susceptible SS2-2. Based on a genetic linkage map incorporating 117 simple sequence repeat markers, QTL analysis revealed two and three QTLs conferring PSD resistance and days to maturity, respectively, in the RIL population. Two QTLs (PSD-6-1 and PSD-10-2) for PSD resistance were identified in the intervals of Satt100-Satt460 and Sat_038-Satt243 on chromosomes 6 and 10, respectively. Two QTLs explained phenotypic variances in PSD resistance of 46.3 and 14.1 %, respectively. At the PSD-6-1 QTL, the PSD-resistant cultivar Taekwangkong contributed the allele with negative effect decreasing the infection rate of PSD and this QTL does not overlap with any previously reported loci for PSD resistance in other soybean genotypes. Among the three QTLs for days to maturity, two (Mat-6-2 and Mat-10-3) were located at positions similar to the PSD-resistance QTLs. The identification of the QTLs linked to both PSD resistance and days to maturity indicates a biological correlation between these two traits. The newly identified QTL for resistance to PSD associated with days to maturity in Taekwangkong will help improve soybean resistance to P. longicolla.

Genetic analysis and QTL mapping of seed coat color in sesame (Sesamum indicum L.)

Seed coat color is an important agronomic trait in sesame, as it is associated with seed biochemical properties, antioxidant content and activity and even disease resistance of sesame. Here, using a high-density linkage map, we analyzed genetic segregation and quantitative trait loci (QTL) for sesame seed coat color in six generations (P1, P2, F1, BC1, BC2 and F2). Results showed that two major genes with additive-dominant-epistatic effects and polygenes with additive-dominant-epistatic effects were responsible for controlling the seed coat color trait. Average heritability of the major genes in the BC1, BC2 and F2 populations was 89.30%, 24.00%, and 91.11% respectively, while the heritability of polygenes was low in the BC1 (5.43%), in BC2 (0.00%) and in F2 (0.89%) populations. A high-density map was constructed using 724 polymorphic markers. 653 SSR, AFLP and RSAMPL loci were anchored in 14 linkage groups (LG) spanning a total of 1,216.00 cM. The average length of each LG was 86.86 cM and the marker density was 1.86 cM per marker interval. Four QTLs for seed coat color, QTL1-1, QTL11-1, QTL11-2 and QTL13-1, whose heritability ranged from 59.33%-69.89%, were detected in F3 populations using CIM and MCIM methods. Alleles at all QTLs from the black-seeded parent tended to increase the seed coat color. Results from QTLs mapping and classical genetic analysis among the P1, P2, F1, BC1, BC2 and F2 populations were comparatively consistent. This first QTL analysis and high-density genetic linkage map for sesame provided a good foundation for further research on sesame genetics and molecular marker-assisted selection (MAS).

Effect of Fusarium virguliforme phytotoxin on soybean gene expression suggests a role in multidimensional defence

Sudden death syndrome (SDS), caused by Fusarium virguliforme, is an important yield-limiting disease of soybean. This soil-borne fungus colonizes soybean roots causing root rot, and also releases a phytotoxin that is translocated to leaves causing interveinal chlorosis and necrosis leading to defoliation and early maturation. The objective of our study was to compare gene expression profiles during the early response of soybean leaves exposed to sterile culture filtrates of F. virguliforme in soybean genotypes with different levels of resistance to SDS. The analysis identified SDS-related defence genes that were induced in the most resistant genotype, but not in the other genotypes. Further functional annotations based on sequence homology suggested that some of the induced genes probably encode proteins involved in cell wall modification, detoxification, defence responses, primary metabolism and membrane transport. Quantitative real-time reverse-transcribed polymerase chain reaction confirmed the differential transcript accumulation of a subset of these genes. In addition, in silico mapping of differentially expressed genes to SDS-resistant quantitative trait loci allowed for the identification of new potential defence genes that could be genetically mapped to the soybean genome, and could be used further in a marker-assisted selection programme. A comparison of the response of soybean to F. virguliforme phytotoxin (Fv toxin) relative to other biotic and abiotic stresses revealed that the resistance response to Fv toxin is quite similar to the response to inoculation with an incompatible Pseudomonas syringae pv. glycinea strain, suggesting that Fv toxin might induce hypersensitive response pathways in soybean leaf tissues in the absence of pathogen in these tissues.

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.

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

BACKGROUND

Soybean (Glycine max (L. Merr.)) resistance to any population of Heterodera glycines (I.), or Fusarium virguliforme (Akoi, O'Donnell, Homma & Lattanzi) required a functional allele at Rhg1/Rfs2. H. glycines, the soybean cyst nematode (SCN) was an ancient, endemic, pest of soybean whereas F. virguliforme causal agent of sudden death syndrome (SDS), was a recent, regional, pest. This study examined the role of a receptor like kinase (RLK) GmRLK18-1 (gene model Glyma_18_02680 at 1,071 kbp on chromosome 18 of the genome sequence) within the Rhg1/Rfs2 locus in causing resistance to SCN and SDS.

RESULTS

A BAC (B73p06) encompassing the Rhg1/Rfs2 locus was sequenced from a resistant cultivar and compared to the sequences of two susceptible cultivars from which 800 SNPs were found. Sequence alignments inferred that the resistance allele was an introgressed region of about 59 kbp at the center of which the GmRLK18-1 was the most polymorphic gene and encoded protein. Analyses were made of plants that were either heterozygous at, or transgenic (and so hemizygous at a new location) with, the resistance allele of GmRLK18-1. Those plants infested with either H. glycines or F. virguliforme showed that the allele for resistance was dominant. In the absence of Rhg4 the GmRLK18-1 was sufficient to confer nearly complete resistance to both root and leaf symptoms of SDS caused by F. virguliforme and provided partial resistance to three different populations of nematodes (mature female cysts were reduced by 30-50%). In the presence of Rhg4 the plants with the transgene were nearly classed as fully resistant to SCN (females reduced to 11% of the susceptible control) as well as SDS. A reduction in the rate of early seedling root development was also shown to be caused by the resistance allele of the GmRLK18-1. Field trials of transgenic plants showed an increase in foliar susceptibility to insect herbivory.

CONCLUSIONS

The inference that soybean has adapted part of an existing pathogen recognition and defense cascade (H.glycines; SCN and insect herbivory) to a new pathogen (F. virguliforme; SDS) has broad implications for crop improvement. Stable resistance to many pathogens might be achieved by manipulation the genes encoding a small number of pathogen recognition proteins.

Transcriptional analysis of soybean root response to Fusarium virguliforme, the causal agent of sudden death syndrome

Sudden death syndrome (SDS) of soybean can be caused by any of four distinct Fusarium species, with F. virguliforme and F. tucumaniae being the main casual agents in North and South America, respectively. Although the fungal tissue is largely confined to the roots, the fungus releases a toxin that is translocated to leaf tissues, in which it causes interveinal chlorosis and necrosis leading to scorching symptoms and possible defoliation. In this study, we report on an Affymetrix analysis measuring transcript abundances in resistant (PI 567.374) and susceptible (Essex) roots upon infection by F. virguliforme, 5 and 7 days postinoculation. Many of the genes with increased expression were common between resistant and susceptible plants (including genes related to programmed cell death, the phenylpropanoid pathway, defense, signal transduction, and transcription factors), but some genotype-specific expression was noted. Changes in small (sm)RNA levels between inoculated and mock-treated samples were also studied and implicate a role for these molecules in this interaction. In total, 2,467 genes were significantly changing in the experiment, with 1,694 changing in response to the pathogen; 93 smRNA and 42 microRNA that have putative soybean gene targets were identified from infected tissue. Comparing genotypes, 247 genes were uniquely modulating in the resistant host, whereas 378 genes were uniquely modulating in the susceptible host. Comparing locations of differentially expressed genes to known resistant quantitative trait loci as well as identifying smRNA that increased while their putative targets decreased (or vice versa) allowed for the narrowing of candidate SDS defense-associated genes.

The genetic control of tolerance to aluminum toxicity in the 'Essex' by 'Forrest' recombinant inbred line population

Aluminum (Al) toxicity to plant roots is a major problem of acidic soils. The main chemical reaction involved is Al hydrolysis. Application of lime or nitrate fertilizers to raise soil pH reduces Al toxicity but not as economically as a plant genotypes with natural tolerance against this stress. Ammonium fertilization of crops and assimilation of ammonium (even that derived from dinitrogen) are particularly acidifying of the root zone. The aims of the present study were to find genotypes of soybean tolerant to aluminum stress and identify QTL underlying that trait. Used were recombinant inbred lines (RILs) derived from the cross of 'Essex' by 'Forrest'. RILs were grown in a greenhouse for 3 weeks and then transferred to hydroponics in a growth chamber. Root lengths (RL) were measured before and 72 h after Al treatment. RL before and after Al treatment were measured and used to calculate root tolerance index (RTI) and relative mean growth (RMG). RILs 1, 85, 40 and 83 had significant (P<0.005) tolerance to Al stress judged by RL after Al, RTI and RMG. Eleven minor but significant marker-trait associations (P<0.05) were detected using one-way ANOVA but only two major loci were significant in composite interval maps (LOD>3.0). The QTL on linkage group F (chromosome 13) was in the interval Satt160-Satt252 with a peak at 24 cM (peak LOD was 3.3). The QTL underlay 31% of trait variation and the Essex allele provided an additional 1.61 cm of root growth over 72 h in the presence of Al. The QTL on linkage group C2 (probably chromosome 4) was in the interval from Satt202 to Satt371 with a peak at 3.2 cM (peak LOD was 14.7). The QTL underlay 34% of trait variation or 1.81 cm of growth over 72 h in the presence of Al. Both loci encompassed genes implicated in citrate metabolism, a method of aluminum detoxification known to vary among soybean cultivars. Two major loci and at least nine minor loci were inferred to underlie tolerance to Al. RILs and markers may be used to select alleles that increase tolerance to soybean against Al stress.

Pyramided QTL underlying tolerance to Phytophthora root rot in mega-environments from soybean cultivars 'Conrad' and 'Hefeng 25'

Phytophthora root rot (PRR) of soybean (Glycine max (L.) Merr.) is the second most important cause of yield loss by disease in North America, surpassed only by soybean cyst nematode (Wrather et al. in Can J Plant Pathol 23:115-121, 2001). Tolerance can provide economically useful disease control, conditioning partial resistance of soybean to PRR. The aims of this study were to identify new quantitative trait loci (QTL) underlying tolerance to PRR, and to evaluate the effects of pyramided or stacked loci on the level of tolerance. A North American cultivar 'Conrad' (tolerant to PRR) was crossed with a northeastern China cultivar 'Hefeng 25' (tolerant to PRR). Through single-seed descent, 140 F2:5 and F2:6 recombinant inbred lines were advanced. A total of 164 simple sequence repeat (SSR) markers were used to construct a genetic linkage map. The percentage of seedling death was measured over 2 years (2007 and 2008) in the field at four naturally infested locations in Canada and China following additional soil infestation and in the greenhouse following inoculation with Phytophthora sojae isolate. A total of eight QTL underlying tolerance to PRR were identified, located in five linkage groups (F, D1b+w, A2, B1, and C2). The phenotypic variation contributed by the loci ranged from 4.24 to 27.98%. QPRR-1 (anchored in the interval of SSR markers Satt325 and Satt343 of LG F), QPRR-2 (anchored in the interval of Satt005 and Satt600 of LG D1b+w), and QPRR-3 (anchored in the interval of Satt579 and Sat_089 of LG D1b+w) derived their beneficial allele from 'Conrad'. They were located at chromosomal locations known to underlie PRR tolerance in diverse germplasm. Five QTL that derived beneficial alleles from 'Hefeng 25' were identified. The QTL (QPRR-1 to QPRR-7) that were detected across at least three environments were selected for loci stacking and to analyze the relationship between number of tolerance loci and disease loss percentage. The accumulation of tolerance loci was positively correlated with decreases in disease loss percentage. The pyramid of loci underlying tolerance to PRR provided germplasm useful for crop improvement by marker-assisted selection and may provide durable cultivar tolerance against the PRR disease.

Iso-lines and inbred-lines confirmed loci that underlie resistance from cultivar 'Hartwig' to three soybean cyst nematode populations

Soybean [Glycine max (L.) Merr.] cultivars varied in their resistance to different populations of the soybean cyst nematode (SCN), Heterodera glycines, called HG Types. The rhg1 locus on linkage group G was necessary for resistance to all HG types. However, the loci for resistance to H. glycines HG Type 1.3- (race 14) and HG Type 1.2.5- (race 2) of the soybean cyst nematode have varied in their reported locations. The aims were to compare the inheritance of resistance to three nematode HG Types in a population segregating for resistance to SCN and to identify the underlying quantitative trait loci (QTL). 'Hartwig', a soybean cultivar resistant to most SCN HG Types, was crossed with the susceptible cultivar 'Flyer'. A total of 92 F5-derived recombinant inbred lines (RILs; or inbred lines) and 144 molecular markers were used for map development. The rhg1 associated QTL found in earlier studies were confirmed and shown to underlie resistance to all three HG Types in RILs (Satt309; HG Type 0, P = 0.0001 R (2) = 22%; Satt275; HG Type 1.3, P = 0.001, R (2) = 14%) and near isogeneic lines (NILs; or iso-lines; Satt309; HG Type 1.2.5-, P = 0.001 R (2) = 24%). A new QTL underlying resistance to HG Type 1.2.5- was detected on LG D2 (Satt574; P = 0.001, R (2) = 11%) among 14 RILs resistant to the other HG types. The locus was confirmed in a small NIL population consisting of 60 plants of ten genotypes (P = 0.04). This QTL (cqSCN-005) is located in an interval previously associated with resistance to both SDS leaf scorch from 'Pyramid' and 'Ripley' (cqSDS-001) and SCN HG Type 1.3- from Hartwig and Pyramid. The QTL detected will allow marker assisted selection for multigenic resistance to complex nematode populations in combination with sudden death syndrome resistance (SDS) and other agronomic traits.

Re-annotation of the physical map of Glycine max for polyploid-like regions by BAC end sequence driven whole genome shotgun read assembly

BACKGROUND

Many of the world's most important food crops have either polyploid genomes or homeologous regions derived from segmental shuffling following polyploid formation. The soybean (Glycine max) genome has been shown to be composed of approximately four thousand short interspersed homeologous regions with 1, 2 or 4 copies per haploid genome by RFLP analysis, microsatellite anchors to BACs and by contigs formed from BAC fingerprints. Despite these similar regions,, the genome has been sequenced by whole genome shotgun sequence (WGS). Here the aim was to use BAC end sequences (BES) derived from three minimum tile paths (MTP) to examine the extent and homogeneity of polyploid-like regions within contigs and the extent of correlation between the polyploid-like regions inferred from fingerprinting and the polyploid-like sequences inferred from WGS matches.

RESULTS

Results show that when sequence divergence was 1-10%, the copy number of homeologous regions could be identified from sequence variation in WGS reads overlapping BES. Homeolog sequence variants (HSVs) were single nucleotide polymorphisms (SNPs; 89%) and single nucleotide indels (SNIs 10%). Larger indels were rare but present (1%). Simulations that had predicted fingerprints of homeologous regions could be separated when divergence exceeded 2% were shown to be false. We show that a 5-10% sequence divergence is necessary to separate homeologs by fingerprinting. BES compared to WGS traces showed polyploid-like regions with less than 1% sequence divergence exist at 2.3% of the locations assayed.

CONCLUSION

The use of HSVs like SNPs and SNIs to characterize BACs wil improve contig building methods. The implications for bioinformatic and functional annotation of polyploid and paleopolyploid genomes show that a combined approach of BAC fingerprint based physical maps, WGS sequence and HSV-based partitioning of BAC clones from homeologous regions to separate contigs will allow reliable de-convolution and positioning of sequence scaffolds (see BES_scaffolds section of SoyGD). This approach will assist genome annotation for paleopolyploid and true polyploid genomes such as soybean and many important cereal and fruit crops.

Separate loci underlie resistance to root infection and leaf scorch during soybean sudden death syndrome

Soybean [Glycine max (L.) Merr.] cultivars show differences in their resistance to both the leaf scorch and root rot of sudden death syndrome (SDS). The syndrome is caused by root colonization by Fusarium virguliforme (ex. F. solani f. sp. glycines). Root susceptibility combined with reduced leaf scorch resistance has been associated with resistance to Heterodera glycines HG Type 1.3.6.7 (race 14) of the soybean cyst nematode (SCN). In contrast, the rhg1 locus underlying resistance to Hg Type 0 was found clustered with three loci for resistance to SDS leaf scorch and one for root infection. The aims of this study were to compare the inheritance of resistance to leaf scorch and root infection in a population that segregated for resistance to SCN and to identify the underlying quantitative trait loci (QTL). "Hartwig", a cultivar partially resistant to SDS leaf scorch, F. virguliforme root infection and SCN HG Type 1.3.6.7 was crossed with the partially susceptible cultivar "Flyer". Ninety-two F5-derived recombinant inbred lines and 144 markers were used for map development. Four QTL found in earlier studies were confirmed. One contributed resistance to leaf scorch on linkage group (LG) C2 (Satt277; P = 0.004, R2 = 15%). Two on LG G underlay root infection at R8 (Satt038; P = 0.0001 R2 = 28.1%; Satt115; P = 0.003, R2 = 12.9%). The marker Satt038 was linked to rhg1 underlying resistance to SCN Hg Type 0. The fourth QTL was on LG D2 underlying resistance to root infection at R6 (Satt574; P = 0.001, R2 = 10%). That QTL was in an interval previously associated with resistance to both SDS leaf scorch and SCN Hg Type 1.3.6.7. The QTL showed repulsion linkage with resistance to SCN that may explain the relative susceptibility to SDS of some SCN resistant cultivars. One additional QTL was discovered on LG G underlying resistance to SDS leaf scorch measured by disease index (Satt130; P = 0.003, R2 = 13%). The loci and markers will provide tagged alleles with which to improve the breeding of cultivars combining resistances to SDS leaf scorch, root infection and SCN HG Type 1.3.6.7.

Microsatellite discovery from BAC end sequences and genetic mapping to anchor the soybean physical and genetic maps

Whole-genome sequencing of the soybean (Glycine max (L.) Merr. 'Williams 82') has made it important to integrate its physical and genetic maps. To facilitate this integration of maps, we screened 3290 microsatellites (SSRs) identified from BAC end sequences of clones comprising the 'Williams 82' physical map. SSRs were screened against 3 mapping populations. We found the AAT and ACT motifs produced the greatest frequency of length polymorphisms, ranging from 17.2% to 32.3% and from 11.8% to 33.3%, respectively. Other useful motifs include the dinucleotide repeats AG, AT, and AG, with frequency of length polymorphisms ranging from 11.2% to 18.4% (AT), 12.4% to 20.6% (AG), and 11.3% to 16.4% (GT). Repeat lengths less than 16 bp were generally less useful than repeat lengths of 40-60 bp. Two hundred and sixty-five SSRs were genetically mapped in at least one population. Of the 265 mapped SSRs, 60 came from BAC singletons not yet placed into contigs of the physical map. One hundred and ten originated in BACs located in contigs for which no genetic map location was previously known. Ninety-five SSRs came from BACs within contigs for which one or more other BACs had already been mapped. For these fingerprinted contigs (FPC) a high percentage of the mapped markers showed inconsistent map locations. A strategy is introduced by which physical and genetic map inconsistencies can be resolved using the preliminary 4x assembly of the whole genome sequence of soybean.

Sexual reproduction in the soybean sudden death syndrome pathogen Fusarium tucumaniae

We investigated the sexual reproductive mode of the two most important etiological agents of soybean sudden death syndrome, Fusarium tucumaniae and Fusarium virguliforme. F. tucumaniae sexual crosses often were highly fertile, making it possible to assign mating type and assess female fertility in 24 South American isolates. These crosses produced red perithecia and oblong-elliptical ascospores, as is typical for sexual members of the F. solani species complex. Genotyping of progeny from three F. tucumaniae crosses confirmed that sexual recombination had occurred. In contrast, pairings among 17 U.S. F. virguliforme isolates never produced perithecia. Inter-species crosses between F. tucumaniae and F. virguliforme, in which infertile perithecia were induced only in one of the two F. tucumaniae mating types, suggest that all U.S. F. virguliforme isolates are of a single mating type. We conclude that the F. tucumaniae life cycle in S. America includes a sexual reproductive mode, and thus this species has greater potential for rapid evolution than the F. virguliforme population in the U.S., which may be exclusively asexual.

The development of BAC-end sequence-based microsatellite markers and placement in the physical and genetic maps of soybean

The composite map of soybean shared among Soybase, LIS and SoyGD (March 2006) contained 3,073 DNA markers in the "Locus" class. Among the markers were 1,019 class I microsatellite markers with 2-3 bp simple sequence repeats (SSRs) of >10 iterations (BARC-SSR markers). However, there were few class II SSRs (2-5 bp repeats with <10 iterations; mostly SIUC-Satt markers). The aims here were to increase the number of classes I and II SSR markers and to integrate bacterial artificial chromosome (BAC) clones onto the soybean physical map using the markers. Used was 10 Mb of BAC-end sequence (BES) derived from 13,473 reads from 7,050 clones constituting minimum tile path 2 of the soybean physical map ( http://www.soybeangenome.siu.edu ; SoyGD). Identified were 1,053 1-6 bp motif, repeat sequences, 333 from class I (>10 repeats) and 720 from class II (<10 repeats). Potential markers were shown on the MTP_SSR track at Gbrowse. Primers were designed as 20-24 bp oligomers that had Tm of 55 +/- 1 C that would generate 100-500 bp amplicons. About 853 useful primer pairs were established. Motifs were not randomly distributed with biases toward AT rich motifs. Strong biases against the GC motif and all tetra-nucleotide repeats were found. The markers discovered were useful. Among the first 135 targeted for use in genetic map improvement about 60% of class II markers and 75% of class I markers were polymorphic among on the parents of four recombinant inbred line (RIL) populations. Many of the BES-based SSRs were located on the soybean genetic map in regions with few BARC-SSR markers. Therefore, BES-based SSRs represent useful tools for genetic map development in soybean. New members of a consortium to map the markers in additional populations are invited.

Leveraging model legume information to find candidate genes for soybean sudden death syndrome using the legume information system

Comparative genomics is an emerging and powerful approach to achieve crop improvement. Using comparative genomics, information from model plant species can accelerate the discovery of genes responsible for disease and pest resistance, tolerance to plant stresses such as drought, and enhanced nutritional value including production of anti-oxidants and anti-cancer compounds. We demonstrate here how to use the Legume Information System for a comparative genomics study, leveraging genomic information from Medicago truncatula (barrel medic), the model legume, to find candidate genes involved with sudden death syndrome (SDS) in Glycine max (soybean). Specifically, genetic maps, physical maps, and annotated tentative consensus and expressed sequence tag (EST) sequences from G. max and M. truncatula can be compared. In addition, the recently published M. truncatula genomic sequences can be used to identify M. truncatula candidate genes in a genomic region syntenic to a quantitative trait loci region for SDS in soybean. Genomic sequences of candidate genes from M. truncatula can then be used to identify ESTs with sequence similarities from soybean for primer design and cloning of potential soybean disease causing alleles.

An updated 'Essex' by 'Forrest' linkage map and first composite interval map of QTL underlying six soybean traits

DNA marker maps based on single populations are the basis for gene, loci and genomic analyses. Individual maps can be integrated to produce composite maps with higher marker densities if shared marker orders are consistent. However, estimates of marker order in composite maps must include sets of markers that were not polymorphic in multiple populations. Often some of the pooled markers were not codominant, or were not correctly scored. The soybean composite map was composed of data from five separate populations based on northern US germplasm but does not yet include 'Essex' by 'Forrest' recombinant inbred line (RIL) population (E x F) or any southern US soybean cultivars. The objectives were, to update the E x F map with codominant markers, to compare marker orders among this map, the Forrest physical map and the composite soybean map and to compare QTL identified by composite interval maps to the earlier interval maps. Two hundred and thirty seven markers were used to construct the core of the E x F map. The majority of marker orders were consistent between the maps. However, 19 putative marker inversions were detected on 12 of 20 linkage groups (LG). Eleven marker distance compressions were also found. The number of inverted markers ranged from 1 to 2 per LG. Thus, marker order inversions may be common in southern compared to northern US germplasm. A total of 61 QTL among 37 measures of six traits were detected by composite interval maps, interval maps and single point analysis. Seventeen of the QTL found in composite intervals had previously been detected among the 29 QTL found in simple interval maps. The genomic locations of the known QTL were more closely delimited. A genome sequencing project to compare Southern and Northern US soybean cultivars would catalog and delimit inverted regions and the associated QTL. Gene introgression in cultivar development programs would be accelerated.

Light is essential for degradation of ribulose-1,5-bisphosphate carboxylase-oxygenase large subunit during sudden death syndrome development in soybean

FUSARIUM SOLANI f. sp. GLYCINES (Fsg) has been reported to produce at least two phytotoxins. Cell-free FSG culture filtrates containing phytotoxins have been shown to develop foliar sudden death syndrome (SDS) in soybean. We have investigated the changes in protein profiles of diseased leaves caused by cell-free FSG culture filtrates prepared from FSG isolates. Two-dimensional sodium dodecyl sulfate-polyacrylamide gel electrophoresis (PAGE) was conducted to investigate the protein profiles of diseased and healthy leaves. An approximately 55 kDa protein was found to be absent in diseased leaves. Matrix-assisted laser desorption-ionization time-of-flight mass spectrometric analyses and a database search revealed that the missing protein is the ribulose 1,5-bisphosphate carboxylase/oxygenase (Rubisco) large subunit, which is involved in carbon assimilation and photorespiration. This result was confirmed by Western blot experiments. We have shown that light is essential for disappearance of the Rubisco large subunit initiated by cell-free FSG culture filtrates. The disappearance of the protein is fairly rapid and occurs within 24 h, presumably due to degradation. Cell-free, FSG culture-induced degradation of the Rubisco large subunit was accompanied by accumulation of reactive oxygen species under light conditions. Terminal deoxynucleotidyl transferase-mediated nick end labelling experiments suggested that programmed cell death was initiated in leaves of seedlings fed with cell-free FSG culture filtrates. These results suggest that, in the presence of light, FSG culture filtrates containing phytotoxins cause degradation of the Rubisco large subunit and accumulation of free radicals and, thereby, initiate programmed cell death leading to foliar SDS development in soybean.

Interaction of Fusarium solani f. sp. glycines and Heterodera glycines in Sudden Death Syndrome of Soybean

ABSTRACT Sudden death syndrome (SDS) of soybean is caused by the soilborne Fusarium solani f. sp. glycines (synonym F. virguliforme). In a sequential approach, two multifactor factorial-design microplot experiments were conducted to investigate the effects of fungal infestation levels and soil moisture on both root necrosis and foliar SDS severity, and the interaction between F. solani f. sp. glycines and Heterodera glycines in fumigated versus nonfumigated soil. In 2003, soybean cv. Spencer was grown in nonfumigated or methyl bromide-fumigated soil and infested with increasing levels of F. solani f. sp. glycines, either under rainfall or irrigated after growth stage V6/R1. In 2004, interactions between F. solani f. sp. glycines and H. glycines were explored in a factorial inoculation design in fumigated or nonfumigated soil, planted to Williams 82 or Cyst-X20-18. In both years, higher levels of foliar SDS severity and root necrosis were found in F. solani f. sp. glycines-infested soils with H. glycines than in soils without the nematode on the soybean cultivars susceptible to both pathogens. Both natural infestations of H. glycines in 2003 and artificially amended populations of H. glycines in 2004 contributed to higher foliar SDS severity. More severe foliar SDS symptoms always were associated with more root necrosis, but elevated levels of root necrosis did not predict severe leaf symptoms. In contrast to the critical role of H. glycines, increasing fungal infestation levels had no significant effects on increasing either foliar SDS symptoms or root necrosis. Effects of moisture regime and fungal infestation levels also were examined in factorial greenhouse and growth chamber experiments. High soil moisture resulted in higher levels of SDS root necrosis. In the greenhouse, root necrosis increased at a higher rate in low soil moisture than the rate in high soil moisture. The two pathogens acted as a complex and the disease development was strongly dependent on high soil moisture.

Quantitative trait loci controlling sulfur containing amino acids, methionine and cysteine, in soybean seeds

Soybean [Glycine max (L.) Merr.] is the single largest source of protein in animal feed. However, a major limitation of soy proteins is their deficiency in sulfur-containing amino acids, methionine (Met) and cysteine (Cys). The objective of this study was to identify quantitative trait loci (QTL) associated with Met and Cys concentration in soybean seed. To achieve this objective, 101 F(6)-derived recombinant inbred lines (RIL) from a population developed from a cross of N87-984-16 x TN93-99 were used. Ground soybean seed samples were analyzed for Met and Cys concentration using a near infrared spectroscopy instrument. Data were analyzed using SAS software and QTL Cartographer. RIL differed (P<0.01) in Met and Cys concentrations, with a range of 5.1-7.3 (g kg(-1) seed dry weight) for Cys and 4.4-8.8 (g kg(-1) seed dry weight) for Met. Heritability estimates on an entry mean basis were 0.14 and 0.57 for Cys and Met, respectively. A total of 94 polymorphic simple sequence repeat molecular genetic markers were screened in the RIL. Single factor ANOVA was used to identify candidate QTL, which were confirmed by composite interval mapping using QTL Cartographer. Four QTL linked to molecular markers Satt235, Satt252, Satt427 and Satt436 distributed on three molecular linkage groups (MLG) D1a, F and G were associated with Cys and three QTL linked to molecular markers Satt252, Satt564 and Satt590 distributed on MLG F, G and M were associated with Met concentration in soybean seed. QTL associated with Met and Cys in soybean seed will provide important information to breeders targeting improvements in the nutritional quality of soybean.

Genomic analysis of a region encompassing QRfs1 and QRfs2: genes that underlie soybean resistance to sudden death syndrome

Candidate genes were identified for two loci, QRfs2 providing resistance to the leaf scorch called soybean (Glycine max (L.) Merr.) sudden death syndrome (SDS) and QRfs1 providing resistance to root infection by the causal pathogen Fusarium solani f.sp. glycines. The 7.5 +/- 0.5 cM region of chromosome 18 (linkage group G) was shown to encompass a cluster of resistance loci using recombination events from 4 near-isogenic line populations and 9 DNA markers. The DNA markers anchored 9 physical map contigs (7 are shown on the soybean Gbrowse, 2 are unpublished), 45 BAC end sequences (41 in Gbrowse), and contiguous DNA sequences of 315, 127, and 110 kbp. Gene density was high at 1 gene per 7 kbp only around the already sequenced regions. Three to 4 gene-rich islands were inferred to be distributed across the entire 7.5 cM or 3.5 Mbp showing that genes are clustered in the soybean genome. Candidate resistance genes were identified and a molecular basis for interactions among the disease resistance genes in the cluster inferred.