Search

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Lactobacillus plantarum is found in various environmental niches such as in the gastrointestinal tract of an animal host or a fermented food. This species isolated from a certain environment is known to possess a variety of properties according to inhabited environment's adaptation. However, a causal relationship of a genetic factor and phenotype affected by a specific environment has not been systematically comprehended. L. plantarum GB-LP3 strain was isolated from Korean traditional fermented vegetable and the whole genome of GB-LP3 was sequenced. Comparative genome analysis of GB-LP3, with other 14 L. plantarum strains, was conducted. In addition, genomic island regions were investigated. The assembled whole GB-LP3 genome contained a single circular chromosome of 3,206,111bp with the GC content of 44.7%. In the phylogenetic tree analysis, GB-LP3 was in the closest distance from ZJ316. The genomes of GB-LP3 and ZJ316 have the high level of synteny. Functional genes that are related to prophage, bacteriocin, and quorum sensing were found through comparative genomic analysis with ZJ316 and investigation of genomic islands. dN/dS analysis identified that the gene coding for phosphonate ABC transporter ATP-binding protein is evolutionarily accelerated in GB-LP3. Our study found that potential candidate genes that are affected by environmental adaptation in Korea traditional fermented vegetable.

Non-alcoholic fatty liver disease (NAFLD) is one of the most frequent causes of liver disease and its prevalence is a serious and growing clinical problem. Caloric restriction (CR) is commonly recommended for improvement of obesity-related diseases such as NAFLD. However, the effects of CR on hepatic metabolism remain unknown. We investigated the effects of CR on metabolic dysfunction in the liver of obese diabetic db/db mice. We found that CR of db/db mice reverted insulin resistance, hepatic steatosis, body weight and adiposity to those of db/m mice. 1H-NMR- and UPLC-QTOF-MS-based metabolite profiling data showed significant metabolic alterations related to lipogenesis, ketogenesis, and inflammation in db/db mice. Moreover, western blot analysis showed that lipogenesis pathway enzymes in the liver of db/db mice were reduced by CR. In addition, CR reversed ketogenesis pathway enzymes and the enhanced autophagy, mitochondrial biogenesis, collagen deposition and endoplasmic reticulum stress in db/db mice. In particular, hepatic inflammation-related proteins including lipocalin-2 in db/db mice were attenuated by CR. Hepatic metabolomic studies yielded multiple pathological mechanisms of NAFLD. Also, these findings showed that CR has a therapeutic effect by attenuating the deleterious effects of obesity and diabetes-induced multiple complications.

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.

We have revealed in a porcine skin injury model that eosinophil recruitment was dose-dependently enhanced by a single high-dose irradiation. In this study, we investigated the underlying mechanism of eosinophil-associated skin fibrosis and the effect of high-dose-per-fraction radiation. The dorsal skin of a mini-pig was divided into two sections containing 4-cm2 fields that were irradiated with 30 Gy in a single fraction or 5 fractions and biopsied regularly over 14 weeks. Eosinophil-related Th2 cytokines such as interleukin (IL)-4, IL-5, and C–C motif chemokine-11 (CCL11/eotaxin) were evaluated by quantitative real-time PCR. RNA-sequencing using 30 Gy-irradiated mouse skin and functional assays in a co-culture system of THP-1 and irradiated-human umbilical vein endothelial cells (HUVECs) were performed to investigate the mechanism of eosinophil-mediated radiation fibrosis. Single high-dose-per-fraction irradiation caused pronounced eosinophil accumulation, increased profibrotic factors collagen and transforming growth factor-β, enhanced production of eosinophil-related cytokines including IL-4, IL-5, CCL11, IL-13, and IL-33, and reduced vessels compared with 5-fraction irradiation. IL-33 notably increased in pig and mouse skin vessels after single high-dose irradiation of 30 Gy, as well as in irradiated HUVECs following 12 Gy. Blocking IL-33 suppressed the migration ability of THP-1 cells and cytokine secretion in a co-culture system of THP-1 cells and irradiated HUVECs. Hence, high-dose-per-fraction irradiation appears to enhance eosinophil-mediated fibrotic responses, and IL-33 may be a key molecule operating in eosinophil-mediated fibrosis in high-dose-per fraction irradiated skin.

Diabetes-induced cognitive decline has been recognized in human patients of type 2 diabetes mellitus and mouse model of obesity, but the underlying mechanisms or therapeutic targets are not clearly identified. We investigated the effect of caloric restriction on diabetes-induced memory deficits and searched a molecular mechanism of caloric restriction-mediated neuroprotection. C57BL/6 mice were fed a high-fat diet for 40 weeks and RNA-seq analysis was performed in the hippocampus of high-fat diet-fed mice. To investigate caloric restriction effect on differential expression of genes, mice were fed high-fat diet for 20 weeks and continued on high-fat diet or subjected to caloric restriction (2 g/day) for 12 weeks. High-fat diet-fed mice exhibited insulin resistance, glial activation, blood–brain barrier leakage, and memory deficits, in that we identified neurogranin, a down-regulated gene in high-fat diet-fed mice using RNA-seq analysis; neurogranin regulates Ca2+/calmodulin-dependent synaptic function. Caloric restriction increased insulin sensitivity, reduced high-fat diet-induced blood–brain barrier leakage and glial activation, and improved memory deficit. Furthermore, caloric restriction reversed high-fat diet-induced expression of neurogranin and the activation of Ca2+/calmodulin-dependent protein kinase II and calpain as well as the downstream effectors. Our results suggest that neurogranin is an important factor of high-fat diet-induced memory deficits on which caloric restriction has a therapeutic effect by regulating neurogranin-associated calcium signaling.

Background Natural and artificial selection following domestication has led to the existence of more than a hundred pig breeds, as well as incredible variation in phenotypic traits. Berkshire pigs are regarded as having superior meat quality compared to other breeds. As the meat production industry seeks selective breeding approaches to improve profitable traits such as meat quality, information about genetic determinants of these traits is in high demand. However, most of the studies have been performed using trained sensory panel analysis without investigating the underlying genetic factors. Here we investigate the relationship between genomic composition and this phenotypic trait by scanning for signatures of positive selection in whole-genome sequencing data. Results We generated genomes of 10 Berkshire pigs at a total of 100.6 coverage depth, using the Illumina Hiseq2000 platform. Along with the genomes of 11 Landrace and 13 Yorkshire pigs, we identified genomic variants of 18.9 million SNVs and 3.4 million Indels in the mapped regions. We identified several associated genes related to lipid metabolism, intramuscular fatty acid deposition, and muscle fiber type which attribute to pork quality (TG, FABP1, AKIRIN2, GLP2R, TGFBR3, JPH3, ICAM2, and ERN1) by applying between population statistical tests (XP-EHH and XP-CLR). A statistical enrichment test was also conducted to detect breed specific genetic variation. In addition, de novo short sequence read assembly strategy identified several candidate genes (SLC25A14, IGF1, PI4KA, CACNA1A) as also contributing to lipid metabolism. Conclusions Results revealed several candidate genes involved in Berkshire meat quality; most of these genes are involved in lipid metabolism and intramuscular fat deposition. These results can provide a basis for future research on the genomic characteristics of Berkshire pigs. Electronic supplementary material The online version of this article (doi:10.1186/s12863-015-0265-1) contains supplementary material, which is available to authorized users.

Adipocytes mainly function as energy storage and endocrine cells. Adipose tissues showed the biological and genetic difference based on their depots. The difference of adipocytes between depots might be influenced by the inherent genetic programing for adipogenesis. We used RNA-seq technique to investigate the transcriptomes in 3 adipose tissues of omental (O), subcutaneous (S) and intramuscular (I) fats in cattle. Sequence reads were obtained from Illumina HiSeq2000 and mapped to the bovine genome using Tophat2. Differentially expressed genes (DEG) between adipose tissues were detected by EdgeR. We identified 5797, 2156, and 5455 DEGs in the comparison between OI, OS, and IS respectively and also found 5657 DEGs in the comparison between the intramuscular and the combined omental and subcutaneous fats (C) (FDR

Since being domesticated about 10,000-12,000 years ago, domestic pigs (Sus scrofa domesticus) have been selected for traits of economic importance, in particular large body size. However, Yucatan miniature pigs have been selected for small body size to withstand high temperature environment and for laboratory use. This renders the Yucatan miniature pig a valuable model for understanding the evolution of body size. We investigate the genetic signature for selection of body size in the Yucatan miniature pig. Phylogenetic distance of Yucatan miniature pig was compared to other large swine breeds (Yorkshire, Landrace, Duroc and wild boar). By estimating the XP-EHH statistic using re-sequencing data derived from 70 pigs, we were able to unravel the signatures of selection of body size. We found that both selections at the level of organism, and at the cellular level have occurred. Selection at the higher levels include feed intake, regulation of body weight and increase in mass while selection at the molecular level includes cell cycle and cell proliferation. Positively selected genes probed by XP-EHH may provide insight into the docile character and innate immunity as well as body size of Yucatan miniature pig.

The Korean Native Chicken (KNC) is an important endemic biological resource in Korea. While numerous studies have been conducted exploring this breed, none have used next-generation sequencing to identify its specific genomic features. We sequenced five strains of KNC and identified 10.9 million SNVs and 1.3 million InDels. Through the analysis, we found that the highly variable region common to all 5 strains had genes like PCHD15, CISD1, PIK3C2A, and NUCB2 that might be related to the phenotypic traits of the chicken such as auditory sense, growth rate and egg traits. In addition, we assembled unaligned reads that could not be mapped to the reference genome. By assembling the unaligned reads, we were able to present genomic sequences characteristic to the KNC. Based on this, we also identified genes related to the olfactory receptors and antigen that are common to all 5 strains. Finally, through the reconstructed mitochondrial genome sequences, we performed phylogenomic analysis and elucidated the maternal origin of the artificially restored KNC. Our results revealed that the KNC has multiple maternal origins which are in agreement with Korea's history of chicken breed imports. The results presented here provide a valuable basis for future research on genomic features of KNC and further understanding of KNC's origin.

Background Whales have captivated the human imagination for millennia. These incredible cetaceans are the only mammals that have adapted to life in the open oceans and have been a source of human food, fuel and tools around the globe. The transition from land to water has led to various aquatic specializations related to hairless skin and ability to regulate their body temperature in cold water. Results We present four common minke whale (Balaenoptera acutorostrata) genomes with depth of ×13 ~ ×17 coverage and perform resequencing technology without a reference sequence. Our results indicated the time to the most recent common ancestors of common minke whales to be about 2.3574 (95% HPD, 1.1521 – 3.9212) million years ago. Further, we found that genes associated with epilation and tooth-development showed signatures of positive selection, supporting the morphological uniqueness of whales. Conclusions This whole-genome sequencing offers a chance to better understand the evolutionary journey of one of the largest mammals on earth. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-1213-1) contains supplementary material, which is available to authorized users.