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 novoassembly, expression profiling, DNA variation discovery, and genotyping. One of the advantages of the NGS systems is the cost-effectiveness to obtain the result of high-throughput DNA sequencing for genome, RNAnome, and miRNAnome studies. 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 and for resequencing the samples with a reference genome DNA sequence. To construct high-quality contig consensus sequences, each DNA fragment read length is important to obtain de novo assembly with long reading sequences of the Roche/454 system. 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, as hybrid assembly for novel genome sequencing would be cost-effective. In some cases, Illumina/Solexa data are used to construct scaffolds through de novo assembly with high coverage depth and large diverse fragment mate paired-end information,even though they are already participating in assembly and have made many contigs. Massive short-length reading data from the Illumina/Solexa system is enough to discover DNA variation, resulting in reducing the cost of DNA sequencing. MAQ and CLC software are useful to both single nucleotide polymorphism discovery and genotyping through a comparison of resequencing data to a reference genome. 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, such as control and exposed tissue. 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. In an in silicomethod, conserved miRNA and novel miRNA discovery is available on massive miRNAnome data in any species. Particularly, the discovered target genes of miRNA could be robust to approach genome biology study.

The differences in the immune response between body lice, Pediculus humanus humanus, and head lice, Pediculus humanus capitis, were regarded as primary factors determining their differential vector competence. To find any differences in genetic components in immune system between body and head lice, whole genome sequences of head lice were determined by both SBS [sequencing by synthesis, Illumina Genome Analyzer (Illumina-GA)] and pyrosequencing (Roche GS FLX), and compared with the reference genome sequences of body lice. The short DNA reads from Illumina-GA (an average mapping depth of 50-fold) were aligned first to the body louse reference genome, to which Roche GS FLX DNA reads (an average depth of 2.5-fold) were subsequently assembled to make up gaps between mapped consensus. Total consensus showed a size of 114 Mb and a coverage of 96% of the published body louse genome sequences. From this head louse genome sequences, a total of 12,651 genes were predicted and used for comparing with the 10,775 genes previously reported from the body louse genome. The homolog analysis identified 873 head louse-specific genes and 422 body lice-specific genes. Comparison of immune response genes between both louse species showed head lice have more number of immune-related genes than body lice. Head lice were determined to possess all of the 107 immune-related genes reported in the previous study (Kim et al., 2011), suggesting that there is no difference in genetic make-up in terms of the 107 immune-related genes between body and head lice.

Natural and artificially induced mutants have provided valuable resources for plant genetic studies and crop improvement. Some variations induced in the process of plant transformation have often been observed in regenerated plants. In this study, we investigated the insertion number of transgene and the flanking sequences of T-DNA in tall-induced line BP23, which was unexpectedly gained in the process of transformation of insect-resistant rice with cryBP1 gene, and also analyzed the whole-genome sequencing by using the NGS technologies to gain a better understanding of the sequence and structural changes between tall line or natural cultivar and rice reference. than others, was confirmed with two copies of foreign gene insertion, which was inserted in one genomic site facing each other between the position 2,430,152~2,430,151 of rice chromosome 12 without any deletion of genomic sequences. Sequencing analysis also revealed that 18bp-unknown sequences were added in the 5′ insertion site of T-DNA. This position in rice genome was confirmed with none of expressed gene sites. By the NGS analysis, we detected 86560 SNPs and 1091/1472 large insertion/deletion (indel) sites (100bp) between BP23 and rice reference, and 84743 SNPs and 1094/1451 large indels between natural cultivar Nagdong and rice reference. The possible mechanisms for the gene mutation, the developmental and tissue expression of the taller height in BP23 line may need to be scrutinized a few more.

Although the overall structure of the chloroplast genome is generally conserved, a number of sequence variations have been identified, which are valuable for plant population and evolutionary studies. Here, we constructed a chloroplast variation map of 30 landrace rice strains of Korean origin, using the Oryza rufipogon chloroplast genome (Genbank: NC_017835) as a reference. Differential distribution of single nucleotide polymorphisms (SNPs) and indels across the rice chloroplast genome is suggestive of a region-specific variation. Population structure clustering revealed the existence of two clear subgroups (indica and japonica) and an admixture group (aus). Phylogenetic analysis of the 30 landrace rice strains and six rice chloroplast references suggested and supported independent evolution of O. sativa indica and japonica. Interestingly, two “aus” type accessions, which were thought to be indica type, shared a closer relationship with the japonica type. One hypothesis is that “Korean aus” was intentionally introduced and may have obtained japonica chloroplasts during cultivation. We also calculated the nucleotide diversity of 30 accessions and compared to six rice chloroplast references, which shown a higher diversity in the indica and aus groups than in the japonica group in lower level substitution diversity.

The generation and analysis of genomic resources information are essential to understand genomic features of crops. Even though medicinal component and its effect of Panax ginseng was well studied, the genomic study has been recently started. The ginseng genome has been known to undergo two rounds of whole genome duplication (WGD), therefore we investigated an evidence of WGD in ginseng draft sequence for understanding current ginseng genome structure. Four paralogous gene-rich genome blocks were found, consisted by eight scaffolds, using about 3.0 Gb whole genome draft sequence and 48,821 unigenes of P. ginseng generated by whole genome shotgun sequencing. The eight scaffold sequences were ordered and connected into four genomic blocks, using zig-zag extension within scaffold sequences recently duplicated. The paralogous scaffold pairs that were recently duplicated showed high sequence conservation in genic and non-genic regions. However, paleo duplicated paralogue scaffold sequences showed little conservation only in genic regions. Finally, a total of 110 paralogous gene pairs and its expression were identified from recently and paleo duplicated scaffold pairs, which were co-linear among four genomic blocks. This study provides the first insight into duplicated genome structure of ginseng and will be a valuable information for further ginseng genomics including improvement of draft sequence quality, chromosome anchoring of scaffolds, and genetic mapping.

GM벼 OsCK는 벼 유래의 OsCK1 유전자를 벼에 형질전환 하여 벼흰잎마름병 및 벼도열병에 대한 저항성을 높게 한벼로 농촌진흥청에서 개발하였다. 형질전환 벡터의 구성은 양쪽 border (LB, RB) 상간에 2개의 MAR 염기서열이 서로 마주보는 형태로 위치하고 있으며, 제초제 저항성 유전자 PAT는 CaMV 35S promoter에 의하여 발현이 유도되고, 목표 유전자인 choline kinase (OsCK)는 actin promoter에 의하여 발현이 조절되며 left border 기준으로 역방향으로 배치되었다. 도입유전자 확인을 위하여 adaptor ligation PCR을 수행하였는데, MAR 영역에 위치하는 제한효소로 GM벼 genomic DNA를 절단한 후 adaptor를 붙였다. 염기서열 분석을 위하여 T-DNA의 양 말단에서 primer를 제작한 후 sequence 분석을 하였다. 분석한 결과, T-DNA의 right border 인근의 MAR sequence가 벼 genome의 10번 염색체 129971번 염기와 연결되어 있음을 확인하였다. Left 영역의 삽입위치는 이후 실시한 Illumina NGS 시스템을 이용하여 확인할 수 있었으며, GM 벼에는 2개의 T-DNA가 도입되었음을 알 수 있고, 첫 번째 T-DNA는 벼 10번 염색체 BAC클론 OSJNBa0014J14의 128947번째 염기와 129970째 염기에 위치하고 벼 genome 염기 1024 bp가 결실됨을 확인하였다. 이 과정에서 첫 번째 T-DNA left border와 첫 번째 MAR sequence 일부(370 bp)가 결실되었고 right border와 두 번째 MAR 영역 199 bp가 결실되었음도 확인하였다. 두 번째 T-DNA는 right border가 결실된 형태로 첫번째 T-DNA의 35S promoter 중간에 삽입되었음을 확인하였다.

Chloroplasts are plant-specific organelles, which have their own genome. Most of the plant chloroplast genomes (CP genome) are highly conserved in terms of its gene contents and genome structures, and they exist in cells with abundant copy numbers. Because of numerous copy numbers, the complete chloroplast sequence assembly pipeline with small amount of whole genome resequencing data, produced by NGS technique, was established in our laboratory. From 14 accessions of cabbage (Brassica oleracea L.) resequencing data produced by Illumina Hi-seq 2000, CP genomes were assembled and compared to each other. 18 sequence variance regions were detected, and 6 HRM(High Resolution Melting curves) markers were developed. Approximately 1 Gb of whole genome sequencing data of 10 Brassica rapa and 2 Brassica napus were also obtained from Institute of Vegetables and Flowers, Chinese Academy of Agricultural Science. With these resequencing data, all CP genomes from these accessions were assembled. Total 27 complete CP genomes of B.oleracea, B.rapa, B.napus, and brassico-raphanus which is a novel allotetraploid species between B.rapa and Raphanus sativus, were compared in sequence level. Phylogenetic analysis based on the comparison revealed that B.rapa could be the maternal species when rapeseeds and brassico-raphanus became allotetraploid species. Additionally, CP genome of B.napus cv.M083 is closer to B.rapa accessions than the other B.napus accessions, thus B.napus could have several different origins.

such as small plant size, small genome, short life cycle, self-fertility, and etc. Moreover, Brachypodium standard line Bd21 had already been sequenced and the data of which now have been available to the public. The development of next-generation sequencing technologies has allowed the discovery of large numbers of genome-wide DNA polymorphisms. Brachypodium standard line Bd21 was exposed to chronic gamma radiation. M2 1376-1 line was less stained with phloroglucinol compared to wild-type, which indicated reduced lingin content. Also, it’s filial generations showed dwarf phenotype. Genomic DNA was extracted from the M3 plant and was used to construct a whole-genome re-sequencing library for using Illumina HiSeq2500. The trimmed reads were aligned to the Brachypodium reference genome sequence (http:// www.phytozome.net). SICKLE, BWA, and Picard were used for accurate variant detection. More than 110 M reads were generated and 96.53% of them were mapped. This represents ~35 fold coverage. As a result, mutant specific SNVs, Insert/Deletions, and non-synonymous mutations were obtained. Moreover, non-synonymous mutations were identified from 5 lignin biosynthesis genes (Bradi1g31320, Bradi3g52350, Bradi5g21550, Bradi3g22980, Bradi5g14720). The obtained results will be incorporated in development of biofuel crops.

Rice is the major food for half of the world population. The nutrition component in rice is critical for improvement of people’s health. Vitamin E serves as important antioxidant by quenching the free radical intermediates and thus protects the cell membrane. Because of the high nutritional value and the benefits of vitamin E in human health, increasing the tocochromanol content of major agricultural crops has long been in the focus of breeding programs and genetic engineering approaches. The key genes involved in tocopherol biosynthesis have been elucidated in Arabidopsis and other model organisms. Quantitative trait locus (QTL) study performed in Arabidopsis suggested that some of these key genes and a few additional loci contribute to natural tocopherol variations. Identifying such genetic variations in rice, enrich our understanding of the genetic mechanisms controlling tocopherol variation, which can be directly applied to rice breeding programs. In this study, we used genome-wide association mapping with high-resolution density SNPs of rice core set to identify natural allelic variations, which contribute to tocopherol increase in rice

Chloroplast (cp) DNA sequence data are a versatile tool for plant identification, barcoding and establishing genetic relationships among plant species. Different chloroplast loci have been utilized to infer evolutionary relationship of plant species. Although the overall structure of the chloroplast genome is generally well conserved, a number of mutations have been observed. Thus, documentation of chloroplast sequence variation has also been an valuable asset in plant population and evolutionary studies for over two decades. Recently, advance in chloroplast genome assembly from whole genome NGS data has become available. In the present study, chloroplast variations among 295 diverse origin accessions were detected based on Oryza rufipogon, which thought to be the progenitor of cultivated rice.. Variation calling was carried out using the whole genome re-sequencing data of those accessions along with the five rice reference cp genomes. Phylogenetic and evolution analysis of the six references and 295 accessions were performed using the whole reference genome sequence and the variation data, respectively. Also, nucleotide polymorphisms of 295 rice accessions were validated by using previously characterized 50 ecotypes. Differential SNP frequency across the rice cp genome suggests a regional dependent preferential high variation occurrence during the evolution of chloroplast

Genetic resources play a great role in crop breeding because of containing a broad array of useful genes. Currently, the harder are rice breeders trying to develop new rice cultivars with the improved traits, they are more often handicapped by the limited availability of germplasm resources. Thus, a desirable core or heuristic (HS) set of germplasm with maximum genetic diversity can be usefully exploited to breakthrough the present and future challenges of the rice breeding. As such we previously developed the rice HS sets of 166 diverse accessions out of a total 24,368 rice germplasms. Here, we report a large-scale analysis of the patterns of genome-wide genetic variations accumulated in the HS as well as Korean rice over the time. We characterized a total of about 11.8 millions of single nucleotide polymorphisms (SNPs) across the rice genome from resequencing a total of 295 rice genomes including 137 HS and 158 KB rice sets, with an average of approximately 10x depth and > 90% coverage. Using about 460,000 high-quality SNPs (HQSNPs), we specified the population structure, confirming our HS set covers all the rice sub-populations. We further traced the relative nucleotide variabilities of HQSNPs and found the level of the diversity was dynamically changing across the KB genome, which reveals the selection history of KB lines in the past and present. In addition, the results of our genome wide association study (GWAS) suggests that our HS can be also a good reservoir of valuable alleles, pinpointing those alleles underlying the important rice agronomical traits. Overall, the resequencing of our HS set re-illuminates the past, present of the germplasm utilization, which will support the Korean rice breeding in the future.

Recently whole genome SNP genotyping has been used to do association analysis and to map a gene of interest. Here we report application of bulked segregant analysis(BSA) using Infinium HD assay with ‘BARC Bean6K_3’, a SNP genotyping beadchip containing 5,399 SNPs for common bean to locate a target gene. We used BSA using Infinium HD assay was performed to find the candidate region of a single dominant rust resistant gene in PI310762, a common bean cultivar. And SSR markers were identified and mapped on the candidate region using F2 population derived from the cross of susceptible Pinto114 x resistant PI310762. BSA revealed the candidate region of the resistant gene is on chromosome 4 where we developed nine SSR markers. Three SSR markers (beanssr1170, beanssr1168, and beanssr1167) of them appeared closely linked to the resistant gene which is located between beanssr1167 at 0.1cM and beanssr1170 at 0.5cM on chromosome 4. This study showed BSA using high-throughput whole genome SNP genotyping is a very fast and efficient method to locate a gene of interest on chromosome.

Progress in next-generation sequencing technologies have enabled discovery of massive amount of genome-wide DNA polymorphisms, single nucleotide polymorphisms (SNPs) and insertion-deletion (InDels), which are an invaluable resource to analyze genetic diversity in a population. We performed whole-genome resequencing of ten Korean rice accessions including six cultivars and four mutant lines. A total of 2,448 million raw reads was generated with 58-fold coverage and uniquely mapped to 87.5% of the Nipponbare as a reference genome. We identified 3,240,025 DNA polymorphisms including 2,867,878 SNPs, 151,845 insertions and 220,302 deletions between the Korean rice accessions and Nipponbare. We observed that in ten Korean rice accessions, the frequency of potential SNPs was estimated to be one per 2.1kb on Nipponbare (382Mb). According to annotation of DNA polymorphisms, 634,617 SNPs were found in gene region, and only 169,738 SNPs were occurred in coding region. Altogether, 86,251 non-synonymous SNPs were located on 76,891 genes. We also examined the cultivar-specific SNPs to select candidate SNPs which would have possibility of being associated with unique phenotype or agronomical trait of each cultivar. It was estimated that the portion of cultivar specific SNPs is 1~12% of the total SNPs. These DNA polymorphisms obtained from our result will provide an invaluable resource to identify molecular markers and genes associated with diverse traits of agronomical importance.

Advances in genome sequencing technologies have aided discovery of millions of genome-wide DNA polymorphisms, single nucleotide polymorphisms (SNPs) and insertion-deletion (InDels), which are an invaluable resource to analyze genetic diversity in a population. We performed whole-genome resequencing of ten Korean rice accessions including six cultivars and four mutant lines. A total of 2,447 million raw reads were generated with over 58x coverage and detected 3,240,025 DNA polymorphisms between the Korean rice accessions and Nipponbare as reference genome. We observed that in ten Korean rice accessions, the frequency of potential SNPs was estimated to be one per 2.1kb on Nipponbare (382Mb). Potential SNPs were classified into two types, homozygous SNP and heterozygous SNP, which approximately 87% of the total was homozygous SNPs from ten accessions and heterozygous SNPs accounted for 13%. According to annotation of DNA polymorphisms, 634,620 SNPs were found in gene region, and only 169,738 SNPs were occurred in coding region. Altogether, 86,251 non-synonymous SNPs were located on 76,891 genes. We also examined genes which had at least one SNP in all ten accessions. It was estimated that the total of 290 genes had one or more non-synonymous SNPs and 25 genes had only synonymous SNPs. These genes were functionally classified based on gene ontology (GO). These DNA polymorphisms obtained from our result will provide an invaluable resource to identify molecular markers and genes associated with diverse traits of agronomical importance.

The application of next generation sequencing technologies allows us to discover the high levels of DNA polymorphism throughout a genome, e.g., single nucleotide polymorphisms (SNPs), and insertions and deletions (InDels). We performed whole-genome resequencing of a Korean rice cultivar (cv. Donganbyeo) and then obtained the sequences of covered 366,042,872 bp (96.63%) with average mapped read depth of 34.17 on 382,788,128 bp of the Japanese cultivar genome (cv. Nipponbare). We characterized the polymorphisms of 173,711 SNPs, 295,334 insertions and 40,642 deletions based on the comparison of both genomes. About 11.5% and 17.8% of the annotated total SNPs were presented in the regions of 1kb upstreams and genes, respectively. The annotated InDels in gene regions were similar with 15.5% insertion (4,588) and 15.9% (5,100) deletions, but not in 1kb upstream regions with 9.0% insertion (2,662) and 14.3% deletions (5,100). In addition, the Korea rice genome sequences were mapped on individual chromosome, resulted that SNPs were shown with different frequencies from each chromosome. The InDels distributions on individual chromosomes exhibited similar pattern as compared to those of SNPs. Some gene families such as NB-ARC (NB-LRR), F-box, RLK (serine/threonine protein kinase) and Zinc-finger (RING) for SNPs occurred the similar pattern with those of Arabidopsis. These results might be useful for better understanding the genome structure and genetic diversity of the Korean rice cultivars.