A bioinformatics approach for identifying transgene insertion sites using whole genome sequencing data.

BACKGROUND: Genetically modified crops (GM crops) have been developed to improve the agricultural traits of modern crop cultivars. Safety assessments of GM crops are of paramount importance in research at developmental stages and before releasing transgenic plants into the marketplace. Sequencing technology is developing rapidly, with higher output and labor efficiencies, and will eventually replace existing methods for the molecular characterization of genetically modified organisms. METHODS: To detect the transgenic insertion locations in the three GM rice gnomes, Illumina sequencing reads are mapped and classified to the rice genome and plasmid sequence. The both mapped reads are classified to characterize the junction site between plant and transgene sequence by sequence alignment. RESULTS: Herein, we present a next generation sequencing (NGS)-based molecular characterization method, using transgenic rice plants SNU-Bt9-5, SNU-Bt9-30, and SNU-Bt9-109. Specifically, using bioinformatics tools, we detected the precise insertion locations and copy numbers of transfer DNA, genetic rearrangements, and the absence of backbone sequences, which were equivalent to results obtained from Southern blot analyses. CONCLUSION: NGS methods have been suggested as an effective means of characterizing and detecting transgenic insertion locations in genomes. Our results demonstrate the use of a combination of NGS technology and bioinformatics approaches that offers cost- and time-effective methods for assessing the safety of transgenic plants.

The involvement of ß-amyrin 28-oxidase (CYP716A52v2) in oleanane-type ginsenoside biosynthesis in Panax ginseng.

Panax species are the most popular medicinal herbs. The root of these plants contains pharmacologically active triterpene saponins, also known as ginsenosides, compounds that are divided into dammarane- and oleanane-type triterpenes. Two CYP716A subfamily genes (CYP716A47 and CYP716A53v2) were recently characterized, encoding an enzyme catalyzing the hydroxylation of dammarane-type triterpenes in Panax ginseng. Herein, we report that one CYP716A subfamily gene (CYP716A52v2) isolated from P. ginseng encodes a ß-amyrin 28-oxidase, which is suggested to modify ß-amyrin into oleanolic acid, a precursor of an oleanane-type saponin (mainly ginsenoside Ro) in P. ginseng. The ectopic expression of both PNY1 and CYP716A52v2 in recombinant yeast resulted in erythrodiol and oleanolic acid production, respectively. In vitro enzymatic activity assays biochemically confirmed that CYP716A52v2 catalyzed the oxidation of ß-amyrin to produce oleanolic acid, and the chemical structure of the oleanolic acid product was confirmed using gas chromatography-mass spectrometry (GC/MS). Transgenic P. ginseng plants were generated via Agrobacterium tumefaciens-mediated transformation: the overexpression of CYP716A52v2 greatly increased the content of oleanane-type ginsenoside (ginsenoside Ro), whereas RNA interference against CYP716A52v2 markedly reduced it. Furthermore, the levels of other dammarene-type ginsenosides were not affected in these transgenic lines. These results indicate that CYP716A52v2 is a ß-amyrin 28-oxidase that plays a key role in the biosynthesis of oleanane-type triterpenes in P. ginseng.

Analysis of genome variants in dwarf soybean lines obtained in F6 derived from cross of normal parents (cultivated and wild soybean).

Plant height is an important component of plant architecture and significantly affects crop breeding practices and yield. We studied DNA variations derived from F5 recombinant inbred lines (RILs) with 96.8% homozygous genotypes. Here, we report DNA variations between the normal and dwarf members of four lines harvested from a single seed parent in an F6 RIL population derived from a cross between Glycine max var. Peking and Glycine soja IT182936. Whole genome sequencing was carried out, and the DNA variations in the whole genome were compared between the normal and dwarf samples. We found a large number of DNA variations in both the dwarf and semi-dwarf lines, with one single nucleotide polymorphism (SNP) per at least 3.68 kb in the dwarf lines and 1 SNP per 11.13 kb of the whole genome. This value is 2.18 times higher than the expected DNA variation in the F6 population. A total of 186 SNPs and 241 SNPs were discovered in the coding regions of the dwarf lines 1282 and 1303, respectively, and we discovered 33 homogeneous nonsynonymous SNPs that occurred at the same loci in each set of dwarf and normal soybean. Of them, five SNPs were in the same positions between lines 1282 and 1303. Our results provide important information for improving our understanding of the genetics of soybean plant height and crop breeding. These polymorphisms could be useful genetic resources for plant breeders, geneticists, and biologists for future molecular biology and breeding projects.

Comparative transcriptome analysis reveals higher expression of stress and defense responsive genes in dwarf soybeans obtained from the crossing of G. max and G. soja.

BACKGROUND: Plant height is an important component of plant architecture and significantly affects crop breeding practices and yield. Dwarfism in plants prevents lodging and therefore it's a desired trait in crops. OBJECTIVE: To find differentially expressed genes to classify and understand the regulation of genes related to plant growth in mutant dwarf soybeans, which appeared in the F5 generation. METHODS: We obtained a few segregated dwarf soybeans in the populations derived from the crossing of Glycine max var. Peking and Glycine soja var. IT182936 in an F5 RIL population. These dwarf soybeans may be useful genetic resources for plant breeders, geneticists and biologists. Using the Illumina high-throughput platform, transcriptomes were generated and compared among normal and dwarf soybeans in triplicate. CONCLUSION: We found complex relationship of the expressed genes to plant growth. There were highly significantly up-/downregulated genes according to the comparison of gene expression in normal and dwarf soybeans. The genes related to disease and stress responses were found to be upregulated in dwarf soybeans. Such over-expression of disease resistance and other immune response genes can be targeted to understand how the immune genes regulate the response of plant growth. In addition, photosynthesis-related genes showed very low expression in dwarf lines. The transcriptome expression and genes classified as related to plant growth may be useful resources to researchers studying plant growth.

RNA-Seq De Novo Assembly and Differential Transcriptome Analysis of Korean Medicinal Herb Cirsium japonicum var. spinossimum.

Cirsium japonicum belongs to the Asteraceae or Compositae family and is a medicinal plant in Asia that has a variety of effects, including tumour inhibition, improved immunity with flavones, and antidiabetic and hepatoprotective effects. Silymarin is synthesized by 4-coumaroyl-CoA via both the flavonoid and phenylpropanoid pathways to produce the immediate precursors taxifolin and coniferyl alcohol. Then, the oxidative radicalization of taxifolin and coniferyl alcohol produces silymarin. We identified the expression of genes related to the synthesis of silymarin in C. japonicum in three different tissues, namely, flowers, leaves and roots, through RNA sequencing. We obtained 51,133 unigenes from transcriptome sequencing by de novo assembly using Trinity v2.1.1, TransDecoder v2.0.1, and CD-HIT v4.6 software. The differentially expressed gene analysis revealed that the expression of genes related to the flavonoid pathway was higher in the flowers, whereas the phenylpropanoid pathway was more highly expressed in the roots. In this study, we established a global transcriptome dataset for C. japonicum. The data shall not only be useful to focus more deeply on the genes related to product medicinal metabolite including flavolignan but also to study the functional genomics for genetic engineering of C. japonicum.

Alteration of Panax ginseng saponin composition by overexpression and RNA interference of the protopanaxadiol 6-hydroxylase gene (CYP716A53v2).

BACKGROUND: The roots of Panax ginseng contain noble tetracyclic triterpenoid saponins derived from dammarenediol-II. Dammarene-type ginsenosides are classified into the protopanaxadiol (PPD) and protopanaxatriol (PPT) groups based on their triterpene aglycone structures. Two cytochrome P450 (CYP) genes (CYP716A47 and CYP716A53v2) are critical for the production of PPD and PPT aglycones, respectively. CYP716A53v2 is a protopanaxadiol 6-hydroxylase that catalyzes PPT production from PPD in P. ginseng. METHODS: We constructed transgenic P. ginseng lines overexpressing or silencing (via RNA interference) the CYP716A53v2 gene and analyzed changes in their ginsenoside profiles. RESULT: Overexpression of CYP716A53v2 led to increased accumulation of CYP716A53v2 mRNA in all transgenic roots compared to nontransgenic roots. Conversely, silencing of CYP716A53v2 mRNA in RNAi transgenic roots resulted in reduced CYP716A53v2 transcription. HPLC analysis revealed that transgenic roots overexpressing CYP716A53v2 contained higher levels of PPT-group ginsenosides (Rg1, Re, and Rf) but lower levels of PPD-group ginsenosides (Rb1, Rc, Rb2, and Rd). By contrast, RNAi transgenic roots contained lower levels of PPT-group compounds and higher levels of PPD-group compounds. CONCLUSION: The production of PPD- and PPT-group ginsenosides can be altered by changing the expression of CYP716A53v2 in transgenic P. ginseng. The biological activities of PPD-group ginsenosides are known to differ from those of the PPT group. Thus, increasing or decreasing the levels of PPT-group ginsenosides in transgenic P. ginseng may yield new medicinal uses for transgenic P. ginseng.