Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

This article is part of the supplement: Selected articles from the Eleventh Asia Pacific Bioinformatics Conference (APBC 2013): Genomics

Open Access Proceedings

Genomic differences between cultivated soybean, G. max and its wild relative G. soja

Trupti Joshi1234, Babu Valliyodan35, Jeng-Hung Wu6, Suk-Ha Lee78, Dong Xu1234 and Henry T Nguyen235*

Author Affiliations

1 Department of Computer Science, University of Missouri, Columbia, MO 65211, USA

2 Christopher S. Bond Life Sciences Center, University of Missouri, Columbia, MO 65211, USA

3 National Center for Soybean Biotechnology, University of Missouri, Columbia, MO 65211, USA

4 Informatics Institute, University of Missouri, Columbia, MO 65211, USA

5 Division of Plant Sciences, University of Missouri, Columbia, MO 65211, USA

6 Department of Medicine, National Yang-Ming University, Taipei, Taiwan, R.O.C

7 Department of Plant Science and Research Institute for Agriculture and Life Sciences, Seoul National University, Seoul 151-921, Korea

8 Plant Genomics and Breeding Institute, Seoul National University, Seoul 151-921, Korea

For all author emails, please log on.

BMC Genomics 2013, 14(Suppl 1):S5  doi:10.1186/1471-2164-14-S1-S5

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/14/S1/S5


Published:21 January 2013

© 2013 Joshi et al.; licensee BioMed Central Ltd.

This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Glycine max is an economically important crop and many different varieties of soybean exist around the world. The first draft sequences and gene models of G. max (domesticated soybean) as well as G. soja (wild soybean), both became available in 2010. This opened the door for comprehensive comparative genomics studies between the two varieties.

Results

We have further analysed the sequences and identified the 425 genes that are unique to G. max and unavailable in G. soja. We further studied the genes with significant number of non-synonymous SNPs in their upstream regions. 12 genes involved in seed development, 3 in oil and 6 in protein concentration are unique to G. max. A significant number of unique genes are seen to overlap with the QTL regions of the three traits including seed, oil and protein. We have also developed a graphical chromosome visualizer as part of the Soybean Knowledge Base (SoyKB) tools for molecular breeding, which was used in the analysis and visualization of overlapping QTL regions for multiple traits with the deletions and SNPs in G. soja.

Conclusions

The comparisons between genome sequences of G. max and G. soja show significant differences between the genomic compositions of the two. The differences also highlight the phenotypic differences between the two in terms of seed development, oil and protein traits. These significant results have been integrated into the SoyKB resource and are publicly available for users to browse at http://soykb.org/GSoja webcite.

Introduction

Soybean [Glycine max (L.) Merr.] is a chief source of protein and oil for human and animal consumption. It has also become a major source of oil for biodiesel production, as well as a general model to study oil production in plants. The closest variety of soybean to G. max is the wild, undomesticated ancestor G. soja. The genome of Glycine soja has been recently sequenced [1] using Illumina Genome Analyzer and has been shown to have 915.4 MB consensus sequences, covering 97.65% of G. max genome sequence [2]. Even though G. soja is considered as the closest relative of the cultivated soybean, G. max, it has significant phenotypic differences. The large phenotypic diversity in soybean is genetically controlled, both qualitative and quantitative aspects. For example, G. soja has mainly tiny, black seeds in contrast to the large yellow seeds of G. max. Some QTL differences between G. max and G. soja have been reported [3-5]. Also, there are very significant differences in the seed oil and protein concentration between wild and domesticated soybeans. Our motivation behind this study is to identify the specific genomic differences that explain these phenotypic changes for the traits of our interest. In G. soja, the seed oil concentration is around 8% when compared to the 25% in G. max.

Several studies suggest that G. soja has important phenotypic characters and specific alleles which are not present in G. max [6]. Major traits of agricultural importance including yield and stress tolerance are polygenic and the presence of these favourable alleles in G. soja will help breeding program to introduce beneficial traits into soybean [7,8]. Earlier protein QTL studies are the best examples to show that G. soja alleles are unique and not present in the G. max lines. We have studied the differences between them with respect to things unique to each type as well as differences which make them different from the oil and protein traits perspective. We have also developed graphical tools to simplify the analysis of QTLs and traits, and overlay insertion and deletion regions of one over the other to highlight the genomic differences. This graphical chromosome visualizer is part of Breeding Program tools in SoyKB [9] web resources and provides extensive analysis and visualization capacity with genomics variation data such as SNPs, GWAS, insertions and deletions being overlaid on different QTL regions from multiple traits.

Materials and methods

Differences between G. max and G. soja

We acquired the sequences for all the 20 chromosomes and the gene model predictions for G. soja [1] as well as G. max [2]. We also acquired the list of inversions and deletions for G.soja. The chromosomal start and end positions for each of the 5794 deletions and 194 inversion positions were mapped against the G. max reference genome to identify the positions differing between the two. We acquired the sequences of 22,927 contigs assembled from sequences novel to G. soja. We also acquired the list of SNPs identified from the MAQ [10] analysis of the sequencing reads for further studies.

Unique to G. max

The 5795 deletion positions were compared against Glycine max gene positions for regions of overlap. The genes were considered to be absent in G. soja but present in G. max if the deleted position of sequences in Glycine soja overlapped with more than 50% of the overlapping annotated gene in Glycine max using Blast [11]. With this filtering criteria applied, we have identified 425 unique genes from the list of higher confidence ~46K gene model predictions from Phytozome [12] that are absent in G. soja and unique in G. max. Additional file 1 shows the entire list of these genes identified to be unique to G. max. When compared against the entire ~75k gene model predictions the number of genes unique to G. max is 2,533.

Additional file 1. Supp Table 1. List of 425 genes found to be unique to G. max when compared to G. soja.

Format: XLSX Size: 16KB Download fileOpen Data

The 425 genes unique to G. max were further analysed to look for overlap with transcription factor predictions from SoyDB [13], functional annotations and the different pathways available in SoyKB [9]. Out of the 425 unique genes, 313 have Pfam domain annotations and 132 can be mapped to various pathways. We identified 12 transcription factors (TFs) amongst this list of unique genes that are present in Glycine max, but absent in Glycine soja (as shown in Table 1). More than 50% of these TFs belong to the family of C2H2 zinc finger protein, one of the biggest family proteins in the eukaryotes. Members of this family show direct role in plant developmental processes including the embryo development [14] in plants. Also, it has been shown that the TF members from Myb and bHLH families are involved in the secondary axis of branching such as vegetative branches, inflorescence branches, or flowers [15].

Table 1. List of transcription factors unique to G. max.

Unique to G. soja

The sequences of 22,927 new contigs unique to G. soja (not present in G. max) were compared against the entire set of CDS sequences for Arabidopsis thaliana from TAIR10 [16], as well as Arabidopsis gene lists associated with seed, fatty acid metabolism, and protein metabolism and biosynthesis for sequence similarity. We used a filtering criteria for this comparison using Blast where hits passing the threshold of identity > = 80%, alignment length > = 100 bps, and e-value < 1e-3 were considered as significant hits. We identified 148, 132, 6, and 2 unique contigs matching significant hits from CDS, related to seed, fatty acid metabolism, and protein metabolism and biosynthesis, respectively. This was used for identifying the associated annotations and we also identified the pathways they were associated with.

SNPs analysis

The total number of all SNPs in G. soja identified from the MAQ analysis is 2,504,985 [1]. The number of non-synonymous SNPs amongst them is 56,948 and synonymous SNPs account for remaining 2,448,037. Out of these 23,379 non-synonymous SNPs in G. soja correspond to 10,770 unique genes annotated with gene IDs ending in ".1" and have one or multiple amino acid substitutions. Polyphen (Polymorphism Phenotyping) [17] and SIFT (Sorting Intolerant from Tolerant) [18] analysis was conducted for these data, to identify the significant non-synonymous SNPs affecting the protein characteristics. SNP positions with normalized probabilities less than 0.05 are predicted to be deleterious, those greater than or equal to 0.05 are predicted to be tolerated by the SIFT analysis.

We considered only the 66,365 Glyma1 gene model predictions with IDs ending in ".1" for further analysis. We compared the positions of these SNPs against the gene model positions and looked for SNPs overlapping within 1,000 bp upstream regions of these genes. Out of these SNPs, 174,208 SNPs were identified to be in regions 1 kb upstream of 46,270 genes (mean: 3.77 SNPs per kb, SD: 3.22). Out of these, 1,177 appear in 279 soybean genes (mean: 4.22, SD: 3.53) that are missing in G. soja but present in G. max, and remaining 173,031 SNPs overlap with 45,991 remaining common genes (mean: 3.77, SD: 3.21). The details are shown in Table 2 below. We further applied a stringent filtering criterion of number of SNPs in 1 kb upstream region > = 12 (p = 0.00467), to select these genes for further analysis. 14 genes (as listed below in Table 3) from 425 genes unique to G. max and 1533 genes from the remaining 45,991 common genes met the criterion and the annotations were studied for each of these genes.

Table 2. Distribution of SNPs positions amongst genes unique to G. max and common to both.

Table 3. 14 genes unique to G. max and having > = 12 SNPs in 1 kb upstream region.

Non-synonymous SNPs (in exons, upstream regions) analysis

The SNPs distribution from Kim et al.[1] was utilized in further analysis. The number of SNPs identified in exon regions is 86,236 and those identified in non-exon regions is 251,021. These were identified to be in 36,861 and 40,282 genes respectively. We calculated the difference between the SNP distribution observed in the 1 kb exon and non-exon regions to look for any distribution biases. We defined,

<a onClick="popup('http://www.biomedcentral.com/1471-2164/14/S1/S5/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2164/14/S1/S5/mathml/M1">View MathML</a>

Thus,

Difference in SNPs in 1 kb region = Average for 1 kb exon region- Average for 1 kb non-exon region.

2,405 genes met the filtering criteria with difference in SNPs in 1 kb region > 0 and # of SNPs in 1 kb exon regions > 10 (1% alteration). These 2,405 genes were further compared against the non-synonymous SNPs list to study and identify genes with higher # of non-synonymous SNPs. The ratio of # of non-synonymous over synonymous SNPs was calculated for this analysis.

<a onClick="popup('http://www.biomedcentral.com/1471-2164/14/S1/S5/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2164/14/S1/S5/mathml/M2">View MathML</a>

139 genes were identified that met the criteria of Ratio > 1. The pathways where these 139 genes were involved in were identified by using Gene Pathway Viewer in SoyKB [9]. Fatty acid biosynthesis, amino acid metabolism and biosynthesis, and protein trafficking and translation pathways are the major pathways these genes fall under.

Upstream region analysis

The 1 kb upstream sequence of 60 genes out of 425 unique genes in G. max with at least one non-synonymous deleterious SNPs and 139 genes with difference> 0, average> 10, and ratio of # of non-synonymous SNPs and # of synonymous SNPs> 1 were extracted and were further analysed using the PLACE database [19] to look for any cis-elements sequences in these regions. The numbers of SNPs in 1 kb upstream for each gene causing changes in the cis-elements were further studied and the detailed information of position and sequence of affected cis-elements was generated.

Homeologous genes analysis

9,807 duplicated homolog pairs from Libault et al. [20] were compared to 425 unique genes and 139 genes with non-synonymous SNPs in G. max as described above. 32 duplicated homolog pairs were identified with one homolog deleted and 18 duplicated homolog pairs were identified with one homolog with non-synonymous/synonymous SNPs ratio > 1. Table 4 and Table 5 show the list of these genes respectively. The functional annotations for all these genes were acquired using multiple gene search tool in SoyKB. Genes associated with storage protein biosynthesis and heat shock protein related chaperon classes were mostly deleted among these 32 homolog pairs. This result indicates that the wild soybean genome has retained a good portion of unique regions, when compared to G. max, and this suggests that G. max genome has undergone deletions during the domestication process.

Table 4. List of 32 genes in G. max with one homolog deleted.

Table 5. List of 18 genes in G. max with one homolog with non-synonymous/synonymous SNPs ratio > 1.

Genes involved in seed development, oil, protein and other functional groups of interest

We analyzed list of genes for various functions of our interest including seed development, oil-related, protein-related and some other functional groups of interest. The analysis was conducted to identify any unique genes in G. max involved in these functions that may have been deleted in G. soja. Figure 1 shows the positions of inversions, deletions, unique genes and TFs in G. max, along with the QTL regions for various traits in a circular genome viewer using Circos [21] tool.

thumbnailFigure 1. Circular genome viewer created using Circos. Circular genome viewer created using Circos, showing of inversions, deletions, unique genes and transcription factors in G. max, along with the QTL regions for various traits. The QTLs for oil trait is shown in orange, protein trait in red and seed trait in pink.

Seed development

A list of 16,365 unique loci involved in seed in Arabidopsis thaliana was extracted based on keyword search at TAIR [16]. These sequences were compared against the soybean cDNA sequences to identify the corresponding homolog using Blast (criteria: identity > = 80%, aligned length > = 100 bps, and e-value < 1e-3). When compared against the 425 soybean genes that are deleted in Glycine soja, 12 genes involved in seed development that are unique to G. max were identified as listed in Table 6.

Table 6. List of 12 genes involved in seed development that are unique to G. max.

Soybean gene Glyma04g12020.1 is annotated as WD40 repeat proteins and these proteins regulate biosynthesis of anthocyanins, proanthocyanidins (PAs), and mucilage in the seed and seed developmental stages [22]. Functional classification of these 12 genes clearly shows their roles in major biochemical processes during seed development in plants. Various molecular functions of these genes such as sugar and carbohydrate metabolism, nitrogen assimilation and metabolism, transporters, trans-membrane, and even the regulation of ubiquitination of proteins [23] are essential to regulate seed developmental processes and seed size in soybean. Above all, the Glyma14g02210.1 encoding enzyme CDP-diacylglycerol synthase is a major enzyme in regulation of fatty acyl fluxes during triacylglycerol (TAG) biosynthesis [24].

Oil-related

The sequences of 694 unique loci of known or genes considered to be involved in acyl lipid metabolism of Arabidopsis thaliana were compared against the soybean cDNA sequences to identify the corresponding homolog using Blast (criteria: identity > = 80% and e-value < 1e-3). When compared against the 425 soybean genes that are deleted in G. soja but present in G. max, 2 genes involved in acyl lipid metabolism were identified that are unique to G. max. 1,261 genes out of 1,266 unique hits to genes involved in acyl lipid metabolism were mapped to genes with ".1" gene models. As shown in Table 7, 3 genes were found in the list of 425 unique genes in G. max and one gene (Glyma10g05750.1) is novel when compared to the previous analysis. The gene Glyma06g39750.1 is shown to overlap with oil, seed-weight and yield related traits QTLs. The gene Glyma10g05750.1 is associated with southern root-knot nematode resistance trait and falls in its QTL. Glyma10g05750.1 is dolichyldiphosphatase which is catalytic and membrane related. This enzyme is reported as cell cycle dependent in human and, dolichol and dolichyl phosphate phosphatase (polyphenols) were characterized from soybean previously [25]. The lipid fractions of the lipo-sugar carriers are essential for the biosynthesis of N-linked glycoproteins. The gene Glyma06g39750.1 encodes for 3-hydroxyacyl-coa dehydrogenase and this enzyme catalyzes the third reaction of the mitochondrial β-oxidation cascade in animals. The gene Glyma04g12260 encodes for the fatty acyl-coenzyme A reductases. Earlier studies has confirmed that the heterologous expression of fatty acyl-coenzyme A reductases FAR1, FAR4, and FAR5 in yeast are active alcohol-forming FARs with distinct, but overlapping, chain length specificities ranging from C18:0 to C24:0 [26]. This report indicates that Arabidopsis FAR1, FAR4, and FAR5 generate the fatty alcohols found in root, seed coat, and wound-induced leaf tissue. The soybean transcript for Phospholipase D (Glyma07g08740) belongs to a lipolytic enzyme subclass which catalyzes the hydrolysis and transesterification of glycerophospholipids at the terminal phosphodiester bond [27].

Table 7. List of genes involved in acyl lipid metabolism that are unique to G. max.

Protein-related

Out of the 425 unique gens is G. max, we have identified several candidates associated with amino acid biosynthesis, amino acid transporters, nitrogen metabolism and glutamate biosynthesis directly correlated with the protein biosynthesis in soybean. One of the major groups of genes belongs to the cellular amino acid biosynthetic processes, which are the basic building blocks and very crucial for protein biosynthesis and metabolism in plants. Some example are, Glyma19g16450 encoding for glutamate synthase, Glyma08g07140 encoding for glutamine amidotransferase class-I, Glyma10g05090 ArgJ family proteins (glutamate ornithine acetyl transaminases). Another group of transcripts predominantly related to the protein metabolism are the vacuolar transporting proteins. Zouher et al. [28] have shown that VSR1, VSR3 and VSR4, but not the remaining VSRs or RMRs, participate in vacuolar sorting of VAC2 in vegetative tissues, and 12S globulins and 2S albumins in seeds, an activity that is essential for seedling germination vigor. Genes encoding for ribosomal protein, eg. Glyma19g06460, are another abundant group came out of the analysis. Ribosomes are the molecular machine for protein biosynthesis and the large (50S) and small (30S) ribosomal subunits contain over 50 distinct ribosomal (r) proteins that interact with the rRNAs and with one another [29]. Transcript Glyma09g36700 encodes for the signal recognition factor. Zhang and Shan [30] reported that the signal recognition particle (SRP) is universally conserved cellular machinery responsible for delivering membrane and secretory proteins to the proper cellular destination. The precise mechanism by which fidelity is achieved by the SRP pathway within the in vivo environment is yet to be understood.

Others functional groups of interest

360 unique loci involved in ubiquitin-dependent catabolic process, 5 unique loci involved in SCF-dependent proteasomal ubiquitin-dependent protein catabolic process, 39 unique loci involved in regulation of translation, 1 unique locus involved in regulation of protein catabolic process, 4 unique loci involved in positive regulation of translation, 8 unique loci involved in negative regulation of translation, 79 unique loci involved in glycoprotein biosynthetic process, and 247 unique loci involved in cellular amino acid biosynthetic process were also mapped to soybean gene models using Blast with best hits with e-value < 1e-3 and 292, 4, 36, 1, 4, 7, 75, and 232 soybean genes were identified, respectively. When compared against the 425 soybean genes that are deleted in G. soja but present in G. max, no genes involved in protein metabolism and regulation were identified. When compared against the 139 significant non-synonymous genes, one gene (Glyma03g28480.1) involved in cellular amino acid biosynthetic process was identified.

QTL, traits and SNPs analysis

The list of non-synonymous SNPs analyzed by Polyphen and Sift was further studied and compared against the different QTLs for various traits such as seed, oil, protein, etc. The genes falling in these QTLs of interest were studied in further details for the difference between the G. max and G. soja. The QTL regions for all the three traits were acquired using Breeding Program and Chromosome Visualizer tools of SoyKB [9]. The total # of all SNPs in G. soja is 2,504,985. Out of these 594,328 overlap with seed trait QTLs; 316,447 overlap with oil trait QTLs; and 380,719 overlap with protein trait QTLs. Each of these traits and QTLs were further analyzed for identifying overlapping genes unique to G. max and deleted in G. soja.

When further analyzed, we identified 315 out of the 425 unique genes in G. max also overlapping with various QTL regions for different traits. Out of these, we identified 97 genes overlapping with the seed trait QTLs, 54 with the Oil trait QTLs and 65 with the Protein trait QTLs in G. max.

Non-synonymous SNPs overlapping QTLs

Analysis was further conducted to identify the overlap on the non-synonymous SNPs with the QTL regions of various traits. The 23,379 non-synonymous SNPs in G. soja corresponding to 10,770 unique genes were utilized for this analysis and compared against various traits and QTL regions. When compared against the unique list of 20,498 genes ending in ".1" and falling in various QTLs for all available traits, we identified non-synonymous SNPs in 20,031 unique genes. Non-synonymous SNPs are present in 3,221, 1,579 and 2,346 genes falling in the QTLs of seed, oil and protein traits, respectively.

SoyKB graphical chromosome visualizer

We performed some analysis and visualization of the QTLs and traits using the SoyKB web resource. The G. soja SNPs, inversions and deletions have been added to the SoyKB resource. The graphical chromosome visualizer is an innovative molecular breeding tool, developed as part of SoyKB Breeding Program Tools, which integrates the genomic variations data such as SNPs and GWAS alongside the traits, QTLs, molecular markers and underlying genes. It displays the data linearly for each chromosome along with other useful information such as QTLs for multiple traits, SNPs from multiple genotypes of soybean including G. soja, as well as insertion, inversion and deletions data from multiple genotypes and G. soja. The chromosome visualizer can be directly accessed via the Tools menu on SoyKB (http://soykb.org/visualizer.php webcite). This tool can be used for highly meaningful analysis for QTLs of specific traits and showing the overlapping genes and regions of deletion and SNPs in comparison to the Williams82 reference genome. The list of genes falling in any QTL regions for single or multiple traits can be saved using the "Save" option on the viewing panel or "Left-clicking" the highlighted QTL region. The panels can be dragged and dropped to rearrange the sequence of visualizing the data. We performed the combined analysis for all the 3 traits. For example, in seed trait QTLs, gene Glyma01g17340.1 is deleted in G. soja and thus unique to G. max as shown in Figure 2. There are 315 such other genes falling in QTLs and being unique to G. max.

thumbnailFigure 2. Chromosome visualizer tool in SoyKB showing the gene Glyma01g17340. Chromosome visualizer tool in SoyKB shows the gene Glyma01g17340.1, which is unique in G. max, overlapping in QTL region for seed trait as well as the deletion region in G. soja.

Multiple traits and QTL analysis

We selected the seed, oil and protein traits to visualize the QTLs for all three of them. We also selected the G. soja SNPs as well as inversion and deletions datasets from the bottom panels for viewing. Figure 3 shows the region of Gm09 chromosome with the overlapping QTLs for all the 3 traits and also the region overlapping with the G. soja deletions region. This clearly shows an example of genes unique to G. max and to be involved important functions related to seed, oil and protein traits. This tool can be used by researchers to identify, confirm as well as generate lists of genes falling in the QTL regions.

thumbnailFigure 3. Chromosome visualizer tool in SoyKB showing the region of chromosome Gm09. Chromosome visualizer tool in SoyKB shows the region of chromosome Gm09 with overlapping QTL regions for multiple traits. Oil trait is shown in pink, protein trait in blue and seed trait in grey.

Combined analysis of seed and oil traits shows a region of deletion in G. soja on chromosome Gm01:19,434,434 - 19,498,249, which overlaps with the QTL regions for both the traits. This involves 7 consecutive genes (Glyma09g16400, Glyma09g16410, Glyma09g16420, Glyma09g16430, Glyma09g16440, Glyma09g16470, and Glyma09g16480), which are deleted in G. soja and are significant in G. max as they all fall in the QTL region for both oil and seed traits (Figure 3), contributing to the differences in the oil and seed development between the two varieties. Combined analysis of oil and protein traits shows a QTL region on chromosome Gm09 that contains several genes deleted in G. soja. Several genes in this region are significant since they overlap with some deleted regions in G. soja and are unique to G. max, in addition to falling in important QTL regions. In Figure 3, the oil trait QTL is highlighted in pink and protein QTL in blue. The functional annotation and experimental data for all of these genes were acquired using the multiple gene search feature of SoyKB.

Discussion/conclusions

With the availability of genome sequence data for both the varieties, it added tremendous valuable data that can be mined in an integrative and innovative manner to draw significant and interesting conclusions that can directly impact breeding of better soybeans. Some analysis has been done on the differences between G. soja and G. max earlier in terms of flowering, carbon metabolism and disease resistance [31]. However, no comparisons have been done in terms of differences in oil, protein and seed related functions and traits, especially at the genomic level. We have mined this data comprehensively for the traits of interest especially oil, seed development and protein. Our analysis shows there are some significant differences between Glycine max and Glycine soja at the genomic level, which corresponds with some functionally important genes for seed, oil and protein traits. These genes can explain some of the phenotypic differences observed between G. soja and G. max especially in terms of the above three traits. More experimental validations will follow including knockout studies of these genes in G. max to see the associated changes in phenotype. It is understood that the genomic differences identified in this study were revealed by comparing genome sequence of single genotype of G. soja and G. max. There are variations that exist within different populations of G. soja and G. max, respectively. Further studies comparing more diverse populations of both G. soja and G. max will be needed and will be extremely valuable in confirming some of these differences identified from a single population study. Nevertheless, the methodology demonstrates strongly how the genomic sequencing data can be combined with QTL and trait information to further narrow unique gene candidates involved in specific functional processes.

As part of the analysis more novel and innovative tools were developed as part of Breeding Program suite of tools in SoyKB, which are beneficial towards analysing and visualizing various traits and genes overlapping with them. This is a resource which is available to all communities freely and so our development can be utilized by others in the soybean area including researchers, scientists, breeders and farmers to their advantage. The overlapping of other genomic variations data along with the above information makes drawing inferences very easy and makes it the most informative tool available thus far.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

TJ conducted the bioinformatics analysis and downstream data analysis for QTL and traits. TJ is also the lead developer for the SoyKB web resource and drafted the initial manuscript. BV provided the expert advice for selection of traits and functional categories to study. Both TJ and BV contributed equally to the work. JW performed the sequence similarity analysis and SNPs comparison analysis. SL collaborated with the team and provided the sequences and associated files for G. soja. HTN and DX provided guidance and valuable insight. All authors have read and approved the manuscript.

Declarations

The funding for the publication of this article came from United Soybean Board (project 8236).

This article has been published as part of BMC Genomics Volume 14 Supplement 1, 2013: Selected articles from the Eleventh Asia Pacific Bioinformatics Conference (APBC 2013): Genomics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/14/S1 webcite.

Acknowledgements

The authors wish to thank Sujitha Murikoli who helped with generating figures in Circos. The development has been supported by Missouri Soybean Merchandising Council (MSMC #306) to DX, HTN; United Soybean Board (project 8236) to HTN and DX; and the National Center for Soybean Biotechnology. SL is supported by the Next Generation BioGreen 21 Program (Code No. PJ00811701), Rural Development Administration, Republic of Korea.

References

  1. Kim MY, Lee S, Van K, Kim TH, Jeong SC, Choi IY, Kim DS, Lee YS, Park D, Ma J, Kim WY, Kim BC, Park S, Lee KA, Kim DH, Kim KH, Shin JH, Jang YE, Kim KD, Liu WX, Chaisan T, Kang YJ, Lee YH, Kim KH, Moon JK, Schmutz J, Jackson SA, Bhak J, Lee SH: Whole-genome sequencing and intensive analysis of the undomesticated soybean (Glycine sojaSieb. and Zucc.) genome.

    Proc Natl Acad Sci USA 2010, 107(51):22032-7. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Schmutz J, Cannon SB, Schlueter J, Ma J, Mitros T, Nelson W, Hyten DL, Song Q, Thelen JJ, Cheng J, Xu D, Hellsten U, May GD, Yu Y, Sakurai T, Umezawa T, Bhattacharyya MK, Sandhu D, Valliyodan B, Lindquist E, Peto M, Grant D, Shu S, Goodstein D, Barry K, Futrell-Griggs M, Abernathy B, Du J, Tian Z, Zhu L, Gill N, Joshi T, Libault M, Sethuraman A, Zhang XC, Shinozaki K, Nguyen HT, Wing RA, Cregan P, Specht J, Grimwood J, Rokhsar D, Stacey G, Shoemaker RC, Jackson SA: Genome sequence of the palaeopolyploid soybean.

    Nature 2010, 463(7278):178-83. PubMed Abstract | Publisher Full Text OpenURL

  3. Liu B, Fujita T, Yan ZH, Sakamoto S, Xu D, Abe J: QTL mapping of domestication-related traits in soybean (Glycine max).

    Ann Bot 2007, 100:1027-1038. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Zhang WK, Wang YJ, Luo GZ, Zhang JS, He CY, Wu XL, Gai JY, Chen SY: QTL mapping of ten agronomic traits on the soybean (Glycine max L. Merr.) genetic map and their association with EST markers.

    Theor Appl Genet 2004, 108:1131-1139. PubMed Abstract | Publisher Full Text OpenURL

  5. Kang ST, Kwak M, Kim HK, Choung MG, Han WY, Baek IY, Kim MY, Van K, Lee SH: Population-specific QTLs and their different epistatic interactions for pod dehiscence in soybean (Glycine max (L.) Merr.).

    Euphytica 2009, 166:15-24. Publisher Full Text OpenURL

  6. Carter TE Jr, Nelson R, Sneller CH, Cui Z: Genetic diversity in soybean. In Soybeans: Improvement, Production and Uses. Edited by Boerma HR, Specht JE. American Society of Agronomy, Madison, WI; 303-416. OpenURL

  7. Tanksley SD, McCouch SR: Seed banks and molecular maps: unlocking genetic potential from the wild.

    Science 1997, 277:1063-1066. PubMed Abstract | Publisher Full Text OpenURL

  8. Li D, Pfeiffer TW, Cornelius PL: Soybean QTL for yield and yield components associated with Glycine soja alleles.

    Crop Sci 2008, 48:571-581. Publisher Full Text OpenURL

  9. Joshi T, Patil K, Fitzpatrick MR, Franklin LD, Yao Q, Cook JR, Wang Z, Libault M, Brechenmacher L, Valliyodan B, Wu X, Cheng J, Stacey G, Nguyen HT, Xu D: Soybean Knowledge Base (SoyKB): a web resource for soybean translational genomics.

    BMC Genomics 2012, 13(Suppl 1):S15. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  10. MAQ[http://maq.sourceforge.net/maq-man.shtml] webcite

  11. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool.

    J Mol Biol 1990, 215(3):403-10. PubMed Abstract | Publisher Full Text OpenURL

  12. Phytozome[http://www.phytozome.net/soybean] webcite

  13. Wang Z, Libault M, Joshi T, Valliyodan B, Nguyen H, Xu D, Stacey G, Cheng J: SoyDB: A Knowledge Database of Soybean Transcription Factors.

    BMC Plant Biology 2010, 10:14. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  14. Nakamura M, Runko AP, Sagerström CG: A novel subfamily of zinc finger genes involved in embryonic development.

    J Cell Biochem 2004, 93(5):887-95. PubMed Abstract | Publisher Full Text OpenURL

  15. Schmitz G, Theres K: Shoot and inflorescence branching.

    Curr Opin Plant Biol 2005, 8(5):506-11. PubMed Abstract | Publisher Full Text OpenURL

  16. Huala E, Dickerman A, Garcia-Hernandez M, Weems D, Reiser L, LaFond F, Hanley D, Kiphart D, Zhuang J, Huang W, Mueller L, Bhattacharyya D, Bhaya D, Sobral B, Beavis B, Somerville C, Rhee SY: The Arabidopsis Information Resource (TAIR): A comprehensive database and web-based information retrieval, analysis, and visualization system for a model plant.

    Nucleic Acids Res 2001, 29(1):102-5. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Ramensky V, Bork P, Sunyaev S: Human non-synonymous SNPs: server and survey.

    Nucleic Acids Res 2002, 30:3894-3900. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. Ng PC, Henikoff S: SIFT: predicting amino acid changes that affect protein function.

    Nucleic Acids Res 2003, 31:3812-3814. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Higo K, Ugawa Y, Iwamoto M, Korenaga T: Plant cis-acting regulatory DNA elements (PLACE) database.

    Nucleic Acids Res 1999, 27(1):297-300. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Libault M, Farmer A, Joshi T, Takahashi K, Langley RJ, Franklin LD, He J, Xu D, May G, Stacey G: An integrated transcriptome atlas of the crop model Glycine max, and its use in comparative analyses in plants.

    Plant J 2010, 63(1):86-99. PubMed Abstract | Publisher Full Text OpenURL

  21. Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA: Circos: an information aesthetic for comparative genomics.

    Genome Res 2009, 19(9):1639-45. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Pang Y, Wenger JP, Saathoff K, Peel GJ, Wen J, Huhman D, Allen SN, Tang Y, Cheng X, Tadege M, Ratet P, Mysore KS, Sumner LW, Marks MD, Dixon RA: A WD40 repeat protein from Medicago truncatula is necessary for tissue-specific anthocyanin and proanthocyanidin biosynthesis but not for trichome development.

    Plant Physiol 2009, 151(3):1114-29. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  23. Park BS, Song JT, Seo HS: Arabidopsis nitrate reductase activity is stimulated by the E3 SUMO ligase AtSIZ1.

    Nat Commun 2011, 2:400. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Hu Z, Ren Z, Lu C: The phosphatidylcholine diacylglycerol cholinephosphotransferase is required for efficient hydroxy fatty acid accumulation in transgenic Arabidopsis.

    Plant Physiol 2012, 158(4):1944-1954. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  25. Ravi K, Rip JW, Carroll KK: Characterization of dolichol and dolichyl phosphate phosphatase from soya beans (Glycine max).

    Biochem J 1983, 213(2):513-518. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  26. Domergue F, Vishwanath SJ, Joubès J, Ono J, Lee JA, Bourdon M, Alhattab R, Lowe C, Pascal S, Lessire R, Rowland O: Three Arabidopsis fatty acyl-coenzyme A reductases, FAR1, FAR4, and FAR5, generate primary fatty alcohols associated with suberin deposition.

    Plant Physiol 2010, 153(4):1539-1554. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  27. Lee J, Welti R, Roth M, Schapaugh WT, Li J, Trick HN: Enhanced seed viability and lipid compositional changes during natural ageing by suppressing phospholipase Dα in soybean seed.

    Plant Biotechnol J 2012, 10(2):164-173. PubMed Abstract | Publisher Full Text OpenURL

  28. Zouhar J, Muñoz A, Rojo E: Functional specialization within the vacuolar sorting receptor family: VSR1, VSR3 and VSR4 sort vacuolar storage cargo in seeds and vegetative tissues.

    Plant J 2010, 64(4):577-588. PubMed Abstract | Publisher Full Text OpenURL

  29. Yutin N, Puigbò P, Koonin EV, Wolf YI: Phylogenomics of prokaryotic ribosomal proteins.

    PLoS One 2012, 7(5):e36972. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  30. Zhang D, Shan SO: Translation elongation regulates substrate selection by the signal recognition particle.

    J Biol Chem 2012, 287(10):7652-7660. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  31. Moon Young Kim, Kyujung Van, Yang Jae Kang, Kil Hyun Kim, Suk-Ha Lee: Tracing soybean domestication history: From nucleotide to genome.

    Breed Sci 2012, 61(5):445-452. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL