Recently, we published a microinjection method for generating transgenic cattle using the DNA transposon system and their analysis by next-generation sequencing (Yum et al. Sci Rep. 2016 Jun 21;6:27185). In that study, we generated transgenic cattle using two different types of DNA transposon system, sleeping beauty (SB) and piggybac (PB), carrying Yellow fluorescent protein with SB (SB-YFP, female) and green fluorescent protein with PB (PB-GFP, male) under the control of the ubiquitous CAG promoter, respectively. The female and male founder cattle have been grown up to date (the female age: 40 months old, the male age: 33 months old) without any health issues. In genomic instability and blood analysis, there was no significant differences between wild type and founder cattle. In the present study, we confirmed germ-line transmission of the transposon-mediated transgene integrations and ubiquitous and persistent expression of transgene in second generation of offspring (F1). The F1 was born without any assistance and expressed GFP in the eyes without UV light. The ubiquitous expression of GFP was detected in skin fibroblast from the ear tissue and confirmed by genomic DNA PCR, which suggest that the transgene from the PB-GFP was successfully transmitted. Unfortunately, no transgene from SB-YFP were identified. To confirm the transgene integration site, the genomic DNA from blood was extracted and performed next-generation sequencing (NGS). The GFP gene was integrated in chromosome 4 (two copies), and 6. As results, a total of two copies of paternal transgene transmitted into the F1. All the integrated position was not related with coding region and there was no significant difference in genomic variants between transgenic and non-transgenic cattle. To our knowledge, this is the first report of germ-line transmission through non-viral transgenic founder cattle. Those transgenic cattle will be valuable resource to many fields of biomedical research and agricultural science.

Lentinula edodes, the popular shiitake mushroom, is one of the most important cultivated edible mushrooms. It is used as a food and for medicinal purposes. Here, we present the 46.1 Mb draft genome of L. edodes, comprising 13,028 predicted gene models. The genome assembly consists of 31 scaffolds. Gene annotation provides key information about various signaling pathways and secondary metabolites. This genomic information should help establish the molecular genetic markers for MAS/MAB and increase our understanding of the genome structure and function.

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.

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.

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.

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.