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.

Genome-wide characterization of long intergenic non-coding RNAs (lincRNAs) provides new insight into viral diseases in honey bees Apis cerana and Apis mellifera.

BACKGROUND: Long non-coding RNAs (lncRNAs) are a class of RNAs that do not encode proteins. Recently, lncRNAs have gained special attention for their roles in various biological process and diseases. RESULTS: In an attempt to identify long intergenic non-coding RNAs (lincRNAs) and their possible involvement in honey bee development and diseases, we analyzed RNA-seq datasets generated from Asian honey bee (Apis cerana) and western honey bee (Apis mellifera). We identified 2470 lincRNAs with an average length of 1011 bp from A. cerana and 1514 lincRNAs with an average length of 790 bp in A. mellifera. Comparative analysis revealed that 5 % of the total lincRNAs derived from both species are unique in each species. Our comparative digital gene expression analysis revealed a high degree of tissue-specific expression among the seven major tissues of honey bee, different from mRNA expression patterns. A total of 863 (57 %) and 464 (18 %) lincRNAs showed tissue-dependent expression in A. mellifera and A. cerana, respectively, most preferentially in ovary and fat body tissues. Importantly, we identified 11 lincRNAs that are specifically regulated upon viral infection in honey bees, and 10 of them appear to play roles during infection with various viruses. CONCLUSIONS: This study provides the first comprehensive set of lincRNAs for honey bees and opens the door to discover lincRNAs associated with biological and hormone signaling pathways as well as various diseases of honey bee.

Uncovering the novel characteristics of Asian honey bee, Apis cerana, by whole genome sequencing.

BACKGROUND: The honey bee is an important model system for increasing understanding of molecular and neural mechanisms underlying social behaviors relevant to the agricultural industry and basic science. The western honey bee, Apis mellifera, has served as a model species, and its genome sequence has been published. In contrast, the genome of the Asian honey bee, Apis cerana, has not yet been sequenced. A. cerana has been raised in Asian countries for thousands of years and has brought considerable economic benefits to the apicultural industry. A cerana has divergent biological traits compared to A. mellifera and it has played a key role in maintaining biodiversity in eastern and southern Asia. Here we report the first whole genome sequence of A. cerana. RESULTS: Using de novo assembly methods, we produced a 238 Mbp draft of the A. cerana genome and generated 10,651 genes. A.cerana-specific genes were analyzed to better understand the novel characteristics of this honey bee species. Seventy-two percent of the A. cerana-specific genes had more than one GO term, and 1,696 enzymes were categorized into 125 pathways. Genes involved in chemoreception and immunity were carefully identified and compared to those from other sequenced insect models. These included 10 gustatory receptors, 119 odorant receptors, 10 ionotropic receptors, and 160 immune-related genes. CONCLUSIONS: This first report of the whole genome sequence of A. cerana provides resources for comparative sociogenomics, especially in the field of social insect communication. These important tools will contribute to a better understanding of the complex behaviors and natural biology of the Asian honey bee and to anticipate its future evolutionary trajectory.

Complex Extract of Polygonatum sibiricum and Nelumbinis semen Improves Menopause Symptoms via Regulation of Estrogen Receptor Beta in an Ovariectomized Rat Model.

Menopause is a hormone-deficiency state that causes facial flushing, vaginal dryness, depression, anxiety, insomnia, obesity, osteoporosis, and cardiovascular disease as ovarian function decreases. Hormone-replacement therapy is mainly used to treat menopause; however, its long-term use is accompanied by side effects such as breast cancer and endometriosis. To identify the effect of a complex extract of Polygonatum sibiricum (PS) and Nelumbinis semen (NS) on improving menopause without side effects, an ovariectomized rat model was established to analyze several menopause symptoms. Compared to single extracts, the complex extract restored vaginal epithelial cell thickness and decreased serotonin concentration by increasing the estrogen receptors ERα (ESR1) and ERß (ESR2), depending on the ratio. Although the complex extract exerted a lower weight-loss effect than the single extracts, improved blood-lipid metabolism was observed after increasing high-density lipoprotein cholesterol levels and decreasing low-density lipoprotein cholesterol and triglyceride levels, and ovariectomy-induced osteoporosis was alleviated by suppressing osteoclast production. Thus, by increasing only ERß expression without regulating ERα expression in the uterus, the complex extract of PS and NS may be a natural treatment for improving menopause symptoms without side effects, such as endometriosis.

Detection of mPing mobilization in transgenic rice plants.

BACKGROUND: Various kinds of transposable elements (TEs) constitute high proportions of eukaryotic genomes. Although most of these TEs are not actively mobile, genome stress can induce mobilization of dormant TEs. Transgenic plants undergo tissue culture and subsequent whole-plant regeneration, which can cause genomic stress and in turn induce mobilization of inactive TEs. OBJECTIVES: To investigate the activation of transposable elements on the genome wide of the GM plant. METHODS: Transposon activities were analyzed in three transgenic rice plants carrying the insect resistance gene Cry1Ac and an herbicide resistance gene by the transposon display technique. These three transgenic plants were derived from a leading Korean rice variety, Illmi. RESULTS: We detected seven mobile activities in the mPing element, which is a MITE family transposon. The identity of the novel fragments in the gel display was confirmed by checking TAA target site duplication via sequence analysis. The genomic integration sites were all on different chromosomes, and the integrations were specific to either one or two T1 transgenic lines, except for one common integration on chromosome 4. One integration was in the 5'-UTR of the Glycerol-3-phosphate acyltransferase 8 gene, two integrations were in introns of expressed genes, and the other four integrations were in intergenic regions. CONCLUSION: Thus, novel mobilization of dormant TEs occurs in transgenic plants, which must be considered in the generation of genetically modified crops (GM crops).

De novo transcriptome sequencing and gene expression profiling with/without B-chromosome plants of Lilium amabile.

Supernumerary B chromosomes were found in Lilium amabile (2n = 2x = 24), an endemic Korean lily that grows in the wild throughout the Korean Peninsula. The extra B chromosomes do not affect the host-plant morphology; therefore, whole transcriptome analysis was performed in 0B and 1B plants to identify differentially expressed genes. A total of 154,810 transcripts were obtained from over 10 Gbp data by de novo assembly. By mapping the raw reads to the de novo transcripts, we identified 7,852 differentially expressed genes (log2FC > |10|), in which 4,059 and 3,794 were up- and down-regulated, respectively, in 1B plants compared to 0B plants. Functional enrichment analysis revealed that various differentially expressed genes were involved in cellular processes including the cell cycle, chromosome breakage and repair, and microtubule formation; all of which may be related to the occurrence and maintenance of B chromosomes. Our data provide insight into transcriptomic changes and evolution of plant B chromosomes and deliver an informative database for future study of B chromosome transcriptomes in the Korean lily.

NGS sequencing reveals that many of the genetic variations in transgenic rice plants match the variations found in natural rice population.

BACKGROUND: As the transformation process can induce mutations in host plants, molecular characterization of the associated genomic changes is important not only for practical food safety but also for understanding the fundamental theories of genome evolution. OBJECTIVES: To investigate a population-scale comparative study of the genome-wide spectrum of sequence variants in the transgenic genome with the variations present in 3000 rice varieties. RESULTS: On average, we identified 19,273 SNPs (including Indels) per transgenic line in which 10,729 SNPs were at the identical locations in the three transgenic rice plants. We found that these variations were predominantly present in specific regions in chromosomes 8 and 10. Majority (88%) of the identified variations were detected at the same genomic locations as those in natural rice population, implying that the transgenic induced mutations had a tendency to be common alleles. CONCLUSION: Genomic variations in transgenic rice plants frequently occurred at the same sites as the major alleles found in the natural rice population, which implies that the sequence variations occur within the limits of a biological system to ensure survival.

Characterization of chloroplast genomes of Alnus rubra and Betula cordifolia, and their use in phylogenetic analyses in Betulaceae.

OBJECTIVE: Betulaceae is a relatively small birch family that comprises about 160 deciduous trees and shrubs. Chloroplast (cp) genome sequencing of Alnus rubra and Betula cordifolia was carried out to elucidate their molecular features and phylogenetic relationship among species in Betulaceae family. METHODS: Chloroplast genome sequencing was carried out using next generation sequencing method. Molecular and genomic features of the two cp genomes were characterized with other cp genomes in Betulaceae. Also, molecular phylogenetic analysis was performed using the whole cp genome sequences. RESULTS: The average cp genome length was 160,136 bp among the Betulaceae species. Base compositions of the cp genomes were skewed toward a high AT ratio, with an average of 63.4%. We identified 117 different genes 83 with protein coding, 4 with ribosomal RNA, and 30 with tRNA. Eighteen genes contained introns which were conserved among the cp genomes of all Betulaceae. We mined 82 SSRs from the cp genomes of A. rubra, A. cordifolia, and A. nana. The SSRs were variable in motif repeat numbers and presence/absence among the cp genomes. CONCLUSION: Chloroplast genome-wide sequence comparison from 11 Betulaceae species and one cp genome of evergreen oak revealed that the patterns of sequence variations were congruent with two subfamily classification Betuloideae (Alnus and Betula) and Corylaceae (Corylus, Ostrya, and Carpinus). Subsequent phylogenetic analysis also supports the sub-classifications of these species.

Karyotype and B chromosome variation in Lilium amabile Palibin.

OBJECTIVE: Lilium amabile Palibin (2n = 2x = 24) is an endemic lily species in Korea. B chromosomes are supernumerary chromosomes and the presence of B chromosome in L. amabile was known by previous researches. The current research was conducted to characterize the genetical and cytological features of the B chromosome plants in L. amabile. METHODS: Karyotype and B chromosome cytotype analyses were carried out among 135 L. amabile accessions that were collected from six geographical locations in Korea using conventional aceto-carmine staining as well as FISH technique with ribosomal RNA gene probes. RESULTS: The karyotype of L. amabile genome consisted of two large metacentric, four intermediate subtelocentric, and six intermediate to small acrocentric chromosomes in which chromosomes 1, 6 and 7 carried the 45S rRNA gene loci and chromosome 3 carried the 5S rRNA gene. There were 4 types of B chromosomes, two large B chromosomes and two small B chromosomes. The ribosomal RNA gene loci were not present in the B chromosomes. The 135 accessions were classified into 13 cytotypes including diploids and different B chromosome aneuploids. Among the aneuploids, the most frequent cytotype was 24 + 1B, which was followed by 24 + 2B, 24 + 1b, 24 + 1B + 2b, 24 + 1B + 4b, and 24 + 2B + 4b. CONCLUSION: The karyotype of L. amabile was consistent with other species in the genus Lilium without polyploids. The B chromosome cytotypes were highly variable and the occurrences of different cytotypes were random among the six populations, implying that the B and b chromosome occurrence was random in each population.

A Comparative Study of Phenolic Antioxidant Activity and Flavonoid Biosynthesis-Related Gene Expression Between Summer and Winter Strawberry Cultivars.

Strawberry (Fragaria ananassa Duch.) possesses good antioxidant properties. Phenolic compounds in strawberries, such as anthocyanins and ellagic acid, mainly act as antioxidants. This study aimed to compare the phenolic content and expression patterns of genes involved in flavonoid biosynthesis between summer and winter strawberry cultivars affected by seasonal variation, degree of ripeness, and genotype. Antioxidant activity and the total content of phenols and flavonoids decreased with fruit ripening. Most notably, summer strawberry cultivars showed higher antioxidant activity than winter cultivars. The expression patterns of flavonoid biosynthetic genes tested were cultivar-dependent and were also affected by ripening. These results help us understand the nutritional and physiological characteristics of selected cultivars and provide a range of information for strawberry consumption.

A novel small-molecule PPI inhibitor targeting integrin αvß3-osteopontin interface blocks bone resorption in vitro and prevents bone loss in mice.

Small molecule-inhibition targeting protein-protein interaction (PPI) is now recognized as an emerging and challenging area in drug design. We developed a novel interactive drug discovery methodology known as Protein Chip technology (ProteoChip) as a cutting-edge PPI assay system applicable for unique PPI-targeting therapeutics integrated with computer-aided drug design (CADD). Here, we describe a novel small molecular PPI inhibitor, IPS-02001, which the blocks integrin αvß3-osteopontin interface a novel PPI inhibitor identified by the interactive methodology of both ProteoChip- and CADD-based PPI assay. IPS-02001 (6,7-Dichloro-2,3,5,8-tetrahydroxy-1,4-naphthoquinone) was screened from different compound libraries (InterBioScreen, Commercial libraries) using an in silico structure-based molecular docking simulation method and a protein chip-based protein-protein interaction assay system. Additionally, integrin αvß3, an adhesion receptor expressed in osteoclasts (OCs), was implicated in the regulation of OC function via regulation of the cytoskeletal organization of OCs. IPS-02001 blocked OC maturation from murine bone marrow-derived macrophages, as well as the resorptive function of OCs. Moreover, treatment with IPS-02001 impaired downstream signaling of integrin αvß3 linked to Pyk2, c-Src, PLCγ2, and Vav3 and disrupted the actin cytoskeleton in mature OCs. Furthermore, IPS-02001 blocked RANKL-induced bone destruction by reducing the number of OCs and protected against ovariectomy-induced bone loss in mice. Thus, IPS-02001 may represent a promising new class of anti-resorptive drugs for treatment of bone diseases associated with increased OC function.

Efficiency to Discovery Transgenic Loci in GM Rice Using Next Generation Sequencing Whole Genome Re-sequencing.

Molecular characterization technology in genetically modified organisms, in addition to how transgenic biotechnologies are developed now require full transparency to assess the risk to living modified and non-modified organisms. Next generation sequencing (NGS) methodology is suggested as an effective means in genome characterization and detection of transgenic insertion locations. In the present study, we applied NGS to insert transgenic loci, specifically the epidermal growth factor (EGF) in genetically modified rice cells. A total of 29.3 Gb (~72× coverage) was sequenced with a 2 × 150 bp paired end method by Illumina HiSeq2500, which was consecutively mapped to the rice genome and T-vector sequence. The compatible pairs of reads were successfully mapped to 10 loci on the rice chromosome and vector sequences were validated to the insertion location by polymerase chain reaction (PCR) amplification. The EGF transgenic site was confirmed only on chromosome 4 by PCR. Results of this study demonstrated the success of NGS data to characterize the rice genome. Bioinformatics analyses must be developed in association with NGS data to identify highly accurate transgenic sites.

Functional characterization of naturally occurring melittin peptide isoforms in two honey bee species, Apis mellifera and Apis cerana.

Insect-derived antimicrobial peptides (AMPs) have diverse effects on antimicrobial properties and pharmacological activities such as anti-inflammation and anticancer properties. Naturally occurring genetic polymorphism have a direct and/or indirect influence on pharmacological effect of AMPs, therefore information on single nucleotide polymorphism (SNP) occurring in natural AMPs provides an important clue to therapeutic applications. Here we identified nucleotide polymorphisms in melittin gene of honey bee populations, which is one of the potent AMP in bee venoms. We found that the novel SNP of melittin gene exists in these two honey bee species, Apis mellifera and Apis cerana. Nine polymorphisms were identified within the coding region of the melittin gene, of which one polymorphism that resulted in serine (Ser) to asparagine (Asp) substitution that can potentially effect on biological activities of melittin peptide. Serine-substituted melittin (Mel-S) showed more cytotoxic effect than asparagine-substituted melittin (Mel-N) against E. coli. Also, Mel-N and Mel-S had different inhibitory effects on the production of inflammatory factors such as IL-6 and TNF-α in BV-2 cells. Moreover, Mel-S showed stronger cytotoxic activities than Mel-N peptide against two human ovarian cancer cell lines. Using carbon nanotube-based transistor, we here characterized that Mel-S interacted with small unilamellar liposomes more strongly than Mel-N. Taken together, our present study demonstrates that there exist different characteristics of the gene frequency and the biological activities of the melittin peptide in two honey bee species, Apis mellifera and A. cerana.