Search

In response to the COVID-19 pandemic, different types of vaccines, such as inactive, live-attenuated, messenger RNA, and protein subunit, have been developed against SARS-CoV2. This circumstance has unintentionally led to heterologous prime-boost vaccination against a single virus in a large human population. Here, we aimed to analyze whether the immunization order of vaccine types influences the efficacy of heterologous prime-boost vaccination, especially mRNA and protein-based vaccines. We developed a new mRNA vaccine expressing hemagglutinin (HA) of influenza using the 3′UTR and 5′UTR of muscle cells (mRNA-HA) and tested its efficacy by heterologous immunization with an HA protein vaccine (protein-HA). The results demonstrated higher IgG2a levels and hemagglutination inhibition titers in mRNA-HA priming/protein-HA boosting (R-P) regimen than that induced by reverse immunization (protein-HA priming/mRNA-HA boosting, P-R). After the virus challenge, the R-P group showed lower virus titers and less inflammation in the lungs than the P-R group. Transcriptome analysis revealed that heterologous prime-boost groups had differentially activated immune response pathways, according to the order of immunization. In summary, our results demonstrate that the sequence of vaccination is critical to sculpt immune responses. This study provides the potential of a heterologous vaccination strategy using mRNA and protein vaccine platforms against viral infection.

High-quality and complete reference genome assemblies are fundamental for the application of genomics to biology, disease, and biodiversity conservation. However, such assemblies are available for only a few non-microbial species1,2,3,4. To address this issue, the international Genome 10K (G10K) consortium5,6 has worked over a five-year period to evaluate and develop cost-effective methods for assembling highly accurate and nearly complete reference genomes. Here we present lessons learned from generating assemblies for 16 species that represent six major vertebrate lineages. We confirm that long-read sequencing technologies are essential for maximizing genome quality, and that unresolved complex repeats and haplotype heterozygosity are major sources of assembly error when not handled correctly. Our assemblies correct substantial errors, add missing sequence in some of the best historical reference genomes, and reveal biological discoveries. These include the identification of many false gene duplications, increases in gene sizes, chromosome rearrangements that are specific to lineages, a repeated independent chromosome breakpoint in bat genomes, and a canonical GC-rich pattern in protein-coding genes and their regulatory regions. Adopting these lessons, we have embarked on the Vertebrate Genomes Project (VGP), an international effort to generate high-quality, complete reference genomes for all of the roughly 70,000 extant vertebrate species and to help to enable a new era of discovery across the life sciences., We thank them for their permission to publish. A.R., S.K., B.P.W. and A.M.P. were supported by the Intramural Research Program of the NHGRI, NIH (1ZIAHG200398). A.R. was also supported by the Korea Health Technology R&D Project through KHIDI, funded by the Ministry of Health & Welfare, Republic of Korea (HI17C2098). S.A.M., I.B. and R.D. were supported by Wellcome Trust grant WT207492; W.C., M. Smith, Z.N., Y.S., J.C., S. Pelan, J.T., A.T., J.W. and Kerstin Howe by WT206194; L.H., F.M., Kevin Howe and P. Flicek by WT108749/Z/15/Z, WT218328/B/19/Z and the European Molecular Biology Laboratory. O.F. and E.D.J. were supported by Howard Hughes Medical Institute and Rockefeller University start-up funds for this project. J.D. and H.A.L. were supported by the Robert and Rosabel Osborne Endowment. M.U.-S. received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement (750747). F.T.-N., J. Hoffman, P. Masterson and K.C. were supported by the Intramural Research Program of the NLM, NIH. C.L., B.J.K., J. Kim and H.K. were supported by the Marine Biotechnology Program of KIMST, funded by the Ministry of Ocean and Fisheries, Republic of Korea (20180430). M.C. was supported by Sloan Research Fellowship (FG-2020-12932). S.C.V. was funded by a Max Planck Research Group award from the Max Planck Society, and a Human Frontiers Science Program (HFSP) Research grant (RGP0058/2016). T.M.L., W.E.J. and the Canada lynx genome were funded by the Maine Department of Inland Fisheries & Wildlife (F11AF01099), including when W.E.J. held a National Research Council Research Associateship Award at the Walter Reed Army Institute of Research (WRAIR). C.B. was supported by the NSF (1457541 and 1456612). D.B. was funded by The University of Queensland (HFSP - RGP0030/2015). D.I. was supported by Science Exchange Inc. (Palo Alto, CA). H.W.D. was supported by NSF grants (OPP-0132032 ICEFISH 2004 Cruise, PLR-1444167 and OPP-1955368) and the Marine Science Center at Northeastern University (416). G.J.P.N. and the thorny skate genome were funded by Lenfest Ocean Program (30884). M.P. was funded by the German Federal Ministry of Education and Research (01IS18026C). M. Malinsky was supported by an EMBO fellowship (ALTF 456-2016). The following authors’ contributions were supported by the NIH: S. Selvaraj (R44HG008118); C.V.M., S.R.F., P.V.L. (R21 DC014432/DC/NIDCD); K.D.M. (R01GM130691); H.C. (5U41HG002371-19); M.D. (U41HG007234); and B.P. (R01HG010485). D.G. was supported by the National Key Research and Development Program of China (2017YFC1201201, 2018YFC0910504 and 2017YFC0907503). F.O.A. was supported by Al-Gannas Qatari Society and The Cultural Village Foundation-Katara, Doha, State of Qatar and Monash University Malaysia. C.T. was supported by The Rockefeller University. M. Hiller was supported by the LOEWE-Centre for Translational Biodiversity Genomics (TBG) funded by the Hessen State Ministry of Higher Education, Research and the Arts (HMWK). H.C. was supported by the NHGRI (5U41HG002371-19). R.H.S.K. was funded by the Max Planck Society with computational resources at the bwUniCluster and BinAC funded by the Ministry of Science, Research and the Arts Baden-Württemberg and the Universities of the State of Baden-Württemberg, Germany (bwHPC-C5). B.V. was supported by the Biomedical Research Council of A*STAR, Singapore. T.M.-B. was funded by the European Research Council under the European Union’s Horizon 2020 research and innovation programme (864203), MINECO/FEDER, UE (BFU2017-86471-P), Unidad de Excelencia María de Maeztu, AEI (CEX2018-000792-M), a Howard Hughes International Early Career award, Obra Social “La Caixa” and Secretaria d’Universitats i Recerca and CERCA Programme del Departament d’Economia i Coneixement de la Generalitat de Catalunya (GRC 2017 SGR 880). E.C.T. was supported by the European Research Council (ERC-2012-StG311000) and an Irish Research Council Laureate Award. M.T.P.G. was supported by an ERC Consolidator Award 681396-Extinction Genomics, and a Danish National Research Foundation Center Grant (DNRF143). T.W. was supported by the NSF (1458652). J. M. Graves was supported by the Australian Research Council (CEO561477). E.W.M. was partially supported by the German Federal Ministry of Education and Research (01IS18026C). Complementary sequencing support for the Anna’s hummingbird and several genomes was provided by Pacific Biosciences, Bionano Genomics, Dovetail Genomics, Arima Genomics, Phase Genomics, 10X Genomics, NRGene, Oxford Nanopore Technologies, Illumina, and DNAnexus. All other sequencing and assembly were conducted at the Rockefeller University, Sanger Institute, and Max Planck Institute Dresden genome labs. Part of this work used the computational resources of the NIH HPC Biowulf cluster (https://hpc.nih.gov). We acknowledge funding from the Wellcome Trust (108749/Z/15/Z) and the European Molecular Biology Laboratory., With funding from the Spanish government through the "Severo Ochoa Centre of Excellence" accreditation (CEX2018-000792-M).

Obesity causes brain injuries with inflammatory and structural changes, leading to neurodegeneration. Although increased circulating lipocalin 2 (LCN2) level has been implicated in neurodegenerative diseases, the precise mechanism of neurodegeneration in obesity is not clear. Here, we investigated whether LCN2-mediated signaling promotes neurodegeneration in the hippocampus of leptin-deficient ob/ob mice, which are characterized by obesity, insulin resistance, systemic inflammation, and neuroinflammation. In particular, there was significant upregulation of both LCN2 and matrix metalloproteinase 9 levels from serum and hippocampus in ob/ob mice. Using RNA-seq analysis, we found that neurodegeneration- sortilin-related receptor 1 (Sorl1) and brain-derived neurotrophic factor (Bdnf) genes were significantly reduced in the hippocampus of ob/ob mice. We additionally found that the endosome-related WD repeat and FYVE-domain-containing 1 (Wdfy1) gene were upregulated in ob/ob mice. In particular, iron overload-related mitochondrial ferritin and nuclear factor erythroid 2-related factor 2 (Nrf2)/heme oxygenase-1 (HO-1) proteins were increased in the hippocampus of ob/ob. Thus, these findings indicate that iron-binding protein LCN2-mediated oxidative stress promotes neurodegeneration in ob/ob mice.

Recent development of long-read sequencing platforms has enabled researchers to explore bacterial community structure through analysis of full-length 16S rRNA gene (∼1,500 bp) or 16S-ITS-23S rRNA operon region (∼4,300 bp), resulting in higher taxonomic resolution than short-read sequencing platforms. Despite the potential of long-read sequencing in metagenomics, resources and protocols for this technology are scarce. Here, we describe MIrROR, the database and analysis tool for metataxonomics using the bacterial 16S-ITS-23S rRNA operon region. We collected 16S-ITS-23S rRNA operon sequences extracted from bacterial genomes from NCBI GenBank and performed curation. A total of 97,781 16S-ITS-23S rRNA operon sequences covering 9,485 species from 43,653 genomes were obtained. For user convenience, we provide an analysis tool based on a mapping strategy that can be used for taxonomic profiling with MIrROR database. To benchmark MIrROR, we compared performance against publicly available databases and tool with mock communities and simulated data sets. Our platform showed promising results in terms of the number of species covered and the accuracy of classification. To encourage active 16S-ITS-23S rRNA operon analysis in the field, BLAST function and taxonomic profiling results with 16S-ITS-23S rRNA operon studies, which have been reported as BioProject on NCBI are provided. MIrROR (http://mirror.egnome.co.kr/) will be a useful platform for researchers who want to perform high-resolution metagenome analysis with a cost-effective sequencer such as MinION from Oxford Nanopore Technologies.

Insights into the evolution of non-model organisms are often limited by the lack of reference genomes. As part of the Vertebrate Genomes Project, we present a new reference genome and a pangenome produced with High-Fidelity long reads for the barn swallow Hirundo rustica. We then generated a reference-free multialignment with other bird genomes to identify genes under selection. Conservation analyses pointed at genes enriched for transcriptional regulation and neurodevelopment. The most conserved gene is CAMK2N2, with a potential role in fear memory formation. In addition, using all publicly available data, we generated a comprehensive catalogue of genetic markers. Genome-wide linkage disequilibrium scans identified potential selection signatures at multiple loci. The top candidate region comprises several genes and includes BDNF, a gene involved in stress response, fear memory formation, and tameness. We propose that the strict association with humans in this species is linked with the evolution of pathways typically under selection in domesticated taxa.

Breast cancer remains a common cause of cancer-related death in women. Therefore, further studies are necessary for the comprehension of breast cancer and the revolution of breast cancer treatment. Cancer is a heterogeneous disease that results from epigenetic alterations in normal cells. Aberrant epigenetic regulation is strongly associated with the development of breast cancer. Current therapeutic approaches target epigenetic alterations rather than genetic mutations due to their reversibility. The formation and maintenance of epigenetic changes depend on specific enzymes, including DNA methyltransferases and histone deacetylases, which are promising targets for epigenetic-based therapy. Epidrugs target different epigenetic alterations, including DNA methylation, histone acetylation, and histone methylation, which can restore normal cellular memory in cancerous diseases. Epigenetic-targeted therapy using epidrugs has anti-tumor effects on malignancies, including breast cancer. This review focuses on the importance of epigenetic regulation and the clinical implications of epidrugs in breast cancer.

We report 37 mitochondrial genome (mitogenome) sequences of Bombyx mori strains (Lepidoptera: Bombycidae) and four of B. mandarina individuals, each preserved and collected, respectively, in South Korea. These mitogenome sequences combined with 45 public data showed a substantial genetic reduction in B. mori strains compared to the presumed ancestor B. mandarina, with the highest diversity detected in the Chinese origin B. mori. Chinese B. mandarina were divided into northern and southern groups, concordant to the Qinling–Huaihe line, and the northern group was placed as an immediate progenitor of monophyletic B. mori strains in phylogenetic analyses, as has previously been detected. However, one individual that was in close proximity to the south Qinling–Huaihe line was exceptional, belonging to the northern group. The enigmatic South Korean population of B. mandarina, which has often been regarded as a closer genetic group to Japan, was most similar to the northern Chinese group, evidencing substantial gene flow between the two regions. Although a substantial genetic divergence is present between B. mandarina in southern China and Japan, a highly supported sister relationship between the two regional populations may suggest the potential origin of Japanese B. mandarina from southern China instead of the Korean peninsula.

Simple Summary Bombyx mandarina (Lepidoptera: Bombycidae), the presumed ancestor of B. mori, has long been a subject of study to illustrate the geographic variation in relationship to origin of B. mori. We report 37 mitochondrial genome (mitogenome) sequences of B. mori strains obtained by whole-genome sequencing and four of B. mandarina individuals obtained by Sanger sequencing from South Korea. These mitogenome sequences were combined with 45 public data to uncover the population genetic and phylogenetic relationships among regional populations of B. mandarina and their relationships to B. mori. A substantial genetic reduction in B. mori strains compared to the B. mandarina was detected, with the highest diversity detected in the Chinese origin B. mori. Chinese B. mandarina were divided into northern and southern groups, largely concordant to the Qinling–Huaihe line, and the northern group was placed as an immediate progenitor of monophyletic B. mori strains in phylogenetic analyses. The enigmatic South Korean population of B. mandarina, which has often been regarded as a closer genetic group to Japan, was most similar to the northern Chinese group, evidencing substantial gene flow between the two regions. This is the first report that the South Korean B. mandarina are closer genetically to the northern Chinese group. Abstract We report 37 mitochondrial genome (mitogenome) sequences of Bombyx mori strains (Lepidoptera: Bombycidae) and four of B. mandarina individuals, each preserved and collected, respectively, in South Korea. These mitogenome sequences combined with 45 public data showed a substantial genetic reduction in B. mori strains compared to the presumed ancestor B. mandarina, with the highest diversity detected in the Chinese origin B. mori. Chinese B. mandarina were divided into northern and southern groups, concordant to the Qinling–Huaihe line, and the northern group was placed as an immediate progenitor of monophyletic B. mori strains in phylogenetic analyses, as has previously been detected. However, one individual that was in close proximity to the south Qinling–Huaihe line was exceptional, belonging to the northern group. The enigmatic South Korean population of B. mandarina, which has often been regarded as a closer genetic group to Japan, was most similar to the northern Chinese group, evidencing substantial gene flow between the two regions. Although a substantial genetic divergence is present between B. mandarina in southern China and Japan, a highly supported sister relationship between the two regional populations may suggest the potential origin of Japanese B. mandarina from southern China instead of the Korean peninsula.

Bifidobacterium longum KACC 91563 secretes family 5 extracellular solute-binding protein via extracellular vesicle. In our previous work, it was demonstrated that the protein effectively alleviated food allergy symptoms via mast cell specific apoptosis, and it has revealed a therapeutic potential of this protein in allergy treatment. In the present study, we cloned the gene encoding extracellular solute-binding protein of the strain into the histidine-tagged pET-28a(+) vector and transformed the resulting plasmid into the Escherichia coli strain BL21 (DE3). The histidine-tagged extracellular solute-binding protein expressed in the transformed cells was purified using Ni-NTA affinity column. To enhance the efficiency of the protein purification, three parameters were optimized; the host bacterial strain, the culturing and induction temperature, and the purification protocol. After the process, two liters of transformed culture produced 7.15 mg of the recombinant proteins. This is the first study describing the production of extracellular solute-binding protein of probiotic bacteria. Establishment of large-scale production strategy for the protein will further contribute to the development of functional foods and potential alternative treatments for allergies.

Assembling high-quality microbial genomes using only cost-effective Nanopore long-read systems such as Flongle is important to accelerate research on the microbial genome and the most critical point for this is the polishing process. In this study, we performed an evaluation based on BUSCO and Prokka gene prediction in terms of microbial genome assembly for eight state-of-the-art Nanopore polishing tools and combinations available. In the evaluation of individual tools, Homopolish, PEPPER, and Medaka demonstrated better results than others. In combination polishing, the second round Homopolish, and the PEPPER × medaka combination also showed better results than others. However, individual tools and combinations have specific limitations on usage and results. Depending on the target organism and the purpose of the downstream research, it is confirmed that there remain some difficulties in perfectly replacing the hybrid polishing carried out by the addition of a short-read. Nevertheless, through continuous improvement of the protein pores, related base-calling algorithms, and polishing tools based on improved error models, a high-quality microbial genome can be achieved using only Nanopore reads without the production of additional short-read data. The polishing strategy proposed in this study is expected to provide useful information for assembling the microbial genome using only Nanopore reads depending on the target microorganism and the purpose of the research.

In this study, complete nucleotide sequences of the mitochondrial genome of a nudibranch species, Hypselodoris festiva (A. Adams, 1861) were determined and characterized. The mitogenome size is 14 880 bp. This is the longest among the known nudibranch mitochondrial genomes. Furthermore, phylogenetic relationship of H. festiva in the Nudibranchia reconstructed due to amino acid sequences of mitochondrial protein coding genes. This is the sixth record for complete mitochondrial genome of the Nudibranchia.

Complete mitochondrial genome of the knobbed rock shell sea snail Menathais tuberosa (Röding, 1798) has been sequenced and phylogenetic relationships evaluated due to mitochondrial protein coding genes. The size of mitochondrial genome for M. tuberosa is 15,294 bp and the nucleotide composition of the mitochondrial genome is 28.4% A, 16.5% C, 17.6% G and 37.5% T. Reconstructed phylogenetic tree of the Neogastropoda showed that M. tuberosa is in the monophyletic Muricidae. This is the first record of complete mitochondrial genome from the genus.

Background Modern sequencing technologies should make the assembly of the relatively small mitochondrial genomes an easy undertaking. However, few tools exist that address mitochondrial assembly directly. Results As part of the Vertebrate Genomes Project (VGP) we develop mitoVGP, a fully automated pipeline for similarity-based identification of mitochondrial reads and de novo assembly of mitochondrial genomes that incorporates both long (> 10 kbp, PacBio or Nanopore) and short (100–300 bp, Illumina) reads. Our pipeline leads to successful complete mitogenome assemblies of 100 vertebrate species of the VGP. We observe that tissue type and library size selection have considerable impact on mitogenome sequencing and assembly. Comparing our assemblies to purportedly complete reference mitogenomes based on short-read sequencing, we identify errors, missing sequences, and incomplete genes in those references, particularly in repetitive regions. Our assemblies also identify novel gene region duplications. The presence of repeats and duplications in over half of the species herein assembled indicates that their occurrence is a principle of mitochondrial structure rather than an exception, shedding new light on mitochondrial genome evolution and organization. Conclusions Our results indicate that even in the “simple” case of vertebrate mitogenomes the completeness of many currently available reference sequences can be further improved, and caution should be exercised before claiming the complete assembly of a mitogenome, particularly from short reads alone.

Additional file 8: Supplementary Note 1. Relationship between mtDNA sequencing, coverage and genome size.

Additional file 12. Review history.

Additional file 1: Fig. S1. Outline of the mitoVGP assembly pipeline. Fig. S2. Assembly success by the availability of long mtDNA reads. Fig. S3. PacBio CLR mitochondrial read counts in different tissues. Fig. S4. mitoVGP assembly results and comparisons. Fig. S5. Correlation between differences in repeat content and assembly length between the mitoVGP versus the Genbank/Refseq assemblies. Fig. S6. Paired comparisons of GC content between the mitoVGP assemblies and their Genbank/RefSeq counterparts. Fig. S7. Heatmap of k-mer-based sequence similarity of repetitive elements. Fig. S8. Read length deviation from the reference in kbp. Fig. S9. Long accurate HiFi reads from the VGP human trio mapped to the child mitogenome assembly.

Modern sequencing technologies should make the assembly of the relatively small mitochondrial genomes an easy undertaking. However, few tools exist that address mitochondrial assembly directly. As part of the Vertebrate Genomes Project (VGP) we have developed mitoVGP, a fully automated pipeline for similarity-based identification of mitochondrial reads and de novo assembly of mitochondrial genomes that incorporates both long (>10 kbp, PacBio or Nanopore) and short (100-300 bp, Illumina) reads. Our pipeline led to successful complete mitogenome assemblies of 100 vertebrate species of the VGP. We have observed that tissue type and library size selection have considerable impact on mitogenome sequencing and assembly. Comparing our assemblies to purportedly complete reference mitogenomes based on short-read sequencing, we have identified errors, missing sequences, and incomplete genes in those references, particularly in repeat regions. Our assemblies have also identified novel gene region duplications, shedding new light on mitochondrial genome evolution and organization.

Continuous monitoring of the present genetic status is essential to preserve the genetic resource of wild populations. In this study, we sequenced regional Pacific abalone Haliotis discus samples from three different locations around the Korean peninsula to assess population structure, utilizing Genotyping-by-Sequencing (GBS) method. Using PstI enzyme for genome reduction, we demonstrated the resultant library represented the whole genome region with even spacing, and as a result 16,603 single nucleotide variants (SNVs) were produced. Genetic diversity and population structure were investigated using several methods, and a strong genetic heterogeneity was observed in the Korean abalone populations. Additionally, by comparison of the variant sets among population groups, we were able to discover 26 Korean abalone population-specific SNVs, potentially associated with phenotype differences. This is the first study demonstrating the feasibility of GBS for population genetic study on H. discus. Our results will provide valuable data for the genetic conservation and management of wild abalone populations in Korea and help future GBS studies on the marine mollusks.

High-quality and complete reference genome assemblies are fundamental for the application of genomics to biology, disease, and biodiversity conservation. However, such assemblies are only available for a few non-microbial species1–4. To address this issue, the international Genome 10K (G10K) consortium5,6 has worked over a five-year period to evaluate and develop cost-effective methods for assembling the most accurate and complete reference genomes to date. Here we summarize these developments, introduce a set of quality standards, and present lessons learned from sequencing and assembling 16 species representing major vertebrate lineages (mammals, birds, reptiles, amphibians, teleost fishes and cartilaginous fishes). We confirm that long-read sequencing technologies are essential for maximizing genome quality and that unresolved complex repeats and haplotype heterozygosity are major sources of error in assemblies. Our new assemblies identify and correct substantial errors in some of the best historical reference genomes. Adopting these lessons, we have embarked on the Vertebrate Genomes Project (VGP), an effort to generate high-quality, complete reference genomes for all ~70,000 extant vertebrate species and help enable a new era of discovery across the life sciences.

Despite the importance of Lactobacillus iners and its unique characteristics for the study of vaginal adaption, its genome and genomic researches for identifying molecular backgrounds of these specific phenotypes are still limited. In this study, the first complete genome of L. iners was constructed using a cost-effective long-read sequencing platform, Flongle from Oxford Nanopore, and comparative genome analysis was conducted using a total of 1,046 strain genomes from 10 vaginal Lactobacillus species. Single-molecule sequencing using Flongle effectively resolved the limitation of the 2nd generation sequencing technologies in dealing with genomic regions of high GC contents, and comparative genome analysis identified three potential core genes (INY, ZnuA, and hsdR) of L. iners which was related to its specific adaption to the vaginal environment. In addition, we performed comparative prophage analysis for 1,046 strain genomes to further identify the species specificity. The number of prophages in L. iners genomes was significantly smaller than other vaginal Lactobacillus species, and one of the specific genes (hsdR) was suggested as the means for defense against bacteriophage. The first complete genome of L. iners and the three specific genes identified in this study will provide useful resources to further expand our knowledge of L. iners and its specific adaption to the vaginal econiche.

Leptin-deficient and leptin-resistant mice manifest obesity, insulin resistance, and left ventricular hypertrophy (LVH); however, LVH’s mechanisms are not fully understood. Cardiac iron dysregulation has been recently implicated in cardiomyopathy. Here we investigated the protective effects of caloric restriction on cardiac remodeling in impaired leptin signaling obese mice. RNA-seq analysis was performed to assess the differential gene expressions in the heart of wild-type and ob/ob mice. In particular, to investigate the roles of caloric restriction on iron homeostasis-related gene expressions, 10-week-old ob/ob and db/db mice were assigned to ad libitum or calorie-restricted diets for 12 weeks. Male ob/ob mice exhibited LVH, cardiac inflammation, and oxidative stress. Using RNA-seq analysis, we identified that an iron uptake-associated gene, transferrin receptor, was upregulated in obese ob/ob mice with LVH. Caloric restriction attenuated myocyte hypertrophy, cardiac inflammation, fibrosis, and oxidative stress in ob/ob and db/db mice. Furthermore, we found that caloric restriction reversed iron homeostasis-related lipocalin 2, divalent metal transporter 1, transferrin receptor, ferritin, ferroportin, and hepcidin expressions in the heart of ob/ob and db/db mice. These findings demonstrate that the cardioprotective effects of caloric restriction result from the cellular regulation of iron homeostasis, thereby decreasing oxidative stress, inflammation, and cardiac remodeling. We suggest that decreasing iron-mediated oxidative stress and inflammation offers new therapeutic approaches for obesity-induced cardiomyopathy.

Lactobacillus plantarum is one of the widely-used probiotics and there have been a large number of advanced researches on the effectiveness of this species. However, the difference between previously reported plantarum strains, and the source of genomic variation among the strains were not clearly specified. In order to understand further on the molecular basis of L. plantarum on Korean traditional fermentation, we isolated the L. plantarum GB-LP4 from Korean fermented vegetable and conducted whole genome assembly. With comparative genomics approach, we identified the candidate genes that are expected to have undergone evolutionary acceleration. These genes have been reported to associate with the maintaining homeostasis, which are generally known to overcome instability in external environment including low pH or high osmotic pressure. Here, our results provide an evolutionary relationship between L. plantarum species and elucidate the candidate genes that play a pivotal role in evolutionary acceleration of GB-LP4 in high osmolarity environment. This study may provide guidance for further studies on L. plantarum.

The Korean native horse (Jeju horse) is one of the most important animals in Korean historical, cultural, and economical viewpoints. In the early 1980s, the Jeju horse was close to extinction. The aim of this study is to explore the phylogenomics of Korean native horse focusing on spatio-temporal dynamics. We determined complete mitochondrial genome sequences for the first Korean native (n = 6) and additional Mongolian (n = 2) horses. Those sequences were analyzed together with 143 published ones using Bayesian coalescent approach as well as three different phylogenetic analysis methods, Bayesian inference, maximum likelihood, and neighbor-joining methods. The phylogenomic trees revealed that the Korean native horses had multiple origins and clustered together with some horses from four European and one Middle Eastern breeds. Our phylogenomic analyses also supported that there was no apparent association between breed or geographic location and the evolution of global horses. Time of the most recent common ancestor of the Korean native horse was approximately 13,200-63,200 years, which was much younger than 0.696 My of modern horses. Additionally, our results showed that all global horse lineages including Korean native horse existed prior to their domestication events occurred in about 6000-10,000 years ago. This is the first study on phylogenomics of the Korean native horse focusing on spatio-temporal dynamics. Our findings increase our understanding of the domestication history of the Korean native horses, and could provide useful information for horse conservation projects as well as for horse genomics, emergence, and the geographical distribution.

Yersinia enterocolitica is a well-known foodborne pathogen causing gastrointestinal infections worldwide. The strain Y. enterocolitica FORC_002 was isolated from the gill of flatfish (plaice) and its genome was sequenced. The genomic DNA consists of 4,837,317 bp with a GC content of 47.1%, and is predicted to contain 4,221 open reading frames, 81 tRNA genes, and 26 rRNA genes. Interestingly, genomic analysis revealed pathogenesis and host immune evasion-associated genes encoding guanylate cyclase (Yst), invasin (Ail and Inv), outer membrane protein (Yops), autotransporter adhesin A (YadA), RTX-like toxins, and a type III secretion system. In particular, guanylate cyclase is a heat-stable enterotoxin causing Yersinia-associated diarrhea, and RTX-like toxins are responsible for attachment to integrin on the target cell for cytotoxic action. This genome can be used to identify virulence factors that can be applied for the development of novel biomarkers for the rapid detection of this pathogen in foods.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Mitochondrial genomes were sequenced from five Raja pulchra individuals, and single-nucleotide polymorphisms (SNPs) were identified by comparing previously announced sequences in this study. Total 117 SNPs were detected and they were present in 2 rRNA genes, 9 tRNA genes, 13 protein coding genes and non-coding region. One deleted polymorphic site, which was located in 16S rRNA gene, was observed in two individuals. Six polymorphic sites were non-synonymous SNPs, which were distributed in ND1, ND2, ATP6 and ND4 gene. Phylogenic analysis validated current taxa. The genome sequences of R. pulchra mitochondria could be comparable information for understanding species divergence and genomic variation among the populations.

Mitochondrial genomes were sequenced from five Raja pulchra individuals, and single-nucleotide polymorphisms (SNPs) were identified by comparing previously announced sequences in this study. Total 117 SNPs were detected and they were present in 2 rRNA genes, 9 tRNA genes, 13 protein coding genes and non-coding region. One deleted polymorphic site, which was located in 16S rRNA gene, was observed in two individuals. Six polymorphic sites were non-synonymous SNPs, which were distributed in ND1, ND2, ATP6 and ND4 gene. Phylogenic analysis validated current taxa. The genome sequences of R. pulchra mitochondria could be comparable information for understanding species divergence and genomic variation among the populations.

Red sea bream, a popular fish resource in Korea and Japan, is being bred in fish farms of the two countries. It is hypothesized that the genomes of red sea bream are influenced by decades of artificial selection. This study investigates the impact of artificial selection on genomes of red sea bream. Whole genome sequencing was conducted for 40 samples of red sea bream either from Ehime, Nagasaki and Tongyeong fish farms or from the wild. Population stratification based on whole genome data was investigated and the genomic regions of fish farm populations under selection were identified using XP-EHH and relative nucleotide diversity. Gene ontology analysis revealed that different functions were enriched in different fish farms. In conclusion, this study highlights the difference between independently cultured red sea bream populations by showing that influence of artificial selection acted upon completely different genes related to different functions including metabolic and developmental processes.

Bifidobacterium longum KACC 91563 secretes family 5 extracellular solute-binding protein via extracellular vesicle. In our previous work, it was demonstrated that the protein effectively alleviated food allergy symptoms via mast cell specific apoptosis, and it has revealed a therapeutic potential of this protein in allergy treatment. In the present study, we cloned the gene encoding extracellular solute-binding protein of the strain into the histidine-tagged pET-28a(+) vector and transformed the resulting plasmid into the Escherichia coli strain BL21 (DE3). The histidine-tagged extracellular solute-binding protein expressed in the transformed cells was purified using Ni-NTA affinity column. To enhance the efficiency of the protein purification, three parameters were optimized; the host bacterial strain, the culturing and induction temperature, and the purification protocol. After the process, two liters of transformed culture produced 7.15 mg of the recombinant proteins. This is the first study describing the production of extracellular solute-binding protein of probiotic bacteria. Establishment of large-scale production strategy for the protein will further contribute to the development of functional foods and potential alternative treatments for allergies.

Acute respiratory virus infectious diseases are a growing health problem, particularly among children and the elderly. Much effort has been made to develop probiotics that prevent influenza virus infections by enhancing innate immunity in the respiratory tract until vaccines are available. Lactobacillus plantarum GB-LP2, isolated from a traditional Korean fermented vegetable, has exhibited preventive effects on influenza virus infection in mice. To identify the molecular basis of this strain, we conducted a whole-genome assembly study. The single circular DNA chromosome of 3,284,304 bp was completely assembled and 3,250 proteinencoding genes were predicted. Evolutionarily accelerated genes related to the phenotypic trait of anti-infective activities for influenza virus were identified. These genes encode three integral membrane proteins, a teichoic acid export ATP-binding protein and a glucosamine - fructose-6-phosphate aminotransferase involved in host innate immunity, the nonspecific DNA-binding protein Dps, which protects bacteria from oxidative damage, and the response regulator of the three-component quorum-sensing regulatory system, which is related to the capacity of adhesion to the surface of the respiratory tract and competition with pathogens. This is the first study to identify the genetic backgrounds of the antiviral activity in L. plantarum strains. These findings provide insight into the anti-infective activities of L. plantarum and the development of preventive probiotics.

As a means of improving Korean Soil Classification System, we have reviewed Australian Soil Classification System by comparing Soil Taxonomy and FAO/WRB Classification System. Australian Soil Classification System is composed of 14 of Order, 87 of Sub-order, 556 of Great-group, 2,451 of Sub-group, and 7,276 of Family. Interestingly, soil order has the Anthroposols which is not classified with Soil Taxonomy, and the classification for some of soils is based on soil texture abruption horizon and soil structure. Seven of 14 soil orders are classified with an old version based on soil color rather than morphological characteristics. The distribution scale of Australian soil order is the largest in Tenosols, and followed by Kandosols, Rudosols, Sodosols and Vertisols in Australia.Key words: Soil Classification, Australia, Soil, Soil Orders

The pinnipeds, which comprise seals, sea lions, and walruses, are a remarkable group of marine animals with unique adaptations to semi-aquatic life. However, their genomes are poorly characterized. In this study, we sequenced and characterized the genomes of three pinnipeds (Phoca largha, Callorhinus ursinus, and Eumetopias jubatus), focusing on site-wise sequence changes. We detected rapidly evolving genes in pinniped lineages and substitutions unique to pinnipeds associated with amphibious sound perception. Phenotypic convergence-related sequence convergences are not common in marine mammals. For example, FASN, KCNA5, and IL17RA contain substitutions specific to pinnipeds, yet are potential candidates of phenotypic convergence (blubber, response to hypoxia, and immunity to pathogens) in all marine mammals. The outcomes of this study will provide insight into targets for future studies of convergent evolution or gene function.

Identifying the microbes present in probiotic products is an important issue in product quality control and public health. The most common methods used to identify genera containing species that produce lactic acid are matrix-assisted laser desorption/ionization–time of flight mass spectrometry (MALDI-TOF MS) and 16S rRNA sequence analysis. However, the high cost of operation, difficulty in distinguishing between similar species, and limitations of the current sequencing technologies have made it difficult to obtain accurate results using these tools. To overcome these problems, a whole-genome shotgun sequencing approach has been developed along with various metagenomic classification tools. Widely used tools include the marker gene and k-mer methods, but their inevitable false-positives (FPs) hampered an accurate analysis. We therefore, designed a coverage-based pipeline to reduce the FP problem and to achieve a more reliable identification of species. The coverage-based pipeline described here not only shows higher accuracy for the detection of species and proportion analysis, based on mapping depth, but can be applied regardless of the sequencing platform. We believe that the coverage-based pipeline described in this study can provide appropriate support for probiotic quality control, addressing current labeling issues.

Background Identification of genetic mechanisms and idiosyncrasies at the breed-level can provide valuable information for potential use in evolutionary studies, medical applications, and breeding of selective traits. Here, we analyzed genomic data collected from 136 Korean Native cattle, known as Hanwoo, using advanced statistical methods. Results Results revealed Hanwoo-specific protein domains which were largely characterized by immunoglobulin function. Furthermore, domain interactions of novel Hanwoo-specific genes reveal additional links to immunity. Novel Hanwoo-specific genes linked to muscle and other functions were identified, including protein domains with functions related to energy, fat storage, and muscle function that may provide insight into the mechanisms behind Hanwoo cattle’s uniquely high percentage of intramuscular fat and fat marbling. Conclusion The identification of Hanwoo-specific genes linked to immunity are potentially useful for future medical research and selective breeding. The significant genomic variations identified here can crucially identify genetic novelties that are arising from useful adaptations. These results will allow future researchers to compare and classify breeds, identify important genetic markers, and develop breeding strategies to further improve significant traits. Electronic supplementary material The online version of this article (10.1186/s12863-018-0623-x) contains supplementary material, which is available to authorized users.

Decapentaplegic (dpp) is a member of the transforming growth factor-β superfamily. Although the dpp gene and related pathways are known to play important roles in insect development, few studies have examined its function in Bombyx mori and Bombyx mandarina. To date, there have been no previous reports on novel splice variants of dpp in silkworm. In the present study, we conducted RT-PCR to examine dpp expression in the mid-gut tissue of B. mandarina and discovered a novel dpp isoform. The isoform sequence was confirmed using sequencing analysis and found to have 333 bp deletion compared to full-length cDNA encoding dpp. The deleted sequence encodes a region of the latency associated peptide (LAP) region of transforming growth factor-β (TGF-β), which may affect the activity and specificity of TGF-β. Using variant calling analyses, we detected 7 candidate single nucleotide variants (SNVs) for different alternative splicing in dpp. This is the first report of a novel splice variant of the dpp gene in B. mandarina and these results provide insight about the domestication process and distinct phenotypic traits of B. mori and B. mandarina.

Haliotis discus hannai, a species of Pacific abalone, is a highly valuable food source throughout Northeast Asia. As H. discus hannai primarily feed on brown algae and largely extract their energy from algal polysaccharides, understanding the mechanisms by which they digest algal polysaccharides is essential for elucidating their energy metabolism. Gut microbes, as well as the host animal, are involved in the process of polysaccharide degradation. To identify algal polysaccharide-digestion mechanisms and their origin, we analyzed the metagenome and metatranscriptome of abalone visceral extracts. Microbial communities were characterized using the 16S rRNA gene sequences in the metagenome and our results differed significantly from those of previous studies using traditional microbiological methods such as bacterial cultivation and cloning. A greater diversity of bacterial taxa was identified here than was previously identified using cultivation methods. Furthermore, the most abundant bacterial taxa also differed from previous studies, which is not common when comparing the results of bacterial culturing with those of molecular methodologies. Based on the metatranscriptome, overall functions were identified and additional analyses were performed on the coding sequences of algal polysaccharide-digestive enzymes, including alginate lyase. Results of the transcriptomic analyses suggest that alginate lyase in the visceral extracts of H. discus hannai was produced by the host itself, not by visceral bacteria. This is the first next-generation sequencing study performed on abalone to characterize the visceral microbiota and the source of the ability to digest algal polysaccharides by analyzing the metagenome and metatranscriptome together.

Table S2. Significantly identified (E- value

Table S1. Summary of sequencing data (DOCX 28 kb)

Staphylococcus aureus is an important foodborne pathogen that causes diverse diseases ranging from minor infections to life-threatening conditions in humans and animals. To further understand its pathogenesis, the genome of the strain S. aureus FORC_001 was isolated from a contaminated food. Its genome consists of 2,886,017 bp double-stranded DNA with a GC content of 32.8%. It is predicted to contain 2,728 open reading frames, 57 tRNAs, and 6 rRNA operons, including 1 additional 5S rRNA gene. Comparative phylogenetic tree analysis of 40 complete S. aureus genome sequences using average nucleotide identity (ANI) revealed that strain FORC_001 belonged to Group I. The closest phylogenetic match was S. aureus MRSA252, according to a whole-genome ANI (99.87%), suggesting that they might share a common ancestor. Comparative genome analysis of FORC_001 and MRSA252 revealed two non-homologous regions: Regions I and II. The presence of various antibiotic resistance genes, including the SCCmec cluster in Region I of MRSA252, suggests that this strain might have acquired the SCCmec cluster to adapt to specific environments containing methicillin. Region II of both genomes contains prophage regions but their DNA sequence identity is very low, suggesting that the prophages might differ. This is the first report of the complete genome sequence of S. aureus isolated from a real foodborne outbreak in South Korea. This report would be helpful to extend our understanding about the genome, general characteristics, and virulence factors of S. aureus for further studies of pathogenesis, rapid detection, and epidemiological investigation in foodborne outbreak.

Background Antheraea yamamai, also known as the Japanese oak silk moth, is a wild species of silk moth. Silk produced by A. yamamai, referred to as tensan silk, shows different characteristics such as thickness, compressive elasticity, and chemical resistance compared with common silk produced from the domesticated silkworm, Bombyx mori. Its unique characteristics have led to its use in many research fields including biotechnology and medical science, and the scientific as well as economic importance of the wild silk moth continues to gradually increase. However, no genomic information for the wild silk moth, including A. yamamai, is currently available. Findings In order to construct the A. yamamai genome, a total of 147G base pairs using Illumina and Pacbio sequencing platforms were generated, providing 210-fold coverage based on the 700-Mb estimated genome size of A. yamamai. The assembled genome of A. yamamai was 656 Mb (>2 kb) with 3675 scaffolds, and the N50 length of assembly was 739 Kb with a 34.07% GC ratio. Identified repeat elements covered 37.33% of the total genome, and the completeness of the constructed genome assembly was estimated to be 96.7% by Benchmarking Universal Single-Copy Orthologs v2 analysis. A total of 15 481 genes were identified using Evidence Modeler based on the gene prediction results obtained from 3 different methods (ab initio, RNA-seq-based, known-gene-based) and manual curation. Conclusions Here we present the genome sequence of A. yamamai, the first genome sequence of the wild silk moth. These results provide valuable genomic information, which will help enrich our understanding of the molecular mechanisms relating to not only specific phenotypes such as wild silk itself but also the genomic evolution of Saturniidae.

Background The history of African indigenous cattle and their adaptation to environmental and human selection pressure is at the root of their remarkable diversity. Characterization of this diversity is an essential step towards understanding the genomic basis of productivity and adaptation to survival under African farming systems. Results We analyze patterns of African cattle genetic variation by sequencing 48 genomes from five indigenous populations and comparing them to the genomes of 53 commercial taurine breeds. We find the highest genetic diversity among African zebu and sanga cattle. Our search for genomic regions under selection reveals signatures of selection for environmental adaptive traits. In particular, we identify signatures of selection including genes and/or pathways controlling anemia and feeding behavior in the trypanotolerant N’Dama, coat color and horn development in Ankole, and heat tolerance and tick resistance across African cattle especially in zebu breeds. Conclusions Our findings unravel at the genome-wide level, the unique adaptive diversity of African cattle while emphasizing the opportunities for sustainable improvement of livestock productivity on the continent. Electronic supplementary material The online version of this article (doi:10.1186/s13059-017-1153-y) contains supplementary material, which is available to authorized users.

Although growth rate is one of the main economic traits of concern in pig production, there is limited knowledge on its epigenetic regulation, such as DNA methylation. In this study, we conducted methyl-CpG binding domain protein-enriched genome sequencing (MBD-seq) to compare genome-wide DNA methylation profile of small intestine and liver tissue between fast- and slow-growing weaning piglets. The genome-wide methylation pattern between the two different growing groups showed similar proportion of CpG (regions of DNA where a cytosine nucleotide occurs next to a guanine nucleotide in the linear sequence) coverage, genomic regions, and gene regions. Differentially methylated regions and genes were also identified for downstream analysis. In canonical pathway analysis using differentially methylated genes, pathways (triacylglycerol pathway, some cell cycle related pathways, and insulin receptor signaling pathway) expected to be related to growth rate were enriched in the two organ tissues. Differentially methylated genes were also organized in gene networks related to the cellular development, growth, and carbohydrate metabolism. Even though further study is required, the result of this study may contribute to the understanding of epigenetic regulation in pig growth.

Korean native cattle, known as Hanwoo, have been raised in the Korean Peninsula since 2000 B.C. for use as a draft animal. However, Hanwoo now have an important position in the Korean livestock industry as a meat source. Therefore, the breeding and selection of Hanwoo are crucial for the industry. Although many researchers have studied the genetic architecture of Hanwoo, no well-established Hanwoo-related databases exist. In order to better understand the genetic contents of Hanwoo, an integrated database is necessary. We constructed a genetic variants database including annotation information. HanwooGDB ( http://hanwoogdb.snu.ac.kr ) provides genetic variants (SNPs, INDELs, CNVs) in the Hanwoo genome produced by Next Generation Sequencing data collected from 23 cattle. The identified SNPs were integrated with SNP chip data and annotation information for checking the concordance of position of each SNP and inferring functional aspects. This database provides browsers to understand and visualize the comprehensive information of these variants and allows users to download data according to their preference from this database without limitation. This database will contribute to genetic research and development of Hanwoo breeding strategies.

Genetics is important for breeding and selection of horses but there is a lack of well-established horse-related browsers or databases. In order to better understand horses, more variants and other integrated information are needed. Thus, we construct a horse genomic variants database including expression and other information. Horse Single Nucleotide Polymorphism and Expression Database (HSDB) (http://snugenome2.snu.ac.kr/HSDB) provides the number of unexplored genomic variants still remaining to be identified in the horse genome including rare variants by using population genome sequences of eighteen horses and RNA-seq of four horses. The identified single nucleotide polymorphisms (SNPs) were confirmed by comparing them with SNP chip data and variants of RNA-seq, which showed a concordance level of 99.02% and 96.6%, respectively. Moreover, the database provides the genomic variants with their corresponding transcriptional profiles from the same individuals to help understand the functional aspects of these variants. The database will contribute to genetic improvement and breeding strategies of Thoroughbreds.

Congenital cataracts can occur as a non-syndromic isolated ocular disease or as a part of genetic syndromes accompanied by a multi-systemic disease. Approximately 50% of all congenital cataract cases have a heterogeneous genetic basis. Here, we describe three generations of a family with an autosomal dominant inheritance pattern and common complex phenotypes, including bilateral congenital cataracts, short stature, macrocephaly, and minor skeletal anomalies. We did not find any chromosomal aberrations or gene copy number abnormalities using conventional genetic tests; accordingly, we conducted whole-exome sequencing (WES) to identify disease-causing genetic alterations in this family. Based on family WES data, we identified a novel BRD4 missense mutation as a candidate causal variant and performed cell-based experiments by ablation of endogenous BRD4 expression in human lens epithelial cells. The protein expression levels of connexin 43, p62, LC3BII, and p53 differed significantly between control cells and cells in which endogenous BRD4 expression was inhibited. We inferred that a BRD4 missense mutation was the likely disease-causing mutation in this family. Our findings may improve the molecular diagnosis of congenital cataracts and support the use of WES to clarify the genetic basis of complex diseases.

Congenital cataracts can occur as a non-syndromic isolated ocular disease or as a part of genetic syndromes accompanied by a multi-systemic disease. Approximately 50% of all congenital cataract cases have a heterogeneous genetic basis. Here, we describe three generations of a family with an autosomal dominant inheritance pattern and common complex phenotypes, including bilateral congenital cataracts, short stature, macrocephaly, and minor skeletal anomalies. We did not find any chromosomal aberrations or gene copy number abnormalities using conventional genetic tests; accordingly, we conducted whole-exome sequencing (WES) to identify disease-causing genetic alterations in this family. Based on family WES data, we identified a novel BRD4 missense mutation as a candidate causal variant and performed cell-based experiments by ablation of endogenous BRD4 expression in human lens epithelial cells. The protein expression levels of connexin 43, p62, LC3BII, and p53 differed significantly between control cells and cells in which endogenous BRD4 expression was inhibited. We inferred that a BRD4 missense mutation was the likely disease-causing mutation in this family. Our findings may improve the molecular diagnosis of congenital cataracts and support the use of WES to clarify the genetic basis of complex diseases.

Supplemental data including Note S1, Tables S1-S7 and Figures S1-S12. (DOCX 14235Â kb)

Background: Abalones are large marine snails in the family Haliotidae and the genus Haliotis belonging to the class Gastropoda of the phylum Mollusca. The family Haliotidae contains only one genus, Haliotis, and this single genus is known to contain several species of abalone. With 18 additional subspecies, the most comprehensive treatment of Haliotidae considers 56 species valid [1]. Abalone is an economically important fishery and aquaculture animal that is considered a highly prized seafood delicacy. The total global supply of abalone has increased 5-fold since the 1970s and farm production increased explosively from 50 mt to 103 464 mt in the past 40 years. Additionally, researchers have recently focused on abalone given their reported tumor suppression effect. However, despite the valuable features of this marine animal, no genomic information is available for the Haliotidae family and related research is still limited. To construct the H. discus hannai genome, a total of 580-G base pairs using Illumina and Pacbio platforms were generated with 322-fold coverage based on the 1.8-Gb estimated genome size of H. discus hannai using flow cytometry. The final genome assembly consisted of 1.86 Gb with 35 450 scaffolds (>2 kb). GC content level was 40.51%, and the N50 length of assembled scaffolds was 211 kb. We identified 29 449 genes using Evidence Modeler based on the gene information from ab initio prediction, protein homology with known genes, and transcriptome evidence of RNA-seq. Here we present the first Haliotidae genome, H. discus hannai, with sequencing data, assembly, and gene annotation information. This will be helpful for resolving the lack of genomic information in the Haliotidae family as well as providing more opportunities for understanding gastropod evolution.