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.     

Uncovering signatures of selection in the soybean genome using SSR diversity near QTLs of agronomic importance

The cultivated soybean [Glycine max (L.) Merr.] is widely considered to descend from the wild soybean (G. soja Sieb. & Zucc.). This study was designed to evaluate the genetic variability and differentiation between G. soja and G. max, and to detect signatures of the selection that may have occurred during the domestication process from G. soja to G. max. A total of 192 G. soja accessions and 104 G. max accessions were genotyped using eight selected simple sequence repeat (SSR) markers assigned to three SSR groups. Four SSRs in group A were not located near any known QTL. Three SSRs in group B were associated with seed protein content, and an SSR in group C was associated with resistance to Sclerotinia stem rot. The number of alleles per locus and the level of genetic variability in G. soja were higher than those in G. max. A total of 122 out of 125 alleles were present in G. soja, but only 59 alleles were detected in G. max. The average gene diversity was 0.74 in G. soja and 0.64 in G. max. Four SSRs near QTLs of agronomic importance showed strong genetic differentiation and shift change in high frequency alleles in groups B and C between G. soja and G. max, revealing selection signatures that may reflect the domestication events and recent selective breeding. With reduced diversity in G. max, some undomesticated genes from G. soja should be prime candidates for introgression to increase the pool of diversity in G. max.     

Study of the biosafety of genetically modified soybean in the center of its origin and diversity in the Far East of the Russian Federation

Wild soybean (Glycine soja Sieb. and Zucc.) is the nearest relative of cultivated soybean (G. max (L.) Merr.). Study of the population genetic structure of wild-growing relatives of genetically modified (GM) plants in the centers of their origin is one of the main procedures before introduction of GM crops in these areas. We studied the genetic variability of nine wild-growing soya populations of Primorskii krai using RAPD analysis. The level of genetic variability of G. soja was considerably higher than that of G. max. We analyzed phylogenetic relationships in the genus Glycine subgenus Soja using RAPD markers. Our data confirm the validity of allocation to G. gracilis of the rank of species.     

Selection of a core set of RILs from Forrest × Williams 82 to develop a framework map in soybean

Soybean BAC-based physical maps provide a useful platform for gene and QTL map-based cloning, EST mapping, marker development, genome sequencing, and comparative genomic research. Soybean physical maps for “Forrest” and “Williams 82” representing the southern and northern US soybean germplasm base, respectively, have been constructed with different fingerprinting methods. These physical maps are complementary for coverage of gaps on the 20 soybean linkage groups. More than 5,000 genetic markers have been anchored onto the Williams 82 physical map, but only a limited number of markers have been anchored to the Forrest physical map. A mapping population of Forrest × Williams 82 made up of 1,025 F₈ recombinant inbred lines (RILs) was used to construct a reference genetic map. A framework map with almost 1,000 genetic markers was constructed using a core set of these RILs. The core set of the population was evaluated with the theoretical population using equality, symmetry and representativeness tests. A high-resolution genetic map will allow integration and utilization of the physical maps to target QTL regions of interest, and to place a larger number of markers into a map in a more efficient way using a core set of RILs.     

Chromosomal features revealed by comparison of genetic maps of Glycine max and Glycine soja

Recombination is a crucial component of evolution and breeding. New combinations of variation on chromosomes are shaped by recombination. Recombination is also involved in chromosomal rearrangements. However, recombination rates vary tremendously among chromosome segments. Genome-wide genetic maps are one of the best tools to study variation of recombination. Here, we describe high density genetic maps of Glycine max and Glycine soja constructed from four segregating populations. The maps were used to identify chromosomal rearrangements and find the highly predictable pattern of cross-overs on the broad scale in soybean. Markers on these genetic maps were used to evaluate assembly quality of the current soybean reference genome sequence. We find a strong inversion candidate larger than 3 Mb based on patterns of cross-overs. We also identify quantitative trait loci (QTL) that control number of cross-overs. This study provides fundamental insights relevant to practical strategy for breeding programs and for pan-genome researches.     

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.     

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.     

Structural variation of mitochondrial genomes sheds light on evolutionary history of soybeans

The architecture and genetic diversity of mitogenome (mtDNA) are largely unknown in cultivated soybean (Glycine max), which is domesticated from the wild progenitor, Glycine soja, 5000 years ago. Here, we de novo assembled the mitogenome of the cultivar ‘Williams 82’ (Wm82_mtDNA) with Illumina PE300 deep sequencing data, and verified it with polymerase chain reaction (PCR) and Southern blot analyses. Wm82_mtDNA maps as two autonomous circular chromosomes (370 871-bp Chr-m1 and 62 661-bp Chr-m2). Its structure is extensively divergent from that of the mono-chromosomal mitogenome reported in the landrace ‘Aiganhuang’ (AGH_mtDNA). Synteny analysis showed that the structural variations (SVs) between two genomes are mainly attributed to ectopic and illegitimate recombination. Moreover, Wm82_mtDNA and AGH_mtDNA each possess six and four specific regions, which are absent in their counterparts and likely result from differential sequence-loss events. Mitogenome SV was further studied in 39 wild and 182 cultivated soybean accessions distributed world-widely with PCR/Southern analyses or a comparable in silico analysis. The results classified both wild and cultivated soybeans into five cytoplasmic groups, named as GSa–GSe and G1–G5; ‘Williams 82’ and ‘Aiganhuang’ belong to G1 and G5, respectively. Notably, except for members in GSe and G5, all accessions carry a bi-chromosomal mitogenome with a common Chr-m2. Phylogenetic analyses based on mtDNA structures and chloroplast gene sequences both inferred that G1–G3, representing >90% of cultigens, likely inherited cytoplasm from the ancestor of domestic soybean, while G4 and G5 likely inherited cytoplasm from wild soybeans carrying GSa- and GSe-like cytoplasm through interspecific hybridization, offering new insights into soybean cultivation history.     

Identification of putative QTL that underlie yield in interspecific soybean backcross populations

Glycine soja, the wild progenitor of soybean, is a potential source of useful genetic variation in soybean improvement. The objective of our study was to map quantitative trait loci (QTL) from G. soja that could improve the crop. Five populations of BC₂F₄-derived lines were developed using the Glycine max cultivar IA2008 as a recurrent parent and the G. soja plant introduction (PI) 468916 as a donor parent. There were between 57 and 112 BC₂F₄-derived lines in each population and a total of 468 lines for the five populations. The lines were evaluated with simple sequence repeat markers and in field tests for yield, maturity, plant height, and lodging. The field testing was done over 2 years and at two locations each year. Marker data were analyzed for linkage and combined with field data to identify QTL. Using an experimentwise significance threshold of P=0.05, four yield QTL were identified across environments on linkage groups C2, E, K, and M. For these yield QTL, the IA2008 marker allele was associated with significantly greater yield than the marker allele from G. soja. In addition, one lodging QTL, four maturity QTL, and five QTL for plant height were identified across environments. Of the 14 QTL identified, eight mapped to regions where QTL with similar effects were previously mapped. Many regions carrying the yield QTL were also significant for other traits, such as plant height and lodging. When the significance threshold was reduced and the data were analyzed with simple linear regression, four QTL with a positive allele for yield from G. soja were mapped. One epistatic interaction between two genetic regions was identified for yield using an experimentwise significance threshold of P=0.05. Additional research is needed to establish whether multiple trait associations are the result of pleiotropy or genetic linkage and to retest QTL with a positive effect from G. soja.Communicated by H.C. Becker     

The salt resistance of wild soybean (Glycine soja Sieb. et Zucc. ZYD 03262) under NaCl stress is mainly determined by Na+ distribution in the plant

Understanding the salt resistance mechanism of wild soybean is important in improving salt tolerance of cultivated soybean. Therefore, we comparatively analyzed effects of NaCl on photosynthesis, antioxidant enzyme activity, and ion distribution in a cultivated (Glycine. max) and a wild (Glycine soja) soybean to study the salt resistance mechanism of the G. soja. The results showed that more Na⁺ was accumulated in the G. soja roots than in the G. max roots, but the Na⁺ in the G. soja leaves was much less than that observed in the G. max leaves. The Na⁺ concentrations in the G. soja leaves were not high enough to affect the photosynthetic apparatus, which was demonstrated by less inhibition of photosynthetic activity, stomatal conductance, carboxylation efficiency in the G. soja leaves than in the G. max leaves after treated with different concentrations of NaCl. Meanwhile, there were no significant changes in intercellular CO₂ concentration, maximum PSII quantum yield, and relative water content in the G. soja leaves after NaCl treatment, while they significantly decreased in the G. max leaves. The non-photochemical quenching and the activities of superoxide dismutase (EC 1.15.1.1) and ascorbate peroxidase (EC 1.11.1.11) in the G. soja leaves increased with the increasing of NaCl concentrations, whereas only the activity of superoxide dismutase increased in G. max leaves. Based on these results, we suggested that the G. soja is able to accumulate higher levels of Na⁺ in its roots, and prevent the transportation of Na⁺ to leaves to protect photosynthetic apparatus from salt damage.     

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.     

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.     

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.     

Inheritance of polymorphic markers generated by DNA amplification fingerprinting and their use as genetic markers in soybean

DNA amplification fingerprinting (DAF) using a high primer-to-template ratio and single, very short arbitrary primers, was used to generate amplified fragment length polymorphic markers (AFLP) in soybean (Glycine max (L.) Merr.). The inheritance of AFLPs was studied using a cross between the ancestral Glycine soja PI468.397 and Glycine max (L.) Merr. line nts382, F1 and F2 progeny. The amplification reaction was carried out with soybean genomic DNA and 8 base long oligounucleotide primers. Silver-stained 5% polyacrylamide gels containing 7 M urea detected from 11 to 28 DAF products with primers of varying GC content (ranging from 50 to 100% GC). Depending on their intensity, AFLPs were classified into three classes. DAF profiles were reproducible for different DNA extractions and gels. Forty AFLPs were detected by 26 primers when comparing G. soja and G. max. Most AFLPs were inherited as dominant Mendelian markers in F1 and F2 populations. However, abnormal inheritance occured with about 25% of polymorphisms. One marker was inherited as a maternal marker, presumably originating from organelle DNA while another showed apparent paternal inheritance. To confirm the nuclear origin and utility of dominant Mendelian markers, three DAF polymorphisms were mapped using a F11 mapping population of recombinant inbred lines from soybean cultivars Minsoy × Noir 1. The study showed that DAF-generated polymorphic markers occur frequently and reliably, that they are inherited as Mendelian dominant loci and that they can be used in genome mapping.     

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.     

Genetic variants in root architecture-related genes in a Glycine soja accession,a potential resource to improve cultivated soybean

Background

Root system architecture is important for water acquisition and nutrient acquisition for all crops. In soybean breeding programs, wild soybean alleles have been used successfully to enhance yield and seed composition traits, but have never been investigated to improve root system architecture. Therefore, in this study, high-density single-feature polymorphic markers and simple sequence repeats were used to map quantitative trait loci (QTLs) governing root system architecture in an inter-specific soybean mapping population developed from a cross between Glycine max and Glycine soja.

Results

Wild and cultivated soybean both contributed alleles towards significant additive large effect QTLs on chromosome 6 and 7 for a longer total root length and root distribution, respectively. Epistatic effect QTLs were also identified for taproot length, average diameter, and root distribution. These root traits will influence the water and nutrient uptake in soybean. Two cell division-related genes (D type cyclin and auxin efflux carrier protein) with insertion/deletion variations might contribute to the shorter root phenotypes observed in G. soja compared with cultivated soybean. Based on the location of the QTLs and sequence information from a second G. soja accession, three genes (slow anion channel associated 1 like, Auxin responsive NEDD8-activating complex and peroxidase), each with a non-synonymous single nucleotide polymorphism mutation were identified, which may also contribute to changes in root architecture in the cultivated soybean. In addition, Apoptosis inhibitor 5-like on chromosome 7 and slow anion channel associated 1-like on chromosome 15 had epistatic interactions for taproot length QTLs in soybean.

Conclusion

Rare alleles from a G. soja accession are expected to enhance our understanding of the genetic components involved in root architecture traits, and could be combined to improve root system and drought adaptation in soybean.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1334-6) contains supplementary material, which is available to authorized users.     

Genetic mapping of a novel gene for soybean aphid resistance in soybean (Glycine max [L.] Merr.) line P203 from China

The soybean aphid (SA: Aphis glycines Matsumura) is a worldwide pest of soybean (Glycine max [L.] Merr.). The objectives of this study were to identify the type of aphid resistance and the resistance phenotype in soybean line ‘P203’, and to map the relative position of the gene involved. Compared with cultivars ‘P746’ and ‘Dongnong 47’, P203 was demonstrated to possess antixenosis resistance. P203 prevented aphids from reproducing in a choice test, but the resistance level decreased significantly in a no-choice test at 11 and 21 days after infestation. Analysis of 273 Dongnong 47/P203 F₂ plants and confirmed using 260 F_2:3 families revealed that a single dominant gene from P203 was positioned between marker loci Sat_377 and Satt409 on chromosome 8. The gene was further mapped to a 1.57 Mb interval flanked by marker loci BARCSOYSSR_08_1451 and BARCSOYSSR_08_1527. We developed five new SSR markers in the target interval and the resistance locus mapped between new markers SSR_08_75 and SSR_08_88 with the genetic distance of 1.1 and 1.0 cM corresponding to a physical distance of 192 kb on the Williams 82 8X draft genome assembly (Glyma1.01). A single serine/threonine protein kinase gene is present in this region, suggesting that the SA resistance mechanism in P203 may be different from those previously reported. Therefore, the resistance gene could very well be novel, and could be valuable in soybean aphid resistance breeding programs.     

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 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.