Identifying Terpenoid Biosynthesis Genes in Euphorbia maculata via Full-Length cDNA Sequencing.

The annual herb Euphorbia maculata L. produces anti-inflammatory and biologically active substances such as triterpenoids, tannins, and polyphenols, and it is used in traditional Chinese medicine. Of these bioactive compounds, terpenoids, also called isoprenoids, are major secondary metabolites in E. maculata. Full-length cDNA sequencing was carried out to characterize the transcripts of terpenoid biosynthesis reference genes and determine the copy numbers of their isoforms using PacBio SMRT sequencing technology. The Illumina short-read sequencing platform was also employed to identify differentially expressed genes (DEGs) in the secondary metabolite pathways from leaves, roots, and stems. PacBio generated 62 million polymerase reads, resulting in 81,433 high-quality reads. From these high-quality reads, we reconstructed a genome of 20,722 genes, in which 20,246 genes (97.8%) did not have paralogs. About 33% of the identified genes had two or more isoforms. DEG analysis revealed that the expression level differed among gene paralogs in the leaf, stem, and root. Whole sets of paralogs and isoforms were identified in the mevalonic acid (MVA), methylerythritol phosphate (MEP), and terpenoid biosynthesis pathways in the E. maculata L. The nucleotide information will be useful for identifying orthologous genes in other terpenoid-producing medicinal plants.

Genome-wide transcriptome profiling of the medicinal plant Zanthoxylum planispinum using a single-molecule direct RNA sequencing approach.

High-throughput RNA sequencing has revolutionized transcriptome-based studies of candidate genes, key pathways and gene regulation in non-model organisms. We analyzed full-length cDNA sequences in Zanthoxylum planispinum (Z. planispinum), a medicinal herb in major parts of East Asia. The full-length mRNA derived from tissues of leaf, early fruit and maturing fruit stage were sequenced using PacBio RSII platform to identify isoform transcriptome. We obtained 51,402 unigenes, with average 1781 bp per gene in 82.473 Mb gene lengths. Among 51,402, 3963 unigenes showed variety of isoform. By selection of one representative gene among each of the various isoforms, we finalized 46,306 unique gene set for this herb. We identified 76 cytochrome P450 (CYP450) and related isoforms that are of the wide diversity in the molecular function and biological process. These transcriptome data of Z. planispinum will provide a good resource to study metabolic engineering for the production of valuable medicinal drugs and phytochemicals.

Ginseng Genome Database: an open-access platform for genomics of Panax ginseng.

BACKGROUND: The ginseng (Panax ginseng C.A. Meyer) is a perennial herbaceous plant that has been used in traditional oriental medicine for thousands of years. Ginsenosides, which have significant pharmacological effects on human health, are the foremost bioactive constituents in this plant. Having realized the importance of this plant to humans, an integrated omics resource becomes indispensable to facilitate genomic research, molecular breeding and pharmacological study of this herb. DESCRIPTION: The first draft genome sequences of P. ginseng cultivar "Chunpoong" were reported recently. Here, using the draft genome, transcriptome, and functional annotation datasets of P. ginseng, we have constructed the Ginseng Genome Database http://ginsengdb.snu.ac.kr /, the first open-access platform to provide comprehensive genomic resources of P. ginseng. The current version of this database provides the most up-to-date draft genome sequence (of approximately 3000 Mbp of scaffold sequences) along with the structural and functional annotations for 59,352 genes and digital expression of genes based on transcriptome data from different tissues, growth stages and treatments. In addition, tools for visualization and the genomic data from various analyses are provided. All data in the database were manually curated and integrated within a user-friendly query page. CONCLUSION: This database provides valuable resources for a range of research fields related to P. ginseng and other species belonging to the Apiales order as well as for plant research communities in general. Ginseng genome database can be accessed at http://ginsengdb.snu.ac.kr /.

Genome and evolution of the shade-requiring medicinal herb Panax ginseng.

Panax ginseng C. A. Meyer, reputed as the king of medicinal herbs, has slow growth, long generation time, low seed production and complicated genome structure that hamper its study. Here, we unveil the genomic architecture of tetraploid P. ginseng by de novo genome assembly, representing 2.98 Gbp with 59 352 annotated genes. Resequencing data indicated that diploid Panax species diverged in association with global warming in Southern Asia, and two North American species evolved via two intercontinental migrations. Two whole genome duplications (WGD) occurred in the family Araliaceae (including Panax) after divergence with the Apiaceae, the more recent one contributing to the ability of P. ginseng to overwinter, enabling it to spread broadly through the Northern Hemisphere. Functional and evolutionary analyses suggest that production of pharmacologically important dammarane-type ginsenosides originated in Panax and are produced largely in shoot tissues and transported to roots; that newly evolved P. ginseng fatty acid desaturases increase freezing tolerance; and that unprecedented retention of chlorophyll a/b binding protein genes enables efficient photosynthesis under low light. A genome-scale metabolic network provides a holistic view of Panax ginsenoside biosynthesis. This study provides valuable resources for improving medicinal values of ginseng either through genomics-assisted breeding or metabolic engineering.

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.

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.

Precision genome engineering with programmable DNA-nicking enzymes.

Zinc finger nucleases (ZFNs) are powerful tools of genome engineering but are limited by their inevitable reliance on error-prone nonhomologous end-joining (NHEJ) repair of DNA double-strand breaks (DSBs), which gives rise to randomly generated, unwanted small insertions or deletions (indels) at both on-target and off-target sites. Here, we present programmable DNA-nicking enzymes (nickases) that produce single-strand breaks (SSBs) or nicks, instead of DSBs, which are repaired by error-free homologous recombination (HR) rather than mutagenic NHEJ. Unlike their corresponding nucleases, zinc finger nickases allow site-specific genome modifications only at the on-target site, without the induction of unwanted indels. We propose that programmable nickases will be of broad utility in research, medicine, and biotechnology, enabling precision genome engineering in any cell or organism.

O- and N-Methyltransferases in benzylisoquinoline alkaloid producing plants.

BACKGROUND: Secondary metabolites such as benzylisoquinoline alkaloids (BIA) have attracted considerable attention because of their pharmacological properties and potential therapeutic applications. Methyltransferases (MTs) can add methyl groups to alkaloid molecules, altering their physicochemical properties and bioactivity, stability, solubility, and recognition by other cellular components. Five types of O-methyltransferases and two types of N-methyltransferases are involved in BIA biosynthesis. OBJECTIVE: Since MTs may be the source for the discovery and development of novel biomedical, agricultural, and industrial compounds, we performed extensive molecular and phylogenetic analyses of O- and N-methyltransferases in BIA-producing plants. METHODS: MTs involved in BIA biosynthesis were isolated from transcriptomes of Berberis koreana and Caulophyllum robustum. We also mined the methyltransferases of Coptis japonica, Papaver somniferum, and Nelumbo nucifera from the National Center for Biotechnology Information protein database. Then, we analyzed the functional motifs and phylogenetic analysis. RESULT: We mined 42 O-methyltransferases and 8 N-methyltransferases from the five BIA-producing plants. Functional motifs for S-adenosyl-L-methionine-dependent methyltransferases were retained in most methyltransferases, except for the three O-methyltransferases from N. nucifera. Phylogenetic analysis revealed that the methyltransferases were grouped into four clades, I, II, III and IV. The clustering patterns in the phylogenetic analysis suggested a monophyletic origin of methyltransferases and gene duplication within species. The coexistence of different O-methyltransferases in the deep branch subclade might support some cases of substrate promiscuity. CONCLUSIONS: Methyltransferases may be a source for the discovery and development of novel biomedical, agricultural, and industrial compounds. Our results contribute to further understanding of their structure and reaction mechanisms, which will require future functional studies.

Characterization of a new high copy Stowaway family MITE, BRAMI-1 in Brassica genome.

BACKGROUND: Miniature inverted-repeat transposable elements (MITEs) are expected to play important roles in evolution of genes and genome in plants, especially in the highly duplicated plant genomes. Various MITE families and their roles in plants have been characterized. However, there have been fewer studies of MITE families and their potential roles in evolution of the recently triplicated Brassica genome. RESULTS: We identified a new MITE family, BRAMI-1, belonging to the Stowaway super-family in the Brassica genome. In silico mapping revealed that 697 members are dispersed throughout the euchromatic regions of the B. rapa pseudo-chromosomes. Among them, 548 members (78.6%) are located in gene-rich regions, less than 3 kb from genes. In addition, we identified 516 and 15 members in the 470 Mb and 15 Mb genomic shotgun sequences currently available for B. oleracea and B. napus, respectively. The resulting estimated copy numbers for the entire genomes were 1440, 1464 and 2490 in B. rapa, B. oleracea and B. napus, respectively. Concurrently, only 70 members of the related Arabidopsis ATTIRTA-1 MITE family were identified in the Arabidopsis genome. Phylogenetic analysis revealed that BRAMI-1 elements proliferated in the Brassica genus after divergence from the Arabidopsis lineage. MITE insertion polymorphism (MIP) was inspected for 50 BRAMI-1 members, revealing high levels of insertion polymorphism between and within species of Brassica that clarify BRAMI-1 activation periods up to the present. Comparative analysis of the 71 genes harbouring the BRAMI-1 elements with their non-insertion paralogs (NIPs) showed that the BRAMI-1 insertions mainly reside in non-coding sequences and that the expression levels of genes with the elements differ from those of their NIPs. CONCLUSION: A Stowaway family MITE, named as BRAMI-1, was gradually amplified and remained present in over than 1400 copies in each of three Brassica species. Overall, 78% of the members were identified in gene-rich regions, and it is assumed that they may contribute to the evolution of duplicated genes in the highly duplicated Brassica genome. The resulting MIPs can serve as a good source of DNA markers for Brassica crops because the insertions are highly dispersed in the gene-rich euchromatin region and are polymorphic between or within species.

Complete mitochondrial genome sequence and identification of a candidate gene responsible for cytoplasmic male sterility in radish (Raphanus sativus L.) containing DCGMS cytoplasm.

A novel cytoplasmic male sterility (CMS) conferred by Dongbu cytoplasmic and genic male-sterility (DCGMS) cytoplasm and its restorer-of-fertility gene (Rfd1) was previously reported in radish (Raphanus sativus L.). Its inheritance of fertility restoration and profiles of mitochondrial DNA (mtDNA)-based molecular markers were reported to be different from those of Ogura CMS, the first reported CMS in radish. The complete mitochondrial genome sequence (239,186 bp; GenBank accession No. KC193578) of DCGMS mitotype is reported in this study. Thirty-four protein-coding genes and three ribosomal RNA genes were identified. Comparative analysis of a mitochondrial genome sequence of DCGMS and previously reported complete sequences of normal and Ogura CMS mitotypes revealed various recombined structures of seventeen syntenic sequence blocks. Short-repeat sequences were identified in almost all junctions between syntenic sequence blocks. Phylogenetic analysis of three radish mitotypes showed that DCGMS was more closely related to the normal mitotype than to the Ogura mitotype. A single 1,551-bp unique region was identified in DCGMS mtDNA sequences and a novel chimeric gene, designated orf463, consisting of 128-bp partial sequences of cox1 gene and 1,261-bp unidentified sequences were found in the unique region. No other genes with a chimeric structure, a major feature of most characterized CMS-associated genes in other plant species, were found in rearranged junctions of syntenic sequence blocks. Like other known CMS-associated mitochondrial genes, the predicted gene product of orf463 contained 12 transmembrane domains. Thus, this gene product might be integrated into the mitochondrial membrane. In total, the results indicate that orf463 is likely to be a casual factor for CMS induction in radish containing the DCGMS cytoplasm.

Whole spectrum of cytochrome P450 genes and molecular responses to water-accommodated fractions exposure in the marine medaka.

Water-accommodated fractions (WAFs) of crude oil include chemicals that are potent toxicants in fish. Increasing oil pollution thus demands a better understanding of molecular mechanisms for detoxification, metabolism, toxicity, and adaptation of fish. Previous studies with fish show modulation of expression of key genes in relation to stress response against WAF exposure, but there is still a lack of studies on responses of cytochrome P450 (CYP) genes and changes in biotransformation upon WAF exposure. In this study, we used the full spectrum of CYP genes of the marine medaka, Oryzias melastigma, to understand their potential mode of action on WAFs-triggered molecular mechanisms. We also analyzed further CYP-involved detoxification and endogenous steroidogenic metabolism after exposure to different concentrations of WAFs over different time courses in the marine medaka. Also, detoxification- and antioxidant-related enzymes' activities were analyzed with different concentrations of WAFs. As a result, the WAF exposure induced CYP-involved detoxification mechanism but reduced CYP-involved steroidogenic metabolism in the marine medaka. These data suggest that whole CYP profiling can be a way of understanding and uncovering the mode of action particularly with respect to emerging chemicals such as WAF exposure with the new finding that WAFs have dual functions on CYP-involved metabolisms.

Isolation and Characterization of the Genes Involved in the Berberine Synthesis Pathway in Asian Blue Cohosh, Caulophyllum robustum.

Caulophyllum robustum, commonly named Asian blue cohosh, is a perennial herb in the family Berberidaceae. It has traditionally been used for folk medicine in China. We isolated berberine from the leaves, stem, roots, and fruits of C. robustum, and this is the first report on berberine in this species. Transcriptome analysis was conducted for the characterization of berberine biosynthesis genes in C. robustum, in which, all the genes for berberine biosynthesis were identified. From 40,094 transcripts, using gene ontology (GO) analysis, 26,750 transcripts were assigned their functions in the categories of biological process, molecular function, and cellular component. In the analysis of genes expressed in different tissues, the numbers of genes in the categories of intrinsic component of membrane and transferase activity were up-regulated in leaves versus stem. The berberine synthesis genes in C. robustum were characterized by phylogenetic analysis with corresponding genes from other berberine-producing species. The co-existence of genes from different plant families in the deepest branch subclade implies that the differentiation of berberine synthesis genes occurred early in the evolution of berberine-producing plants. Furthermore, the copy number increment of the berberine synthesis genes was detected at the species level.

Genome sequence of Pectobacterium carotovorum subsp. carotovorum strain PCC21, a pathogen causing soft rot in Chinese cabbage.

Pectobacterium carotovorum is a plant-pathogenic enterobacterium responsible for soft rot in various commercially important plants. Here we report the complete genome sequence and automatic annotation of strain PCC21.

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 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.

Complete genome sequence of Mycobacterium intracellulare strain ATCC 13950(T).

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-64, belonging to the INT1 genotype.

Here, we report the complete genome sequence of the Mycobacterium intracellulare clinical strain MOTT-64, previously grouped into the INT1 genotype among five genotypes of M. intracellulare. This genome sequence will serve as a valuable reference for understanding the disparity in the virulence and epidemiologic traits among M. intracellulare genotypes.

Comparative mapping, genomic structure, and expression analysis of eight pseudo-response regulator genes in Brassica rapa.

Circadian clocks regulate plant growth and development in response to environmental factors. In this function, clocks influence the adaptation of species to changes in location or climate. Circadian-clock genes have been subject of intense study in models such as Arabidopsis thaliana but the results may not necessarily reflect clock functions in species with polyploid genomes, such as Brassica species, that include multiple copies of clock-related genes. The triplicate genome of Brassica rapa retains high sequence-level co-linearity with Arabidopsis genomes. In B. rapa we had previously identified five orthologs of the five known Arabidopsis pseudo-response regulator (PRR) genes that are key regulators of the circadian clock in this species. Three of these B. rapa genes, BrPRR1, BrPPR5, and BrPPR7, are present in two copies each in the B. rapa genome, for a total of eight B. rapa PRR (BrPRR) orthologs. We have now determined sequences and expression characteristics of the eight BrPRR genes and mapped their positions in the B. rapa genome. Although both members of each paralogous pair exhibited the same expression pattern, some variation in their gene structures was apparent. The BrPRR genes are tightly linked to several flowering genes. The knowledge about genome location, copy number variation and structural diversity of these B. rapa clock genes will improve our understanding of clock-related functions in this important crop. This will facilitate the development of Brassica crops for optimal growth in new environments and under changing conditions.

The polychaete, Perinereis nuntia ESTs and its use to uncover potential biomarker genes for molecular ecotoxicological studies.

The polychaete, Perinereis nuntia, has been used as an indicator species to assess the environmental condition of benthic communities in coastal marine environments. Recently, high-throughput sequencing technology has been proven to be a useful method for analyzing expressed sequence tags (ESTs) in non-model species. Thus, we have obtained extensive cDNA information by the pyrosequencing method, to utilize the polychaete species as a test organism for sediment quality monitoring studies. From the total RNA of P. nuntia, cDNA was reversely synthesized and randomly sequenced using a GS-FLX sequencer. In the assembly stage 1, 40,379 transcripts (13,666 contigs and 26,713 singletons) were acquired and showed 47% hitting rate compared with the GenBank non-redundant (NR) amino acid sequence database using BLASTX. After the stage-2 assembly, 21,657 transcripts were identified and showed 28% hitting rate. Finally, we obtained 6 064 unigenes that corresponded to the GenBank NR amino acid sequence database using BLASTX. Of the transcripts obtained in this species, we found a number of stress- and cell defense-related genes (e.g. heat shock protein family, antioxidant-related genes, cytochrome P450 genes) that are potentially useful for sediment monitoring at the molecular level, indicating that the pyrosequencing method is an effective approach to uncover gene families of potential biomarker genes simultaneously, and thus make transcriptomic studies possible. To confirm the usefulness of those potential biomarker genes, we analyzed the comparative profiling of P. nuntia mRNA transcripts between the samples collected from the polychaete aquaculture farm and the southern coast field of South Korea. In this paper, we summarize the expressed cDNA information of P. nuntia and discussed its potential use in environmental genomics and ecotoxicological studies for uncovering the potential molecular mechanisms of environmental stresses and chemical toxicity to the indicator species, P. nuntia in marine sediments.