SNPs discovery and CAPS marker conversion in soybean

To discover new SNPs and develop an easy assay method in soybean, we compared the high-throughput pyrosequencing ESTs with whole genome sequences in different soybean varieties and identified 3899 SNPs. Transitions were found to be much more frequent than transversions in these SNPs. We found that SNPs were widely distributed in the soybean genome, targeting numerous genes involved in various physiological and biochemical processes influencing important agronomic traits. A set of 16 SNPs were validated in nine soybean varieties, and seven SNPs were converted into CAPS. From functional gene association analysis, the marker CAPS282 on the 3′-UTR of gene Glyma07g03490 was identified as associated with 100-seed weight in soybean. The SNP discovery and CAPS markers conversion system developed in this study is fast and cost effective, and holds great promise for molecular-assisted breeding of soybean.     

Evaluation of cytotoxicity and immune modulatory activities of soyasaponin Ab: An in vitro and in vivo study

To improve the immune efficacy of protein subunit vaccines, novel adjuvants are needed to elicit a suitable protective immune response and to promote long term immunologic memory. In this work, soyasaponin Ab, a major constituent among group A soyasaponins in soybeans was purified and prepared from soy hypocotyls. The immunomodulatory effects of soyasaponin Ab both in vitro and in vivo were investigated, and its pro-immunomodulatory molecular mechanism was also studied. For in vitro assays, with mouse macrophage cell line RAW264.7 as the studying model, both cytotoxicity and immune stimulatory activity were investigated to evaluate the potential of soyasaponin Ab as the vaccine adjuvant. The results indicated that soyasaponin Ab could be significantly safer than Quillaja saponins (QS). Soyasaponin Ab showed no toxicities over the tested concentration ranges compared to QS. Soyasaponin Ab was proved able to promote releases of inflammatory cytokines like TNFα and IL-1β in a dose-dependent manner. Furthermore, NF-κB signalling was also activated by soyasaponin Ab effectively. In addition, with TLR4 gene expression of RAW264.7 cell inhibited by RNA interference, immune stimulatory effects by soyasaponin Ab dropped down significantly. On the other hand, the in vivo experiment results showed that anti-ovalbumin (OVA) IgG, IgG1, IgG2a, IgG2b were significantly enhanced by the soyasaponin Ab and QS groups (p < 0.05 or p < 0.01). The results suggested that compared to QS, soyasaponin Ab may represent a viable candidate for effective vaccine adjuvant. TLR4 receptor dependent pathway may be involved in immune stimulatory effects of soyasaponin Ab.     

Divergence of flowering genes in soybean

Soybean genome sequences were blasted with Arabidopsis thaliana regulatory genes involved in photoperiod-dependent flowering. This approach enabled the identification of 118 genes involved in the flowering pathway. Two genome sequences of cultivated (Williams 82) and wild (IT182932) soybeans were employed to survey functional DNA variations in the flowering-related homologs. Forty genes exhibiting nonsynonymous substitutions between G. max and G. soja were catalogued. In addition, 22 genes were found to co-localize with QTLs for six traits including flowering time, first flower, pod maturity, beginning of pod, reproductive period, and seed filling period. Among the genes overlapping the QTL regions, two LHY/CCA1 genes, GI and SFR6 contained amino acid changes. The recently duplicated sequence regions of the soybean genome were used as additional criteria for the speculation of the putative function of the homologs. Two duplicated regions showed redundancy of both flowering-related genes and QTLs. ID 12398025, which contains the homeologous regions between chr 7 and chr 16, was redundant for the LHY/CCA1 and SPA1 homologs and the QTLs. Retaining of the CRY1 gene and the pod maturity QTLs were observed in the duplicated region of ID 23546507 on chr 4 and chr 6. Functional DNA variation of the LHY/CCA1 gene (Glyma07g05410) was present in a counterpart of the duplicated region on chr 7, while the gene (Glyma16g01980) present in the other portion of the duplicated region on chr 16 did not show a functional sequence change. The gene list catalogued in this study provides primary insight for understanding the regulation of flowering time and maturity in soybean.     

The origin and functional transition of P34

P34, a storage protein and major soybean allergen, has undergone a functional transition from a cysteine peptidase to a syringolide receptor. An exploration of the evolutionary mechanism of this functional transition is made. To identify homologous genes of P34, syntenic network was constructed using syntenic relationships from the Plant Genome Duplication Database. The collected homologous genes, along with SPE31, a highly homologous protein to P34 from the seeds of Pachyrhizus erosus, were used to construct a phylogenetic tree. The results show that multiple gene duplications, exon shuffling and following granulin domain loss and some critical point mutations are associated with the functional transition. Although some tests suggested the existence of positive selection, the possibility that random fixation under relaxation of purifying selection results in the functional transition is also supported. In addition, the genes Glyma08g12340 and Medtr8g086470 may belong to a new group within the papain family.     

Fine mapping of a <Emphasis Type="Italic">Phytophthora</Emphasis>-resistance gene <Emphasis Type="Italic">RpsWY</Emphasis> in soybean (<Emphasis Type="Italic">Glycine max</Emphasis> L.) by high-throughput genome-wide sequencing

Key message

Using a combination of phenotypic screening, genetic and statistical analyses, and high-throughput genome-wide sequencing, we have finely mapped a dominant Phytophthora resistance gene in soybean cultivar Wayao.

Abstract

Phytophthora root rot (PRR) caused by Phytophthora sojae is one of the most important soil-borne diseases in many soybean-production regions in the world. Identification of resistant gene(s) and incorporating them into elite varieties are an effective way for breeding to prevent soybean from being harmed by this disease. Two soybean populations of 191 F₂ individuals and 196 F_7:8 recombinant inbred lines (RILs) were developed to map Rps gene by crossing a susceptible cultivar Huachun 2 with the resistant cultivar Wayao. Genetic analysis of the F₂ population indicated that PRR resistance in Wayao was controlled by a single dominant gene, temporarily named RpsWY, which was mapped on chromosome 3. A high-density genetic linkage bin map was constructed using 3469 recombination bins of the RILs to explore the candidate genes by the high-throughput genome-wide sequencing. The results of genotypic analysis showed that the RpsWY gene was located in bin 401 between 4466230 and 4502773 bp on chromosome 3 through line 71 and 100 of the RILs. Four predicted genes (Glyma03g04350, Glyma03g04360, Glyma03g04370, and Glyma03g04380) were found at the narrowed region of 36.5 kb in bin 401. These results suggest that the high-throughput genome-wide resequencing is an effective method to fine map PRR candidate genes.

     

Genetic dissection of soybean partial resistance to sclerotinia stem rot through genome wide association study and high throughout single nucleotide polymorphisms

Sclerotinia stem rot (SSR) is a disease of soybean [Glycine max (L.) Merr] that causes severe yield losses. We studied 185 representative soybean accessions to evaluate partial SSR resistance and sequenced these by the specific-locus amplified fragment sequencing method. In total, 22,048 single-nucleotide polymorphisms (SNPs), with minor allele frequencies (MAF) ≥5% and missing data <3%, were developed and applied to genome-wide association study of SSR responsiveness and assess linkage disequilibrium (LD) level for candidate gene selection. We identified 18 association signals related to SSR partial resistance. Among them, six overlapped the regions of previous quantitative trait loci, and twelve were novel. We identified 243 candidate genes located in the 200 kb genomic region of these peak SNPs. Based on quantitative real-time polymerase chain reaction and haplotype analysis, Glyma.03G196000 and Glyma.20G095100, encoding pentatricopeptide repeat proteins, might be important factors in the resistance response of soybean to SSR.     

GmPep914, an eight-amino acid peptide isolated from soybean leaves, activates defense-related genes

Only a handful of endogenous peptide defense signals have been isolated from plants. Herein, we report a novel peptide from soybean (Glycine max) leaves that is capable of alkalinizing the media of soybean suspension cells, a response that is generally associated with defense peptides. The peptide, DHPRGGNY, was synthesized and found to be active at 0.25 nM and requiring only 5 to 10 min to obtain a maximal pH change. The peptide is located on the carboxy-terminal end of a 52-amino acid precursor protein (Glyma12g00990) deduced from the soybean genome project. A search of the soybean databank revealed a homolog (Glyma09g36370) that contained a similar peptide, DLPRGGNY, which was synthesized and shown to have identical activity. The peptides, designated GmPep914 (DHPRGGNY) and GmPep890 (DLPRGGNY), were capable of inducing the expression of both Glyma12g00990 (GmPROPEP914) and Glyma09g36370 (GmPROPEP890) in cultured soybean suspension cells within 1 h. Both peptides induced the expression of defense genes, including CYP93A1, a cytochrome P450 gene involved in phytoalexin synthesis, chitinaseb1-1, a chitinase involved in pathogen defense, and Glycine max chalcone synthase1 (Gmachs1), chalcone synthase, involved in phytoalexin production. Both GmPROPEP914 and GmPROPEP890 were highly expressed in the roots, relative to the aerial portions of the plant. However, treatment of the aerial portion of soybean plants with hormones involved in elicitation of defense responses revealed a significant increase in expression levels of GmPROPEP914 and GmPROPEP890. A search of gene databases revealed homologous sequences in other members of the Fabales and also in the closely related Cucurbitales but not in any other order of plants.     

Fine-mapping and identification of a novel locus <Emphasis Type="Italic">Rsc15</Emphasis> underlying soybean resistance to <Emphasis Type="Italic">Soybean mosaic virus</Emphasis>

Key message

Rsc15, a novel locus underlying soybean resistance to SMV, was fine mapped to a 95-kb region on chromosome 6. The Rsc15- mediated resistance is likely attributed to the gene GmPEX14 , the relative expression of which was highly correlated with the accumulation of H ₂ O ₂ along with the activities of POD and CAT during the early stages of SMV infection in RN-9.

Abstract

Soybean mosaic virus (SMV) causes severe yield losses and seed quality deterioration in soybean [Glycine max (L.) Merr.] worldwide. A series of single dominant SMV resistance genes have been identified on respective soybean chromosomes 2, 13 and 14, while one novel locus, Rsc15, underlying resistance to the virulent SMV strain SC15 from soybean cultivar RN-9 has been recently mapped to a 14.6-cM region on chromosome 6. However, candidate gene has not yet been identified within this region. In the present study, we aimed to fine map the Rsc15 region and identify candidate gene(s) for this invaluable locus. High-resolution fine-mapping revealed that the Rsc15 gene was located in a 95-kb genomic region which was flanked by the two simple sequence repeat (SSR) markers SSR_06_17 and BARCSOYSSR_06_0835. Allelic sequence comparison and expression profile analysis of candidate genes inferred that the gene Glyma.06g182600 (designated as GmPEX14) was the best candidate gene attributing for the resistance of Rsc15, and that genes encoding receptor-like kinase (RLK) (i.e., Glyma.06g175100 and Glyma.06g184400) and serine/threonine kinase (STK) (i.e., Glyma.06g182900 and Glyma.06g183500) were also potential candidates. High correlations were established between the relative expression level of GmPEX14 and the hydrogen peroxide (H₂O₂) concentration and activities of catalase (CAT) and peroxidase (POD) during the early stages of SMV-SC15 infection in RN-9. The results of the present study will be useful in marker-assisted breeding for SMV resistance and will lead to further understanding of the molecular mechanisms of host resistance against SMV.

     

Direct Determination of Copper,Lead and Cadmium in Tea Infusions by Flameless Atomic Absorption Spectrometry

Two kinds of “group A saponin,” Aa and Ab, are present as the main constituent in soybean seed. The saponin composition and content in F₁ and F₂ seeds derived from the cross parents of Aa and Ab types were analyzed. The “group A saponin” was of Aa–Ab type in all the F₁ seeds, and the ratio of Aa type : Aa–Ab type : Ab type was 1:2:1 in the F₂ seeds. From these results, it appears that Aa and Ab were controlled by codominant allelic alternatives at a single locus. An investigation of the saponin composition of the seed hypocotyls of 18 wild lines revealed some lines in which “group A saponin” was absent.     

Rapid genotyping of soybean cultivars using high throughput sequencing

Soybean (Glycine max) breeding involves improving commercially grown varieties by introgressing important agronomic traits from poor yielding accessions and/or wild relatives of soybean while minimizing the associated yield drag. Molecular markers associated with these traits are instrumental in increasing the efficiency of producing such crosses and Single Nucleotide Polymorphisms (SNPs) are particularly well suited for this task, owing to high density in the non-genic regions and thus increased likelihood of finding a tightly linked marker to a given trait. A rapid method to develop SNP markers that can differentiate specific loci between any two parents in soybean is thus highly desirable. In this study we investigate such a protocol for developing SNP markers between multiple soybean accessions and the reference Williams 82 genome. To restrict sampling frequency reduced representation libraries (RRLs) of genomic DNA were generated by restriction digestion followed by library construction. We chose to sequence four accessions Dowling (PI 548663), Dwight (PI 597386), Komata (PI200492) and PI 594538A for their agronomic importance as well as Williams 82 as a control.MseI was chosen to digest genomic DNA based on predictions that it will cut sparingly in the mathematically defined high-copy-number regions of the genome. All RRLs were sequenced on the Illumina genome analyzer. Reads were aligned to the Glyma1 reference assembly and SNP calls made from the alignments. We identified from 4294 to 14550 SNPs between the four accessions and the Williams 82 reference. In addition a small number of SNPs (1142) were found by aligning Williams 82 reads to the reference assembly (Glyma1) suggesting limited genetic variation within the Williams 82 line. The SNP data allowed us to estimate genetic diversity between the four lines and Williams 82. Restriction digestion of soybean genomic DNA with MseI followed by high throughput sequencing provides a rapid and reproducible method for generating SNP markers.     

不同Bt蛋白对棉铃虫存活率和生长发育的影响

在我国Bt棉主要以Cry1Ab或Cry1Ac为主,其他新型Bt基因未被转入棉花中用来控制害虫,然而大面积种植单价Bt基因的棉花,将可能会大大增加靶标害虫对该类型Bt棉花抗性频率,因此研究其他新型Bt蛋白对靶标害虫的控制作用显得十分必要。采用蛋白混入人工饲料的生物测定方法,在室内测定了6种Bt蛋白对棉铃虫初孵幼虫的毒力,比较了浓度为1．0μg· g-1时不同Bt蛋白对棉铃虫幼虫生长发育的影响。毒力测定结果表明,不同Bt蛋白对棉铃虫初孵幼虫的毒力不同,LC50值由低到高依次为Cry1Ab 0．065μg· g-1、Cry1Ac 0．074μg· g-1、Cry2Ab 0．133μg· g-1、Cry2Aa 11．670μg· g-1、Cry1Ah 13．010μg· g-1和Cry1Ca＞20μg· g-1。生长发育测定结果表明,Cry1Ab和Cry1Ac对棉铃虫幼虫的生长发育影响最大,Cry2Ab次之;Cry1Ah和Cry2Aa对1龄幼虫的校正死亡率和体重抑制率差别不大,但对2龄幼虫的差异较大,Cry1Ah处理2龄幼虫后体重和生长发育参数与Cry2Ab接近,而Cry1Ca对棉铃虫幼虫生长发育几乎没影响。Cry1Ah、Cry2Aa和Cry2Ab的毒力不如Cry1Ac和Cry1Ab,但仍可以作为控制棉铃虫幼虫的替代策略。     

Single nucleotide polymorphism markers for rapid detection of the <Emphasis Type="Italic">Rsv4</Emphasis> locus for soybean mosaic virus resistance in diverse germplasm

Soybean mosaic virus (SMV) causes a substantial decrease in soybean yield and reduction of seed quality. The most effective management strategy to control the virus is the deployment of host resistance. Seven SMV strains and three independent multi-allelic loci for SMV resistance have been identified previously. The goal of this research was to detect single nucleotide polymorphisms (SNPs) associated with SMV resistance at the Rsv4 locus. Ten soybean accessions, with confirmed resistance genes, were used for sequencing the candidate gene Glyma.02g121400. Alignment of these sequences revealed three SNPs displaying 100% consistency for genotypes carrying the Rsv4 gene. These SNPs were applied for a rapid screen of diverse soybean germplasm using the Sequenom iPLEX Gold platform, phenotyped with SMV-G1 and G7 strains to determine phenotype and classified into several groups carrying the proposed R-gene. The population of V94-5152 (Rsv4) × Lee 68 (rsv) was screened using novel SNPs to create a genetic map with improved resolution to determine the location of the Rsv4. To observe the recombination frequencies within the population, three additional SNPs on both sides of the Glyma.02g121400 gene were added. A linkage map revealed a distance of 3.6 cM between the Rsv4 locus and the closest SNP, thus shifting the putative Rsv4 region downstream on chromosome 2. With this region, five candidate genes have been proposed. The genomic position of the discovered SNPs, linked to the Rsv4, could increase screening precision and accelerate breeding efforts to develop multi-strain-resistant crops.     

A candidate male-fertility female-fertility gene tagged by the soybean endogenous transposon, Tgm9

In soybean, the W4 gene encoding dihydroflavonol-4-reductase controls anthocyanin pigment biosynthesis in flowers. The mutant allele, w4-m, is characterized by variegated flowers and was evolved from the insertion of an endogenous transposable element, Tgm9, in intron II of the W4 gene. In the w4-m mutant line, reversion of the unstable allele from variegated to normal purple flower in revertants would indicate Tgm9’s excision accompanied by its insertion into a second locus. We identified a male-sterile, female-sterile mutant from such germinal revertant bearing purple flowers. The objectives of our investigation were to map the sterility locus, identify candidate genes for the male-fertile, female-fertile phenotype, and then determine if sterility is associated with the insertion of Tgm9 in the sterility locus. We used bulked segregant analysis to map the locus to molecular linkage group J (chromosome 16). Fine mapping enabled us to flank the locus to a 62-kb region that contains only five predicted genes. One of the genes in that region, Glyma16g07850.1, codes for a helicase. A rice homolog of this gene has been shown to control crossing over and fertility phenotype. Thus, Glyma16g07850.1 is most likely the gene regulating the male and female fertility phenotype in soybean. DNA blot analysis of the segregating individuals for Tgm9 showed perfect association between sterility and the presence of the transposon. Most likely, the sterility mutation was caused by the insertion of Tgm9. The transposable element should facilitate identification of the male- and female-fertility gene. Characterization of the fertility gene will provide vital molecular insight on the reproductive biology of soybean and other plants.     

The Sg-1 glycosyltransferase locus regulates structural diversity of triterpenoid saponins of soybean

Triterpene saponins are a diverse group of biologically functional products in plants. Saponins usually are glycosylated, which gives rise to a wide diversity of structures and functions. In the group A saponins of soybean (Glycine max), differences in the terminal sugar species located on the C-22 sugar chain of an aglycone core, soyasapogenol A, were observed to be under genetic control. Further genetic analyses and mapping revealed that the structural diversity of glycosylation was determined by multiple alleles of a single locus, Sg-1, and led to identification of a UDP-sugar-dependent glycosyltransferase gene (Glyma07g38460). Although their sequences are highly similar and both glycosylate the nonacetylated saponin A0-αg, the Sg-1(a) allele encodes the xylosyltransferase UGT73F4, whereas Sg-1(b) encodes the glucosyltransferase UGT73F2. Homology models and site-directed mutagenesis analyses showed that Ser-138 in Sg-1(a) and Gly-138 in Sg-1(b) proteins are crucial residues for their respective sugar donor specificities. Transgenic complementation tests followed by recombinant enzyme assays in vitro demonstrated that sg-1(0) is a loss-of-function allele of Sg-1. Considering that the terminal sugar species in the group A saponins are responsible for the strong bitterness and astringent aftertastes of soybean seeds, our findings herein provide useful tools to improve commercial properties of soybean products.     

Fine mapping and analyses of <Emphasis Type="Italic">R</Emphasis><Subscript><Emphasis Type="Italic">SC8</Emphasis></Subscript> resistance candidate genes to soybean mosaic virus in soybean

Soybean mosaic virus (SMV) in soybean [Glycine max (L.) Merr.] is a destructive foliar disease in soybean-producing countries worldwide. In this study, F₂, F_2:3, and F_7:11 recombinant inbred lines populations derived from Kefeng No.1 × Nannong 1138-2 were used to study inheritance and linkage mapping of the SMV strain SC8 resistance gene in Kefeng No.1. Results indicated that a single dominant gene (designated R _SC8 ) controls resistance, which is located on chromosome 2 (MLG D1b). A mixed segregating population was developed by selfing two heterozygous plants (RHL153-1 and RHL153-2) at four markers adjacent to the locus and used in fine mapping R _SC8 . In addition, two genomic-simple sequence repeats (SSR) markers BARCSOYSSR_02_0610 and BARCSOYSSR_02_0616 were identified that flank the two sides of R _SC8 . Sequence analysis of the soybean genome indicated that the interval between the two genomic-SSR markers is 200 kb. QRT-PCR analysis of the candidate genes determined that five genes (Glyma02g13310, 13320, 13400, 13460, and 13470) are likely involved in soybean SMV resistance. These results will have utility in cloning, transferring, and pyramiding of the R _SC8 through marker-assisted selection in soybean breeding programs.     

Complete genome sequence of Bacillus thuringiensis BR145, a strain with insecticidal activity against Lepidoptera pests

Bacillus thuringiensis BR145 isolated from a soybean field in Southern Brazil showed toxicity against two important insect pests from soybean crop, Helicoverpa armigera, and Chrysodeixis includens, with LC₅₀ 0.294 µg.cm^-2 and 0.277 µg.cm^-2, respectively. We analyzed the genome of this strain through sequences obtained by Next Generation DNA Sequencing and de novo assembly. The analysis of the genome revealed insecticidal genes cry1Aa, cry1Ab, cry1Ac, cry1Ia, cry2Ab, cyt1, and vip3Aa, suggesting the use of this strain in new strategies of biological control.     

RNA sequencing analysis of salt tolerance in soybean (Glycine max)

Salt stress causes foliar chlorosis and scorch, plant stunting, and eventually yield reduction in soybean. There are differential responses, namely tolerance (excluder) and intolerance (includer), among soybean germplasm. However, the genetic and physiological mechanisms for salt tolerance is complex and not clear yet. Based on the results from the screening of the RA-452 x Osage mapping population, two F_4:6 lines with extreme responses, most tolerant and most sensitive, were selected for a time-course gene expression study in which the 250 mM NaCl treatment was initially imposed at the V1 stage and continued for 24 h (hrs). Total RNA was isolated from the leaves harvested at 0, 6, 12, 24 h after the initiation of salt treatment, respectively. The RNA-Seq analysis was conducted to compare the salt tolerant genotype with salt sensitive genotype at each time point using RNA-Seq pipeline method. A total of 2374, 998, 1746, and 630 differentially expressed genes (DEGs) between salt-tolerant line and salt-sensitive line, were found at 0, 6, 12, and 24 h, respectively. The expression patterns of 154 common DEGs among all the time points were investigated, of which, six common DEGs were upregulated and seven common DEGs were downregulated in salt-tolerant line. Moreover, 13 common DEGs were dramatically expressed at all the time points. Based on Log2 (fold change) of expression level of salt-tolerant line to salt-sensitive line and gene annotation, Glyma.02G228100, Glyma.03G226000, Glyma.03G031000, Glyma.03G031400, Glyma.04G180300, Glyma.04G180400, Glyma.05 g204600, Glyma.08G189600, Glyma.13G042200, and Glyma.17G173200, were considered to be the key potential genes involving in the salt-tolerance mechanism in the soybean salt-tolerant line.     

Quantitative trait loci controlling aluminum tolerance in soybean: candidate gene and single nucleotide polymorphism marker discovery

Aluminum (Al) toxicity is an important abiotic stress that affects soybean production in acidic soils throughout the world. Development of Al-tolerant cultivars is an efficient and environmentally friendly solution to the problem. A previous report identified quantitative trait loci (QTL) for Al tolerance inherited from PI 416937, using restriction fragment length polymorphism markers, in a population of Young × PI 416937. The population was genotyped with 162 simple sequence repeats to enhance the power of QTL detection and enable the selection of candidate genes for functional marker development. Two QTL that explained 54 % of the phenotypic variation in root extension under Al stress conditions (HIAL) were refined on chromosomes (chr) Gm08 and Gm16. Three QTL located on chr Gm08, Gm16 and Gm19 explained 59 % of the phenotypic variation in root extension as a percent of control (PC). Two major QTL, designated qAL_HIAL_08 and qAL_PC_08, controlling HIAL and PC, respectively, were mapped to the same genomic region on chr Gm08 and inherited their favorable allele from PI 416937. These QTL explained 45 and 41 % of phenotypic variation in HIAL and PC, respectively. Six homologues for citrate synthase (CS) genes were found in the soybean genome sequence at chr Gm02, Gm08, Gm14, Gm15, and Gm18. Sixteen single nucleotide polymorphisms (SNPs) were identified in the CS homologue on chr Gm08. A SimpleProbe assay of Glyma08g42400-SNP was developed for the major QTL on chr Gm08. The SNPs identified from this region could be used for marker-assisted selection of Al tolerance.