Beom Soon Choi

Genome sequence of the hot pepper provides insights into the evolution of pungency in Capsicum species

Hot pepper (Capsicum annuum), one of the oldest domesticated crops in the Americas, is the most widely grown spice crop in the world. We report whole-genome sequencing and assembly of the hot pepper (Mexican landrace of Capsicum annuum cv. CM334) at 186.6× coverage. We also report resequencing of two cultivated peppers and de novo sequencing of the wild species Capsicum chinense. The genome size of the hot pepper was approximately fourfold larger than that of its close relative tomato, and the genome showed an accumulation of Gypsy and Caulimoviridae family elements. Integrative genomic and transcriptomic analyses suggested that change in gene expression and neofunctionalization of capsaicin synthase have shaped capsaicinoid biosynthesis. We found differential molecular patterns of ripening regulators and ethylene synthesis in hot pepper and tomato. The reference genome will serve as a platform for improving the nutritional and medicinal values of Capsicum species.

Complete chloroplast and ribosomal sequences for 30 accessions elucidate evolution of Oryza AA genome species.

Cytoplasmic chloroplast (cp) genomes and nuclear ribosomal DNA (nR) are the primary sequences used to understand plant diversity and evolution. We introduce a high-throughput method to simultaneously obtain complete cp and nR sequences using Illumina platform whole-genome sequence. We applied the method to 30 rice specimens belonging to nine Oryza species. Concurrent phylogenomic analysis using cp and nR of several of specimens of the same Oryza AA genome species provides insight into the evolution and domestication of cultivated rice, clarifying three ambiguous but important issues in the evolution of wild Oryza species. First, cp-based trees clearly classify each lineage but can be biased by inter-subspecies cross-hybridization events during speciation. Second, O. glumaepatula, a South American wild rice, includes two cytoplasm types, one of which is derived from a recent interspecies hybridization with O. longistminata. Third, the Australian O. rufipogan-type rice is a perennial form of O. meridionalis.

Development of Reproducible EST-derived SSR Markers and Assessment of Genetic Diversity in Panax ginseng Cultivars and Related Species

Little is known about the genetics or genomics of Panax ginseng. In this study, we developed 70 expressed sequence tag-derived polymorphic simple sequence repeat markers by trials of 140 primer pairs. All of the 70 markers showed reproducible polymorphism among four Panax speciesand 19 of them were polymorphic in six P. ginseng cultivars. These markers segregated 1:2:1 manner of Mendelian inheritance in an F2 population of a cross between two P. ginseng cultivars, ‘Yunpoong’ and ‘Chunpoong’, indicating that these are reproducible and inheritable mappable markers. A phylogenetic analysis using the genotype data showed three distinctive groups: a P. ginseng-P. japonicus clade, P. notoginseng and P. quinquefolius, with similarity coefficients of 0.70. P. japonicus was intermingled with P. ginseng cultivars, indicating that both species have similar genetic backgrounds. P. ginseng cultivars were subdivided into three minor groups: an independent cultivar ‘Chunpoong’, a subgroup with three accessions including two cultivars, ‘Gumpoong’ and ‘Yunpoong’ and one landrace ‘Hwangsook’ and another subgroup with two accessions including one cultivar, ‘Gopoong’ and one landrace ‘Jakyung’. Each primer pair produced 1 to 4 bands, indicating that the ginseng genome has a highly replicated paleopolyploid genome structure.

Major repeat components covering one-third of the ginseng (Panax ginseng C.A. Meyer) genome and evidence for allotetraploidy

Ginseng (Panax ginseng) is a famous medicinal herb, but the composition and structure of its genome are largely unknown. Here we characterized the major repeat components and inspected their distribution in the ginseng genome. By analyzing three repeat-rich bacterial artificial chromosome (BAC) sequences from ginseng, we identified complex insertion patterns of 34 long terminal repeat retrotransposons (LTR-RTs) and 11 LTR-RT derivatives accounting for more than 80% of the BAC sequences. The LTR-RTs were classified into three Ty3/gypsy (PgDel, PgTat and PgAthila) and two Ty1/Copia (PgTork and PgOryco) families. Mapping of 30-Gbp Illumina whole-genome shotgun reads to the BAC sequences revealed that these five LTR-RT families occupy at least 34% of the ginseng genome. The Ty3/Gypsy families were predominant, comprising 74 and 33% of the BAC sequences and the genome, respectively. In particular, the PgDel family accounted for 29% of the genome and presumably played major roles in enlargement of the size of the ginseng genome. Fluorescence in situ hybridization (FISH) revealed that the PgDel1 elements are distributed throughout the chromosomes along dispersed heterochromatic regions except for ribosomal DNA blocks. The intensity of the PgDel2 FISH signals was biased toward 24 out of 48 chromosomes. Unique gene probes showed two pairs of signals with different locations, one pair in subtelomeric regions on PgDel2-rich chromosomes and the other in interstitial regions on PgDel2-poor chromosomes, demonstrating allotetraploidy in ginseng. Our findings promote understanding of the evolution of the ginseng genome and of that of related species in the Araliaceae.

Genotyping-by-sequencing map permits identification of clubroot resistance QTLs and revision of the reference genome assembly in cabbage (Brassica oleracea L.)

Clubroot is a devastating disease caused by Plasmodiophora brassicae and results in severe losses of yield and quality in Brassica crops. Many clubroot resistance genes and markers are available in Brassica rapa but less is known in Brassica oleracea. Here, we applied the genotyping-by-sequencing (GBS) technique to construct a high-resolution genetic map and identify clubroot resistance (CR) genes. A total of 43,821 SNPs were identified using GBS data for two parental lines, one resistant and one susceptible lines to clubroot, and 18,187 of them showed >5× coverage in the GBS data. Among those, 4,103 were credibly genotyped for all 78 F2 individual plants. These markers were clustered into nine linkage groups spanning 879.9 cM with an average interval of 1.15 cM. Quantitative trait loci (QTLs) survey based on three rounds of clubroot resistance tests using F2 : 3 progenies revealed two and single major QTLs for Race 2 and Race 9 of P. brassicae, respectively. The QTLs show similar locations to the previously reported CR loci for Race 4 in B. oleracea but are in different positions from any of the CR loci found in B. rapa. We utilized two reference genome sequences in this study. The high-resolution genetic map developed herein allowed us to reposition 37 and 2 misanchored scaffolds in the 02–12 and TO1000DH genome sequences, respectively. Our data also support additional positioning of two unanchored 3.3 Mb scaffolds into the 02–12 genome sequence.

Complete Genome Sequence of the Fenitrothion-Degrading Burkholderia sp. Strain YI23

Burkholderia species are ubiquitous in soil environments. Many Burkholderia species isolated from various environments have the potential to biodegrade man-made chemicals. Burkholderia sp. strain YI23 was isolated from a golf course soil and identified as a fenitrothion-degrading bacterium. In this study, we report the complete genome sequence of Burkholderia sp. strain YI23.

Complete Genome Sequence of Mycobacterium intracellulare Strain ATCC 13950T

Here we report the first complete genome sequence of Mycobacterium intracellulare ATCC 13950(T), a Mycobacterium avium complex (MAC) strain. This genome sequence will serve as a valuable reference for understanding the epidemiologic, biological, and pathogenic aspects of the disparity between MAC members.

Complete Genome Sequence of Mycobacterium intracellulare Clinical Strain MOTT-02

Here, we report the first complete genome sequence of the Mycobacterium intracellulare clinical strain MOTT-02, which was previously grouped in the INT2 genotype of M. intracellulare. This genome sequence will serve as a valuable reference for improving the understanding of the disparity in the virulence and epidemiologic traits between M. intracellulare genotypes.

Complete Genome Sequence of Mycobacterium intracellulare Clinical Strain MOTT-36Y, Belonging to the INT5 Genotype

Here we report the complete genome sequence of the Mycobacterium intracellulare clinical strain MOTT-36Y, previously grouped into the INT5 genotype among the 5 genotypes of M. intracellulare. This genome sequence will serve as a valuable reference for understanding the disparity in virulence and epidemiologic traits between M. intracellulare-related strains.

Survey of the Applications of NGS to Whole-Genome Sequencing and Expression Profiling

Recently, the technologies of DNA sequence variation and gene expression profiling have been used widely as approaches in the expertise of genome biology and genetics. The application to genome study has been particularly developed with the introduction of the next-generation DNA sequencer (NGS) Roche/454 and Illumina/Solexa systems, along with bioinformation analysis technologies of whole-genome de novo assembly, expression profiling, DNA variation discovery, and genotyping. Both massive whole-genome shotgun paired-end sequencing and mate paired-end sequencing data are important steps for constructing de novo assembly of novel genome sequencing data. It is necessary to have DNA sequence information from a multiplatform NGS with at least 2× and 30× depth sequence of genome coverage using Roche/454 and Illumina/Solexa, respectively, for effective an way of de novo assembly. Massive short-length reading data from the Illumina/Solexa system is enough to discover DNA variation, resulting in reducing the cost of DNA sequencing. Whole-genome expression profile data are useful to approach genome system biology with quantification of expressed RNAs from a whole-genome transcriptome, depending on the tissue samples. The hybrid mRNA sequences from Rohce/454 and Illumina/Solexa are more powerful to find novel genes through de novo assembly in any whole-genome sequenced species. The 20× and 50× coverage of the estimated transcriptome sequences using Roche/454 and Illumina/Solexa, respectively, is effective to create novel expressed reference sequences. However, only an average 30× coverage of a transcriptome with short read sequences of Illumina/Solexa is enough to check expression quantification, compared to the reference expressed sequence tag sequence.