Synthetic ShK-like Peptide from the Jellyfish Nemopilema nomurai Has Human Voltage-Gated Potassium-Channel-Blocking Activity.

We identified a new human voltage-gated potassium channel blocker, NnK-1, in the jellyfish Nemopilema nomurai based on its genomic information. The gene sequence encoding NnK-1 contains 5408 base pairs, with five introns and six exons. The coding sequence of the NnK-1 precursor is 894 nucleotides long and encodes 297 amino acids containing five presumptive ShK-like peptides. An electrophysiological assay demonstrated that the fifth peptide, NnK-1, which was chemically synthesized, is an effective blocker of hKv1.3, hKv1.4, and hKv1.5. Multiple-sequence alignment with cnidarian Shk-like peptides, which have Kv1.3-blocking activity, revealed that three residues (3Asp, 25Lys, and 34Thr) of NnK-1, together with six cysteine residues, were conserved. Therefore, we hypothesized that these three residues are crucial for the binding of the toxin to voltage-gated potassium channels. This notion was confirmed by an electrophysiological assay with a synthetic peptide (NnK-1 mu) where these three peptides were substituted with 3Glu, 25Arg, and 34Met. In conclusion, we successfully identified and characterized a new voltage-gated potassium channel blocker in jellyfish that interacts with three different voltage-gated potassium channels. A peptide that interacts with multiple voltage-gated potassium channels has many therapeutic applications in various physiological and pathophysiological contexts.

The complete mitochondrial genome of Rana coreana (Anura: Ranidae).

Rana coreana is a brown frog species native to the Korean Peninsula. We characterized the complete mitochondrial genome of the species. The mitochondrial genome sequence of R. coreana is 22,262 bp and comprises 13 protein-coding genes, two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and two control regions (CRs). The CR duplication and gene organization were identical to those observed in Rana kunyuensis and Rana amurensis. A total of 13 protein-coding genes were used to examine the phylogenetic relationships between this species and the genus Rana. R. coreana living on the Korean Peninsula, formed a cluster with R. kunyuensis and R. amurensis, with R. coreana showing the closest phylogenetic affinity for R. kunyuensis.

PharmaKoVariome database for supporting genetic testing.

Pharmacogenomics (PGx) provides information about routine precision medicine, based on the patient's genotype. However, many of the available information about human allele frequencies, and about clinical drug-gene interactions, is based on American and European populations. PharmaKoVariome database was constructed to support genetic testing for safe prescription and drug development. It consolidated and stored 2507 diseases, 11 459 drugs and 61 627 drug-target or druggable genes from public databases. PharmaKoVariome precomputed ethnic-specific abundant variants for approximately 120 M single-nucleotide variants of drug-target or druggable genes. A user can search by gene symbol, drug name, disease and reference SNP ID number (rsID) to statistically analyse the frequency of ethnical variations, such as odds ratio and P-values for related genes. In an example study, we observed five Korean-enriched variants in the CYP2B6 and CYP2D6 genes, one of which (rs1065852) is known to be incapable of metabolizing drug. It is also shown that 4-6% of North and East Asians have risk factors for drugs metabolized by the CYP2D6 gene. Therefore, PharmaKoVariome is a useful database for pharmaceutical or diagnostic companies for developing diagnostic technologies that can be applied in the Asian PGx industry. Database URL: http://www.pharmakovariome.com/.

Comparative Genome and Evolution Analyses of an Endangered Stony Coral Species Dendrophyllia cribrosa Near Dokdo Islands in the East Sea.

Stony corals often harbor intracellular photosynthetic dinoflagellate algae that receive dissolved inorganic nutrients. However, Dendrophyllia cribrosa is a nonsymbiotic stony coral distributed in the western Pacific. We assembled a chromosome-level D. cribrosa genome using PacBio and Hi-C technologies. The final assembly was 625 Mb, distributed on 14 chromosomes, and contained 30,493 protein-coding genes. The Benchmarking Universal Single-Copy Orthologs analysis revealed a percentage of 96.8 of the metazoan genome. A comparative phylogenetic analysis revealed that D. cribrosa, which lacks symbionts, evolved to acquire cellular energy by expanding genes related to acyl-CoA metabolism and carbohydrate transporters. This species also has expanded immune-related genes involved in the receptor protein tyrosine kinase signaling pathway. In addition, we observed a specific expansion of calcification genes, such as coral acid-rich proteins and carbonic anhydrase, in D. cribrosa. This high-quality reference genome and comparative analysis provides insights into the ecology and evolution of nonsymbiotic stony corals.

Chromosome-Level Genome Assembly of the Butter Clam Saxidomus purpuratus.

Herein, we provide the first whole-genome sequence of the purple butter clam (Saxidomus purpuratus), an economically important bivalve shellfish. Specifically, we sequenced and de novo assembled the genome of Sa. purpuratus based on PromethION long reads and Hi-C data. The 978-Mb genome of Sa. purpuratus comprises 19 chromosomes with 36,591 predicted protein-coding genes. The N50 length of Sa. purpuratus genome is 52 Mb, showing the highest continuous assembly among bivalve genomes. The Benchmarking by Universal Single-Copy Orthologs assessment indicated that 95.07% of complete metazoan universal single-copy orthologs (n = 954) were present in the assembly. Approximately 51% of Sa. purpuratus genome comprises repetitive sequences. Based on the high-quality Sa. purpuratus genome, we resolved half of the immune-associated genes, namely, scavenger receptor (SR) proteins, which are collinear to those in the closely related Cyclina sinensis genome. This finding suggested a high degree of conservation among immune-associated genes. Twenty-two (19%) SR proteins are tandemly duplicated in Sa. purpuratus genome, suggesting putative convergence evolution. Overall, Sa. purpuratus genome provides a new resource for the discovery of economically important traits and immune-response genes.

A chromosome-scale genome assembly and annotation of the spring orchid (Cymbidium goeringii).

Cymbidium goeringii, commonly known as the spring orchid, has long been favoured for horticultural purposes in Asian countries. It is a popular orchid with much demand for improvement and development for its valuable varieties. Until now, its reference genome has not been published despite its popularity and conservation efforts. Here, we report the de novo assembly of the C. goeringii genome, which is the largest among the orchids published to date, using a strategy that combines short- and long-read sequencing and chromosome conformation capture (Hi-C) information. The total length of all scaffolds is 3.99 Gb, with an N50 scaffold size of 178.2 Mb. A total of 29,556 protein-coding genes were annotated and 3.55 Gb (88.87% of genome) repetitive sequences were identified. We constructed pseudomolecular chromosomes using Hi-C, incorporating 89.4% of the scaffolds in 20 chromosomes. We identified 220 expanded and 106 contracted genes families in C. goeringii after divergence from its close relative. We also identified new gene families, resistance gene analogues and changes within the MADS-box genes, which control a diverse set of developmental processes during orchid evolution. Our high quality chromosomal-level assembly of C. goeringii can provide a platform for elucidating the genomic evolution of orchids, mining functional genes for agronomic traits and for developing molecular markers for accelerated breeding as well as accelerating conservation efforts.

Whole-genome, transcriptome, and methylome analyses provide insights into the evolution of platycoside biosynthesis in Platycodon grandiflorus, a medicinal plant.

Triterpenoid saponins (TSs) are common plant defense phytochemicals with potential pharmaceutical properties. Platycodon grandiflorus (Campanulaceae) has been traditionally used to treat bronchitis and asthma in East Asia. The oleanane-type TSs, platycosides, are a major component of the P. grandiflorus root extract. Recent studies show that platycosides exhibit anti-inflammatory, antiobesity, anticancer, antiviral, and antiallergy properties. However, the evolutionary history of platycoside biosynthesis genes remains unknown. In this study, we sequenced the genome of P. grandiflorus and investigated the genes involved in platycoside biosynthesis. The draft genome of P. grandiflorus is 680.1 Mb long and contains 40,017 protein-coding genes. Genomic analysis revealed that the CYP716 family genes play a major role in platycoside oxidation. The CYP716 gene family of P. grandiflorus was much larger than that of other Asterid species. Orthologous gene annotation also revealed the expansion of ß-amyrin synthases (bASs) in P. grandiflorus, which was confirmed by tissue-specific gene expression. In these expanded gene families, we identified key genes showing preferential expression in roots and association with platycoside biosynthesis. In addition, whole-genome bisulfite sequencing showed that CYP716 and bAS genes are hypomethylated in P. grandiflorus, suggesting that epigenetic modification of these two gene families affects platycoside biosynthesis. Thus whole-genome, transcriptome, and methylome data of P. grandiflorus provide novel insights into the regulation of platycoside biosynthesis by CYP716 and bAS gene families.

The Origin and Composition of Korean Ethnicity Analyzed by Ancient and Present-Day Genome Sequences.

Koreans are thought to be an ethnic group of admixed northern and southern subgroups. However, the exact genetic origins of these two remain unclear. In addition, the past admixture is presumed to have taken place on the Korean peninsula, but there is no genomic scale analysis exploring the origin, composition, admixture, or the past migration of Koreans. Here, 88 Korean genomes compared with 91 other present-day populations showed two major genetic components of East Siberia and Southeast Asia. Additional paleogenomic analysis with 115 ancient genomes from Pleistocene hunter-gatherers to Iron Age farmers showed a gradual admixture of Tianyuan (40 ka) and Devil's gate (8 ka) ancestries throughout East Asia and East Siberia up until the Neolithic era. Afterward, the current genetic foundation of Koreans may have been established through a rapid admixture with ancient Southern Chinese populations associated with Iron Age Cambodians. We speculate that this admixing trend initially occurred mostly outside the Korean peninsula followed by continuous spread and localization in Korea, corresponding to the general admixture trend of East Asia. Over 70% of extant Korean genetic diversity is explained to be derived from such a recent population expansion and admixture from the South.

Whole Genome Analysis of the Red-Crowned Crane Provides Insight into Avian Longevity.

The red-crowned crane (Grus japonensis) is an endangered, large-bodied crane native to East Asia. It is a traditional symbol of longevity and its long lifespan has been confirmed both in captivity and in the wild. Lifespan in birds is known to be positively correlated with body size and negatively correlated with metabolic rate, though the genetic mechanisms for the red-crowned crane's long lifespan have not previously been investigated. Using whole genome sequencing and comparative evolutionary analyses against the grey-crowned crane and other avian genomes, including the long-lived common ostrich, we identified redcrowned crane candidate genes with known associations with longevity. Among these are positively selected genes in metabolism and immunity pathways (NDUFA5, NDUFA8, NUDT12, SOD3, CTH , RPA1, PHAX, HNMT , HS2ST1 , PPCDC , PSTK CD8B, GP9, IL-9R, and PTPRC). Our analyses provide genetic evidence for low metabolic rate and longevity, accompanied by possible convergent adaptation signatures among distantly related large and long-lived birds. Finally, we identified low genetic diversity in the red-crowned crane, consistent with its listing as an endangered species, and this genome should provide a useful genetic resource for future conservation studies of this rare and iconic species.

The complete mitochondrial genome of a Dokdo shrimp, Lebbeus groenlandicus.

Lebbeus groenlandicus is a shrimp species indigenous to the Dokdo islands in the East Sea of Korea. We report the 17,399 bp mitochondrial genome (mitogenome) of the species that consists of 13 protein-coding genes, 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (rRNAs), and a control region (CR). A maximum-likelihood tree, constructed with 18 prawn and 45 shrimp mitogenomes, confirmed that L. groenlandicus occupies the most basal position within the Caridea infra-order and is closely related to Pandalidae shrimps.

Intracellular Ca(2+) and K(+) concentration in Brassica oleracea leaf induces differential expression of transporter and stress-related genes.

BACKGROUND: One of the most important members of the genus Brassica, cabbage, requires a relatively high level of calcium for normal growth (Plant Cell Environ 7: 397-405, 1984; Plant Physiol 60: 854-856, 1977). Localized Ca(2+) deficiency in cabbage leaves causes tip-burn, bringing about serious economic losses (Euphytica 9:203-208, 1960; Ann Bot 43:363-372, 1979; Sci Hortic 14:131-138, 1981). Although it has been known that the occurrence of tip-burn is related to Ca(2+) deficiency, there is limited information on the underlying mechanisms of tip-burn or the relationship between Ca(2+) and tip-burn incidence. To obtain more information on the genetic control of tip-burn symptoms, we focused on the identification of genes differentially expressed in response to increasing intracellular Ca(2+) and K(+) concentrations in B. oleracea lines derived from tip-burn susceptible, tip-burn resistant cabbages (B. oleracea var. capitata), and kale (B. oleracea var. acephala). RESULTS: We compared the levels of major macronutrient cations, including Ca(2+) and K(+), in three leaf segments, the leaf apex (LA), middle of leaf (LM), and leaf base (LB), of tip-burn susceptible, tip-burn resistant cabbages, and kale. Ca(2+) and K(+) concentrations were highest in kale, followed by tip-burn resistant and then tip-burn susceptible cabbages. These cations generally accumulated to a greater extent in the LB than in the LA. Transcriptome analysis identified 58,096 loci as putative non-redundant genes in the three leaf segments of the three B. oleracea lines and showed significant changes in expression of 27,876 loci based on Ca(2+) and K(+) levels. Among these, 1844 loci were identified as tip-burn related phenotype-specific genes. Tip-burn resistant cabbage and kale-specific genes were largely related to stress and transport activity based on GO annotation. Tip-burn resistant cabbage and kale plants showed phenotypes clearly indicative of heat-shock, freezing, and drought stress tolerance compared to tip-burn susceptible cabbages, demonstrating a correlation between intracellular Ca(2+) and K(+) concentrations and tolerance of abiotic stress with differential gene expression. We selected 165 genes that were up- or down-regulated in response to increasing Ca(2+) and K(+) concentrations in the three leaf segments of the three plant lines. Gene ontology enrichment analysis indicated that these genes participated in regulatory metabolic processes or stress responses. CONCLUSIONS: Our results indicate that the genes involved in regulatory metabolic processes or stress responses were differentially expressed in response to increasing Ca(2+) and K(+) concentrations in the B. oleracea leaf. Our transcriptome data and the genes identified may serve as a starting point for understanding the mechanisms underlying essential macronutrient deficiencies in plants, as well as the features of tip-burn in cabbage and other Brassica species.

Functional innovations of three chronological mesohexaploid Brassica rapa genomes.

BACKGROUND: The Brassicaceae family is an exemplary model for studying plant polyploidy. The Brassicaceae knowledge-base includes the well-annotated Arabidopsis thaliana reference sequence; well-established evidence for three rounds of whole genome duplication (WGD); and the conservation of genomic structure, with 24 conserved genomic blocks (GBs). The recently released Brassica rapa draft genome provides an ideal opportunity to update our knowledge of the conserved genomic structures in Brassica, and to study evolutionary innovations of the mesohexaploid plant, B. rapa. RESULTS: Three chronological B. rapa genomes (recent, young, and old) were reconstructed with sequence divergences, revealing a trace of recursive WGD events. A total of 636 fast evolving genes were unevenly distributed throughout the recent and young genomes. The representative Gene Ontology (GO) terms for these genes were 'stress response' and 'development' both through a change in protein modification or signaling, rather than by enhancing signal recognition. In retention patterns analysis, 98% of B. rapa genes were retained as collinear gene pairs; 77% of those were singly-retained in recent or young genomes resulting from death of the ancestral copies, while others were multi-retained as long retention genes. GO enrichments indicated that single retention genes mainly function in the interpretation of genetic information, whereas, multi-retention genes were biased toward signal response, especially regarding development and defense. In the recent genome, 13,302, 5,790, and 20 gene pairs were multi-retained following Brassica whole genome triplication (WGT) events with 2, 3, and 4 homoeologous copies, respectively. Enriched GO-slim terms from B. rapa homomoelogues imply that a major effect of the B. rapa WGT may have been to acquire environmental adaptability or to change the course of development. These homoeologues seem to more frequently undergo subfunctionalization with spatial expression patterns compared with other possible events including nonfunctionalization and neofunctionalization. CONCLUSION: We refined Brassicaceae GB information using the latest genomic resources, and distinguished three chronologically ordered B. rapa genomes. B. rapa genes were categorized into fast evolving, single- and multi-retention genes, and long retention genes by their substitution rates and retention patterns. Representative functions of the categorized genes were elucidated, providing better understanding of B. rapa evolution and the Brassica genus.

Small RNA and transcriptome deep sequencing proffers insight into floral gene regulation in Rosa cultivars.

BACKGROUND: Roses (Rosa sp.), which belong to the family Rosaceae, are the most economically important ornamental plants--making up 30% of the floriculture market. However, given high demand for roses, rose breeding programs are limited in molecular resources which can greatly enhance and speed breeding efforts. A better understanding of important genes that contribute to important floral development and desired phenotypes will lead to improved rose cultivars. For this study, we analyzed rose miRNAs and the rose flower transcriptome in order to generate a database to expound upon current knowledge regarding regulation of important floral characteristics. A rose genetic database will enable comprehensive analysis of gene expression and regulation via miRNA among different Rosa cultivars. RESULTS: We produced more than 0.5 million reads from expressed sequences, totalling more than 110 million bp. From these, we generated 35,657, 31,434, 34,725, and 39,722 flower unigenes from Rosa hybrid: 'Vital', 'Maroussia', and 'Sympathy' and Rosa rugosa Thunb., respectively. The unigenes were assigned functional annotations, domains, metabolic pathways, Gene Ontology (GO) terms, Plant Ontology (PO) terms, and MIPS Functional Catalogue (FunCat) terms. Rose flower transcripts were compared with genes from whole genome sequences of Rosaceae members (apple, strawberry, and peach) and grape. We also produced approximately 40 million small RNA reads from flower tissue for Rosa, representing 267 unique miRNA tags. Among identified miRNAs, 25 of them were novel and 242 of them were conserved miRNAs. Statistical analyses of miRNA profiles revealed both shared and species-specific miRNAs, which presumably effect flower development and phenotypes. CONCLUSIONS: In this study, we constructed a Rose miRNA and transcriptome database, and we analyzed the miRNAs and transcriptome generated from the flower tissues of four Rosa cultivars. The database provides a comprehensive genetic resource which can be used to better understand rose flower development and to identify candidate genes for important phenotypes.

A genome-wide comparison of NB-LRR type of resistance gene analogs (RGA) in the plant kingdom.

Plants express resistance (R) genes to recognize invaders and prevent the spread of pathogens. To analyze nucleotide binding site, leucine-rich repeat (NB-LRR) genes, we constructed a fast pipeline to predict and classify the R gene analogs (RGAs) by applying in-house matrices. With predicted ~37,000 RGAs, we can directly compare RGA contents across entire plant lineages, from green algae to flowering plants. We focused on the highly divergent NBLRRs in land plants following the emergence of mosses. We identified entire loss of Toll/Interleukin-1 receptor, NBLRR (TNL) in Poaceae family of monocots and interestingly from Mimulus guttatus (a dicot), which leads to the possibility of species-specific TNL loss in other sequenced flowering plants. Using RGA maps, we have elucidated a positive correlation between the cluster sizes of NB-LRRs and their numbers. The cluster members were observed to consist of the same class of NB-LRRs or their variants, which were probably generated from a single locus for an R gene. Our website ( http://sol.kribb.re.kr/PRGA/ ), called plant resistance gene analog (PRGA), provides useful information, such as RGA annotations, tools for predicting RGAs, and analyzing domain profiles. Therefore, PRGA provides new insights into R-gene evolution and is useful in applying RGA as markers in breeding and or systematic studies.

ASPMF: a new approach for identifying alternative splicing isoforms using peptide mass fingerprinting.

Alternative splicing is generally accepted as a mechanism that explains the discrepancy between the number of genes and proteins. We used peptide mass fingerprinting with a theoretical database and scoring method to discover and identify alternative splicing isoforms. Our theoretical database was built using published alternative splicing databases such as ECgene, H-DBAS, and TISA. According to our theoretical database of 190,529 isoforms, 37% of human genes have multiple isoforms. The isoforms produced from a gene partially share common peptide fragments because they have common exons, making it difficult to distinguish isoforms. Therefore, we developed a new method that effectively distinguishes a true isoform among multiple isoforms in a gene. In order to evaluate our algorithm, we made test sets for 4226 protein isoforms extracted from our theoretical database randomly. Consequently, 94% of true isoforms were identified by our scoring algorithm.

A combined approach for the classification of G protein-coupled receptors and its application to detect GPCR splice variants.

G protein-coupled receptors (GPCRs) constitute the largest family of cell surface receptors and play a central role in cellular signaling pathways. The importance of GPCRs has led to their becoming the targets of more than 50% of prescription drugs. However, drug compounds that do not differentiate between receptor subtypes can have considerable side effects and efficacy problems. An accurate classification of GPCRs can solve the side effect problems and raise the efficacy of drugs. Here, we introduce an approach that combines a fingerprint method, statistical profiles and physicochemical properties of transmembrane (TM) domains for a highly accurate classification of the receptors. The approach allows both the recognition and classification for GPCRs at the subfamily and subtype level, and allows the identification of splice variants. We found that the approach demonstrates an overall accuracy of 97.88% for subfamily classification, and 94.57% for subtype classification.

Clinical manifestations and diagnosis of extrapulmonary tuberculosis.

Since the diagnosis of extrapulmonary tuberculosis (EPT) is largely depended on the physician's suspicion in respect of the disease, we believed that it would be worthwhile to scrutinize the clinical characteristics of EPT. Thus, here we present retrospectively evaluated clinical manifestations of patients who were diagnosed as EPT cases in a tertiary referral care hospital. Medical records of 312 patients, diagnosed as having EPT at Yongdong Severance hospital from January 1997 to December 1999, were reviewed retrospectively. In total 312 patients, 149 (47.8%) males and 163 (52.2%) females aged from 13 years to 87 years, were included into this study. The most common site of the involvement was pleura (35.6%). The patients complained of localized symptoms (72.4%) more frequently than systemic symptoms (52.2%). The most common symptom was pain at the infected site (48.1%). Leukocytosis, anemia, and elevated erythrocyte sedimentation rate (ESR) and C-reactive protein (CRP) were found in 12.8%, 50.3%, 79.3% and 63.1% of the patients, respectively. Twenty-four percent of the patients had underlying medical illnesses such as, diabetes mellitus or liver cirrhosis, or were over 60 years old. In 67.3% of patients, tuberculosis was suspected at the initial visit. However, tuberculosis was microbiologically proven in only 23.7% of the patients. The time interval from the symptom onset to the diagnosis varied, with the mean duration of the period 96 days. Pulmonary parenchymal abnormal lesions were found in 133 patients (42.6%) on chest radiographs. EPT has a wide spectrum of clinical manifestations, so it is difficult to diagnose it. Based on our studies, only 11.2% of the patients were confirmed as EPT. So it is important that the physician who first examines the patient should have a high degree of suspicion based on the chest radiography, localized or systemic symptoms and several laboratory parameters reviewed in this study.