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.     

Restriction fragment length polymorphism diversity in soybean

Summary Fifty-eight soybean accessions from the genus Glycine, subgenus Soja, were surveyed with 17 restriction fragment length polymorphism (RFLP) genetic markers to assess the level of molecular diversity and to evaluate the usefulness of previously identified RFLP markers. In general, only low levels of molecular diversity were observed: 2 of the 17 markers exhibited three alleles per locus, whereas all others had only two alleles. Thirty-five percent of the markers had rare alleles present in only 1 or 2 of the 58 accessions. Molecular diversity was least among cultivated soybeans and greatest between accessions of different soybean species such as Glycine max (L.) Merr. and G. soja Sieb. and Zucc. Principal component analysis was useful in reducing the multidimensional genotype data set and identifying genetic relationships.     

Genetic diversity and peculiarity of annual wild soybean (<Emphasis Type="Italic">G. soja</Emphasis> Sieb. et Zucc.) from various eco-regions in China

Annual wild soybean (Glycine soja Sieb. et Zucc.) is believed to be a potential gene source for future soybean improvement in coping with the world climate change for food security. To evaluate the wild soybean genetic diversity and differentiation, we analyzed allelic profiles at 60 simple-sequence repeat (SSR) loci and variation of eight morph-biological traits of a representative sample with 196 accessions from the natural growing area in China. For comparison, a representative sample with 200 landraces of Chinese cultivated soybean was included in this study. The SSR loci produced 1,067 alleles (17.8 per locus) with a mean gene diversity of 0.857 in the wild sample, which indicated the genetic diversity of G. soja was much higher than that of its cultivated counterpart (total 826 alleles, 13.7 per locus, mean gene diversity 0.727). After domestication, the genetic diversity of the cultigens decreased, with its 65.5% alleles inherited from the wild soybean, while 34.5% alleles newly emerged. AMOVA analysis showed that significant variance did exist among Northeast China, Huang-Huai-Hai Valleys and Southern China subpopulations. UPGMA cluster analysis indicated very significant association between the geographic grouping and genetic clustering, which demonstrated the geographic differentiation of the wild population had its relevant genetic bases. In comparison with the other two subpopulations, the Southern China subpopulation showed the highest allelic richness, diversity index and largest number of specific-present alleles, which suggests Southern China should be the major center of diversity for annual wild soybean. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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.     

The comparison of RFLP,RAPD, AFLP and SSR (microsatellite) markers for germplasm analysis

The utility of RFLP (restriction fragment length polymorphism), RAPD (random-amplified polymorphic DNA), AFLP (amplified fragment length polymorphism) and SSR (simple sequence repeat, microsatellite) markers in soybean germplasm analysis was determined by evaluating information content (expected heterozygosity), number of loci simultaneously analyzed per experiment (multiplex ratio) and effectiveness in assessing relationships between accessions. SSR markers have the highest expected heterozygosity (0.60), while AFLP markers have the highest effective multiplex ratio (19). A single parameter, defined as the marker index, which is the product of expected heterozygosity and multiplex ratio, may be used to evaluate overall utility of a marker system. A comparison of genetic similarity matrices revealed that, if the comparison involved both cultivated (Glycine max) and wild soybean (Glycine soja) accessions, estimates based on RFLPs, AFLPs and SSRs are highly correlated, indicating congruence between these assays. However, correlations of RAPD marker data with those obtained using other marker systems were lower. This is because RAPDs produce higher estimates of interspecific similarities. If the comparisons involvedG. max only, then overall correlations between marker systems are significantly lower. WithinG. max, RAPD and AFLP similarity estimates are more closely correlated than those involving other marker systems.Abbreviations RFLP restriction fragment length plymorphism - RAPD random-amplified polymorphic DNA - AFLP amplified fragment length polymorphism - SSR simple sequence repeat - PCR polymerase chain reaction - TBE Tris-borate-EDTA buffer - MI marker index - SENA sum of effective numbers of alleles     

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.     

QTL affecting fitness of hybrids between wild and cultivated soybeans in experimental fields

The objective of this study was to identify quantitative trait loci (QTL) affecting fitness of hybrids between wild soybean (Glycine soja) and cultivated soybean (Glycine max). Seed dormancy and seed number, both of which are important for fitness, were evaluated by testing artificial hybrids of G. soja × G. max in a multiple‐site field trial. Generally, the fitness of the F₁ hybrids and hybrid derivatives from self‐pollination was lower than that of G. soja due to loss of seed dormancy, whereas the fitness of hybrid derivatives with higher proportions of G. soja genetic background was comparable with that of G. soja. These differences were genetically dissected into QTL for each population. Three QTLs for seed dormancy and one QTL for total seed number were detected in the F₂ progenies of two diverse cross combinations. At those four QTLs, the G. max alleles reduced seed number and severely reduced seed survival during the winter, suggesting that major genes acquired during soybean adaptation to cultivation have a selective disadvantage in natural habitats. In progenies with a higher proportion of G. soja genetic background, the genetic effects of the G. max alleles were not expressed as phenotypes because the G. soja alleles were dominant over the G. max alleles. Considering the highly inbreeding nature of these species, most hybrid derivatives would disappear quickly in early self‐pollinating generations in natural habitats because of the low fitness of plants carrying G. max alleles.     

Molecular characterization and genetic diversity analysis of soybean (Glycine max (L.) Merr.) germplasm accessions in India

Molecular characterization and genetic diversity among 82 soybean accessions was carried out by using 44 simple sequence repeat (SSR) markers. Of the 44 SSR markers used, 40 markers were found polymorphic among 82 soybean accessions. These 40 polymorphic markers produced a total of 119 alleles, of which five were unique alleles and four alleles were rare. The allele number for each SSR locus varied between two to four with an average of 2.97 alleles per marker. Polymorphic information content values of SSRs ranged from 0.101 to 0.742 with an average of 0.477. Jaccard’s similarity coefficient was employed to study the molecular diversity of 82 soybean accessions. The pairwise genetic similarity among 82 soybean accessions varied from 0.28 to 0.90. The dendrogram constructed based on genetic similarities among 82 soybean accessions identified three major clusters. The majority of genotypes including four improved cultivars were grouped in a single subcluster IIIa of cluster III, indicating high genetic resemblance among soybean germplasm collection in India.

Electronic supplementary material

The online version of this article (doi:10.1007/s12298-014-0266-y) contains supplementary material, which is available to authorized users.     

An attempt to elucidate the origin of cultivated soybean via comparison of nucleotide sequences encoding glycinin B4 polypeptide of cultivated soybean,Glycine max,and its presumed wild progenitor,Glycine soja

Summary Nucleotide sequences of cDNAs encoding soybean glycinin B₄ polypeptide were compared for three soybean cultivars and two introductions of wild soybean, G. soja. For three G. max cultivars, only two nucleotide substitutions were found, while G. max and G. soja nucleotide sequences had four substitutions. These data give added proof that G. max originated from G. soja. On the other hand, the time required for the accumulation of four nucleotide substitutions (calculated from the parameters of 11S globulin molecular evolution) appeared to be longer than the duration of the soybean domestication period.     

Evaluation of the geographical pattern of genetic diversity of Glycine soja and Glycine max based on four single copy nuclear gene loci: For conservation of soybean germplasm

Genetic diversity and its geographical patterns play a very important role in species conservation and exploitation. Here, nucleotide polymorphism patterns of four single copy nuclear gene loci in wild (Glycine soja) and cultivated soybean (Glycine max) populations from different geographical regions as well as their demographic history were analyzed. The results showed that: (1) Southern subpopulation has the highest, while central subpopulation revealed the lowest genetic diversity among three Chinese G. soja subpopulations. (2) Northern Chinese G. max subpopulation depicted higher genetic diversity than other two Chinese, Korean, Japanese and American G. max subpopulations. (3) Significant genetic differentiation (P < 0.001) was observed among Chinese G. soja subpopulations from three ecological zones. There was also a significant genetic differentiation（P < 0.01）between three Chinese and Japanese subpopulations of G. max. (4) The demographic dynamics revealed that effective population size of G. soja is expanding, while it was constant in G. max. G. soja is a useful germplasm resource to widen the genetic diversity of G. max. This study suggests that native populations of G. soja from different geo-ecological regions should be protected to conserve the genetic diversity.     

Evaluation of Genetic Diversity in Chinese Wild Apple Species Along with Apple Cultivars Using SSR Markers

China, one of the primary centers of genetic diversity for the genus Malus, is very rich in wild apple germplasm. In this study, genetic diversity in 29 Malus accessions, including 12 accessions from 7 Chinese Malus species, 4 Chinese landraces, and 13 introduced apple cultivars, was assessed using a set of 19 single-locus simple sequence repeat (SSR) markers distributed across all 17 linkage groups of the apple genome. The number of alleles detected at each locus ranged from 2 to 11, with an average of 5.3 per SSR marker. In some accessions, 16 unique alleles were identified. Ten out of these 16 unique alleles (62.5%) were detected exclusively in wild species, indicating that these Chinese wild apple species have considerable genetic diversity and can be used in breeding programs to increase the genetic diversity of apple cultivars. Using 19 SSRs, an unweighted pair-group method with arithmetic average cluster analysis was conducted, and the resulting dendrogram revealed that all cultivars, except for E??peMeBckoe, were clustered together in the same group. The Russian cultivar E??peMeBckoe was closely related to the Chinese crabapple Baihaitang (M. prunifolia), with a high similarity coefficient value of 0.94. Of the two M. sieversii accessions used, one accession showed a close relationship to apple cultivars, while the other accession was closely related to wild apple species, suggesting the presence of a wider genetic diversity in Chinese M. sieversii species. The influence of SSR marker selection on genetic diversity analysis in this Malus collection was also discussed.     

Microsatellite diversity in tetraploid <Emphasis Type="Italic">Gossypium</Emphasis> germplasm: assembling a highly informative genotyping set of cotton SSRs

A series of 320 mapped simple sequence repeats (SSRs) have been used to screen the allelic diversity of tetraploid Gossypium species. Fourty-seven genotypes were analyzed representing (i) the wide spectrum of diversity of the cultivated pool and of the primitive landraces of species G. hirsutum (‘marie-galante’, ‘punctatum’, ‘richmondi’, ‘morrilli’, ‘palmeri’, and ‘latifolium’, and ‘yucatanense’), and (ii) species G. barbadense, G. darwinii and G. tomentosum. The polymorphism of 201 SSR loci revealed 1128 allelic variants ranging from 3 to 17 per locus. Neighbor-joining (NJ) method based on genetic dissimilarities produced groupings consistent with the assignments of accessions both at species and at race level. Our data confirmed the proximity of the Galapagos endemic species G. darwinii to species G. barbadense. Within species G. hirsutum, and as compared to the other 6 races, race yucatanense appeared as the most distant from cultivated genotypes. Race yucatanense also exhibited the highest number of unique alleles. The important informative heterogeneity of the 201 SSR loci was exploited to select the most polymorphic ones that were assembled into three series of genome-wide (i.e. each homoeologous AD chromosome pair being equally represented) and mutliplexable (× 3) SSRs. Using one of these ‘genotyping set’, consisting of 39 SSRs (one 3-plex for each of the 13 AD chromosomes pairs) or 45 loci, we were able to assess the relationships between accessions and the topology in the genetic diversity sampled. Such genotyping set of highly informative SSR markers assembled in PCR-multiplex, while increasing genotyping throughput, will be applicable for molecular genetic diversity studies of large germplasm collections. Electronic Supplementary Material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Exploration of presence/absence variation and corresponding polymorphic markers in soybean genome

This study was designed to reveal the genome‐wide distribution of presence/absence variation (PAV) and to establish a database of polymorphic PAV markers in soybean. The 33 soybean whole‐genome sequences were compared to each other with that of Williams 82 as a reference genome. A total of 33,127 PAVs were detected and 28,912 PAV markers with their primer sequences were designed as the database NJAUSoyPAV_1.0. The PAVs scattered on whole genome while only 518 (1.8%) overlapped with simple sequence repeats (SSRs) in BARCSOYSSR_1.0 database. In a random sample of 800 PAVs, 713 (89.13%) showed polymorphism among the 12 differential genotypes. Using 126 PAVs and 108 SSRs to test a Chinese soybean germplasm collection composed of 828 Glycine soja Sieb. et Zucc. and Glycine max (L.) Merr. accessions, the per locus allele number and its variation appeared less in PAVs than in SSRs. The distinctness among alleles/bands of PCR (polymerase chain reaction) products showed better in PAVs than in SSRs, potential in accurate marker‐assisted allele selection. The association mapping results showed SSR + PAV was more powerful than any single marker systems. The NJAUSoyPAV_1.0 database has enriched the source of PCR markers, and may fit the materials with a range of per locus allele numbers, if jointly used with SSR markers.     

Nucleotide diversity patterns of three divergent soybean populations: evidences for population‐dependent linkage disequilibrium and taxonomic status of Glycine gracilis

The level of linkage disequilibrium (LD) is a major factor to determine DNA polymorphism pattern of a population and to construct high‐resolution maps useful in localizing and gene cloning of complicated traits. Here, we investigated LD level of three soybean populations with different genetic backgrounds and taxonomic status of G. gracilis by comparing the DNA polymorphism patterns of four high‐diversity single‐copy nuclear genes. A total of 152, 22, and 77 accessions of G. soja, G. gracilis, and G. max were observed. The results indicated that G. max retained only 75.3 (π) and 39% (θ) of the nucleotide polymorphism found in G. soja. Four gene loci evolved according to neutrality in both G. max and G. gracilis populations, and three gene loci evolved according to neutrality in G. soja population by Tajima's and Fu and Li's test. However, one gene locus deviated from neutrality by Fu and Li's test in the G. soja population. Further, medial level of LD (average r² = 0.2426) was found in intragene in G. max and G. gracilis populations, but unexpected low level of LD (r² ≤ 0.0539) was found in G. soja population. Significant genetic differentiation was detected between G. max and G. soja populations and also between G. max and G. gracilis populations; however, nonsignificant genetic differentiation was found between G. gracilis and G. soja populations. The results suggest that LD level depends on genetic background of soybean population, and implicit that G. gracilis should be regarded as the variant of G. soja, not as an independent species.     

Identification of <Emphasis Type="Italic">GmSALT3</Emphasis> haplotypes and development of molecular markers based on their diversity associated with salt tolerance in soybean

Soybean [Glycine max (L.) Merr.] is a major agricultural crop and generally known as a salt-sensitive crop. In a previous study, GmSALT3 was identified as a salt tolerance gene in soybean, and its nine haplotypes (H1 to H9) were reported in Chinese soybean accessions. In the present study, we aimed to identify new haplotypes of GmSALT3 in soybean accessions and to develop molecular markers for selection of salt-tolerant and -sensitive accessions. To do so, we examined genomic variations in the GmSALT3 coding region of 216 accessions of G. max and G. soja from Korea, China, and Japan. As a result, 40 different haplotypes, including three known haplotypes (H1, H2, and H5), were identified. Performing salt tolerance tests for the haplotypes, we were able to classify them into salt-tolerant (8) and salt-sensitive (32) categories. We also found that more variations in the haplotype composition of G. soja accessions exist than in that of G. max accessions. Quantitative expression analysis showed that almost all of the salt-tolerant haplotypes had much higher levels of GmSALT3 expression than the salt-sensitive haplotypes did. Finally, we developed molecular markers and applied them to screen salt tolerance of soybean accessions. The molecular markers performed well with an accuracy of 98.8% in identifying phenotypes of soybean accessions.     

Assessment of molecular diversity and evolutionary relationship of kenaf (Hibiscus cannabinus L.), roselle (H. sabdariffa L.) and their wild relatives

Kenaf (Hibiscus cannabinus L.) and roselle (H. sabdariffa L.) are valuable fibre crop species with diverse end use. Phylogenetic relationship of 73 accessions of kenaf, roselle and their wild relatives from 15 countries was assessed using 44 inter-simple sequence repeat (ISSR) and jute (Corchorus olitorius L.) specific simple sequence repeats (SSR) markers. A total of 113 alleles were identified of which 61.95 % were polymorphic. Jute specific SSR markers exhibited high polymorphism and resolving power in kenaf, although ISSR markers exhibited higher resolving power than SSR markers. Number of polymorphic alleles varied from 1 to 5 for ISSR and 1 to 6 for SSR markers. Cultivated species exhibited higher allele polymorphism (57 %) than the wild species (35 %), but the improved cultivars exhibited lower genetic diversity compared to germplasm accessions. Accessions with common genetic lineage and geographical distribution clustered together. Indian kenaf varieties were distinct from cultivars bred in other countries and shared more genetic homology with African accessions. High genetic diversity was observed in the Indian (J = 0.35–0.74) and exotic kenaf germplasm collections (J = 0.38–0.79), suggesting kenaf might have been introduced in India from Africa through Central Asia during early domestication. Genetic similarity-based cluster analysis was in close accordance with taxonomic classification of Hibiscus.     

Analysis of Genetic Diversity in the North Eastern Himalayan Maize Landraces using Microsatellite Markers

Population DNA fingerprinting of 48 selected North Eastern Himalayan (NEH) landrace accessions was undertaken using 41 polymorphic fluorescent dye-labelled microsatellite/Simple Sequence Repeat (SSR) markers, using a DNA Sequencer. The analysis revealed a large number of SSR alleles (576), with high mean number of alleles per locus (13.8), and Polymorphism Information Content (PIC) of 0.63, reflecting the level of diversity in the NEH accessions and the informativeness of the SSR markers. The study also led to identification of 135 unique alleles, differentiating 44 out of the 48 accessions. Five highly frequent (major) SSR alleles (umc1545 _80bp, phi062 _162bp, umc1367 _159bp, umc2250 _152bp and phi112 _152bp) were detected indicating that chromosomal regions harbouring these S SR alleles might not be selectively neutral. Analysis of population genetic parameters, including Wright’s F statistics, revealed high level of genetic differentiation, very low levels of inbreeding, and restricted gene flow between the NEH landraces. AMOVA (Analysis of Molecular Variance) showed that 67 per cent of the total variation in the accessions could be attributed to within-population diversity, and the rest between the accessions. Cluster analysis of SSR data using Rogers’ genetic distance and UPGMA, showed significant genetic diversity among the landraces from Sikkim. This is the first detailed study of SSR allele frequency-based analysis of genetic diversity in the NEH maize landraces of India.     

Microsatellite marker-based diversity and population genetic analysis of selected lowland and mid-altitude maize landrace accessions of India

Maize (Zea mays L.) harbours significant genetic diversity not only in its centre of origin (Mexico) but also in several countries worldwide, including India, in the form of landraces. In this study, DNA fingerprinting of 48 landrace accessions from diverse regions of India was undertaken using 42 fluorescent dye-labeled Simple Sequence Repeat (SSR) markers, followed by allele resolution using DNA sequencer and analysis of molecular diversity within and among these landraces. The study revealed a large number of alleles (550), with high mean number of alleles per locus (13.1), and Polymorphism Information Content (PIC) of 0.60, reflecting the level of diversity in the landrace accessions. Besides identification of 174 unique alleles in 44 accessions, six highly frequent SSR alleles were detected at six loci (phi014, phi090, phi112, umc1367, phi062 and umc1266) with individual frequencies greater than 0.75, indicating that chromosomal regions harboring these SSR alleles are not selectively neutral. F statistics revealed very high genetic differentiation, population subdivision and varying levels of inbreeding in the landraces. Analysis of Molecular Variance showed that 63 % of the total variation in the accessions could be attributed to within-population diversity, and 37 % represented between population diversity. Cluster analysis of SSR data using Nei’s genetic distance and UPGMA revealed considerable genetic diversity in these populations, although no clear separation of accessions was observed based on their geographic origin.