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.     

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.

     

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     

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.     

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.     

Two Erysiphe species associated with recent outbreak of soybean powdery mildew: results of molecular phylogenetic analysis based on nuclear rDNA sequences

 Serious outbreaks of powdery mildew by a fungus belonging to the mitosporic genus Oidium subgenus Pseudoidium have been reported on soybean (Glycine max) in a wide area of eastern Asia since 1998. The taxonomic and phylogenetic placement of the causal fungus has not yet been determined because of lack of the perfect stage. We found ascomata having mycelioid appendages on a single leaf of soybean infested by powdery mildew. Molecular phylogenetic analysis was conducted based on a total of 14 sequences of the rDNA internal transcribed spacer (ITS) region from 13 soybean and wild soybean (Glycine soja) materials collected in Japan, Korea, Vietnam, and the United States, combined with 47 sequence data obtained from the DNA databases. It was revealed that two Erysiphe species were associated with the outbreak of soybean powdery mildew. There was 16% difference between the two species in genetic divergence of the ITS sequence. One species with perfect stage has an ITS sequence identical to that of Erysiphe glycines on Amphicarpaea and is identified as Erysiphe glycines based on the ITS sequence and morphology of ascomata. The second species, without the perfect stage, is likely to be Erysiphe diffusa (= Microsphaera diffusa), known as the fungus causing soybean powdery mildew in the United States, because the ITS sequences are identical to those from materials collected in the United States. However, we need materials having ascomata of E. diffusa to confirm the species name. Received: March 15, 2002 / Accepted: May 22, 2002     

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.     

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 characterization of a leaf margin necrosis mutant in wild annual soybean (Glycine soja)

A leaf margin necrosis mutant was observed in a wild annual soybean (Glycine soja Sieb. & Zucc.) population from South Korea. Genetic studies showed that it was controlled by a single recessive nuclear gene, designatedlmn. TheLmn locus segregated independently of theAp, Dial, Dia3, Idh2, Pgi1 andTi isozyme loci.     

5S ribosomal gene variation in the soybean and its progenitor

Summary The soybean, Glycine max and its wild progenitor, Glycine soja, have been surveyed for repeat length variation for the nuclearly encoded 5S ribosomal RNA genes. There is little variation among the 33 accessions assayed, with a common repeat length of 345 bases being typical of both taxa. A 334 base size variant was encountered in individuals from two populations of G. soja from China. The low level of variability is in marked contrast to the variation observed within and between the species of the perennial subgenus Glycine.     

Genetic spatial clustering: significant implications for conservation of wild soybean (Glycine soja: Fabaceae)

Knowledge of spatial patterns of genetic variation within populations of wild relative species has significant implications with respect to sampling strategies for ex situ and in situ conservation. To study spatial genetic structure of wild soybean (Glycine soja Sieb. et Zucc.) at the fine scale, three natural populations in northern China were analyzed using inter-simple sequence repeat (ISSR) fingerprints for estimating kinship coefficients. A regression analysis of kinship coefficients against spatial distances revealed that individuals occurring close together tended to be more genetically related. The Sp statistic further indicated a comparable spatial pattern among the three wild soybean populations with similar Sp values (mean = 0.0734, varied from 0.0645 to 0.0943) detected across the three populations. Genetic patches were on average ca. 20 m in size, and the effective neighborhood sizes varied between 10 and 15 m. The spatial genetic structure evident in the wild soybean populations may be attributed to the restricted seed dispersal and predominant inbreeding mating system of this species. The detection of family structure in the populations of wild soybean has a significant implication for the effective conservation of the important genetic resources.     

The genetic diversity and population structure of wild soybean evaluated by chloroplast and nuclear gene sequences

Glycine soja, also called wild soybean, is the wild ancestor of domesticated soybean (Glycine max), and one of the world's major cultivated crops. Wild soybean is a valuable resource for the breeding of cultivated soybean and harbors useful genes or agronomic traits. To use and conserve this valuable resource, we conducted a study to evaluate the genetic diversity and population structure of wild soybean using the sequencing data of two nuclear loci (AF105221 and PhyB) and one chloroplast locus (trnQ-rps16) of more than 600 individuals representing 53 populations throughout the natural distribution range. The results showed that most of the variation was found within the populations and groups, but significant genetic differentiation was also detected among different eco-geographical groups. Correlations between genetic and geographical distance at all the loci were consistent with the isolation by distance gene flow model. G. soja exhibited the highest genetic diversity in middle and downstream of Yangzi River (MDYR) region, followed by North East China (NEC), and was the lowest in North West China (NWC). We concluded that both in situ and ex situ conservation strategies required for wild soybean populations, especially which are native to MDYR and NEC regions.     

Postglacial range shift and demographic expansion of the marine intertidal snail Batillaria attramentaria

To address the impacts of past climate changes, particularly since the last glacial period, on the history of the distribution and demography of marine species, we investigated the evolutionary and demographic responses of the intertidal batillariid gastropod, Batillaria attramentaria, to these changes, using the snail as a model species in the northwest Pacific. We applied phylogeographic and divergence population genetic approaches to mitochondrial COI sequences from B. attramentaria. To cover much of its distributional range, 197 individuals collected throughout Korea and 507 publically available sequences (mostly from Japan) were used. Finally, a Bayesian skyline plot (BSP) method was applied to reconstruct the demographic history of this species. We found four differentiated geographic groups around Korea, confirming the presence of two distinct, geographically subdivided haplogroups on the Japanese coastlines along the bifurcated routes of the warm Tsushima and Kuroshio Currents. These two haplogroups were estimated to have begun to split approximately 400,000 years ago. Population divergence analysis supported the hypothesis that the Yellow Sea was populated by a northward range expansion of a small fraction of founders that split from a southern ancestral population since the last glacial maximum (LGM: 26,000–19,000 years ago), when the southern area became re‐submerged. BSP analyses on six geographically and genetically defined groups in Korea and Japan consistently demonstrated that each group has exponentially increased approximately since the LGM. This study resolved the phylogeography of B. attramentaria as a series of events connected over space and time; while paleoceanographic conditions determining the connectivity of neighboring seas in East Asia are responsible for the vicariance of this species, the postglacial sea‐level rise and warming temperatures have played a crucial role in rapid range shifts and broad demographic expansions of its populations.     

野生大豆耐盐性研究进展

野生大豆对于拓宽大豆种质遗传基础和丰富大豆种质基因库具有重要意义.该文从野生大豆的资源概况及优良性状、耐盐机理和利用野生大豆提高栽培大豆耐盐性等三个方面,对近年来国内外有关野生大豆耐盐性的解剖结构、生理基础、分子生物学基础等方面的研究进展进行了系统综述,并提出野生大豆通过茎叶表皮上的"腺体"及对Na+和Cl-的排斥性,实现对盐渍环境的颉颃作用.较强的抗氧化能力、大豆异黄酮代谢和耐盐基因也是其适应盐渍环境的重要原因.今后应对野生大豆耐盐机理的遗传学基础进行深入研究,并通过种群保护以保障野生大豆的发掘鉴定和创新利用.     

Population structure of the wild soybean (Glycine soja) in China: implications from microsatellite analyses

Background and Aims Wild soybean (Glycine soja), a native species of East Asia, is the closest wild relative of the cultivated soybean (G. max) and supplies valuable genetic resources for cultivar breeding. Analyses of the genetic variation and population structure of wild soybean are fundamental for effective conservation studies and utilization of this valuable genetic resource. Methods In this study, 40 wild soybean populations from China were genotyped with 20 microsatellites to investigate the natural population structure and genetic diversity. These results were integrated with previous microsatellite analyses for 231 representative individuals from East Asia to investigate the genetic relationships of wild soybeans from China. Key Results Analysis of molecular variance (AMOVA) revealed that 43·92 % of the molecular variance occurred within populations, although relatively low genetic diversity was detected for natural wild soybean populations. Most of the populations exhibited significant effects of a genetic bottleneck. Principal co-ordinate analysis, construction of a Neighbor-Joining tree and Bayesian clustering indicated two main genotypic clusters of wild soybean from China. The wild soybean populations, which are distributed in north-east and south China, separated by the Huang-Huai Valley, displayed similar genotypes, whereas those populations from the Huang-Huai Valley were different. Conclusions The previously unknown population structure of the natural populations of wild soybean distributed throughout China was determined. Two evolutionarily significant units were defined and further analysed by combining genetic diversity and structure analyses from Chinese populations with representative samples from Eastern Asia. The study suggests that during the glacial period there may have been an expansion route between south-east and north-east China, via the temperate forests in the East China Sea Land Bridge, which resulted in similar genotypes of wild soybean populations from these regions. Genetic diversity and bottleneck analysis supports that both extensive collection of germplasm resources and habitat management strategies should be undertaken for effective conservation studies of these important wild soybean resources.     

A new allele of flower color gene <Emphasis Type="Italic">W1</Emphasis> encoding flavonoid 3'5'-hydroxylase is responsible for light purple flowers in wild soybean <Emphasis Type="Italic">Glycine soja</Emphasis>

Background  

Glycine soja is a wild relative of soybean that has purple flowers. No flower color variant of Glycine soja has been found in the natural habitat.     

Genetic differentiation in relation to seed weights in wild soybean species (Glycine soja Sieb. & Zucc.)

Seed weight is one of the most important botanical and phylogenetic characteristics. The study objective was to understand whether there is genetic difference in different seed weights of wild soybean (Glycine soja Sieb. & Zucc.). A total of 563 wild soybean samples, which belonged separately to genebank germplasm accessions (220 samples), one regional population samples (293 plants) and one natural population (150 plants), were analyzed using microsatellite markers. Of four size classes, the smallest seed size type had the highest coefficient of variation in seed weight; small and large seed types had relatively great genetic differences. In the national genebank germplasm accessions, genetic diversity gradually decreased from quantitatively dominant small and middling seed types to less frequent large seed types. In the regional and natural populations, generally, small to middling seed sizes had higher genetic diversity than the smallest and larger seed sizes. Cluster analysis revealed genetic differences in seed size traits. The semi-wild type (Glycine gracilis Skvortzow) was the most genetically differentiated from other seed sizes. However, it was also clearly shown that the phylogenic genetic differentiation among seed sizes was less than the genetic differentiation among geographical habitat populations in the wild soybean species.     

QTL associated with horizontal resistance to soybean cyst nematode in Glycine soja PI464925B

Soybean cyst nematode (Heterodera glycines Ichinohe; SCN) is the primary disease responsible for yield loss of soybean [Glycine max (L.) Merr.]. Resistant cultivars are an effective management tool; however, the sources currently available have common resistant genes. Glycine soja Sieb. and Zucc., the wild ancestor of domesticated soybean, represents a diverse germplasm pool with known SCN resistance. The objectives of this research were to: (1) determine the genetic variation and inheritance of SCN resistance in a G. max (‘S08-80’) × G. soja (PI464925B) F _4:5 recombinant inbred line (RIL) population; and (2) identify and evaluate quantitative trait loci (QTL) associated with SCN resistance. Transgressive segregation for resistance was observed, although neither parent was resistant to the Chatham and Ruthven SCN isolates. Broad sense heritability was 0.81 for the Ruthven and 0.91 for the Chatham isolate. Root dry weight was a significant covariate that influenced cyst counts. One RIL [female index (FI) = 5.2 ± 1.11] was identified as resistant to the Chatham isolate (FI < 10). Seventeen and three RILs infected with Chatham and Ruthven isolates, respectively, had mean adjusted cyst counts of zero. Unique and novel QTL, which derived resistance from G. soja, were identified on linkage groups I, K, and O, and individually explained 8, 7 and 5% (LOD = 2.1–2.7) of the total phenotypic variation, respectively. Significant epistatic interactions were found between pairs of SSR markers that individually may or may not have been associated with SCN resistance, which explained between 10 and 15% of the total phenotypic variation. Best-fit regression models explained 21 and 31% of the total phenotypic variation in the RIL population to the Chatham and Ruthven isolates, respectively. The results of this study help to improve the understanding of the genetic control of SCN resistance in soybean caused by minor genes resulting in horizontal resistance. The incorporation of the novel resistance QTL from G. soja could increase the durability of SCN-resistance in soybean cultivars, especially if major gene resistance breaks down.