Kil Hyun Kim

Whole-genome sequencing and intensive analysis of the undomesticated soybean (Glycine soja Sieb. and Zucc.) genome

The genome of soybean (Glycine max), a commercially important crop, has recently been sequenced and is one of six crop species to have been sequenced. Here we report the genome sequence of G. soja, the undomesticated ancestor of G. max (in particular, G. soja var. IT182932). The 48.8-Gb Illumina Genome Analyzer (Illumina-GA) short DNA reads were aligned to the G. max reference genome and a consensus was determined for G. soja. This consensus sequence spanned 915.4 Mb, representing a coverage of 97.65% of the G. max published genome sequence and an average mapping depth of 43-fold. The nucleotide sequence of the G. soja genome, which contains 2.5 Mb of substituted bases and 406 kb of small insertions/deletions relative to G. max, is ∼0.31% different from that of G. max. In addition to the mapped 915.4-Mb consensus sequence, 32.4 Mb of large deletions and 8.3 Mb of novel sequence contigs in the G. soja genome were also detected. Nucleotide variants of G. soja versus G. max confirmed by Roche Genome Sequencer FLX sequencing showed a 99.99% concordance in single-nucleotide polymorphism and a 98.82% agreement in insertion/deletion calls on Illumina-GA reads. Data presented in this study suggest that the G. soja/G. max complex may be at least 0.27 million y old, appearing before the relatively recent event of domestication (6,000∼9,000 y ago). This suggests that soybean domestication is complicated and that more in-depth study of population genetics is needed. In any case, genome comparison of domesticated and undomesticated forms of soybean can facilitate its improvement.

Genome sequence of mungbean and insights into evolution within Vigna species

Mungbean (Vigna radiata) is a fast-growing, warm-season legume crop that is primarily cultivated in developing countries of Asia. Here we construct a draft genome sequence of mungbean to facilitate genome research into the subgenus Ceratotropis, which includes several important dietary legumes in Asia, and to enable a better understanding of the evolution of leguminous species. Based on the de novo assembly of additional wild mungbean species, the divergence of what was eventually domesticated and the sampled wild mungbean species appears to have predated domestication. Moreover, the de novo assembly of a tetraploid Vigna species (V. reflexo-pilosa var. glabra) provides genomic evidence of a recent allopolyploid event. The species tree is constructed using de novo RNA-seq assemblies of 22 accessions of 18 Vigna species and protein sets of Glycine max. The present assembly of V. radiata var. radiata will facilitate genome research and accelerate molecular breeding of the subgenus Ceratotropis.

Genome-wide mapping of NBS-LRR genes and their association with disease resistance in soybean

BackgroundR genes are a key component of genetic interactions between plants and biotrophic bacteria and are known to regulate resistance against bacterial invasion. The most common R proteins contain a nucleotide-binding site and a leucine-rich repeat (NBS-LRR) domain. Some NBS-LRR genes in the soybean genome have also been reported to function in disease resistance. In this study, the number of NBS-LRR genes was found to correlate with the number of disease resistance quantitative trait loci (QTL) that flank these genes in each chromosome. NBS-LRR genes co-localized with disease resistance QTL. The study also addressed the functional redundancy of disease resistance on recently duplicated regions that harbor NBS-LRR genes and NBS-LRR gene expression in the bacterial leaf pustule (BLP)-induced soybean transcriptome.ResultsA total of 319 genes were determined to be putative NBS-LRR genes in the soybean genome. The number of NBS-LRR genes on each chromosome was highly correlated with the number of disease resistance QTL in the 2-Mb flanking regions of NBS-LRR genes. In addition, the recently duplicated regions contained duplicated NBS-LRR genes and duplicated disease resistance QTL, and possessed either an uneven or even number of NBS-LRR genes on each side. The significant difference in NBS-LRR gene expression between a resistant near-isogenic line (NIL) and a susceptible NIL after inoculation of Xanthomonas axonopodis pv. glycines supports the conjecture that NBS-LRR genes have disease resistance functions in the soybean genome.ConclusionsThe number of NBS-LRR genes and disease resistance QTL in the 2-Mb flanking regions of each chromosome was significantly correlated, and several recently duplicated regions that contain NBS-LRR genes harbored disease resistance QTL for both sides. In addition, NBS-LRR gene expression was significantly different between the BLP-resistant NIL and the BLP-susceptible NIL in response to bacterial infection. From these observations, NBS-LRR genes are suggested to contribute to disease resistance in soybean. Moreover, we propose models for how NBS-LRR genes were duplicated, and apply Ks values for each NBS-LRR gene cluster.

RNA-Seq Analysis of a Soybean Near-Isogenic Line Carrying Bacterial Leaf Pustule-Resistant and -Susceptible Alleles

Bacterial leaf pustule (BLP) disease is caused by Xanthomonas axonopodis pv. glycines (Xag). To investigate the plant basal defence mechanisms induced in response to Xag, differential gene expression in near-isogenic lines (NILs) of BLP-susceptible and BLP-resistant soybean was analysed by RNA-Seq. Of a total of 46 367 genes that were mapped to soybean genome reference sequences, 1978 and 783 genes were found to be up- and down-regulated, respectively, in the BLP-resistant NIL relative to the BLP-susceptible NIL at 0, 6, and 12h after inoculation (hai). Clustering analysis revealed that these genes could be grouped into 10 clusters with different expression patterns. Functional annotation based on gene ontology (GO) categories was carried out. Among the putative soybean defence response genes identified (GO:0006952), 134 exhibited significant differences in expression between the BLP-resistant and -susceptible NILs. In particular, pathogen-associated molecular pattern (PAMP) and damage-associated molecular pattern (DAMP) receptors and the genes induced by these receptors were highly expressed at 0 hai in the BLP-resistant NIL. Additionally, pathogenesis-related (PR)-1 and -14 were highly expressed at 0 hai, and PR-3, -6, and -12 were highly expressed at 12 hai. There were also significant differences in the expression of the core JA-signalling components MYC2 and JASMONATE ZIM-motif. These results indicate that powerful basal defence mechanisms involved in the recognition of PAMPs or DAMPs and a high level of accumulation of defence-related gene products may contribute to BLP resistance in soybean.

Draft genome sequence of adzuki bean, Vigna angularis

Adzuki bean (Vigna angularis var. angularis) is a dietary legume crop in East Asia. The presumed progenitor (Vigna angularis var. nipponensis) is widely found in East Asia, suggesting speciation and domestication in these temperate climate regions. Here, we report a draft genome sequence of adzuki bean. The genome assembly covers 75% of the estimated genome and was mapped to 11 pseudo-chromosomes. Gene prediction revealed 26,857 high confidence protein-coding genes evidenced by RNAseq of different tissues. Comparative gene expression analysis with V. radiata showed that the tissue specificity of orthologous genes was highly conserved. Additional re-sequencing of wild adzuki bean, V. angularis var. nipponensis, and V. nepalensis, was performed to analyze the variations between cultivated and wild adzuki bean. The determined divergence time of adzuki bean and the wild species predated archaeology-based domestication time. The present genome assembly will accelerate the genomics-assisted breeding of adzuki bean.

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.

Tracing soybean domestication history: From nucleotide to genome

Since the genome sequences of wild species may provide key information about the genetic elements involved in speciation and domestication, the undomesticated soybean (Glycine soja Sieb. and Zucc.), a wild relative of the current cultivated soybean (G. max), was sequenced. In contrast to the current hypothesis of soybean domestication, which holds that the current cultivated soybean was domesticated from G. soja, our previous work has suggested that soybean was domesticated from the G. soja/G. max complex that diverged from a common ancestor of these two species of Glycine. In this review, many structural genomic differences between the two genomes are described and a total of 705 genes are identified as structural variations (SVs) between G. max and G. soja. After protein families database of alignments and hidden Markov models IDs and gene ontology terms were assigned, many interesting genes are discussed in detail using four domestication related traits, such as flowering time, transcriptional factors, carbon metabolism and disease resistance. Soybean domestication history is explored by studying these SVs in genes. Analysis of SVs in genes at the population-level may clarify the domestication history of soybean.

Genome-wide comparative analysis of flowering genes between Arabidopsis and mungbean

Mungbean is one of the major crops grown in South, East and Southeast Asia because of a high quality of amino acid profile; however, its asynchronous flowering time makes difficult to harvest at a time. Synchronization of flowering time is important to reduce labor costs for harvesting. With the availability of next generation sequencing data of mungbean, we approached a strategy of comparative genomics to identify mungbean homologous counterpart of A. thaliana genes that are known to be involved in flowering pathways, followed by a comparative soybean quantitative trait loci (QTL) analysis of the putative mungbean flowering-related genes. Co-localization of mungbean QTL associated with days to first flower day (FLD) was also identified using the EST-SSR markers from a previous study. Additionally, based on the mungbean transcriptome data with a distinct flowering stage of R2, FPKM (Fragments Per Kilobase per Million mapped reads) expression analysis of all the genes found in paralogous synteny blocks was conducted to examine expression patterns of the genes that have undergone a whole genome duplication event. Our results indicate that the paralogous flowering genes along with other genes within a same synteny block have evolved together at the macro-synteny scale. This study provides insights into mungbean flowering genes, in which they can be used as tools in order to improve flowering synchronization and to increase yield.

Genetic diversity and population structure of wild soybean ( Glycine soja Sieb. and Zucc.) accessions in Korea

Genetic variation in wild soybean ( Glycine soja Sieb. and Zucc.) is a valuable resource for crop improvement efforts. Soybean is believed to have originated from China, Korea, and Japan, but little is known about the diversity or evolution of Korean wild soybean. Therefore, in this study, we evaluated the genetic diversity and population structure of 733 G. soja accessions collected in Korea using 21 simple sequence repeat (SSR) markers. The SSR loci produced 539 alleles (25.7 per locus) with a mean genetic diversity of 0.882 in these accessions. Rare alleles, those with a frequency of less than 5%, represented 75% of the total number. This collection was divided into two populations based on the principal coordinate analysis. Accessions from population 1 were distributed throughout the country, whereas most of the accessions from population 2 were distributed on the western side of the Taebaek and Sobaek mountains. The Korean G. soja collection evaluated in this study should provide useful background information for allele mining approach and breeding programmes to introgress alleles into the cultivated soybean ( G. max (L). Merr.) from wild soybean.

Optimization of a Virus-Induced Gene Silencing System with Soybean yellow common mosaic virus for Gene Function Studies in Soybeans

Virus-induced gene silencing (VIGS) is an effective tool for the study of soybean gene function. Successful VIGS depends on the interaction between virus spread and plant growth, which can be influenced by environmental conditions. Recently, we developed a new VIGS system derived from the Soybean yellow common mosaic virus (SYCMV). Here, we investigated several environmental and developmental factors to improve the efficiency of a SYCMV-based VIGS system to optimize the functional analysis of the soybean. Following SYCMV: Glycine max-phytoene desaturase (GmPDS) infiltration, we investigated the effect of photoperiod, inoculation time, concentration of Agrobacterium inoculm, and growth temperature on VIGS efficiency. In addition, the relative expression of GmPDS between non-silenced and silenced plants was measured by qRT-PCR. We found that gene silencing efficiency was highest at a photoperiod of 16/8 h (light/dark) at a growth temperature of approximately 27°C following syringe infiltration to unrolled unifoliolate leaves in cotyledon stage with a final SYCMV:GmPDS optimal density (OD)600 of 2.0. Using this optimized protocol, we achieved high efficiency of GmPDS-silencing in various soybean germplasms including cultivated and wild soybeans. We also confirmed that VIGS occurred in the entire plant, including the root, stem, leaves, and flowers, and could transmit GmPDS to other soybean germplasms via mechanical inoculation. This optimized protocol using a SYCMV-based VIGS system in the soybean should provide a fast and effective method to elucidate gene functions and for use in large-scale screening experiments.