Kyung Do 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.

The genome sequences of Arachis duranensis and Arachis ipaensis , the diploid ancestors of cultivated peanut

Cultivated peanut (Arachis hypogaea) is an allotetraploid with closely related subgenomes of a total size of ∼2.7 Gb. This makes the assembly of chromosomal pseudomolecules very challenging. As a foundation to understanding the genome of cultivated peanut, we report the genome sequences of its diploid ancestors (Arachis duranensis and Arachis ipaensis). We show that these genomes are similar to cultivated peanuts A and B subgenomes and use them to identify candidate disease resistance genes, to guide tetraploid transcript assemblies and to detect genetic exchange between cultivated peanuts subgenomes. On the basis of remarkably high DNA identity of the A. ipaensis genome and the B subgenome of cultivated peanut and biogeographic evidence, we conclude that A. ipaensis may be a direct descendant of the same population that contributed the B subgenome to cultivated peanut.

Widespread natural variation of DNA methylation within angiosperms

BackgroundDNA methylation is an important feature of plant epigenomes, involved in the formation of heterochromatin and affecting gene expression. Extensive variation of DNA methylation patterns within a species has been uncovered from studies of natural variation. However, the extent to which DNA methylation varies between flowering plant species is still unclear. To understand the variation in genomic patterning of DNA methylation across flowering plant species, we compared single base resolution DNA methylomes of 34 diverse angiosperm species.ResultsBy analyzing whole-genome bisulfite sequencing data in a phylogenetic context, it becomes clear that there is extensive variation throughout angiosperms in gene body DNA methylation, euchromatic silencing of transposons and repeats, as well as silencing of heterochromatic transposons. The Brassicaceae have reduced CHG methylation levels and also reduced or loss of CG gene body methylation. The Poaceae are characterized by a lack or reduction of heterochromatic CHH methylation and enrichment of CHH methylation in genic regions. Furthermore, low levels of CHH methylation are observed in a number of species, especially in clonally propagated species.ConclusionsThese results reveal the extent of variation in DNA methylation in angiosperms and show that DNA methylation patterns are broadly a reflection of the evolutionary and life histories of plant species.

A comparative epigenomic analysis of polyploidy-derived genes in soybean and common bean

Reference methylomes provide insights into the evolutionary role of DNA methylation in paleopolyploid genomes. Soybean (Glycine max) and common bean (Phaseolus vulgaris) share a paleopolyploidy (whole-genome duplication [WGD]) event, approximately 56.5 million years ago, followed by a genus Glycine-specific polyploidy, approximately 10 million years ago. Cytosine methylation is an epigenetic mark that plays an important role in the regulation of genes and transposable elements (TEs); however, the role of DNA methylation in the fate/evolution of genes following polyploidy and speciation has not been fully explored. Whole-genome bisulfite sequencing was used to produce nucleotide resolution methylomes for soybean and common bean. We found that, in soybean, CG body-methylated genes were abundant in WGD genes, which were, on average, more highly expressed than single-copy genes and had slower evolutionary rates than unmethylated genes, suggesting that WGD genes evolve more slowly than single-copy genes. CG body-methylated genes were also enriched in shared single-copy genes (single copy in both species) that may be responsible for the broad and high expression patterns of this class of genes. In addition, diverged methylation patterns in non-CG contexts between paralogs were due mostly to TEs in or near genes, suggesting a role for TEs and non-CG methylation in regulating gene expression post polyploidy. Reference methylomes for both soybean and common bean were constructed, providing resources for investigating epigenetic variation in legume crops. Also, the analysis of methylation patterns of duplicated and single-copy genes has provided insights into the functional consequences of polyploidy and epigenetic regulation in plant genomes.

Dynamic Rearrangements Determine Genome Organization and Useful Traits in Soybean

Soybean (Glycine max) is a paleopolyploid whose genome has gone through at least two rounds of polyploidy and subsequent diploidization events. Several studies have investigated the changes in genome structure produced by the relatively recent polyploidy event, but little is known about the ancient polyploidy due to the high frequency of gene loss after duplication. Our previous study, regarding a region responsible for bacterial leaf pustule, reported two homeologous Rxp regions produced by the recent whole-genome duplication event. In this study, we identified the full set of four homeologous Rxp regions (ranging from 1.96 to 4.60 Mb) derived from both the recent and ancient polyploidy events, and this supports the quadruplicated structure of the soybean genome. Among the predicted genes on chromosome 17 (linkage group D2), 71% of them were conserved in a recently duplicated region, while 21% and 24% of duplicated genes were retained in two homeologous regions formed by the ancient polyploidy. Furthermore, comparative analysis showed a 2:1 relationship between soybean and Medicago truncatula, since M. truncatula did not undergo the recent polyploidy event that soybean did. Unlike soybean, M. truncatula homeologous regions were highly fractionated and their synteny did not exist, revealing different rates of diploidization process between the two species. Our data show that extensive synteny remained in the four homeologous regions in soybean, even though the soybean genome experienced dynamic genome rearrangements following paleopolyploidy events. Moreover, multiple Rxp quantitative trait loci on different soybean chromosomes actually comprise homeologous regions produced by two rounds of polyploidy events.

Transcriptomic changes due to water deficit define a general soybean response and accession-specific pathways for drought avoidance

BackgroundAmong abiotic stresses, drought is the most common reducer of crop yields. The slow-wilting soybean genotype PI 416937 is somewhat robust to water deficit and has been used previously to map the trait in a bi-parental population. Since drought stress response is a complex biological process, whole genome transcriptome analysis was performed to obtain a deeper understanding of the drought response in soybean.ResultsContrasting data from PI 416937 and the cultivar ‘Benning’, we developed a classification system to identify genes that were either responding to water-deficit in both genotypes or that had a genotype x environment (GxE) response. In spite of very different wilting phenotypes, 90% of classifiable genes had either constant expression in both genotypes (33%) or very similar response profiles (E genes, 57%). By further classifying E genes based on expression profiles, we were able to discern the functional specificity of transcriptional responses at particular stages of water-deficit, noting both the well-known reduction in photosynthesis genes as well as the less understood up-regulation of the protein transport pathway. Two percent of classifiable genes had a well-defined GxE response, many of which are located within slow-wilting QTLs. We consider these strong candidates for possible causal genes underlying PI 416937’s unique drought avoidance strategy.ConclusionsThere is a general and functionally significant transcriptional response to water deficit that involves not only known pathways, such as down-regulation of photosynthesis, but also up-regulation of protein transport and chromatin remodeling. Genes that show a genotypic difference are more likely to show an environmental response than genes that are constant between genotypes. In this study, at least five genes that clearly exhibited a genotype x environment response fell within known QTL and are very good candidates for further research into slow-wilting.

Landscape and evolutionary dynamics of terminal repeat retrotransposons in miniature in plant genomes

BackgroundTerminal repeat retrotransposons in miniature (TRIMs) are a unique group of small long terminal repeat retrotransposons that are difficult to identify. Thus far, only a few TRIMs have been characterized in the euphyllophytes, and their evolutionary and biological significance as well as their transposition mechanisms are poorly understood.ResultsUsing a combination of de novo and homology-based methods, we annotate TRIMs in 48 plant genome sequences, spanning land plants to algae. The TRIMs are grouped into 156 families including 145 that were previously undefined. Notably, we identify the first TRIMs in a lycophyte and non-vascular plants. The majority of the TRIM families are highly conserved and shared within and between plant families. Unlike other long terminal repeat retrotransposons, TRIMs are enriched in or near genes; they are also targeted by sRNAs between 21 and 24 nucleotides in length, and are frequently found in CG body-methylated genes. Importantly, we also identify putative autonomous retrotransposons and very recent transpositions of a TRIM element in Oryza sativa.ConclusionsWe perform the most comprehensive analysis of TRIM transposons thus far and report that TRIMs are ubiquitous across plant genomes. Our results show that TRIMs are more frequently associated with large and CG body-methylated genes that have undergone strong purifying selection. Our findings also indicate that TRIMs are likely derived from internal deletions of large long terminal repeat retrotransposons. Finally, our data and methodology are important resources for the characterization and evolutionary and genomic studies of long terminal repeat retrotransposons in other genomes.

A new approach for annotation of transposable elements using small RNA mapping

Transposable elements (TEs) are mobile genomic DNA sequences found in most organisms. They so densely populate the genomes of many eukaryotic species that they are often the major constituents. With the rapid generation of many plant genome sequencing projects over the past few decades, there is an urgent need for improved TE annotation as a prerequisite for genome-wide studies. Analogous to the use of RNA-seq for gene annotation, we propose a new method for de novo TE annotation that uses as a guide 24 nt-siRNAs that are a part of TE silencing pathways. We use this new approach, called TASR (for Transposon Annotation using Small RNAs), for de novo annotation of TEs in Arabidopsis, rice and soybean and demonstrate that this strategy can be successfully applied for de novo TE annotation in plants. Executable PERL is available for download from: http://tasr-pipeline.sourceforge.net/

Divergent cytosine DNA methylation patterns in single-cell, soybean root hairs.

Chromatin modifications, such as cytosine methylation of DNA, play a significant role in mediating gene expression in plants, which affects growth, development, and cell differentiation. As root hairs are single-cell extensions of the root epidermis and the primary organs for water uptake and nutrients, we sought to use root hairs as a single-cell model system to measure the impact of environmental stress. We measured changes in cytosine DNA methylation in single-cell root hairs as compared with multicellular stripped roots, as well as in response to heat stress. Differentially methylated regions (DMRs) in each methylation context showed very distinct methylation patterns between cell types and in response to heat stress. Intriguingly, at normal temperature, root hairs were more hypermethylated than were stripped roots. However, in response to heat stress, both root hairs and stripped roots showed hypomethylation in each context, especially in the CHH context. Moreover, expression analysis of mRNA from similar tissues and treatments identified some associations between DMRs, genes and transposons. Taken together, the data indicate that changes in DNA methylation are directly or indirectly associated with expression of genes and transposons within the context of either specific tissues/cells or stress (heat).

Genomic and epigenomic immunity in common bean: the unusual features of NB-LRR gene family

Abstract In plants, a key class of genes comprising most of disease resistance (R) genes encodes Nucleotide-binding leucine-rich repeat (NL) proteins. Access to common bean (Phaseolus vulgaris) genome sequence provides unparalleled insight into the organization and evolution of this large gene family (∼400 NL) in this important crop. As observed in other plant species, most common bean NL are organized in cluster of genes. However, a particularity of common bean is that these clusters are often located in subtelomeric regions close to terminal knobs containing the satellite DNA khipu. Phylogenetically related NL are spread between different chromosome ends, suggesting frequent exchanges between non-homologous chromosomes. NL peculiar location, in proximity to heterochromatic regions, led us to study their DNA methylation status using a whole-genome cytosine methylation map. In common bean, NL genes displayed an unusual body methylation pattern since half of them are methylated in the three contexts, reminiscent of the DNA methylation pattern of repeated sequences. Moreover, 90 NL were also abundantly targeted by 24 nt siRNA, with 90% corresponding to methylated NL genes. This suggests the existence of a transcriptional gene silencing mechanism of NL through the RdDM (RNA-directed DNA methylation) pathway in common bean that has not been described in other plant species.