Introgression between wild and cultivated soybeans of Japan revealed by RFLP analysis for chloroplast DNAs

Wild soybeans collected in Japan were surveyed for RFLPs of chloroplast DNA. Three haplotypes were detected in RFLPs with a cpDNA clone which contains a LSC region adjacent to the left member of IR. Most of the plants tested possessed haplotype III, and a few plants, collected mostly in southern Japan, had haplotype II. Haplotype I, which is the predominant form in modern cultivars, was detected at six sites from four widely separated regions. Our results indicate that haplotype III is predominant in wild soybean of Japan. Some of the plants having haplotype I were phenotypically intermediate between wild and cultivated soybeans, while the others possessed a seed morphology and plant architecture typical of ordinary wild soybean. The plants having haplotype I appear to be either derivatives of hybridization between wild and cultivated soybeans or relics of a direct progenitor of soybean cultivars with the haplotype I chloroplast genome.     

Environmental versus geographical effects on genomic variation in wild soybean (Glycine soja) across its native range in northeast Asia

A fundamental goal in evolutionary biology is to understand how various evolutionary factors interact to affect the population structure of diverse species, especially those of ecological and/or agricultural importance such as wild soybean (Glycine soja). G. soja, from which domesticated soybeans (Glycine max) were derived, is widely distributed throughout diverse habitats in East Asia (Russia, Japan, Korea, and China). Here, we utilize over 39,000 single nucleotide polymorphisms genotyped in 99 ecotypes of wild soybean sampled across their native geographic range in northeast Asia, to understand population structure and the relative contribution of environment versus geography to population differentiation in this species. A STRUCTURE analysis identified four genetic groups that largely corresponded to the geographic regions of central China, northern China, Korea, and Japan, with high levels of admixture between genetic groups. A canonical correlation and redundancy analysis showed that environmental factors contributed 23.6% to population differentiation, much more than that for geographic factors (6.6%). Precipitation variables largely explained divergence of the groups along longitudinal axes, whereas temperature variables contributed more to latitudinal divergence. This study provides a foundation for further understanding of the genetic basis of climatic adaptation in this ecologically and agriculturally important species.     

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.     

The origin and fate of morphological intermediates between wild and cultivated soybeans in their natural habitats in Japan

The spread of transgenes into the genome of wild soybean is a concern when transgenic and wild soybeans are planted sympatrically. The objectives of this study were to investigate the origin and fate of morphological intermediates between wild and cultivated soybeans in their natural habitats in Japan. Twenty nuclear microsatellite and two chloroplast dCAPS markers were used to evaluate genetic variation of 468 wild, 17 intermediate, and 12 cultivated soybean samples collected from six sites between 2003 and 2006. Allelic differentiation of microsatellite markers between wild and cultivated soybeans was sufficient to detect their hybrids. Based on levels of observed heterozygosity, intermediate soybean plants were from two generations: either F₁ or an early segregating generation. Genetic admixture analysis and parentage assignment analysis revealed that the parents of all intermediate soybean plants could be assigned to a particular wild soybean plant and late‐maturing cultivar. The chloroplast DNA haplotypes revealed that all intermediate soybean plants originated from gene flow from cultivated to wild soybeans at all sites. Based on monitoring at both the phenotypic and molecular levels, hybrids quickly disappeared from natural habitats, and secondary gene flow from these plants to wild soybean was not detected. Thus, while gene flow from transgenic soybean into wild soybean can occur, gene introgression appears to be rare in natural habitats in Japan. This is the first report on the detection of gene flow from cultivated to wild soybean at the molecular level.     

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.     

A single origin and moderate bottleneck during domestication of soybean (Glycine max): implications from microsatellites and nucleotide sequences

Background and Aims

It is essential to illuminate the evolutionary history of crop domestication in order to understand further the origin and development of modern cultivation and agronomy; however, despite being one of the most important crops, the domestication origin and bottleneck of soybean (Glycine max) are poorly understood. In the present study, microsatellites and nucleotide sequences were employed to elucidate the domestication genetics of soybean.

Methods

The genomes of 79 landrace soybeans (endemic cultivated soybeans) and 231 wild soybeans (G. soja) that represented the species-wide distribution of wild soybean in East Asia were scanned with 56 microsatellites to identify the genetic structure and domestication origin of soybean. To understand better the domestication bottleneck, four nucleotide sequences were selected to simulate the domestication bottleneck.

Key Results

Model-based analysis revealed that most of the landrace genotypes were assigned to the inferred wild soybean cluster of south China, South Korea and Japan. Phylogeny for wild and landrace soybeans showed that all landrace soybeans formed a single cluster supporting a monophyletic origin of all the cultivars. The populations of the nearest branches which were basal to the cultivar lineage were wild soybeans from south China. The coalescent simulation detected a bottleneck severity of K′ = 2 during soybean domestication, which could be explained by a foundation population of 6000 individuals if domestication duration lasted 3000 years.

Conclusions

As a result of integrating geographic distribution with microsatellite genotype assignment and phylogeny between landrace and wild soybeans, a single origin of soybean in south China is proposed. The coalescent simulation revealed a moderate genetic bottleneck with an effective wild soybean population used for domestication estimated to be ≈2 % of the total number of ancestral wild soybeans. Wild soybeans in Asia, especially in south China contain tremendous genetic resources for cultivar improvement.     

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.     

The possible origin of thick stem in Chinese wild soybean (Glycine soja)

Thicker, erect stem and enlarged seeds are characteristic of the domestication of cultivated soybeans (Glycine max) from its progenitor, wild soybean (G. soja). Wild soybeans have different stem thicknesses but the thick stem as defined here appears in a small number of small-seeded wild soybeans (≤2.0 g/100-seeds) in China. However, little attention has been paid to this phenomenon in considering the origin of thick stem in wild soybean. Here, we addressed this question through the study of a mixed field of wild, semi-wild and cultivated soybeans. Thick-stemmed samples had lower sensitivity to light period, higher mean genetic diversity (H _e = 0.090, H = 0.535) and higher mean multilocus outcrossing rate (t _m = 9.77 %), while thin-stemmed plants were the opposite (H _e = 0.029, H = 0.416) and lower mean outcrossing rate (t _m = 5.88 %). F statistics calculations indicated that there was genetic differentiation between the thin and thick stems. UPGM cluster analysis showed that not only were thick-stemmed wild soybeans genetically different from thin-stemmed ones, but they were also genetically closer to semi-wild soybean, to varying degrees completely dependent on seed size. These data strongly implied that the plants with thick stems had more complicated genetic backgrounds than the thin-stemmed ones, and that they were related to cultivated soybeans. This study suggests that if plants have distinctly thick stems (an average 2.5-fold thicker than other thin-stemmed plants) or stems similar to semi-wild plants and/or near to local soybeans in a natural wild population adjacent to farmlands, such plants could be cultivar-introgressive offspring.     

Genome-wide signatures of the geographic expansion and breeding of soybean

Soybean is a leguminous crop that provides oil and protein. Exploring the genomic signatures of soybean evolution is crucial for breeding varieties with improved adaptability to environmental extremes. We analyzed the genome sequences of 2,214 soybeans and proposed a soybean evolutionary route, i.e., the expansion of annual wild soybean (Glycine soja Sieb. & Zucc.) from southern China and its domestication in central China, followed by the expansion and local breeding selection of its landraces (G. max (L.) Merr.). We observed that the genetic introgression in soybean landraces was mostly derived from sympatric rather than allopatric wild populations during the geographic expansion. Soybean expansion and breeding were accompanied by the positive selection of flowering time genes, including GmSPA3c. Our study sheds light on the evolutionary history of soybean and provides valuable genetic resources for its future breeding.

     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Study of interactions between wild soybean subpopulations (<Emphasis Type="Italic">Glycine soja</Emphasis>) in the valley of the Tsukanovka river in the south of Russia’s Far East

Comparative analysis of the genetic structure of natural and anthropogenic populations of G. soja gives significant information about the formation of different populations and allows for the developing of measures for the preservation of the unique natural gene bank of wild soybean, which is a species closely related to cultivated soybean. In this study, ISSR-markers were used to carry out a comparative analysis of the genetic structure of natural and anthropogenic subpopulations of G. soja for studying the possible mutual influence of subpopulations of anthropogenic and natural phytocenosis on the formation of their genetic diversity and studying the genetic structure of natural subpopulations of wild soybean in the contact places between the two types of cenoses. As a result, the characteristics that describe the genetic diversity of the studied populations have been identified, and the important role of interaction between subpopulations of different phytocenoses in the formation of the spatial genetic structure of the population in the Tsukanovka river valley have been demonstrated.     

Ribosomal gene variation in soybean (Glycine) and its relatives

Summary The genes encoding the 18S25S ribosomal RNA gene repeat in soybean (Glycine max) and its relatives in the genus Glycine are surveyed for variation in repeat length and restriction enzyme site locations. Within the wild species of subgenus Glycine, considerable differences in repeat size occur, with a maximum observed in G. falcata. Repeat length and site polymorphisms occur in several species, but within individual plants only single repeat types are observed. The rDNA of the cultivated soybean and its wild progenitor, G. soja are identical at the level of this study, and no variation is found in over 40 accessions of the two species. Data from rDNA mapping studies are congruent with those of previous biosystematic studies, and in some instances give evidence of divergences not seen with other approaches.     

Distribution of restriction site polymorphism within the chloroplast genome of the genus Glycine,subgenus Soja

Summary Restriction fragment length polymorphisms (RFLPs) have been used to detect intragenic sequence diversity in Glycine subgenus soja chloroplast DNA. The distribution of these RFLPs allow Glycine max and G. soja accessions to be grouped according to cytoplasmic genetic relatedness. DNA clones from mung bean chloroplast DNA were used to locate the RFLPs to specific regions of the chloroplast genome. In the course of the experiments, several previously unobserved RFLPs were also identified. At least six molecular changes were detected, including both restriction site loss or gain and insertion/deletion events. Three of the fragment polymorphisms detected are due to changes in the juncture region between one inverted repeat region and the large single-copy region. Probes detecting polymorphisms in three representative soybean genotypes were used to screen additional cultivars and Plant Introductions. The distribution of RFLP patterns in these accessions were consistent with the patterns of previously described cytoplasmic groupings, with the exception of one accession, which formed a new plastome group.