Construction and comparison of three reference‐quality genome assemblies for soybean

We report reference‐quality genome assemblies and annotations for two accessions of soybean (Glycine max) and for one accession of Glycine soja, the closest wild relative of G. max. The G. max assemblies provided are for widely used US cultivars: the northern line Williams 82 (Wm82) and the southern line Lee. The Wm82 assembly improves the prior published assembly, and the Lee and G. soja assemblies are new for these accessions. Comparisons among the three accessions show generally high structural conservation, but nucleotide difference of 1.7 single‐nucleotide polymorphisms (snps) per kb between Wm82 and Lee, and 4.7 snps per kb between these lines and G. soja. snp distributions and comparisons with genotypes of the Lee and Wm82 parents highlight patterns of introgression and haplotype structure. Comparisons against the US germplasm collection show placement of the sequenced accessions relative to global soybean diversity. Analysis of a pan‐gene collection shows generally high conservation, with variation occurring primarily in genomically clustered gene families. We found approximately 40–42 inversions per chromosome between either Lee or Wm82v4 and G. soja, and approximately 32 inversions per chromosome between Wm82 and Lee. We also investigated five domestication loci. For each locus, we found two different alleles with functional differences between G. soja and the two domesticated accessions. The genome assemblies for multiple cultivated accessions and for the closest wild ancestor of soybean provides a valuable set of resources for identifying causal variants that underlie traits for the domestication and improvement of soybean, serving as a basis for future research and crop improvement efforts for this important crop species.     

Composite and clinal distribution of Glycine soja in Japan revealed by RFLP analysis of mitochondrial DNA

 Wild soybean (Glycine soja Sieb. et Zucc.), regarded as the progenitor of cultivated soybean [G. max (L.) Merr.], is widely distributed in East Asia. We have collected 1097 G. soja plants from all over Japan and analyzed restriction fragment length polymorphisms (RFLPs) of mitochondrial DNA (mtDNA) in them. Based on the RFLPs detected by gel-blot analysis, using coxII and atp6 as probes, the collected plants were divided into 18 groups. Five mtDNA types accounted for 94% of the plants examined. The geographic distribution of mtDNA types revealed that, in many regions, wild soybeans grown in Japan consisted of a mixture of plants with different types of mtDNA, occasionally even within sites. Some of the mtDNA types showed marked geographic clines among the regions. Additionally, some wild soybeans possessed mtDNA types that were identical to those widely detected in cultivated soybeans. Our results suggest that the analysis of mtDNA could resolve the maternal lineage among plants of the genus Glycine subgenus Soja. Received: 16 June 1997/Accepted: 5 August 1997     

Development,validation and genetic analysis of a large soybean SNP genotyping array

Cultivated soybean (Glycine max) suffers from a narrow germplasm relative to other crop species, probably because of under‐use of wild soybean (Glycine soja) as a breeding resource. Use of a single nucleotide polymorphism (SNP) genotyping array is a promising method for dissecting cultivated and wild germplasms to identify important adaptive genes through high‐density genetic mapping and genome‐wide association studies. Here we describe a large soybean SNP array for use in diversity analyses, linkage mapping and genome‐wide association analyses. More than four million high‐quality SNPs identified from high‐depth genome re‐sequencing of 16 soybean accessions and low‐depth genome re‐sequencing of 31 soybean accessions were used to select 180 961 SNPs for creation of the Axiom^® SoyaSNP array. Validation analysis for a set of 222 diverse soybean lines showed that 170 223 markers were of good quality for genotyping. Phylogenetic and allele frequency analyses of the validation set data indicated that accessions showing an intermediate morphology between cultivated and wild soybeans collected in Korea were natural hybrids. More than 90 unanchored scaffolds in the current soybean reference sequence were assigned to chromosomes using this array. Finally, dense average spacing and preferential distribution of the SNPs in gene‐rich chromosomal regions suggest that this array may be suitable for genome‐wide association studies of soybean germplasm. Taken together, these results suggest that use of this array may be a powerful method for soybean genetic analyses relating to many aspects of soybean breeding.     

Diversity of chloroplast DNA SSRs in wild and cultivated soybeans: evidence for multiple origins of cultivated soybean

Soybean [ Glycine max (L.) Merr.] is one of the major crops in the world and was domesticated from a wild progenitor, Glycine soja Sieb. & Zucc., in East Asia. In order to address the questions concerning the evolution and maternal lineage of soybean, we surveyed the variation in chloroplast DNA simple sequence repeats (cpSSR) of 326 wild and cultivated soybean accessions that were collected from various Asian countries. Twenty-three variants were detected at six cpSSRs in the accessions tested. All of the variants were found in wild soybean, whereas only 14 variants existed in the cultigen. Combining the variants at the six cpSSRs gave 52 haplotypes in the former and eight haplotypes in the latter. Both analyses indicated a considerably higher genetic diversity in the wild soybean. Around 75% of the cultivated accessions tested possessed a common haplotype (no. 49), which was detected in only seven wild accessions, six from southern Japan and one from southern China. The predominant haplotype in the cultigen may therefore have originated from a rare haplotype of the wild soybean that is presently distributed in the southern areas of Japan and China. The remaining seven haplotypes in the cultigen were distributed regionally, and except for three rare haplotypes, largely overlapped with the distributions of wild accessions with the same respective haplotypes. Our results strongly suggest that the cultivated soybeans with different cpDNA haplotypes originated independently in different regions from different wild gene pools and/or hybrid swarms between cultivated and wild forms.     

栽培大豆和野生大豆线粒体基因组密码子使用偏性的比较分析

为分析栽培大豆和野生大豆线粒体基因组的密码子使用特征差异,该文以其线粒体基因组编码序列为研究对象,比较其密码子偏性形成的影响因素和演化过程。结果表明:(1)栽培大豆和野生大豆线粒体基因组编码区的GC含量分别为44.56%和44.58%,说明栽培大豆和野生大豆线粒体编码基因均富含A/T碱基。(2)栽培大豆和野生大豆线粒体基因组密码子第1位、第2位GC含量平均值与第3位GC含量的相关性均呈极显著水平,说明突变在其密码子偏性形成中的作用不可忽略; PR2-plot分析显示,在同义密码子第3位碱基的使用频率上,嘌呤低于嘧啶; Nc-plot分析中Nc比值位于-0.1～0.2区间的基因数占总基因数的95%以上;突变和选择等多重因素共同作用影响了大豆线粒体基因组编码序列密码子使用偏性的形成。(3)有20、21个密码子分别被确定为栽培大豆和野生大豆线粒体基因组编码序列的最优密码子,其中除丝氨酸TCC密码子外均以A或T结尾。综上结果认为,栽培大豆线粒体密码子偏性的形成受选择的影响要高于野生大豆,这可能是栽培大豆由野生大豆经长期人工栽培驯化的结果。     

Differential responses of molecular mechanisms and physiochemical characters in wild and cultivated soybeans against invasion by the pathogenic Fusarium oxysporum Schltdl

Cultivated soybean (Glycine max) was derived from the wild soybean (Glycine soja), which has genetic resources that can be critically important for improving plant stress resistance. However, little information is available pertaining to the molecular and physiochemical comparison between the cultivated and wild soybeans in response to the pathogenic Fusarium oxysporum Schltdl. In this study, we first used comparative phenotypic and paraffin section analyses to indicate that wild soybean is indeed more resistant to F. oxysporum than cultivated soybean. Genome‐wide RNA‐sequencing approach was then used to elucidate the genetic mechanisms underlying the differential physiological and biochemical responses of the cultivated soybean, and its relative, to F. oxysporum. A greater number of genes related to cell wall synthesis and hormone metabolism were significantly altered in wild soybean than in cultivated soybean under F. oxysporum infection. Accordingly, a higher accumulation of lignins was observed in wild soybean than cultivated soybean under F. oxysporum infection. Collectively, these results indicated that secondary metabolites and plant hormones may play a vital role in differentiating the response between cultivated and wild soybeans against the pathogen. These important findings may provide future direction to breeding programs to improve resistance to F. oxysporum in the elite soybean cultivars by taking advantage of the genetic resources within wild soybean germplasm.     

Fine-scale phylogenetic structure and major events in the history of the current wild soybean (Glycine soja) and taxonomic assignment of semi-wild type (Glycine gracilis Skvortz.) within the Chinese subgenus Soja

Wild and cultivated species of soybeans have coexisted for 5000 years in China. Despite this long history, there is very little information on the genetic relationship of Glycine soja and G. max. To gain insight into the major events in the history of the subgenus Soja, we examined 20 simple sequence repeat (SSR) markers of a large number of accessions (910). The results showed no significant differences between wild and semi-wild soybeans in genetic diversity but significant differences between G. soja and G. max. Ancestry and cluster analyses revealed that semi-wild soybeans should belong to the wild category and not to G. max. Our results also showed that differentiation had occurred not only among G. soja, G. gracilis, and G. max but also within G. soja and within G. gracilis. Glycine soja had 3 clear genetic categories: typical small-seeded (≤2.0 g 100-seed weight), dual-origin middle-seeded (2.0-2.5 g), and large-seeded plants (2.51-3.0 g). These last were genetically close to G. gracilis, their defining some traits having been acquired mainly by introgression from soybeans. Small-seeded G. gracilis (3.01-3.5 g) were genetically different from larger seeded ones (from 3.51 to 4.0 to over 10 g). Seed size predominated over seed coat color in evolutionary degree. Typical and large-seeded G. soja were found to have 0.7% and 12% introgressive cultivar genes, respectively. The genetic boundary of G. gracilis was at the range of 2.51-3.0 g of G. soja. In the great majority of wild accessions, traits such as white flowers, gray pubescences, no-seed bloom, and colored seed coats were likely introgressive from domesticated soybeans.     

Single nucleotide mutation leading to an amino acid substitution in the variant <Emphasis Type="Italic">Tik</Emphasis> soybean Kunitz trypsin inhibitor (SKTI) identified in Chinese wild soybean (<Emphasis Type="Italic">Glycine soja</Emphasis> Sieb. &; Zucc.)

The wild soybean (Glycine soja), which is the progenitor of cultivated soybean (Glycine max), is expected to offer more information about genetic variability and more useful mutants for evolutionary research and breeding applications. Here, a total of 1,600 wild soybean samples from China were investigated for genetic variation with regard to the soybean Kunitz trypsin inhibitor (SKTI). A new mutant SKTI, Tik, was identified. It was found to be a Tia-derived codominant allele caused by a transversion point mutation from C to G at nucleotide +171, leading to an alteration of one codon (AAC → AAG) and a corresponding amino acid substitution (Asn → Lys) at the ninth residue. Upon examination of this variant and others previously found in wild soybeans, it became clear that SKTI has undergone high-level evolutionary differentiation. There were more abundant polymorphisms in the wild than in the cultivated soybean.     

Population genetic structure of Japanese wild soybean (Glycine soja) based on microsatellite variation

The research objectives were to determine aspects of the population dynamics relevant to effective monitoring of gene flow in the soybean crop complex in Japan. Using 20 microsatellite primers, 616 individuals from 77 wild soybean (Glycine soja) populations were analysed. All samples were of small seed size (< 0.03 g), were directly collected in the field and came from all parts of Japan where wild soybeans grow, except Hokkaido. Japanese wild soybean showed significant reduction in observed heterozygosity, low outcrossing rate (mean 3.4%) and strong genetic differentiation among populations. However, the individual assignment test revealed evidence of rare long-distance seed dispersal (> 10 km) events among populations, and spatial autocorrelation analysis revealed that populations within a radius of 100 km showed a close genetic relationship to one another. When analysis of graphical ordination was applied to compare the microsatellite variation of wild soybean with that of 53 widely grown Japanese varieties of cultivated soybean (Glycine max), the primary factor of genetic differentiation was based on differences between wild and cultivated soybeans and the secondary factor was geographical differentiation of wild soybean populations. Admixture analysis revealed that 6.8% of individuals appear to show introgression from cultivated soybeans. These results indicated that population genetic structure of Japanese wild soybean is (i) strongly affected by the founder effect due to seed dispersal and inbreeding strategy, (ii) generally well differentiated from cultivated soybean, but (iii) introgression from cultivated soybean occurs. The implications of the results for the release of transgenic soybeans where wild soybeans grow are discussed.     

Population Structure and Domestication Revealed by High-Depth Resequencing of Korean Cultivated and Wild Soybean Genomes

Despite the importance of soybean as a major crop, genome-wide variation and evolution of cultivated soybeans are largely unknown. Here, we catalogued genome variation in an annual soybean population by high-depth resequencing of 10 cultivated and 6 wild accessions and obtained 3.87 million high-quality single-nucleotide polymorphisms (SNPs) after excluding the sites with missing data in any accession. Nuclear genome phylogeny supported a single origin for the cultivated soybeans. We identified 10-fold longer linkage disequilibrium (LD) in the wild soybean relative to wild maize and rice. Despite the small population size, the long LD and large SNP data allowed us to identify 206 candidate domestication regions with significantly lower diversity in the cultivated, but not in the wild, soybeans. Some of the genes in these candidate regions were associated with soybean homologues of canonical domestication genes. However, several examples, which are likely specific to soybean or eudicot crop plants, were also observed. Consequently, the variation data identified in this study should be valuable for breeding and for identifying agronomically important genes in soybeans. However, the long LD of wild soybeans may hinder pinpointing causal gene(s) in the candidate regions.     

The genetic diversity of annual wild soybeans grown in China

Annual wild soybeans (Glycine soja), the ancestors of cultivated soybeans (G. max), are important sources of major genes for resistance to pests, diseases and environmental stresses. The study of their genetic diversity is invaluable for efficient utilization, conservation and management of germplasm collections. In this paper, the number of accessions, the variation of traits, the genetic diversity indexes (Shannon index) and the coefficient of variation were employed to study the geographical distribution of accessions, genetic diversity of characters and genetic diversity centers of annual wild soybean by statistical analysis of the database from the National Germplasm Evaluation Program of China. Most annual wild soybeans are distributed in Northeast China, and the number of accessions decreases from the Northeast to other directions in China. The genetic diversity indexes (Shannon index) were 0.49, 0.74, 0.02, 0.55, 1.45, 2.41, 1.27 and 1.89 for flower color, sootiness of seed coat, cotyledon color, pubescence color, hilum color, leaf shape, stem type and seed color, respectively. Coefficients of variation were 7.1%, 28.7%, 76.43% and 18.2% for protein content, oil content, 100-seed weight and days to maturity, respectively. Three genetic diversity centers, the Northeast, the Yellow River Valley and the Southeast Coasts of China, are proposed based on the geographical distribution of the number of accessions, genetic diversity and the multivariate variation coefficient. Based on these results and Vavilov’s theory of crop origination, two opposing possible models for the formation of the three centers are proposed, either these centers are independent of each other and the annual wild soybeans in these centers originated separately, or the Northeast center was the primary center for annual wild soybeans in China, while the Yellow River Valley center was derived from this primary center and served as the origin for the Southeast Coast center. Received: 25 June 2000 / Accepted: 18 October 2000     

Sequence variation of non-coding regions of chloroplast DNA of soybean and related wild species and its implications for the evolution of different chloroplast haplotypes

Soybean chloroplast DNAs (cpDNAs) are classified into three types (I, II and III) based on RFLP profiles. Type I is mainly observed in cultivated soybean (Glycine max), while type II and type III are frequently found in both cultivated and wild soybean (Glycine soja), although type III is predominant in wild soybean. In order to evaluate the diversity of cpDNA and to determine the phylogenetic relationship among different chloroplast types, we sequenced nine non-coding regions of cpDNA for seven cultivated and 12 wild soybean accessions with different cpDNA types. Eleven single-base substitutions and a deletion of five bases were detected in a total of 3849 bases identified. Five mutations distinguished the accessions with types I and II from those with type III, and seven were found in the accessions with type III, independently of their taxa. Four species of the subgenus Glycine shared bases that were identical to those with types I and II at two of the five mutation sites and shared bases that were identical to those with type III at the remaining three sites. Therefore, the different cpDNA types may not have originated monophyletically, but rather may have differentiated from a common ancestor in different evolutionary directions. A neighbor-joining tree resulting from the sequence data revealed that the subgenus Soja connected with Glycine microphylla which formed a distinct clade from Clycine clandestina and the tetraploid cytotypes of Glycine tabacina and Glycine tomentella. Several informative length mutations of 54 to 202 bases, due to insertions or deletions, were also detected among the species of the genus Glycine. Received: 16 December 1999 / Accepted: 12 February 2000     

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.     

Molecular sequence variations of the lipoxygenase-2 gene in soybean

Soybean lipoxygenase genes comprise a multi-gene family, with the seed lipoxygenase isozymes LOX1, LOX2, and LOX3 present in soybean seeds. Among these, the LOX2 isozyme is primarily responsible for the “beany” flavor of most soybean seeds. The variety, Jinpumkong 2, having null alleles (lx1, lx2, and lx3) lacks the three seed lipoxygenases; so, sequence variations between the lipoxygenase-2 genes of Pureunkong (Lx2) and Jinpumkong 2 (lx2) cultivars were examined. One indel, four single nucleotide polymorphisms (SNPs), a 175-bp fragment in the 5′-flanking sequence, and a missense mutation within the coding region were found in Jinpumkong 2. The distribution of the sequence variations was investigated among 90 recombinant inbred lines (RILs) derived from a cross of Pureunkong × Jinpumkong 2 and in 480 germplasm accessions with various origins and maturity groups. Evidence for a genetic bottleneck was observed: the 175-bp fragment was rare in Glycine max, but present in the majority of the G. soja accessions. Furthermore, the 175-bp fragment was not detected in the 5′ upstream region of the Lx2 gene on chromosome (Chr) 13 in Williams 82; instead, a similar 175-bp fragment was positioned in the homeologous region on Chr 15. The findings indicated that the novel fragment identified was originally present in the Lx2 region prior to the recent genome duplication in soybean, but became rare in the G. max gene pool. The missense mutation of the conserved histidine residue of the lx2 allele was developed into a single nucleotide-amplified polymorphism (SNAP) marker. The missense mutation showed a perfect correlation with the LOX2-lacking phenotype, so the SNAP marker is expected to facilitate breeding of soybean cultivars which lack the LOX2 isozyme.     

Genetic divergence between North American ancestral soybean lines and introductions with resistance to soybean cyst nematode revealed by chloroplast haplotype

Domesticated soybean [Glycine max (L.) Merr.] is a major crop with an established ancestral relationship to wild soybean (Glycine soja Sieb. & Zucc.) native to Asia. Soybean genetic diversity can be assessed at different levels by identification of polymorphic alleles at genetic loci, in either the plastid or nuclear genomes. The objective of this study was to evaluate genetic diversity based on chloroplast haplotypes for soybean genotypes present in the USDA germplasm resource collection. Shared chloroplast haplotypes represent broad groups of genetic relatedness. Previous work categorized three-quarters of the cultivated soybeans from Asia into a single haplotype group. Our results confirmed the close relationship of North American soybean ancestors and G. max plant introductions previously identified as representing potential sources of soybean genetic variation with the finding that these genotypes belonged to a single chloroplast haplotype group. Genetic diversity was identified in soybean genotypes determined to have a high density of single nucleotide polymorphisms and in a screen of accessions with resistance to soybean cyst nematode. Characterization of soybean plant introduction lines into chloroplast haplotype group may be an important initial step in evaluating the appropriate use of particular soybean genotypes.     

A platform for soybean molecular breeding: the utilization of core collections for food security

Soybean is an important crop not only for human consumption but also for its addition of nitrogen to the soil during crop rotation. China has the largest collection of cultivated soybeans (Glycine max) and wild soybeans (Glycine soja) all over the world. The platform of soybean core, mini core and integrated applied core collections has been developed in the past decade based on systematic researches which included the sampling strategies, statistical methods, phenotypic data and SSR markers. Meanwhile, intergrated applied core collections including accessions with single or integrated favorite traits are being developed in order to meet the demand of soybean breeding. These kinds of core collections provide powerful materials for evaluation of germplasm, identification of trait-specific accessions, gene discovery, allele mining, genomic study, maker development, and molecular breeding. Some successful cases have proved the usefulness and efficiency of this platform. The platform is helpful for enhancing utilization of soybean genetic resources in sustainable crop improvement for food security. The efficient utilization of this platform in the future is relying on accurate phenotyping methods, abundant functional markers, high-throughput genotyping platforms, and effective breeding programs.     

Nucleotide diversity of the upstream region of the putative MADS-box gene controlling soybean maturity

MADS-box genes are involved in plant reproductive development. However, the role of gene nucleotide diversity in soybean flowering and maturity remains unknown. Therefore, in this study, the distribution of DNA polymorphisms in the putative MADS-box gene located near the quantitative trait loci (QTL) for flowering time and maturity was targeted for association analysis using Glycine max (cultivated soybean) and Glycine soja (wild soybean). Sixteen single nucleotide polymorphisms identified in the upstream region of the putative MADS-box gene around QTL Pod mat 13-7 and Fflr 4-2 on chromosome 7 were found to be highly associated with maturity in soybean. The genetic diversity between cultivated soybeans and the wild relative was comparable, although the early maturity group (EMG) was less diverse than the late maturity group (LMG) of the cultivated soybean. Population size changes of the MADS-box gene in this soybean germplasm appeared to result from non-random selection. A selective pressure seemed to act on this gene in the EMG, while the LMG and G. soja were in genetic equilibrium. Neutrality tests and the constructed neighbor-joining tree indicate that the EMG of G. max has experienced strong artificial selection for its domestication and genetic improvement.     

盐生野大豆的异黄酮积累及其生态学意义

以自然生长在盐碱地上的野大豆(Glycine soja)和不耐盐的栽培大豆(G. max)为材料,测定了它们在不同盐度条件下叶片、根部和种子的异黄酮含量,并测定了它们叶片的L-苯丙氨酸含量和苯丙氨酸裂解酶(PAL)活性,还测定了它们根部的结瘤量和固氮酶活性。通过两者比较,分析了它们的大豆异黄酮代谢和盐渍环境的关系。结果表明:盐渍处理不抑制盐生野大豆PAL酶的活性,其大豆异黄酮大量积累;相反,盐渍处理明显抑制栽培大豆PAL酶活性,其大豆异黄酮含量减少,而大豆异黄酮合成前体L-苯丙氨酸积累。结果还显示:在盐渍条件下,盐生野大豆根部异黄酮积累的同时,其根瘤结瘤量较多,且固氮酶活性也较高;而栽培大豆随着其根部异黄酮的减少,其根瘤结瘤量大大减少,且固氮活性大大下降。野大豆和栽培大豆的这些差别说明:盐生野大豆积累大豆异黄酮有其生态学意义,这很可能是野大豆通过异黄酮次生代谢途径适应盐渍环境的一种重要机制。     

Amplified fragment length polymorphism (AFLP) in soybean: species diversity, inheritance, and near-isogenic line analysis

Amplified fragment length polymorphism (AFLP) analysis is a PCR-based technique capable of detecting more than 50 independent loci in a single PCR reaction. The objectives of the present study were to: (1) assess the extent of AFLP variation in cultivated (Gycine max L. Merr.) and wild soybean (G. soja Siebold & Zucc.), (2) determine genetic relationships among soybean accessions using AFLP data, and (3) evaluate the usefulness of AFLPs as genetic markers. Fifteen AFLP primer pairs detected a total of 759 AFLP fragments in a sample of 23 accessions of wild and cultivated soybean, with an average of 51 fragments produced per primer pair per accession. Two-hundred and seventy four fragments (36% of the total observed) were polymorphic, among which 127 (17%) were polymorphic in G. max and 237 (31%) were polymorphic in G. soja. F₂ segregation analysis of six AFLP fragments indicated that they segregate as stable Mendelian loci. The number of polymorphic loci detected per AFLP primer pair in a sample of 23 accessions ranged from 9 to 27. The AFLP phenotypic diversity values were greater in wild than in cultivated soybean. Cluster and principal component analyses using AFLP data clearly separated G. max and G. soja accessions. Within the G. max group, adapted soybean cultivars were tightly clustered, illustrating the relatively low genetic diversity present in cultivated soybean. AFLP analysis of four soybean near-isogenic lines (NILs) identified three AFLP markers putatively linked to a virus resistance gene from two sources. The capacity of AFLP analysis to detect thousands of independent genetic loci with minimal cost and time requirements makes them an ideal marker for a wide array of genetic investigations.