BACKGROUND
the cultivated soybean  is an economically important crop that grown all over the world.
with an average of ~38% protein and ~18% oil content in seeds, soybean provides 69% of dietary protein and 30% of vegetable oil consumption worldwide .
modern soybean cultivars were originally domesticated from its wild progenitor  more than 3000 years ago, which was an endemic species in china  <cit> .
since then, a variety of morphological and physiological changes except for reproductive isolation have occurred that distinguish soybean cultivars from their wild ancestor.
wild soybeans possess much higher adaptability to various natural environments such as drought and salt stress, whereas cultivated soybeans exhibit a bush-type growth habit with large seeds, variable seed coat colors and a stout primary stem.
wild soybeans also differ in the extent of photosynthesis capacity, pod dehiscence and number from cultivated soybeans  <cit> .
heritable changes occurred during plant domestication are being revealed by gene mapping and genomic analyses  <cit> .
the availability of soybean genome and high throughput sequencing technologies provides excellent opportunity to excavate the domestication events and phenotypic diversification at the genome level  <cit> .
re-sequenced soybeans representing wild and cultivated accessions revealed the nature and extent of genetic diversity in both populations  <cit> .
another research reported a reservoir of genes that were affected by early domestication and modern genetic improvement  <cit> .
besides, several domestication-related traits have been studied and proposed to be controlled by a small number of genes or several major qtls  <cit> .
however, more analyses are needed to delimit the regions of these qtls and the footprints of domestication for further gene mapping.
from an evolutionary perspective, if a mutation happens to be beneficial to the species, it will spread to the population immediately by selection  <cit> .
during crop domestication, strong selective pressure caused traits of interests to be fixed in a founder population in quite a short time  <cit> .
probably, advantageous mutations underlying traits of interests will be subject to fixation in the population.
these fixation events differ from those in natural populations, because artificial selection usually acted on alleles that were likely neutral or nearly neutral before domestication.
thus, understanding nucleotide fixation driven by artificial selection is indispensable to complete the picture of soybean evolution.
in this research, the published soybean sequencing data were collected to identify single nucleotide variations , based on which we detected the genomic regions affected by artificial selection during domestication and improvement.
in these footprints, nucleotide fixations that happened in all cultivars were potentially caused by artificial selection, and the genes with these nucleotides were further analyzed, and some of these genes were associated with agronomic traits through functional annotation and qtl meta-analysis.
this kind of investigation will provide clues to understand the differentiation of wild and cultivated soybeans.
besides, fundamental practical information will be obtained for future enhancement of cultivars through traditional breeding and transgenic methods.
RESULTS
estimation of single nucleotide variations among soybean populations
recently, a set of diverse soybean individuals was sequenced and reported based on ngs platforms  <cit> .
these soybeans, representing wild and cultivars that mainly consist of landrace and modern elite accessions in east asia, were selected based on intensive molecular and phenotypic analysis to maximally reflect the genetic diversity of soybeans .
it provides us an important resource to depict the genetic diversity of wild and cultivated populations, and to detect the footprints of domestication events.
thus, we downloaded all the short reads of sequencing soybeans from ncbi short read archive under accession numbers sra <dig>  sra <dig>  srp <dig>  and erp <dig>  these reads were aligned to the soybean reference genome glycine max  with soap <dig>  <cit> , and were subsequently used to detect snvs with soapsnp pipeline  <cit> .
a total of  <dig> , <dig> snvs were identified across all accessions, of which  <dig> , <dig> and  <dig> , <dig> appear in wild and cultivars, respectively.
previous reports with the same pipeline have shown that the snv calling accuracy is 95-99%, with false-positive and false-negative rates to be ~2% and ~3%, respectively  <cit> .
to estimate the coverage of these snvs in the whole soybean germplasm, we employed a random sampling approach to investigate the accumulation of snvs detected in different accessions .
the end of the snv curve tends to be flat, which indicates that the snvs identified here probably reach saturation in soybean germplasm.
it is sufficient for as few as  <dig> accessions to detect 95% of all snvs in different populations.
for cultivated soybeans, only  <dig> individuals can achieve 95% of snvs.
approximately  <dig>  million snvs would reach saturation in cultivars, which are far less than those in wild soybeans.
in previous work  <cit> , lam et al reported  <dig>  million snvs in  <dig> soybeans, while we discovered  <dig> , <dig> more in the same individuals by a larger population.
a large number of rare snvs and those with low allele frequency were missed in former analysis due to strict filtering conditions and a small number of individuals .
although some very rare snvs still remain to be discovered, we have identified a substantial majority of the common snvs in soybeans.figure  <dig>
detection of single nucleotide variations in sequencing soybean accessions.
accumulated snv coverage in cultivars, wild and all accessions;  distribution of missing snvs in previous report by lam et al.
soybean has suffered several genetic bottlenecks, such as early domestication producing lots of asian landrace, the introduction of few landraces to north america, and modern extensive breeding activities  <cit> .
subsequently, different level of genetic diversity was reduced during these human-mediated events.
more snvs were identified in wild than in cultivated accessions.
two common statistics used to measure nucleotide diversity are the pairwise divergence per nucleotide θπ  <cit>  and watterson estimator θw  <cit>  that corrected for sample size.
whole-genome analysis using these parameters shows a higher level of genetic diversity in wild populations .
estimated by θπ, the average diversity within wild, landrace and elite cultivars are  <dig>  × 10- <dig>   <dig>  × 10- <dig>  and  <dig>  × 10- <dig> per nucleotide, respectively.
considering the cultivars consist of landrace and elites, the average θπ is  <dig> × 10- <dig> in cultivated population.
it is notable that the cultivars have retained only  <dig> % of the sequence diversity present in wild soybeans, which is lower than previous estimation  <cit> .
the genetic diversity was reduced by  <dig> % in early domestication and further reduced by  <dig> % in genetic improvement.figure  <dig>
analysis of genetic diversity and phylogenetic relationship among soybean accessions.
reduction of genetic diversity from wild, to landrace and then to elite soybeans;  a neighbor-joining tree;  principal component analysis of soybeans.
the reduction of genetic diversity eroded by artificial selection could also be reflected by phylogenetic tree  and principle component analysis .
the wild soybeans shattered in a loose 3-dimension space, while cultivated soybeans formed a relatively tight cluster distinct from the wild individuals.
within the cluster, however, the landraces were not clearly separated from elite cultivars.
some landraces mixed with wild group in our analysis, indicating the early domestication process probably accompanied with considerable gene flow with the wild ancestors.
in addition to artificial selection, the genetic erosion can also reflect the narrow genetic base of cultivated soybeans  <cit> .
analysis of representative wild and cultivated soybeans provides us a comprehensive insight into such evolutionary events that affected population dynamics of soybeans.
detecting artificial selection and nucleotide fixation in soybeans
the signal of artificial selection could be detected by the loss of genetic diversity, which shaped selective sweeps around beneficial alleles on the genomes  <cit> .
to further elucidate the effects of domestication, we detected the genomic regions showing artificial selection signals by genetic bottleneck model  <cit>  and population branch statistics  <cit> .
the sequenced accessions except c <dig> and c <dig> were grouped into wild and cultivated population to detect selection signals in early domestication process.
using a sliding window approach, we calculated the distribution of θπ and tajima’s d  <cit>  in wild and cultivated populations along the genome.
regions with significantly lower θπ  and lower tajima’s d  in cultivars than that in wild accessions were treated as putative candidates that were affected by early domestication .
however, signals of very recent natural selection could be easily omitted using the above bottleneck model.
to detect signatures that shaped in modern crop improvement, we employed an effective method known as population branch statistics.
taking wild soybeans as control, we recalculated the divergence index fst  <cit>  in a sliding window along the genome, based on which we detected significant signals  to infer selective footprints from landraces to elite cultivars .
this approach had been shown to be effective in identifying recent artificial selection considering the very short time of modern breeding practice  <cit> .
a total of  <dig> regions comprising  <dig>  mb genome sequences and  <dig> regions with a length of  <dig>  mb were affected by early domestication and genetic improvement, respectively.
based on the latest annotation,  <dig>  genes with  <dig>  transcripts were involved in early domestication, whereas  <dig>  genes with  <dig>  transcripts were affected in subsequent improvement.figure  <dig>
footprints of artificial selection during  early domestication and  modern improvement.
during the human-mediated breeding process, the strongly selected advantageous mutations could become fixed as these mutations increase in frequency in a population  <cit> .
a selective sweep is shaped when a selected mutation goes to fixation, because it reduces variability in the neighboring region where neutral variants are segregating  <cit> .
a nucleotide fixation locus was defined when a snv has a unique genotype in one population while it exhibits polymorphic genotypes in the others.
to better understand how genes were affected by domestication events, we primarily focused on those with nucleotide fixation in the selective footprints.
we calculated the likelihood of genotypes of each individual and then we allocated the allele type with the maximum likelihood back to each individual as the consensus genotype.
after calibration,  <dig>  nucleotide fixations were identified in the selective regions in cultivars, which could be potentially caused by artificial selection.
compared with the genome-wide distribution, nucleotide fixations happened more frequently in the candidate regions of artificial selection .
nucleotide fixation accumulated substantially in cultivars and happened unevenly along chromosomes , indicating that some chromosomes were more susceptible to be affected by artificial selection.
nucleotide fixation also explains the reduction of genetic diversity in cultivated crops compared with their wild ancestors.
we analyzed the allele frequency of snvs in wild soybeans that were fixed in cultivars, as it represents the initial status of these nucleotide fixations before domestication.
the frequency spectrum shows that these snvs were almost neutral at the beginning of domestication .
since non-synonymous substitutions may result in a change in functions, they are subject to natural or artificial selection  <cit> .
of the nucleotide fixation happened in early domestication,  <dig>  located in coding sequences and  <dig>  of them caused non-synonymous substitutions in  <dig>  genes, which altered the amino acid sequences of the proteins.
for those loci fixed in modern improvement,  <dig>  located in coding sequences with  <dig> non-synonymous in  <dig> genes.
apparently, more nucleotide fixations were introduced to cultivars during domestication than those during improvement.figure  <dig>
the distribution of nucleotide fixation over the genome versus in the selective regions.
the window size was set to be 20 kb.
a central question in analyzing the genetic variations in a given population is to explore whether the population has different substructures  <cit> .
when analyzing the nucleotide fixations by pca and phylogenetic tree, two distinct clusters shaped between the cultivars and wild soybeans .
some noise always exists in inferring phylogenetic relationships among individuals, especially when they are subject to introgressive hybridization  <cit> .
cultivars tightly joined together without noise, supporting the hypothesis of a single rather than multiple evolutionary origins in soybean domestication  <cit> .
nucleotide fixation in wild soybeans
in the process of nucleotide substitution, the fixation of a mutation could spread through the population by random genetic drift or extreme natural selection  <cit> .
in the regions affected by artificial selection,  <dig>  nucleotide fixations happened in wild soybeans, which located in  <dig> transcripts corresponding to  <dig> genes.
nucleotide fixation happened more frequently in cultivars compared with wild soybeans.
to some degree, artificial selection could have promoted the occurrence of fixation events.
however, genetic bottlenecks caused by domestication often results in a smaller effective population size of cultivars than that of wild soybeans, which would also contribute to an elevated level of nucleotide fixation.
genes affected by nucleotide fixations were involved in kinds of biological activities as described in the kyoto encyclopedia of genes and genomes  database .
the ability of resistance to pathogen in wild soybeans is much broader than that in cultivated soybeans  <cit> .
interestingly, glyma20g <dig> gene is an ortholog of the disease resistance gene rpm <dig>  which was detected and characterized using molecular genetic approach in arabidopsis  <cit> .
in soybeans, the rpm <dig> gene was recently reported being under purifying selection  <cit> .
it serves as an example that natural selection in the wild population also caused nucleotide fixations, although its strength was less than artificial selection.
agronomic traits affected by selective nucleotide fixation
during domestication, artificial selection is thought to have extremely strong selective pressure on ancestral population for desired phenotypes  <cit> .
the strong selection exerted by human led to an excessive amount of nucleotide fixations during domestication.
artificial selection during soybean domestication has modified a number of traits including seed size, seed color, plant height and prostrate habitat, shaping the domestication syndrome  <cit> .
to analyze the effects of nucleotide fixation during artificial selection, we focused on genes within qtls responsible for domestication-related traits , such as oil content, pod number, lodging, plant height, etc.
meta-analysis of these qtls revealed that  <dig> of them responsible for  <dig> traits and  <dig> for  <dig> traits were affected by nucleotide fixation in early domestication and modern improvement, respectively .
total qtl regions were narrowed down from  <dig>  mb to  <dig>  mb assisted by selective signals.
analysis of related genes, as well as their orthologs through comparative genomics, could provide information on their potential functions under artificial selection.
as an agriculturally important trait, grain filling makes a significant contribution to grain weight  <cit> .
the gene grain incomplete filling  <dig>  was detected to be responsible and associated with this domestication syndrome  <cit> .
it was reported to encode a cell-wall invertase required for carbon partitioning during early grain filling in rice.
a selective gene glyma03g <dig> with nucleotide fixation in domestication is an ortholog of gif <dig> and this gene was involved in the carbohydrate metabolism pathway by searching kegg .
besides, this gene was covered by the qtls responsible for lodging and pod number.
it indicates that glyma03g <dig> is a potential candidate gene, which could be used for further soybean breeding.
flower and pod numbers per plant are important agronomic traits for grain yield in soybean.
to detect the genes involved in flower and pod numbers will help to understand the genetic basis of soybean yield  <cit> .
two genes, glyma07g <dig> and glyma07g <dig>  with nucleotide fixation introduced in improvement, are orthologs of comt <dig> gene encoding caffeic acid 3-o-methyltransferase .
it differentially expressed in hair cells of growing pod, the possible location of vanillin biosynthesis  <cit> .
another five selective genes with nucleotide fixation mediated by domestication and improvement encode a kind of protein responsible for the transportation of inositol.
these genes were covered by qtls responsible for seed-coat color, protein and pod number.
previous study showed that the total number of mature pods considerably higher due to the application of inositol, indicating the positive effect in pod number  <cit> .
it suggested that deficiency of lignin biosynthesis resulted in growth reduction and dwarfing  <cit> .
the gene glyma13g <dig> is linked to marker sat <dig> that associate with seed weight.
as an orthologs of nifu gene, it is required for full activation of nitrogenase catalytic components  <cit> .
nifu protein has been suggested to either mobilize the fe necessary for nitrogenase fe-s cluster formation or provide an intermediate fe-s cluster assembly site  <cit> .
in addition, the gene was reported to be related to seed weight  <cit> .
as nitrogen fixation is imperative in soybean growth, glyma13g <dig> gene could also be a putative candidate gene responsible for seed weight through activating biological nitrogen fixation.
the flowering of soybean represents the transition from a vegetative state to a reproductive state, making a contribution to the yield.
meta-analysis of qtls identified  <dig> selective genes with non-synonymous nucleotide fixation responsible for flower number in soybean.
carbon fixation in the process of photosynthesis is pivotal to soybean production.
seven selective genes with nucleotide fixation were involved in photosynthesis or photosystem.
besides, two selective genes glyma03g <dig> and glyma19g <dig> with nucleotide fixation were identified as orthologs of luminidependens, which is involved in the timing of flowering in arabidopsis  <cit> .
interestingly,  <dig> and  <dig> selective genes with nucleotide fixation in domestication and improvement, respectively, were annotated to be, or related with transcription factors.
analysis of all the genes subject to artificial selection with agrigo  <cit>  also told an accumulation of transcription factors by fisher’s exact test and the permutation test .
most of the genes cloned to date that responsible for domestication related traits in crops were proved to be transcription factors, such as teosinte branched  <dig>   <cit> , shattering   <cit> , six-rowed spike   <cit> , etc.
it is probably because the human mediated domestication history was momentary compared with the long natural evolution; changing the transcription factors probably the easiest way happened to affect the agricultural traits of interest.
however, putative candidate genes underlying these domestication-targeted phenotypes have diverse functions, which need to be validated by further experiments.
plant-pathogen interaction affected by artificial selection
domestication caused complex morphological and physiological changes in soybeans.
annotated by the kegg and agrigo database, selective genes were associated with different biological functions, among which, plant-pathogen interaction, sequence-specific dna binding, phenylpropanoid biosynthesis, starch and sucrose metabolism are over-represented categories .
plant-pathogen interactions are conducted between a pathogen and the host plant.
in nature, plants are generally resistant to most invading pathogens due to innate ability to recognize them through successful defenses.
when an exception happens, a pathogen would cause diseases in its host  <cit> .
however, pathogens could also cause diseases if they have evolved to evade detection or suppress host defense mechanisms, or both.
the effects of plant-pathogen interactions are of particular relevance during early domestication events on agricultural systems  <cit> .
thus, understanding the genetic basis of why a certain pathogen causes disease in its host plant instead of others has long intrigued and motivated plant pathologists.figure  <dig>
functional annotation of selective genes with nucleotide fixation introduced in early domestication and modern improvement.
a total of  <dig> selective genes with nucleotide fixation were involved in plant-pathogen interactions .
of them, two selective genes glyma14g <dig> and glyma08g <dig> with nucleotide fixation are orthologs of rps <dig> gene.
the disease resistance gene rps <dig> was isolated using positional cloning and further screen for susceptible mutant  <cit> .
the rps <dig> protein contains two characteristics of a large family of plant r genes: a nucleotide-binding site and a leucine-rich repeat region  <cit> .
it is consistent with previous report that rps <dig> locus exhibit selection signals by examining a worldwide sample of  <dig> arabidopsis accessions, and the n-terminal part of the leucine-rich repeat region was a probable target of selection  <cit> .
cyclic nucleotide-gated ion  channels are ion channels that function in the pathogen signaling cascade by facilitating ca2+ uptake into the cytosol  <cit> .
two selective genes with nucleotide fixation were detected to encode cng channels.
the topology of their proteins was predicted using tmhmm, which is based on a hidden markov model  <cit> .
the two genes encode transmembrane proteins with nucleotide fixation located outside the membrane .
besides, eight selective genes are orthologs of transmembrane receptor kinase fls <dig>  which acts as pathogen-associated molecular pattern signals in triggering the innate immune response  <cit> .
in addition, the category of terpene synthase activity was also enriched with six selective genes involved in .
terpenes are one of the most important defensive plant compounds against herbivores and pathogens  <cit> .
recently, a new monoterpene synthase gene gmnes was identified and characterized in soybean  <cit> .
its transcription was up-regulated in soybean when infested with cotton leafworm.
our analysis indicates the gene was possibly affected by artificial selection during soybean domestication.
discussion
nucleotide fixation was crucial in soybean divergence
domestication led to significant morphological divergence between cultivated and wild soybeans.
wild soybean exhibits, for example, twining and vine stem, severer shattering, impermeable seed coats, pod cracking sensitivity, small seeds, and low oleic acid, all of which are seldom observed in cultivars  <cit> .
deciphering how cultivated soybean have been transformed from its wild ancestor is advantageous both from genetic and evolutionary perspectives.
with the available sequencing data, we comprehensively estimated the saturation number of snvs in soybean germplasm and detected a set of candidate genes showing artificial selection signals.
to some degree, analysis of artificial selection and nucleotide fixation unravels the mystery of soybean domestication and subsequent improvement.
based on nucleotide fixation, our analysis supports a single evolutionary origin of domesticated soybean.
during domestication, only lines with certain agriculturally important traits were selected, resulting in a genome-wide reduction of genetic diversity or so-called selective sweep in cultivated crops  <cit> .
one possible explanation for the reduction is that an excess of nucleotide fixation happened in cultivars compared to wild soybeans.
meta-analysis of qtls responsible for domestication related traits and the selective genes provided insights into the role of nucleotide fixation played in morphological differentiation between wild and cultivated soybeans.
using comparative genomics, an amount of genes was found to be orthologs of those whose function was validated and responsible for corresponding traits in other plants.
nucleotide fixation happened in those genes responsible for agronomically important traits.
although traditional linkage and association mapping were used to dissect these traits, they failed to detect genetic changes caused by domestication and improvement  <cit> .
our analysis here provides valuable information for further qtl mapping and will facilitate molecular assisted selection in soybean breeding practice.
artificial selection accelerates nucleotide fixation
domestication was an evolutionary process where the characters of interests were selected, such as loss of seed dispersal, higher yield and increasing abiotic resistance.
the detection of selective loci during crop domestication contributes to modern breeding efforts and the opportunity to improve genomic selection models  <cit> .
recently, genome-wide scans based on genetic bottlenecks have been successfully applied to detect footprints of selection in plants by surveying both natural and cultivated species  <cit> .
artificial selection of a beneficial mutation will lead to an elevated frequency in a population.
eventually, allele frequencies will be skewed and nucleotide fixation happened after plant domestication.
our analysis focused on to what degree nucleotide fixation was caused by artificial selection during soybean domestication.
more nucleotide fixation happened in cultivars than those in wild soybeans, indicating that artificial selection was much stronger than natural selection.
however, the effective population size of cultivated soybeans was substantially reduced during domestication  <cit> , which could make a nucleotide seem to be fixed in cultivars.
that mainly explains why nucleotide fixations were observed in cultivars across the soybean genome.
considering nucleotide fixation accumulated in footprints of domestication and improvement, artificial selection probably accelerated the occurrence of fixation in soybean breeding activities.
even thought, some of them could be also caused by the shrinking population size, especially when different haplotypes shaped in those selective sweeps.
these fixations are extremely hard to be distinguished in current samples.
morphological transition can be achieved by a mutation at a single locus  <cit> , and artificial selection can rapidly change domestication targeted phenotypes within  <dig> generations  <cit> .
domestication could be a rapid instead of a slow or gradual process, given strong selective pressures and a suitable genetic architecture.
this was supported by the severe reduction of genetic diversity and large selective sweeps.
in the process of domestication, any mutations detrimental to the traits of interests were eliminated immediately, whereas those advantageous ones were strongly selected, diffused and eventually fixed in a population.
the environments wild soybeans grow in are various and usually harsh, resulting in diversifying selection instead of strong directional selection.
what’s more, selection intensity imposed by natural selection was disparate in diverse habitats.
these reasons also explain why artificial selection was much stronger than natural selection in crop domestication.
evolutionary perspective of nucleotide fixation
a long-term goal of crop genomics is to determine to what extent artificial selection impacts genomic variation patterns within and between populations.
there are both genetic and statistical approaches to detect signals of hitchhiking caused selective sweeps  <cit> .
the hitchhiking effect is contingent on the nature of genetic variations and how selection acts on them.
generally, there are at least three evolutionary routes by which a novel mutation may fix: drift to fixation for nearly neutral mutation; rapidly sweep to fixation, so-called hard sweep for beneficial mutation; and soft sweep to fixation for those initially neutral but later become beneficial for some reason.
affected by artificial selection, a pre-exist mutation which became beneficial during domestication rapidly increased in frequency toward nucleotide fixation, as what we found in our analysis.
when traits of interests during domestication were determined by multiple adaptive mutations at the same locus, artificial selection usually generates soft rather hard selective sweeps.
many studies focus on hard sweeps in which only a single adaptive haplotype was skewed to fixation in the population  <cit> , whereas multiple adaptive haplotypes formed simultaneously in a soft sweep.
lots of nucleotide fixations happened within quantitative traits, indicating the corresponding traits of interests were incrementally changed at various causal loci.
as a consequence, these sweeps related with artificial selection are likely to be both soft and incomplete.
in soybean, some traits related to yield were selected, such as seed weight, seed blooming and prostrate habit, for which these are usually major qtls responsible.
nevertheless, during intensive breeding human pursuits quality related traits such as protein content and lipid content, for which there are lots of small effect qtls responsible.
analysis of nucleotide fixation indicates that more soft selective sweeps happened in extensive breeding than in early domestication in soybean, which still needs further investigation.
CONCLUSIONS
we integrated the available sequencing accessions to describe a whole picture of soybean genetic diversity, artificial selection and concomitant nucleotide fixation.
there are approximately  <dig>  million snvs in soybean germplasm, of which about  <dig>  million reserved in cultivars.
the genetic diversity was reduced by  <dig> % in early domestication and subsequently reduced by  <dig> % in genetic improvement.
a total of  <dig>  and  <dig>  genes were involved in early domestication and subsequent improvement, respectively.
both processes introduced about  <dig>  million nucleotide fixations, which contributed to the divergence of wild and cultivated soybeans.
artificial selection probably accelerated the occurrence of nucleotide fixation, which affected some agronomic traits, as well as related biological pathways such as plant-pathogen interaction.

 the cultivated soybean is an economically important crop that grown all over the world. with an average of ~38% protein and ~18% oil content in seeds, soybean provides 69% of dietary protein and 30% of vegetable oil consumption worldwide. modern soybean cultivars were originally domesticated from its wild progenitor more than 3000 years ago, which was an endemic species in china. since then, a variety of morphological and physiological changes except for reproductive isolation have occurred that distinguish soybean cultivars from their wild ancestor. wild soybeans possess much higher adaptability to various natural environments such as drought and salt stress, whereas cultivated soybeans exhibit a bush-type growth habit with large seeds, variable seed coat colors and a stout primary stem. wild soybeans also differ in the extent of photosynthesis capacity, pod dehiscence and number from cultivated soybeans. heritable changes occurred during plant domestication are being revealed by gene mapping and genomic analyses. the availability of soybean genome and high throughput sequencing technologies provides excellent opportunity to excavate the domestication events and phenotypic diversification at the genome level. re-sequenced soybeans representing wild and cultivated accessions revealed the nature and extent of genetic diversity in both populations. another research reported a reservoir of genes that were affected by early domestication and modern genetic improvement. besides, several domestication-related traits have been studied and proposed to be controlled by a small number of genes or several major qtls. however, more analyses are needed to delimit the regions of these qtls and the footprints of domestication for further gene mapping. from an evolutionary perspective, if a mutation happens to be beneficial to the species, it will spread to the population immediately by selection. during crop domestication, strong selective pressure caused traits of interests to be fixed in a founder population in quite a short time. probably, advantageous mutations underlying traits of interests will be subject to fixation in the population. these fixation events differ from those in natural populations, because artificial selection usually acted on alleles that were likely neutral or nearly neutral before domestication. thus, understanding nucleotide fixation driven by artificial selection is indispensable to complete the picture of soybean evolution. in this research, the published soybean sequencing data were collected to identify single nucleotide variations, based on which we detected the genomic regions affected by artificial selection during domestication and improvement. in these footprints, nucleotide fixations that happened in all cultivars were potentially caused by artificial selection, and the genes with these nucleotides were further analyzed, and some of these genes were associated with agronomic traits through functional annotation and qtl meta-analysis. this kind of investigation will provide clues to understand the differentiation of wild and cultivated soybeans. besides, fundamental practical information will be obtained for future enhancement of cultivars through traditional breeding and transgenic methods.  estimation of single nucleotide variations among soybean populations recently, a set of diverse soybean individuals was sequenced and reported based on ngs platforms. these soybeans, representing wild and cultivars that mainly consist of landrace and modern elite accessions in east asia, were selected based on intensive molecular and phenotypic analysis to maximally reflect the genetic diversity of soybeans. it provides us an important resource to depict the genetic diversity of wild and cultivated populations, and to detect the footprints of domestication events. thus, we downloaded all the short reads of sequencing soybeans from ncbi short read archive under accession numbers sra sra srp and erp these reads were aligned to the soybean reference genome glycine max with soap, and were subsequently used to detect snvs with soapsnp pipeline. a total of, snvs were identified across all accessions, of which, and, appear in wild and cultivars, respectively. previous reports with the same pipeline have shown that the snv calling accuracy is 95-99%, with false-positive and false-negative rates to be ~2% and ~3%, respectively. to estimate the coverage of these snvs in the whole soybean germplasm, we employed a random sampling approach to investigate the accumulation of snvs detected in different accessions. the end of the snv curve tends to be flat, which indicates that the snvs identified here probably reach saturation in soybean germplasm. it is sufficient for as few as accessions to detect 95% of all snvs in different populations. for cultivated soybeans, only individuals can achieve 95% of snvs. approximately million snvs would reach saturation in cultivars, which are far less than those in wild soybeans. in previous work, lam et al reported million snvs in soybeans, while we discovered, more in the same individuals by a larger population. a large number of rare snvs and those with low allele frequency were missed in former analysis due to strict filtering conditions and a small number of individuals. although some very rare snvs still remain to be discovered, we have identified a substantial majority of the common snvs in soybeans.figure detection of single nucleotide variations in sequencing soybean accessions. accumulated snv coverage in cultivars, wild and all accessions; distribution of missing snvs in previous report by lam et al. soybean has suffered several genetic bottlenecks, such as early domestication producing lots of asian landrace, the introduction of few landraces to north america, and modern extensive breeding activities. subsequently, different level of genetic diversity was reduced during these human-mediated events. more snvs were identified in wild than in cultivated accessions. two common statistics used to measure nucleotide diversity are the pairwise divergence per nucleotide θπ and watterson estimator θw that corrected for sample size. whole-genome analysis using these parameters shows a higher level of genetic diversity in wild populations. estimated by θπ, the average diversity within wild, landrace and elite cultivars are  × 10-  × 10- and  × 10- per nucleotide, respectively. considering the cultivars consist of landrace and elites, the average θπ is × 10- in cultivated population. it is notable that the cultivars have retained only % of the sequence diversity present in wild soybeans, which is lower than previous estimation. the genetic diversity was reduced by % in early domestication and further reduced by % in genetic improvement.figure analysis of genetic diversity and phylogenetic relationship among soybean accessions. reduction of genetic diversity from wild, to landrace and then to elite soybeans; a neighbor-joining tree; principal component analysis of soybeans. the reduction of genetic diversity eroded by artificial selection could also be reflected by phylogenetic tree and principle component analysis. the wild soybeans shattered in a loose 3-dimension space, while cultivated soybeans formed a relatively tight cluster distinct from the wild individuals. within the cluster, however, the landraces were not clearly separated from elite cultivars. some landraces mixed with wild group in our analysis, indicating the early domestication process probably accompanied with considerable gene flow with the wild ancestors. in addition to artificial selection, the genetic erosion can also reflect the narrow genetic base of cultivated soybeans. analysis of representative wild and cultivated soybeans provides us a comprehensive insight into such evolutionary events that affected population dynamics of soybeans. detecting artificial selection and nucleotide fixation in soybeans the signal of artificial selection could be detected by the loss of genetic diversity, which shaped selective sweeps around beneficial alleles on the genomes. to further elucidate the effects of domestication, we detected the genomic regions showing artificial selection signals by genetic bottleneck model and population branch statistics. the sequenced accessions except c and c were grouped into wild and cultivated population to detect selection signals in early domestication process. using a sliding window approach, we calculated the distribution of θπ and tajima’s d in wild and cultivated populations along the genome. regions with significantly lower θπ and lower tajima’s d in cultivars than that in wild accessions were treated as putative candidates that were affected by early domestication. however, signals of very recent natural selection could be easily omitted using the above bottleneck model. to detect signatures that shaped in modern crop improvement, we employed an effective method known as population branch statistics. taking wild soybeans as control, we recalculated the divergence index fst in a sliding window along the genome, based on which we detected significant signals to infer selective footprints from landraces to elite cultivars. this approach had been shown to be effective in identifying recent artificial selection considering the very short time of modern breeding practice. a total of regions comprising  mb genome sequences and regions with a length of  mb were affected by early domestication and genetic improvement, respectively. based on the latest annotation, genes with transcripts were involved in early domestication, whereas genes with transcripts were affected in subsequent improvement.figure footprints of artificial selection during early domestication and modern improvement. during the human-mediated breeding process, the strongly selected advantageous mutations could become fixed as these mutations increase in frequency in a population. a selective sweep is shaped when a selected mutation goes to fixation, because it reduces variability in the neighboring region where neutral variants are segregating. a nucleotide fixation locus was defined when a snv has a unique genotype in one population while it exhibits polymorphic genotypes in the others. to better understand how genes were affected by domestication events, we primarily focused on those with nucleotide fixation in the selective footprints. we calculated the likelihood of genotypes of each individual and then we allocated the allele type with the maximum likelihood back to each individual as the consensus genotype. after calibration, nucleotide fixations were identified in the selective regions in cultivars, which could be potentially caused by artificial selection. compared with the genome-wide distribution, nucleotide fixations happened more frequently in the candidate regions of artificial selection. nucleotide fixation accumulated substantially in cultivars and happened unevenly along chromosomes, indicating that some chromosomes were more susceptible to be affected by artificial selection. nucleotide fixation also explains the reduction of genetic diversity in cultivated crops compared with their wild ancestors. we analyzed the allele frequency of snvs in wild soybeans that were fixed in cultivars, as it represents the initial status of these nucleotide fixations before domestication. the frequency spectrum shows that these snvs were almost neutral at the beginning of domestication. since non-synonymous substitutions may result in a change in functions, they are subject to natural or artificial selection. of the nucleotide fixation happened in early domestication, located in coding sequences and of them caused non-synonymous substitutions in genes, which altered the amino acid sequences of the proteins. for those loci fixed in modern improvement, located in coding sequences with non-synonymous in genes. apparently, more nucleotide fixations were introduced to cultivars during domestication than those during improvement.figure the distribution of nucleotide fixation over the genome versus in the selective regions. the window size was set to be 20 kb. a central question in analyzing the genetic variations in a given population is to explore whether the population has different substructures. when analyzing the nucleotide fixations by pca and phylogenetic tree, two distinct clusters shaped between the cultivars and wild soybeans. some noise always exists in inferring phylogenetic relationships among individuals, especially when they are subject to introgressive hybridization. cultivars tightly joined together without noise, supporting the hypothesis of a single rather than multiple evolutionary origins in soybean domestication. nucleotide fixation in wild soybeans in the process of nucleotide substitution, the fixation of a mutation could spread through the population by random genetic drift or extreme natural selection. in the regions affected by artificial selection, nucleotide fixations happened in wild soybeans, which located in transcripts corresponding to genes. nucleotide fixation happened more frequently in cultivars compared with wild soybeans. to some degree, artificial selection could have promoted the occurrence of fixation events. however, genetic bottlenecks caused by domestication often results in a smaller effective population size of cultivars than that of wild soybeans, which would also contribute to an elevated level of nucleotide fixation. genes affected by nucleotide fixations were involved in kinds of biological activities as described in the kyoto encyclopedia of genes and genomes database. the ability of resistance to pathogen in wild soybeans is much broader than that in cultivated soybeans. interestingly, glyma20g gene is an ortholog of the disease resistance gene rpm which was detected and characterized using molecular genetic approach in arabidopsis. in soybeans, the rpm gene was recently reported being under purifying selection. it serves as an example that natural selection in the wild population also caused nucleotide fixations, although its strength was less than artificial selection. agronomic traits affected by selective nucleotide fixation during domestication, artificial selection is thought to have extremely strong selective pressure on ancestral population for desired phenotypes. the strong selection exerted by human led to an excessive amount of nucleotide fixations during domestication. artificial selection during soybean domestication has modified a number of traits including seed size, seed color, plant height and prostrate habitat, shaping the domestication syndrome. to analyze the effects of nucleotide fixation during artificial selection, we focused on genes within qtls responsible for domestication-related traits, such as oil content, pod number, lodging, plant height, etc. meta-analysis of these qtls revealed that of them responsible for traits and for traits were affected by nucleotide fixation in early domestication and modern improvement, respectively. total qtl regions were narrowed down from  mb to  mb assisted by selective signals. analysis of related genes, as well as their orthologs through comparative genomics, could provide information on their potential functions under artificial selection. as an agriculturally important trait, grain filling makes a significant contribution to grain weight. the gene grain incomplete filling was detected to be responsible and associated with this domestication syndrome. it was reported to encode a cell-wall invertase required for carbon partitioning during early grain filling in rice. a selective gene glyma03g with nucleotide fixation in domestication is an ortholog of gif and this gene was involved in the carbohydrate metabolism pathway by searching kegg. besides, this gene was covered by the qtls responsible for lodging and pod number. it indicates that glyma03g is a potential candidate gene, which could be used for further soybean breeding. flower and pod numbers per plant are important agronomic traits for grain yield in soybean. to detect the genes involved in flower and pod numbers will help to understand the genetic basis of soybean yield. two genes, glyma07g and glyma07g with nucleotide fixation introduced in improvement, are orthologs of comt gene encoding caffeic acid 3-o-methyltransferase. it differentially expressed in hair cells of growing pod, the possible location of vanillin biosynthesis. another five selective genes with nucleotide fixation mediated by domestication and improvement encode a kind of protein responsible for the transportation of inositol. these genes were covered by qtls responsible for seed-coat color, protein and pod number. previous study showed that the total number of mature pods considerably higher due to the application of inositol, indicating the positive effect in pod number. it suggested that deficiency of lignin biosynthesis resulted in growth reduction and dwarfing. the gene glyma13g is linked to marker sat that associate with seed weight. as an orthologs of nifu gene, it is required for full activation of nitrogenase catalytic components. nifu protein has been suggested to either mobilize the fe necessary for nitrogenase fe-s cluster formation or provide an intermediate fe-s cluster assembly site. in addition, the gene was reported to be related to seed weight. as nitrogen fixation is imperative in soybean growth, glyma13g gene could also be a putative candidate gene responsible for seed weight through activating biological nitrogen fixation. the flowering of soybean represents the transition from a vegetative state to a reproductive state, making a contribution to the yield. meta-analysis of qtls identified selective genes with non-synonymous nucleotide fixation responsible for flower number in soybean. carbon fixation in the process of photosynthesis is pivotal to soybean production. seven selective genes with nucleotide fixation were involved in photosynthesis or photosystem. besides, two selective genes glyma03g and glyma19g with nucleotide fixation were identified as orthologs of luminidependens, which is involved in the timing of flowering in arabidopsis. interestingly, and selective genes with nucleotide fixation in domestication and improvement, respectively, were annotated to be, or related with transcription factors. analysis of all the genes subject to artificial selection with agrigo also told an accumulation of transcription factors by fisher’s exact test and the permutation test. most of the genes cloned to date that responsible for domestication related traits in crops were proved to be transcription factors, such as teosinte branched, shattering, six-rowed spike, etc. it is probably because the human mediated domestication history was momentary compared with the long natural evolution; changing the transcription factors probably the easiest way happened to affect the agricultural traits of interest. however, putative candidate genes underlying these domestication-targeted phenotypes have diverse functions, which need to be validated by further experiments. plant-pathogen interaction affected by artificial selection domestication caused complex morphological and physiological changes in soybeans. annotated by the kegg and agrigo database, selective genes were associated with different biological functions, among which, plant-pathogen interaction, sequence-specific dna binding, phenylpropanoid biosynthesis, starch and sucrose metabolism are over-represented categories. plant-pathogen interactions are conducted between a pathogen and the host plant. in nature, plants are generally resistant to most invading pathogens due to innate ability to recognize them through successful defenses. when an exception happens, a pathogen would cause diseases in its host. however, pathogens could also cause diseases if they have evolved to evade detection or suppress host defense mechanisms, or both. the effects of plant-pathogen interactions are of particular relevance during early domestication events on agricultural systems. thus, understanding the genetic basis of why a certain pathogen causes disease in its host plant instead of others has long intrigued and motivated plant pathologists.figure functional annotation of selective genes with nucleotide fixation introduced in early domestication and modern improvement. a total of selective genes with nucleotide fixation were involved in plant-pathogen interactions. of them, two selective genes glyma14g and glyma08g with nucleotide fixation are orthologs of rps gene. the disease resistance gene rps was isolated using positional cloning and further screen for susceptible mutant. the rps protein contains two characteristics of a large family of plant r genes: a nucleotide-binding site and a leucine-rich repeat region. it is consistent with previous report that rps locus exhibit selection signals by examining a worldwide sample of arabidopsis accessions, and the n-terminal part of the leucine-rich repeat region was a probable target of selection. cyclic nucleotide-gated ion channels are ion channels that function in the pathogen signaling cascade by facilitating ca2+ uptake into the cytosol. two selective genes with nucleotide fixation were detected to encode cng channels. the topology of their proteins was predicted using tmhmm, which is based on a hidden markov model. the two genes encode transmembrane proteins with nucleotide fixation located outside the membrane. besides, eight selective genes are orthologs of transmembrane receptor kinase fls which acts as pathogen-associated molecular pattern signals in triggering the innate immune response. in addition, the category of terpene synthase activity was also enriched with six selective genes involved in. terpenes are one of the most important defensive plant compounds against herbivores and pathogens. recently, a new monoterpene synthase gene gmnes was identified and characterized in soybean. its transcription was up-regulated in soybean when infested with cotton leafworm. our analysis indicates the gene was possibly affected by artificial selection during soybean domestication. discussion nucleotide fixation was crucial in soybean divergence domestication led to significant morphological divergence between cultivated and wild soybeans. wild soybean exhibits, for example, twining and vine stem, severer shattering, impermeable seed coats, pod cracking sensitivity, small seeds, and low oleic acid, all of which are seldom observed in cultivars. deciphering how cultivated soybean have been transformed from its wild ancestor is advantageous both from genetic and evolutionary perspectives. with the available sequencing data, we comprehensively estimated the saturation number of snvs in soybean germplasm and detected a set of candidate genes showing artificial selection signals. to some degree, analysis of artificial selection and nucleotide fixation unravels the mystery of soybean domestication and subsequent improvement. based on nucleotide fixation, our analysis supports a single evolutionary origin of domesticated soybean. during domestication, only lines with certain agriculturally important traits were selected, resulting in a genome-wide reduction of genetic diversity or so-called selective sweep in cultivated crops. one possible explanation for the reduction is that an excess of nucleotide fixation happened in cultivars compared to wild soybeans. meta-analysis of qtls responsible for domestication related traits and the selective genes provided insights into the role of nucleotide fixation played in morphological differentiation between wild and cultivated soybeans. using comparative genomics, an amount of genes was found to be orthologs of those whose function was validated and responsible for corresponding traits in other plants. nucleotide fixation happened in those genes responsible for agronomically important traits. although traditional linkage and association mapping were used to dissect these traits, they failed to detect genetic changes caused by domestication and improvement. our analysis here provides valuable information for further qtl mapping and will facilitate molecular assisted selection in soybean breeding practice. artificial selection accelerates nucleotide fixation domestication was an evolutionary process where the characters of interests were selected, such as loss of seed dispersal, higher yield and increasing abiotic resistance. the detection of selective loci during crop domestication contributes to modern breeding efforts and the opportunity to improve genomic selection models. recently, genome-wide scans based on genetic bottlenecks have been successfully applied to detect footprints of selection in plants by surveying both natural and cultivated species. artificial selection of a beneficial mutation will lead to an elevated frequency in a population. eventually, allele frequencies will be skewed and nucleotide fixation happened after plant domestication. our analysis focused on to what degree nucleotide fixation was caused by artificial selection during soybean domestication. more nucleotide fixation happened in cultivars than those in wild soybeans, indicating that artificial selection was much stronger than natural selection. however, the effective population size of cultivated soybeans was substantially reduced during domestication, which could make a nucleotide seem to be fixed in cultivars. that mainly explains why nucleotide fixations were observed in cultivars across the soybean genome. considering nucleotide fixation accumulated in footprints of domestication and improvement, artificial selection probably accelerated the occurrence of fixation in soybean breeding activities. even thought, some of them could be also caused by the shrinking population size, especially when different haplotypes shaped in those selective sweeps. these fixations are extremely hard to be distinguished in current samples. morphological transition can be achieved by a mutation at a single locus, and artificial selection can rapidly change domestication targeted phenotypes within generations. domestication could be a rapid instead of a slow or gradual process, given strong selective pressures and a suitable genetic architecture. this was supported by the severe reduction of genetic diversity and large selective sweeps. in the process of domestication, any mutations detrimental to the traits of interests were eliminated immediately, whereas those advantageous ones were strongly selected, diffused and eventually fixed in a population. the environments wild soybeans grow in are various and usually harsh, resulting in diversifying selection instead of strong directional selection. what’s more, selection intensity imposed by natural selection was disparate in diverse habitats. these reasons also explain why artificial selection was much stronger than natural selection in crop domestication. evolutionary perspective of nucleotide fixation a long-term goal of crop genomics is to determine to what extent artificial selection impacts genomic variation patterns within and between populations. there are both genetic and statistical approaches to detect signals of hitchhiking caused selective sweeps. the hitchhiking effect is contingent on the nature of genetic variations and how selection acts on them. generally, there are at least three evolutionary routes by which a novel mutation may fix: drift to fixation for nearly neutral mutation; rapidly sweep to fixation, so-called hard sweep for beneficial mutation; and soft sweep to fixation for those initially neutral but later become beneficial for some reason. affected by artificial selection, a pre-exist mutation which became beneficial during domestication rapidly increased in frequency toward nucleotide fixation, as what we found in our analysis. when traits of interests during domestication were determined by multiple adaptive mutations at the same locus, artificial selection usually generates soft rather hard selective sweeps. many studies focus on hard sweeps in which only a single adaptive haplotype was skewed to fixation in the population, whereas multiple adaptive haplotypes formed simultaneously in a soft sweep. lots of nucleotide fixations happened within quantitative traits, indicating the corresponding traits of interests were incrementally changed at various causal loci. as a consequence, these sweeps related with artificial selection are likely to be both soft and incomplete. in soybean, some traits related to yield were selected, such as seed weight, seed blooming and prostrate habit, for which these are usually major qtls responsible. nevertheless, during intensive breeding human pursuits quality related traits such as protein content and lipid content, for which there are lots of small effect qtls responsible. analysis of nucleotide fixation indicates that more soft selective sweeps happened in extensive breeding than in early domestication in soybean, which still needs further investigation.  we integrated the available sequencing accessions to describe a whole picture of soybean genetic diversity, artificial selection and concomitant nucleotide fixation. there are approximately million snvs in soybean germplasm, of which about million reserved in cultivars. the genetic diversity was reduced by % in early domestication and subsequently reduced by % in genetic improvement. a total of and genes were involved in early domestication and subsequent improvement, respectively. both processes introduced about million nucleotide fixations, which contributed to the divergence of wild and cultivated soybeans. artificial selection probably accelerated the occurrence of nucleotide fixation, which affected some agronomic traits, as well as related biological pathways such as plant-pathogen interaction. 
analysis of available sequencing accessions estimates that ~ million single nucleotide variations reach saturation in cultivars, and then ~ million in soybean germplasm. artificial selection during soybean domestication and improvement results in substantial phenotypic divergence between wild and cultivated soybeans. nucleotide fixation mediated by artificial selection affected diverse molecular functions and biological reactions that associated with soybean morphological and physiological changes. soybean domesticationgenetic improvementartificial selectionnucleotide fixationplant-pathogen interactionissue-copyright-statement© the author 2015 strong selective pressure on beneficial phenotypes could cause nucleotide fixations in the founder population of soybean cultivars in quite a short time. meta-analysis of reported quantitative trait loci and selective signals with nucleotide fixation identifies a series of putative candidate genes responsible for agronomically important traits. of them, plant-pathogen interactions are of particular relevance as selective nucleotide fixations happened in disease resistance genes, cyclic nucleotide-gated ion channels and terpene synthases. selective sweeps defined by loss of genetic diversity reveal and genes were involved in domestication and subsequent improvement, respectively. both processes introduced ~ million nucleotide fixations, which contributed to the divergence of wild and cultivated soybeans. our analysis provides insights into the impacts of nucleotide fixation during soybean domestication and improvement, which would facilitate future qtl mapping and molecular breeding practice. plant domestication involves complex morphological and physiological modification of wild species to meet human needs.
 plant domestication involves complex morphological and physiological modification of wild species to meet human needs. artificial selection during soybean domestication and improvement results in substantial phenotypic divergence between wild and cultivated soybeans. strong selective pressure on beneficial phenotypes could cause nucleotide fixations in the founder population of soybean cultivars in quite a short time.  analysis of available sequencing accessions estimates that ~ million single nucleotide variations reach saturation in cultivars, and then ~ million in soybean germplasm. selective sweeps defined by loss of genetic diversity reveal and genes were involved in domestication and subsequent improvement, respectively. both processes introduced ~ million nucleotide fixations, which contributed to the divergence of wild and cultivated soybeans. meta-analysis of reported quantitative trait loci and selective signals with nucleotide fixation identifies a series of putative candidate genes responsible for agronomically important traits. nucleotide fixation mediated by artificial selection affected diverse molecular functions and biological reactions that associated with soybean morphological and physiological changes. of them, plant-pathogen interactions are of particular relevance as selective nucleotide fixations happened in disease resistance genes, cyclic nucleotide-gated ion channels and terpene synthases.  our analysis provides insights into the impacts of nucleotide fixation during soybean domestication and improvement, which would facilitate future qtl mapping and molecular breeding practice. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords soybean domesticationgenetic improvementartificial selectionnucleotide fixationplant-pathogen interactionissue-copyright-statement© the author 2015 