A Genome-Wide Association Study of Protein, Oil, and Amino Acid Content in Wild Soybean (Glycine soja).

Soybean (Glycine max L.) is a globally important source of plant proteins, oils, and amino acids for both humans and livestock. Wild soybean (Glycine soja Sieb. and Zucc.), the ancestor of cultivated soybean, could be a useful genetic source for increasing these components in soybean crops. In this study, 96,432 single-nucleotide polymorphisms (SNPs) across 203 wild soybean accessions from the 180K Axiom® Soya SNP array were investigated using an association analysis. Protein and oil content exhibited a highly significant negative correlation, while the 17 amino acids exhibited a highly significant positive correlation with each other. A genome-wide association study (GWAS) was conducted on the protein, oil, and amino acid content using the 203 wild soybean accessions. A total of 44 significant SNPs were associated with protein, oil, and amino acid content. Glyma.11g015500 and Glyma.20g050300, which contained SNPs detected from the GWAS, were selected as novel candidate genes for the protein and oil content, respectively. In addition, Glyma.01g053200 and Glyma.03g239700 were selected as novel candidate genes for nine of the amino acids (Ala, Asp, Glu, Gly, Leu, Lys, Pro, Ser, and Thr). The identification of the SNP markers related to protein, oil, and amino acid content reported in the present study is expected to help improve the quality of selective breeding programs for soybeans.

Genetic Diversity of Korean Wild Soybean Core Collections and Genome-Wide Association Study for Days to Flowering.

The utilization of wild soybean germplasms in breeding programs increases genetic diversity, and they contain the rare alleles of traits of interest. Understanding the genetic diversity of wild germplasms is essential for determining effective strategies that can improve the economic traits of soybeans. Undesirable traits make it challenging to cultivate wild soybeans. This study aimed to construct a core subset of 1467 wild soybean accessions of the total population and analyze their genetic diversity to understand their genetic variations. Genome-wild association studies were conducted to detect the genetic loci underlying the time to flowering for a core subset collection, and they revealed the allelic variation in E genes for predicting maturity using the available resequencing data of wild soybean. Based on principal component and cluster analyses, 408 wild soybean accessions in the core collection covered the total population and were explained by 3 clusters representing the collection regions, namely, Korea, China, and Japan. Most of the wild soybean collections in this study had the E1e2E3 genotype according to association mapping and a resequencing analysis. Korean wild soybean core collections can provide helpful genetic resources to identify new flowering and maturity genes near the E gene loci and genetic materials for developing new cultivars, facilitating the introgression of genes of interest from wild soybean.

The discovery of novel SNPs associated with group A soyasaponin biosynthesis from Korea soybean core collection.

Soyasaponin is a type of glycoside such as steroids, steroidal alkaloids or triterpenes, which enhance the body immunity. In order to efficiently identify genes and markers related to the soyasaponin, we used a 180K Axiom® SoyaSNP array and whole genome resequencing data from the Korean soybean core collection. As a result of conducting GWAS for group A soyasaponin (Aa and Ab derivatives), 16 significant common markers associated with Aa and Ab derivatives were mapped to chromosome 7, and three candidate genes including Glyma.07g254600 were detected. The functional haplotypes for candidate genes showed that Aa and Ab contents were mainly determined by alleles of AX-90322128, the marker of Glyma.07g254600. In addition, 14 novel SNPs variants closely associated with Aa and Ab derivatives were discovered for Glyma.07g254600. Therefore, the results of this study that identified soyasaponin-associated markers and useful genes utilizing various genomic information could provide insight into functional soybean breeding.

Genome assembly of the popular Korean soybean cultivar Hwangkeum.

Massive resequencing efforts have been undertaken to catalog allelic variants in major crop species including soybean, but the scope of the information for genetic variation often depends on short sequence reads mapped to the extant reference genome. Additional de novo assembled genome sequences provide a unique opportunity to explore a dispensable genome fraction in the pan-genome of a species. Here, we report the de novo assembly and annotation of Hwangkeum, a popular soybean cultivar in Korea. The assembly was constructed using PromethION nanopore sequencing data and two genetic maps and was then error-corrected using Illumina short-reads and PacBio SMRT reads. The 933.12 Mb assembly was annotated as containing 79,870 transcripts for 58,550 genes using RNA-Seq data and the public soybean annotation set. Comparison of the Hwangkeum assembly with the Williams 82 soybean reference genome sequence (Wm82.a2.v1) revealed 1.8 million single-nucleotide polymorphisms, 0.5 million indels, and 25 thousand putative structural variants. However, there was no natural megabase-scale chromosomal rearrangement. Incidentally, by adding two novel subfamilies, we found that soybean contains four clearly separated subfamilies of centromeric satellite repeats. Analyses of satellite repeats and gene content suggested that the Hwangkeum assembly is a high-quality assembly. This was further supported by comparison of the marker arrangement of anthocyanin biosynthesis genes and of gene arrangement at the Rsv3 locus. Therefore, the results indicate that the de novo assembly of Hwangkeum is a valuable additional reference genome resource for characterizing traits for the improvement of this important crop species.

Fine-mapping and candidate gene analysis for the foxglove aphid resistance gene Raso2 from wild soybean PI 366121.

KEY MESSAGE: The foxglove aphid resistance gene Raso2 from PI 366121 was fine-mapped to 77 Kb region, and one candidate gene was identified. The foxglove aphid (FA: Aulacorthum solani Kaltenbach) is an important insect pest that causes serious yield losses in soybean. The FA resistance gene Raso2 from wild soybean PI 366121 was previously mapped to a 13 cM interval on soybean chromosome 7. However, fine-mapping of Raso2 was needed to improve the effectiveness of marker-assisted selection (MAS) and to eventually clone it. The objectives of this study were to fine-map Raso2 from PI 366121 using Axiom® 180 K SoyaSNP array, to confirm the resistance and inheritance of Raso2 in a different background, and to identify candidate gene(s). The 105 F4:8 recombinant inbred lines were used to fine-map the gene and to test antibiosis and antixenosis of Raso2 to FA. These efforts resulted in the mapping of Raso2 on 1 cM interval which corresponds to 77 Kb containing eight annotated genes based on the Williams 82 reference genome assembly (Wm82.a2.v1). Interestingly, all nonsynonymous substitutions were in Glyma.07g077700 which encodes the disease resistance protein containing LRR domain and expression of the gene in PI 366121 was significantly higher than that in Williams 82. In addition, distinct SNPs within Glyma.07g077700 that can distinguish PI 366121 and diverse FA-susceptible soybeans were identified. We also confirmed that Raso2 presented the resistance to FA and the Mendelian inheritance for single dominant gene in a different background. The results of this study would provide fundamental information on MAS for development of FA-resistant cultivars as well as functional study and cloning of the candidate gene in soybean.

The patterns of deleterious mutations during the domestication of soybean.

Globally, soybean is a major protein and oil crop. Enhancing our understanding of the soybean domestication and improvement process helps boost genomics-assisted breeding efforts. Here we present a genome-wide variation map of 10.6 million single-nucleotide polymorphisms and 1.4 million indels for 781 soybean individuals which includes 418 domesticated (Glycine max), 345 wild (Glycine soja), and 18 natural hybrid (G. max/G. soja) accessions. We describe the enhanced detection of 183 domestication-selective sweeps and the patterns of putative deleterious mutations during domestication and improvement. This predominantly selfing species shows 7.1% reduction of overall deleterious mutations in domesticated soybean relative to wild soybean and a further 1.4% reduction from landrace to improved accessions. The detected domestication-selective sweeps also show reduced levels of deleterious alleles. Importantly, genotype imputation with this resource increases the mapping resolution of genome-wide association studies for seed protein and oil traits in a soybean diversity panel.

Elucidation of genomic organizations of transgenic soybean plants through de novo genome assembly with short paired-end reads.

Elucidation of the genomic organizations of transgene insertion sites is essential for the genetic studies of transgenic plants. Herein, we establish an analysis pipeline that identifies the transgene insertion sites as well as the presence of vector backbones, through de novo genome assembly with high-throughput sequencing data in two transgenic soybean lines, AtYUCCA6-#5 and 35S-UGT72E3/2-#7. Sequencing data of approximately 28× and 29× genome coverages for each line generated by high-throughput sequencing were de novo assembled. The databases generated from the de novo assembled sequences were used to search contigs that contained putative insertion sites and their flanking sequences (integration sites) of transgene fragments using transgenic vector sequences as queries. The predicted integration site sequences, which are located at three annotated genes that might regulate plant development or confer disease resistance, were then confirmed by local alignment against the soybean reference genome and PCR amplification. As results, we revealed the precise transgene-flanking sequences and sequence rearrangements at insertion sites in both the transgenic lines, as well as the aberrant insertion of a transgene fragment. Consequently, relative to experimental or enrichment technologies, our approach is straightforward and time-effective, providing an alternative method for the identification of insertion sites in transgenic plants.

Environmental risk assessment of glufosinate-resistant soybean by pollen-mediated gene flow under field conditions in the region of the genetic origin.

Pollen-mediated gene flow of genetically modified crops to their wild relatives can facilitate the spread of transgenes into the ecosystem and alter the fitness of the consequential progeny. A two-year field study was conducted to quantify the gene flow from glufosinate-ammonium resistant (GR) soybean (Glycinemax) to its wild relative, wild soybean (G. soja), and assess the potential weed risk of hybrids resulting from the gene flow during their entire life cycle under field conditions in Korea, where wild soybean is the natural inhabitant. Pollen-mediated gene flow from GR soybeans to wild soybeans ranged from 0.292% (mixed planting) to 0.027% at 8 m distance. The log-logistic model described the gene flow rate with increasing distance from GR soybean to wild soybean; the estimated effective isolation distance for 0.01% gene flow between GR and wild soybeans was 37.7 m. The F1 and F2 hybrids exhibited the intermediate characteristics of their parental soybeans in their vegetative and reproductive stages. Canopy height and stem length of hybrids were close to those of wild soybean, which shows an indeterminate growth; the numbers of flowers, pods, and seeds per hybrid plant were close to those of wild soybean and significantly higher than those of GR soybean. Seed longevity of F2 hybrid plants was also intermediate but significantly greater than that of GR soybean due to high seed dormancy. Our results suggest that transgenes of the GR soybean might disperse into wild populations and persist in the agroecosystem of the genetic origin regions due to the pollen-mediated gene flow and the relatively high fitness of the hybrid progeny.

Genome Resequencing, Improvement of Variant Calling, and Population Genomic Analyses Provide Insights into the Seedlessness in the Genus Vitis.

The seedlessness of grape derived from stenospermocarpy is one of the most prized traits of table or raisin grapes. It is controlled by a complex genetic system containing one dominant gene and multiple recessive genes. Here, we collected dense variation data from high-depth resequencing data of seeded, seedless, and wild relative grape genomes sequenced to > 37x mean depth. Variant calls were made using a modified variant calling pipeline that was suitable for highly diverse interspecific grape accessions. The modified pipeline enabled us to call several million more variants than the commonly recommended pipeline. The quality was validated by Sanger sequencing data and subsequently supported by the genetic population structure and the phylogenetic tree constructed using the obtained variation data, results of which were generally consistent with known pedigree and taxonomic classifications. Variation data enabled us to confirm a dominant gene and identify recessive loci for seedlessness. Incidentally, we found that grape cultivar Rizamat contains an ancestral chromosomal region of the dominant gene in Sultanina, a predominant seedlessness donor cultivar. Furthermore, we predicted new candidate causal genes including Vitvi01g00455, Vitvi08g01528, and Vitvi18g01237 associated with the recessive seedless-regulating loci, which showed high homology with genes that regulate seed development in Arabidopsis This study provides fundamental insights relevant to variant calling from genome resequencing data of diverse interspecific hybrid germplasms such as grape and will accelerate future efforts aimed at crop improvement.

Performance of hybrids between abiotic stress-tolerant transgenic rice and its weedy relatives under water-stressed conditions.

Gene transfer from transgenic crops to their weedy relatives may introduce undesired ecological consequences that can increase the fitness and invasiveness of weedy populations. Here, we examined the rate of gene flow from abiotic stress-tolerant transgenic rice that over-express AtCYP78A7, a gene encoding cytochrome P450 protein, to six weedy rice accessions and compared the phenotypic performance and drought tolerance of their hybrids over generations. The rate of transgene flow from AtCYP78A7-overexpressing transgenic to weedy rice varied between 0% and 0.0396%. F1 hybrids containing AtCYP78A7 were significantly taller and heavier, but the percentage of ripened grains, grain numbers and weight per plant were significantly lower than their transgenic and weedy parents. The homozygous and hemizygous F2 progeny showed higher tolerance to drought stress than the nullizygous F2 progeny, as indicated by leaf rolling scores. Shoot growth of nullizygous F3 progeny was significantly greater than weedy rice under water-deficient conditions in a rainout shelter, however, that of homozygous F3 progeny was similar to weedy rice, indicating the cost of continuous expression of transgene. Our findings imply that gene flow from AtCYP78A7-overexpressing transgenic to weedy rice might increase drought tolerance as shown in the pot experiment, however, increased fitness under stressed conditions in the field were not observed for hybrid progeny containing transgenes.

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

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

Korean soybean core collection: Genotypic and phenotypic diversity population structure and genome-wide association study.

A core collection is a subset that represents genetic diversity of the total collection. Soybean (Glycine max (L.) Merr.) is one of major food and feed crops. It is the world's most cultivated annual herbaceous legume. Constructing a core collection for soybean could play a pivotal role in conserving and utilizing its genetic variability for research and breeding programs. To construct and evaluate a Korean soybean core collection, genotypic and phenotypic data as well as population structure, were analyzed. The Korean soybean core collection consisted of 430 accessions selected from 2,872 collections based on Affymetrix Axiom® 180k SoyaSNP array data. The core collection represented 99% of genotypic diversity of the total collection. Analysis of population structure clustered the core collection into five subpopulations. Accessions from South Korea and North Korea were distributed across five subpopulations. Analysis of molecular variance indicated that only 2.01% of genetic variation could be explained by geographic origins while 16.18% of genetic variation was accounted for by subpopulations. Genome-wide association study (GWAS) for days to flowering, flower color, pubescent color, and growth habit confirmed that the core collection had the same genetic diversity for tested traits as the total collection. The Korean soybean core collection was constructed based on genotypic information of the 180k SNP data. Size and phenotypic diversity of the core collection accounted for approximately 14.9% and 18.1% of the total collection, respectively. GWAS of core and total collections successfully confirmed loci associated with tested traits. Consequently, the present study showed that the Korean soybean core collection could provide fundamental and practical material and information for both soybean genetic research and breeding programs.

GmBRC1 is a Candidate Gene for Branching in Soybean (Glycine max (L.) Merrill).

Branch number is one of the main factors affecting the yield of soybean (Glycine max (L.)). In this study, we conducted a genome-wide association study combined with linkage analysis for the identification of a candidate gene controlling soybean branching. Five quantitative trait nucleotides (QTNs) were associated with branch numbers in a soybean core collection. Among these QTNs, a linkage disequilibrium (LD) block qtnBR6-1 spanning 20 genes was found to overlap a previously identified major quantitative trait locus qBR6-1. To validate and narrow down qtnBR6-1, we developed a set of near-isogenic lines (NILs) harboring high-branching (HB) and low-branching (LB) alleles of qBR6-1, with 99.96% isogenicity and different branch numbers. A cluster of single nucleotide polymorphisms (SNPs) segregating between NIL-HB and NIL-LB was located within the qtnBR6-1 LD block. Among the five genes showing differential expression between NIL-HB and NIL-LB, BRANCHED1 (BRC1; Glyma.06G210600) was down-regulated in the shoot apex of NIL-HB, and one missense mutation and two SNPs upstream of BRC1 were associated with branch numbers in 59 additional soybean accessions. BRC1 encodes TEOSINTE-BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTORS 1 and 2 transcription factor and functions as a regulatory repressor of branching. On the basis of these results, we propose BRC1 as a candidate gene for branching in soybean.

Genetic diversity patterns and domestication origin of soybean.

KEY MESSAGE: Genotyping data of a comprehensive Korean soybean collection obtained using a large SNP array were used to clarify global distribution patterns of soybean and address the evolutionary history of soybean. Understanding diversity and evolution of a crop is an essential step to implement a strategy to expand its germplasm base for crop improvement research. Accessions intensively collected from Korea, which is a small but central region in the distribution geography of soybean, were genotyped to provide sufficient data to underpin population genetic questions. After removing natural hybrids and duplicated or redundant accessions, we obtained a non-redundant set comprising 1957 domesticated and 1079 wild accessions to perform population structure analyses. Our analysis demonstrates that while wild soybean germplasm will require additional sampling from diverse indigenous areas to expand the germplasm base, the current domesticated soybean germplasm is saturated in terms of genetic diversity. We then showed that our genome-wide polymorphism map enabled us to detect genetic loci underlying flower color, seed-coat color, and domestication syndrome. A representative soybean set consisting of 194 accessions was divided into one domesticated subpopulation and four wild subpopulations that could be traced back to their geographic collection areas. Population genomics analyses suggested that the monophyletic group of domesticated soybeans was likely originated at a Japanese region. The results were further substantiated by a phylogenetic tree constructed from domestication-associated single nucleotide polymorphisms identified in this study.

Differences in the metabolic profiles and antioxidant activities of wild and cultivated black soybeans evaluated by correlation analysis.

Wild soybeans are considered a potential resource for soybean domestication and an important source of genetic diversity for soybean crop improvement. Understanding metabolite-caused bioactivity differences between cultivated and wild soybeans is essential for designing a soybean with enhanced nutritional traits. In this study, the non-targeted metabolic profiling of 26 soybean varieties, 15 wild black soybeans (WBS) and 11 cultivated black soybeans (CBS), using liquid chromatography-mass spectrometry (LC-MS) in combination with multivariate analysis revealed significant differences in 25 differential metabolites. Among these, the soyasaponins Ab and Bb were found to be characteristic metabolites expressed more substantially in CBS than in WBS. Three different antioxidant assays and correlation analysis identified major and minor antioxidants that contributed to WBS having an antioxidant activity 4- to 8-fold stronger than that of CBS. Epicatechin, procyanidin B2, and cyanidin-3-O-glucoside were identified by both association analysis and the online LC-ABTS radical scavenging assay as being major antioxidants.

GenoCore: A simple and fast algorithm for core subset selection from large genotype datasets.

Selecting core subsets from plant genotype datasets is important for enhancing cost-effectiveness and to shorten the time required for analyses of genome-wide association studies (GWAS), and genomics-assisted breeding of crop species, etc. Recently, a large number of genetic markers (>100,000 single nucleotide polymorphisms) have been identified from high-density single nucleotide polymorphism (SNP) arrays and next-generation sequencing (NGS) data. However, there is no software available for picking out the efficient and consistent core subset from such a huge dataset. It is necessary to develop software that can extract genetically important samples in a population with coherence. We here present a new program, GenoCore, which can find quickly and efficiently the core subset representing the entire population. We introduce simple measures of coverage and diversity scores, which reflect genotype errors and genetic variations, and can help to select a sample rapidly and accurately for crop genotype dataset. Comparison of our method to other core collection software using example datasets are performed to validate the performance according to genetic distance, diversity, coverage, required system resources, and the number of selected samples. GenoCore selects the smallest, most consistent, and most representative core collection from all samples, using less memory with more efficient scores, and shows greater genetic coverage compared to the other software tested. GenoCore was written in R language, and can be accessed online with an example dataset and test results at https://github.com/lovemun/Genocore.

Genetic analysis of the Lf1 gene that controls leaflet number in soybean.

KEY MESSAGE: A high-resolution genetic map that was constructed for the Lf1 -residing region will provide valuable information for map-based cloning and genetic improvement efforts in soybean. Changes in leaf architecture as photosynthesis factories remain a major challenge for the improvement of crop productivity. Unlike most soybeans, which have compound leaves comprising three leaflets, the soybean Lf1 mutant has a high frequency of compound leaves with five leaflets in a partially dominant manner. Here, we generated a fine genetic map to determine the genetic basis of this multifoliolate leaf trait. A five-leaflet variant Dusam was found in a recently collected landrace cultivar. Phenotypic data were collected from the F2 population of a cross between the Dusam and three-leaflet cultivar V94-5152. The mapping results generated using public markers indicated that the five-leaflet determining gene in Dusam is an allele of the previously studied Lf1 gene on chromosome 8. A high-resolution map delimited the genomic region controlling the leaflet number trait to a sequence length of 49 kb. AP2 domain-containing Glyma.08g281900 annotated in this 49 kb region appeared to be a strong candidate for the Lf1-encoding gene, as members of the AP2-type transcription factor family regulate lateral organ development. Dusam additionally exhibits visually distinct phenotypes for shattering and seed-coat cracking traits. However, the two traits were clearly unlinked to the Lf1 gene in our mapping population. Interestingly, the mapping results suggest that the Lf1 gene most likely exerts a pleiotropic effect on the number of seeds per pod. Thus, our results provide a strong foundation towards the cloning of this compound leaf development gene and marker-assisted selection of the seeds per pod trait.

Identification of haplotypes at the Rsv4 genomic region in soybean associated with durable resistance to soybean mosaic virus.

KEY MESSAGE: Discovery of new germplasm sources and identification of haplotypes for the durable Soybean mosaic virus resistance gene, Rsv 4, provide novel resources for map-based cloning and genetic improvement efforts in soybean. The Soybean mosaic virus (SMV) resistance locus Rsv4 is of interest because it provides a durable type of resistance in soybean [Glycine max (L.) Merr.]. To better understand its molecular basis, we used a population of 309 BC3F2 individuals to fine-map Rsv4 to a ~120 kb interval and leveraged this genetic information in a second study to identify accessions 'Haman' and 'Ilpumgeomjeong' as new sources of Rsv4. These two accessions along with three other Rsv4 and 14 rsv4 accessions were used to examine the patterns of nucleotide diversity at the Rsv4 region based on high-depth resequencing data. Through a targeted association analysis of these 19 accessions within the ~120 kb interval, a cluster of four intergenic single-nucleotide polymorphisms (SNPs) was found to perfectly associate with SMV resistance. Interestingly, this ~120 kb interval did not contain any genes similar to previously characterized dominant disease resistance genes. Therefore, a haplotype analysis was used to further resolve the association signal to a ~94 kb region, which also resulted in the identification of at least two Rsv4 haplotypes. A haplotype phylogenetic analysis of this region suggests that the Rsv4 locus in G. max is recently introgressed from G. soja. This integrated study provides a strong foundation for efforts focused on the cloning of this durable virus resistance gene and marker-assisted selection of Rsv4-mediated SMV resistance in soybean breeding programs.

Development, validation and genetic analysis of a large soybean SNP genotyping array.

Cultivated soybean (Glycine max) suffers from a narrow germplasm relative to other crop species, probably because of under-use of wild soybean (Glycine soja) as a breeding resource. Use of a single nucleotide polymorphism (SNP) genotyping array is a promising method for dissecting cultivated and wild germplasms to identify important adaptive genes through high-density genetic mapping and genome-wide association studies. Here we describe a large soybean SNP array for use in diversity analyses, linkage mapping and genome-wide association analyses. More than four million high-quality SNPs identified from high-depth genome re-sequencing of 16 soybean accessions and low-depth genome re-sequencing of 31 soybean accessions were used to select 180,961 SNPs for creation of the Axiom(®) SoyaSNP array. Validation analysis for a set of 222 diverse soybean lines showed that 170,223 markers were of good quality for genotyping. Phylogenetic and allele frequency analyses of the validation set data indicated that accessions showing an intermediate morphology between cultivated and wild soybeans collected in Korea were natural hybrids. More than 90 unanchored scaffolds in the current soybean reference sequence were assigned to chromosomes using this array. Finally, dense average spacing and preferential distribution of the SNPs in gene-rich chromosomal regions suggest that this array may be suitable for genome-wide association studies of soybean germplasm. Taken together, these results suggest that use of this array may be a powerful method for soybean genetic analyses relating to many aspects of soybean breeding.

Cross-family translational genomics of abiotic stress-responsive genes between Arabidopsis and Medicago truncatula.

Cross-species translation of genomic information may play a pivotal role in applying biological knowledge gained from relatively simple model system to other less studied, but related, genomes. The information of abiotic stress (ABS)-responsive genes in Arabidopsis was identified and translated into the legume model system, Medicago truncatula. Various data resources, such as TAIR/AtGI DB, expression profiles and literatures, were used to build a genome-wide list of ABS genes. tBlastX/BlastP similarity search tools and manual inspection of alignments were used to identify orthologous genes between the two genomes. A total of 1,377 genes were finally collected and classified into 18 functional criteria of gene ontology (GO). The data analysis according to the expression cues showed that there was substantial level of interaction among three major types (i.e., drought, salinity and cold stress) of abiotic stresses. In an attempt to translate the ABS genes between these two species, genomic locations for each gene were mapped using an in-house-developed comparative analysis platform. The comparative analysis revealed that fragmental colinearity, represented by only 37 synteny blocks, existed between Arabidopsis and M. truncatula. Based on the combination of E-value and alignment remarks, estimated translation rate was 60.2% for this cross-family translation. As a prelude of the functional comparative genomic approaches, in-silico gene network/interactome analyses were conducted to predict key components in the ABS responses, and one of the sub-networks was integrated with corresponding comparative map. The results demonstrated that core members of the sub-network were well aligned with previously reported ABS regulatory networks. Taken together, the results indicate that network-based integrative approaches of comparative and functional genomics are important to interpret and translate genomic information for complex traits such as abiotic stresses.