Search

Abstract NA4/NA6, an intermediate degradation product of β-agarase, is a high value-added product with anticancer, anti-obesity, and anti-diabetic effects. Therefore, a method that enables the efficient production of NA4/NA6 would be useful from economic and medical perspectives. In this study, we aimed to generate a Streptomyces coelicolor A3(2) mutant M22-2C43 that produces NA4/NA6 as a final product; this method serves as a more efficient alternative to the enzymatic conversion of β-agarase for the generation of these products. The M22-2C43 strain was generated through two rounds of mutagenesis and screening for increased β-agarase activity and effective production of NA4/NA6. We assembled the complete genomes of two mutants, M22 and M22-2C43, which were identified following a two-round screening. Large and small genetic changes were found in these two mutants, including the loss of two plasmids present in wild-type S. coelicolor A3(2) and chromosome circularization of mutant M22-2C43. These findings suggest that mutant M22-2C43 can produce NA4/NA6 as a degradation product due to functional inactivation of the dagB gene through a point mutation (G474A), ultimately preventing further degradation of NA4/NA6 to NA2. To our knowledge, this is the first report of a microbial strain that can effectively produce NA4/NA6 as the main degradation product of β-agarase, opening the door for the use of this species for the large-scale production of this valuable product.

Abstract In response to the COVID-19 pandemic, different types of vaccines, such as inactive, live-attenuated, messenger RNA (mRNA), and protein subunit, have been developed against SARS-CoV-2. This has unintentionally created a unique scenario where heterologous prime-boost vaccination against a single virus has been administered to 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 encoding the hemagglutinin (HA) glycoprotein of the influenza virus 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 the mRNA-HA priming/protein-HA boosting (R-P) regimen than those induced by reverse immunization (protein-HA priming/mRNA-HA boosting, P-R). After the viral challenge, the R-P group showed lower virus loads and less inflammation in the lungs than the P-R group did. Transcriptome analysis revealed that the 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 direct desired immune responses. This study demonstrates the potential of a heterologous vaccination strategy using mRNA and protein vaccine platforms against viral infection.

Dermacoccus barathri is the first reported pathogen within the Dermacoccus genus to cause a catheter-related bloodstream infection, which occurred in 2015. In this study, the complete genome assembly of Dermacoccus barathri was constructed, and the complete genome of Dermacoccus barathri FBCC-B549 consists of a single chromosome (3,137,745 bp) without plasmids. The constructed genome of D. barathri was compared with those of two closely related species within the Dermacoccus genus. D. barathri exhibited a pattern similar to Dermacoccus abyssi in terms of gene clusters and synteny analysis. Contrary to previous studies, biosynthetic gene cluster (BGC) analysis for predicting secondary metabolites revealed the presence of the LAP biosynthesis pathway in the complete genome of D. barathri, predicting the potential synthesis of the secondary metabolite plantazolicin. Furthermore, an analysis to investigate the potential pathogenicity of D. barathri did not reveal any antibiotic resistance genes; however, nine virulence factors were identified in the Virulence Factor Database (VFDB). According to these matching results in the VFDB, despite identifying a few factors involved in biofilm formation, further research is required to determine the actual impact of D. barathri on pathogenicity. The complete genome of D. barathri is expected to serve as a valuable resource for future studies on D. barathri, which currently lack sufficient genomic sequence information.

Most pancreatic cancers are pancreatic ductal adenocarcinomas. This is an extremely lethal disease with poor prognosis and almost no treatment choices. Considering the profound role of the pancreas in the human body, malfunction of this organ can significantly affect quality of life. Although multiple metabolic pathways are altered in cancer cells, certain metabolic gene signatures may be critical for immunotherapy. The reprogrammed metabolism of glucose, amino acids, and lipids can nourish the tumor microenvironment (TME). Previous studies have also shown that reprogrammed metabolism influences immune responses. Tumor-infiltrating immune cells in the TME can adapt their metabolism to blunt the immune system, leading to immunosuppression and tumor progression. The identification of metabolism-related genes (MRGs) associated with immune reactions in pancreatic cancer may lead to improved treatments. This review highlights the characteristics of MRGs in pancreatic cancer and suggests that enhanced anti-cancer therapies could be used to overcome resistance to immunotherapy.

The fall armyworm (FAW), Spodoptera frugiperda (Lepidoptera: Noctuidae), is native to tropical and subtropical regions of the Western Hemisphere, but is now regularly appearing in crop fields across South Korea, particularly in corn fields. Therefore, it is crucial to promptly and accurately identify the presence of FAW in crop fields to effectively eradicate it as a regulated quarantine species. We developed a loop-mediated isothermal amplification (LAMP) assay, which allows for rapid in-filed identification. To develop the LAMP assay, we selected FAW-specific genomic regions from the whole-genome sequences of one FAW and 13 other lepidopteran species and validated five primer sets that consistently produced positive reactions in ten FAW samples collected from eight different locations in four countries. The assay successfully identified FAW in a maximum of 45 min, starting from crude DNA extraction (~15 min) to diagnosis (30 min) from the following samples, which were deposited outdoors for 30 days: a 1st-instar larva, an adult leg, an adult antenna, and 1/16 and 1/8 of an adult thorax. The five assays can be used selectively or in combination to cross-check and provide further confidence in the in-field diagnosis of FAW.

Measurement(s) Larval period • Cocoon yield • Single cocoon weight • Egg color • Cocoon shape • Number of cocoons per liter • Voltinism • Moltinism • Cocoon color• Larva pattern Technology Type(s) Clock Device • scale • Eye • amount Factor Type(s) Artificial feed • Target Season Sample Characteristic - Organism Bombyx mori Sample Characteristic - Environment house Sample Characteristic - Location South Korea

Summary: Insights into the evolution of non-model organisms are limited by the lack of reference genomes of high accuracy, completeness, and contiguity. Here, we present a chromosome-level, karyotype-validated reference genome and pangenome for the barn swallow (Hirundo rustica). We complement these resources with a reference-free multialignment of the reference genome with other bird genomes and with the most comprehensive catalog of genetic markers for the barn swallow. We identify potentially conserved and accelerated genes using the multialignment and estimate genome-wide linkage disequilibrium using the catalog. We use the pangenome to infer core and accessory genes and to detect variants using it as a reference. Overall, these resources will foster population genomics studies in the barn swallow, enable detection of candidate genes in comparative genomics studies, and help reduce bias toward a single reference genome.

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

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

ABSTRACT 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. IMPORTANCE Metabarcoding is a powerful tool to investigate community diversity in an economic and efficient way by amplifying a specific gene marker region. With the advancement of long-read sequencing technologies, the field of metabarcoding has entered a new phase. The technologies have brought a need for development in several areas, including new markers that long-read can cover, database for the markers, tools that reflect long-read characteristics, and compatibility with downstream analysis tools. By constructing MIrROR, we met the need for a database and tools for the 16S-ITS-23S rRNA operon region, which has recently been shown to have sufficient resolution at the species level. Bacterial community analysis using the 16S-ITS-23S rRNA operon region with MIrROR will provide new insights from various research fields.

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.

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.

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

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.

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.

Ribosomal RNA is an indispensable molecule in living organisms that plays an essential role in protein synthesis. Especially in bacteria, 16S, 23S, and 5S rRNAs are usually co-transcribed as operons. Despite the positive effects of rRNA co-transcription on growth and reproduction rate, a recent study revealed that bacteria with unlinked rRNA operons are more widespread than expected. However, it is still unclear why the rRNA operon is broken. Here, we explored rRNA operon linkage status in 15,898 bacterial genomes and investigated whether they have common features or lifestyles; 574 genomes were found to have unlinked rRNA operons and tended to be phylogenetically conserved. Most of them were symbionts and showed enhanced symbiotic genomic features such as reduced genome size and high adenine–thymine (AT) content. In an eggNOG-mapper analysis, they were also found to have significantly fewer genes than rRNA operon-linked bacteria in the “transcription” and “energy production and conversion in metabolism” categories. These genomes also tend to decrease RNases related to the synthesis of ribosomes and tRNA processing. Based on these results, the disruption of the rRNA operon seems to be one of the tendencies associated with the characteristics of bacteria requiring a low dynamic range.

Recent comparative genomics studies have suggested that horizontal gene transfer (HGT) is one of the major processes in bacterial evolution. In this study, HGT events of 64 Chlamydia strains were investigated based on the pipeline employed in HGTree database constructed in our recent study. Tree reconciliation method was applied in order to calculate feasible HGT events. Following initial detection and an evaluation procedure, evidence of the HGT was identified in 548 gene families including 42 gene families transferred from outside of Chlamydiae phylum with high reliability. The donor species of inter-phylum HGT consists of 12 different bacterial and archaeal phyla, suggesting that Chlamydia might have even more various host range than in previous reports. In addition, each species of Chlamydia showed varying preference towards HGT, and genes engaged in HGT within Chlamydia and between other species showed different functional distribution. Also, examination of individual gene flows of niche-specific genes suggested that many of such genes are transferred mainly within Chlamydia genus. Our results uncovered novel features of HGT acting on Chlamydia genome evolution, and it would be also strong evidence that HGT is an ongoing process for intracellular pathogens. We expect that the results provide more insight into lineage- and niche-specific adaptations regarding their infectivity and pathogenicity.

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.

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.

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.

[This corrects the article DOI: 10.1371/journal.pone.0169226.].

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.

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